The variational two-electron reduced-density-matrix (2-RDM) method allows for the computation of accurate ground-state energies and 2-RDMs of atoms and molecules without the explicit construction of an N-electron wave function. While previous work on variational 2-RDM theory has focused on calculating full configuration-interaction energies, this work presents the first application toward approximating multiconfiguration self-consistent-field (MCSCF) energies via low-rank restrictions on the 1- and 2-RDMs. The 2-RDM method with two- or three-particle N-representability conditions reduces the exponential active-space scaling of MCSCF methods to a polynomial scaling. Because the first-order algorithm [Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)] represents each form of the 1- and 2-RDMs by a matrix factorization, the RDMs are readily defined to have a low rank rather than a full rank by setting the matrix factors to be rectangular rather than square. Results for the potential energy surfaces of hydrogen fluoride, water, and the nitrogen molecule show that the low-rank 2-RDM method yields accurate approximations to the MCSCF energies. We also compute the energies along the symmetric stretch of a 20-atom hydrogen chain where traditional MCSCF calculations, requiring more than 17x10(9) determinants in the active space, could not be performed.
The Pauli exclusion principle requires the spectrum of the occupation numbers of the one-electron reduced density matrix (1-RDM) to be bounded by one and zero. However, for a 1-RDM from a wave function, there exist additional conditions on the spectrum of occupation numbers, known as pure
N-representability conditions or generalized Pauli conditions. For atoms and molecules, we measure through a Euclidean-distance metric the proximity of the 1-RDM spectrum to the facets of the convex set (polytope) generated by the generalized Pauli conditions. For the ground state of any spin symmetry, as long as time-reversal symmetry is considered in the definition of the polytope, we find that the 1-RDM's spectrum is pinned to the boundary of the polytope. In contrast, for excited states, we find that the 1-RDM spectrum is not pinned. Proximity of the 1-RDM to the boundary of the polytope provides a measurement and classification of electron correlation and entanglement within the quantum system. For comparison, this distance to the boundary of the generalized Pauli conditions is also compared to the distance to the polytope of the traditional Pauli conditions, and the distance to the nearest 1-RDM spectrum from a Slater determinant. We explain the difference in pinning in the ground- and excited-state 1-RDMs through a connection to the
N-representability conditions of the two-electron reduced density matrix.
Nuclear dynamics is incorporated into an efficient density matrix formalism of electronic dynamics which has been applied to molecular systems containing thousands of atoms. The formalism for the combined dynamics of electrons and nuclei is derived from the Dirac–Frenkel variational principle. The single electron reduced density matrices and the Glauber coherent states are used for the electronic and nuclear degrees of freedom, respectively. The new formalism is applicable to simulate the dynamics of large molecular systems. As an illustration of its validity, the formalism is employed to calculate the electron and nuclei dynamics of hydrogen molecules.