A novel interpretation of reduced density matrix and cumulant for electronic structure theories

On the notion of strong correlation in electronic structure theory

Strong correlation has been said to have many faces, and appears to have many synonyms of questionable suitability. In this work we aim not to define the term once and for all, but to highlight one possibility that is both rigorously defined and physically transparent, and remains so in reference to molecules and quantum lattice models. We survey both molecular examples - hydrogen systems (Hn, n = 2, 4, 6), Be2, H-He-H, and benzene - and the half-filled Hubbard model over a range of correlation regimes. Various quantities are examined including the extent of spin symmetry breaking in correlated single-reference wave functions, energetic ratios inspired by the Hubbard model and the Virial theorem, and metrics derived from the one- and two-electron reduced density matrices (RDMs). The trace and the square norm of the cumulant of the two-electron reduced density matrix capture what may well be defined as strong correlation. Accordingly, strong correlation is understood as a statistical dependence between two electrons, and is distinct from the concepts of "correlation energy" and more general than entanglement quantities that require a partitioning of a quantum system into distinguishable subspaces. This work enables us to build a bridge between a rigorous and quantifiable regime of strong electron correlation and more familiar chemical concepts such as anti-aromaticity in the context of Baird's rule.

Intruder-free cumulant-truncated driven similarity renormalization group second-order multireference perturbation theory

Accurate multireference electronic structure calculations are important for constructing potential energy surfaces. Still, even in the case of low-scaling methods, their routine use is limited by the steep growth of the computational and storage costs as the active space grows. This is primarily due to the occurrence of three- and higher-body density matrices or, equivalently, their cumulants. This work examines the effect of various cumulant truncation schemes on the accuracy of the driven similarity renormalization group second-order multireference perturbation theory. We test four different levels of three-body reduced density cumulant truncations that set different classes of cumulant elements to zero. Our test cases include the singlet-triplet gap of CH2, the potential energy curves of the XΣg+1 and AΣu+3 states of N2, and the singlet-triplet splittings of oligoacenes. Our results show that both relative and absolute errors introduced by these cumulant truncations can be as small as 0.5 kcal mol-1 or less. At the same time, the amount of memory required is reduced from O(NA6) to O(NA5), where NA is the number of active orbitals. No additional regularization is needed to prevent the intruder state problem in the cumulant-truncated second-order driven similarity renormalization group multireference perturbation theory methods.

Geometrical picture of the electron-electron correlation at the large-D limit

In electronic structure calculations, the correlation energy is defined as the difference between the mean field and the exact solution of the non relativistic Schrödinger equation. Such an error in the different calculations is not directly observable as there is no simple quantum mechanical operator, apart from correlation functions, that correspond to such quantity. Here, we use the dimensional scaling approach, in which the electrons are localized at the large-dimensional scaled space, to describe a geometric picture of the electronic correlation. Both, the mean field, and the exact solutions at the large-D limit have distinct geometries. Thus, the difference might be used to describe the correlation effect. Moreover, correlations can be also described and quantified by the entanglement between the electrons, which is a strong correlation without a classical analog. Entanglement is directly observable and it is one of the most striking properties of quantum mechanics and bounded by the area law for local gapped Hamiltonians of interacting many-body systems. This study opens up the possibility of presenting a geometrical picture of the electron-electron correlations and might give a bound on the correlation energy. The results at the large-D limit and at D = 3 indicate the feasibility of using the geometrical picture to get a bound on the electron-electron correlations.

Reduced Density Matrix Cumulants: The Combinatorics of Size-Consistency and Generalized Normal Ordering

Reduced density matrix cumulants play key roles in the theory of both reduced density matrices and multiconfigurational normal ordering. We present a new, simpler generating function for reduced density matrix cumulants that is formally identical with equating the coupled cluster and configuration interaction ansätze. This is shown to be a general mechanism to convert between a multiplicatively separable quantity and an additively separable quantity, as defined by a set of axioms. It is shown that both the cumulants of probability theory and the reduced density matrices are entirely combinatorial constructions, where the differences can be associated with changes in the notion of "multiplicative separability" for expectation values of random variables compared to reduced density matrices. We compare our generating function to that of previous works and criticize previous claims of probabilistic significance of the reduced density matrix cumulants. Finally, we present a simple proof of the generalized normal ordering formalism to explore the role of reduced density matrix cumulants therein. While the formalism can be used without cumulants, the combinatorial structure of expressing RDMs in terms of cumulants is the same combinatorial structure on cumulants that allows for a simple extended generalized Wick's theorem.

Signature of van der Waals interactions in the cumulant density matrix

Here we propose and implement a universal signature of the van der Waals interactions based on the cumulant part of the two-electron reduced density matrix (2-RDM). Due to the connected property of the cumulant, we can use it to detect the van der Waals interactions between two molecular moieties. In particular, we use the squared Frobenius norm of the cumulant of the 2-RDM, which has been previously shown to provide a size-extensive measure of the electron correlation. As two moieties are separated to infinity, the cumulant Frobenius norm exhibits an r-6 decay to its asymptotic limit, providing a density-based measure of the van der Waals interaction. We study this signature of van der Waals forces in a collection of small molecules of varying geometries. These computations agree with experimental trends of known literature values.

Density cumulant functional theory from a unitary transformation: N-representability, three-particle correlation effects, and application to O4(+)

A new approach to density cumulant functional theory is developed that derives density cumulant N-representability conditions from an approximate Fock space unitary transformation. We present explicit equations for the third- and fourth-order two-particle cumulant N-representability, as well as the second-order contributions that depend on the connected three-particle density cumulant. These conditions are used to formulate the ODC-13 method and the non-iterative (λ3) correction that employ an incomplete description of the fourth-order two-particle cumulant N-representability and the second-order three-particle correlation effects, respectively. We perform an analysis of the ODC-13 N-representability description for the dissociation of H2 and apply the ODC-13 method and the (λ3) correction to diatomic molecules with multiple bond character and the symmetry-breaking tetraoxygen cation (O4(+)). For the O4(+) molecule, the vibrational frequencies of the ODC-13(λ3) method do not exhibit spatial symmetry breaking and are in a good agreement with the recent infrared photodissociation experiment. We report the O4(+) equilibrium structure, harmonic frequencies, and dissociation energy computed using ODC-13(λ3) with a diffuse, core-correlated aug-cc-pCVTZ basis set.

Orbital-optimized density cumulant functional theory

In density cumulant functional theory (DCFT) the electronic energy is evaluated from the one-particle density matrix and two-particle density cumulant, circumventing the computation of the wavefunction. To achieve this, the one-particle density matrix is decomposed exactly into the mean-field (idempotent) and correlation components. While the latter can be entirely derived from the density cumulant, the former must be obtained by choosing a specific set of orbitals. In the original DCFT formulation [W. Kutzelnigg, J. Chem. Phys. 125, 171101 (2006)] the orbitals were determined by diagonalizing the effective Fock operator, which introduces partial orbital relaxation. Here we present a new orbital-optimized formulation of DCFT where the energy is variationally minimized with respect to orbital rotations. This introduces important energy contributions and significantly improves the description of the dynamic correlation. In addition, it greatly simplifies the computation of analytic gradients, for which expressions are also presented. We offer a perturbative analysis of the new orbital stationarity conditions and benchmark their performance for a variety of chemical systems.

Density cumulant functional theory: the DC-12 method, an improved description of the one-particle density matrix

Density cumulant functional theory (DCFT) is a theory that, in principle, can compute energies and properties exactly without a wavefunction. To accomplish this, the energy is expressed as an exact, known functional of the one-particle density matrix and two-particle density cumulant. The correlation contribution to the one-particle density matrix is obtained from the cumulant, to eliminate redundancy in the equations. The previous formulation of DCFT introduced this relationship in an approximate way, to obtain tractable equations. In this research, it is demonstrated that the correlation contribution to the one-particle density matrix can be extracted exactly from the cumulant, with minimal computational overhead and no increase in the asymptotic cost of the theory. We present numerical results, showing the improvements resulting from this reformulation (DC-12), and offer a perturbative analysis of the new equations to compare them to their predecessors.

Analytic gradients for density cumulant functional theory: the DCFT-06 model

Density cumulant functional theory (DCFT) is one of a number of nascent electron correlation methods that are derived from reduced density matrices and cumulants thereof, instead of the wavefunction. Deriving properties from the density cumulant naturally yields methods that are size extensive and size consistent. In this work, we derive expressions for the analytic gradient, with respect to an external perturbation, for the DCFT-06 variant of density cumulant functional theory. Despite the fact that the DCFT-06 energy functional is stationary with respect to the density cumulant, the analytic gradients of the energy require the solution of perturbation-independent equations for both orbital and cumulant response. These two sets of linear response equations are coupled in nature and are solved iteratively with the solution of orbital and cumulant response equations each macroiteration, exhibiting rapid convergence. The gradients are implemented and benchmarked against coupled cluster theory with single and double excitations (CCSD) and CCSD with perturbative triple excitations [CCSD(T)], as well as accurate empirically corrected experimental data, for a test set comprising 15 small molecules. For most of the test cases, results from DCFT-06 are closer to CCSD(T) and empirical data than those from CCSD. Although the total energy and analytic gradient have the same asymptotic scaling, the present experience shows that the computational cost of the gradient is significantly lower.