The Journal of Chemical Physics

Cumulant expansion of the reduced density matrices

Werner Kutzelnigg, Debashis Mukherjee

Published Feb 8, 1999• 291 citations

Abstract

k -particle cumulants λk (for 2⩽k⩽n) corresponding to the k-particle reduced density matrices γk for an n-fermion system are defined via a generating function. The two-particle cumulant λ2 describes two-particle correlations (excluding exchange), λ3 genuine three-particle correlations etc. The properties of these cumulants are analyzed. Conditions for vanishing of certain λk are formulated. Necessary and sufficient for λ2=0 is the well-known idempotency condition γ2=γ for γ≡γ1. For λ3=0 to hold, a general necessary condition is Tr{2γ3−3γ2+γ}=0, for three special forms of the wave function (arbitrary two-electron state, antisymmetrized product of strongly orthogonal geminals on antisymmetrized geminal power wave function of extreme type) 2γ3−3γ2+γ=0 turns out to be necessary and sufficient. For a multiconfiguration self-consistent field wave function the only nonvanishing matrix elements of the cumulants are those where all labels refer to active (partially occupied) spin orbitals. Spin-free cumulants Λk corresponding to the spin-free reduced density matrices Γk are also defined and analyzed. The main interest in the density cumulants is in connection with the recently formulated normal ordering and the corresponding Wick theorem for arbitrary reference functions, but they are also useful for an analysis of electron correlation.

DOI: 10.1063/1.478189

Related papers

Spinfree formulation of reduced density matrices, density cumulants and generalised normal ordering

Werner Kutzelnigg, K. R. Shamasundar, Debashis Mukherjee

The k-particle reduced density matrices γ k for a singlet state or a spin-orientation averaged ensemble are invariant under the spin rotation group SU(2), i.e. they are independent of spin, even if they are expressed in terms of spin orbitals. The γ k contain the same information as their spinfree counterparts Γ k expressed in terms of spinfree orbitals. However, the relation between the two quantities is, by no means, elementary. The summation over spin is particularly demanding and surprisingly complicated for the k-particle density cumulants λ k and for the formulation of a generalised normal ordering with respect to an arbitrary reference state, including the generalised Wick theorem. Plausible extrapolations from small k turned out to be unjustified and some published expressions for k = 4 need to be revised. A compact general result, valid for all k has not yet been found.

Molecular PhysicsFeb 10, 201054 citations

Irreducible Brillouin conditions and contracted Schrödinger equations for n -electron systems. I. The equations satisfied by the density cumulants

Debashis Mukherjee, Werner Kutzelnigg

Two alternative conditions for the stationarity of the energy expectation value with respect to k-particle excitations are the k-particle Brillouin conditions BCk and the k-particle contracted Schrödinger equations, CSEk. These conditions express the k-particle density matrices γk in terms of density matrices of higher particle rank. The latter can be eliminated if one expresses the γk in terms of their cumulants λk, but this is not sufficient to make the BCk or CSEk separable (extensive), i.e., they are not expressible in terms of only connected diagrams. However, in a formulation based on the recently introduced general normal ordering with respect to arbitrary wave functions, the irreducible counterparts IBCk and ICSEk of the BCk and CSEk can be defined. They are easily evaluated explicitly in terms of the generalized Wick theorem for arbitrary wave functions, and they lead to equations for the direct construction of the cumulants λk, which are additively separable quantities and which scale linearly with the system size. The IBCk or the ICSEk are necessary conditions for γ and the λk to represent an exact n-fermionic eigenstate of the given Hamiltonian. To specify the desired state, additional conditions must be satisfied as well, e.g., the partial trace relations which relate λ2 to γ and γ2. The particle number and the total spin must be specified and n-representability conditions enter implicitly. While the nondiagonal elements of γ and the λk are determined by the IBCk or the ICSEk, the additional conditions mainly serve to fix the diagonal elements. A hierarchy of k-particle approximations is defined. It is based on the fact that the expansion in terms of cumulants λk can be truncated at any particle rank, which would not be possible for the density matrices γk. For closed-shell states the one-particle approximation agrees with Hartree–Fock.

The Journal of Chemical PhysicsFeb 1, 2001106 citations

Cumulant decomposition of reduced density matrices, multireference normal ordering, and Wicks theorem: A spin-free approach

K. R. Shamasundar

We propose a spin-free approach to the cumulant decomposition of reduced density matrices of singlet and spin-rotation or SU(2) invariant ensemble of nonsinglet states as in [W. Kutzelnigg and D. Mukherjee, J. Chem. Phys. 110, 2800 (1999); 116, 4787 (2002)]. We provide a simple recursive procedure to obtain expressions which relate different spin components of spin-orbital reduced density matrices and cumulants of such states to the spin-free counterparts. These results are used to set up a spin-summation procedure to arrive at a definition of spin-free cumulants of any order. Alternatively, an analytic formula for the spin-free form resulting from a spin summation involving product of two spin-orbital cumulants is derived and its utility in spin-free cumulant decomposition of reduced density matrices is demonstrated. This leads to suitable definitions of spin-free analog of multireference normal ordering and the associated Wicks theorem. The results of this formulation are expected to be useful in investigations of spin-free multireference internally contracted coupled-cluster methods where cumulant approximations to the active reduced density matrices are employed.

The Journal of Chemical PhysicsNov 4, 200949 citations

Irreducible Brillouin conditions and contracted Schrödinger equations for n-electron systems. II. Spin-free formulation

Werner Kutzelnigg, Debashis Mukherjee

Recently [W. Kutzelnigg and D. Mukherjee, Chem. Phys. Lett. 317, 567 (2000); D. Mukherjee and W. Kutzelnigg, J. Chem. Phys 114, 2047 (2001)] the irreducible k-particle Brillouin conditions IBCk and the irreducible k-particle contracted Schrödinger equations ICSEk were derived. These permit the definition of a hierarchy of k-particle approximations for the direct calculation of the cumulants λk of the k-particle density matrices. Now, the spin-free form of these conditions, appropriate for a spin-free Hamiltonian, is given. This is particularly useful for open-shell states. The definition of the cumulants of the reduced densities has to be generalized for these anyway, making use of irreducible tensor operators with respect to SU2. There are two alternative definitions of spin-free cumulants, of which the one in terms of spin-free reduced density matrices appears to be preferable. Alternatively to the straight spin-free formulation, we also present a theory in terms of spin-free operators adapted to the symmetric group. Partial trace relations that relate the elements of the cumulants of different particle rank are derived. There are partial trace relations for “exchange elements,” which are determined by the total spin quantum number S. From these relations the individual exchange elements of Λ2 can be obtained in special cases. This allows a simple formulation of the stationarity conditions for open-shell states of any spin multiplicity.

The Journal of Chemical PhysicsMar 22, 200256 citations

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

Jonathon P. Misiewicz, Justin M. Turney, Henry F. Schaefer

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.

Journal of Chemical Theory and ComputationAug 31, 202019 citations

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

Liguo Kong, Edward F. Valeev

We propose a novel interpretation of the reduced density matrix (RDM) and its cumulant that combines linear and exponential parametrizations of the wavefunction. Any n-particle RDM can be written as a weighted average of "configuration interaction" amplitudes. The corresponding n-particle cumulant is represented in terms of two types of contributions: "connected" (statistical averages of substitution amplitudes) and "disconnected" (cross-correlations of substitution amplitudes). A diagonal element of n-RDM represents the average occupation number of the orbital n-tuple. The diagonal elements of 2- and 3-cumulants take particularly elegant forms in the natural spin-orbital basis: they represent the covariances (correlated fluctuations) of the occupation numbers of the orbital pair and triples, respectively. Thus, the diagonal elements of the cumulants quantify the correlation between the orbital occupation numbers. Our interpretation is used to examine the weak to strong correlation transition in the "two electrons in two orbitals" problem.

The Journal of Chemical PhysicsJun 2, 201128 citations