Communication: Coupled cluster and many-body perturbation theory for fractional charges and spins

Optimal Tuning Perspective of Range-Separated Double Hybrid Functionals

We study the optimal tuning of the free parameters in range-separated double hybrid functionals, based on enforcing the exact conditions of piecewise linearity and spin constancy. We find that introducing the range separation in both the exchange and the correlation terms allows for the minimization of both fractional charge and fractional spin errors for singlet atoms. The optimal set of parameters is system specific, underlining the importance of the tuning procedure. We test the performance of the resulting optimally tuned functionals for the dissociation curves of diatomic molecules. We find that they recover the correct dissociation curve for the one-electron system, H2+, and improve the dissociation curves of many-electron molecules such as H2 and Li2, but they also yield a nonphysical maximum and only converge to the correct dissociation limit at very large distances.

Regularized second-order correlation methods for extended systems

Second-order Møller-Plesset perturbation theory (MP2) constitutes the simplest form of many-body wavefunction theory and often provides a good compromise between efficiency and accuracy. There are, however, well-known limitations to this approach. In particular, MP2 is known to fail or diverge for some prototypical condensed matter systems like the homogeneous electron gas (HEG) and to overestimate dispersion-driven interactions in strongly polarizable systems. In this paper, we explore how the issues of MP2 for metallic, polarizable, and strongly correlated periodic systems can be ameliorated through regularization. To this end, two regularized second-order methods (including a new, size-extensive Brillouin-Wigner approach) are applied to the HEG, the one-dimensional Hubbard model, and the graphene-water interaction. We find that regularization consistently leads to improvements over the MP2 baseline and that different regularizers are appropriate for the various systems.

Pure non-local machine-learned density functional theory for electron correlation

Density-functional theory (DFT) is a rigorous and (in principle) exact framework for the description of the ground state properties of atoms, molecules and solids based on their electron density. While computationally efficient density-functional approximations (DFAs) have become essential tools in computational chemistry, their (semi-)local treatment of electron correlation has a number of well-known pathologies, e.g. related to electron self-interaction. Here, we present a type of machine-learning (ML) based DFA (termed Kernel Density Functional Approximation, KDFA) that is pure, non-local and transferable, and can be efficiently trained with fully quantitative reference methods. The functionals retain the mean-field computational cost of common DFAs and are shown to be applicable to non-covalent, ionic and covalent interactions, as well as across different system sizes. We demonstrate their remarkable possibilities by computing the free energy surface for the protonated water dimer at hitherto unfeasible gold-standard coupled cluster quality on a single commodity workstation.

Semilocal density functionals, while computationally efficient, do not account for non-local correlation. Here, the authors propose a machine-learning approach to DFT that leads to non-local and transferable functionals applicable to non-covalent, ionic and covalent interactions across system of different sizes.

Capturing static and dynamic correlation with ΔNO-MP2 and ΔNO-CCSD

The ΔNO method for static correlation is combined with second-order Møller-Plesset perturbation theory (MP2) and coupled-cluster singles and doubles (CCSD) to account for dynamic correlation. The MP2 and CCSD expressions are adapted from finite-temperature CCSD, which includes orbital occupancies and vacancies, and expanded orbital summations. Correlation is partitioned with the aid of damping factors incorporated into the MP2 and CCSD residual equations. Potential energy curves for a selection of diatomics are in good agreement with extrapolated full configuration interaction results and on par with conventional multireference approaches.

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.

Non-empirical, low-cost recovery of exact conditions with model-Hamiltonian inspired expressions in jmDFT

Density functional theory (DFT) is widely applied to both molecules and materials, but well known energetic delocalization and static correlation errors in practical exchange-correlation approximations limit quantitative accuracy. Common methods that correct energetic delocalization errors, such as the Hubbard U correction in DFT+U or Hartree-Fock exchange in global hybrids, do so at the cost of worsening static correlation errors. We recently introduced an alternate approach [Bajaj et al., J. Chem. Phys. 147, 191101 (2017)] known as judiciously modified DFT (jmDFT), wherein the deviation from exact behavior of semilocal functionals over both fractional spin and charge, i.e., the so-called flat plane, was used to motivate functional forms of second order analytic corrections. In this work, we introduce fully nonempirical expressions for all four coefficients in a DFT+U+J-inspired form of jmDFT, where all coefficients are obtained only from energies and eigenvalues of the integer-electron systems. We show good agreement for U and J coefficients obtained nonempirically as compared with the results of numerical fitting in a jmDFT U+J/J' correction. Incorporating the fully nonempirical jmDFT correction reduces and even eliminates the fractional spin error at the same time as eliminating the energetic delocalization error. We show that this approach extends beyond s-electron systems to higher angular momentum cases including p- and d-electrons. Finally, we diagnose some shortcomings of the current jmDFT approach that limit its ability to improve upon DFT results for cases such as weakly bound anions due to poor underlying semilocal functional behavior.

Finite Temperature Coupled Cluster Theories for Extended Systems

At zero temperature, coupled cluster theories are widely used to predict total energies, ground state expectation values, and even excited states for molecules and extended systems. However, for systems with a small band gap, such as metals, the zero-temperature approximation does not necessarily hold. Thermal effects may even give rise to interesting chemistry on metal surfaces. Most approaches to temperature dependent electronic properties employ finite temperature perturbation theory in the Matsubara frequency formulation. Computations require a large number of Matsubara frequencies to yield sufficiently accurate results, especially at low temperatures. This work, and independently the work of White and Chan J. Chem. Theory Comput. 2018 , DOI: 10.1021/acs.jctc.8b00773 , proposes a coupled cluster implementation directly in the imaginary time domain on the compact interval [0, β], closely related to the thermal cluster cumulant approach of Sanyal et al. [ Chem. Phys. Lett. 1992 , 192 , 55 - 61 ] , Sanyal et al. [ Phys. Rev. E 1993 , 48 , 3373 - 3389 ], and Mandal et al. [ Int. J. Mod. Phys. B 2003 , 17 , 5367 - 5377 ]. Here, the arising imaginary time dependent coupled cluster amplitude integral equations are solved in the linearized direct ring doubles approximation, also referred to as Tamm-Dancoff approximation with second order (linearized) screened exchange. In this framework, the transition from finite to zero temperature is uniform and comes at no extra costs, allowing to go to temperatures as low as room temperature. In this approximation, correlation grand potentials are calculated over a wide range of temperatures for solid lithium, a metallic system, and for solid silicon, a semiconductor.

Shrinking Self-Interaction Errors with the Fermi-Löwdin Orbital Self-Interaction-Corrected Density Functional Approximation

The self-interaction error (SIE) is one of the major drawbacks of practical exchange-correlation functionals for Kohn-Sham density functional theory. Despite this, the use of methods that explicitly remove SIE from approximate density functionals is scarce in the literature due to their relatively high computational cost and lack of consistent improvement over standard modern functionals. In this article we assess the performance of a novel approach recently proposed by Pederson, Ruzsinszky, and Perdew [ J. Chem. Phys. 2014, 140, 121103] for performing self-interaction free calculations in density functional theory based on Fermi orbitals. To this end, we employ test sets consisting of reaction energies that are considered particularly sensitive to SIE. We found that the parameter-free Fermi-Löwdin orbital self-interaction correction method combined with the standard local spin density approximation (LSDA) and Perdew-Burke-Ernzerhof (PBE) functionals gives a much better estimate of reaction energies compared to their parent LSDA and PBE functionals for most of the reactions in these two sets. They also perform on par with the global PBE0 and range-separated LC-ωPBE hybrids, which partially eliminate the SIE by including Hartree-Fock exchange. This shows the potential of the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method for practical density functional calculations without SIE.