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:1707.01022 · The Journal of Chemical Physics

Method for making 2-electron response reduced density matrices approximately N-representable

Caitlin Lanssens, Paul W. Ayers, Dimitri Van Neck, Stijn De Baerdemacker, Klaas Gunst, Patrick Bultinck

Published Feb 22, 2018• 14 citations
arXivHTMLPDF

Abstract

In methods like geminal-based approaches or coupled cluster that are solved using the projected Schrödinger equation, direct computation of the 2-electron reduced density matrix (2-RDM) is impractical and one falls back to a 2-RDM based on response theory. However, the 2-RDMs from response theory are not N-representable. That is, the response 2-RDM does not correspond to an actual physical N-electron wave function. We present a new algorithm for making these non-N-representable 2-RDMs approximately N-representable, i.e., it has the right symmetry and normalization and it fulfills the P-, Q-, and G-conditions. Next to an algorithm which can be applied to any 2-RDM, we have also developed a 2-RDM optimization procedure specifically for seniority-zero 2-RDMs. We aim to find the 2-RDM with the right properties which is the closest (in the sense of the Frobenius norm) to the non-N-representable 2-RDM by minimizing the square norm of the difference between this initial response 2-RDM and the targeted 2-RDM under the constraint that the trace is normalized and the 2-RDM, Q-matrix, and G-matrix are positive semidefinite, i.e., their eigenvalues are non-negative. Our method is suitable for fixing non-N-representable 2-RDMs which are close to being N-representable. Through the N-representability optimization algorithm we add a small correction to the initial 2-RDM such that it fulfills the most important N-representability conditions.

DOI
10.1063/1.4994618

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

Quantum Chemistry without Wave Functions: Two-Electron Reduced Density Matrices

David A. Mazziotti

For 50 years, progress toward the direct calculation of the ground-state two-electron reduced density matrix (2-RDM) was stymied from an inability to constrain the 2-RDM to represent an N-electron wave function. Recent advances in theory and optimization have realized the direct calculation of the 2-RDM. A variational 2-RDM procedure, using first-order semidefinite programming, has been shown to capture multireference correlation effects important at nonequilibrium geometries [Mazziotti, Phys. Rev. Lett. 2004, 93, No. 213001]. This method emerged from research on a nonvariational calculation of the 2-RDM by the contracted Schrödinger equation. Both approaches will be discussed and illustrated.

Accounts of Chemical Research•Feb 9, 2006•128 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

Spin and symmetry adaptation of the variational two-electron reduced-density-matrix method

Gergely Gidofalvi, David A. Mazziotti

The variational two-electron reduced-density-matrix (2-RDM) method computes the ground-state energy and 2-RDM of an atom or molecule without calculation of the many-electron wave function. Recently, the computational efficiency of the 2-RDM method has been significantly enhanced through the use of a first-order algorithm for semidefinite programming [Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)]. In this paper we develop a spin- and symmetry-adapted formulation of the method that further improves its efficiency by incorporating both the spin and spatial symmetries of many-electron atoms and molecules. While previous work on density-matrix symmetry focused on only one form of the 2-RDM, the variational method employs three different forms of the 2-RDM, known as the $D$, $Q$, and $G$ matrices, to restrict the 2-RDM to be approximately $N$-representable, that is representable by an $N$-electron wave function. We apply spin symmetries to the three forms of the 2-RDM, each of which breaks into four diagonal spin-blocks, namely one singlet and three triplet blocks. If the molecules have point-group symmetry, each of the 2-RDMs may be further subdivided into smaller diagonal blocks according to the spatial symmetry of the basis functions. The subdivision of the 2-RDMs into diagonal blocks generates significant computational savings in both floating-point operations and memory storage. Calculations illustrate the computational savings. Spin adaptation also enforces the correct expectation value of the ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{S}}^{2}$ operator, which in earlier work is applied as a separate constraint.

Physical Review A•Nov 14, 2005•69 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

Rank restriction for the variational calculation of two-electron reduced density matrices of many-electron atoms and molecules

Kasra Naftchi-Ardebili, Nathania W. Hau, David A. Mazziotti

Variational minimization of the ground-state energy as a function of the two-electron reduced density matrix (2-RDM), constrained by necessary $N$-representability conditions, provides a polynomial-scaling approach to studying strongly correlated molecules without computing the many-electron wave function. Here we introduce a route to enhancing necessary conditions for $N$ representability through rank restriction of the 2-RDM. Rather than adding computationally more expensive $N$-representability conditions, we directly enhance the accuracy of two-particle (2-positivity) conditions through rank restriction, which removes degrees of freedom in the 2-RDM that are not sufficiently constrained. We select the rank of the particle-hole 2-RDM by deriving the ranks associated with model wave functions, including both mean-field and antisymmetrized geminal power (AGP) wave functions. Because the 2-positivity conditions are exact for quantum systems with AGP ground states, the rank of the particle-hole 2-RDM from the AGP ansatz provides a minimum for its value in variational 2-RDM calculations of general quantum systems. To implement the rank-restricted conditions, we extend a first-order algorithm for large-scale semidefinite programming. The rank-restricted conditions significantly improve the accuracy of the energies; for example, the percentages of correlation energies recovered for HF, CO, and N${}_{2}$ improve from 115.2%, 121.7%, and 121.5% without rank restriction to 97.8%, 101.1%, and 100.0% with rank restriction. Similar results are found at both equilibrium and nonequilibrium geometries. While more accurate, the rank-restricted $N$-representability conditions are less expensive computationally than the full-rank conditions.

arXiv:1111.1978•Physical Review A•Nov 8, 2011•14 citations
In methods like geminal-based approaches or coupled cluster that are solved using the projected Schrödinger equation, direct computation of the 2-electron reduced density matrix (2-RDM) is impractical and one falls back to a 2-RDM based on response theory. However, the 2-RDMs from response theory are not N-representable. That is, the response 2-RDM does not correspond to an actual physical N-electron wave function. We present a new algorithm for making these non-N-representable 2-RDMs approximately N-representable, i.e., it has the right symmetry and normalization and it fulfills the P-, Q-, and G-conditions. Next to an algorithm which can be applied to any 2-RDM, we have also developed a 2-RDM optimization procedure specifically for seniority-zero 2-RDMs. We aim to find the 2-RDM with the right properties which is the closest (in the sense of the Frobenius norm) to the non-N-representable 2-RDM by minimizing the square norm of the difference between this initial response 2-RDM and the targeted 2-RDM under the constraint that the trace is normalized and the 2-RDM, Q-matrix, and G-matrix are positive semidefinite, i.e., their eigenvalues are non-negative. Our method is suitable for fixing non-N-representable 2-RDMs which are close to being N-representable. Through the N-representability optimization algorithm we add a small correction to the initial 2-RDM such that it fulfills the most important N-representability conditions.
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.
For 50 years, progress toward the direct calculation of the ground-state two-electron reduced density matrix (2-RDM) was stymied from an inability to constrain the 2-RDM to represent an N-electron wave function. Recent advances in theory and optimization have realized the direct calculation of the 2-RDM. A variational 2-RDM procedure, using first-order semidefinite programming, has been shown to capture multireference correlation effects important at nonequilibrium geometries [Mazziotti, Phys. Rev. Lett. 2004, 93, No. 213001]. This method emerged from research on a nonvariational calculation of the 2-RDM by the contracted Schrödinger equation. Both approaches will be discussed and illustrated.
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.
The variational two-electron reduced-density-matrix (2-RDM) method computes the ground-state energy and 2-RDM of an atom or molecule without calculation of the many-electron wave function. Recently, the computational efficiency of the 2-RDM method has been significantly enhanced through the use of a first-order algorithm for semidefinite programming [Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)]. In this paper we develop a spin- and symmetry-adapted formulation of the method that further improves its efficiency by incorporating both the spin and spatial symmetries of many-electron atoms and molecules. While previous work on density-matrix symmetry focused on only one form of the 2-RDM, the variational method employs three different forms of the 2-RDM, known as the DDD, QQQ, and GGG matrices, to restrict the 2-RDM to be approximately NNN-representable, that is representable by an NNN-electron wave function. We apply spin symmetries to the three forms of the 2-RDM, each of which breaks into four diagonal spin-blocks, namely one singlet and three triplet blocks. If the molecules have point-group symmetry, each of the 2-RDMs may be further subdivided into smaller diagonal blocks according to the spatial symmetry of the basis functions. The subdivision of the 2-RDMs into diagonal blocks generates significant computational savings in both floating-point operations and memory storage. Calculations illustrate the computational savings. Spin adaptation also enforces the correct expectation value of the ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{S}}^{2}$ operator, which in earlier work is applied as a separate constraint.
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.
Variational minimization of the ground-state energy as a function of the two-electron reduced density matrix (2-RDM), constrained by necessary NNN-representability conditions, provides a polynomial-scaling approach to studying strongly correlated molecules without computing the many-electron wave function. Here we introduce a route to enhancing necessary conditions for NNN representability through rank restriction of the 2-RDM. Rather than adding computationally more expensive NNN-representability conditions, we directly enhance the accuracy of two-particle (2-positivity) conditions through rank restriction, which removes degrees of freedom in the 2-RDM that are not sufficiently constrained. We select the rank of the particle-hole 2-RDM by deriving the ranks associated with model wave functions, including both mean-field and antisymmetrized geminal power (AGP) wave functions. Because the 2-positivity conditions are exact for quantum systems with AGP ground states, the rank of the particle-hole 2-RDM from the AGP ansatz provides a minimum for its value in variational 2-RDM calculations of general quantum systems. To implement the rank-restricted conditions, we extend a first-order algorithm for large-scale semidefinite programming. The rank-restricted conditions significantly improve the accuracy of the energies; for example, the percentages of correlation energies recovered for HF, CO, and N2{}_{2}2​ improve from 115.2%, 121.7%, and 121.5% without rank restriction to 97.8%, 101.1%, and 100.0% with rank restriction. Similar results are found at both equilibrium and nonequilibrium geometries. While more accurate, the rank-restricted NNN-representability conditions are less expensive computationally than the full-rank conditions.