We consider, for systems of indistinguishable fermions, approximate reconstruction of the three- and four-particle reduced density matrices (RDMs) from the one- and two-particle RDMs, $\mathit{\ensuremath{\gamma}}$ and $\mathit{\ensuremath{\Gamma}}.$ Our ansatz for reconstructing the four-particle RDM is the linear combination $a(\mathit{\ensuremath{\Gamma}}\ensuremath{\wedge}\mathit{\ensuremath{\Gamma}})+b(\mathit{\ensuremath{\gamma}}\ensuremath{\wedge}\mathit{\ensuremath{\gamma}}\ensuremath{\wedge}\mathit{\ensuremath{\Gamma}})+c(\mathit{\ensuremath{\gamma}}\ensuremath{\wedge}\mathit{\ensuremath{\gamma}}\ensuremath{\wedge}\mathit{\ensuremath{\gamma}}\ensuremath{\wedge}\mathit{\ensuremath{\gamma}}),$ where ``\ensuremath{\wedge}'' denotes the antisymmetrized (Grassmann) product. This is a generalization of reconstruction functionals employed recently to perform direct RDM calculations without wave functions via the contracted Schr\"odinger equation. Here we consider relationships between the parameters a, b, and c that are required in order for the reconstruction functionals to respect the hierarchy of contraction relations between RDMs. To this end we establish several general theorems concerning contractions of antisymmetrized tensor products of $\mathit{\ensuremath{\gamma}},$ $\mathit{\ensuremath{\Gamma}},$ and various products thereof. The accuracy of proposed reconstruction functionals is evaluated using accurate density matrices for the ground state of Be.
We consider, for systems of indistinguishable fermions, approximate reconstruction of the three- and four-particle reduced density matrices (RDMs) from the one- and two-particle RDMs, \ensuremathγ and \ensuremathΓ. Our ansatz for reconstructing the four-particle RDM is the linear combination a(\ensuremathΓ\ensuremath∧\ensuremathΓ)+b(\ensuremathγ\ensuremath∧\ensuremathγ\ensuremath∧\ensuremathΓ)+c(\ensuremathγ\ensuremath∧\ensuremathγ\ensuremath∧\ensuremathγ\ensuremath∧\ensuremathγ), where ``\ensuremath{\wedge}'' denotes the antisymmetrized (Grassmann) product. This is a generalization of reconstruction functionals employed recently to perform direct RDM calculations without wave functions via the contracted Schr"odinger equation. Here we consider relationships between the parameters a, b, and c that are required in order for the reconstruction functionals to respect the hierarchy of contraction relations between RDMs. To this end we establish several general theorems concerning contractions of antisymmetrized tensor products of \ensuremathγ,\ensuremathΓ, and various products thereof. The accuracy of proposed reconstruction functionals is evaluated using accurate density matrices for the ground state of Be.
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.
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%.
Preface. I: Properties of Reduced Density Matrices. 1. RDMs: How Did We Get Here? A.J. Coleman. 2. Some Theorems on Uniqueness and Reconstruction of Higher-Order Density Matrices M. Rosina. 3. Cumulant Expansions of Reduced Densities, Reduced Density Matrices, and Green's Functions P. Ziesche. 4. On Calculating Approximate and Exact Density Matrices R. Erdahl, B. Jin. II: The Contracted Schroedinger Equation. 5. Density Equation Theory in Chemical Physics H. Nakatsuji. 6. Critical Questions Concerning Iterative Solution of the Contracted Schroedinger Equation C. Valdemoro, L.M. Tel, E. Perez-Romero 7. Cumulants and the Contracted Schroedinger Equation D.A. Mazziotti.. III: Density Matrix Functional Theory. 8. Natural Orbital Functional Theory S. Goedecker, C.J. Umrigar. 9. The Pair Density in Approximate Density Functional Theory: The Hidden Agent N.T. Maitra, K. Burke. 10. Functional N-representability in Density Matrix and Density Functional Theory: An Illustration for Hooke's Atom E.V. Ludena, V. Karasiev, A. Artemiev, D. Gomez.. IV: Electron Intracule and Extracule Densities. 11. Intracule and Extracule Densities: Historical Perspectives and Future Prospects E. Valderrama, J.M. Ugalde, R.J. Boyd. 12. Topology of Electron Correlation J. Cioslowski. 13. Electron-Pair Densities of Atoms T. Koga. Index.
The solution of N-representability for the two-particle reduced density matrix (2-RDM) is intimately connected through the contracted Schr\"odinger equation with the problem of reconstructing the 3- and 4-RDM's from the 2-RDM [D. A. Mazziotti, Phys. Rev. A 57, 4219 (1998)]. We derive an explicit formula for building the 3-RDM from the 2-RDM which significantly improves the accuracy of the 3-RDM cumulant expansion. Calculations indicate that the improvement is better than that achieved by the correction of Nakatsuji and Yasuda (NY) [Phys. Rev. Lett. 76, 1039 (1996); Phys. Rev. A 56, 2648 (1997)]. We also present a more efficient formulation of the NY correction. Both the present and NY reconstructions are shown to satisfy the particle-hole duality. Comparison of the reconstruction methods is made through several molecules LiH, BeH2,BH3,H2O, and CH4, as well as systems with more general pairwise correlation.