Self-consistent double-hybrid density-functional theory using the optimized-effective-potential method

A Critical Evaluation of the Hybrid KS DFT Functionals Based on the KS Exchange-Correlation Potentials

We have developed a critical methodology for the evaluation of the quality of hybrid exchange-correlation (XC) density functional approximations (DFAs) based on very fundamental quantities, i.e., Kohn-Sham (KS) XC potentials, self-consistent electron densities, first ionization potentials (IPs), and total energies. Since the XC potentials, the primary objects in the current study, are not directly accessible for the hybrids, we calculate them by inverting the KS electron densities. Utilizing this methodology, we tested 155 hybrid DFAs available in the LIBXC library using FCI and CCSD(T) methods as a reference. We have found that a group of functionals produces very decent XC potentials, mainly those with a large mixture of Hartree-Fock exchange. Moreover, the value of IP strongly depends on the XC potential quality. On the other hand, we show that the XC energy is dominated by functional-driven error, which in some cases leads to substantial errors in electronic densities. The study shows new directions for constructing more accurate XC functionals within the KS-DFT framework.

Chemical Potentials and the One-Electron Hamiltonian of the Second-Order Perturbation Theory from the Functional Derivative Approach

We develop a functional derivative approach to calculate the chemical potentials of second-order perturbation theory (MP2). In the functional derivative approach, the correlation part of the MP2 chemical potential, which is the derivative of the MP2 correlation energy with respect to the occupation number of frontier orbitals, is obtained from the chain rule via the noninteracting Green's function. First, the MP2 correlation energy is expressed in terms of the noninteracting Green's function, and its functional derivative to the noninteracting Green's function is the second-order self-energy. Then, the derivative of the noninteracting Green's function to the occupation number is obtained by including the orbital relaxation effect. We show that the MP2 chemical potentials obtained from the functional derivative approach agree with that obtained from the finite difference approach. The one-electron Hamiltonian, defined as the derivative of the MP2 energy with respect to the one particle density matrix, is also derived using the functional derivative approach, which can be used in the self-consistent calculations of MP2 and double-hybrid density functionals. The developed functional derivative approach is promising for calculating the chemical potentials and the one-electron Hamiltonian of approximate functionals and many-body perturbation approaches dependent explicitly on the noninteracting Green's function.

Vertical Ionization Potentials and Electron Affinities at the Double-Hybrid Density Functional Level

The double-hybrid (DH) time-dependent density functional theory is extended to vertical ionization potentials (VIPs) and electron affinities (VEAs). Utilizing the density fitting approximation, efficient implementations are presented for the genuine DH ansatz relying on the perturbative second-order correction, while an iterative analogue is also elaborated using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH approach. The favorable computational requirements of the present schemes are discussed in detail. The performance of the recently proposed spin-component-scaled and spin-opposite-scaled (SOS) range-separated (RS) and long-range corrected (LC) DH functionals is comprehensively assessed, while popular hybrid and global DH approaches are also discussed. For the benchmark calculations, up-to-date test sets are selected with high-level coupled-cluster references. Our results show that the ADC(2)-based SOS-RS-PBE-P86 approach is the most accurate and robust functional. This method consistently outperforms the excellent SOS-ADC(2) approach for VIPs, although the results are somewhat less satisfactory for VEAs. Among the genuine DH functionals, the SOS-ωPBEPP86 approach is also recommended for describing ionization processes, but its performance is even less reliable for electron-attached states. In addition, surprisingly good results are attained by the LC hybrid ωB97X-D functional, where the corresponding occupied (unoccupied) orbital energies are retrieved as VIPs (VEAs) within the present formalism.

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.

Improving Results by Improving Densities: Density-Corrected Density Functional Theory

Density functional theory (DFT) calculations have become widespread in both chemistry and materials, because they usually provide useful accuracy at much lower computational cost than wavefunction-based methods. All practical DFT calculations require an approximation to the unknown exchange-correlation energy, which is then used self-consistently in the Kohn-Sham scheme to produce an approximate energy from an approximate density. Density-corrected DFT is simply the study of the relative contributions to the total energy error. In the vast majority of DFT calculations, the error due to the approximate density is negligible. But with certain classes of functionals applied to certain classes of problems, the density error is sufficiently large as to contribute to the energy noticeably, and its removal leads to much better results. These problems include reaction barriers, torsional barriers involving π-conjugation, halogen bonds, radicals and anions, most stretched bonds, etc. In all such cases, use of a more accurate density significantly improves performance, and often the simple expedient of using the Hartree-Fock density is enough. This Perspective explains what DC-DFT is, where it is likely to improve results, and how DC-DFT can produce more accurate functionals. We also outline challenges and prospects for the field.

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.

Screened range-separated hybrid by balancing the compact and slowly varying density regimes: Satisfaction of local density linear response

Due to their quantitative accuracy and ability to solve several difficulties, screened range-separated hybrid exchange-correlation functionals are now a standard approach for ab initio simulation of condensed matter systems. However, the screened range-separated hybrid functionals proposed so far are biased either toward compact or slowly varying densities. In this paper, we propose a screened range-separated hybrid functional, named HSEint, which can well describe these density regimes, achieving good accuracy for both molecular and solid-state systems. The semilocal part of the proposed functional is based on the PBEint generalized gradient approximation [E. Fabiano et al., Phys. Rev. B 82, 113104 (2010)], constructed for hybrid interfaces. To improve the functional performance, we employ exact or nearly exact constraints in the construction of range-separated hybrid functional, such as recovering of the local density linear response and semiclassical atom linear response.

Investigation of the Exchange-Correlation Potentials of Functionals Based on the Adiabatic Connection Interpolation

We have studied the correlation potentials produced by various adiabatic connection models (ACMs) for several atoms and molecules. The results have been compared to accurate reference potentials (coupled cluster and quantum Monte Carlo results) as well as to state-of-the-art ab initio DFT approaches. We have found that all the ACMs yield correlation potentials that exhibit a correct behavior, quite resembling scaled second-order Görling-Levy (GL2) potentials and including most of the physically meaningful features of the accurate reference data. The behavior and contribution of the strong-interaction limit potentials have also been investigated and discussed.

Wavefunction-like Correlation Model for Use in Hybrid Density Functionals

We present Unsöld-W12 (UW12), an approximation to the correlation energy of molecules that is an explicit functional of the single-particle reduced-density matrix. The approximation resembles one part of modern explicitly correlated second-order Møller-Plesset (MP2) theory and is intended as an alternative to MP2 in double-hybrid exchange-correlation functionals. Orbital optimization with UW12 is straightforward, and the UW12 energy is evaluated without a double summation over unoccupied orbitals, leading to a faster basis-set convergence than is seen in double-hybrid functionals. We suggest a one-parameter hybrid exchange-correlation functional XCH-BLYP-UW12. XCH-BLYP-UW12 is similar to double-hybrid functionals, but contains UW12 correlation instead of MP2 correlation. We find that XCH-BLYP-UW12 is more accurate than the existing double-hybrid functional B2-PLYP for small-molecule main-group reaction barrier heights and has roughly the same accuracy as the existing hybrid functional B3LYP for atomization energies.

Combining Wave Function Methods with Density Functional Theory for Excited States

We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.

Generalized Optimized Effective Potential for Orbital Functionals and Self-Consistent Calculation of Random Phase Approximations

A new self-consistent procedure for calculating the total energy with an orbital-dependent density functional approximation (DFA), the generalized optimized effective potential (GOEP), is developed in the present work. The GOEP is a nonlocal Hermitian potential that delivers the sets of occupied and virtual orbitals and minimizes the total energy. The GOEP optimization leads to the same minimum as does the orbital optimization. The GOEP method is promising as an effective optimization approach for orbital-dependent functionals, as demonstrated for the self-consistent calculations of the random phase approximation (RPA) to the correlation functionals in the particle-hole (ph) and particle-particle (pp) channels. The results show that the accuracy in describing the weakly interacting van der Waals systems is significantly improved in the self-consistent calculations. In particular, the important single excitations contribution in non-self-consistent RPA calculations can be captured self-consistently through the GOEP optimization, leading to orbital renormalization, without using the single excitations in the energy functional.