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:2601.21088

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

Thomas Blommel, M. Rey Lambert, Michael Alexander Kurniawan, Annabelle Canestraight, Vojtěch Vlček

Published Jan 28, 2026• 0 citations
arXivHTMLPDF

Abstract

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established non-equilibrium Green's function approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal Green's functions, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz (GKBA), and exact diagonalization for small systems using second order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the Green's functions, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich non-equilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and bandgap renormalization, which are largely washed out in self-consistent approaches, yet are present in the exact solutions. Our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

DOI
10.48550/arxiv.2601.21088

Related papers

Real-time Dyson expansion for the nonequilibrium $GW$ approximation: Correlated spectra and dependence on the reference propagator

Erik Schroedter, Jan‐Philip Joost, M. Bönitz +1 more

Electronic spectra provide direct insight into the excitations and correlations of condensed matter systems. Their description requires electron correlations beyond mean field. In equilibrium, the $GW$ approximation has become the method of choice for many materials. Extending this approximation to nonequilibrium, however, is challenging. Full two-time $GW$ simulations within the nonequilibrium Green's functions approach scale at least cubically with propagation time, whereas improved scaling schemes such as time-local (adiabatic) approximations or the generalized Kadanoff--Baym ansatz usually retain only mean-field character in the spectra. The recently introduced real-time Dyson expansion recovers dynamical correlations in the spectrum at time-linear cost, but has so far been restricted to the second-order Born approximation with mean-field reference propagators. Here we formulate and implement the RT-DE for the nonequilibrium $GW$ self-energy with reference propagators of general form whose off-diagonal evolution is time local. In particular, we assess Hartree--Fock propagators, propagators with statically screened exchange with nonequilibrium screening, and a correlated propagator based on the Hartree--Fock GKBA. Benchmarks against exact diagonalization for a driven two-band lattice model with long-range interactions show that the mean-field and statically screened references yield the most accurate spectra, including satellite structures absent at mean-field level, whereas the GKBA-based approach best captures scattering-induced occupation dynamics, but artificially broadens and splits spectral peaks. When applied to large systems that are beyond the reach of exact methods, the scheme resolves the excitonic replica of the valence band and satellites identified as exciton shake-up, as well as their reshaping with increasing excitation density.

arXiv:2609.08932•Sep 8, 2026•0 citations

A Real-time Dyson Expansion Scheme: Efficient Inclusion of Dynamical Correlations in Non-equilibrium Spectral Properties

Cian C. Reeves, Vojtěch Vlček

Time-resolved photoemission spectroscopy is the key technique to probe the real-time non-equilibrium dynamics of electronic states. Theoretical predictions of the time dependent spectral function for realistic systems is however, a challenge. Employing the Kadanoff-Baym equations to find this quantity results in a cubic scaling in the total number of time steps, quickly becoming prohibitive and often fail quantitatively and even qualitatively. In comparison, mean-field methods have more favorable numerical scaling both in the number of time steps and in the complexity associated with the cost of evolving for a single time step, however they miss key spectral properties such as emergent spectral features. Here we present a scheme that allows for the inclusion of dynamical correlations to the spectral function while maintaining the same scaling in the number of time steps as for mean-field approaches, while capturing the emergent physics. Further, the scheme can be efficiently implemented on top of equilibrium real-time many-body perturbation theory schemes and codes. We see excellent agreement with exact results for test systems. Furthermore we exemplify the method on a periodic system and demonstrate clear evidence that our proposed scheme produces complex spectral features including excitonic band replicas, features that are not observed using static mean-field approaches.

arXiv:2403.07155•Mar 11, 2024•1 citation

Real-Time Dyson-Expansion Scheme: Efficient Inclusion of Dynamical Correlations in Nonequilibrium Spectral Properties

Cian C. Reeves, Vojtěch Vlček

Time-resolved photoemission spectroscopy is the key technique to probe the real-time nonequilibrium dynamics of electronic states. Theoretical predictions of the time dependent spectral function for realistic systems is however, a challenge. Employing the Kadanoff-Baym equations to find this quantity results in a cubic scaling in the total number of time steps, quickly becoming prohibitive and often fail quantitatively and even qualitatively. In comparison, mean-field methods have more favorable numerical scaling both in the number of time steps and in the complexity associated with the cost of evolving for a single time step, however they miss key spectral properties such as emergent spectral features. Here we present a scheme that allows for the inclusion of dynamical correlations to the spectral function while maintaining the same scaling in the number of time steps as for mean-field approaches, while capturing the emergent physics. Further, the scheme can be efficiently implemented on top of equilibrium real-time many-body perturbation theory schemes and codes. We see excellent agreement with exact results for test systems. Furthermore, we exemplify the method on a periodic system and demonstrate clear evidence that our proposed scheme produces complex spectral features including excitonic band replicas, features that are not observed using static mean-field approaches.

Physical Review Letters•Nov 26, 2024•14 citations

Spectral properties of correlated quantum wires and carbon nanotubes within the Generalized Kadanoff-Baym Ansatz

F. Cosco, N. Walter Talarico, Riku Tuovinen +1 more

We investigate the spectral properties of an open interacting system by solving the Generalized Kadanoff-Baym Ansatz (GKBA) master equation for the single-particle density matrix, namely the time-diagonal lesser Green's function. To benchmark its validity, we compare the solution obtained within the GKBA with the solution of the Dyson equation (equivalently the full Kadanoff-Baym equations). In both approaches, we treat the interaction within the self-consistent second-order Born approximation, whereas the GKBA still retains the retarded propagator calculated at the Hartree-Fock level. We consider the case of two leads connected through a central correlated region where particles can interact and exploit the stationary particle current at the boundary of the junction as a probe of the spectral features of the system. In this work, as an example, we take the central region to be a one-dimensional quantum wire and a two-dimensional carbon nanotube and show that the solution of the GKBA master equation well captures their spectral features. Our result demonstrates that, even when the propagator used is at the Hartree-Fock level, the GBKA solution retains the main spectral features of the self-energy used.

arXiv:2007.08901•Jul 17, 2020•5 citations

Hubbard nanoclusters far from equilibrium

S. Hermanns, Niclas Schlünzen, M. Bönitz

The Hubbard model is a prototype for strongly correlated many-particle systems, including electrons in condensed matter and molecules, as well as for fermions or bosons in optical lattices. While the equilibrium properties of these systems have been studied in detail, the nonequilibrium dynamics following a strong nonperturbative excitation only recently came into the focus of experiments and theory. It is of particular interest how the dynamics depend on the coupling strength and on the particle number and whether there exist universal features in the time evolution. Here, we present results for the dynamics of finite Hubbard clusters based on a self-consistent nonequilibrium Green functions (NEGF) approach invoking the generalized Kadanoff-Baym ansatz (GKBA). We discuss the conserving properties of the GKBA with Hartree-Fock propagators in detail and present a generalized form of the energy conservation criterion of Baym and Kadanoff for NEGF. Furthermore, we demonstrate that the HF-GKBA cures some artifacts of prior two-time NEGF simulations. Besides, this approach substantially speeds up the numerical calculations and thus presents the capability to study comparatively large systems and to extend the analysis to long times allowing for an accurate computation of the excitation spectrum via time propagation. Our data obtained within the second Born approximation compare favorably with exact diagonalization results (available for up to 13 particles) and are expected to have predictive capability for substantially larger systems in the weak-coupling limit.

arXiv:1402.7300•Physical Review B•Sep 5, 2014•80 citations

Unified simulation framework for correlated driven-dissipative quantum dynamics

Thomas Blommel, Enrico Perfetto, Gianluca Stefanucci +1 more

Time-resolved photoemission spectroscopy provides a unique and direct way to explore the real-time nonequilibrium dynamics of electrons and holes. The formal theory of spectral function evolution requires inclusion of electronic correlations and dissipation, which are challenging due to the associated long simulation timescales which translate to a high computational cost. Recent methodological developments, namely, the real-time Dyson expansion, as well as theoretical developments of many-body perturbation theory for dissipative systems, have allowed for the study of driven-dissipative interacting quantum systems. In this work, we implement the hitherto unrealized study of driven dynamically correlated systems and utilize these methods and perturbative expansions. We illustrate the combined formalism on a prototypical two-band semiconductor model with long-range density-density Coulomb interactions. We show that the intraband thermalization of conduction band electrons induces nontrivial time-dependent changes in the system's band structure and a time-evolving band-gap renormalization (with a reduction by up to 8.4%). We show that the qualitative features are preserved for a variety of parameters, discuss the corresponding spectral dynamics, and provide an outlook on the introduced simulation framework, which enables treating electron-electron scattering and dissipation effects on equal footing.

Physical review. B./Physical review. B•Sep 22, 2025•5 citations
We present a systematic comparison of the real-time Dyson expansion (RTDE) with established non-equilibrium Green's function approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal Green's functions, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz (GKBA), and exact diagonalization for small systems using second order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the Green's functions, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich non-equilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and bandgap renormalization, which are largely washed out in self-consistent approaches, yet are present in the exact solutions. Our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.
Real-time Dyson expansion for the nonequilibrium GWGWGW approximation: Correlated spectra and dependence on the reference propagator
Electronic spectra provide direct insight into the excitations and correlations of condensed matter systems. Their description requires electron correlations beyond mean field. In equilibrium, the GWGWGW approximation has become the method of choice for many materials. Extending this approximation to nonequilibrium, however, is challenging. Full two-time GWGWGW simulations within the nonequilibrium Green's functions approach scale at least cubically with propagation time, whereas improved scaling schemes such as time-local (adiabatic) approximations or the generalized Kadanoff--Baym ansatz usually retain only mean-field character in the spectra. The recently introduced real-time Dyson expansion recovers dynamical correlations in the spectrum at time-linear cost, but has so far been restricted to the second-order Born approximation with mean-field reference propagators. Here we formulate and implement the RT-DE for the nonequilibrium GWGWGW self-energy with reference propagators of general form whose off-diagonal evolution is time local. In particular, we assess Hartree--Fock propagators, propagators with statically screened exchange with nonequilibrium screening, and a correlated propagator based on the Hartree--Fock GKBA. Benchmarks against exact diagonalization for a driven two-band lattice model with long-range interactions show that the mean-field and statically screened references yield the most accurate spectra, including satellite structures absent at mean-field level, whereas the GKBA-based approach best captures scattering-induced occupation dynamics, but artificially broadens and splits spectral peaks. When applied to large systems that are beyond the reach of exact methods, the scheme resolves the excitonic replica of the valence band and satellites identified as exciton shake-up, as well as their reshaping with increasing excitation density.
Time-resolved photoemission spectroscopy is the key technique to probe the real-time non-equilibrium dynamics of electronic states. Theoretical predictions of the time dependent spectral function for realistic systems is however, a challenge. Employing the Kadanoff-Baym equations to find this quantity results in a cubic scaling in the total number of time steps, quickly becoming prohibitive and often fail quantitatively and even qualitatively. In comparison, mean-field methods have more favorable numerical scaling both in the number of time steps and in the complexity associated with the cost of evolving for a single time step, however they miss key spectral properties such as emergent spectral features. Here we present a scheme that allows for the inclusion of dynamical correlations to the spectral function while maintaining the same scaling in the number of time steps as for mean-field approaches, while capturing the emergent physics. Further, the scheme can be efficiently implemented on top of equilibrium real-time many-body perturbation theory schemes and codes. We see excellent agreement with exact results for test systems. Furthermore we exemplify the method on a periodic system and demonstrate clear evidence that our proposed scheme produces complex spectral features including excitonic band replicas, features that are not observed using static mean-field approaches.
Time-resolved photoemission spectroscopy is the key technique to probe the real-time nonequilibrium dynamics of electronic states. Theoretical predictions of the time dependent spectral function for realistic systems is however, a challenge. Employing the Kadanoff-Baym equations to find this quantity results in a cubic scaling in the total number of time steps, quickly becoming prohibitive and often fail quantitatively and even qualitatively. In comparison, mean-field methods have more favorable numerical scaling both in the number of time steps and in the complexity associated with the cost of evolving for a single time step, however they miss key spectral properties such as emergent spectral features. Here we present a scheme that allows for the inclusion of dynamical correlations to the spectral function while maintaining the same scaling in the number of time steps as for mean-field approaches, while capturing the emergent physics. Further, the scheme can be efficiently implemented on top of equilibrium real-time many-body perturbation theory schemes and codes. We see excellent agreement with exact results for test systems. Furthermore, we exemplify the method on a periodic system and demonstrate clear evidence that our proposed scheme produces complex spectral features including excitonic band replicas, features that are not observed using static mean-field approaches.
We investigate the spectral properties of an open interacting system by solving the Generalized Kadanoff-Baym Ansatz (GKBA) master equation for the single-particle density matrix, namely the time-diagonal lesser Green's function. To benchmark its validity, we compare the solution obtained within the GKBA with the solution of the Dyson equation (equivalently the full Kadanoff-Baym equations). In both approaches, we treat the interaction within the self-consistent second-order Born approximation, whereas the GKBA still retains the retarded propagator calculated at the Hartree-Fock level. We consider the case of two leads connected through a central correlated region where particles can interact and exploit the stationary particle current at the boundary of the junction as a probe of the spectral features of the system. In this work, as an example, we take the central region to be a one-dimensional quantum wire and a two-dimensional carbon nanotube and show that the solution of the GKBA master equation well captures their spectral features. Our result demonstrates that, even when the propagator used is at the Hartree-Fock level, the GBKA solution retains the main spectral features of the self-energy used.
The Hubbard model is a prototype for strongly correlated many-particle systems, including electrons in condensed matter and molecules, as well as for fermions or bosons in optical lattices. While the equilibrium properties of these systems have been studied in detail, the nonequilibrium dynamics following a strong nonperturbative excitation only recently came into the focus of experiments and theory. It is of particular interest how the dynamics depend on the coupling strength and on the particle number and whether there exist universal features in the time evolution. Here, we present results for the dynamics of finite Hubbard clusters based on a self-consistent nonequilibrium Green functions (NEGF) approach invoking the generalized Kadanoff-Baym ansatz (GKBA). We discuss the conserving properties of the GKBA with Hartree-Fock propagators in detail and present a generalized form of the energy conservation criterion of Baym and Kadanoff for NEGF. Furthermore, we demonstrate that the HF-GKBA cures some artifacts of prior two-time NEGF simulations. Besides, this approach substantially speeds up the numerical calculations and thus presents the capability to study comparatively large systems and to extend the analysis to long times allowing for an accurate computation of the excitation spectrum via time propagation. Our data obtained within the second Born approximation compare favorably with exact diagonalization results (available for up to 13 particles) and are expected to have predictive capability for substantially larger systems in the weak-coupling limit.
Time-resolved photoemission spectroscopy provides a unique and direct way to explore the real-time nonequilibrium dynamics of electrons and holes. The formal theory of spectral function evolution requires inclusion of electronic correlations and dissipation, which are challenging due to the associated long simulation timescales which translate to a high computational cost. Recent methodological developments, namely, the real-time Dyson expansion, as well as theoretical developments of many-body perturbation theory for dissipative systems, have allowed for the study of driven-dissipative interacting quantum systems. In this work, we implement the hitherto unrealized study of driven dynamically correlated systems and utilize these methods and perturbative expansions. We illustrate the combined formalism on a prototypical two-band semiconductor model with long-range density-density Coulomb interactions. We show that the intraband thermalization of conduction band electrons induces nontrivial time-dependent changes in the system's band structure and a time-evolving band-gap renormalization (with a reduction by up to 8.4%). We show that the qualitative features are preserved for a variety of parameters, discuss the corresponding spectral dynamics, and provide an outlook on the introduced simulation framework, which enables treating electron-electron scattering and dissipation effects on equal footing.