The ab initio density functional theory applied for spin-polarized calculations

Semilocal Meta-GGA Exchange-Correlation Approximation from Adiabatic Connection Formalism: Extent and Limitations

The incorporation of a strong-interaction regime within the approximate semilocal exchange-correlation functionals still remains a very challenging task for density functional theory. One of the promising attempts in this direction is the recently proposed adiabatic connection semilocal correlation (ACSC) approach [Constantin, L. A.; Phys. Rev. B 2019, 99, 085117] allowing one to construct the correlation energy functionals by interpolation of the high and low-density limits for the given semilocal approximation. The current study extends the ACSC method to the meta-generalized gradient approximations (meta-GGA) level of theory, providing some new insights in this context. As an example, we construct the correlation energy functional on the basis of the high- and low-density limits of the Tao-Perdew-Staroverov-Scuseria (TPSS) functional. Arose in this way, the TPSS-ACSC functional is one-electron self-interaction free and accurate for the strictly correlated and quasi-two-dimensional regimes. Based on simple examples, we show the advantages and disadvantages of ACSC semilocal functionals and provide some new guidelines for future developments in this context.

Self-Consistent Implementation of Kohn-Sham Adiabatic Connection Models with Improved Treatment of the Strong-Interaction Limit

Adiabatic connection models (ACMs), which interpolate between the limits of weak and strong interaction, are powerful tools to build accurate exchange–correlation functionals. If the exact weak-interaction expansion from the second-order perturbation theory is included, a self-consistent implementation of these functionals is challenging and still absent in the literature. In this work, we fill this gap by presenting a fully self-consistent-field (SCF) implementation of some popular ACM functionals. While using second-order perturbation theory at weak interactions, we have also introduced new generalized gradient approximations (GGAs), beyond the usual point-charge-plus-continuum model, for the first two leading terms at strong interactions, which are crucial to ensure robustness and reliability. We then assess the SCF–ACM functionals for molecular systems and for prototypical strong-correlation problems. We find that they perform well for both the total energy and the electronic density and that the impact of SCF orbitals is directly connected to the accuracy of the ACM functional form. For the H₂ dissociation, the SCF–ACM functionals yield significant improvements with respect to standard functionals also thanks to the use of the new GGAs for the strong-coupling functionals.

Solid-state performance of a meta-GGA screened hybrid density functional constructed from Pauli kinetic enhancement factor dependent semilocal exchange hole

The semilocal form of the exchange hole is highly useful in developing non-local range-separated hybrid density functionals for finite and extended systems. The way to construct the conventional exact exchange hole model is based on either the Taylor series expansion or the reverse engineering technique from the corresponding exchange energy functional. Although the latter technique is quite popular in context of generalized gradient approximation (GGA) functionals, the same for the meta-GGA functionals is not so much explored. Thus, in this study, we propose a reverse-engineered semilocal exchange hole of a meta-GGA functional, that depends only on the meta-GGA ingredient α (also known as the Pauli kinetic energy enhancement factor). The model is used subsequently to design the short-range-separated meta-GGA hybrid density functional. We show that the present method can be successfully applied for several challenging problems in the context of solids, especially for which the GGA based hybrid fails drastically. This assessment proves that the present functional is quite useful for materials sciences. Finally, we also use this method for several molecular test cases, where the results are also as comparative as its base semilocal functional.

Accurate density functional made more versatile

We propose a one-electron self-interaction-free correlation energy functional compatible with the order-of-limit problem-free Tao-Mo (TM) semilocal functional (regTM) [J. Tao and Y. Mo, Phys. Rev. Lett. 117, 073001 (2016) and Patra et al., J. Chem. Phys. 153, 184112 (2020)] to be used for general purpose condensed matter physics and quantum chemistry. The assessment of the proposed functional for large classes of condensed matter and chemical systems shows its improvement in most cases compared to the TM functional, e.g., when applied to the relative energy difference of MnO₂ polymorphs. In this respect, the present exchange-correction functional, which incorporates the TM technique of the exchange hole model combined with the slowly varying density correction, can achieve broad applicability, being able to solve difficult solid-state problems.

Random-Phase Approximation in Many-Body Noncovalent Systems: Methane in a Dodecahedral Water Cage

The many-body expansion (MBE) of energies of molecular clusters or solids offers a way to detect and analyze errors of theoretical methods that could go unnoticed if only the total energy of the system was considered. In this regard, the interaction between the methane molecule and its enclosing dodecahedral water cage, CH₄···(H₂O)₂₀, is a stringent test for approximate methods, including density functional theory (DFT) approximations. Hybrid and semilocal DFT approximations behave erratically for this system, with three- and four-body nonadditive terms having neither the correct sign nor magnitude. Here, we analyze to what extent these qualitative errors in different MBE contributions are conveyed to post-Kohn-Sham random-phase approximation (RPA), which uses approximate Kohn-Sham orbitals as its input. The results reveal a correlation between the quality of the DFT input states and the RPA results. Moreover, the renormalized singles energy (RSE) corrections play a crucial role in all orders of the many-body expansion. For dimers, RSE corrects the RPA underbinding for every tested Kohn-Sham model: generalized-gradient approximation (GGA), meta-GGA, (meta-)GGA hybrids, as well as the optimized effective potential at the correlated level. Remarkably, the inclusion of singles in RPA can also correct the wrong signs of three- and four-body nonadditive energies as well as mitigate the excessive higher-order contributions to the many-body expansion. The RPA errors are dominated by the contributions of compact clusters. As a workable method for large systems, we propose to replace those compact contributions with CCSD(T) energies and to sum up the remaining many-body contributions up to infinity with supermolecular or periodic RPA. As a demonstration of this approach, we show that for RPA(PBE0)+RSE it suffices to apply CCSD(T) to dimers and 30 compact, hydrogen-bonded trimers to get the methane-water cage interaction energy to within 1.6% of the reference value.

Generalizing Double-Hybrid Density Functionals: Impact of Higher-Order Perturbation Terms

Connections between the Görling-Levy (GL) perturbation theory and the parameters of double-hybrid (DH) density functional are established via adiabatic connection formalism. Moreover, we present a more general DH density functional theory, where the higher-order perturbation terms beyond the second-order GL2 one, such as GL3 and GL4, also contribute. It is shown that a class of DH functionals including previously proposed ones can be formed using the present construction. Based on the proposed formalism, we assess the performance of higher-order DH and long-range corrected DH formed on the Perdew-Burke-Ernzerhof (PBE) semilocal functional and second-order GL2 correlation energy. The underlying construction of DH functionals based on the generalized many-body perturbation approaches is physically appealing in terms of the development of the non-local forms using more accurate and sophisticated semilocal functionals.

Modified Interaction-Strength Interpolation Method as an Important Step toward Self-Consistent Calculations

The modified point charge plus continuum (mPC) model [ConstantinL. A.; Phys. Rev. B2019, 99, 085117] solves the important failures of the original counterpart, namely, the divergences when the reduced gradient of the density is large, such as in the tail of the density and in quasi-dimensional density regimes. The mPC allows us to define a modified interaction-strength interpolation (mISI) method inheriting these good features, which are important steps toward the full self-consistent treatment. Here, we provide an assessment of mISI for molecular systems (i.e., considering thermochemistry properties, correlation energies, vertical ionization potentials, and several noncovalent interactions), harmonium atoms, and functional derivatives in the strong-interaction limit. For all our tests, mISI provides a systematic improvement over the original ISI method. Semilocal approximations of the second-order Görling–Levy (GL2) perturbation theory are also considered in the mISI method, showing considerable worsening of the results. Possible further development of mISI is briefly discussed.

Unveiling the Physics Behind Hybrid Functionals

We show that accurate exchange-correlation hybrid functionals give very physically optimized effective-correlation potentials, capable of correctly describing the quantum oscillations of atoms and molecules. Based on this analysis and on understanding the error cancellation between semilocal exchange and correlation functionals, we propose a very simple, semilocal correlation potential model compatible with the exact exchange of density functional theory, which performs remarkably well for charge densities and orbital energies.