A More Accurate Generalized Gradient Approximation for Solids

Zhigang Wu,Ronald E. Cohen
DOI: https://doi.org/10.1103/PhysRevB.73.235116
2006-01-25
Abstract:We present a new nonempirical density functional generalized gradient approximation (GGA) that gives significant improvements for lattice constants, crystal structures, and metal surface energies over the most popular Perdew-Burke-Ernzerhof (PBE) GGA. The new functional is based on a diffuse radial cutoff for the exchange-hole in real space, and the analytic gradient expansion of the exchange energy for small gradients. There are no adjustable parameters, the constraining conditions of PBE are maintained, and the functional is easily implemented in existing codes.
Materials Science
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to improve the computational accuracy of the lattice constant, crystal structure, and metal surface energy of solid materials. Specifically, the author proposes a new non - empirical generalized gradient approximation (GGA) method, aiming to significantly improve the performance of the existing most popular Perdew - Burke - Ernzerhof (PBE) GGA method in these aspects. ### Background problems 1. **Lattice constant and crystal structure**: - LSD (local spin density approximation) usually underestimates the lattice constant \(a_0\) by 1 - 3%, while GGA methods such as PBE tend to overestimate the lattice constant by 1 - 2%. - For some materials, such as ferroelectrics, this deviation may lead to serious prediction errors, especially in volume - sensitive properties. 2. **Metal surface energy**: - The prediction of metal surface energy by PBE has deviations, especially when dealing with extended systems such as metal surfaces, and its exchange - correlation hole may exhibit a diffusive long - tail behavior. ### The new method proposed in the paper The author proposes a new GGA method, which is improved based on the following points: 1. **Diffusive radial cutoff**: - A diffusive radial cutoff factor is used to describe the exchange - correlation hole instead of the traditional sharp cutoff. This can describe the change in electron density in solids more accurately. 2. **Analytical gradient expansion**: - An analytical gradient expansion under a small gradient is introduced to describe the exchange energy more precisely. 3. **No adjustable parameters**: - The new method does not contain any adjustable parameters, ensuring the universality and reliability of the method. 4. **Maintaining constraint conditions**: - The constraint conditions of PBE are maintained, making the new method easy to implement in existing codes. ### Mathematical expression The core of the new method lies in defining a new exchange enhancement factor \(F_X\), whose form is as follows: \[ E_{\text{GGA}}^X=\int n\epsilon_{\text{unif}}^X(n)F_X(s)d^3r \] where \(s\) is the normalized gradient, defined as: \[ s = \frac{|\nabla n|}{[2(3\pi^2)^{1/3}n^{4/3}]} \] The new expression of \(F_X\) is: \[ x=\frac{10}{81}s^2+(\mu - \frac{10}{81})s^2\exp(-s^2)+\ln(1 + cs^4) \] where the parameter \(c\) is set to recover the fourth - order parameter at small \(s\), and its specific value is: \[ c=\frac{146}{2025}(\frac{2}{3})^2-\frac{73}{405}\cdot\frac{2}{3}+(\mu - \frac{10}{81}) = 0.0079325 \] ### Results and verification By calculating the equilibrium lattice constant, binding energy, and metal surface energy of various solid materials, the author proves that the new method has significant improvements in these aspects. In particular, for ferroelectric materials, the new method can predict volume, strain, and atomic displacement more accurately. ### Summary The main contribution of this paper is to provide a new GGA method, which significantly improves the computational accuracy of the properties of solid materials, especially for materials that require high - precision prediction, such as ferroelectrics. This method not only improves the prediction of lattice constant and crystal structure but also performs well in metal surface energy and other aspects.