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 two-electron reduced density matrix theory for many-electron atoms and molecules: Implementation of the spin- and symmetry-adapted T2 condition through first-order semidefinite programming

David A. Mazziotti

Published Sep 30, 2005• 74 citations
Publisher

Abstract

The energy and properties of a many-electron atom or molecule may be directly computed from a variational optimization of a two-electron reduced density matrix (2RDM) that is constrained to represent many-electron quantum systems. In this paper we implement a variational 2RDM method with a representability constraint, known as the ${T}_{2}$ condition. The optimization of the 2RDM is performed with a first-order algorithm for semidefinite programming [D. A. Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)] which, because of its lower computational cost in comparison to second-order methods, allows the treatment of larger basis sets. We also derive and implement a spin- and symmetry-adapted formulation of the ${T}_{2}$ condition that significantly decreases the size of the largest block in the ${T}_{2}$ matrix. The ${T}_{2}$ condition, originally derived by Erdahl [Int. J. Quantum Chem. 13, 697 (1978)], was recently applied via a second-order algorithm to atoms and molecules [Z. Zhao et al., J. Chem. Phys. 120, 2095 (2004)]. While these calculations were restricted to molecules at equilibrium geometries in minimal basis sets, we apply the 2RDM method with the ${T}_{2}$ condition to compute the electronic energies of molecules in both minimal and nonminimal basis sets at equilibrium as well as nonequilibrium geometries. Accurate potential energies curves are produced for BH, HF, and ${\mathrm{N}}_{2}$. Results are compared with the 2RDM method without the ${T}_{2}$ condition as well as several wave-function methods.

DOI
10.1103/physreva.72.032510

Related papers

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

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

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 applied to the electronic structure of few-electron quantum dots

Adam E. Rothman, David A. Mazziotti

Variational two-electron reduced-density-matrix (2RDM) theory is applied to computing energy spectra and properties of few-electron quantum dots. The model system is a two-dimensional electron gas with a central confinement potential. For each orbital angular momentum $J$, the energy and 2RDM are computed by the variational 2RDM method in which the energy is minimized as a functional of the 2RDM. In the minimization, which is performed by semidefinite programming, the 2RDM is constrained to represent a $N$-electron wave function with angular momentum $J$ by $N$- and $J$-representability conditions [D. A. Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)]. Advantages of the variational 2RDM method include (i) lower bounds on the energies for all $J$ values, (ii) calculation of approximate 2RDMs in polynomial time without many-electron wave functions, (iii) exploitation of angular symmetry in the sparse block-diagonal structure of the 2RDM, (iv) accurate description of multireference correlation (entanglement) effects, and (v) direct calculation of one- and two-electron properties from the 2RDM. With the 2RDM we directly compute pair-correlation functions, radial charge densities, and average radial electron displacements in the quantum dot. Energies and properties are compared to those from solving the $N$-electron Schr\"odinger equation by large-scale, exact diagonalization. It is found that the accuracy of the variational 2RDM approach is sensitive to the total orbital angular momentum and the symmetry of the final wave function. For quantum dots of high symmetry, the variational algorithm isolates a highly accurate solution that recovers the correlation energy within a few percent.

Physical Review A•Sep 9, 2008•27 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
The energy and properties of a many-electron atom or molecule may be directly computed from a variational optimization of a two-electron reduced density matrix (2RDM) that is constrained to represent many-electron quantum systems. In this paper we implement a variational 2RDM method with a representability constraint, known as the T2{T}_{2}T2​ condition. The optimization of the 2RDM is performed with a first-order algorithm for semidefinite programming [D. A. Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)] which, because of its lower computational cost in comparison to second-order methods, allows the treatment of larger basis sets. We also derive and implement a spin- and symmetry-adapted formulation of the T2{T}_{2}T2​ condition that significantly decreases the size of the largest block in the T2{T}_{2}T2​ matrix. The T2{T}_{2}T2​ condition, originally derived by Erdahl [Int. J. Quantum Chem. 13, 697 (1978)], was recently applied via a second-order algorithm to atoms and molecules [Z. Zhao et al., J. Chem. Phys. 120, 2095 (2004)]. While these calculations were restricted to molecules at equilibrium geometries in minimal basis sets, we apply the 2RDM method with the T2{T}_{2}T2​ condition to compute the electronic energies of molecules in both minimal and nonminimal basis sets at equilibrium as well as nonequilibrium geometries. Accurate potential energies curves are produced for BH, HF, and N2{\mathrm{N}}_{2}N2​. Results are compared with the 2RDM method without the T2{T}_{2}T2​ condition as well as several wave-function methods.
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.
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.
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.
Variational two-electron reduced-density-matrix (2RDM) theory is applied to computing energy spectra and properties of few-electron quantum dots. The model system is a two-dimensional electron gas with a central confinement potential. For each orbital angular momentum JJJ, the energy and 2RDM are computed by the variational 2RDM method in which the energy is minimized as a functional of the 2RDM. In the minimization, which is performed by semidefinite programming, the 2RDM is constrained to represent a NNN-electron wave function with angular momentum JJJ by NNN- and JJJ-representability conditions [D. A. Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)]. Advantages of the variational 2RDM method include (i) lower bounds on the energies for all JJJ values, (ii) calculation of approximate 2RDMs in polynomial time without many-electron wave functions, (iii) exploitation of angular symmetry in the sparse block-diagonal structure of the 2RDM, (iv) accurate description of multireference correlation (entanglement) effects, and (v) direct calculation of one- and two-electron properties from the 2RDM. With the 2RDM we directly compute pair-correlation functions, radial charge densities, and average radial electron displacements in the quantum dot. Energies and properties are compared to those from solving the NNN-electron Schr\"odinger equation by large-scale, exact diagonalization. It is found that the accuracy of the variational 2RDM approach is sensitive to the total orbital angular momentum and the symmetry of the final wave function. For quantum dots of high symmetry, the variational algorithm isolates a highly accurate solution that recovers the correlation energy within a few percent.
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.