Towards Efficient Orbital-Dependent Density Functionals for Weak and Strong Correlation

How accurate can Kohn-Sham density functional be for both main-group and transition metal reactions

Achieving chemical accuracy in describing reactions involving both main-group elements and transition metals poses a substantial challenge for density functional approximations (DFAs), primarily due to the significantly different behaviors for electrons moving in the s,p-orbitals or in the d,f-orbitals. MOR41, a representative dataset of transition metal chemistry, has highlighted the PWPB95-D3(BJ) functional, a B2PLYP-type doubly hybrid (bDH) approximation equipped with an empirical dispersion correction, as the leading functional thus far (Dohm et al., J Chem Theory Comput 2018;14: 2596-2608). However, this functional is not among the top bDH methods for main-group chemistry (Goerigk et al., Phys Chem Chem Phys. 2017;19: 32184). Conversely, bDH methods such as DSD-BLYP-D3, proficient in main-group chemistry, often falter for transition metal chemistry. Herein, taking advantage of the home-made Rust-based Electronic-Structure Toolkits, we examine a suite of XYG3-type doubly hybrid (xDH) methods. We confirm that the trade-off in descriptive accuracy between main-group and transition metal systems persists within the realm of perturbation theory (PT2)-based xDH methods. Notably, however, our study ushers in a pivotal advance with the recently proposed renormalized xDH method, R-xDH7-SCC15. This method not only distinguishes itself among the elite methods for main-group chemistry, but also achieves an unprecedented accuracy for the MOR41 dataset, outperforming all other reported DFAs. The efficacy of R-xDH7-SCC15 stems from the successful integration of a renormalized PT2 correlation model (rPT2) and a machine-learning strong-correlation correction (SCC15), marking a significant step forward in the realm of computational chemistry.

Toward Efficient and Unified Treatment of Static and Dynamic Correlations in Generalized Kohn-Sham Density Functional Theory

Accurate description of the static correlation poses a persistent challenge in electronic structure theory, particularly when it has to be concurrently considered with the dynamic correlation. We develop here a method in the generalized Kohn-Sham density functional theory (DFT) framework, named R-xDH7-SCC15, which achieves an unprecedented accuracy in capturing the static correlation, while maintaining a good description of the dynamic correlation on par with the state-of-the-art DFT and wave function theory methods, all grounded in the same single-reference black-box methodology. Central to R-xDH7-SCC15 is a general-purpose static correlation correction (SCC) model applied to the renormalized XYG3-type doubly hybrid method (R-xDH7). The SCC model development involves a hybrid machine learning strategy that integrates symbolic regression with nonlinear parameter optimization, aiming to achieve a balance between generalization capability, numerical accuracy, and interpretability. Extensive benchmark studies confirm the robustness and broad applicability of R-xDH7-SCC15 across a diverse array of main-group chemical scenarios. Notably, it displays exceptional aptitude in accurately characterizing intricate reaction kinetics and dynamic processes in regions distant from equilibrium, where the influence of static correlation is most profound. Its capability to consistently and efficiently predict the whole energy profiles, activation barriers, and reaction pathways within a user-friendly "black-box" framework represents an important advance in the field of electronic structure theory.

Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy

Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree-Fock reference is also), and orbital invariant. Our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

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.

Spin-Opposite-Scaled Range-Separated Exchange Double-Hybrid Models (SOS-RSX-DHs): Marriage Between DH and RSX/SOS-RSX Is Not Always a Happy Match

The range-separated version of double-hybrid density functional theory (DH-DFT) with a remarkable efficiency for both ground-state and excited-state characteristics has recently come into spotlight. In this work, based on theoretical arguments, several variants of spin-opposite-scaled range-separated exchange double-hybrid models (SOS-RSX-DHs) have been proposed and validated. More specifically, we first extend the RSX-DHs to design some other related models. Next, the SOS version of the resulting approximations is constructed and thoroughly evaluated using standard benchmark compilations of various properties. It is shown that although there are properties for which the RSX-DH and SOS-RSX-DH frameworks are rival, there are still some problems particularly prone to the self-interaction error issues where our proposed models seem to be beneficial. Furthermore, some of the presented models devoid of any additional corrections can also surpass the recently proposed approximations from different rungs of "Jacob's Ladder". Nonetheless, perusing the results of different methods and detailed comparisons with the predecessors discloses that all things may not necessarily be well with the RSX and SOS-RSX schemes, where the parent DHs as well as their SOS counterparts can still come into play.

Doubly Hybrid Functionals Close to Chemical Accuracy for Both Finite and Extended Systems: Implementation and Test of XYG3 and XYGJ-OS

While being widely used to understand the chemical reactions in heterogeneous catalysis or other multidisciplinary systems, a great challenge that semilocal and hybrid density functional approximations (DFAs) are facing is to deliver a uniformly accurate description for both finite and extended systems. Herein, we perform reliable and well-converged periodic calculations of two doubly hybrid approximations (DHAs), XYG3 and XYGJ-OS, and demonstrate that the good accuracy of DHAs achieved for molecules is transferable to the semiconductors and insulators. Such an accuracy is not only for energetic properties but also for the first- and second-order response properties, which is general for different kinds of chemical environments, including simple cubic bulks, perovskite-type transition metal oxides like TiO2, and heterogeneous systems like CO adsorption on the NaCl(100) surface. The present finding has strengthened the predictive power of DFT, which not only will inspire the future development of the top-rung DFAs but also will boost their applications in multidisciplinary studies with high accuracy and efficiency.

Pure non-local machine-learned density functional theory for electron correlation

Density-functional theory (DFT) is a rigorous and (in principle) exact framework for the description of the ground state properties of atoms, molecules and solids based on their electron density. While computationally efficient density-functional approximations (DFAs) have become essential tools in computational chemistry, their (semi-)local treatment of electron correlation has a number of well-known pathologies, e.g. related to electron self-interaction. Here, we present a type of machine-learning (ML) based DFA (termed Kernel Density Functional Approximation, KDFA) that is pure, non-local and transferable, and can be efficiently trained with fully quantitative reference methods. The functionals retain the mean-field computational cost of common DFAs and are shown to be applicable to non-covalent, ionic and covalent interactions, as well as across different system sizes. We demonstrate their remarkable possibilities by computing the free energy surface for the protonated water dimer at hitherto unfeasible gold-standard coupled cluster quality on a single commodity workstation.

Semilocal density functionals, while computationally efficient, do not account for non-local correlation. Here, the authors propose a machine-learning approach to DFT that leads to non-local and transferable functionals applicable to non-covalent, ionic and covalent interactions across system of different sizes.

The correlation factor approach: Combining density functional and wave function theory

Several of the limitations of approximate exchange-correlation functionals within Kohn-Sham density functional theory can be eliminated by extending the single-determinant reference system to a multi-determinant one. Here, we employ the correlation factor ansatz to combine multi-configurational, self-consistent field (MCSCF) with approximate density functionals. In the proposed correlation factor approach, the exchange-correlation hole ρ_XC(r, u), a function of the reference point r and the electron-electron separation u, is written as a product of the correlation factor f_C(r, u) and an exchange plus static-correlation hole ρ_XS(r, u), i.e., ρ_XC ^CFXS(r, u) = f_C(r, u)ρ_XS(r, u). ρ_XS(r, u) is constructed to reproduce the exchange-correlation energy of an MCSCF reference wave function. The correlation factor f_C(r, u) is designed to account for dynamic correlation effects that are absent in ρ_XS(r, u). The resulting approximation to the exchange-correlation energy, which we refer to as CFXStatic, is free of empirical parameters, and it combines the qualitatively correct description of the electronic structure obtainable with MCSCF with the advantages of approximate density functionals in accounting for dynamic correlation.

Simultaneous Attenuation of Both Self-Interaction Error and Nondynamic Correlation Error in Density Functional Theory: A Spin-Pair Distinctive Adiabatic-Connection Approximation

We present a spin-pair distinctive algorithm in the context of adiabatic-connection fluctuation-dissipation (ACFD) theorem, which enables to quantify the self-interaction error (SIE) and the nondynamic/strong correlation error (NCE) in the direct random-phase approximation (dRPA). Using this knowledge, we propose a spin-component scaled dRPA (scsRPA) correlation model with simultaneous attenuation of both the SIE and the NCE. Along with the exact exchange, scsRPA is shown to present a comprehensive improvement over dRPA, as well as the well-established PBE and PBE0 functionals, for bonding energies of pronounced multireference characters and transition-metal complexes of strongly correlated systems, while consistently providing an accurate description for reaction energies, reaction barriers, and noncovalent bond interactions of weakly correlated systems.

Range-Separated Double-Hybrid Functional from Nonempirical Constraints

On the basis of our previous developments in the field of nonempirical double hybrids, we present here a new exchange-correlation functional based on a range-separated model for the exchange part and integrating a nonlocal perturbative correction to the electron correlation contribution. Named RSX-QIDH, the functional is free from any kind of empirical parametrization. Its range-separation parameter is set to recover the total energy of the hydrogen atom, thus eliminating the self-interaction error for this one-electron system. Subsequent tests on some relevant benchmark data sets confirm that the self-interaction error is particularly low for RSX-QIDH. This new functional provides also correct dissociation profiles for charged rare-gas dimers and very accurate ionization potentials directly from Kohn-Sham orbital energies. Above all, these good results are not obtained at the expense of other properties. Indeed, further tests on standard benchmarks show that RSX-QIDH is competitive with the more empirical ωB97X-2 double hybrid and outperforms the parent LC-PBE long-range corrected hybrid, thus underlining the important role of the nonlocal perturbative correlation.

Doubly hybrid density functionals that correctly describe both density and energy for atoms

Recently, it was argued [Medvedev MG, et al. (2017) Science 355:49-52] that the development of density functional approximations (DFAs) is "straying from the path toward the exact functional." The exact functional should yield both exact energy and density for a system of interest; nevertheless, they found that many heavily fitted functionals for molecular energies actually lead to poor electron densities of atoms. They also observed a trend that, for the nonempirical and few-parameter functionals, densities can be improved as one climbs up the first four rungs of the Jacob's ladder of DFAs. The XYG3 type of doubly hybrid functionals (xDHs) represents a less-empirical and fewer-parameter functional on the top fifth rung, in which both the Hartree-Fock-like exchange and the second-order perturbative (MP2-like) correlation are hybridized with the low rung functionals. Here, we show that xDHs can well describe both density and energy for the same atomic set of Medvedev et al., showing that the latter trend can well be extended to the top fifth rung.

Simple Fully Nonlocal Density Functionals for Electronic Repulsion Energy

From a simplified version of the mathematical structure of the strong coupling limit of the exact exchange-correlation functional, we construct an approximation for the electronic repulsion energy at physical coupling strength, which is fully nonlocal. This functional is self-interaction free and yields energy densities within the definition of the electrostatic potential of the exchange-correlation hole that are locally accurate and have the correct asymptotic behavior. The model is able to capture strong correlation effects that arise from chemical bond dissociation, without relying on error cancellation. These features, which are usually missed by standard density functional theory (DFT) functionals, are captured by the highly nonlocal structure, which goes beyond the "Jacob's ladder" framework for functional construction, by using integrals of the density as the key ingredient. Possible routes for obtaining the full exchange-correlation functional by recovering the missing kinetic component of the correlation energy are also implemented and discussed.