We present a constructive solution to the N-representability problem: a full characterization of the conditions for constraining the two-electron reduced density matrix to represent an N-electron density matrix. Previously known conditions, while rigorous, were incomplete. Here, we derive a hierarchy of constraints built upon (i) the bipolar theorem and (ii) tensor decompositions of model Hamiltonians. Existing conditions D, Q, G, T1, and T2, known classical conditions, and new conditions appear naturally. Subsets of the conditions are amenable to polynomial-time computations of strongly correlated systems.
We recently presented a constructive solution to the
N-representability problem of the two-electron reduced density matrix (2-RDM)---a systematic approach to constructing complete conditions to ensure that the 2-RDM represents a realistic
N-electron quantum system [D. A. Mazziotti, Phys. Rev. Lett. (to be published)]. In this paper we provide additional details and derive further
N-representability conditions on the 2-RDM that follow from the constructive solution. The resulting conditions can be classified into a hierarchy of constraints, known as the
(2,q)-positivity conditions, where the
q indicates their derivation from the non-negativity of
q-body operators. In addition to the known T1 and T2 conditions, we derive another class of (2,3)-positivity conditions. We also derive 3 classes of (2,4)-positivity conditions, 6 classes of (2,5)-positivity conditions, and 24 classes of (2,6)-positivity conditions. The constraints obtained can be divided into two general types: (i) lifting conditions, that is, conditions which arise from lifting lower
(2,q)-positivity conditions to higher
(2,q+1)-positivity conditions, and (ii) pure conditions, that is, conditions which cannot be derived from a simple lifting of the lower conditions. All of the lifting conditions and the pure
(2,q)-positivity conditions for
q>3 require tensor decompositions of the coefficients in the model Hamiltonians. Subsets of the derived
N-representability conditions can be employed with the previously known conditions to achieve polynomially scaling calculations of ground-state energies and 2-RDMs of many-electron quantum systems even in the presence of strong electron correlation.
We derive necessary conditions for any two-fermion reduced density matrix (2-RDM) to be representable by a pure
N-fermion density matrix $\mathrm{\ensuremath{\Psi}}{\mathrm{\ensuremath{\Psi}}}^{*}$ where $\mathrm{\ensuremath{\Psi}}$ is the wave function. These pure
N-representability conditions of the 2-RDM are important because they provide stringent constraints beyond those from the Pauli and the generalized Pauli constraints on the structure of many-fermion 2-RDMs and their wave functions. The pure 2-RDM conditions are derived as generalized Pauli constraints on effective one-fermion reduced density matrices (1-RDMs) generated by the removal or addition of a fermion from the wave function. Computationally, we show for four-electron molecules that the derived pure
N-representability conditions are nontrivially active for exact ground-state 2-RDMs and that they provide significant restrictions beyond the D, Q, and G ensemble
N-representability conditions. Constraints on higher-order
p-RDMs where
p>2 are derived in a similar fashion. The constraints have potentially significant applications to computing strongly correlated many-fermion states with enhanced accuracy and decreased computational complexity.
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.
For a quantum-mechanical system of
N identical fermions, the
N-representability problem is the problem of recognizing whether, for a given pth-order reduced density matrix ${\ensuremath{\Gamma}}^{(p)}(12\ensuremath{\cdots}p|{1}^{\ensuremath{'}}{2}^{\ensuremath{'}}\ensuremath{\cdots}{p}^{\ensuremath{'}})$, there exists an antisymmetric
N-particle wave function $\ensuremath{\Psi}(12\ensuremath{\cdots}N)$ such that ${\ensuremath{\Gamma}}^{(p)}(12\ensuremath{\cdots}p|{1}^{\ensuremath{'}}{2}^{\ensuremath{'}}\ensuremath{\cdots}{p}^{\ensuremath{'}})=({N}{p})\ensuremath{\int}\ensuremath{\Psi}(12\ensuremath{\cdots}N)\ifmmode\times\else\texttimes\fi{}{\ensuremath{\Psi}}^{*}({1}^{\ensuremath{'}}{2}^{\ensuremath{'}}\ensuremath{\cdots}{p}^{\ensuremath{'}},p+1\ensuremath{\cdots}N){d}_{{\ensuremath{\tau}}_{p+1}}\ensuremath{\cdots}{d}_{{\ensuremath{\tau}}_{N}}$. It is shown that if the Hamiltonian of a system is time-reversal invariant, and the number of particles,
N, is even, the necessary and sufficient condition that an approximate first-order density matrix corresponding to a nondegenerate energy eigenstate be
N-representable is that its natural spin-orbital occupation numbers be equal in pairs.
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.
The direct variational optimization of the two-electron reduced density matrix (2RDM) can provide a reference-independent description of the electronic structure of many-electron systems that naturally capture strong or nondynamic correlation effects. Such variational 2RDM approaches can often provide a highly accurate description of strong electron correlation, provided that the 2RDMs satisfy at least partial three-particle N-representability conditions (e.g., the T2 condition). However, recent benchmark calculations on hydrogen clusters [N. H. Stair and F. A. Evangelista, J. Chem. Phys. 153, 104108 (2020)] suggest that even the T2 condition leads to unacceptably inaccurate results in the case of two- and three-dimensional clusters. We demonstrate that these failures persist under the application of full three-particle N-representability conditions (3POS). A variety of correlation metrics are explored in order to identify regimes under which 3POS calculations become unreliable, and we find that the relative squared magnitudes of the cumulant three- and two-particle reduced density matrices correlate reasonably well with the energy error in these systems. However, calculations on other molecular systems reveal that this metric is not a universal indicator for the reliability of the reduced-density-matrix theory with 3POS conditions.