Contracted Schrödinger equation: Determining quantum energies and two-particle density matrices without wave functions

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.