Partial self-energy parts of Kadanoff-Baym

Takeo Nishigori

Michael Potthoff

Segura Oscar Moreno, Y. Pavlyukh, Riku Tuovinen

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.

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.

Akihiro Nishiyama

We present a theoretical study of thermalization of quantum fields with Kadanoff-Baym (KB) equation. First we introduce a field theoretical technique known as 2 Particle Irreducible (2PI) effective action and review derivation of equation of motion for mean fields and fluctuations, then derive Kadanoff-Baym equation in scalar and gauge theories. In order to analyze thermalization processes we introduce relativistic entropy current based on the first order gradient expansion of Kadanoff-Baym equation. Then we show that in taking into account next-to-leading order (NLO) self-energy of the coupling expansion in φ4 theory (in symmetric phase <φ^> = 0), it is possible to prove the H-theorem within the 1st order gradient expansion. We also show that it is possible to prove the H-theorem for KB equation with NLO self-energy of 1/N expansion in O(N) theory (in symmetric phase). Furthermore we suggest a possibility for H-theorem to be satisfied in KB equation with leading-order (LO) self-energy of the coupling expansion in non-Abelian gauge theory. In d+1 dimensions (d＞2) LO particle number changing processes, such as 0-to-3, 1-to-2 and 2-to-1, which are prohibited in the normal Boltzmann approach with on-shell particles, may contribute to entropy production in KB approach although we must deal with gauge dependence and infrared divergence to complete the proof. Next we perform numerical simulations in spatially homogeneous configurations to investigate thermalization properties of the system in 1+1 and 2+1 dimensions for φ4 theory with NLO self-energy of the coupling expansion and 1+1 dimensions for O(N) theory with NLO self-energy of 1/N expansion. We also estimate the time evolution of the kinetic entropy in the case of scalar φ4 and O(N) theory. Then we find that at later times X0≫1/m, where m represents mass of particles, the kinetic entropy increases monotonically and approaches the equilibrium value, although the limited time interval at the earlier times X0～1/m invalidates the use of it. No thermalization occurs in the simulations with the normal Boltzmann equation with 2-to-2 collision term in 1+1 dimensions, while it occurs in the simulations with KB equation in 1+1 dimensions. It is due to two types of offshell effects: memory effects and spectral functions with the decay width.