For the description of ground-state correlation phenomena an accurate mapping of many-body quantum mechanics onto four particles is developed. The energy for a quantum system with no more than two-particle interactions may be expressed in terms of a two-particle reduced density matrix (2-RDM), but variational optimization of the 2-RDM requires that it corresponds to an N-particle wave function. We derive N-representability conditions on the 2-RDM that guarantee the validity of the uncertainty relations for all operators with two-particle interactions. One of these conditions is shown to be necessary and sufficient to make the RDM solutions of the dispersion condition equivalent to those from the contracted Schr\"odinger equation (CSE) [Mazziotti, Phys. Rev. A 57, 4219 (1998)]. In general, the CSE is a stronger N-representability condition than the dispersion condition because the CSE implies the dispersion condition as well as additional N-representability constraints from the Hellmann-Feynman theorem. Energy minimization subject to the representability constraints is performed for a boson model with 10, 30, and 75 particles. Even when traditional wave-function methods fail at large perturbations, the present method yields correlation energies within 2%.
For the description of ground-state correlation phenomena an accurate mapping of many-body quantum mechanics onto four particles is developed. The energy for a quantum system with no more than two-particle interactions may be expressed in terms of a two-particle reduced density matrix (2-RDM), but variational optimization of the 2-RDM requires that it corresponds to an N-particle wave function. We derive N-representability conditions on the 2-RDM that guarantee the validity of the uncertainty relations for all operators with two-particle interactions. One of these conditions is shown to be necessary and sufficient to make the RDM solutions of the dispersion condition equivalent to those from the contracted Schr"odinger equation (CSE) [Mazziotti, Phys. Rev. A 57, 4219 (1998)]. In general, the CSE is a stronger N-representability condition than the dispersion condition because the CSE implies the dispersion condition as well as additional N-representability constraints from the Hellmann-Feynman theorem. Energy minimization subject to the representability constraints is performed for a boson model with 10, 30, and 75 particles. Even when traditional wave-function methods fail at large perturbations, the present method yields correlation energies within 2%.
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.
A many-body theory for interacting bosons is developed within the framework of minimizing the ground-state energy with respect to the two-particle reduced-density matrix (2-RDM) subject to N-representability conditions. The N-representability conditions, which ensure that the 2-RDM may be derived from an N-particle wave function, are imposed through a hierarchy of positivity conditions where the p-positivity conditions restrict the metric matrices for p∕2-body operators to be positive semidefinite. Using two-positivity, we minimize the ground-state energies of 5−−10000 harmonically interacting bosons in a harmonic external potential. The energies and 2-RDMs obtained are in agreement with the exact solution except for round-off errors, which implies that for this class of boson interactions two-positivity conditions alone yield exact results for any interaction strength. The ground-state energies obtained at strong interactions are more accurate than many-body perturbative techniques by many orders of magnitude.
The contracted Schr\"odinger equation (CSE) technique through its direct determination of the two-particle reduced density matrix (2RDM) without the wave function may offer a fresh alternative to traditional many-body quantum calculations. Without additional information the CSE, also known as the density equation, cannot be solved for the 2RDM because it also requires a knowledge of the 4RDM. We provide theoretical foundations through a reconstruction theorem for recent attempts at generating higher RDMs from the 2RDM to remove the indeterminacy of the CSE. With Grassmann algebra a more concise representation for Valdemoro's reconstruction functionals [F. Colmenero, C. Perez del Valle, and C. Valdemoro, Phys. Rev. A 47, 971 (1993)] is presented. From the perspective of the particle-hole equivalence we obtain Nakatsuji and Yasuda's correction for the 4RDM formula [H. Nakatsuji and K. Yasuda, Phys. Rev. Lett. 76, 1039 (1996)] as well as a corrective approach for the 3RDM functional. A different reconstruction strategy, the ensemble representability method (ERM), is introduced to build the 3- and 4-RDMs by enforcing four-ensemble representability and contraction conditions. We derive the CSE in second quantization without Valdemoro's matrix contraction mapping and offer the first proof of Nakatsuji's theorem for the second-quantized CSE. Both the functional and ERM reconstruction strategies are employed with the CSE to solve for the energies and the 2RDMs of a quasispin model without wave functions. We elucidate the iterative solution of the CSE through an analogy with the power method for eigenvalue equations. Resulting energies of the CSE methods are comparable to single-double configuration-interaction (SDCI) energies, and the 2RDMs are more accurate by an order of magnitude than those from SDCI. While the CSE has been applied to systems with 14 electrons, we present results for as many as 40 particles. Results indicate that the 2RDM remains accurate as the number of particles increases. We also report a direct determination of excited-state 2RDMs through the CSE. By circumventing the wave function, the CSE presents new possibilities for treating electron correlation.
Here we present a many-body theory based on a solution of the N-representability problem in which the ground-state two-particle reduced density matrix (2-RDM) is determined directly without the many-particle wave function. We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM, which are required for its derivation from an N-particle density matrix, known as N-representability conditions. The approach produces a complete hierarchy of 2-RDM constraints that do not depend explicitly upon the higher RDMs or the wave function. By using the two-particle part of a unitary decomposition of higher-order constraint matrices, we can solve the energy minimization by semidefinite programming in a form where the low-rank structure of these matrices can be potentially exploited. We illustrate by computing the ground-state electronic energy and properties of the H_{8} ring.
Variationally calculating the ground state of a many-electron quantum system using only the two-electron reduced-density-matrix (2-RDM) requires N-representability conditions that constrain the 2-RDM to correspond to an N-electron wave function. A systematic hierarchy of N-representability conditions, known as p-positivity conditions, has been developed [D. A. Mazziotti and R. M. Erdahl, Phys. Rev. A 63, 042113 (2001)], and many-electron atoms and molecules in nonminimal basis sets have been solved with useful accuracy by a variational 2-RDM method with 2-positivity conditions [D. A. Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)]. This paper considers two forms of partial 3-positivity conditions, the lifting conditions and the T1∕T2 conditions, to further enhance the accuracy of the 2-RDM methods without the computational cost of full 3-positivity conditions. Variational 2-RDM methods with different N-representability constraints including 2-positivity conditions, the two types of partial 3-positivity conditions, as well as the complete 3-positivity conditions are applied to compute the ground state of the Lipkin spin model. The energies and 2-RDMs are compared to the results from full and truncated configuration interaction, many-body perturbation theory, and couple cluster theory with single and double excitations. Implications of using partial 3-positivity for variational 2-RDM calculations of many-electron atoms and molecules will be discussed.
Atomic and molecular ground-state energies are variationally determined by constraining the two-particle reduced density matrix (2-RDM) to satisfy positivity conditions. Because each positivity condition corresponds to correcting the ground-state energies for a class of Hamiltonians with two-particle interactions, these conditions collectively provide a new approach to many-body theory that, unlike perturbation theory, can capture significantly correlated phenomena including the multireference effects of potential-energy surfaces. The D, Q, and G conditions for the 2-RDM are extended through generalized lifting operators inspired from the formal solution of N-representability. These lifted conditions agree with the hierarchy of positivity conditions presented by Mazziotti and Erdahl [Phys. Rev. A 63, 042113 (2001)]. The connection between positivity and the formal solution explains how constraining higher RDMs to be positive semidefinite improves the N representability of the 2-RDM and suggests using pieces of higher positivity conditions that computationally scale like the D condition. With the D, Q, and G conditions as well as pieces of higher positivity the electronic energies for Be, LiH, H2O, and BH are computed through a primal-dual interior-point algorithm for positive semidefinite programming. The variational method produces potential-energy surfaces that are highly accurate even far from the equilibrium geometry where single-reference perturbation-based methods often fail to produce realistic energies.