Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAssessing the Performance of Recent Density Functionals for Bulk Solids
142
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We assess the performance of recent density functionals for the exchange-correlation energy of a nonmolecular solid, by applying accurate calculations with the GAUSSIAN, BAND, and VASP codes to a test set of 24 solid metals and non-metals. The functionals tested are the modified Perdew-Burke-Ernzerhof generalized gradient approximation (PBEsol GGA), the second-order GGA (SOGGA), and the Armiento-Mattsson 2005 (AM05) GGA. For completeness, we also test more-standard functionals: the local density approximation, the original PBE GGA, and the Tao-Perdew-Staroverov-Scuseria (TPSS) meta-GGA. We find that the recent density functionals for solids reach a high accuracy for bulk properties (lattice constant and bulk modulus). For the cohesive energy, PBE is better than PBEsol overall, as expected, but PBEsol is actually better for the alkali metals and alkali halides. For fair comparison of calculated and experimental results, we consider the zero-point phonon and finite-temperature effects ignored by many workers. We show how Gaussian basis sets and inaccurate experimental reference data may affect the rating of the quality of the functionals. The results show that PBEsol and AM05 perform somewhat differently from each other for alkali metal, alkaline earth metal and alkali halide crystals (where the maximum value of the reduced density gradient is about 2), but perform very similarly for most of the other solids (where it is often about 1). Our explanation for this is consistent with the importance of exchange-correlation nonlocality in regions of core-valence overlap.
rate research
Read More
One of the standard generalized-gradient approximations (GGAs) in use in modern electronic-structure theory, PBE, and a recently proposed modification designed specifically for solids, PBEsol, are identified as particular members of a family of functionals taking their parameters from different properties of homogeneous or inhomogeneous electron liquids. Three further members of this family are constructed and tested, together with the original PBE and PBEsol, for atoms, molecules and solids. We find that PBE, in spite of its popularity in solid-state physics and quantum chemistry, is not always the best performing member of the family, and that PBEsol, in spite of having been constructed specifically for solids, is not the best for solids. The performance of GGAs for finite systems is found to sensitively depend on the choice of constraints steaming from infinite systems. Guidelines both for users and for developers of density functionals emerge from this work.
Kinetic energy (KE) approximations are key elements in orbital-free density functional theory. To date, the use of non-local functionals, possibly employing system dependent parameters, has been considered mandatory in order to obtain satisfactory accuracy for different solid-state systems, whereas semilocal approximations are generally regarded as unfit to this aim. Here, we show that instead properly constructed semilocal approximations, the Pauli-Gaussian (PG) KE functionals, especially at the Laplacian-level of theory, can indeed achieve similar accuracy as non-local functionals and can be accurate for both metals and semiconductors, without the need of system-dependent parameters.
The crystalline structure of ground-state orthorhombic SrRuO$_3$ is reproduced by applying hybrid density functional theory scheme to the functionals based on the revised generalized-gradient approximations for solid-state calculations. The amount of Hartree-Fock (HF) exchange energy is varied in the range of $5-20%$ in order to systematically ascertain the optimum value of HF mixing which in turn ensures the best correspondence to the experimental measurements. Such investigation allows to expand the set of tools that could be used for the efficient theoretical modelling of, for example, only recently stabilized phases of SrRuO$_3$.
Eleven density functionals are compared with regard to their performance for the lattice constants of solids. We consider standard functionals, such as the local-density approximation and the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA), as well as variations of PBE GGA, such as PBEsol and similar functionals, PBE-type functionals employing a tighter Lieb-Oxford bound, and combinations thereof. Several of these variations are proposed here for the first time. On a test set of 60 solids we perform a system-by-system analysis for selected functionals and a full statistical analysis for all of them. The impact of restoring the gradient expansion and of tightening the Lieb-Oxford bound is discussed, and confronted with previous results obtained from other codes, functionals or test sets. No functional is uniformly good for all investigated systems, but surprisingly, and pleasingly, the simplest possible modifications to PBE turn out to have the most beneficial effect on its performance. The atomization energy of molecules was also considered and on a testing set of six molecules, we found that the PBE functional is clearly the best, the others leading to strong overbinding.
A recent study of Mejia-Rodriguez and Trickey [Phys. Rev. A 96, 052512 (2017)] showed that the deorbitalization procedure (replacing the exact Kohn-Sham kinetic-energy density by an approximate orbital-free expression) applied to exchange-correlation functionals of the meta-generalized gradient approximation (MGGA) can lead to important changes in the results for molecular properties. For the present work, the deorbitalization of MGGA functionals is further investigated by considering various properties of solids. It is shown that depending on the MGGA, common orbital-free approximations to the kinetic-energy density can be sufficiently accurate for the lattice constant, bulk modulus, and cohesive energy. For the band gap, calculated with the modified Becke-Johnson MGGA potential, the deorbitalization has a larger impact on the results.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
