Two-electron reduced density matrix method for the electronic structure of correlated 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 $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.