Two-electron reduced density matrices (2-RDMs) of many-electron molecules are directly determined without calculation of their wave functions by solving the anti-Hermitian contracted Schrödinger equation. Approximation of the 3-RDM in the anti-Hermitian contracted Schrödinger equation by a corrected cumulant expansion [Mazziotti, Phys. Rev. A 60, 3618 (1999)] permits the direct calculation of the energy and 2-RDM with many high-order correlation effects included. The method is illustrated for the molecules BeH2, H2O, NH3, CH4, and CO as well as the dissociation of BH. Correlation energies are obtained within 95%-100% of full-configuration interaction, and 2-RDMs very nearly satisfy known N-representability conditions.
Two-electron reduced density matrices (2-RDMs) of many-electron molecules are directly determined without calculation of their wave functions by solving the anti-Hermitian contracted Schrödinger equation. Approximation of the 3-RDM in the anti-Hermitian contracted Schrödinger equation by a corrected cumulant expansion [Mazziotti, Phys. Rev. A 60, 3618 (1999)] permits the direct calculation of the energy and 2-RDM with many high-order correlation effects included. The method is illustrated for the molecules BeH2, H2O, NH3, CH4, and CO as well as the dissociation of BH. Correlation energies are obtained within 95%-100% of full-configuration interaction, and 2-RDMs very nearly satisfy known N-representability conditions.
A recent advance in the theory of the contracted Schr\"odinger equation (CSE), in which only the anti-Hermitian part of the equation is solved, permits the direct determination of ground-state two-electron reduced density matrices (2-RDM's) that yield 95%--100% of the correlation energy of atoms and molecules [D. A. Mazziotti, Phys. Rev. Lett. 97, 143002 (2006)]. Here we discuss in detail the anti-Hermitian contracted Schr\"odinger equation (ACSE) and its comparison to the CSE with regard to cumulant reconstruction of RDM's, the role of Nakatsuji's theorem, and the structure of the wave function. The ACSE is also formulated in the Heisenberg representation and related to canonical diagonalization. The solution of the ACSE is illustrated with a variety of molecules including H2O, CH2, NH4+, HF, and N2, and potential energy and dipole-moment surfaces are computed for boron hydride in a polarized double-$\ensuremath{\zeta}$ basis set. The computed 2-RDM's very closely satisfy known N-representability conditions.
Two-electron reduced density matrices (2-RDMs) have recently been directly calculated by solving the anti-Hermitian contracted Schr\"odinger equation (ACSE) to obtain 95--100 of the ground-state correlation energy of atoms and molecules with the accuracy increasing with the size of the one-electron basis set [Mazziotti, Phys. Rev. Lett. 97, 143002 (2006).] In this paper, the ACSE method is extended to treat strong multireference correlation effects that are often important at nonequilibrium molecular geometries. While previous ACSE calculations have employed an initial 2-RDM from the Hartree-Fock method, we initialize the solution of the ACSE with a 2-RDM guess from a multiconfiguration self-consistent field calculation. Applications are made to multireference correlation in the potential energy surfaces of the molecules HF, H2O, and C2 in polarized valence double-zeta basis sets while N2 is treated in polarized valence double- and triple-zeta basis sets. Accurate ground-state energies and 1-RDM occupation numbers are obtained at both equilibrium and nonequilibrium geometries where the energies are within a few millihartrees of those from full configuration interaction.
For 50 years, progress toward the direct calculation of the ground-state two-electron reduced density matrix (2-RDM) was stymied from an inability to constrain the 2-RDM to represent an N-electron wave function. Recent advances in theory and optimization have realized the direct calculation of the 2-RDM. A variational 2-RDM procedure, using first-order semidefinite programming, has been shown to capture multireference correlation effects important at nonequilibrium geometries [Mazziotti, Phys. Rev. Lett. 2004, 93, No. 213001]. This method emerged from research on a nonvariational calculation of the 2-RDM by the contracted Schrödinger equation. Both approaches will be discussed and illustrated.
Correlation energies and reduced density matrices (RDMs) of atoms and molecules are directly computed by solving the 1,3-contracted Schr\"odinger equation (1,3-CSE). The solution of the 1,3-CSE synthesizes two optimization strategies recently employed for the direct determination of the 2-RDM: (i) variational minimization of the energy with respect to a 2-RDM constrained by positivity conditions [D. A. Mazziotti, Phys. Rev. A 65, 062511 (2002)] and (ii) the contracted power method for solving the 2,4-CSE [D. A. Mazziotti, J. Chem. Phys. 116, 1239 (2002)]. While both the 3- and the 4-RDMs in the 2,4-CSE are reconstructed from the 2-RDM by cumulant expansions, similar techniques cannot be directly applied to the 1,3-CSE because constructing the 2-RDM from the 1-RDM with cumulant theory does not improve upon the mean-field approximation. We, however, establish a unique mapping from the 1-RDM to the 2-RDM by searching for the 2-RDM, constrained by contraction and N-representability conditions, which minimizes the energy. The 2-RDM constrained search is practically implemented through recent advances in positive semidefinite programming. With the variational reconstruction of the 2-RDM and a cumulant reconstruction of the 3-RDM, the 1,3-CSE may be solved via a contracted power method for the ground-state energy and RDMs. The initial RDMs, it is shown, need not be N representable for the contracted power method to converge; this allows us to choose the original RDMs from a variational calculation with approximate N-representability conditions on the 2-RDM. Application of the 1,3-CSE algorithm to atoms and molecules yields highly accurate correlation energies both near and far from equilibrium geometries.