Local-hybrid functional based on the correlation length

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.

Paramagnetic Nuclear Magnetic Resonance Shifts for Triplet Systems and Beyond with Modern Relativistic Density Functional Methods

An efficient framework for the calculation of paramagnetic NMR (pNMR) shifts within exact two-component (X2C) theory and (current-dependent) density functional theory (DFT) up to the class of local hybrid functionals (LHFs) is presented. Generally, pNMR shifts for systems with more than one unpaired electron depend on the orbital shielding contribution and a temperature-dependent term. The latter includes zero-field splitting (ZFS), hyperfine coupling (HFC), and the g-tensor. For consistency, we calculate these three tensors at the same level of theory, i.e., using scalar-relativistic X2C augmented with spin-orbit perturbation theory. Results for pNMR chemical shifts of transition-metal complexes reveal that this X2C-DFT framework can yield good results for both the shifts and the individual tensor contributions of metallocenes and related systems, especially if the HFC constant is large. For small HFC constants, the relative error is often large, and sometimes the sign may be off. 4d and 5d complexes with more complicated structures demonstrate the limitations of a fully DFT-based approach. Additionally, a Co-based complex with a very large ZFS and pronounced multireference character is not well described. Here, a hybrid DFT-multireference framework is necessary for accurate results. Our results show that X2C is sufficient to describe relativistic effects and computationally cheaper than a fully relativistic approach. Thus, it allows use of large basis sets for converged HFCs. Overall, current-dependent meta-generalized gradient approximations and LHFs show some potential; however, the currently available functionals leave a lot to be desired, and the predictive power is limited.

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.

Zero-field splitting parameters within exact two-component theory and modern density functional theory using seminumerical integration

An efficient implementation of zero-field splitting parameters based on the work of Schmitt et al. [J. Chem. Phys. 134, 194113 (2011)] is presented. Seminumerical integration techniques are used for the two-electron spin-dipole contribution and the response equations of the spin-orbit perturbation. The original formulation is further generalized. First, it is extended to meta-generalized gradient approximations and local hybrid functionals. For these functional classes, the response of the paramagnetic current density is considered in the coupled-perturbed Kohn-Sham equations for the spin-orbit perturbation term. Second, the spin-orbit perturbation is formulated within relativistic exact two-component theory and the screened nuclear spin-orbit (SNSO) approximation. The accuracy of the implementation is demonstrated for transition-metal and diatomic main-group compounds. The efficiency is assessed for Mn and Mo complexes. Here, it is found that coarse integration grids for the seminumerical schemes lead to drastic speedups while introducing clearly negligible errors. In addition, the SNSO approximation substantially reduces the computational demands and leads to very similar results as the spin-orbit mean field Ansatz.

Reducing Exact Two-Component Theory for NMR Couplings to a One-Component Approach: Efficiency and Accuracy

The self-consistent and complex spin-orbit exact two-component (X2C) formalism for NMR spin-spin coupling constants [ J. Chem. Theory Comput. 17, 2021, 3874-3994] is reduced to a scalar one-component ansatz. This way, the first-order response term can be partitioned into the Fermi-contact (FC) and spin-dipole (SD) interactions as well as the paramagnetic spin-orbit (PSO) contribution. The FC+SD terms are real and symmetric, while the PSO term is purely imaginary and antisymmetric. The relativistic one-component approach is combined with a modern density functional treatment up to local hybrid functionals including the response of the current density. Computational demands are reduced by factors of 8-24 as shown for a large tin compound consisting of 137 atoms. Limitations of the current ansatz are critically assessed for Sn, Pb, Pd, and Pt compounds, i.e. the one-component treatment is not sufficient for tin compounds featuring a few heavy halogen atoms.

Paramagnetic NMR Shielding Tensors Based on Scalar Exact Two-Component and Spin-Orbit Perturbation Theory

The temperature-dependent Fermi-contact and pseudocontact terms are important contributions to the paramagnetic NMR shielding tensor. Herein, we augment the scalar-relativistic (local) exact two-component (X2C) framework with spin-orbit perturbation theory including the screened nuclear spin-orbit correction for the EPR hyperfine coupling and g tensor to compute these temperature-dependent terms. The accuracy of this perturbative ansatz is assessed with the self-consistent spin-orbit two-component and four-component treatments serving as reference. This shows that the Fermi-contact and pseudocontact interaction is sufficiently described for paramagnetic NMR shifts; however, larger deviations are found for the EPR spectra and the principle components of the EPR properties of heavy elements. The impact of the perturbative treatment is further compared to that of the density functional approximation and the basis set. Large-scale calculations are routinely possible with the multipole-accelerated resolution of the identity approximation and the seminumerical exchange approximation, as shown for [CeTi₆O₃(OiPr)₉(salicylate)₆].

A Local Hybrid Exchange Functional Approximation from First Principles

Local hybrid functionals are a more flexible class of density functional approximations allowing for a position-dependent admixture of exact exchange. This additional flexibility, however, comes with a more involved mathematical form and a more complicated design. A common denominator for previously constructed local hybrid funtionals is usage of thermochemical benchmark data to construct these functionals. Herein, we design a local hybrid functional without relying on benchmark data. Instead, we construct it in a more ab initio manner, following the principles of modern meta-generalized gradient approximations and considering theoretical constrains. To achieve this, we make use of the density matrix expansion and a local mixing function based on an approximate correlation length. The accuracy of the developed density functional approximation is assessed for thermochemistry, excitation energies, polarizabilities, magnetizabilities, NMR spin-spincoupling constants, NMR shieldings and shifts, as well as EPR g-tensors and hyperfine coupling constants. Here, the new exchange functional shows a robust performance and is especially well suited for atomization energies, barrier heights, excitation energies, NMR coupling constants, and EPR properties, whereas it looses some ground for the NMR shifts.Therefore, the designed functional is a major step forwards for functionals that have been designed from first principles.

Communication: Impact of the current density on paramagnetic NMR properties

Meta-generalized gradient approximations (meta-GGAs) and local hybrid functionals generally depend on the kinetic energy density τ. For magnetic properties, this necessitates generalizations to ensure gauge invariance. In most implementations, τ is generalized by incorporating the external magnetic field. However, this introduces artifacts in the response of the density matrix and does not satisfy the iso-orbital constraint. Here, we extend previous approaches based on the current density to paramagnetic NMR shieldings and EPR g-tensors. The impact is assessed for main-group compounds and transition-metal complexes considering 25 density functional approximations. It is shown that the current density leads to substantial improvements-especially for the popular Minnesota and SCAN functional families. Thus, we strongly recommend to use the current density generalized τ in paramagnetic NMR and EPR calculations with meta-GGAs.

Hybrid functionals with local range separation: Accurate atomization energies and reaction barrier heights

Range-separated hybrid approximations to the exchange-correlation density functional mix exact and semi-local exchange in a position-dependent manner. In their conventional form, the range separation is controlled by a constant parameter. Turning this constant into a density functional leads to a locally space-dependent range-separation function and thus a more powerful and flexible range-separation approach. In this work, we explore the self-consistent implementation of a local range-separated hybrid, taking into account a one-electron self-interaction correction and the behavior under uniform density scaling. We discuss different forms of the local range-separation function that depend on the electron density, its gradient, and the kinetic energy density. For test sets of atomization energies, reaction barrier heights, and total energies of atoms, we demonstrate that our best model is a clear improvement over common global range-separated hybrid functionals and can compete with density functionals that contain multiple empirical parameters. Promising results for equilibrium bond lengths, harmonic vibrational frequencies, and vertical ionization potentials further underline the potential and flexibility of our approach.

NMR Coupling Constants Based on the Bethe-Salpeter Equation in the GW Approximation

We present the first steps to extend the Green's function GW method and the Bethe-Salpeter equation (BSE) to molecular response properties such as nuclear magnetic resonance (NMR) indirect spin-spin coupling constants. We discuss both a nonrelativistic one-component and a quasi-relativistic two-component formalism. The latter describes scalar-relativistic and spin-orbit effects and allows us to study heavy-element systems with reasonable accuracy. Efficiency is maintained by the application of the resolution of the identity approximation throughout. The performance is demonstrated using conventional central processing units (CPUs) and modern graphics processing units (GPUs) for molecules involving several thousand basis functions. Our results show that a large amount of Hartree-Fock exchange is vital to provide a sufficient Kohn-Sham starting point to compute the GW quasi-particle energies. As the GW-BSE approach is generally less accurate for triplet excitations or related properties such as the Fermi-contact interaction, the admixture of the Kohn-Sham correlation kernel through the contracted BSE (cBSE) method improves the results for NMR coupling constants. This leads to remarkable results when combined with the eigenvalue-only self-consistent variant (evGW) and Becke's half and half functional (BH&HLYP) or the CAM-QTP family. The developed methodology is used to calculate the Karplus curve of tin molecules, illustrating its applicability to extended chemically relevant molecules. Here, the GW-cBSE method improves upon the chosen BH&HLYP Kohn-Sham starting points.

Systematic Evaluation of Modern Density Functional Methods for the Computation of NMR Shifts of 3d Transition-Metal Nuclei

A wide range of density functionals from all rungs of Jacob's ladder have been evaluated systematically for a set of experimental 3d transition-metal NMR shifts of 70 complexes encompassing 12 × ⁴⁹Ti, 10 × ⁵¹V, 10 × ⁵³Cr, 11 × ⁵⁵Mn, 9 × ⁵⁷Fe, 9 × ⁵⁹Co, and 9 × ⁶¹Ni shift values, as well as a diverse range of electronic structure characteristics. The overall 39 functionals evaluated include one LDA, eight GGAs, seven meta-GGAs (including their current-density-functional─CDFT─versions), nine global hybrids, four range-separated hybrids, eight local hybrids, and two double hybrids, and we also include Hartree-Fock and MP2 calculations. While recent evaluations of the same functionals for a very large coupled-cluster-based benchmark of main-group shieldings and shifts achieved in some cases aggregate percentage mean absolute errors clearly below 2%, the best results for the present 3d-nuclei set are in the range between 4 and 5%. Strikingly, the overall best-performing functionals are the recently implemented CDFT versions of two meta-GGAs, namely cM06-L (4.0%) and cVSXC (4.3%), followed by cLH14t-calPBE (4.9%), B3LYP (5.0%), and cLH07t-SVWN (5.1%), i.e., the previously best-performing global hybrid and two local hybrids. A number of further functionals achieve aggregate deviations in the range 5-6%. Range-separated hybrids offer no particular advantage over global hybrids. Due to the overall poor performance of Hartree-Fock theory for all systems except the titanium complexes, MP2 and double-hybrid functionals are unsuitable for these 3d-nucleus shifts and provide large errors. Global hybrid functionals with larger EXX admixtures, such as BHLYP or M06-2X, also perform poorly, and some other highly parametrized global hybrids also are unsuitable. For many functionals depending on local kinetic energy τ, their CDFT variants perform much better than their "non-CDFT" versions. This holds notably also for the above-mentioned M06-L and VSXC, while the effect is small for τ-dependent local hybrids and can even be somewhat detrimental to the agreement with experiment for a few other cases. The separation between well-performing and more poorly performing functionals is mainly determined by their results for the most critical nuclei ⁵⁵Mn, ⁵⁷Fe, and ⁵⁹Co. Here either moderate exact-exchange admixtures or CDFT versions of meta-GGAs are beneficial for the accuracy. The overall deviations of the better-performing global or local hybrids are then typically dominated by the ⁵³Cr shifts, where triplet instabilities appear to disfavor exact-exchange admixture. Further detailed analyses help to pinpoint specific nuclei and specific types of complexes that are challenges for a given functional.

Hyperfine Coupling Constants in Local Exact Two-Component Theory

We present a highly efficient implementation of the electron-nucleus hyperfine coupling matrix within the one-electron exact two-component (X2C) theory. The complete derivative of the X2C Hamiltonian is formed, that is, the derivatives of the unitary decoupling transformation are considered. This requires the solution of the response and Sylvester equations, consequently increasing the computational costs. Therefore, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The finite nucleus model is employed for both the scalar potential and the vector potential. Two-electron picture-change effects are modeled with the (modified) screened nuclear spin-orbit approach. Our implementation is fully integral direct and OpenMP-parallelized. An extensive benchmark study regarding the Hamiltonian, the basis set, and the density functional approximation is carried out for a set of 12-17 transition-metal compounds. The error introduced by DLU is negligible, and the DLU-X2C Hamiltonian accurately reproduces its four-component "fully" relativistic parent results. Functionals with a large amount of Hartree-Fock exchange such as CAM-QTP-02 and ωB97X-D are generally favorable. The pure density functional r²SCAN performs remarkably and even outperforms the common hybrid functionals TPSSh and CAM-B3LYP. Fully uncontracted basis sets or contracted quadruple-ζ bases are required for accurate results. The capability of our implementation is demonstrated for [Pt(C₆Cl₅)₄]^- with more than 4700 primitive basis functions and four rare-earth single-molecule magnets: [La(OAr*)₃]^-, [Lu(NR₂)₃]^-, [Lu(OAr*)₃]^-, and [TbPc₂]^-. Here, the results with the spin-orbit DLU-X2C Hamiltonian are in an excellent agreement with the experimental findings of all Pt, La, Lu, and Tb molecules.

Extended Benchmark Set of Main-Group Nuclear Shielding Constants and NMR Chemical Shifts and Its Use to Evaluate Modern DFT Methods

An extended theoretical benchmark set, NS372, for light main-group nuclear shieldings and NMR shifts has been constructed based on high-level GIAO-CCSD(T)/pcSseg-3//CCSD(T)/cc-pVQZ reference data. After removal of the large static-correlation cases O₃, F₃^-, and BH from the statistical evaluations for the ¹⁷O, ¹⁹F, and ¹¹B subsets, the benchmark comprises overall 372 shielding values in 117 molecules with a wide range of electronic-structure situations, containing 124 ¹H, 14 ¹¹B, 93 ¹³C, 43 ¹⁵N, 31 ¹⁷O, 47 ¹⁹F, 14 ³¹P, and 6 ³³S shielding constants. The CCSD(T)/pcSseg-3 data are shown to be close to the basis-set and method limit and thus provide an excellent benchmark to evaluate more approximate methods, such as density functional approaches. This dataset has been used to evaluate Hartree-Fock (HF) and MP2, and a wide range of exchange-correlation functionals from local density approximation (LDA) to generalized gradient approximations (GGAs) and meta-GGAs (focusing on their current-density functional implementations), as well as global hybrid, range-separated hybrid, local hybrid, and double-hybrid functionals. Starting with absolute shielding constants, the DSD-PBEP86 double hybrid is confirmed to provide the highest accuracy, with an aggregate relative mean absolute error (rel. MAE) of only 0.9%, followed by MP2 (1.1%). MP2 and double hybrids only show larger errors for a few systems with the largest static-correlation effects. The double-hybrid B2GP-PLYP, the two local hybrids cLH12ct-SsirPW92 and cLH12ct-SsifPW92, and the current-density functional meta-GGA cB97M-V follow closely behind (all 1.5%), as do some further functionals, cLH20t and cMN15-L (both 1.6%), as well as B2PLYP and KT3 (both 2.0%). Functionals on the lower rungs of the usual ladder offer the advantage of lower computational cost and access to larger molecules. Closer examination also reveals the best-performing methods for individual nuclei in the test set. Different ways of treating τ-dependent functionals are evaluated. When moving from absolute shielding constants to chemical shifts, some of the methods can benefit from systematic error compensation, and the overall error range somewhat narrows. Further methods now achieve the 2% threshold of relative MAEs, including functionals based on TPSS (TPSSh, cmPSTS).

Assessment of hybrid functionals for singlet and triplet excitations: Why do some local hybrid functionals perform so well for triplet excitation energies?

The performance of various hybrid density functionals is assessed for 105 singlet and 105 corresponding triplet vertical excitation energies from the QUEST database. The overall lowest mean absolute error is obtained with the local hybrid (LH) functional LH12ct-SsirPW92 with individual errors of 0.11 eV (0.11 eV) for singlet (triplet) n → π* excitations and 0.29 eV (0.17 eV) for π → π* excitations. This is slightly better than with the overall best performing global hybrid M06-2X [n → π*: 0.13 eV (0.17 eV), π → π*: 0.30 eV (0.20 eV)], while most other global and range-separated hybrids and some LHs suffer from the "triplet problem" of time-dependent density functional theory. This is exemplified by correlating the errors for singlet and triplet excitations on a state-by-state basis. The excellent performance of LHs based on a common local mixing function, i.e., an LMF constructed from the spin-summed rather than the spin-resolved semilocal quantities, is systematically investigated by the introduction of a spin-channel interpolation scheme that allows us to continuously modulate the fraction of opposite-spin terms used in the LMF. The correlation of triplet and singlet errors is systematically improved for the n → π* excitations when larger fractions of the opposite-spin-channel are used in the LMF, whereas this effect is limited for the π → π* excitations. This strongly supports a previously made hypothesis that attributes the excellent performance of LHs based on a common LMF to cross-spin-channel nondynamical correlation terms.

NMR Indirect Spin-Spin Coupling Constants in a Modern Quasi-Relativistic Density Functional Framework

A quasi-relativistic implementation of NMR indirect spin-spin coupling constants is presented. The exact two-component (X2C) Hamiltonian and its diagonal local approximation to the unitary decoupling transformation (DLU) are utilized together with the (modified) screened nuclear spin-orbit approach. In a restricted kinetic balance, the finite nucleus model is available for both the scalar and vector potentials. The implementation supports density functionals up to the fourth rung of Jacob's ladder, i.e., (range-separated) hybrid and local hybrid functionals based on a seminumerical ansatz. We assess the quality of our quasi-relativistic X2C approach by comparison with "fully" relativistic four-component results for small main-group molecules and alkynyl compounds. The mean absolute error introduced by the DLU scheme is less than 0.05 × 10¹⁹ T J^-2 of the reduced coupling constant for the small main-group molecules and 0.5 Hz for the alkynyl compounds. Thus, the error is significantly smaller than finite nucleus size effects for heavy elements. The basis set convergence and the impact of different density functional approximations are further studied. We propose a simple scheme to develop segmented-contracted relativistic all-electron basis sets for NMR spin-spin couplings. Our implementation allows us to perform calculations of extended molecules with reasonable computational effort, which is illustrated for the ¹J(¹¹⁹Sn, ³¹P) coupling constant of a low-valent tin phosphinidenide complex. The corresponding results are in good agreement with the experimental findings.

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.

Assessing the Accuracy of Local Hybrid Density Functional Approximations for Molecular Response Properties

A comprehensive overview of the performance of local hybrid functionals for molecular properties like excited states, ionization potentials within the GW framework, polarizabilities, magnetizabilities, NMR chemical shifts, and NMR spin-spin coupling constants is presented. We apply the generalization of the kinetic energy, τ, with the paramagnetic current density to all magnetic properties and the excitation energies from time-dependent density functional theory. This restores gauge invariance for these properties. Different ansätze for local mixing functions such as the iso-orbital indicator, the correlation length, the Görling-Levy second-order limit, and the spin polarization are compared. For the latter, we propose a modified version of the corresponding hyper-generalized gradient approximation functional of Perdew, Staroverov, Tao, and Scuseria (PSTS) [Phys. Rev. A 2008, 78, 052513] to allow for a numerically stable evaluation of the exchange-correlation kernel and hyperkernel. The PSTS functional leads to a very consistent improvement compared to the related TPSSh functional. It is further shown that the "best" choice of the local mixing function depends on the studied property and molecular class. While functionals based on the iso-orbital indicator lead to rather accurate excitation energies and ionization energies, the results are less impressive for NMR properties, for which a considerable dependence on the considered molecular test set and nuclei is observed. Johnson's local hybrid functional based on the correlation length yields remarkable results for NMR shifts of compounds featuring heavy elements and also for the excitation energies of organic compounds.

A comprehensive analysis of the history of DFT based on the bibliometric method RPYS

This bibliometric study aims at providing a comprehensive analysis of the history of density functional theory (DFT) from a perspective of chemistry by using reference publication year spectroscopy (RPYS). 114,138 publications with their 4,412,152 non-distinct cited references are analyzed. The RPYS analysis revealed three different groups of seminal papers which researchers in DFT have drawn from: (i) some long-known experimental studies from the 19th century about physical and chemical phenomena were referenced rather frequently in contemporary DFT publications. (ii) Fundamental quantum-chemical papers from the time period 1900–1950 which predate DFT form another group of seminal papers. (iii) Finally, various very frequently employed DFT approximations, basis sets, and other techniques (e.g., implicit descriptions of solvents) constitute another group of seminal papers. The earliest cited reference we found was published in 1806. The references to papers published in the 19th century mainly served the purpose of referring to long-known physical and chemical phenomena which were used to test if DFT approximations deliver correct results (e.g., Van der Waals interactions). The foundational papers of DFT by Hohenberg and Kohn as well as Kohn and Sham do not seem to be affected by obliteration by incorporation as they appear as pronounced peaks in our RPYS analysis. Since the 1990s, only very few pronounced peaks occur as most years were referenced nearly equally often. Exceptions are 1993 and 1996 due to seminal papers by Axel Becke, John P. Perdew and co-workers, and Georg Kresse and co-workers.

Separation of dynamic and nondynamic correlation

The account of electron correlation and its efficient separation into dynamic and nondynamic parts plays a key role in the development of computational methods. In this paper we suggest a physically-sound matrix formulation to split electron correlation into dynamic and nondynamic parts using the two-particle cumulant matrix and a measure of the deviation from idempotency of the first-order density matrix. These matrices are applied to a two-electron model, giving rise to a simplified electron correlation index that (i) depends only on natural orbitals and their occupancies, (ii) can be straightforwardly decomposed into orbital contributions and (iii) splits into dynamic and nondynamic correlation parts that (iv) admit a local version. These expressions are shown to account for dynamic and nondynamic correlation in a variety of systems containing different electron correlation regimes, thus providing the first separation of dynamic and nondynamic correlation using solely natural orbital occupancies.

Evolution of DFT studies in view of a scientometric perspective

Background

This bibliometric study aims to analyze the publications in which density functional theory (DFT) plays a major role. The bibliometric analysis is performed on the full publication volume of 114,138 publications as well as sub-sets defined in terms of six different types of compounds and nine different research topics. Also, a compound analysis is presented that shows how many compounds with specific elements are known to be calculated with DFT. This analysis is done for each element from hydrogen to nobelium.

Results

We find that hydrogen, carbon, nitrogen, and oxygen occur most often in compounds calculated with DFT in terms of absolute numbers, but a relative perspective shows that DFT calculations were performed rather often in comparison with experiments for rare gas elements, many actinides, some transition metals, and polonium.

Conclusions

The annual publication volume of DFT literature continues to grow steadily. The number of publications doubles approximately every 5–6 years while a doubling of publication volume every 11 years is observed for the CAplus database (14 years if patents are excluded). Calculations of the structure and energy of compounds dominate the DFT literature.

New approaches for the calibration of exchange-energy densities in local hybrid functionals

A new derivation of semi-local calibration functions is proposed, to ameliorate the gauge problem of exchange-energy densities in local hybrid functionals.