Practical auxiliary basis implementation of Rung 3.5 functionals

Unification of Perdew-Zunger self-interaction correction, DFT+U, and Rung 3.5 density functionals

This Communication presents a unified derivation of three different approximations used in density functional theory (DFT): the Perdew-Zunger self-interaction correction (PZSIC), the Hubbard correction DFT+U, and the Rung 3.5 density functionals. All three approximations can be derived by introducing electron self-interaction into the Kohn-Sham (KS) reference system of noninteracting electrons. The derivation uses the Adiabatic Projection formalism: one projects the electron-electron interaction operator onto certain states, introduces the projected operator into the reference system, and defines a density functional for the remainder. Projecting onto individual localized KS orbitals recovers our previous derivation of the PZSIC [B. G. Janesko, J. Phys. Chem. Lett. 13, 5698-5702 (2022)]. Projecting onto localized atom-centered orbitals recovers a variant of DFT+U. Projecting onto localized states at each point in space recovers Rung 3.5 approaches. New results include an "atomic state PZSIC" that does not require localizing the KS orbitals, a demonstration that typical Hubbard U parameters reproduce a scaled-down PZSIC, and a Rung 3.5 variant of DFT+U that does not require choosing atom-dependent states.

Replacing hybrid density functional theory: motivation and recent advances

Density functional theory (DFT) is the most widely-used electronic structure approximation across chemistry, physics, and materials science. Every year, thousands of papers report hybrid DFT simulations of chemical structures, mechanisms, and spectra. Unfortunately, hybrid DFT's accuracy is ultimately limited by tradeoffs between over-delocalization and under-binding. This review summarizes these tradeoffs, and introduces six modern attempts to go beyond them while maintaining hybrid DFT's relatively low computational cost: DFT+U, self-interaction corrections, localized orbital scaling corrections, local hybrid functionals, real-space nondynamical correlation, and our rung-3.5 approach. The review concludes with practical suggestions for DFT users to identify and mitigate these tradeoffs' impact on their simulations.

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.

A General Route to Include Pauli Repulsion and Quantum Dispersion Effects in QM/MM Approaches

A methodology to account for nonelectrostatic interactions in Quantum Mechanical (QM)/Molecular Mechanics (MM) approaches is developed. Formulations for Pauli repulsion and dispersion energy, explicitly depending on the QM density, are derived. Such expressions are based on the definition of an auxiliary density on the MM portion and the Tkatchenko-Scheffler (TS) approach, respectively. The developed method is general enough to be applied to any QM/MM method and partition, provided an accurate tuning of a small number of parameters is obtained. The coupling of the method with both nonpolarizable and fully polarizable QM/fluctuating charge (FQ) approaches is reported and applied. A suitable parametrization for the aqueous solution, so that its most representative features are well reproduced, is outlined. Then, the obtained parametrization and method are applied to calculate the nonelectrostatic (repulsion and dispersion) interaction energy of nicotine in aqueous solution.

Comparative Study of Nonhybrid Density Functional Approximations for the Prediction of 3d Transition Metal Thermochemistry

The utility of several nonhybrid density functional approximations (DFAs) is considered for the prediction of gas phase enthalpies of formation for a large set of 3d transition metal-containing molecules. Nonhybrid DFAs can model thermochemical values for 3d transition metal-containing molecules with accuracy comparable to that of hybrid functionals. The GAM-generalized gradient approximation (GGA); the TPSS, M06-L, and MN15-L meta-GGAs; and the Rung 3.5 PBE+ΠLDA(s) DFAs all give root-mean-square deviations below that of the widely used B3LYP hybrid. Modern nonhybrid DFAs continue to show utility for transition metal thermochemistry.

Practical Density Functionals beyond the Overdelocalization-Underbinding Zero-Sum Game

Density functional theory (DFT) uses a density functional approximation (DFA) to add electron correlation to mean-field electronic structure calculations. Standard strategies (generalized gradient approximations GGAs, meta-GGAs, hybrids, etc.) for building DFAs, no matter whether based on exact constraints or empirical parametrization, all face a zero-sum game between overdelocalization (fractional charge error, FC) and underestimation of covalent bonding (fractional spin error, FS). This work presents an alternative strategy. Practical "Rung 3.5" ingredients are used to implement insights from hyper-GGA DFAs that reduce both FS and FC errors. Prototypes of this strategy qualitatively improve FS and FC error over 40 years of standard DFAs while maintaining low cost and practical evaluation of properties. Numerical results ranging from transition metal thermochemistry to absorbance peaks and excited-state geometry optimizations highlight this strategy's promise and indicate areas requiring further development.

Testing Exact Upper Bounds to Exact Exchange

The exact exchange energy and its energy density are useful but computationally expensive ingredients in density functional approximations for Kohn-Sham density functional theory. We present detailed tests of some exact nonempirical upper bounds to exact exchange. These "Rung 3.5" upper bounds contract the Kohn-Sham one-particle density matrix with model density matrices used to construct semilocal model exchange holes and invoke the Cauchy-Schwarz inequality. The contraction automatically eliminates the computationally expensive long-range component of the exact exchange hole. Numerical tests show that the exchange upper bounds underestimate total exchange energies while predicting other properties with accuracy approaching standard hybrid approximations.

Reducing density-driven error without exact exchange

PBE calculations, performed non-self-consistently on densities evaluated with Rung 3.5 density functionals, give improved performance for hydrogen transfer reaction barriers.

Quantifying solvated electrons' delocalization

The electron delocalization range EDR(r;u_av) (left) captures the spin density (right) of an electron delocalized over u_av = 5.77 Å on the surface of an (H₂O)₂₀⁻ cluster.