A methodology is introduced for the systematic treatment of electron-electron correlations in solids and other interacting quantum
N-particle systems. The method is developed within the framework of electronic structure theory (band theory) but, in contrast to conventional approaches, which are based on the single-particle picture, it is formulated within a many-particle picture in which
n particles in
d-dimensional phase space are treated as a single particle in a phase space of
nd dimensions. In this phase space, interparticle interactions appear as external potentials allowing the treatment of the system of particles through the use of single-particle methods, while at the same time allowing a systematic, direct, and nonperturbative treatment of interparticle interactions. The method makes use of the invariance of the Hamiltonian describing an interacting-particle system under partitioning into subsystems of
n particles. This treatment leads to exact results in the limit $n\ensuremath{\rightarrow}N$. Based on such partitioning, we propose a generalization of density functional theory and an appropriately defined local density approximation to treat the interactions between the
n-particle units in a system of
N> n particles. This approach yields
n-particle correlated densities and
n-particle states which can be used in an analysis of the electronic properties of materials, such as total energy, phase stability, electronic transport, and others. We use the formal construct of multiple-scattering theory to develop the method for the calculation of the two-particle electronic structure of a solid and the corresponding total energy of the ground state. We also illustrate some of the properties of the method in terms of a Hubbard model Hamiltonian on a linear ring. Various features of the method and further possible applications are presented in a discussion section.
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.
This manuscript introduces a methodology (within the Born-Oppenheimer picture) to compute electronic ground-state properties of molecules and solids/surfaces with fractionally occupied components. Given a user-defined division of the molecule into subsystems, our theory uses an auxiliary global Hamiltonian that is defined as the sum of subsystem Hamiltonians, plus the spatial integral of a second-quantized local operator that allows the electrons to be transferred between subsystems. This electron transfer operator depends on a local potential that can be determined using density functional approximations and/or other techniques such as machine learning. The present framework employs superpositions of tensor-product wave functions, which can satisfy size consistency and avoid spurious fractional charges at large bond distances. The electronic population of each subsystem is in general a positive real number and is obtained from wave-function amplitudes, which are calculated by means of ground-state matrix diagonalization (or matrix propagation in the time-dependent case). Our method can provide pathways to explore charge-transfer effects in environments where dividing the molecule into subsystems is convenient and to develop computationally affordable electronic structure algorithms.
In materials with strong electron correlation the proper treatment of local atomic physics described by orbital occupations is crucial. Reduced density-matrix functional theory is a natural extension of density functional theory for systems that are dominated by orbital physics. We review the current state of reduced density-matrix functional theory (RDMFT). For atomic structure relaxations or ab-initio molecular dynamics the combination of density functional theory (DFT) and dynamical mean-field theory (DMFT) possesses a number of disadvantages, like the cumbersome evaluation of forces. We therefore describe a method, DFT+RDMFT, that combines many-particle effects based on reduced density-matrix functional theory with a density functional-like framework. A recent development is the construction of density-matrix functionals directly from many-particle theory such as methods from quantum chemistry or many-particle Green’s functions. We present the underlying exact theorems and describe current progress towards quantitative functionals.
A method for calculating the Kohn-Sham exchange-correlation potential v(XC)(r) from a given electronic wave function is devised and implemented. It requires on input one- and two-electron reduced density matrices and involves construction of the generalized Fock matrix. The method is free from numerical limitations and basis-set artifacts of conventional schemes for constructing v(XC)(r) in which the potential is recovered from a given electron density, and is simpler than various many-body techniques. The chief significance of this development is that it allows one to directly probe the functional derivative of the true exchange-correlation energy functional and to rigorously test and improve various density-functional approximations.
The methods of quantum chemistry and solid state theory to solve the many-body problem are reviewed. We start with the definitions of reduced density matrices, their properties (contraction sum rules, spectral resolutions, cumulant expansion,
N-representability), and their determining equations (contracted Schrödinger equations) and we summarize recent extensions and generalizations of the traditional quantum chemical methods, of the density functional theory, and of the quasi-particle theory: from finite to extended systems (incremental method), from density to density matrix (density matrix functional theory), from weak to strong correlation (dynamical mean field theory), from homogeneous (Kimball-Overhauser approach) to inhomogeneous and finite systems. Measures of the correlation strength are discussed. The cumulant two-body reduced density matrix proves to be a key quantity. Its spectral resolution contains geminals, being possibly the solutions of an approximate effective two-body equation, and the idea is sketched of how its contraction sum rule can be used for a variational treatment.