Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1207.0541 · Physical Review A

Significant conditions for the two-electron reduced density matrix from the constructive solution of N representability

David A. Mazziotti

Published Jun 7, 2012• 30 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreva.85.062507

Related papers

Variational two-electron reduced-density-matrix theory: Partial 3-positivity conditions for N -representability

Jeff R. Hammond, David A. Mazziotti

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 ${T}_{1}∕{T}_{2}$ 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.

Physical Review A•Jun 21, 2005•39 citations

Structure of Fermionic Density Matrices: Complete N -Representability Conditions

David A. Mazziotti

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.

arXiv:1112.5866•Physical Review Letters•Jun 27, 2012•221 citations

Variational minimization of atomic and molecular ground-state energies via the two-particle reduced density matrix

David A. Mazziotti

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, ${\mathrm{H}}_{2}\mathrm{O},$ 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.

Physical Review A•Jun 18, 2002•190 citations

Variational reduced-density-matrix method using three-particle N -representability conditions with application to many-electron molecules

David A. Mazziotti

Molecular two-electron reduced density matrices (2-RDMs) are computed variationally without the many-electron wave function by constraining the 2-RDM with a set of three-particle $N$-representability conditions known as three-positivity conditions. These constraints restrict four distinct three-particle probability distributions, which can be defined for any $N$-particle system, to be nonnegative. The variational calculation of the 2-RDM with full three-positivity conditions is implemented with the first-order semidefinite programming algorithm [D.A. Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)]. We also derive and implement a generalization of the ${T}_{2}$ representability condition, which is a subset of the three-positivity conditions. Energies and 2-RDMs are computed for several molecules as well as the nitrogen molecule at both equilibrium and nonequilibrium geometries. The ground-state energies for nitrogen are within $0.0015\phantom{\rule{0.3em}{0ex}}\mathrm{a.u.}$ of the values from full configuration interaction at all internuclear distances.

Physical Review A•Sep 1, 2006•123 citations

Pure- N -representability conditions of two-fermion reduced density matrices

David A. Mazziotti

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.

Physical Review A•Sep 29, 2016•67 citations

Quantum Many-Body Theory from a Solution of the N -Representability Problem

David A. Mazziotti

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.

arXiv:2304.08570•Physical Review Letters•Apr 12, 2023•41 citations
We recently presented a constructive solution to the NNN-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 NNN-electron quantum system [D. A. Mazziotti, Phys. Rev. Lett. (to be published)]. In this paper we provide additional details and derive further NNN-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)(2,q)(2,q)-positivity conditions, where the qqq indicates their derivation from the non-negativity of qqq-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)(2,q)(2,q)-positivity conditions to higher (2,q+1)(2,q+1)(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)(2,q)(2,q)-positivity conditions for q>3q>3q>3 require tensor decompositions of the coefficients in the model Hamiltonians. Subsets of the derived NNN-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.
Variationally calculating the ground state of a many-electron quantum system using only the two-electron reduced-density-matrix (2-RDM) requires NNN-representability conditions that constrain the 2-RDM to correspond to an NNN-electron wave function. A systematic hierarchy of NNN-representability conditions, known as ppp-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{T}_{1}∕{T}_{2}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 NNN-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.
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.
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,{\mathrm{H}}_{2}\mathrm{O},H2​O, 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.
Molecular two-electron reduced density matrices (2-RDMs) are computed variationally without the many-electron wave function by constraining the 2-RDM with a set of three-particle NNN-representability conditions known as three-positivity conditions. These constraints restrict four distinct three-particle probability distributions, which can be defined for any NNN-particle system, to be nonnegative. The variational calculation of the 2-RDM with full three-positivity conditions is implemented with the first-order semidefinite programming algorithm [D.A. Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)]. We also derive and implement a generalization of the T2{T}_{2}T2​ representability condition, which is a subset of the three-positivity conditions. Energies and 2-RDMs are computed for several molecules as well as the nitrogen molecule at both equilibrium and nonequilibrium geometries. The ground-state energies for nitrogen are within 0.0015a.u.0.0015\phantom{\rule{0.3em}{0ex}}\mathrm{a.u.}0.0015a.u. of the values from full configuration interaction at all internuclear distances.
We derive necessary conditions for any two-fermion reduced density matrix (2-RDM) to be representable by a pure NNN-fermion density matrix $\mathrm{\ensuremath{\Psi}}{\mathrm{\ensuremath{\Psi}}}^{*}$ where $\mathrm{\ensuremath{\Psi}}$ is the wave function. These pure NNN-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 NNN-representability conditions are nontrivially active for exact ground-state 2-RDMs and that they provide significant restrictions beyond the D, Q, and G ensemble NNN-representability conditions. Constraints on higher-order ppp-RDMs where p>2p>2p>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.