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

Chemical Physics Letters

Recursively generated linear constraints for variational two-particle reduced-density-matrix theory

Tamás Juhász, Neil Shenvi, David A. Mazziotti

Published Aug 2, 2007• 6 citations
Publisher
DOI
10.1016/j.cplett.2007.07.067

Related papers

Variational determination of the two-particle reduced density matrix on the GPU

Ward Poelmans, Brecht Verstichel, Dimitri Van Neck +1 more

•Jan 1, 2013•0 citations

Simple Hamiltonians which exhibit drastic failures by variational determination of the two-particle reduced density matrix with some well known N-representability conditions

真秀 中田, Maho Nakata, Bastiaan J. Braams +5 more

•Jan 1, 2006•0 citations

Reexamining the variational two-particle reduced density matrix for nuclear systems

J. G. Li, N. Michel, Zuo Wei +1 more

In most nuclear many-body methods, observables are calculated using many-body wave functions explicitly. The variational two-particle reduced density matrix method is one of the few exceptions to the rule. Ground-state energies of both closed-shell and open-shell nuclear systems can indeed be evaluated by minimizing a constrained linear functional of the two-particle reduced density matrix. However, it has virtually never been used in nuclear theory, because nuclear ground states were found to be well overbound, contrary to those of atoms and molecules. Consequently, we introduced new constraints in the nuclear variational two-particle reduced density matrix method, developed recently for atomic and molecular systems. Our calculations then show that this approach can provide a proper description of nuclear systems where only valence neutrons are included. For the nuclear systems where both neutrons and protons are active, however, the energies obtained with the variational two-particle reduced density matrix method are still overbound. The possible reasons for the noticed discrepancies and solutions to this problem will be discussed.

arXiv:2106.04103•Physical Review C•Jun 23, 2021•2 citations

Two algorithms for the lower bound method of reduced density matrix theory

Robert Erdahl

Reports on Mathematical Physics•Apr 1, 1979•79 citations

Active-Space N -Representability Constraints for Variational Two-Particle Reduced Density Matrix Calculations

Neil Shenvi, Artur F. Izmaylov

The ground-state energy of a system of fermions can be calculated by minimizing a linear functional of the two-particle reduced density matrix (2-RDM) if an accurate set of N-representability conditions is applied. In this Letter we introduce a class of linear N-representability conditions based on exact calculations on a reduced active space. Unlike wave-function-based approaches, the 2-RDM methodology allows us to combine information from calculations on different active spaces. By adding active-space constraints, we can iteratively improve our estimate for the ground-state energy. Applying our methodology to a 1D Hubbard model yields a significant improvement over traditional 2-positivity constraints with the same computational scaling.

Physical Review Letters•Nov 18, 2010•60 citations

Spurious zero-dispersion continua in full $(2,q)$-positivity: A linear hierarchy-depth lower bound from an exactly solvable counting operator

Maho Nakata

Variational two-particle reduced-density-matrix methods impose hierarchies of semidefinite constraints. We ask whether a fixed level of such a hierarchy can recover even an elementary many-fermion spectrum, using the dispersion test, in which a value is accepted as an eigenvalue when the minimized dispersion vanishes. For an exactly solvable operator that counts fermions in one half of the spin orbitals, whose spectrum consists of integers, we determine the complete set of accepted values at every level of the full hierarchy of particle-hole semidefinite constraints. Besides the integers, this set contains an entire interval of noninteger values unless the hierarchy level grows linearly with the particle number, and we give the exact threshold. The proof constructs explicit pseudo reduced density matrices that are diagonal in the occupation basis, obtained by continuing the occupation probabilities of a fixed-particle-number state to a noninteger particle number. Classical semidefinite constraints on occupation correlations, supplemented by nonnegative occupation-pattern probabilities, suffice to establish fermionic positivity. For this counting model the two relaxations have the same dispersion zero set, not necessarily the same feasible states, and a converse argument shows that no other false eigenvalues occur. The limitation is specific to the hierarchy: two-body inequalities that encode the integer spectrum would remove every false value but are not generated at a fixed level. We also verify that the full semidefinite constraints and their complete projection onto two-particle reduced-density-matrix space are equivalent in finite dimensions, so that the result applies equally to the polar formulation.

quant-pharXiv:2610.09698•Oct 8, 2026•0 citations
In most nuclear many-body methods, observables are calculated using many-body wave functions explicitly. The variational two-particle reduced density matrix method is one of the few exceptions to the rule. Ground-state energies of both closed-shell and open-shell nuclear systems can indeed be evaluated by minimizing a constrained linear functional of the two-particle reduced density matrix. However, it has virtually never been used in nuclear theory, because nuclear ground states were found to be well overbound, contrary to those of atoms and molecules. Consequently, we introduced new constraints in the nuclear variational two-particle reduced density matrix method, developed recently for atomic and molecular systems. Our calculations then show that this approach can provide a proper description of nuclear systems where only valence neutrons are included. For the nuclear systems where both neutrons and protons are active, however, the energies obtained with the variational two-particle reduced density matrix method are still overbound. The possible reasons for the noticed discrepancies and solutions to this problem will be discussed.
The ground-state energy of a system of fermions can be calculated by minimizing a linear functional of the two-particle reduced density matrix (2-RDM) if an accurate set of N-representability conditions is applied. In this Letter we introduce a class of linear N-representability conditions based on exact calculations on a reduced active space. Unlike wave-function-based approaches, the 2-RDM methodology allows us to combine information from calculations on different active spaces. By adding active-space constraints, we can iteratively improve our estimate for the ground-state energy. Applying our methodology to a 1D Hubbard model yields a significant improvement over traditional 2-positivity constraints with the same computational scaling.
Spurious zero-dispersion continua in full (2,q)(2,q)(2,q)-positivity: A linear hierarchy-depth lower bound from an exactly solvable counting operator
Variational two-particle reduced-density-matrix methods impose hierarchies of semidefinite constraints. We ask whether a fixed level of such a hierarchy can recover even an elementary many-fermion spectrum, using the dispersion test, in which a value is accepted as an eigenvalue when the minimized dispersion vanishes. For an exactly solvable operator that counts fermions in one half of the spin orbitals, whose spectrum consists of integers, we determine the complete set of accepted values at every level of the full hierarchy of particle-hole semidefinite constraints. Besides the integers, this set contains an entire interval of noninteger values unless the hierarchy level grows linearly with the particle number, and we give the exact threshold. The proof constructs explicit pseudo reduced density matrices that are diagonal in the occupation basis, obtained by continuing the occupation probabilities of a fixed-particle-number state to a noninteger particle number. Classical semidefinite constraints on occupation correlations, supplemented by nonnegative occupation-pattern probabilities, suffice to establish fermionic positivity. For this counting model the two relaxations have the same dispersion zero set, not necessarily the same feasible states, and a converse argument shows that no other false eigenvalues occur. The limitation is specific to the hierarchy: two-body inequalities that encode the integer spectrum would remove every false value but are not generated at a fixed level. We also verify that the full semidefinite constraints and their complete projection onto two-particle reduced-density-matrix space are equivalent in finite dimensions, so that the result applies equally to the polar formulation.