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:1609.04151 · Physical review. B./Physical review. B

Post-quench dynamics and suppression of thermalization in an open half-filled Hubbard layer

Igor V. Blinov, Pedro Ribeiro, Alexey N. Rubtsov

Published Jan 20, 2017• 4 citations
arXivHTMLPDF

Abstract

We study the time evolution of a half-filled Hubbard layer coupled to a magnon bath after a quench of the Hubbard interaction. Qualitatively different regimes, regarding the asymptotic long-time dynamics, are identified and characterized within the mean-field approximation. In the absence of the bath, the dynamics of the closed system is similar to that of a quenched BCS condensate. Though the presence of the bath introduces an additional relaxation mechanism, our numerical results and analytical arguments show that equilibration with the bath is not necessarily attained within the approximations used. Instead, nonequilibrium states, similar to the ones observed in the closed system, can emerge at long times as a consequence of the competition between the intrinsic relaxation mechanism (Landau damping, for example), and the bath-induced dissipation.

DOI
10.1103/physrevb.95.024309

Related papers

Transient phases and dynamical transitions in the post-quench evolution of the generalized Bose-Anderson model

Dmitry V. Chichinadze, Pedro Ribeiro, Yulia E. Shchadilova +1 more

The exact description of the time evolution of open correlated quantum systems remains one of the major challenges of condensed matter theory, especially for asymptotic long times where most numerical methods fail. Here, the post-quench dynamics of the $N$-component Bose-Anderson impurity model is studied in the $N\ensuremath{\rightarrow}\ensuremath{\infty}$ limit. The equilibrium phase diagram is similar to that of the Bose-Hubbard model in that it contains local versions of Mott and Bose phases. Using a numerically exact procedure, we are able to study the real-time evolution including asymptotic long-time regimes. The formation of long-lived transient phases is observed for quench paths crossing foreign phases. For quenches inside the local Bose condensed phase, a dynamical phase transition is reported that separates the evolution towards a new equilibrium state and a regime characterized at large times by a persistent phase rotation of the order parameter. We explain how such nondecaying modes can exist in the presence of a dissipative bath. We discuss the extension of our results to the experimental relevant finite-$N$ case and their implication for the existence of nondecaying modes in generic quantum systems in the presence of a bath.

arXiv:1604.01995•Physical review. B./Physical review. B•Aug 1, 2016•7 citations

Dynamical transition in interaction quenches of the one-dimensional Hubbard model

Simone A. Hamerla, Götz Silvester Uhrig

We show that the nonequilibrium time evolution after interaction quenches in the one-dimensional, integrable Hubbard model exhibits a dynamical transition in the half-filled case. This transition ceases to exist upon doping. Our study is based on systematically extended equations of motion. Thus it is controlled for small and moderate times; no relaxation effects are neglected. Remarkable similarities to the quench dynamics in the infinite-dimensional Hubbard model are found, suggesting that dynamical transitions are a general feature of quenches in such models.

arXiv:1302.4109•Physical Review B•Feb 28, 2013•33 citations

Thermalization Dynamics in the Two-Dimensional Hubbard Model with Neural-Network Quantum States

Alessandro Sinibaldi, Luciano Loris Viteritti, Riccardo Rende +2 more

Thermalization in strongly correlated fermionic systems remains a central open problem in quantum many-body physics. In this work, we investigate the real-time dynamics and the approach to thermalization in the two-dimensional Hubbard model, a paradigmatic framework for correlated electrons, relevant to high-temperature superconductivity and ultracold quantum simulation. Focusing on the half-filled square lattice, we monitor the time evolution of the double occupancy following a quench in the on-site interaction $U$, and assess whether its long-time value is captured by a canonical thermal ensemble. We employ time-dependent variational Monte Carlo methods combined with transformer-based Neural-Network Quantum States to accurately describe the nonequilibrium dynamics of fermions, especially for the behavior at long times, thereby accessing regimes that were previously inaccessible to numerical simulations. Our results reveal two dynamical behaviors: for weak to intermediate interactions, the double occupancy rapidly approaches the thermal prediction, consistent with ergodic evolution; beyond a critical interaction $U_{C}$, the dynamics remains distinct from the thermal expectation on the timescales investigated, revealing signatures of a prethermal plateau delaying fast relaxation. These results establish numerical simulation as a powerful tool to probe nonequilibrium quantum phenomena in correlated fermionic matter.

cond-mat.str-elCross-listedarXiv:2606.05293•Jun 5, 2026•0 citations

Thermalization after an Interaction Quench in the Hubbard Model

Martin Eckstein, Marcus Kollar, Philipp Werner

We use nonequilibrium dynamical mean-field theory to study the time evolution of the fermionic Hubbard model after an interaction quench. Both in the weak-coupling and in the strong-coupling regime the system is trapped in quasistationary states on intermediate time scales. These two regimes are separated by a sharp crossover at U(c)dyn=0.8 in units of the bandwidth, where fast thermalization occurs. Our results indicate a dynamical phase transition which should be observable in experiments on trapped fermionic atoms.

arXiv:0904.0976•Physical Review Letters•Jul 29, 2009•460 citations

Time and spatially resolved quench of the fermionic Hubbard model showing restricted equilibration

Florian Goth, Fakher F. Assaad

We investigate the quench of half-filled one-dimensional and two-dimensional fermionic Hubbard models to models without Coulomb interaction. Since the time propagation is Gaussian, we can use a variety of time-dependent quantum Monte Carlo methods to tackle this problem without generating a dynamical sign problem. Using a continuous time quantum Monte Carlo method, we achieve a system size of 128 sites in one dimension, and using a Blankenbecler-Scalapino-Sugar type algorithm, we were able to simulate $20\ifmmode\times\else\texttimes\fi{}20$ square lattices. Applying these methods to study the dynamics after the quench, we observe that the final state of the system can be reasonably well described by a thermal single-particle density matrix that takes the initial single-particle conservation laws into account. The characteristic decay toward this limit is found to be oscillatory with an additional power law decay that depends on the dimensionality. This numerically exact result is shown to compare favorably to mean-field approximations as well as to perturbation theory. Furthermore, we observe the information propagation in the one-dimensional case in the charge-charge and spin-spin correlations and find that it is linear with a velocity of roughly $v\ensuremath{\approx}4$ in units of the hopping amplitude.

arXiv:1108.2703•Physical Review B•Feb 27, 2012•29 citations

Nonequilibrium dynamics from BCS to the bosonic limit

G. Seibold, J. Lorenzana

Within the time-dependent Gutzwiller approximation applied to the negative-$U$ Hubbard model, we investigate the dynamics of a superconductor after an interaction quench for different values of the final interaction from weak to strong coupling. The equilibrium BCS-BEC crossover becomes a sharp transition in the out-of-equilibrium dynamics between weak- and strong-coupling dynamical phases. Two different frequencies (${\mathrm{\ensuremath{\Omega}}}_{J}<{\mathrm{\ensuremath{\Omega}}}_{U}$) dominate the order parameter dynamics. In the weak coupling phase, ${\mathrm{\ensuremath{\Omega}}}_{J}$ follows twice the asymptotic average value of the gap (${\mathrm{\ensuremath{\Omega}}}_{J}\ensuremath{\approx}2{\mathrm{\ensuremath{\Delta}}}_{\ensuremath{\infty}}$) and ${\mathrm{\ensuremath{\Omega}}}_{U}$ is much larger but hardy visible in the anomalous density. At long times the BCS dynamics is recovered but surprisingly differences remain during the transient phase. In the strong coupling phase, the dynamics decouples from the asymptotic value of the order parameter except at exactly half filling and close to the dynamical transition where ${\mathrm{\ensuremath{\Omega}}}_{U}\ensuremath{\approx}2{\mathrm{\ensuremath{\Delta}}}_{\ensuremath{\infty}}$. The out-of-equilibrium transient spectral density and optical conductivity are presented and discussed in relation to pump probe experiments. Both ${\mathrm{\ensuremath{\Omega}}}_{J}$ and ${\mathrm{\ensuremath{\Omega}}}_{U}$ give rise to a complex structure of self-driven slow Rabi oscillations which are visible in the nonequilibrium optical conductivity where also sidebands appear due to the modulation of the double occupancy by superconducting amplitude oscillations. Analogous results apply to CDW and SDW systems. Our results show that in systems with long coherence times, pump-probe experiments allow us to characterize the regime (BCS vs preformed pairs) through very specific out-of-equilibrium fingerprints.

arXiv:2004.02454•Physical review. B./Physical review. B•Oct 5, 2020•12 citations
We study the time evolution of a half-filled Hubbard layer coupled to a magnon bath after a quench of the Hubbard interaction. Qualitatively different regimes, regarding the asymptotic long-time dynamics, are identified and characterized within the mean-field approximation. In the absence of the bath, the dynamics of the closed system is similar to that of a quenched BCS condensate. Though the presence of the bath introduces an additional relaxation mechanism, our numerical results and analytical arguments show that equilibration with the bath is not necessarily attained within the approximations used. Instead, nonequilibrium states, similar to the ones observed in the closed system, can emerge at long times as a consequence of the competition between the intrinsic relaxation mechanism (Landau damping, for example), and the bath-induced dissipation.
The exact description of the time evolution of open correlated quantum systems remains one of the major challenges of condensed matter theory, especially for asymptotic long times where most numerical methods fail. Here, the post-quench dynamics of the NNN-component Bose-Anderson impurity model is studied in the $N\ensuremath{\rightarrow}\ensuremath{\infty}$ limit. The equilibrium phase diagram is similar to that of the Bose-Hubbard model in that it contains local versions of Mott and Bose phases. Using a numerically exact procedure, we are able to study the real-time evolution including asymptotic long-time regimes. The formation of long-lived transient phases is observed for quench paths crossing foreign phases. For quenches inside the local Bose condensed phase, a dynamical phase transition is reported that separates the evolution towards a new equilibrium state and a regime characterized at large times by a persistent phase rotation of the order parameter. We explain how such nondecaying modes can exist in the presence of a dissipative bath. We discuss the extension of our results to the experimental relevant finite-NNN case and their implication for the existence of nondecaying modes in generic quantum systems in the presence of a bath.
We show that the nonequilibrium time evolution after interaction quenches in the one-dimensional, integrable Hubbard model exhibits a dynamical transition in the half-filled case. This transition ceases to exist upon doping. Our study is based on systematically extended equations of motion. Thus it is controlled for small and moderate times; no relaxation effects are neglected. Remarkable similarities to the quench dynamics in the infinite-dimensional Hubbard model are found, suggesting that dynamical transitions are a general feature of quenches in such models.
Thermalization in strongly correlated fermionic systems remains a central open problem in quantum many-body physics. In this work, we investigate the real-time dynamics and the approach to thermalization in the two-dimensional Hubbard model, a paradigmatic framework for correlated electrons, relevant to high-temperature superconductivity and ultracold quantum simulation. Focusing on the half-filled square lattice, we monitor the time evolution of the double occupancy following a quench in the on-site interaction UUU, and assess whether its long-time value is captured by a canonical thermal ensemble. We employ time-dependent variational Monte Carlo methods combined with transformer-based Neural-Network Quantum States to accurately describe the nonequilibrium dynamics of fermions, especially for the behavior at long times, thereby accessing regimes that were previously inaccessible to numerical simulations. Our results reveal two dynamical behaviors: for weak to intermediate interactions, the double occupancy rapidly approaches the thermal prediction, consistent with ergodic evolution; beyond a critical interaction UCU_{C}UC​, the dynamics remains distinct from the thermal expectation on the timescales investigated, revealing signatures of a prethermal plateau delaying fast relaxation. These results establish numerical simulation as a powerful tool to probe nonequilibrium quantum phenomena in correlated fermionic matter.
We use nonequilibrium dynamical mean-field theory to study the time evolution of the fermionic Hubbard model after an interaction quench. Both in the weak-coupling and in the strong-coupling regime the system is trapped in quasistationary states on intermediate time scales. These two regimes are separated by a sharp crossover at U(c)dyn=0.8 in units of the bandwidth, where fast thermalization occurs. Our results indicate a dynamical phase transition which should be observable in experiments on trapped fermionic atoms.
We investigate the quench of half-filled one-dimensional and two-dimensional fermionic Hubbard models to models without Coulomb interaction. Since the time propagation is Gaussian, we can use a variety of time-dependent quantum Monte Carlo methods to tackle this problem without generating a dynamical sign problem. Using a continuous time quantum Monte Carlo method, we achieve a system size of 128 sites in one dimension, and using a Blankenbecler-Scalapino-Sugar type algorithm, we were able to simulate $20\ifmmode\times\else\texttimes\fi{}20$ square lattices. Applying these methods to study the dynamics after the quench, we observe that the final state of the system can be reasonably well described by a thermal single-particle density matrix that takes the initial single-particle conservation laws into account. The characteristic decay toward this limit is found to be oscillatory with an additional power law decay that depends on the dimensionality. This numerically exact result is shown to compare favorably to mean-field approximations as well as to perturbation theory. Furthermore, we observe the information propagation in the one-dimensional case in the charge-charge and spin-spin correlations and find that it is linear with a velocity of roughly $v\ensuremath{\approx}4$ in units of the hopping amplitude.
Within the time-dependent Gutzwiller approximation applied to the negative-UUU Hubbard model, we investigate the dynamics of a superconductor after an interaction quench for different values of the final interaction from weak to strong coupling. The equilibrium BCS-BEC crossover becomes a sharp transition in the out-of-equilibrium dynamics between weak- and strong-coupling dynamical phases. Two different frequencies (${\mathrm{\ensuremath{\Omega}}}_{J}<{\mathrm{\ensuremath{\Omega}}}_{U}$) dominate the order parameter dynamics. In the weak coupling phase, ${\mathrm{\ensuremath{\Omega}}}_{J}$ follows twice the asymptotic average value of the gap (${\mathrm{\ensuremath{\Omega}}}_{J}\ensuremath{\approx}2{\mathrm{\ensuremath{\Delta}}}_{\ensuremath{\infty}}$) and ${\mathrm{\ensuremath{\Omega}}}_{U}$ is much larger but hardy visible in the anomalous density. At long times the BCS dynamics is recovered but surprisingly differences remain during the transient phase. In the strong coupling phase, the dynamics decouples from the asymptotic value of the order parameter except at exactly half filling and close to the dynamical transition where ${\mathrm{\ensuremath{\Omega}}}_{U}\ensuremath{\approx}2{\mathrm{\ensuremath{\Delta}}}_{\ensuremath{\infty}}$. The out-of-equilibrium transient spectral density and optical conductivity are presented and discussed in relation to pump probe experiments. Both ${\mathrm{\ensuremath{\Omega}}}_{J}$ and ${\mathrm{\ensuremath{\Omega}}}_{U}$ give rise to a complex structure of self-driven slow Rabi oscillations which are visible in the nonequilibrium optical conductivity where also sidebands appear due to the modulation of the double occupancy by superconducting amplitude oscillations. Analogous results apply to CDW and SDW systems. Our results show that in systems with long coherence times, pump-probe experiments allow us to characterize the regime (BCS vs preformed pairs) through very specific out-of-equilibrium fingerprints.