Full Implementation, Optimization, and Evaluation of a Range-Separated Local Hybrid Functional with Wide Accuracy for Ground and Excited States

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.

Combining Local Range Separation and Local Hybrids: A Step in the Quest for Obtaining Good Energies and Eigenvalues from One Functional

Some of the most successful exchange-correlation approximations in density functional theory are "hybrids", i.e., they rely on combining semilocal density functionals with exact nonlocal Fock exchange. In recent years, two classes of hybrid functionals have emerged as particularly promising: range-separated hybrids on the one hand, and local hybrids on the other hand. These functionals offer the hope to overcome a long-standing "observable dilemma", i.e., the fact that density functionals typically yield either a good description of binding energies, as obtained, e.g., in global and local hybrids, or physically interpretable eigenvalues, as obtained, e.g., in optimally tuned range-separated hybrids. Obtaining both of these characteristics from one and the same functional with the same set of parameters has been a long-standing challenge. We here discuss combining the concepts of local range separation and local hybrids as part of a constraint-guided quest for functionals that overcome the observable dilemma.

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.

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.

Visualizing and characterizing excited states from time-dependent density functional theory

Time-dependent density functional theory (TD-DFT) is the most widely-used electronic structure method for excited states, due to a favorable combination of low cost and semi-quantitative accuracy in many contexts, even if there are well recognized limitations. This Perspective describes various ways in which excited states from TD-DFT calculations can be visualized and analyzed, both qualitatively and quantitatively. This includes not just orbitals and densities but also well-defined statistical measures of electron-hole separation and of Frenkel-type exciton delocalization. Emphasis is placed on mathematical connections between methods that have often been discussed separately. Particular attention is paid to charge-transfer diagnostics, which provide indicators of when TD-DFT may not be trustworthy due to its categorical failure to describe long-range electron transfer. Measures of exciton size and charge separation that are directly connected to the underlying transition density are recommended over more ad hoc metrics for quantifying charge-transfer character.

Revisiting Gauge-Independent Kinetic Energy Densities in Meta-GGAs and Local Hybrid Calculations of Magnetizabilities

In a recent study [J. Chem. Theory Comput. 2021, 17, 1457-1468], some of us examined the accuracy of magnetizabilities calculated with density functionals representing the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA (mGGA), as well as global hybrid (GH) and range-separated (RS) hybrid functionals by assessment against accurate reference values obtained with coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)]. Our study was later extended to local hybrid (LH) functionals by Holzer et al. [J. Chem. Theory Comput. 2021, 17, 2928-2947]; in this work, we examine a larger selection of LH functionals, also including range-separated LH (RSLH) functionals and strong-correlation LH (scLH) functionals. Holzer et al. also studied the importance of the physically correct handling of the magnetic gauge dependence of the kinetic energy density (τ) in mGGA calculations by comparing the Maximoff-Scuseria formulation of τ used in our aforementioned study to the more physical current-density extension derived by Dobson. In this work, we also revisit this comparison with a larger selection of mGGA functionals. We find that the newly tested LH, RSLH, and scLH functionals outperform all of the functionals considered in the previous studies. The various LH functionals afford the seven lowest mean absolute errors while also showing remarkably small standard deviations and mean errors. Most strikingly, the best two functionals are scLHs that also perform remarkably well in cases with significant multiconfigurational character, such as the ozone molecule, which is traditionally excluded from statistical error evaluations due to its large errors with common density functionals.

Toward a correct treatment of core properties with local hybrid functionals

In local hybrid functionals (LHs), a local mixing function (LMF) determines the position-dependent exact-exchange admixture. We report new LHs that focus on an improvement of the LMF in the core region while retaining or partly improving upon the high accuracy in the valence region exhibited by the LH20t functional. The suggested new pt-LMFs are based on a Padé form and modify the previously used ratio between von Weizsäcker and Kohn-Sham local kinetic energies by different powers of the density to enable flexibly improved approximations to the correct high-density and iso-orbital limits relevant for the innermost core region. Using TDDFT calculations for a set of K-shell core excitations of second- and third-period systems including accurate state-of-the-art relativistic orbital corrections, the core part of the LMF is optimized, while the valence part is optimized as previously reported for test sets of atomization energies and reaction barriers (Haasler et al., J Chem Theory Comput 2020, 16, 5645). The LHs are completed by a calibration function that minimizes spurious nondynamical correlation effects caused by the gauge ambiguities of exchange-energy densities, as well as by B95c meta-GGA correlation. The resulting LH23pt functional relates to the previous LH20t functional but specifically improves upon the core region.

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

Extending recent developments on strong-correlation (sc) corrections to local hybrid functionals to the recent accurate ωLH22t range-separated local hybrid, a series of highly flexible strong-correlation-corrected range-separated local hybrids (scRSLHs) has been constructed and evaluated. This has required the position-dependent reduction of both short- and long-range exact-exchange admixtures in regions of space characterized by strong static correlations. Using damping procedures provides scRSLHs that retain largely the excellent performance of ωLH22t for weakly correlated situations and, in particular, for accurate quasiparticle energies of a wide variety of systems while reducing dramatically static-correlation errors, e.g., in stretched-bond situations. An additional correction to the local mixing function to reduce delocalization errors in abnormal open-shell situations provides further improvements in thermochemical and kinetic parameters, making scRSLH functionals such as ωLH23tdE or ωLH23tdP promising tools for complex molecular or condensed-phase systems, where low fractional-charge and fractional-spin errors are simultaneously important. The proposed rung 4 functionals thereby largely escape the usual zero-sum game between these two quantities and are expected to open new areas of accurate computations by Kohn-Sham DFT. At the same time, they require essentially no extra computational effort over the underlying ωLH22t functional, which means that their use is only moderately more demanding than that of global, local, or range-separated hybrid functionals.

Learning from the 4-(dimethylamino)benzonitrile twist: Two-parameter range-separated local hybrid functional with high accuracy for triplet and charge-transfer excitations

The recent ωLH22t range-separated local hybrid (RSLH) is shown to provide outstanding accuracy for the notorious benchmark problem of the two lowest excited-state potential energy curves for the amino group twist in 4-(dimethylamino)benzonitrile (DMABN). However, the design of ωLH22t as a general-purpose functional resulted in less convincing performance for triplet excitations, which is an important advantage of previous LHs. Furthermore, ωLH22t uses 8 empirical parameters to achieve broad accuracy. In this work, the RSLH ωLH23ct-sir is constructed with minimal empiricism by optimizing its local mixing function prefactor and range-separation parameter for only 8 excitation energies. ωLH23ct-sir maintains the excellent performance of ωLH22t for the DMABN twist and charge-transfer benchmarks but significantly improves the errors for triplet excitation energies (0.17 vs 0.24 eV). Additional test calculations for the AE6BH6 thermochemistry test set and large dipole moment and static polarizability test sets confirm that the focus on excitation energies in the optimization of ωLH23ct-sir has not caused any dramatic errors for ground-state properties. Although ωLH23ct-sir cannot replace ωLH22t as a general-purpose functional, it is preferable for problems requiring a universally good description of localized and charge-transfer excitations of both singlet and triplet multiplicity. Current limitations on the application of ωLH23ct-sir and other RSLHs to the study of singlet-triplet gaps of emitters for thermally activated delayed fluorescence are discussed. This work also includes the first systematic analysis of the influence of the local mixing function prefactor and the range-separation parameter in an RSLH on different types of excitations.

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Meta-generalized gradient approximations in time dependent generalized Kohn-Sham theory: Importance of the current density correction

We revisit the use of Meta-Generalized Gradient Approximations (mGGAs) in time-dependent density functional theory, reviewing conceptual questions and solving the generalized Kohn-Sham equations by real-time propagation. After discussing the technical aspects of using mGGAs in combination with pseudopotentials and comparing real-space and basis set results, we focus on investigating the importance of the current-density based gauge invariance correction. For the two modern mGGAs that we investigate in this work, TASK and r2SCAN, we observe that for some systems, the current density correction leads to negligible changes, but for others, it changes excitation energies by up to 40% and more than 0.8 eV. In the cases that we study, the agreement with the reference data is improved by the current density correction.

Correlation Consistent Basis Sets for Explicitly Correlated Theory: The Transition Metals

We present correlation consistent basis sets for explicitly correlated (F12) calculations, denoted VnZ(-PP)-F12-wis (n = D,T), for the d-block elements. The cc-pVDZ-F12-wis basis set is contracted to [8s7p5d2f] for the 3d-block, while its ECP counterpart for the 4d and 5d-blocks, cc-pVDZ-PP-F12-wis, is contracted to [6s6p5d2f]. The corresponding contracted sizes for cc-pVTZ(-PP)-F12-wis are [9s8p6d3f2g] for the 3d-block elements and [7s7p6d3f2g] for the 4d and 5d-block elements. Our VnZ(-PP)-F12-wis basis sets are evaluated on challenging test sets for metal-organic barrier heights (MOBH35) and group-11 metal clusters (CUAGAU-2). In F12 calculations, they are found to be about as close to the complete basis set limit as the combination of standard cc-pVnZ-F12 on main-group elements with the standard aug-cc-pV(n+1)Z(-PP) basis sets on the transition metal(s). While our basis sets are somewhat more compact than aug-cc-pV(n+1)Z(-PP), the CPU time benefit is negligible for catalytic complexes that contain only one or two transition metals among dozens of main-group elements; however, it is somewhat more significant for metal clusters.