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

The Journal of Physical Chemistry C

Strong Electron Correlation in Materials from Pair-Interacting Model Hamiltonians

Nicholas C. Rubin, David A. Mazziotti

Published Mar 26, 2015• 10 citations
Publisher

Abstract

Strong electron correlation in materials is explored within a class of model Hamiltonians that treat only pair interactions between electrons. The model is unique among typical spin Hamiltonians in that it does not have an effective mean-field reference wave function. The ground-state wave functions from all Hamiltonians in the model have the same one-electron reduced density matrix (1-RDM); consequently, one-electron theories such as the Hartree–Fock and density functional theories are inapplicable. In contrast, the ground-state two-electron reduced density matrix (2-RDM) has a one-to-one mapping to the ground-state wave function. For a range of lattices including linear, ladder, and square topologies, we variationally compute the 2-RDM subject to constraints, known as N -representability conditions, that are necessary for the 2-RDM to represent an N -electron ensemble density matrix. We find that for all model Hamiltonians the 2-RDM is accurately computed as long as the D, Q, and G N -representability conditions are supplemented with the T 2 condition. Energies, orbital and geminal occupation numbers, and correlation functions are computed. The model has applications to superconductivity as well as more general pairing phenomena in electronic systems. Effective methods are needed to treat strong electron correlation arising in the study of materials including the study of high-temperature superconductivity and phase transitions.

DOI
10.1021/jp5130266

Related papers

THE VARIATIONAL TWO-ELECTRON REDUCED-DENSITY-MATRIX METHOD FOR EXTENDED SYSTEMS

Nicholas C. Rubin

In this thesis we develop the variational two-electron reduced-density-matrix method for extended systems. Extended systems are represented in two ways: i) lattice models describing the dominant valence electronic structure with periodic boundaries to account for their extended nature and ii) a crystalline-orbital basis built from atomic orbitals using the generalization of molecular orbital theory to polymers. The first part of this thesis (Ch. 3-4) examines the performance of the variational 2-RDM method on lattice systems with tunable electron correlation. The first of these systems is the classic Hubbard model with linear and ladder lattice topologies. Because electron correlation functions, such as charge- and spin-ordering, are linear functions of the 2-RDM, the difference in electronic structure between one- and quasi-one-dimensional systems is accurately characterized. The second model contains only two-body interactions and is unique among typical spin models in that it does not have a mean-field reference wave function. The ground state wave functions from all Hamiltonians in the model have the same 1-electron reduced density matrix; consequently, one-electron theories are largely inapplicable. The superconducting eta-pairing ground states make the model a unique tool for demonstrating the necessary N-representability in highly correlated environments. The second part of this thesis (Ch. 5-6) develops a formalism for modeling materials by solving the full Schrodinger equation. Crystalline-orbital Hartree-Fock provides a set of orbitals and integral tensors for the variational 2-RDM method. We demonstrate that time-reversal symmetry, which is implicitly included in position space electronic structure calculations, must be explicitly included as an N-representability constraint on the 2-RDM when using a momentum space basis. The necessity of these equality constraints is demonstrated by the accurate recovery of the binding energy of two polymers and the symmetry of their natural orbital occupations.

•Jan 1, 2016•0 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

Two-electron reduced density matrix method for the electronic structure of correlated materials

Anna O. Schouten, Simon Ewing, David A. Mazziotti

The last several decades have seen significant advances in the theoretical modeling of materials within the fields of solid-state physics and materials science, but many methods commonly applied to this problem struggle to capture strong electron correlation accurately. Recent widespread interest in quantum materials---where strong correlation plays a crucial role in the quantum effects governing their behavior---further highlights the need for theoretical methods capable of rigorously treating such correlation. Here, we present a periodic generalization of variational two-electron reduced density matrix (2-RDM) theory, a bootstrapping-type method that minimizes the ground-state energy as a functional of the 2-RDM without relying on the wave function. The 2-RDM is computed directly by semidefinite programming with $N$-representability conditions, ensuring accurate treatment of strongly correlated electronic systems. By exploiting translational symmetry, we significantly reduce computational scaling, enabling applications to realistic materials-scale systems. Additionally, we introduce an alternative to conventional energy band structures: natural-orbital occupation-number bands, which, being independent of mean-field assumptions, offer deeper insights into electron correlation effects. We demonstrate the effectiveness of this approach by applying the theory to hydrogen chains, molybdenum disulfide, and nickel oxide, showing that natural-orbital occupation bands correctly capture electronic character in regimes where density functional theory fails. This work represents a major step toward accurately describing the electronic structure of quantum materials using reduced density matrices rather than wave functions.

Physical review. B./Physical review. B•Nov 12, 2025•5 citations

Bootstrapping the Electronic Structure of Quantum Materials

Anna O. Schouten, Simon Ewing, David A. Mazziotti

The last several decades have seen significant advances in the theoretical modeling of materials within the fields of solid-state physics and materials science, but many methods commonly applied to this problem struggle to capture strong electron correlation accurately. Recent widespread interest in quantum materials -- where strong correlation plays a crucial role in the quantum effects governing their behavior -- further highlights the need for theoretical methods capable of rigorously treating such correlation. Here, we present a periodic generalization of variational two-electron reduced density matrix (2-RDM) theory, a bootstrapping-type method that minimizes the ground-state energy as a functional of the 2-RDM without relying on the wavefunction. The 2-RDM is computed directly by semidefinite programming with $N$-representability conditions, ensuring accurate treatment of strongly correlated electronic systems. By exploiting translational symmetry, we significantly reduce computational scaling, enabling applications to realistic materials-scale systems. Additionally, we introduce an alternative to conventional energy band structures: natural-orbital occupation-number bands, which, being independent of mean-field assumptions, offer deeper insights into electron correlation effects. We demonstrate the effectiveness of this approach by applying the theory to hydrogen chains, molybdenum disulfide, and nickel oxide, showing that natural-orbital occupation bands correctly capture electronic character in regimes where density functional theory fails. This work represents a major step toward accurately describing the electronic structure of quantum materials using reduced density matrices rather than wavefunctions.

arXiv:2504.02861•Mar 31, 2025•1 citation

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

Strong Correlation in Molecular Periodic Systems from a Variational Reduced Density Matrix Theory

Simon Ewing, David A. Mazziotti

Molecular periodic systems have the challenge of not only representing an infinite, albeit periodically repeating, molecule but also treating the strong electron correlation that emerges in such systems. Many post-Hartree-Fock methods that are efficient in moderately correlated molecules often fail to describe strongly correlated electrons and their properties in periodic systems. Here we develop and apply an electronic structure approach to periodic systems based on the variational calculation of the two-electron reduced density matrix (2-RDM). The variational 2-RDM method calculates a lower bound to the ground-state energy of a molecular periodic system without computation or storage of its many-electron wave function. The 2-RDM is constrained by $N$-representability conditions that are necessary for it to represent at least one $N$-electron density matrix---conditions that allow it to treat strongly correlated systems. Two periodic systems are treated to demonstrate the methodology: (i) a metallic hydrogen chain and (ii) a strongly correlated acene chain. In the hydrogen chains the 2-RDM theory describes the strongly correlated Mott metal-to-insulator transition while in acene chains it describes the emergence of polyradical character with increasing chain length. The methodology is applicable to treating strong electron correlation in more general periodic molecules and materials.

•Aug 6, 2020•2 citations
Strong electron correlation in materials is explored within a class of model Hamiltonians that treat only pair interactions between electrons. The model is unique among typical spin Hamiltonians in that it does not have an effective mean-field reference wave function. The ground-state wave functions from all Hamiltonians in the model have the same one-electron reduced density matrix (1-RDM); consequently, one-electron theories such as the Hartree–Fock and density functional theories are inapplicable. In contrast, the ground-state two-electron reduced density matrix (2-RDM) has a one-to-one mapping to the ground-state wave function. For a range of lattices including linear, ladder, and square topologies, we variationally compute the 2-RDM subject to constraints, known as N -representability conditions, that are necessary for the 2-RDM to represent an N -electron ensemble density matrix. We find that for all model Hamiltonians the 2-RDM is accurately computed as long as the D, Q, and G N -representability conditions are supplemented with the T 2 condition. Energies, orbital and geminal occupation numbers, and correlation functions are computed. The model has applications to superconductivity as well as more general pairing phenomena in electronic systems. Effective methods are needed to treat strong electron correlation arising in the study of materials including the study of high-temperature superconductivity and phase transitions.
In this thesis we develop the variational two-electron reduced-density-matrix method for extended systems. Extended systems are represented in two ways: i) lattice models describing the dominant valence electronic structure with periodic boundaries to account for their extended nature and ii) a crystalline-orbital basis built from atomic orbitals using the generalization of molecular orbital theory to polymers. The first part of this thesis (Ch. 3-4) examines the performance of the variational 2-RDM method on lattice systems with tunable electron correlation. The first of these systems is the classic Hubbard model with linear and ladder lattice topologies. Because electron correlation functions, such as charge- and spin-ordering, are linear functions of the 2-RDM, the difference in electronic structure between one- and quasi-one-dimensional systems is accurately characterized. The second model contains only two-body interactions and is unique among typical spin models in that it does not have a mean-field reference wave function. The ground state wave functions from all Hamiltonians in the model have the same 1-electron reduced density matrix; consequently, one-electron theories are largely inapplicable. The superconducting eta-pairing ground states make the model a unique tool for demonstrating the necessary N-representability in highly correlated environments. The second part of this thesis (Ch. 5-6) develops a formalism for modeling materials by solving the full Schrodinger equation. Crystalline-orbital Hartree-Fock provides a set of orbitals and integral tensors for the variational 2-RDM method. We demonstrate that time-reversal symmetry, which is implicitly included in position space electronic structure calculations, must be explicitly included as an N-representability constraint on the 2-RDM when using a momentum space basis. The necessity of these equality constraints is demonstrated by the accurate recovery of the binding energy of two polymers and the symmetry of their natural orbital occupations.
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 last several decades have seen significant advances in the theoretical modeling of materials within the fields of solid-state physics and materials science, but many methods commonly applied to this problem struggle to capture strong electron correlation accurately. Recent widespread interest in quantum materials---where strong correlation plays a crucial role in the quantum effects governing their behavior---further highlights the need for theoretical methods capable of rigorously treating such correlation. Here, we present a periodic generalization of variational two-electron reduced density matrix (2-RDM) theory, a bootstrapping-type method that minimizes the ground-state energy as a functional of the 2-RDM without relying on the wave function. The 2-RDM is computed directly by semidefinite programming with NNN-representability conditions, ensuring accurate treatment of strongly correlated electronic systems. By exploiting translational symmetry, we significantly reduce computational scaling, enabling applications to realistic materials-scale systems. Additionally, we introduce an alternative to conventional energy band structures: natural-orbital occupation-number bands, which, being independent of mean-field assumptions, offer deeper insights into electron correlation effects. We demonstrate the effectiveness of this approach by applying the theory to hydrogen chains, molybdenum disulfide, and nickel oxide, showing that natural-orbital occupation bands correctly capture electronic character in regimes where density functional theory fails. This work represents a major step toward accurately describing the electronic structure of quantum materials using reduced density matrices rather than wave functions.
The last several decades have seen significant advances in the theoretical modeling of materials within the fields of solid-state physics and materials science, but many methods commonly applied to this problem struggle to capture strong electron correlation accurately. Recent widespread interest in quantum materials -- where strong correlation plays a crucial role in the quantum effects governing their behavior -- further highlights the need for theoretical methods capable of rigorously treating such correlation. Here, we present a periodic generalization of variational two-electron reduced density matrix (2-RDM) theory, a bootstrapping-type method that minimizes the ground-state energy as a functional of the 2-RDM without relying on the wavefunction. The 2-RDM is computed directly by semidefinite programming with NNN-representability conditions, ensuring accurate treatment of strongly correlated electronic systems. By exploiting translational symmetry, we significantly reduce computational scaling, enabling applications to realistic materials-scale systems. Additionally, we introduce an alternative to conventional energy band structures: natural-orbital occupation-number bands, which, being independent of mean-field assumptions, offer deeper insights into electron correlation effects. We demonstrate the effectiveness of this approach by applying the theory to hydrogen chains, molybdenum disulfide, and nickel oxide, showing that natural-orbital occupation bands correctly capture electronic character in regimes where density functional theory fails. This work represents a major step toward accurately describing the electronic structure of quantum materials using reduced density matrices rather than wavefunctions.
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.
Molecular periodic systems have the challenge of not only representing an infinite, albeit periodically repeating, molecule but also treating the strong electron correlation that emerges in such systems. Many post-Hartree-Fock methods that are efficient in moderately correlated molecules often fail to describe strongly correlated electrons and their properties in periodic systems. Here we develop and apply an electronic structure approach to periodic systems based on the variational calculation of the two-electron reduced density matrix (2-RDM). The variational 2-RDM method calculates a lower bound to the ground-state energy of a molecular periodic system without computation or storage of its many-electron wave function. The 2-RDM is constrained by NNN-representability conditions that are necessary for it to represent at least one NNN-electron density matrix---conditions that allow it to treat strongly correlated systems. Two periodic systems are treated to demonstrate the methodology: (i) a metallic hydrogen chain and (ii) a strongly correlated acene chain. In the hydrogen chains the 2-RDM theory describes the strongly correlated Mott metal-to-insulator transition while in acene chains it describes the emergence of polyradical character with increasing chain length. The methodology is applicable to treating strong electron correlation in more general periodic molecules and materials.