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

Stability Bounds for the Generalized Kadanoff‐Baym Ansatz in the Holstein Dimer

Segura Oscar Moreno, Y. Pavlyukh, Riku Tuovinen

Published Jan 1, 2026• 0 citations
Publisher

Abstract

Predicting real-time dynamics in correlated systems is demanding: exact two-time Green's function methods are accurate but often too costly, while the Generalized Kadanoff-Baym Ansatz (GKBA) offers time-linear propagation at the risk of uncontrolled behavior. We examine when and why GKBA fails in a minimal yet informative setting, the Holstein dimer that describes electron–phonon coupling. Using a conserving, fully self-consistent electron–phonon self-energy, we map out parameter regions where GKBA dynamics is stable and where it becomes unstable. We trace the onset of these failures to qualitative changes in the model's ground-state solutions obtained from the full nonequilibrium Green's function theory, thereby providing practical stability bounds for GKBA time evolution. We further show that coupling the dimer to electronic leads can damp and, in part, cure these instabilities. The results supply simple diagnostics and guidelines for reliable GKBA simulations of electron–phonon dynamics.

DOI
10.1002/ctpp.70096&quot

Related papers

Stability Bounds for the Generalized Kadanoff‐Baym Ansatz in the Holstein Dimer

Oscar Moreno Segura, Y. Pavlyukh, Riku Tuovinen

ABSTRACT Predicting real‐time dynamics in correlated systems is demanding: exact two‐time Green's function methods are accurate but often too costly, while the Generalized Kadanoff‐Baym Ansatz (GKBA) offers time‐linear propagation at the risk of uncontrolled behavior. We examine when and why GKBA fails in a minimal yet informative setting, the Holstein dimer that describes electron–phonon coupling. Using a conserving, fully self‐consistent electron–phonon self‐energy, we map out parameter regions where GKBA dynamics is stable and where it becomes unstable. We trace the onset of these failures to qualitative changes in the model's ground‐state solutions obtained from the full nonequilibrium Green's function theory, thereby providing practical stability bounds for GKBA time evolution. We further show that coupling the dimer to electronic leads can damp and, in part, cure these instabilities. The results supply simple diagnostics and guidelines for reliable GKBA simulations of electron–phonon dynamics.

arXiv:2511.15582•Contributions to Plasma Physics•Mar 10, 2026•1 citation

Many-body Green’s function theory for electron-phonon interactions: The Kadanoff-Baym approach to spectral properties of the Holstein dimer

Niko Säkkinen, Yang Peng, Heiko Appel +1 more

We present a Kadanoff-Baym formalism to study time-dependent phenomena for systems of interacting electrons and phonons in the framework of many-body perturbation theory. The formalism takes correctly into account effects of the initial preparation of an equilibrium state and allows for an explicit time-dependence of both the electronic and phononic degrees of freedom. The method is applied to investigate the charge neutral and non-neutral excitation spectra of a homogeneous, two-site, two-electron Holstein model. This is an extension of a previous study of the ground state properties in the Hartree (H), partially self-consistent Born (Gd) and fully self-consistent Born (GD) approximations published in Säkkinen et al. [J. Chem. Phys. 143, 234101 (2015)]. Here, the homogeneous ground state solution is shown to become unstable for a sufficiently strong interaction while a symmetry-broken ground state solution is shown to be stable in the Hartree approximation. Signatures of this instability are observed for the partially self-consistent Born approximation but are not found for the fully self-consistent Born approximation. By understanding the stability properties, we are able to study the linear response regime by calculating the density-density response function by time-propagation. This amounts to a solution of the Bethe-Salpeter equation with a sophisticated kernel. The results indicate that none of the approximations is able to describe the response function during or beyond the bipolaronic crossover for the parameters investigated. Overall, we provide an extensive discussion on when the approximations are valid and how they fail to describe the studied exact properties of the chosen model system.

arXiv:1507.04726•The Journal of Chemical Physics•Dec 15, 2015•36 citations

Nonequilibrium Dynamics of the Hubbard Dimer

Y. Pavlyukh

Electron dynamics in a two‐sites Hubbard model is studied using the nonequilibrium Green's function approach. The study is motivated by the empirical observation that a full solution of the integro‐differential Kadanoff–Baym equation (KBE) is more stable and often accompanied by artificial damping [Marc Puig von Friesen, C. Verdozzi, and C.‐O. Almbladh (2009)] than its time‐linear reformulations relying on the generalized Kadanoff–Baym ansatz (GKBA). Additionally, for conserving theories, numerical simulations suggest that KBE produces natural occupations bounded by one and zero in agreement with the Pauli exclusion principle, whereas, in some regimes, GKBA‐based theories violate this principle. As the first step for understanding these issues, the electron dynamics arising in the adiabatic switching scenario is studied. Many‐body approximations are classified according to the channel of the Bethe–Salpeter equation in which electronic correlations are explicitly treated. They give rise to the so‐called second Born, T‐matrix, and GW approximations. In each of these cases, the model is reduced to a system of ordinary differential equations, which resemble equations of motion for a driven harmonic oscillator with time‐dependent frequencies. A more complete treatment of electronic correlations is achieved by combining different correlation channels, with parquet theory serving as a starting point.

physica status solidi (b)•Feb 14, 2024•6 citations

Nonequilibrium Dynamics of the Hubbard Dimer

Y. Pavlyukh

Electron dynamics in a two-sites Hubbard model is studied using the nonequilibrium Green's function approach. The study is motivated by the empirical observation that a full solution of the integro-differential Kadanoff-Baym equation (KBE) is more stable and often accompanied by artificial damping [Marc Puig von Friesen, C. Verdozzi, and C.-O. Almbladh (2009)] than its time-linear reformulations relying on the generalized Kadanoff-Baym ansatz (GKBA). Additionally, for conserving theories, numerical simulations suggest that KBE produces natural occupations bounded by one and zero in agreement with the Pauli exclusion principle, whereas, in some regimes, GKBA-based theories violate this principle. As the first step for understanding these issues, the electron dynamics arising in the adiabatic switching scenario is studied. Many-body approximations are classified according to the channel of the Bethe-Salpeter equation in which electronic correlations are explicitly treated. They give rise to the so-called second Born, $T$-matrix and $GW$ approximations. In each of these cases, the model is reduced to a system of ordinary differential equations, which resemble equations of motion for a driven harmonic oscillator with time-dependent frequencies. A more complete treatment of electronic correlations is achieved by combining different correlation channels, with parquet theory serving as a starting point.

arXiv:2401.11301•Jan 20, 2024•0 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

The non-equilibrium Green function approach to inhomogeneous quantum many-body systems using the generalized Kadanoff–Baym ansatz

S Hermanns, K Balzer, M. Bönitz

In the non-equilibrium Green function calculations, the use of the generalized Kadanoff–Baym ansatz (GKBA) allows for a simple approximate reconstruction of the two-time Green function from its time-diagonal value. With this, a drastic reduction of the computational needs is achieved in time-dependent calculations, making longer time propagation possible and more complex systems accessible. This paper gives credit to the GKBA that was introduced 25 years ago. After a detailed derivation of the GKBA, we recall its application to homogeneous systems and show how to extend it to strongly correlated, inhomogeneous systems. As a proof of concept, we present the results for a two-electron quantum well, where the correct treatment of the correlated electron dynamics is crucial for a correct description of the equilibrium and dynamic properties.

arXiv:1205.4427•Physica Scripta•Nov 1, 2012•45 citations
Predicting real-time dynamics in correlated systems is demanding: exact two-time Green's function methods are accurate but often too costly, while the Generalized Kadanoff-Baym Ansatz (GKBA) offers time-linear propagation at the risk of uncontrolled behavior. We examine when and why GKBA fails in a minimal yet informative setting, the Holstein dimer that describes electron–phonon coupling. Using a conserving, fully self-consistent electron–phonon self-energy, we map out parameter regions where GKBA dynamics is stable and where it becomes unstable. We trace the onset of these failures to qualitative changes in the model's ground-state solutions obtained from the full nonequilibrium Green's function theory, thereby providing practical stability bounds for GKBA time evolution. We further show that coupling the dimer to electronic leads can damp and, in part, cure these instabilities. The results supply simple diagnostics and guidelines for reliable GKBA simulations of electron–phonon dynamics.
ABSTRACT Predicting real‐time dynamics in correlated systems is demanding: exact two‐time Green's function methods are accurate but often too costly, while the Generalized Kadanoff‐Baym Ansatz (GKBA) offers time‐linear propagation at the risk of uncontrolled behavior. We examine when and why GKBA fails in a minimal yet informative setting, the Holstein dimer that describes electron–phonon coupling. Using a conserving, fully self‐consistent electron–phonon self‐energy, we map out parameter regions where GKBA dynamics is stable and where it becomes unstable. We trace the onset of these failures to qualitative changes in the model's ground‐state solutions obtained from the full nonequilibrium Green's function theory, thereby providing practical stability bounds for GKBA time evolution. We further show that coupling the dimer to electronic leads can damp and, in part, cure these instabilities. The results supply simple diagnostics and guidelines for reliable GKBA simulations of electron–phonon dynamics.
We present a Kadanoff-Baym formalism to study time-dependent phenomena for systems of interacting electrons and phonons in the framework of many-body perturbation theory. The formalism takes correctly into account effects of the initial preparation of an equilibrium state and allows for an explicit time-dependence of both the electronic and phononic degrees of freedom. The method is applied to investigate the charge neutral and non-neutral excitation spectra of a homogeneous, two-site, two-electron Holstein model. This is an extension of a previous study of the ground state properties in the Hartree (H), partially self-consistent Born (Gd) and fully self-consistent Born (GD) approximations published in Säkkinen et al. [J. Chem. Phys. 143, 234101 (2015)]. Here, the homogeneous ground state solution is shown to become unstable for a sufficiently strong interaction while a symmetry-broken ground state solution is shown to be stable in the Hartree approximation. Signatures of this instability are observed for the partially self-consistent Born approximation but are not found for the fully self-consistent Born approximation. By understanding the stability properties, we are able to study the linear response regime by calculating the density-density response function by time-propagation. This amounts to a solution of the Bethe-Salpeter equation with a sophisticated kernel. The results indicate that none of the approximations is able to describe the response function during or beyond the bipolaronic crossover for the parameters investigated. Overall, we provide an extensive discussion on when the approximations are valid and how they fail to describe the studied exact properties of the chosen model system.
Electron dynamics in a two‐sites Hubbard model is studied using the nonequilibrium Green's function approach. The study is motivated by the empirical observation that a full solution of the integro‐differential Kadanoff–Baym equation (KBE) is more stable and often accompanied by artificial damping [Marc Puig von Friesen, C. Verdozzi, and C.‐O. Almbladh (2009)] than its time‐linear reformulations relying on the generalized Kadanoff–Baym ansatz (GKBA). Additionally, for conserving theories, numerical simulations suggest that KBE produces natural occupations bounded by one and zero in agreement with the Pauli exclusion principle, whereas, in some regimes, GKBA‐based theories violate this principle. As the first step for understanding these issues, the electron dynamics arising in the adiabatic switching scenario is studied. Many‐body approximations are classified according to the channel of the Bethe–Salpeter equation in which electronic correlations are explicitly treated. They give rise to the so‐called second Born, T‐matrix, and GW approximations. In each of these cases, the model is reduced to a system of ordinary differential equations, which resemble equations of motion for a driven harmonic oscillator with time‐dependent frequencies. A more complete treatment of electronic correlations is achieved by combining different correlation channels, with parquet theory serving as a starting point.
Electron dynamics in a two-sites Hubbard model is studied using the nonequilibrium Green's function approach. The study is motivated by the empirical observation that a full solution of the integro-differential Kadanoff-Baym equation (KBE) is more stable and often accompanied by artificial damping [Marc Puig von Friesen, C. Verdozzi, and C.-O. Almbladh (2009)] than its time-linear reformulations relying on the generalized Kadanoff-Baym ansatz (GKBA). Additionally, for conserving theories, numerical simulations suggest that KBE produces natural occupations bounded by one and zero in agreement with the Pauli exclusion principle, whereas, in some regimes, GKBA-based theories violate this principle. As the first step for understanding these issues, the electron dynamics arising in the adiabatic switching scenario is studied. Many-body approximations are classified according to the channel of the Bethe-Salpeter equation in which electronic correlations are explicitly treated. They give rise to the so-called second Born, TTT-matrix and GWGWGW approximations. In each of these cases, the model is reduced to a system of ordinary differential equations, which resemble equations of motion for a driven harmonic oscillator with time-dependent frequencies. A more complete treatment of electronic correlations is achieved by combining different correlation channels, with parquet theory serving as a starting point.
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.
In the non-equilibrium Green function calculations, the use of the generalized Kadanoff–Baym ansatz (GKBA) allows for a simple approximate reconstruction of the two-time Green function from its time-diagonal value. With this, a drastic reduction of the computational needs is achieved in time-dependent calculations, making longer time propagation possible and more complex systems accessible. This paper gives credit to the GKBA that was introduced 25 years ago. After a detailed derivation of the GKBA, we recall its application to homogeneous systems and show how to extend it to strongly correlated, inhomogeneous systems. As a proof of concept, we present the results for a two-electron quantum well, where the correct treatment of the correlated electron dynamics is crucial for a correct description of the equilibrium and dynamic properties.