Abstract Local correlation methods rely on the assumption that electron correlation is nearsighted. In this work, we develop a method to alleviate this assumption. This new method is demonstrated by calculating the random phase approximation (RPA) correlation energies in several one‐dimensional model systems. In this new method, the first step is to approximately decompose the RPA correlation energy to the nearsighted and farsighted components based on the wavelength decomposition of electron correlation developed by Langreth and Perdew. The short‐wavelength (SW) component of the RPA correlation energy is then considered to be nearsighted, and the long‐wavelength (LW) component of the RPA correlation energy is considered to be farsighted. The SW RPA correlation energy is calculated using a recently developed local correlation method: the embedded cluster density approximation (ECDA). The LW RPA correlation energy is calculated globally based on the system's Kohn‐Sham orbitals. This new method is termed λ‐ECDA, where λ indicates the wavelength decomposition. The performance of λ‐ECDA is examined on a one‐dimensional model system: a H24 chain, in which the RPA correlation energy is highly nonlocal. In this model system, a softened Coulomb interaction is used to describe the electron‐electron and electron‐ion interactions, and slightly stronger nuclear charges (1.2e) are assigned to the pseudo‐H atoms. Bond stretching energies, RPA correlation potentials, and Kohn‐Sham eigenvalues predicted by λ‐ECDA are in good agreement with the benchmarks when the clusters are made reasonably large. We find that the LW RPA correlation energy is critical for obtaining accurate prediction of the RPA correlation potential, even though the LW RPA correlation energy contributes to only a few percent of the total RPA correlation energy.
Abstract Random phase approximation (RPA) electron correlation methods have gained in popularity in the recent years. A number of RPA correlation energy variants emerged in the Kohn–Sham DFT framework and in the theory of strongly orthogonal geminals. Foundations of most RPA approaches trace back to an exact expression for two‐electron reduced density matrix written in terms of one‐electron density matrix and dynamic one‐electron response functions, originally presented in the seminal paper of McLachlan and Ball (Rev. Mod. Phys. 1964, 36, 844). The aim of this article is to give a pedagogical introduction to possible approaches for describing electron correlation energy based on the McLachlan and Ball formula. The focus of the presentation is to formulate electron‐interaction energy expressions as functions of reduced density matrices. On one hand, it provides a common umbrella for RPA approximations proposed for uncorrelated (Hartree–Fock or Kohn–Sham) as well as partially correlated (strongly orthogonal geminals) references. On the other hand, such presentation may stimulate new developments in density matrix functional theory.
The exact static and time-dependent Kohn-Sham (KS) exchange-correlation potential is extremely challenging to approximate as it is a local multiplicative potential that depends on the electron density everywhere in the system. The KS approach can be generalized by allowing part of the potential to be spatially nonlocal. We take this nonlocal part to be that of unrestricted Hartree-Fock theory. The additional local correlation potential in principle ensures that the single-particle density exactly equals the many-body density. In our case, the local correlation potential is predominantly nearsighted in its dependence on the density and hence an (adiabatic) local-density approximation to this potential yields accurate ground-state properties and real-time densities for one-dimensional test systems.
We use the exact strong-interaction limit of the Hohenberg-Kohn energy density functional to approximate the exchange-correlation energy of the restricted Kohn-Sham scheme. Our approximation corresponds to a highly nonlocal density functional whose functional derivative can be easily constructed, thus transforming exactly, in a physically transparent way, an important part of the electron-electron interaction into an effective local one-body potential. We test our approach on quasi-one-dimensional systems, showing that it captures essential features of strong correlation that restricted Kohn-Sham calculations using the currently available approximations cannot describe.
We present an embedding approach to treat local electron correlation effects in periodic environments. In a single consistent framework, our plane wave based scheme embeds a local high-level correlation calculation [here, Coupled Cluster (CC) theory], employing localized orbitals, into a low-level correlation calculation [here, the direct Random Phase Approximation (RPA)]. This choice allows for an accurate and efficient treatment of long-range dispersion effects. Accelerated convergence with respect to the local fragment size can be observed if the low-level and high-level long-range dispersions are quantitatively similar, as is the case for CC in RPA. To demonstrate the capabilities of the introduced embedding approach, we calculate adsorption energies of molecules on a surface and in a chabazite crystal cage, as well as the formation energy of a lattice impurity in a solid at the level of highly accurate many-electron perturbation theories. The absorption energy of a methane molecule in a zeolite chabazite is converged with an error well below 20 meV at the CC level. As our largest periodic benchmark system, we apply our scheme to the adsorption of a water molecule on titania in a supercell containing more than 1000 electrons.
Perhaps the most important approximations to the electronic structure problem in quantum chemistry are those based on coupled cluster and density functional theories. Coupled cluster theory has been called the ``gold standard'' of quantum chemistry due to the high accuracy that it achieves for weakly correlated systems. Kohn--Sham density functionals based on semilocal approximations are, without a doubt, the most widely used methods in chemistry and material science because of their high accuracy/cost ratio. The root of the success of coupled cluster and density functionals is their ability to efficiently describe the dynamic part of the electron correlation. However, both traditional coupled cluster and density functional approximations may fail catastrophically when substantial static correlation is present. This severely limits the applicability of these methods to a plethora of important chemical and physical problems such as, e.g., the description of bond breaking, transition states, transition metal-, lanthanide- and actinide-containing compounds, and superconductivity. In an attempt to tackle this problem, nonstandard (single-reference) coupled cluster-based techniques that aim to describe static correlation have been recently developed: pair coupled cluster doubles (pCCD) and singlet-paired coupled cluster doubles (CCD0). The ability to describe static correlation in pCCD and CCD0 comes, however, at the expense of important amounts of dynamic correlation so that the high accuracy of standard coupled cluster becomes unattainable. Thus, the reliable and efficient description of static and dynamic correlation in a simultaneous manner remains an open problem for quantum chemistry and many-body theory in general. In this thesis, different ways to combine pCCD and CCD0 with density functionals in order to describe static and dynamic correlation simultaneously (and efficiently) are explored. The combination of wavefunction and density functional methods has a long history in quantum chemistry (practical implementations have appeared in the literature since the 1970s). However, this kind of techniques have not achieved widespread use due to problems such as double counting of correlation and the symmetry dilemma---the fact that wavefunction methods respect the symmetries of Hamiltonian, while modern functionals are designed to work with broken symmetry densities. Here, particular mathematical features of pCCD and CCD0 are exploited to avoid these problems in an efficient manner. The two resulting families of approximations, denoted as pCCD+DFT and CCD0+DFT, are shown to be able to describe static and dynamic correlation in standard benchmark calculations. Furthermore, it is also shown that CCD0+DFT lends itself to combination with correlation from the direct random phase approximation (dRPA). Inclusion of dRPA in the long-range via the technique of range-separation allows for the description of dispersion correlation, the remaining part of the correlation. Thus, when combined with the dRPA, CCD0+DFT can account for all three-types of electron correlation that are necessary to accurately describe molecular systems. Lastly, applications of CCD0+DFT to actinide chemistry are considered in this work. The accuracy of CCD0+DFT for predicting equilibrium geometries and vibrational frequencies of actinide molecules and ions is assessed and compared to that of well-established quantum chemical methods. For this purpose, the
f0 actinyl series (UO
22+, NpO
23+, PuO
24+), the isoelectronic NUN, and Thorium (ThO, ThO
2+) and Nobelium (NoO, NoO
2) oxides are studied. It is shown that the CCD0+DFT description of these species agrees with available experimental data and is comparable with the results given by the highest-level calculations that are possible for such heavy compounds while being, at least, an order of magnitude lower in computational cost.