A note on the accuracy of KS-DFT densities

Ionization energies of metallocenes: a coupled cluster study of cobaltocene

Open-shell transition metal chemistry presents challenges to contemporary electronic structure methods, based on either density functional or wavefunction theory. While CCSD(T) is the well-trusted gold standard for maingroup thermochemistry, the accuracy and robustness of the method is less clear for open-shell transition metal chemistry, requiring benchmarking of CCSD(T)-based protocols against either higher-level theory or experiment. Ionization energies (IEs) of metallocenes provide an interesting test case with metallocenes being common redox reagents as well as playing roles as redox mediators and cocatalysts in redox catalysis. Using highly accurate ZEKE-MATI experimental measurements of gas phase adiabatic (5.3275 ± 0.0006 eV) and vertical (5.4424 ± 0.0006 eV) ionization energies of cobaltocene, we systematically assessed the accuracy of the local coupled-cluster method DLPNO-CCSD(T) with respect to geometry, reference determinant, basis set size and extrapolation schemes, PNO cut-off and extrapolation, local triples approximation, relativistic effects and core-valence correlation. We show that PNO errors are controllable via the recently introduced PNO extrapolation schemes and that the expensive iterative triples (T₁) contribution can be made more manageable by calculating it as a smaller-basis/smaller PNO-cutoff correction. The reference determinant turns out to be a critical aspect in these calculations with the HF determinant resulting in large DLPNO-CCSD(T) errors, likely due to the qualitatively flawed molecular orbital spectrum. The BP86 functional on the other hand was found to provide reference orbitals giving small DLPNO-CCSD(T) errors, likely due to more realistic orbitals as suggested by the more consistent MO spectrum compared to HF. A protocol including complete basis set extrapolations with correlation-consistent basis sets, complete PNO space extrapolations, iterative triples- and core-valence correlation corrections was found to give errors of -0.07 eV and -0.03 eV for adiabatic- and vertical-IE of cobaltocene, respectively, giving close to chemical accuracy for both properties. A computationally efficient DLPNO-CCSD(T) protocol was devised and tested against adiabatic ionization energies of 6 different metallocenes (V, Cr, Mn, Fe, Co, Ni). For the other metallocenes, the iterative triples (T₁) and PNO extrapolation contributions turn out to be even more important. The results give errors close to the experimental uncertainty, similar to recent auxiliary-field quantum Monte Carlo results. The quality of the reference determinant orbitals is identified as the main source of uncertainty in CCSD(T) calculations of metallocenes.

Sensitivity of coupled cluster electronic properties on the reference determinant: Can Kohn-Sham orbitals be more beneficial than Hartree-Fock orbitals?

Coupled cluster calculations are traditionally performed over Hartree-Fock reference orbitals (HF-CC methodology). However, in the literature it has been repeatedly raised whether the use of a Kohn-Sham reference (KS-CC methodology) might result in improved performance relative to HF-CC. In the present study, we re-examine the relation of HF-CC and KS-CC methods by comparing the results of widely applied truncated CC calculations (CCSD, CCSD(T), CCSDT) to the limit of full configuration interaction (FCI), which serves as an undebatable reference point of accuracy. Based on a series of CC calculations on diatoms and transition metal complexes, we demonstrate that no systematic improvement of coupled cluster electronic energies, densities and chemical reaction energies is expected when changing from HF to a KS reference. Nevertheless, fortuitous error cancellations might occasionally result in illusory improvement compared to HF-CC. Altogether, the application of KS-CC is not advantageous over HF-CC, but it is also not unreasonable as the choice of reference has negligible influence on the results at sufficiently high CC levels. In addition, KS-CC can be a particularly useful alternative if difficulties are encountered in HF or HF-CC convergence. It is also notable that KS-CC results are found to be practically independent of the chosen density functional, which implies that almost any KS-CC method can be used in place of HF-CC.

Benchmarking isotropic hyperfine coupling constants using (QTP) DFT functionals and coupled cluster theory

Significant effort has been devoted to benchmarking isotropic hyperfine coupling constants for both wavefunction and density-based approaches in recent years, as accurate theoretical predictions aid the fitting of experimental model Hamiltonians. However, literature examining the predictive quality of a Density Functional Theory (DFT) functional abiding by the Bartlett IP condition is absent. In an attempt to rectify this, we report isotropic hyperfine coupling constant predictions of 24 commonly used DFT functionals on a total of 56 radicals, with the intent of exploring the successes and failures of the Quantum Theory Project (QTP) line of DFT functionals (i.e., CAM-QTP00, CAM-QTP01, CAM-QTP02, and QTP17) for this property. Included in this benchmark study are both small and large organic radicals as well as transition metal complexes, all of which have been studied to some extent in prior work. Subsequent coupled-cluster singles and doubles (CCSD) and CCSD withperturbative triples [CCSD(T)] calculations on small and large organic radicals show modest improvement as compared to prior work and offer an additional avenue for evaluation of DFT functional performance. We find that the QTP17 and CAM-QTP00 functionals consistently underperform, despite being parameterized to satisfy an IP eigenvalue condition primarily focused on inner shell electrons. On the other hand, the CAM-QTP01 functional is the most accurate functional in both organic radical datasets. Furthermore, both CAM-QTP01 and CAM-QTP02 are the most accurate functionals tested on the transition metal dataset. A significant portion of functionals were found to have comparable errors (within 5–15 MHz), but the hybrid class of DFT functionals maintains a consistently optimal balance between accuracy and precision across all datasets.

New venues in electron density analysis

We provide a comprehensive overview of the chemical information within the electron density: how to extract information, but also how to obtain and how to assess the quality of the...

Electronic Energy and Local Property Errors at QTAIM Critical Points while Climbing Perdew's Ladder of Density-Functional Approximations

We investigate the relationships between electron-density and electronic-energy errors produced by modern exchange-correlation density-functional approximations belonging to all of the rungs of Perdew's ladder. To this aim, a panel of relevant (semi)local properties evaluated at critical points of the electron-density field (as defined within the framework of Bader's atoms-in-molecules theory) are computed on a large selection of molecular systems involved in thermodynamic, kinetic, and noncovalent interaction chemical databases using density functionals developed in a nonempirical and minimally and highly parametrized fashion. The comparison of their density- and energy-based performance, also discussed in terms of density-driven errors, casts light on the strengths and weaknesses of the most recent and efficient density-functional approximations.

On the Use of Normalized Metrics for Density Sensitivity Analysis in DFT

In the past years, there has been a discussion about how the errors in density functional theory might be related to errors in the self-consistent densities obtained from different density functional approximations. This, in turn, brings up the discussion about the different ways in which we can measure such errors and develop metrics that assess the sensitivity of calculated energies to changes in the density. It is important to realize that there cannot be a unique metric in order to look at this density sensitivity, simultaneously needing size-extensive and size-intensive metrics. In this study, we report two metrics that are widely applicable to any density functional approximation. We also show how they can be used to classify different chemical systems of interest with respect to their sensitivity to small variations in the density.

Representing Exact Electron Densities by a Single Slater Determinant in Finite Basis Sets

The premise for Kohn-Sham density functional theory is that the exact electron density can be generated by a set of orbitals in a single Slater determinant. While this is ensured in a complete basis set, it has been shown that it cannot hold in small basis sets. The present work probes how accurately a reference electron density of the full-CI type can be reproduced by a set of orbitals in a single Slater determinant, as a function of the basis set used for the fitting electron density. The key finding is that the fitting error may be significant for basis sets of double- or triple-ζ quality. It is also shown that it is important that the fitting basis set includes the same basis functions as used for generating the reference electron density. The main limitation in a given basis set is the lack of higher order polarization functions. The error for practical purposes becomes insignificant for basis sets of quadruple-ζ or better quality, and this should be the choice when assessing the accuracy of exchange-correlation functionals by comparing electron densities to accurate reference results generated by wave function methods. The methodology in the present work can be used to transform an electron density from a multideterminant wave function into a set of orbitals in a single Slater determinant, and this may be useful for developing and testing new exchange-correlation functionals.

How accurate is density functional theory in predicting spin density? An insight from the prediction of hyperfine coupling constants

Electron paramagnetic resonance (EPR) spectroscopy has been proven to be an important technique for studying paramagnetic systems. Probably, the most accessible EPR parameter and the one that provides a significant amount of information about molecular structure and spin density is the hyperfine coupling constant (HFCC). Hence, accurate quantum-chemical modeling of HFCCs is frequently essential to the adequate interpretation of EPR spectra. It requires the precise spin density, which is the difference between the densities of α- and β-electrons, and thus, its quality is expected to reflect the quality of the total electron density. The question of which approximate exchange-correlation density functional yields sufficiently accurate HFCCs, and thus, the spin density remains open. To assess the performance of well-established density functionals for calculating HFCCs, we used a series of 26 small paramagnetic species and compared the obtained results to the CCSD reference values. The performance of DFT was also tested on EPR-studied o-semiquinone radical interacting with water molecules and Mg²⁺ cation. The HFCCs were additionally calculated by the DLPNO-CCSD method, and this wave function-based technique was found superior to all functionals we tested. Although some functionals were found, on average, to be fairly efficient, we found that the most accurate functional is system-dependent, and therefore, the DLPNO-CCSD method should be preferred for theoretical investigations of the HFCCs and spin density.

Adventures in DFT by a wavefunction theorist

The attraction density functional theory (DFT) has for electronic structure theory is that it is easier to do computationally than ab initio, correlated wavefunction methods, due to its effective one-particle structure. On the contrary, ab initio theorists insist on the ability to converge to the right answer in appropriate limits, but this requires a treatment of the reduced two-particle density matrix. DFT avoids that by appealing to an "existence" theorem (not a constructive one) that all its effects are subsummed into a DFT functional of the one-particle density. However, the existence of thousands of DFT functionals emphasizes that there is no satisfactory way to systematically improve the Kohn-Sham (KS) version as most changes in parameterization or formulation seldom lead to a new functional that is genuinely better than others. Some researchers in the DFT community try to address this issue by imposing conditions rigorously derived from exact DFT considerations, but to date, no one has shown how this route will ever lead to converged results even for the ground state, much less for all the other electronic states obtained from time-dependent DFT that are critically important for chemistry. On the contrary, coupled-cluster (CC) theory and its equation-of-motion extensions provide rigorous results for both that KS-DFT methods are attempting to emulate. How to use them and their exact formal properties to tie CC theory to an effective one-particle form is the target of this perspective. This route addresses the devil's triangle of KS-DFT problems: the one-particle spectrum, self-interaction, and the integer discontinuity.

Localizing electron density errors in density functional theory

The accuracy of different density functional approximations is assessed through the use of quantum chemical topology on molecular electron densities. In particular, three simple yet ever-important systems are studied: N₂, CO and ethane. Our results exemplify how real-space descriptors can help understand the sources of errors in density functional theory, avoiding unwanted error compensation present in simplified statistical metrics. Errors in "well-built" functionals are shown to be concentrated in chemically meaningful regions of space, and hence they are predictable. Conversely, strongly parametrized functionals show isotropic errors that cannot be traced back to chemically transferable units. Moreover, we will show that energetic corrections are mapped back into improvements in the density in chemically meaningful regions. These results point at the relevance of real-space perspectives when parametrizing or relating energy and density errors.

The maximum occupancy condition for the localized property-optimized orbitals

It is shown analytically that the Chemist's Localized Property-optimized Orbitals (CLPOs), which are the localized orbitals obtainable from the results of ab initio calculations by the open-source program JANPA (http://janpa.sourceforge.net/) according to the recently proposed optimal property partitioning condition, form the Lewis structure with nearly maximum possible total electron occupancy. The conditions required for this additional optimality to hold are discussed. In particular, when a single-determinant wavefunction is used to describe the molecular system without a noticeable electron delocalization, CLPOs derived from this wavefunction approximately optimize the same target quantity as the Natural Bond Orbitals (NBOs), establishing in this way the link between the two sets of localized orbitals. The performance of CLPO and NBO methods is compared by using a dataset containing 7101 small molecules, and the relevant methodological features of both methods are discussed.

Kinetic-energy-based error quantification in Kohn–Sham density functional theory

We present a basis-independent metric to assess the quality of the electron density obtained from Kohn–Sham (KS) density functional theory (DFT).

Predictive Simulations of Ionization Energies of Solvated Halide Ions with Relativistic Embedded Equation of Motion Coupled Cluster Theory

A subsystem approach for obtaining electron binding energies in the valence region is presented and applied to the case of halide ions (X^{-},X=F-At) in water. This approach is based on electronic structure calculations combining the relativistic equation-of-motion coupled cluster method for electron detachment and density functional theory via the frozen density embedding approach, using structures from classical molecular dynamics with polarizable force fields for discrete systems (in our study, droplets containing the anion and 50 water molecules). Our results indicate that one can accurately capture both the large solvent effect observed for the halides and the splitting of their ionization signals due to the increasingly large spin-orbit coupling of the p_{3/2}-p_{1/2} manifold across the series, at an affordable computational cost. Furthermore, owing to the quantum mechanical treatment of both solute and solvent electron binding energies of semiquantitative quality are also obtained for (bulk) water as by-products of the calculations for the halogens (in droplets).