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

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

Simon Ewing, David A. Mazziotti

Published Aug 6, 2020• 2 citations
OpenAlex

Abstract

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.

Related papers

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

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

Strong Electron Correlation in Materials from Pair-Interacting Model Hamiltonians

Nicholas C. Rubin, David A. Mazziotti

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.

The Journal of Physical Chemistry C•Mar 26, 2015•10 citations

Strong Correlation in Acene Sheets from the Active-Space Variational Two-Electron Reduced Density Matrix Method: Effects of Symmetry and Size

Kenley M. Pelzer, Loren Greenman, Gergely Gidofalvi +1 more

Polyaromatic hydrocarbons (PAHs) are a class of organic molecules with importance in several branches of science, including medicine, combustion chemistry, and materials science. The delocalized π-orbital systems in PAHs require highly accurate electronic structure methods to capture strong electron correlation. Treating correlation in PAHs has been challenging because (i) traditional wave function methods for strong correlation have not been applicable since they scale exponentially in the number of strongly correlated orbitals, and (ii) alternative methods such as the density-matrix renormalization group and variational two-electron reduced density matrix (2-RDM) methods have not been applied beyond linear acene chains. In this paper we extend the earlier results from active-space variational 2-RDM theory [Gidofalvi, G.; Mazziotti, D. A. J. Chem. Phys. 2008, 129, 134108] to the more general two-dimensional arrangement of rings--acene sheets--to study the relationship between geometry and electron correlation in PAHs. The acene-sheet calculations, if performed with conventional wave function methods, would require wave function expansions with as many as 1.5 × 10(17) configuration state functions. To measure electron correlation, we employ several RDM-based metrics: (i) natural-orbital occupation numbers, (ii) the 1-RDM von Neumann entropy, (iii) the correlation energy per carbon atom, and (iv) the squared Frobenius norm of the cumulant 2-RDM. The results confirm a trend of increasing polyradical character with increasing molecular size previously observed in linear PAHs and reveal a corresponding trend in two-dimensional (arch-shaped) PAHs. Furthermore, in PAHs of similar size they show significant variations in correlation with geometry. PAHs with the strictly linear geometry (chains) exhibit more electron correlation than PAHs with nonlinear geometries (sheets).

The Journal of Physical Chemistry A•May 12, 2011•98 citations
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.
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.
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.
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.
Polyaromatic hydrocarbons (PAHs) are a class of organic molecules with importance in several branches of science, including medicine, combustion chemistry, and materials science. The delocalized π-orbital systems in PAHs require highly accurate electronic structure methods to capture strong electron correlation. Treating correlation in PAHs has been challenging because (i) traditional wave function methods for strong correlation have not been applicable since they scale exponentially in the number of strongly correlated orbitals, and (ii) alternative methods such as the density-matrix renormalization group and variational two-electron reduced density matrix (2-RDM) methods have not been applied beyond linear acene chains. In this paper we extend the earlier results from active-space variational 2-RDM theory [Gidofalvi, G.; Mazziotti, D. A. J. Chem. Phys. 2008, 129, 134108] to the more general two-dimensional arrangement of rings--acene sheets--to study the relationship between geometry and electron correlation in PAHs. The acene-sheet calculations, if performed with conventional wave function methods, would require wave function expansions with as many as 1.5 × 10(17) configuration state functions. To measure electron correlation, we employ several RDM-based metrics: (i) natural-orbital occupation numbers, (ii) the 1-RDM von Neumann entropy, (iii) the correlation energy per carbon atom, and (iv) the squared Frobenius norm of the cumulant 2-RDM. The results confirm a trend of increasing polyradical character with increasing molecular size previously observed in linear PAHs and reveal a corresponding trend in two-dimensional (arch-shaped) PAHs. Furthermore, in PAHs of similar size they show significant variations in correlation with geometry. PAHs with the strictly linear geometry (chains) exhibit more electron correlation than PAHs with nonlinear geometries (sheets).