Cumulant reconstruction functionals

Using cumulant reconstruction functionals A3[Ai, A2] and A4[Ai, A2], one can certainly derive closed, nonlinear equations for the elements of Ai and A2, which could be solved using an iterative procedure that does not exploit the reconstruction functionals at each iteration. Of the RDM reconstruction functionals derived to date, several [7, 8, 11] utilize the cumulant decompositions in Eqs. (25c) and (25d) to obtain the unconnected portions of D3 and D4 exactly (in terms of the lower-order RDMs), then use many-body perturbation theory to estimate the connected parts A3 and A4 in terms of Aj and A2, the latter essentially serving as a renormalized pair interaction. Reconstruction functionals of this type are equally useful in solving ICSE(l) and ICSE(2), but the reconstruction functionals introduced by Valdemoro and co-workers [25, 26] cannot be used to solve the ICSEs because they contain no connected terms in D3 or D4 (and thus no contributions to A3 or A4). 

The reconstruction functionals, derived in the previous section through the particle-hole duality, may also be produced through the theory of cumulants [21,22,24,26,39,55-57]. We begin by constructing a functional whose derivatives with respect to probe variables generate the reduced density matrices in second quantization. Because we require that additional derivatives increase the number of second quantization operators, we are led to the following exponential form ... 

Rosina s theorem states that for an unspecified Hamiltonian with no more than two-particle interactions the ground-state 2-RDM alone has sufficient information to build the higher ROMs and the exact wavefunction [20, 51]. Cumulants allow us to divide the reconstruction functional into two parts (i) an unconnected part that may be written as antisymmetrized products of the lower RDMs, and (ii) a connected part that cannot be expressed as products of the lower RDMs. As shown in the previous section, cumulant theory alone generates all of the unconnected terms in RDM reconstruction, but cumulants do not directly indicate how to compute the connected portions of the 3- and 4-RDMs from the 2-RDM. In this section we discuss a systematic approximation of the connected (or cumulant) 3-RDM [24, 26]. 

The theory of cumulants allows us to partition an RDM into contributions that scale differently with the number N of particles. Because aU of the particles are connected by interactions, the cumulant RDMs scale linearly with the number N of particles. The unconnected terms in the p-RDM reconstruction formulas scale between N and W according to the number of connected RDMs in the wedge product. For example, the term scales as NP since all p particles are statistically independent of each other. By examining the scaling of terms with N in the contraction of higher reconstruction functionals, we may derive an important set of relations for the connected RDMs. 

Much of the recent literature on RDM reconstruction functionals is couched in terms of cumulant decompositions [13, 27-38]. Insofar as the p-RDM represents a quantum mechanical probability distribution for p-electron subsystems of an M-electron supersystem, the RDM cumulant formalism bears much similarity to the cumulant formalism of classical statistical mechanics, as formalized long ago by by Kubo [39]. (Quantum mechanics introduces important differences, however, as we shall discuss.) Within the cumulant formalism, the p-RDM is decomposed into connected and unconnected contributions, with the latter obtained in a known way from the lower-order -RDMs, q < p. The connected part defines the pth-order RDM cumulant (p-RDMC). In contrast to the p-RDM, the p-RDMC is an extensive quantity, meaning that it is additively separable in the case of a composite system composed of noninteracting subsystems. (The p-RDM is multiphcatively separable in such cases [28, 32]). The implication is that the RDMCs, and the connected equations that they satisfy, behave correctly in the limit of noninteracting subsystems by construction, whereas a 2-RDM obtained by approximate solution of the CSE may fail to preserve extensivity, or in other words may not be size-consistent [40, 42]. 

The CSE allows us to recast A-representability as a reconstruction problem. If we knew how to build from the 2-RDM to the 4-RDM, the CSE in Eq. (12) furnishes us with enough equations to solve iteratively for the 2-RDM. Two approaches for reconstruction have been explored in previous work on the CSE (i) the explicit representation of the 3- and 4-RDMs as functionals of the 2-RDM [17, 18, 20, 21, 29], and (ii) the construction of a family of higher 4-RDMs from the 2-RDM by imposing ensemble representability conditions [20]. After justifying reconstmction from the 2-RDM by Rosina s theorem, we develop in Sections III.B and III.C the functional approach to the CSE from two different perspectives—the particle-hole duality and the theory of cumulants. 

The 2-RDM is automatically antisymmetric, but it may require an adjustment of the trace to correct the normalization. The functionals in Table I from cumulant theory allow us to approximate the 3- and the 4-RDMs from the 2-RDM and, hence, to iterate with the contracted power method. Because of the approximate reconstruction the contracted power method does not yield energies that are strictly above the exact energy. As in the full power method the updated 2-RDM in Eq. (116) moves toward the eigenstate whose eigenvalue has the largest magnitude. 

Importantly, the anti-Hermitian CSE may be evaluated through second order of a renormalized perturbation theory even when the cumulant 3-RDM is neglected in the reconstruction. The anti-Hermitian part of the CSE [27, 31, 63] is the stationary condition for two-body unitary transformations of the A-particle wave-function [31, 32], and hence the two-body unitary transformations may easily be evaluated with the anti-Hermitian CSE and RDM reconstruction without the many-electron Schrodinger equation. The contracted Schrodinger equation in conjunction with the concepts of reconstruction and purification provides a new, important approach to computing the 2-RDM directly without the many-electron wavefunction. 

Under this definition, the E function is characterized by an experimental 5 factor, which can be estimated experimentally (Saad et al, 2001). Note that the cumulative B factor applied in the final reconstruction is a composite of experimental and computational causes. The computational B factor is attributable to additional blurring effects such as inaccuracy in determining the orientation of particles, which could also be described by a Gaussian function type. The dampening of the image contrast by the... 

The defocus values for electron micrographs can be readily estimated on the basis of the CTF rings visible in the incoherently averaged Fourier transforms of individual particle images (Zhou et at., 1994, 1996). This method has become a routine practice universally adapted for the initial evaluation and determination of CTF parameters as defined in Eqs. (2) and (3). So far, the determination of the cumulative B factors for 3-D reconstruction has been somewhat ad hoc. The cumulative B factor used in these studies is determined either by trial and error with the initial value derived from previous results, or from the incoherently averaged Fourier transforms of particle images. Different approaches have been adopted to make corrections for the CTF and the function of the micrographs. They differ in the steps where these corrections are made and in whether or how the E function is corrected. 

---