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

Journal of Statistical Mechanics Theory and Experiment

Variational determination of the two-particle reduced density matrix within the doubly occupied configuration interaction space: exploiting translational and reflection invariance

Gustavo E. Massaccesi, A Rubio-García, P. Capuzzi, Elías Ríos, Ofelia B. Oña, J. Dukelsky, Luis Laín, Alicia Torre, Diego R. Alcoba

Published Jan 1, 2021• 14 citations
Publisher

Abstract

Abstract This work incorporates translational and reflection symmetry reductions to the variational determination of the two-particle reduced density matrix (2-RDM) corresponding to the ground state of N -particle systems, within the doubly occupied configuration interaction (DOCI) space. By exploiting these symmetries within this lower-bound variational methodology it is possible to treat larger systems than those previously studied. The 2-RDM matrix elements are calculated by imposing up to four-particle N -representability constraint conditions using standard semidefinite programing algorithms. The method is applied to the one- and two-dimensional XXZ spin 1/2 model of quantum magnetism. Several observables including the energy and the spin–spin correlation functions are obtained to assess the physical content of the variationally determined 2-RDM. Comparison with quantum-Monte Carlo and matrix product state simulations shows that in most cases only requiring up to three-particle positivity conditions is enough to correctly describe the ground-state properties of these one- and two-dimensional models.

DOI
10.1088/1742-5468/abd940

Related papers

Variational reduced density matrix method in the doubly-occupied configuration interaction space using four-particle N-representability conditions: Application to the XXZ model of quantum magnetism

Álvaro Rubio-García, J. Dukelsky, Diego R. Alcoba +5 more

This work deals with the variational determination of the two-particle reduced density matrix (2-RDM) and the energy corresponding to the ground state of N-particle systems within the doubly occupied configuration interaction (DOCI) space. Here, we impose for the first time up to four-particle N-representability constraint conditions in the variational determination of the 2-RDM matrix elements using the standard semidefinite programming algorithms. The energies and 2-RDMs obtained from this treatment and the corresponding computational costs are compared with those arisen from previously reported less restrictive variational methods [D. R. Alcoba et al., J. Chem. Phys. 149, 194105 (2018)] as well as with the exact DOCI values. We apply the different approximations to the one-dimensional XXZ model of quantum magnetism, which has a rich phase diagram with one critical phase and constitutes a stringent test for the method. The numerical results show the usefulness of our treatment to achieve a high degree of accuracy.

The Journal of Chemical Physics•Oct 15, 2019•20 citations

Variational reduced density matrix method in the doubly occupied configuration interaction space using three-particle N-representability conditions

Diego R. Alcoba, P. Capuzzi, Álvaro Rubio-García +5 more

Ground-state energies and two-particle reduced density matrices (2-RDMs) corresponding to N-particle systems are computed variationally within the doubly occupied configuration interaction (DOCI) space by constraining the 2-RDM to satisfy a complete set of three-particle N-representability conditions known as three-positivity conditions. These conditions are derived and implemented in the variational calculation of the 2-RDM with standard semidefinite programming algorithms. Ground state energies and 2-RDMs are computed for N2, CO, CN−, and NO+ molecules at both equilibrium and nonequilibrium geometries as well as for pairing models at different repulsive interaction strengths. The results from the full three-positivity conditions are compared with those from the exact DOCI method and with approximated 2-RDM variational ones obtained within two-positivity and two-positivity plus a subset of three-positivity conditions, as recently reported [D. R. Alcoba et al., J. Chem. Phys. 148, 024105 (2018) and A. Rubio-García et al., J. Chem. Theory Comput. 14, 4183 (2018)]. The accuracy of these numerical determinations and their low computational cost demonstrate the usefulness of the three-particle variational constraints within the DOCI framework.

The Journal of Chemical Physics•Nov 20, 2018•20 citations

Direct variational determination of the two-electron reduced density matrix for doubly occupied-configuration-interaction wave functions: The influence of three-index N-representability conditions

Diego R. Alcoba, Alicia Torre, Luis Laín +7 more

This work proposes the variational determination of two-electron reduced density matrices corresponding to the ground state of N-electron systems within the doubly occupied-configuration-interaction methodology. The P, Q, and G two-index N-representability conditions have been extended to the T1 and T2 (T2') three-index ones and the resulting optimization problem has been addressed using a standard semidefinite program. We report results obtained from the doubly occupied-configuration-interaction method, from the two-index constraint variational procedure and from the two- and three-index constraint variational treatment. The discussion of these results along with a study of the computational cost demanded shows the usefulness of our proposal.

The Journal of Chemical Physics•Jan 9, 2018•37 citations

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

Jeff R. Hammond, David A. Mazziotti

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.

Physical Review A•Jun 13, 2006•41 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

Variational determination of ground and excited-state two-electron reduced density matrices in the doubly occupied configuration space: A dispersion operator approach

Diego R. Alcoba, Ofelia B. Oña, Luis Laín +6 more

This work implements a variational determination of the elements of two-electron reduced density matrices corresponding to the ground and excited states of N-electron interacting systems based on the dispersion operator technique. The procedure extends the previously reported proposal [Nakata et al., J. Chem. Phys. 125, 244109 (2006)] to two-particle interaction Hamiltonians and N-representability conditions for the two-, three-, and four-particle reduced density matrices in the doubly occupied configuration interaction space. The treatment has been applied to describe electronic spectra using two benchmark exactly solvable pairing models: reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev Hamiltonians. The dispersion operator combined with N-representability conditions up to the four-particle reduced density matrices provides excellent results.

The Journal of Chemical Physics•Jun 10, 2021•12 citations
Abstract This work incorporates translational and reflection symmetry reductions to the variational determination of the two-particle reduced density matrix (2-RDM) corresponding to the ground state of N -particle systems, within the doubly occupied configuration interaction (DOCI) space. By exploiting these symmetries within this lower-bound variational methodology it is possible to treat larger systems than those previously studied. The 2-RDM matrix elements are calculated by imposing up to four-particle N -representability constraint conditions using standard semidefinite programing algorithms. The method is applied to the one- and two-dimensional XXZ spin 1/2 model of quantum magnetism. Several observables including the energy and the spin–spin correlation functions are obtained to assess the physical content of the variationally determined 2-RDM. Comparison with quantum-Monte Carlo and matrix product state simulations shows that in most cases only requiring up to three-particle positivity conditions is enough to correctly describe the ground-state properties of these one- and two-dimensional models.
This work deals with the variational determination of the two-particle reduced density matrix (2-RDM) and the energy corresponding to the ground state of N-particle systems within the doubly occupied configuration interaction (DOCI) space. Here, we impose for the first time up to four-particle N-representability constraint conditions in the variational determination of the 2-RDM matrix elements using the standard semidefinite programming algorithms. The energies and 2-RDMs obtained from this treatment and the corresponding computational costs are compared with those arisen from previously reported less restrictive variational methods [D. R. Alcoba et al., J. Chem. Phys. 149, 194105 (2018)] as well as with the exact DOCI values. We apply the different approximations to the one-dimensional XXZ model of quantum magnetism, which has a rich phase diagram with one critical phase and constitutes a stringent test for the method. The numerical results show the usefulness of our treatment to achieve a high degree of accuracy.
Ground-state energies and two-particle reduced density matrices (2-RDMs) corresponding to N-particle systems are computed variationally within the doubly occupied configuration interaction (DOCI) space by constraining the 2-RDM to satisfy a complete set of three-particle N-representability conditions known as three-positivity conditions. These conditions are derived and implemented in the variational calculation of the 2-RDM with standard semidefinite programming algorithms. Ground state energies and 2-RDMs are computed for N2, CO, CN−, and NO+ molecules at both equilibrium and nonequilibrium geometries as well as for pairing models at different repulsive interaction strengths. The results from the full three-positivity conditions are compared with those from the exact DOCI method and with approximated 2-RDM variational ones obtained within two-positivity and two-positivity plus a subset of three-positivity conditions, as recently reported [D. R. Alcoba et al., J. Chem. Phys. 148, 024105 (2018) and A. Rubio-García et al., J. Chem. Theory Comput. 14, 4183 (2018)]. The accuracy of these numerical determinations and their low computational cost demonstrate the usefulness of the three-particle variational constraints within the DOCI framework.
This work proposes the variational determination of two-electron reduced density matrices corresponding to the ground state of N-electron systems within the doubly occupied-configuration-interaction methodology. The P, Q, and G two-index N-representability conditions have been extended to the T1 and T2 (T2') three-index ones and the resulting optimization problem has been addressed using a standard semidefinite program. We report results obtained from the doubly occupied-configuration-interaction method, from the two-index constraint variational procedure and from the two- and three-index constraint variational treatment. The discussion of these results along with a study of the computational cost demanded shows the usefulness of our proposal.
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.
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.
This work implements a variational determination of the elements of two-electron reduced density matrices corresponding to the ground and excited states of N-electron interacting systems based on the dispersion operator technique. The procedure extends the previously reported proposal [Nakata et al., J. Chem. Phys. 125, 244109 (2006)] to two-particle interaction Hamiltonians and N-representability conditions for the two-, three-, and four-particle reduced density matrices in the doubly occupied configuration interaction space. The treatment has been applied to describe electronic spectra using two benchmark exactly solvable pairing models: reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev Hamiltonians. The dispersion operator combined with N-representability conditions up to the four-particle reduced density matrices provides excellent results.