Range-Separated Local Hybrid Functionals with Small Fractional-Charge and Fractional-Spin Errors: Escaping the Zero-Sum Game of DFT Functionals

Chemically Accurate Singlet-Triplet Gaps of Arylcarbenes from Local Hybrid Density Functionals

Singlet-triplet (ST) gaps are key descriptors of carbenes, because their properties and reactivity are strongly spin-dependent. However, the theoretical prediction of ST gaps is challenging and generally thought to require elaborate correlated wave function methods or double-hybrid density functionals. By evaluating two recent test sets of arylcarbenes (AC12 and AC18), we show that local hybrid functionals based on the "common t" local mixing function (LMF) model achieve mean absolute errors below 1 kcal/mol at a computational cost only slightly higher than that of global hybrid functionals. An analysis of correlation contributions to the ST gaps suggests that the accuracy of the common t-LMF model is mainly due to an improved description of nondynamical correlation which, unlike exchange, is not additive in each spin-channel. Although spin-nonadditivity can be achieved using the local spin polarization alone, using the "common", i.e., spin-unresolved, iso-orbital indicator t for constructing the LMF is found to be critical for consistent accuracy in ST gaps of arylcarbenes. The results support the view of LHs as vehicles to improve the description of nondynamical correlation rather than sophisticated exchange mixing approaches.

Toward the Next Generation of Density Functionals: Escaping the Zero-Sum Game by Using the Exact-Exchange Energy Density

ConspectusKohn-Sham density functional theory (KS DFT) is arguably the most widely applied electronic-structure method with tens of thousands of publications each year in a wide variety of fields. Its importance and usefulness can thus hardly be overstated. The central quantity that determines the accuracy of KS DFT calculations is the exchange-correlation functional. Its exact form is unknown, or better "unknowable", and therefore the derivation of ever more accurate yet efficiently applicable approximate functionals is the "holy grail" in the field. In this context, the simultaneous minimization of so-called delocalization errors and static correlation errors is the greatest challenge that needs to be overcome as we move toward more accurate yet computationally efficient methods. In many cases, an improvement on one of these two aspects (also often termed fractional-charge and fractional-spin errors, respectively) generates a deterioration in the other one. Here we report on recent notable progress in escaping this so-called "zero-sum-game" by constructing new functionals based on the exact-exchange energy density. In particular, local hybrid and range-separated local hybrid functionals are discussed that incorporate additional terms that deal with static correlation as well as with delocalization errors. Taking hints from other coordinate-space models of nondynamical and strong electron correlations (the B13 and KP16/B13 models), position-dependent functions that cover these aspects in real space have been devised and incorporated into the local-mixing functions determining the position-dependence of exact-exchange admixture of local hybrids as well as into the treatment of range separation in range-separated local hybrids. While initial functionals followed closely the B13 and KP16/B13 frameworks, meanwhile simpler real-space functions based on ratios of semilocal and exact-exchange energy densities have been found, providing a basis for relatively simple and numerically convenient functionals. Notably, the correction terms can either increase or decrease exact-exchange admixture locally in real space (and in interelectronic-distance space), leading even to regions with negative admixture in cases of particularly strong static correlations. Efficient implementations into a fast computer code (Turbomole) using seminumerical integration techniques make such local hybrid and range-separated local hybrid functionals promising new tools for complicated composite systems in many research areas, where simultaneously small delocalization errors and static correlation errors are crucial. First real-world application examples of the new functionals are provided, including stretched bonds, symmetry-breaking and hyperfine coupling in open-shell transition-metal complexes, as well as a reduction of static correlation errors in the computation of nuclear shieldings and magnetizabilities. The newest versions of range-separated local hybrids (e.g., ωLH23tdE) retain the excellent frontier-orbital energies and correct asymptotic exchange-correlation potential of the underlying ωLH22t functional while improving substantially on strong-correlation cases. The form of these functionals can be further linked to the performance of the recent impactful deep-neural-network "black-box" functional DM21, which itself may be viewed as a range-separated local hybrid.

Multiconfigurational Correlation at DFT + U Cost: On-Site Electron-Electron Interactions Yield a Block-Localized Configuration Interaction Hamiltonian

This work presents a first-principles wavefunction-in-DFT approach based on the Hubbard density functional theory (DFT) + U method. This approach begins with the standard DFT reference system of noninteracting electrons and introduces an electron-electron interaction projected onto DFT+U-type atomic states. The reference system's configuration interaction Hamiltonian is block-localized to these states and can be expressed in terms of state occupation numbers, state self-energies (which correspond to unscreened Hubbard U values), and the promotion energies of doubly excited Slater determinants. Simple approximations for the promotion energies provide multiconfigurational correlation energies without requiring explicit orbital localization/transform. Numerical results for fractionally occupied chromium atom, bonded chromium dimer, dissociating covalent bonds, and large active spaces show that the approach provides beyond-zero-sum accuracy at computational cost comparable to standard DFT+U.

Implementation and First Evaluation of Strong-Correlation-Corrected Local Hybrid Functionals for the Calculation of NMR Shieldings and Shifts

Local hybrid functionals containing strong-correlation factors (scLHs) and range-separated local hybrids (RSLHs) have been integrated into an efficient coupled-perturbed Kohn-Sham implementation for the calculation of nuclear shielding constants. Several scLHs and the ωLH22t RSLH have then been evaluated for the first time for the extended NS372 benchmark set of main-group shieldings and shifts and the TM70 benchmark of 3d transition-metal shifts. The effects of the strong-correlation corrections have been analyzed with respect to the spatial distribution of the sc-factors, which locally diminish exact-exchange admixture at certain regions in a molecule. The scLH22t, scLH23t-mBR, and scLH23t-mBR-P functionals, which contain a "damped" strong-correlation factor to retain the excellent performance of the underlying LH20t functional for weakly correlated situations, tend to make smaller corrections to shieldings and shifts than the "undamped" scLH22ta functional. While the latter functional can also deteriorate agreement with the reference data in certain weakly correlated cases, it provides overall better performance, in particular for systems where static correlation is appreciable. This pertains only to a minority of systems in the NS372 main-group test set but to many more systems in the TM70 transition-metal test set, in particular for high-oxidation-state complexes, e.g., Cr(+VI) complexes and other systems with stretched bonds. Another undamped scLH, the simpler LDA-based scLH21ct-SVWN-m, also tends to provide significant improvements in many cases. The differences between the functionals and species can be rationalized on the basis of one-dimensional plots of the strong-correlation factors, augmented by isosurface plots of the fractional orbital density (FOD). Position-dependent exact-exchange admixture is thus shown to provide substantial flexibility in treating response properties like NMR shifts for both weakly and strongly correlated systems.

Spin-Symmetry Breaking and Hyperfine Couplings in Transition-Metal Complexes Revisited Using Density Functionals Based on the Exact-Exchange Energy Density

A small set of mononuclear manganese complexes evaluated previously for their Mn hyperfine couplings (HFCs) has been analyzed using density functionals based on the exact-exchange energy density─in particular, the spin symmetry breaking (SSB) found previously when using hybrid functionals. Employing various strong-correlation corrected local hybrids (scLHs) and strong-correlation corrected range-separated local hybrids (scRSLHs) with or without additional corrections to their local mixing functions (LMFs) to mitigate delocalization errors (DE), the SSB and the associated dipolar HFCs of [Mn(CN)4]2-, MnO3, [Mn(CN)4N]-, and [Mn(CN)5NO]2- (the latter with cluster embedding) have been examined. Both strong-correlation (sc)-correction and DE-correction terms help to diminish SSB and correct the dipolar HFCs. The DE corrections are more effective, and the effects of the sc corrections depend on their damping factors. Interestingly, the DE-corrections reduce valence-shell spin polarization (VSSP) and thus SSB by locally enhancing exact-exchange (EXX) admixture near the metal center and thereby diminishing spin-density delocalization onto the ligand atoms. In contrast, sc corrections diminish EXX admixture locally, mostly on specific ligand atoms. This then reduces VSSP and SSB as well. The performance of scLHs and scRSLHs for the isotropic Mn HFCs has also been analyzed, with particular attention to core-shell spin-polarization contributions. Further sc-corrected functionals, such as the KP16/B13 construction and the DM21 deep-neural-network functional, have been examined.