Comparison of performance of van der Waals-corrected exchange-correlation functionals for interlayer interaction in graphene and hexagonal boron nitride

Irina V. Lebedeva,Alexander V. Lebedev,Andrey M. Popov,Andrey A. Knizhnik

DOI: https://doi.org/10.1016/j.commatsci.2016.11.011

2018-08-21

Abstract:Exchange-correlation functionals with corrections for van der Waals interactions (PBE-D2, PBE-D3, PBE-D3(BJ), PBE-TS, optPBE-vdW and vdW-DF2) are tested for graphene and hexagonal boron nitride, both in the form of bulk and bilayer. The characteristics of the potential energy surface, such as the barrier to relative sliding of the layers and magnitude of corrugation, and physically measurable properties associated with relative in-plane and out-of-plane motion of the layers including the shear modulus and modulus for axial compression, shear mode frequency and frequency of out-of-plane vibrations are considered. The PBE-D3(BJ) functional gives the best results for the stackings of hexagonal boron nitride and graphite that are known to be ground-state from the experimental studies. However, it fails to describe the order of metastable states of boron nitride in energy. The PBE-D3 and vdW-DF2 functionals, which reproduce this order correctly, are identified as the optimal choice for general studies. The vdW-DF2 functional is preferred for evaluation of the modulus for axial compression and frequency of out-of-plane vibrations, while the PBE-D3 functional is somewhat more accurate in calculations of the shear modulus and shear mode frequency. The best description of the latter properties, however, is achieved also using the vdW-DF2 functional combined with consideration of the experimental interlayer distance. In the specific case of graphene, the PBE-D2 functional works very well and can be further improved by adjustment of the parameters.

Mesoscale and Nanoscale Physics

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to evaluate and compare the performance of different van der Waals (vdW) -corrected exchange - correlation functionals in describing the inter - layer interactions of graphene and hexagonal boron nitride (h - BN). Specifically, the study focuses on the following aspects:

1. **Inter - layer distance and binding energy**: These are the basic properties of materials and directly affect the manifestation of other physical properties. Experimentally, the inter - layer distances of graphene and hexagonal boron nitride are known very precisely. Therefore, accurate prediction of these values is crucial for validating theoretical models.

2. **Potential - energy - surface characteristics**: Including relative sliding barriers, surface - roughness degrees, etc. These characteristics determine the relative motion behavior between layers, such as shear modulus and axial compressive modulus, as well as related vibration frequencies (e.g., in - plane and out - of - plane vibration modes). These parameters are very important for understanding the response of materials under mechanical stress.

3. **Energy order of stacking configurations**: Especially for hexagonal boron nitride, understanding the energy differences between different stacking methods (such as AA', AB1', etc.) helps to determine the most stable structure and explain the experimentally observed phenomena.

4. **Consistency with experimental data**: By comparing the calculation results with the existing experimental measurements, the effectiveness and reliability of each vdW correction method can be evaluated. This is not limited to the specific properties mentioned above, but also includes any information that can be obtained directly or indirectly through experiments.

In summary, this paper aims to find the best scheme for describing the long - range vdW interactions in two - dimensional materials such as graphene and hexagonal boron nitride by comprehensively testing multiple vdW - corrected density - functional - theory (DFT) methods. In addition, it also explores which functionals can provide more accurate results in specific application scenarios, such as functionals that perform excellently in calculating certain mechanical properties.