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

Efficient non-equilibrium Green's function simulations of correlated electron-boson systems

Daniel Karlsson, Robert van Leeuwen, Y. Pavlyukh, Enrico Perfetto, Gianluca Stefanucci

Published Jun 26, 2020• 0 citations
OpenAlex

Abstract

This work contains a computationally efficient Green's function scheme for treating correlated electron-boson systems out of equilibrium. The numerical effort scales linearly with the simulation time, thus competing with that of the fastest quantum methods available. On a more fundamental level the proposed scheme fulfills energy balance as well as other fundamental conservation laws, and it can be systematically improved using standard many-body perturbation theory. We demonstrate its accuracy and stability through benchmarks against exact results in the Holstein model. However, our formulation equally applies to molecules and solids for first-principles transient calculations with quantum bosons, such as photons and phonons.

Related papers

Fast Green’s Function Method for Ultrafast Electron-Boson Dynamics

Daniel Karlsson, Robert van Leeuwen, Y. Pavlyukh +2 more

The interaction of electrons with quantized phonons and photons underlies the ultrafast dynamics of systems ranging from molecules to solids, and it gives rise to a plethora of physical phenomena experimentally accessible using time-resolved techniques. Green's function methods offer an invaluable interpretation tool since scattering mechanisms of growing complexity can be selectively incorporated in the theory. Currently, however, real-time Green's function simulations are either prohibitively expensive due to the cubic scaling with the propagation time or do neglect the feedback of electrons on the bosons, thus violating energy conservation. We put forward a computationally efficient Green's function scheme which overcomes both limitations. The numerical effort scales linearly with the propagation time while the simultaneous dressing of electrons and bosons guarantees the fulfillment of all fundamental conservation laws. We present a real-time study of the phonon-driven relaxation dynamics in an optically excited narrow band-gap insulator, highlighting the nonthermal behavior of the phononic degrees of freedom. Our formulation paves the way to first-principles simulations of electron-boson systems with unprecedented long propagation times.

arXiv:2006.14965•Physical Review Letters•Jul 15, 2021•73 citations

Time-linear scaling nonequilibrium Green's function methods for real-time simulations of interacting electrons and bosons. I. Formalism

Y. Pavlyukh, Enrico Perfetto, Daniel Karlsson +2 more

Simulations of interacting electrons and bosons out of equilibrium, starting from first principles and aiming at realistic multiscale scenarios, is a grand theoretical challenge. Here, using the formalism of nonequilibrium Green's functions and relying in a crucial way on the recently discovered time-linear formulation of the Kadanoff-Baym equations, we present a versatile toolbox for the simulation of correlated electron-boson dynamics. A large class of methods are available, from the Ehrenfest to the dressed $GD$ for the treatment of electron-boson interactions in combination with perturbative, i.e., Hartree-Fock and second-Born, or nonperturbative, i.e., $GW$ and $T$ matrices either without or with exchange effects, for the treatment of the Coulomb interaction. In all cases the numerical scaling is linear in time and the equations of motion satisfy all fundamental conservation laws.

arXiv:2111.06698•Physical review. B./Physical review. B•Mar 28, 2022•46 citations

Time-Linear Quantum Transport Simulations with Correlated Nonequilibrium Green’s Functions

Riku Tuovinen, Y. Pavlyukh, Enrico Perfetto +1 more

We present a time-linear scaling method to simulate open and correlated quantum systems out of equilibrium. The method inherits from many-body perturbation theory the possibility to choose selectively the most relevant scattering processes in the dynamics, thereby paving the way to the real-time characterization of correlated ultrafast phenomena in quantum transport. The open system dynamics is described in terms of an "embedding correlator" from which the time-dependent current can be calculated using the Meir-Wingreen formula. We show how to efficiently implement our approach through a simple grafting into recently proposed time-linear Green's function methods for closed systems. Electron-electron and electron-phonon interactions can be treated on equal footing while preserving all fundamental conservation laws.

arXiv:2211.15635•Physical Review Letters•Jun 16, 2023•32 citations

Accelerating Nonequilibrium Green Functions Simulations: The G1–G2 Scheme and Beyond

M. Bönitz, Jan‐Philip Joost, Christopher Makait +3 more

The theory of nonequilibrium Green functions (NEGF) has seen a rapid development over the recent three decades. Applications include diverse correlated many‐body systems in and out of equilibrium. Very good agreement with experiments and available exact theoretical results could be demonstrated if the proper selfenergy approximations were used. However, full two‐time NEGF simulations are computationally costly, as they suffer from a cubic scaling of the computation time with the simulation duration. Recently the G1–G2 scheme that exactly reformulates the generalized Kadanoff–Baym ansatz with Hartree–Fock propagators (HF‐GKBA) into time‐local equations is introduced, which achieves time‐linear scaling and allows for a dramatic speedup and extension of the simulations (Schluenzen et al. Phys. Rev. Lett. 2020, 124, 076601). Remarkably, this scaling is achieved quickly, and also for high‐level selfenergies, including the nonequilibrium GW and T‐matrix approximations (Joost et al. Phys. Rev. B 2020, 101, 245101). Even the dynamically screened ladder approximation is now feasible (Joost et al. Phys. Rev. B 2022, 105, 165155), and also applications to electron‐boson systems are demonstrated. Herein, an overview on recent results that are achieved with the G1–G2 scheme is presented. Problems and open questions are discussed and further ideas of how to overcome the current limitations of the scheme and present are presented. The G1–G2 scheme is illustrated by presenting applying it to the excitation dynamics of Hubbard clusters, to optical excitation of graphene, and to charge transfer during stopping of ions by correlated materials.

arXiv:2312.15030•physica status solidi (b)•Apr 25, 2024•15 citations

In- and out-of-equilibrium {\em ab initio} theory of electrons and phonons

Gianluca Stefanucci, Robert van Leeuwen, Enrico Perfetto

We lay down the {\em ab initio} many-body quantum theory of electrons and phonons in- and out-of-equilibrium at any temperature. We begin by addressing a fundamental issue concerning the {\em ab initio} Hamiltonian in the harmonic approximation, which we show must be determined {\em self-consistently} to avoid inconsistencies. After identifying the most suitable partitioning into a ``noninteracting'' and an ``interacting'' part we embark on the Green's function diagrammatic analysis. We single out key diagrammatic structures to carry on the expansion in terms of dressed propagators and screened interaction. The final outcome is the finite-temperature nonequilibrium extension of the Hedin equations, featuring the appearance of the time-local Ehrenfest diagram in the electronic self-energy. The Hedin equations have limited practical utility for real-time simulations of driven systems. We leverage the versatility of diagrammatic expansion to generate a closed system of integro-differential equations for the Green's functions and nuclear displacements. These are the Kadanoff-Baym equations for electrons and phonons. Another advantage of the diagrammatic derivation is the ability to use conserving approximations, which ensure the satisfaction of all fundamental conservation laws during the time evolution. As an example we show that the adiabatic Born-Oppenheimer approximation is not conserving whereas its dynamical extension is conserving provided that the electrons are treated in the Fan-Migdal approximation with a dynamically screened electron-phonon coupling. We also derive the formal solution of the Kadanoff-Baym equations in the long time limit and at the steady state. The expansion of the phononic Green's function around the quasi-phonon energies points to a possible correlation-induced splitting of the phonon dispersion in materials with no time-reversal invariance.

arXiv:2303.02102•Mar 3, 2023•1 citation

Time-linear scaling nonequilibrium Green's function method for real-time simulations of interacting electrons and bosons. II. Dynamics of polarons and doublons

Y. Pavlyukh, Enrico Perfetto, Daniel Karlsson +2 more

Nonequilibrium dynamics of the open chain Holstein-Hubbard model is studied using the linear time-scaling $\mathrm{GKBA}+\mathrm{ODE}$ scheme developed in Pavlyukh et al. [Phys. Rev. B 105, 125134 (2022)]. We focus on the set of parameters relevant for photovoltaic materials, i.e., a pair of electrons interacting with phonons at the crossover between the adiabatic and antiadiabatic regimes and at moderately large electron-electron interaction. By comparing with exact solutions for two corner cases, we demonstrate the accuracy of the $T$ matrix (in the $pp$ channel) and the second-order Fan ($GD$) approximations for the treatment of electronic ($e\text{\ensuremath{-}}e$) and electron-phonon ($e$-ph) correlations, respectively. The feedback of electrons on phonons is consistently included and is shown to be mandatory for the total energy conservation. When two interactions are simultaneously present, our simulations offer a glimpse into the dynamics of doublons and polarons unveiling the formation, propagation and decay of these quasiparticles, energy redistribution between them and self-trapping of electrons.

Physical review. B./Physical review. B•Mar 28, 2022•27 citations
This work contains a computationally efficient Green's function scheme for treating correlated electron-boson systems out of equilibrium. The numerical effort scales linearly with the simulation time, thus competing with that of the fastest quantum methods available. On a more fundamental level the proposed scheme fulfills energy balance as well as other fundamental conservation laws, and it can be systematically improved using standard many-body perturbation theory. We demonstrate its accuracy and stability through benchmarks against exact results in the Holstein model. However, our formulation equally applies to molecules and solids for first-principles transient calculations with quantum bosons, such as photons and phonons.
The interaction of electrons with quantized phonons and photons underlies the ultrafast dynamics of systems ranging from molecules to solids, and it gives rise to a plethora of physical phenomena experimentally accessible using time-resolved techniques. Green's function methods offer an invaluable interpretation tool since scattering mechanisms of growing complexity can be selectively incorporated in the theory. Currently, however, real-time Green's function simulations are either prohibitively expensive due to the cubic scaling with the propagation time or do neglect the feedback of electrons on the bosons, thus violating energy conservation. We put forward a computationally efficient Green's function scheme which overcomes both limitations. The numerical effort scales linearly with the propagation time while the simultaneous dressing of electrons and bosons guarantees the fulfillment of all fundamental conservation laws. We present a real-time study of the phonon-driven relaxation dynamics in an optically excited narrow band-gap insulator, highlighting the nonthermal behavior of the phononic degrees of freedom. Our formulation paves the way to first-principles simulations of electron-boson systems with unprecedented long propagation times.
Simulations of interacting electrons and bosons out of equilibrium, starting from first principles and aiming at realistic multiscale scenarios, is a grand theoretical challenge. Here, using the formalism of nonequilibrium Green's functions and relying in a crucial way on the recently discovered time-linear formulation of the Kadanoff-Baym equations, we present a versatile toolbox for the simulation of correlated electron-boson dynamics. A large class of methods are available, from the Ehrenfest to the dressed GDGDGD for the treatment of electron-boson interactions in combination with perturbative, i.e., Hartree-Fock and second-Born, or nonperturbative, i.e., GWGWGW and TTT matrices either without or with exchange effects, for the treatment of the Coulomb interaction. In all cases the numerical scaling is linear in time and the equations of motion satisfy all fundamental conservation laws.
We present a time-linear scaling method to simulate open and correlated quantum systems out of equilibrium. The method inherits from many-body perturbation theory the possibility to choose selectively the most relevant scattering processes in the dynamics, thereby paving the way to the real-time characterization of correlated ultrafast phenomena in quantum transport. The open system dynamics is described in terms of an "embedding correlator" from which the time-dependent current can be calculated using the Meir-Wingreen formula. We show how to efficiently implement our approach through a simple grafting into recently proposed time-linear Green's function methods for closed systems. Electron-electron and electron-phonon interactions can be treated on equal footing while preserving all fundamental conservation laws.
The theory of nonequilibrium Green functions (NEGF) has seen a rapid development over the recent three decades. Applications include diverse correlated many‐body systems in and out of equilibrium. Very good agreement with experiments and available exact theoretical results could be demonstrated if the proper selfenergy approximations were used. However, full two‐time NEGF simulations are computationally costly, as they suffer from a cubic scaling of the computation time with the simulation duration. Recently the G1–G2 scheme that exactly reformulates the generalized Kadanoff–Baym ansatz with Hartree–Fock propagators (HF‐GKBA) into time‐local equations is introduced, which achieves time‐linear scaling and allows for a dramatic speedup and extension of the simulations (Schluenzen et al. Phys. Rev. Lett. 2020, 124, 076601). Remarkably, this scaling is achieved quickly, and also for high‐level selfenergies, including the nonequilibrium GW and T‐matrix approximations (Joost et al. Phys. Rev. B 2020, 101, 245101). Even the dynamically screened ladder approximation is now feasible (Joost et al. Phys. Rev. B 2022, 105, 165155), and also applications to electron‐boson systems are demonstrated. Herein, an overview on recent results that are achieved with the G1–G2 scheme is presented. Problems and open questions are discussed and further ideas of how to overcome the current limitations of the scheme and present are presented. The G1–G2 scheme is illustrated by presenting applying it to the excitation dynamics of Hubbard clusters, to optical excitation of graphene, and to charge transfer during stopping of ions by correlated materials.
We lay down the {\em ab initio} many-body quantum theory of electrons and phonons in- and out-of-equilibrium at any temperature. We begin by addressing a fundamental issue concerning the {\em ab initio} Hamiltonian in the harmonic approximation, which we show must be determined {\em self-consistently} to avoid inconsistencies. After identifying the most suitable partitioning into a ``noninteracting'' and an ``interacting'' part we embark on the Green's function diagrammatic analysis. We single out key diagrammatic structures to carry on the expansion in terms of dressed propagators and screened interaction. The final outcome is the finite-temperature nonequilibrium extension of the Hedin equations, featuring the appearance of the time-local Ehrenfest diagram in the electronic self-energy. The Hedin equations have limited practical utility for real-time simulations of driven systems. We leverage the versatility of diagrammatic expansion to generate a closed system of integro-differential equations for the Green's functions and nuclear displacements. These are the Kadanoff-Baym equations for electrons and phonons. Another advantage of the diagrammatic derivation is the ability to use conserving approximations, which ensure the satisfaction of all fundamental conservation laws during the time evolution. As an example we show that the adiabatic Born-Oppenheimer approximation is not conserving whereas its dynamical extension is conserving provided that the electrons are treated in the Fan-Migdal approximation with a dynamically screened electron-phonon coupling. We also derive the formal solution of the Kadanoff-Baym equations in the long time limit and at the steady state. The expansion of the phononic Green's function around the quasi-phonon energies points to a possible correlation-induced splitting of the phonon dispersion in materials with no time-reversal invariance.
Nonequilibrium dynamics of the open chain Holstein-Hubbard model is studied using the linear time-scaling GKBA+ODE\mathrm{GKBA}+\mathrm{ODE}GKBA+ODE scheme developed in Pavlyukh et al. [Phys. Rev. B 105, 125134 (2022)]. We focus on the set of parameters relevant for photovoltaic materials, i.e., a pair of electrons interacting with phonons at the crossover between the adiabatic and antiadiabatic regimes and at moderately large electron-electron interaction. By comparing with exact solutions for two corner cases, we demonstrate the accuracy of the TTT matrix (in the pppppp channel) and the second-order Fan (GDGDGD) approximations for the treatment of electronic ($e\text{\ensuremath{-}}e$) and electron-phonon (eee-ph) correlations, respectively. The feedback of electrons on phonons is consistently included and is shown to be mandatory for the total energy conservation. When two interactions are simultaneously present, our simulations offer a glimpse into the dynamics of doublons and polarons unveiling the formation, propagation and decay of these quasiparticles, energy redistribution between them and self-trapping of electrons.