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

arXiv:1411.0495 · Physical Review A

Propagating two-particle reduced density matrices without wave functions

Fabian Lackner, Iva Březinová, Takeshi Sato, Kenichi L. Ishikawa, Joachim Burgdörfer

Published Feb 11, 2015• 54 citations
arXivHTMLPDF

Abstract

Describing time-dependent many-body systems where correlation effects play an important role remains a major theoretical challenge. In this paper we develop a time-dependent many-body theory that is based on the two-particle reduced density matrix (2-RDM). We present a closed equation of motion for the 2-RDM by developing a reconstruction functional for the three-particle reduced density matrix (3-RDM) that preserves norm, energy, and spin symmetries during time propagation. We show that approximately enforcing $N$-representability during time evolution is essential for achieving stable solutions. As a prototypical test case which features long-range Coulomb interactions we employ the one-dimensional model for lithium hydride (LiH) in strong infrared laser fields. We probe both one-particle observables such as the time-dependent dipole moment and two-particle observables such as the pair density and mean electron-electron interaction energy. Our results are in very good agreement with numerically exact solutions for the $N$-electron wave function obtained from the multiconfigurational time-dependent Hartree-Fock method.

DOI
10.1103/physreva.91.023412

Related papers

The time-dependent two-particle reduced density matrix method: LiH in strong laser fields

Fabian Lackner, Iva Březinová, Takeshi Sato +2 more

We present a time-dependent many-body theory that is based on the propagation of the two-particle reduced density matrix (2-RDM) [1]. We apply our method to simulate the non-linear response of one-dimensional lithium hydride in femtosecond mid-IR laser pulses.

Journal of Physics Conference Series•Sep 7, 2015•2 citations

Exact solutions to the quantum many-body problem using the geminal density matrix

Nicholas J. Cox

It is virtually impossible to directly solve the Schrödinger equation for a many-electron wave function due to the exponential growth in degrees of freedom with increasing particle number. The two-body reduced density matrix (2-RDM) formalism reduces this coordinate dependence to that of four particles irrespective of the wave function's dimensionality, providing a promising path to solve the many-body problem. Unfortunately, errors arise in this approach because the 2-RDM cannot practically be constrained to guarantee that it corresponds to a valid wave function. Here we approach this so-called $N$-representability problem by expanding the 2-RDM in a complete basis of two-electron wave functions and studying the matrix formed by the expansion coefficients. This quantity, which we call the geminal density matrix (GDM), is found to evolve in time by a unitary transformation that preserves $N$-representability. This evolution law enables us to calculate eigenstates of strongly correlated systems by a fictitious adiabatic evolution in which the electron-electron interaction is slowly switched on. We show how this technique is used to diagonalize atomic Hamiltonians, finding that the problem reduces to the solution of $\sim N(N-1)/2$ two-electron eigenstates of the Helium atom on a grid of electron-electron interaction scaling factors.

arXiv:2112.11400•Dec 21, 2021•0 citations

Uncertainty relations and reduced density matrices: Mapping many-body quantum mechanics onto four particles

David A. Mazziotti, Robert Erdahl

For the description of ground-state correlation phenomena an accurate mapping of many-body quantum mechanics onto four particles is developed. The energy for a quantum system with no more than two-particle interactions may be expressed in terms of a two-particle reduced density matrix (2-RDM), but variational optimization of the 2-RDM requires that it corresponds to an N-particle wave function. We derive N-representability conditions on the 2-RDM that guarantee the validity of the uncertainty relations for all operators with two-particle interactions. One of these conditions is shown to be necessary and sufficient to make the RDM solutions of the dispersion condition equivalent to those from the contracted Schr\"odinger equation (CSE) [Mazziotti, Phys. Rev. A 57, 4219 (1998)]. In general, the CSE is a stronger N-representability condition than the dispersion condition because the CSE implies the dispersion condition as well as additional N-representability constraints from the Hellmann-Feynman theorem. Energy minimization subject to the representability constraints is performed for a boson model with 10, 30, and 75 particles. Even when traditional wave-function methods fail at large perturbations, the present method yields correlation energies within 2%.

Physical Review A•Mar 21, 2001•187 citations

Time-Dependent Quantum Many-Body Systems: Linear Response, Electronic Transport, and Reduced Density Matrices

Heiko Appel

Time-dependent density functional theory (TDDFT) provides a successful approach to calculate excitation energies of atomic and molecular systems. In part I of this work we present a double-pole approximation (DPA) to the response equations of TDDFT. The double-pole approximation provides an exact description of systems with two strongly coupled excitations which are isolated from the rest of the spectrum. In contrast to the traditional single- pole approximation of TDDFT the DPA also yields corrections to the Kohn-Sham oscillator strengths. Several critical pole separations can be identified, e.g. we find that the pole coupling can cause transitions to vanish entirely from the optical spectrum. We also demonstrate how to invert the double-pole solution which allows us to predict matrix elements of the exchange- correlation kernel fxc from experimental input. This can serve as benchmark for the construction of future approximations for the kernel fxc. Reduced density matrix functional theory (RDMFT) has emerged recently as promising candidate to treat strongly correlated electronic many-body systems beyond traditional density functional theory (DFT). The research within RDMFT was so far focussed on the static theory. In this work we attempt some first steps towards a time-dependent generalization of RDMFT. In part II we derive equations of motion for natural orbitals and occupation numbers. Using the equation of motion for the occupation numbers we show that an adiabatic extension of presently known ground-state functionals of static RDMFT always leads to occupation numbers which are constant in time. From the stationary conditions of the equations of motion for the N-body correlations (correlated parts of the N-body matrices) we derive a new class of ground-state functionals which can be used in static RDMFT. Applications are presented for a one-dimensional model system where the time-dependent many-body Schrödinger equation can be propagated numerically. We use optimal control theory to find optimized laser pulses for transitions in a model for atomic Helium. From the numerically exact correlated wavefunction we extract the exact time evolution of natural orbitals and occupation numbers for (i) laser-driven Helium and (ii) electron-ion scattering. Part III of this work considers time-dependent quantum transport within TDDFT. We present an algorithm for the calculation of extended eigenstates of single- particle Hamiltonians which is especially tailored to a finite-difference discretization of the Schrödinger equation. We consider the propagation of finite mesoscopic systems and demonstrate the limitations of such an approach. To overcome the shortcomings of a description of quantum transport in terms of finite systems we develop a time-propagation scheme for extended states which utilizes a mixed basis representation. Our discretization scheme allows to treat central device and lead regions on the same footing thus preventing artificial reflections at grid boundaries.

•Jan 1, 2007•10 citations

Quantum Many-Body Theory from a Solution of the N -Representability Problem

David A. Mazziotti

Here we present a many-body theory based on a solution of the N-representability problem in which the ground-state two-particle reduced density matrix (2-RDM) is determined directly without the many-particle wave function. We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM, which are required for its derivation from an N-particle density matrix, known as N-representability conditions. The approach produces a complete hierarchy of 2-RDM constraints that do not depend explicitly upon the higher RDMs or the wave function. By using the two-particle part of a unitary decomposition of higher-order constraint matrices, we can solve the energy minimization by semidefinite programming in a form where the low-rank structure of these matrices can be potentially exploited. We illustrate by computing the ground-state electronic energy and properties of the H_{8} ring.

arXiv:2304.08570•Physical Review Letters•Apr 12, 2023•41 citations

Capturing reduced-order quantum many-body dynamics out of equilibrium via neural ordinary differential equations

Patrick Egenlauf, Iva Březinová, Sabine Andergassen +1 more

Out-of-equilibrium quantum many-body systems exhibit rapid correlation buildup that underlies many emerging phenomena. Exact wave-function methods to describe this scale exponentially with particle number; simpler mean-field approaches neglect essential two-particle correlations. The time-dependent two-particle reduced density matrix (TD2RDM) formalism offers a middle ground by propagating the two-particle reduced density matrix (2RDM) and closing the BBGKY hierarchy with a reconstruction of the three-particle cumulant. But the validity and existence of time-local reconstruction functionals ignoring memory effects remain unclear across different dynamical regimes. We show that a neural ODE model trained on exact 2RDM data (no dimensionality reduction) can reproduce its dynamics without any explicit three-particle information -- but only in parameter regions where the Pearson correlation between the two- and three-particle cumulants is large. In the anti-correlated or uncorrelated regime, the neural ODE fails, indicating that no simple time-local functional of the instantaneous two-particle cumulant can capture the evolution. The magnitude of the time-averaged three-particle-correlation buildup appears to be the primary predictor of success: For a moderate correlation buildup, both neural ODE predictions and existing TD2RDM reconstructions are accurate, whereas stronger values lead to systematic breakdowns. These findings pinpoint the need for memory-dependent kernels in the three-particle cumulant reconstruction for the latter regime. Our results place the neural ODE as a model-agnostic diagnostic tool that maps the regime of applicability of cumulant expansion methods and guides the development of non-local closure schemes. More broadly, the ability to learn high-dimensional RDM dynamics from limited data opens a pathway to fast, data-driven simulation of correlated quantum matter.

cs.LGCross-listedarXiv:2512.13913•Machine Learning Science and Technology•Dec 17, 2025•0 citations
Describing time-dependent many-body systems where correlation effects play an important role remains a major theoretical challenge. In this paper we develop a time-dependent many-body theory that is based on the two-particle reduced density matrix (2-RDM). We present a closed equation of motion for the 2-RDM by developing a reconstruction functional for the three-particle reduced density matrix (3-RDM) that preserves norm, energy, and spin symmetries during time propagation. We show that approximately enforcing NNN-representability during time evolution is essential for achieving stable solutions. As a prototypical test case which features long-range Coulomb interactions we employ the one-dimensional model for lithium hydride (LiH) in strong infrared laser fields. We probe both one-particle observables such as the time-dependent dipole moment and two-particle observables such as the pair density and mean electron-electron interaction energy. Our results are in very good agreement with numerically exact solutions for the NNN-electron wave function obtained from the multiconfigurational time-dependent Hartree-Fock method.
It is virtually impossible to directly solve the Schrödinger equation for a many-electron wave function due to the exponential growth in degrees of freedom with increasing particle number. The two-body reduced density matrix (2-RDM) formalism reduces this coordinate dependence to that of four particles irrespective of the wave function's dimensionality, providing a promising path to solve the many-body problem. Unfortunately, errors arise in this approach because the 2-RDM cannot practically be constrained to guarantee that it corresponds to a valid wave function. Here we approach this so-called NNN-representability problem by expanding the 2-RDM in a complete basis of two-electron wave functions and studying the matrix formed by the expansion coefficients. This quantity, which we call the geminal density matrix (GDM), is found to evolve in time by a unitary transformation that preserves NNN-representability. This evolution law enables us to calculate eigenstates of strongly correlated systems by a fictitious adiabatic evolution in which the electron-electron interaction is slowly switched on. We show how this technique is used to diagonalize atomic Hamiltonians, finding that the problem reduces to the solution of ∼N(N−1)/2\sim N(N-1)/2∼N(N−1)/2 two-electron eigenstates of the Helium atom on a grid of electron-electron interaction scaling factors.
For the description of ground-state correlation phenomena an accurate mapping of many-body quantum mechanics onto four particles is developed. The energy for a quantum system with no more than two-particle interactions may be expressed in terms of a two-particle reduced density matrix (2-RDM), but variational optimization of the 2-RDM requires that it corresponds to an N-particle wave function. We derive N-representability conditions on the 2-RDM that guarantee the validity of the uncertainty relations for all operators with two-particle interactions. One of these conditions is shown to be necessary and sufficient to make the RDM solutions of the dispersion condition equivalent to those from the contracted Schr\"odinger equation (CSE) [Mazziotti, Phys. Rev. A 57, 4219 (1998)]. In general, the CSE is a stronger N-representability condition than the dispersion condition because the CSE implies the dispersion condition as well as additional N-representability constraints from the Hellmann-Feynman theorem. Energy minimization subject to the representability constraints is performed for a boson model with 10, 30, and 75 particles. Even when traditional wave-function methods fail at large perturbations, the present method yields correlation energies within 2%.
Time-dependent density functional theory (TDDFT) provides a successful approach to calculate excitation energies of atomic and molecular systems. In part I of this work we present a double-pole approximation (DPA) to the response equations of TDDFT. The double-pole approximation provides an exact description of systems with two strongly coupled excitations which are isolated from the rest of the spectrum. In contrast to the traditional single- pole approximation of TDDFT the DPA also yields corrections to the Kohn-Sham oscillator strengths. Several critical pole separations can be identified, e.g. we find that the pole coupling can cause transitions to vanish entirely from the optical spectrum. We also demonstrate how to invert the double-pole solution which allows us to predict matrix elements of the exchange- correlation kernel fxc from experimental input. This can serve as benchmark for the construction of future approximations for the kernel fxc. Reduced density matrix functional theory (RDMFT) has emerged recently as promising candidate to treat strongly correlated electronic many-body systems beyond traditional density functional theory (DFT). The research within RDMFT was so far focussed on the static theory. In this work we attempt some first steps towards a time-dependent generalization of RDMFT. In part II we derive equations of motion for natural orbitals and occupation numbers. Using the equation of motion for the occupation numbers we show that an adiabatic extension of presently known ground-state functionals of static RDMFT always leads to occupation numbers which are constant in time. From the stationary conditions of the equations of motion for the N-body correlations (correlated parts of the N-body matrices) we derive a new class of ground-state functionals which can be used in static RDMFT. Applications are presented for a one-dimensional model system where the time-dependent many-body Schrödinger equation can be propagated numerically. We use optimal control theory to find optimized laser pulses for transitions in a model for atomic Helium. From the numerically exact correlated wavefunction we extract the exact time evolution of natural orbitals and occupation numbers for (i) laser-driven Helium and (ii) electron-ion scattering. Part III of this work considers time-dependent quantum transport within TDDFT. We present an algorithm for the calculation of extended eigenstates of single- particle Hamiltonians which is especially tailored to a finite-difference discretization of the Schrödinger equation. We consider the propagation of finite mesoscopic systems and demonstrate the limitations of such an approach. To overcome the shortcomings of a description of quantum transport in terms of finite systems we develop a time-propagation scheme for extended states which utilizes a mixed basis representation. Our discretization scheme allows to treat central device and lead regions on the same footing thus preventing artificial reflections at grid boundaries.
Here we present a many-body theory based on a solution of the N-representability problem in which the ground-state two-particle reduced density matrix (2-RDM) is determined directly without the many-particle wave function. We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM, which are required for its derivation from an N-particle density matrix, known as N-representability conditions. The approach produces a complete hierarchy of 2-RDM constraints that do not depend explicitly upon the higher RDMs or the wave function. By using the two-particle part of a unitary decomposition of higher-order constraint matrices, we can solve the energy minimization by semidefinite programming in a form where the low-rank structure of these matrices can be potentially exploited. We illustrate by computing the ground-state electronic energy and properties of the H_{8} ring.
Out-of-equilibrium quantum many-body systems exhibit rapid correlation buildup that underlies many emerging phenomena. Exact wave-function methods to describe this scale exponentially with particle number; simpler mean-field approaches neglect essential two-particle correlations. The time-dependent two-particle reduced density matrix (TD2RDM) formalism offers a middle ground by propagating the two-particle reduced density matrix (2RDM) and closing the BBGKY hierarchy with a reconstruction of the three-particle cumulant. But the validity and existence of time-local reconstruction functionals ignoring memory effects remain unclear across different dynamical regimes. We show that a neural ODE model trained on exact 2RDM data (no dimensionality reduction) can reproduce its dynamics without any explicit three-particle information -- but only in parameter regions where the Pearson correlation between the two- and three-particle cumulants is large. In the anti-correlated or uncorrelated regime, the neural ODE fails, indicating that no simple time-local functional of the instantaneous two-particle cumulant can capture the evolution. The magnitude of the time-averaged three-particle-correlation buildup appears to be the primary predictor of success: For a moderate correlation buildup, both neural ODE predictions and existing TD2RDM reconstructions are accurate, whereas stronger values lead to systematic breakdowns. These findings pinpoint the need for memory-dependent kernels in the three-particle cumulant reconstruction for the latter regime. Our results place the neural ODE as a model-agnostic diagnostic tool that maps the regime of applicability of cumulant expansion methods and guides the development of non-local closure schemes. More broadly, the ability to learn high-dimensional RDM dynamics from limited data opens a pathway to fast, data-driven simulation of correlated quantum matter.