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:1602.03723 · International Journal of Quantum Chemistry

Density‐matrix propagation driven by semiclassical correlation

Peter D. T. A. Elliott, Neepa T. Maitra

Published Feb 1, 2016• 6 citations
arXivHTMLPDF

Abstract

Methods based on propagation of the one‐body reduced density‐matrix hold much promise for the simulation of correlated many‐electron dynamics far from equilibrium, but difficulties with finding good approximations for the interaction term in its equation of motion have so far impeded their application. These difficulties include the violation of fundamental physical principles such as energy conservation, positivity conditions on the density, or unchanging natural orbital occupation numbers. We review some of the recent efforts to confront these problems, and explore a semiclassical approximation for electron correlation coupled to time‐dependent Hartree‐Fock propagation. We find that this approach captures changing occupation numbers, and excitations to doubly‐excited states, improving over TDHF and adiabatic approximations in density‐matrix propagation. However, it does not guarantee N‐representability of the density‐matrix, consequently resulting sometimes in violation of positivity conditions, even though a purely semiclassical treatment preserves these conditions. © 2016 Wiley Periodicals, Inc.

DOI
10.1002/qua.25087

Related papers

Challenges in truncating the hierarchy of time-dependent reduced density matrices equations

A. Akbari, M. J. Hashemi, Ángel Rubio +2 more

In this work, we analyze the Born, Bogoliubov, Green, Kirkwood, and Yvon (BBGKY) hierarchy of equations for describing the full time evolution of a many-body fermionic system in terms of its reduced density matrices (at all orders). We provide an exhaustive study of the challenges and open problems linked to the truncation of such a hierarchy of equations to make them practically applicable. We restrict our analysis to the coupled evolution of the one- and two-body reduced density matrices, where higher-order correlation effects are embodied into the approximation used to close the equations. We prove that within this approach, the number of electrons and total energy are conserved, regardless of the employed approximation. Further, we demonstrate that although most of the truncation schemes available in the literature give acceptable ground-state energy, when applied to describe driven electron dynamics, they exhibit undesirable and unphysical behavior, e.g., violation and even divergence in local electronic density, both in weakly and strongly correlated regimes. We illustrate and analyze these problems within the few-site Hubbard model. The model can be solved exactly and provides a unique reference for our detailed study of electron dynamics for different values of interaction, different initial conditions, and the large set of approximations considered here. Moreover, we study the role of compatibility between two hierarchical equations and positive semidefiniteness of reduced density matrices in the instability of electron dynamics. We show that even if the used approximation holds the compatibility, electron dynamics can still diverge when positive definitiveness is violated. We propose some partial solutions of such problems and point out the main paths for future work in order to make this approach applicable for the description of the correlated electron dynamics in complex systems.

arXiv:1204.4395•Physical Review B•Jun 12, 2012•48 citations

Semiclassical Electron Correlation in Density-Matrix Time Propagation

Arun K. Rajam, Izabela Raczkowska, Neepa T. Maitra

Lack of memory (locality in time) is a major limitation of almost all present time-dependent density-functional approximations. By using semiclassical dynamics to compute correlation effects in the time propagation of the density matrix, we incorporate memory, including initial-state dependence, as well as changing occupation numbers, and predict more observables in strong-field applications.

arXiv:1005.1341•Physical Review Letters•Sep 7, 2010•24 citations

Gaussian phase-space representation of fermion dynamics: Beyond the time-dependent Hartree-Fock approximation

Saar Rahav, Shaul Mukamel

A Gaussian operator representation for the many-body density matrix of fermionic systems, developed by Corney and Drummond [Phys. Rev. Lett. 93, 260401 (2004)], is used to derive approximate decoupling schemes for their dynamics. In this approach the reduced single electron density-matrix elements serve as stochastic variables which satisfy an exact Fokker-Planck equation. The number of variables scales as $\ensuremath{\sim}{N}^{2}$ rather than $\ensuremath{\sim}\text{exp}(N)$ with the basis set size, and the time-dependent Hartree-Fock approximation (TDHF) is recovered in the ``classical'' limit. An approximate closed set of equations of motion for the one- and two-particle reduced density matrices provides a direct generalization of the TDHF.

arXiv:0902.3227•Physical Review B•Apr 10, 2009•10 citations

Density matrix methods for semiconductor Coulomb dynamics

James W. Dufty, Chang Sub Kim, M. Bönitz +1 more

Current experiments on semiconductor devices using femtosecond lasers provide new theoretical challenges for the description of charge carrier dynamics. Among the new features of such experiments are states driven very far from equilibrium and probes on time scales short compared to scattering and other characteristic material relaxation times. Standard many-body methods must be modified and extended to accommodate these features. We propose that the quantum hierarchy for reduced density operators is an ideal formulation of such initial value problems and describe how the dominant effects of exchange and charge correlations can be accounted for in a simple and physically transparent closure of the hierarchy of equations. The transformations, approximations, and interpretation can be accomplished independent of any particular matrix representation. Decomposition into kinetic equations for band occupation densities and polarization densities follows in a straightforward way after the many-body problem has been brought under control. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 65: 929–940, 1997

International Journal of Quantum Chemistry•Jan 1, 1997•5 citations

Adiabatic approximation in time-dependent reduced-density-matrix functional theory

Ryan Requist, Oleg Pankratov

With the aim of describing real-time electron dynamics, we introduce an adiabatic approximation for the equation of motion of the one-body reduced density matrix (one-matrix). The eigenvalues of the one-matrix, which represent the occupation numbers of single-particle orbitals, are obtained from the constrained minimization of the instantaneous ground-state energy functional rather than from their dynamical equations. The performance of the approximation vis-\`a-vis nonadiabatic effects is assessed in real-time simulations of a two-site Hubbard model. Due to Landau-Zener-type transitions, the system evolves into a nonstationary state with persistent oscillations in the observables. The amplitude of the oscillations displays a strongly nonmonotonic dependence on the strength of the electron-electron interaction and the rate of variation of the external potential. We interpret an associated resonance behavior in the phase of the oscillations in terms of ``scattering'' with spectator energy levels. To clarify the motivation for the minimization condition, we derive a sequence of energy functionals ${E}_{v}^{(n)}$, for which the corresponding sequence of minimizing one-matrices is asymptotic to the exact one-matrix in the adiabatic limit.

arXiv:0911.0945•Physical Review A•Apr 26, 2010•46 citations

Methods for electronic‐structure calculations: Overview from a reduced‐density‐matrix point of view

P. Ziesche, Ferenc Tasnádi

Abstract 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, and 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 quasiparticle 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, possibly the solutions of an approximate effective two‐body equation, and we sketch the idea of how its contraction sum rule can be used for a variational treatment. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004

International Journal of Quantum Chemistry•Jan 1, 2004•2 citations
Methods based on propagation of the one‐body reduced density‐matrix hold much promise for the simulation of correlated many‐electron dynamics far from equilibrium, but difficulties with finding good approximations for the interaction term in its equation of motion have so far impeded their application. These difficulties include the violation of fundamental physical principles such as energy conservation, positivity conditions on the density, or unchanging natural orbital occupation numbers. We review some of the recent efforts to confront these problems, and explore a semiclassical approximation for electron correlation coupled to time‐dependent Hartree‐Fock propagation. We find that this approach captures changing occupation numbers, and excitations to doubly‐excited states, improving over TDHF and adiabatic approximations in density‐matrix propagation. However, it does not guarantee N‐representability of the density‐matrix, consequently resulting sometimes in violation of positivity conditions, even though a purely semiclassical treatment preserves these conditions. © 2016 Wiley Periodicals, Inc.
In this work, we analyze the Born, Bogoliubov, Green, Kirkwood, and Yvon (BBGKY) hierarchy of equations for describing the full time evolution of a many-body fermionic system in terms of its reduced density matrices (at all orders). We provide an exhaustive study of the challenges and open problems linked to the truncation of such a hierarchy of equations to make them practically applicable. We restrict our analysis to the coupled evolution of the one- and two-body reduced density matrices, where higher-order correlation effects are embodied into the approximation used to close the equations. We prove that within this approach, the number of electrons and total energy are conserved, regardless of the employed approximation. Further, we demonstrate that although most of the truncation schemes available in the literature give acceptable ground-state energy, when applied to describe driven electron dynamics, they exhibit undesirable and unphysical behavior, e.g., violation and even divergence in local electronic density, both in weakly and strongly correlated regimes. We illustrate and analyze these problems within the few-site Hubbard model. The model can be solved exactly and provides a unique reference for our detailed study of electron dynamics for different values of interaction, different initial conditions, and the large set of approximations considered here. Moreover, we study the role of compatibility between two hierarchical equations and positive semidefiniteness of reduced density matrices in the instability of electron dynamics. We show that even if the used approximation holds the compatibility, electron dynamics can still diverge when positive definitiveness is violated. We propose some partial solutions of such problems and point out the main paths for future work in order to make this approach applicable for the description of the correlated electron dynamics in complex systems.
A Gaussian operator representation for the many-body density matrix of fermionic systems, developed by Corney and Drummond [Phys. Rev. Lett. 93, 260401 (2004)], is used to derive approximate decoupling schemes for their dynamics. In this approach the reduced single electron density-matrix elements serve as stochastic variables which satisfy an exact Fokker-Planck equation. The number of variables scales as $\ensuremath{\sim}{N}^{2}$ rather than $\ensuremath{\sim}\text{exp}(N)$ with the basis set size, and the time-dependent Hartree-Fock approximation (TDHF) is recovered in the ``classical'' limit. An approximate closed set of equations of motion for the one- and two-particle reduced density matrices provides a direct generalization of the TDHF.
Current experiments on semiconductor devices using femtosecond lasers provide new theoretical challenges for the description of charge carrier dynamics. Among the new features of such experiments are states driven very far from equilibrium and probes on time scales short compared to scattering and other characteristic material relaxation times. Standard many-body methods must be modified and extended to accommodate these features. We propose that the quantum hierarchy for reduced density operators is an ideal formulation of such initial value problems and describe how the dominant effects of exchange and charge correlations can be accounted for in a simple and physically transparent closure of the hierarchy of equations. The transformations, approximations, and interpretation can be accomplished independent of any particular matrix representation. Decomposition into kinetic equations for band occupation densities and polarization densities follows in a straightforward way after the many-body problem has been brought under control. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 65: 929–940, 1997
With the aim of describing real-time electron dynamics, we introduce an adiabatic approximation for the equation of motion of the one-body reduced density matrix (one-matrix). The eigenvalues of the one-matrix, which represent the occupation numbers of single-particle orbitals, are obtained from the constrained minimization of the instantaneous ground-state energy functional rather than from their dynamical equations. The performance of the approximation vis-\`a-vis nonadiabatic effects is assessed in real-time simulations of a two-site Hubbard model. Due to Landau-Zener-type transitions, the system evolves into a nonstationary state with persistent oscillations in the observables. The amplitude of the oscillations displays a strongly nonmonotonic dependence on the strength of the electron-electron interaction and the rate of variation of the external potential. We interpret an associated resonance behavior in the phase of the oscillations in terms of ``scattering'' with spectator energy levels. To clarify the motivation for the minimization condition, we derive a sequence of energy functionals Ev(n){E}_{v}^{(n)}Ev(n)​, for which the corresponding sequence of minimizing one-matrices is asymptotic to the exact one-matrix in the adiabatic limit.
Abstract 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, and 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 quasiparticle 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, possibly the solutions of an approximate effective two‐body equation, and we sketch the idea of how its contraction sum rule can be used for a variational treatment. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004