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:1407.8288 · Physical Review B

Interaction quench in the Holstein model: Thermalization crossover from electron- to phonon-dominated relaxation

Yuta Murakami, Philipp Werner, Naoto Tsuji, Hideo Aoki

Published Jan 20, 2015• 70 citations
arXivHTMLPDF

Abstract

We study the relaxation of the Holstein model after a sudden switch-on of the interaction by means of the nonequilibrium dynamical mean field theory, with the self-consistent Migdal approximation as an impurity solver. We show that there exists a qualitative change in the thermalization dynamics as the interaction is varied in the weak-coupling regime. On the weaker interaction side of this crossover, the phonon oscillations are damped more rapidly than the electron thermalization time scale, as determined from the relaxation of the electron momentum distribution function. On the stronger interaction side, the relaxation of the electrons becomes faster than the phonon damping. In this regime, despite long-lived phonon oscillations, a thermalized momentum distribution is realized temporarily. The origin of the ``thermalization crossover'' found here is traced back to different behaviors of the electron and phonon self-energies as a function of the electron-phonon coupling. In addition, the importance of the phonon dynamics is demonstrated by comparing the self-consistent Migdal results with those obtained with a simpler Hartree-Fock impurity solver that neglects the phonon self-energy. The latter scheme does not properly describe the evolution and thermalization of isolated electron-phonon systems.

DOI
10.1103/physrevb.91.045128

Related papers

Thermalization Fronts in the Hubbard-Holstein Model

Antonio Picano, Marco Schirò

We investigate the nonequilibrium dynamics of the weak-coupling Hubbard-Holstein model after a sudden switch-on of the electron-phonon interaction within nonequilibrium dynamical mean-field theory (DMFT). Using the self-consistent Migdal approximation for the electron-phonon coupling together with second-order perturbation theory for the electron-electron interaction, we show that the relaxation dynamics exhibits a crossover between electron-dominated and phonon-dominated regimes, extending to finite Hubbard interaction the scenario previously identified in the Holstein model. To investigate the microscopic buildup of the thermal state, we analyze the dynamics within the Step-by-Step DMFT framework. In the plane of real time and DMFT iteration number, thermalization is marked by a sharp propagating front. This front appears in electronic observables already for weak quenches within the simulated time window, whereas the phononic sector exhibits a visible front only at sufficiently strong coupling. Thus, at weak coupling the local dispersionless phonons show a delayed onset of front formation, while near and beyond the crossover the front develops on comparable timescales in both the electronic and phononic sectors. Whenever both fronts are resolved, they propagate with the same velocity, showing that thermalization spreads coherently through the coupled electron-phonon system.

arXiv:2604.10775•Apr 12, 2026•0 citations

Phonon-enhanced relaxation and excitation in the Holstein-Hubbard model

Philipp Werner, Martin Eckstein

We study quenches of the interaction and electron-phonon coupling parameter in the Hubbard-Holstein model, using nonequilibrium dynamical mean field theory. The calculations are based on a generalized Lang-Firsov scheme for time-dependent interactions or externally driven phonons, and an approximate strong-coupling impurity solver. The interaction quench calculations reveal the phonon-assisted decay of excess doublons, while the quenches of the electron-phonon coupling lead to persistent oscillations of the phonons and to a phonon-enhanced doublon production.

arXiv:1307.8096•Physical Review B•Oct 7, 2013•65 citations

Efficient Dynamical Mean Field Simulation of the Holstein-Hubbard Model

Philipp Werner, Andrew J. Millis

We present a method for solving impurity models with electron-phonon coupling, which treats the phonons efficiently and without approximations. The algorithm is applied to the Holstein-Hubbard model in the dynamical mean field approximation, where it allows access to strong interactions, very low temperatures, and arbitrary fillings. We show that a renormalized Migdal-Eliashberg theory provides a reasonable description of the phonon contribution to the electronic self-energy in strongly doped systems, but fails if the quasiparticle energy becomes of order of the phonon frequency.

arXiv:cond-mat/0701730•Physical Review Letters•Oct 4, 2007•180 citations

Integrable-to-thermalizing crossover in nonequilibrium superconductors

Andrey Grankin, Victor Galitski

Motivated by the experiment by M. Budden {\em et al.} [Nature Physics {\bf 17}, 611 (2021)], who observed signatures of long-lived photo-induced superconductivity, we develop an accurate analytical/computational approach to non-equilibrium superconductivity following a quench. We consider the BCS-Holstein model, which includes both integrable local electron-electron interactions and integrability-breaking electron-phonon coupling. We develop Keldysh-Eliashberg theory on the Kadanoff-Baym contour, which enables us to describe non-equilibrium dynamics of the superconductor. We consider a quench in interactions, which results in a dynamic transition from the initial superconducting state to a normal thermal state in the end of the evolution. It is shown that the dynamics contain two stages: The early-time integrable behavior, involving coherent oscillations of the superconducting order parameter, crosses over to the late-time ergodic dynamics exhibiting a thermal decay into an equilibrium state. In the former regime, our computational approach both reproduces exact analytical results on the integrable dynamics of the order parameter and generalizes those to the case of an initial thermal state. The method also succeeds for the first time in describing both integrable-to-thermalizing crossover and the late-time thermal decay, which is shown to be consistent with the time-dependent Ginzburg-Landau theory (with the exponential decay time dependent on the density of quasiparticle excitations). We observe the electron distribution function approaching the Fermi-Dirac thermal distribution at final stages. The details of two-time non-equilibrium dynamics depend on the density of quasiparticles in the initial state and the integrability-breaking parameters, which under certain conditions may result in a long-lived transient superconductivity consistent with experiment.

arXiv:2312.13391•Physical review. B./Physical review. B•Dec 12, 2025•1 citation

Comparative study of nonequilibrium insulator-to-metal transitions in electron-phonon systems

Sharareh Sayyad, Rok Žitko, Hugo U. R. Strand +2 more

We study equilibrium and nonequilibrium properties of electron-phonon systems described by the Hubbard-Holstein model using dynamical mean-field theory. In equilibrium, we benchmark the results for impurity solvers based on the one-crossing approximation and slave-rotor approximation against non-perturbative numerical renormalization group reference data. We also examine how well the low-energy properties of the electron-boson coupled systems can be reproduced by an effective static electron-electron interaction. The one-crossing and slave-rotor approximations are then used to simulate insulator-to-metal transitions induced by a sudden switch-on of the electron-phonon interaction. The slave-rotor results suggest the existence of a critical electron-phonon coupling above which the system is transiently trapped in a nonthermal metallic state with coherent quasiparticles. The same quench protocol in the one-crossing approximation results in a bad metallic state.

arXiv:1810.01598•Physical review. B./Physical review. B•Jan 9, 2019•15 citations

Stochastic semiclassical theory for nonequilibrium electron-phonon coupled systems

Antonio Picano, Francesco Grandi, Philipp Werner +1 more

We discuss a semiclassical approach to solve the quantum impurity model within nonequilibrium dynamical mean-field theory for electron-lattice models. The effect of electronic fluctuations on the phonon is kept beyond Ehrenfest dynamics, leading to a stochastic phonon evolution with damping and noise terms that are self-consistently determined by the electronic correlation functions in the fluctuating phonon field. Together with a solution of the electronic model based on a nonperturbative quantum Boltzmann equation, the approach can be used to address the coupled dynamics of the electrons and the lattice during photoinduced phase transitions. Results for the Anderson-Holstein model are benchmarked against numerically exact quantum Monte Carlo data. We find good agreement for the phonon distribution function at temperatures comparable to the charge ordering temperature. The general formulation can be extended to models with electron-electron interactions or multiorbital systems.

arXiv:2209.00428•Physical review. B./Physical review. B•Jul 10, 2023•11 citations
We study the relaxation of the Holstein model after a sudden switch-on of the interaction by means of the nonequilibrium dynamical mean field theory, with the self-consistent Migdal approximation as an impurity solver. We show that there exists a qualitative change in the thermalization dynamics as the interaction is varied in the weak-coupling regime. On the weaker interaction side of this crossover, the phonon oscillations are damped more rapidly than the electron thermalization time scale, as determined from the relaxation of the electron momentum distribution function. On the stronger interaction side, the relaxation of the electrons becomes faster than the phonon damping. In this regime, despite long-lived phonon oscillations, a thermalized momentum distribution is realized temporarily. The origin of the ``thermalization crossover'' found here is traced back to different behaviors of the electron and phonon self-energies as a function of the electron-phonon coupling. In addition, the importance of the phonon dynamics is demonstrated by comparing the self-consistent Migdal results with those obtained with a simpler Hartree-Fock impurity solver that neglects the phonon self-energy. The latter scheme does not properly describe the evolution and thermalization of isolated electron-phonon systems.
We investigate the nonequilibrium dynamics of the weak-coupling Hubbard-Holstein model after a sudden switch-on of the electron-phonon interaction within nonequilibrium dynamical mean-field theory (DMFT). Using the self-consistent Migdal approximation for the electron-phonon coupling together with second-order perturbation theory for the electron-electron interaction, we show that the relaxation dynamics exhibits a crossover between electron-dominated and phonon-dominated regimes, extending to finite Hubbard interaction the scenario previously identified in the Holstein model. To investigate the microscopic buildup of the thermal state, we analyze the dynamics within the Step-by-Step DMFT framework. In the plane of real time and DMFT iteration number, thermalization is marked by a sharp propagating front. This front appears in electronic observables already for weak quenches within the simulated time window, whereas the phononic sector exhibits a visible front only at sufficiently strong coupling. Thus, at weak coupling the local dispersionless phonons show a delayed onset of front formation, while near and beyond the crossover the front develops on comparable timescales in both the electronic and phononic sectors. Whenever both fronts are resolved, they propagate with the same velocity, showing that thermalization spreads coherently through the coupled electron-phonon system.
We study quenches of the interaction and electron-phonon coupling parameter in the Hubbard-Holstein model, using nonequilibrium dynamical mean field theory. The calculations are based on a generalized Lang-Firsov scheme for time-dependent interactions or externally driven phonons, and an approximate strong-coupling impurity solver. The interaction quench calculations reveal the phonon-assisted decay of excess doublons, while the quenches of the electron-phonon coupling lead to persistent oscillations of the phonons and to a phonon-enhanced doublon production.
We present a method for solving impurity models with electron-phonon coupling, which treats the phonons efficiently and without approximations. The algorithm is applied to the Holstein-Hubbard model in the dynamical mean field approximation, where it allows access to strong interactions, very low temperatures, and arbitrary fillings. We show that a renormalized Migdal-Eliashberg theory provides a reasonable description of the phonon contribution to the electronic self-energy in strongly doped systems, but fails if the quasiparticle energy becomes of order of the phonon frequency.
Motivated by the experiment by M. Budden {\em et al.} [Nature Physics {\bf 17}, 611 (2021)], who observed signatures of long-lived photo-induced superconductivity, we develop an accurate analytical/computational approach to non-equilibrium superconductivity following a quench. We consider the BCS-Holstein model, which includes both integrable local electron-electron interactions and integrability-breaking electron-phonon coupling. We develop Keldysh-Eliashberg theory on the Kadanoff-Baym contour, which enables us to describe non-equilibrium dynamics of the superconductor. We consider a quench in interactions, which results in a dynamic transition from the initial superconducting state to a normal thermal state in the end of the evolution. It is shown that the dynamics contain two stages: The early-time integrable behavior, involving coherent oscillations of the superconducting order parameter, crosses over to the late-time ergodic dynamics exhibiting a thermal decay into an equilibrium state. In the former regime, our computational approach both reproduces exact analytical results on the integrable dynamics of the order parameter and generalizes those to the case of an initial thermal state. The method also succeeds for the first time in describing both integrable-to-thermalizing crossover and the late-time thermal decay, which is shown to be consistent with the time-dependent Ginzburg-Landau theory (with the exponential decay time dependent on the density of quasiparticle excitations). We observe the electron distribution function approaching the Fermi-Dirac thermal distribution at final stages. The details of two-time non-equilibrium dynamics depend on the density of quasiparticles in the initial state and the integrability-breaking parameters, which under certain conditions may result in a long-lived transient superconductivity consistent with experiment.
We study equilibrium and nonequilibrium properties of electron-phonon systems described by the Hubbard-Holstein model using dynamical mean-field theory. In equilibrium, we benchmark the results for impurity solvers based on the one-crossing approximation and slave-rotor approximation against non-perturbative numerical renormalization group reference data. We also examine how well the low-energy properties of the electron-boson coupled systems can be reproduced by an effective static electron-electron interaction. The one-crossing and slave-rotor approximations are then used to simulate insulator-to-metal transitions induced by a sudden switch-on of the electron-phonon interaction. The slave-rotor results suggest the existence of a critical electron-phonon coupling above which the system is transiently trapped in a nonthermal metallic state with coherent quasiparticles. The same quench protocol in the one-crossing approximation results in a bad metallic state.
We discuss a semiclassical approach to solve the quantum impurity model within nonequilibrium dynamical mean-field theory for electron-lattice models. The effect of electronic fluctuations on the phonon is kept beyond Ehrenfest dynamics, leading to a stochastic phonon evolution with damping and noise terms that are self-consistently determined by the electronic correlation functions in the fluctuating phonon field. Together with a solution of the electronic model based on a nonperturbative quantum Boltzmann equation, the approach can be used to address the coupled dynamics of the electrons and the lattice during photoinduced phase transitions. Results for the Anderson-Holstein model are benchmarked against numerically exact quantum Monte Carlo data. We find good agreement for the phonon distribution function at temperatures comparable to the charge ordering temperature. The general formulation can be extended to models with electron-electron interactions or multiorbital systems.