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

Physical Review A

Variational reduced-density-matrix calculation of the one-dimensional Hubbard model

Jeff R. Hammond, David A. Mazziotti

Published Jun 13, 2006• 41 citations
Publisher

Abstract

Variational reduced-density-matrix theory is applied to calculating the ground-state energy and two-electron reduced density matrices (2-RDMs) of the one-dimensional Hubbard model for a range of interaction strengths. The 2-RDM is constrained to represent an $N$-particle wave function by two sets of $N$-representability conditions, known as the 2- and partial 3-positivity conditions. Variational optimization of the energy with the 2-RDM constrained by $N$-representability conditions is performed using a first-order semidefinite-programming algorithm that was developed for treating atoms and molecules [D. A. Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)]. Accurate energies for a broad range of interaction strengths indicate that the variational 2-RDM method is a valuable tool for studying strongly correlated electrons.

DOI
10.1103/physreva.73.062505

Related papers

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

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

Variational calculation of second-order reduced density matrices by strong N-representability conditions and an accurate semidefinite programming solver

Maho Nakata, Bastiaan J. Braams, Katsuki Fujisawa +4 more

The reduced density matrix (RDM) method, which is a variational calculation based on the second-order reduced density matrix, is applied to the ground state energies and the dipole moments for 57 different states of atoms, molecules, and to the ground state energies and the elements of 2-RDM for the Hubbard model. We explore the well-known N-representability conditions (P, Q, and G) together with the more recent and much stronger T1 and T2(') conditions. T2(') condition was recently rederived and it implies T2 condition. Using these N-representability conditions, we can usually calculate correlation energies in percentage ranging from 100% to 101%, whose accuracy is similar to CCSD(T) and even better for high spin states or anion systems where CCSD(T) fails. Highly accurate calculations are carried out by handling equality constraints and/or developing multiple precision arithmetic in the semidefinite programming (SDP) solver. Results show that handling equality constraints correctly improves the accuracy from 0.1 to 0.6 mhartree. Additionally, improvements by replacing T2 condition with T2(') condition are typically of 0.1-0.5 mhartree. The newly developed multiple precision arithmetic version of SDP solver calculates extraordinary accurate energies for the one dimensional Hubbard model and Be atom. It gives at least 16 significant digits for energies, where double precision calculations gives only two to eight digits. It also provides physically meaningful results for the Hubbard model in the high correlation limit.

The Journal of Chemical Physics•Apr 28, 2008•75 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

First-order semidefinite programming for the direct determination of two-electron reduced density matrices with application to many-electron atoms and molecules

David A. Mazziotti

Direct variational calculation of two-electron reduced density matrices (2-RDMs) for many-electron atoms and molecules in nonminimal basis sets has recently been achieved through the use of first-order semidefinite programming [D. A. Mazziotti, Phys. Rev. Lett. (in press)]. With semidefinite programming, the electronic ground-state energy of a molecule is minimized with respect to the 2-RDM subject to N-representability constraints known as positivity conditions. Here we present a detailed account of the first-order algorithm for semidefinite programming and its comparison with the primal-dual interior-point algorithms employed in earlier variational 2-RDM calculations. The first-order semidefinite-programming algorithm, computations show, offers an orders-of-magnitude reduction in floating-point operations and storage in comparison with previous implementations. We also examine the ability of the positivity conditions to treat strong correlation and multireference effects through an analysis of the Hamiltonians for which the conditions are exact. Calculations are performed in nonminimal basis sets for a variety of atoms and molecules and the potential-energy curves for CO and H(2)O.

The Journal of Chemical Physics•Nov 30, 2004•87 citations

First-order semidefinite programming for the two-electron treatment of many-electron atoms and molecules

David A. Mazziotti

The ground-state energy and properties of any many-electron atom or molecule may be rigorously computed by variationally computing the two-electron reduced density matrix rather than the many-electron wavefunction. While early attempts fifty years ago to compute the ground-state 2-RDM directly were stymied because the 2-RDM must be constrained to represent an N-electron wavefunction, recent advances in theory and optimization have made direct computation of the 2-RDM possible. The constraints in the variational calculation of the 2-RDM require a special optimization known as a semidefinite programming. Development of first-order semidefinite programming for the 2-RDM method has reduced the computational costs of the calculation by orders of magnitude [Mazziotti, Phys. Rev. Lett. 93 (2004) 213001]. The variational 2-RDM approach is effective at capturing multi-reference correlation effects that are especially important at non-equilibrium molecular geometries. Recent work on 2-RDM methods will be reviewed and illustrated with particular emphasis on the importance of advances in large-scale semidefinite programming.

ESAIM Mathematical Modelling and Numerical Analysis•Mar 1, 2007•35 citations
Variational reduced-density-matrix theory is applied to calculating the ground-state energy and two-electron reduced density matrices (2-RDMs) of the one-dimensional Hubbard model for a range of interaction strengths. The 2-RDM is constrained to represent an NNN-particle wave function by two sets of NNN-representability conditions, known as the 2- and partial 3-positivity conditions. Variational optimization of the energy with the 2-RDM constrained by NNN-representability conditions is performed using a first-order semidefinite-programming algorithm that was developed for treating atoms and molecules [D. A. Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)]. Accurate energies for a broad range of interaction strengths indicate that the variational 2-RDM method is a valuable tool for studying strongly correlated electrons.
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.
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.
The reduced density matrix (RDM) method, which is a variational calculation based on the second-order reduced density matrix, is applied to the ground state energies and the dipole moments for 57 different states of atoms, molecules, and to the ground state energies and the elements of 2-RDM for the Hubbard model. We explore the well-known N-representability conditions (P, Q, and G) together with the more recent and much stronger T1 and T2(') conditions. T2(') condition was recently rederived and it implies T2 condition. Using these N-representability conditions, we can usually calculate correlation energies in percentage ranging from 100% to 101%, whose accuracy is similar to CCSD(T) and even better for high spin states or anion systems where CCSD(T) fails. Highly accurate calculations are carried out by handling equality constraints and/or developing multiple precision arithmetic in the semidefinite programming (SDP) solver. Results show that handling equality constraints correctly improves the accuracy from 0.1 to 0.6 mhartree. Additionally, improvements by replacing T2 condition with T2(') condition are typically of 0.1-0.5 mhartree. The newly developed multiple precision arithmetic version of SDP solver calculates extraordinary accurate energies for the one dimensional Hubbard model and Be atom. It gives at least 16 significant digits for energies, where double precision calculations gives only two to eight digits. It also provides physically meaningful results for the Hubbard model in the high correlation limit.
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.
Direct variational calculation of two-electron reduced density matrices (2-RDMs) for many-electron atoms and molecules in nonminimal basis sets has recently been achieved through the use of first-order semidefinite programming [D. A. Mazziotti, Phys. Rev. Lett. (in press)]. With semidefinite programming, the electronic ground-state energy of a molecule is minimized with respect to the 2-RDM subject to N-representability constraints known as positivity conditions. Here we present a detailed account of the first-order algorithm for semidefinite programming and its comparison with the primal-dual interior-point algorithms employed in earlier variational 2-RDM calculations. The first-order semidefinite-programming algorithm, computations show, offers an orders-of-magnitude reduction in floating-point operations and storage in comparison with previous implementations. We also examine the ability of the positivity conditions to treat strong correlation and multireference effects through an analysis of the Hamiltonians for which the conditions are exact. Calculations are performed in nonminimal basis sets for a variety of atoms and molecules and the potential-energy curves for CO and H(2)O.
The ground-state energy and properties of any many-electron atom or molecule may be rigorously computed by variationally computing the two-electron reduced density matrix rather than the many-electron wavefunction. While early attempts fifty years ago to compute the ground-state 2-RDM directly were stymied because the 2-RDM must be constrained to represent an N-electron wavefunction, recent advances in theory and optimization have made direct computation of the 2-RDM possible. The constraints in the variational calculation of the 2-RDM require a special optimization known as a semidefinite programming. Development of first-order semidefinite programming for the 2-RDM method has reduced the computational costs of the calculation by orders of magnitude [Mazziotti, Phys. Rev. Lett. 93 (2004) 213001]. The variational 2-RDM approach is effective at capturing multi-reference correlation effects that are especially important at non-equilibrium molecular geometries. Recent work on 2-RDM methods will be reviewed and illustrated with particular emphasis on the importance of advances in large-scale semidefinite programming.