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

Nonequilibrium dynamics of superconductivity in the attractive Hubbard model

Gia-Wei Chern, Kipton M. Barros

Published Jan 30, 2019• 11 citations
arXivHTMLPDF

Abstract

We present a framework of semiclassical superconductivity (SC) dynamics that properly includes effects of spatial fluctuations for the attractive Hubbard model. We consider both coherent and adiabatic limits. To model the coherent SC dynamics, we develop a real-space von Neumann equation based on the time-dependent Hartree-Fock-Bogoliubov theory. Applying our method to interaction quenches in the negative-$U$ Hubbard model, we show that the relaxation of SC order at weak coupling is dominated by Landau-damping. At strong coupling, we find a two-stage relaxation of the pairing field: a collapse of the synchronized oscillation of Cooper pairs due to spatial inhomogeneity, followed by a slow relaxation to a quasistationary state. SC dynamics in the adiabatic limit is described by a quantum Landau-Lifshitz equation with Ginzburg-Landau relaxation. Numerical simulations of the pump-probe process show that long-time recovery of the pairing field is dominated by defect dynamics. Our results demonstrate the important role of spatial fluctuations in both limits.

DOI
10.1103/physrevb.99.035162

Related papers

Dynamic inhomogeneities and phase separation after quantum quenches in strongly correlated systems

Gia-Wei Chern

We present a Gutzwiller von Neumann dynamics (GvND) method for simulating equilibrium and nonequilibrium phenomena in strongly correlated electron systems. Our approach is a real-space formulation of the time-dependent Gutzwiller approximation method. Applying the GvND method to simulate interaction quenches in the Peierls-Hubbard model, we demonstrate the amplification of initial inhomogeneities, which in turn results in the collapse of quench-induced synchronized oscillation. Moreover, we find a dynamical phase transition separating two quasi-stationary regimes with rather distinct spatial distributions of physical quantities after the collapsed oscillation. In particular, in the strong-coupling regime, the system exhibits a dynamic phase separation in the quasi-stationary state. Our results thus underscore the importance of spatial fluctuations in the nonequilibrium dynamics of strongly correlated systems.

arXiv:1711.07972•Nov 21, 2017•0 citations

Spatial instability of a homogeneous superconductor after an interaction quench

Medha Roy Choudhury, Sankha Subhra Bakshi, Pinaki Majumdar +1 more

We investigate the quantum dynamics of a conventional superconductor, described by the attrac- tive Hubbard model, driven far from equilibrium by quenching the Hubbard attraction. We evolve the full density matrix of the system in time using the Heisenberg equation of motion, tracking the spatiotemporal dynamics of the pairing and density fields. We use very weakly inhomogeneous initial conditions to probe the stability of the homogeneous dynamics to spatial perturbations. At initial times after the quench the pairing and density fields remain homogeneous in space while fluctuating in time. However, depending on the quench strength the spatially uniform dynamics can become unstable beyond a threshold time, with emergence of spatial textures. The onset time for the instability depends on the strength of the quench and the initial inhomogeneity. At long times the system evolves towards a state with significantly reduced spatial heterogeneity. The time dependence of the global superconducting order parameter suggests two broad categories of re- sponse to a quench. One of these involves strongly inhomogeneous pairing and bond currents in the nonequilibrium superconductor.

arXiv:2609.13426•Sep 11, 2026•0 citations

Electronic and fluctuation dynamics following a quench to the superconducting phase

Christopher Stahl, Martin Eckstein

We investigate the dynamics of superconducting fluctuations in the attractive three-dimensional Hubbard model after a quench from the disordered phase to the ordered regime. While the long-time evolution is well understood in terms of dissipative time-dependent Ginzburg-Landau models with unstable potentials, early times are more demanding due to the inseparable dynamics of the pairing fluctuations and the electronic quasiparticles. Our simulation using the time-dependent fluctuation exchange approximation treats both degrees of freedom on the same footing and reveals a nonthermal electronic regime causing a nonmonotonous growth of the fluctuations. This feature is not directly captured by the Ginzburg-Landau theory but nevertheless remains observable beyond the thermalization time of the electrons. We further explore how the growth of the order parameter fluctuations leads to an opening of a pseudogap in the electronic spectrum and identify Andreev reflections as the dominant mechanism behind the gap opening.

arXiv:2007.10305•Physical review. B./Physical review. B•Jan 13, 2021•17 citations

Charge and pairing dynamics in the attractive Hubbard model: Mode coupling and the validity of linear-response theory

Jörg Bünemann, G. Seibold

Pump-probe experiments have turned out as a powerful tool in order to study the dynamics of competing orders in a large variety of materials. The corresponding analysis of the data often relies on standard linear-response theory generalized to nonequilibrium situations. Here we examine the validity of such an approach for the charge and pairing response of systems with charge-density wave and (or) superconducting (SC) order. Our investigations are based on the attractive Hubbard model which we study within the time-dependent Hartree-Fock approximation. In particular, we calculate the quench and pump-probe dynamics for SC and charge order parameters in order to analyze the frequency spectra and the coupling of the probe field to the specific excitations. Our calculations reveal that the ``linear-response assumption'' is justified for small to moderate nonequilibrium situations (i.e., pump pulses) in the case of a purely charge-ordered ground state. However, the pump-probe dynamics on top of a superconducting ground state is determined by phase and amplitude modes which get coupled far from the equilibrium state indicating the failure of the linear-response assumption.

arXiv:1708.05984•Physical review. B./Physical review. B•Dec 26, 2017•7 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

From sudden quench to adiabatic dynamics in the attractive Hubbard model

Giacomo Mazza

We study the crossover between the sudden quench limit and the adiabatic dynamics of superconducting states in the attractive Hubbard model. We focus on the dynamics induced by the change of the attractive interaction during a finite ramp time. The ramp time is varied in order to track the evolution of the dynamical phase diagram from the sudden quench to the equilibrium limit. Two different dynamical regimes are realized for quenches towards weak and strong coupling interactions. At weak coupling the dynamics depends only on the energy injected into the system, whereas a dynamics retaining memory of the initial state takes place at strong coupling. We show that this is related to a sharp transition between a weak and a strong coupling quench dynamical regime, which defines the boundaries beyond which a dynamics independent from the initial state is recovered. Comparing the dynamics in the superconducting and nonsuperconducting phases, we argue that this is due to the lack of an adiabatic connection to the equilibrium ground state for nonequilibrium superconducting states in the strong coupling quench regime.

arXiv:1708.01096•Physical review. B./Physical review. B•Nov 7, 2017•24 citations
We present a framework of semiclassical superconductivity (SC) dynamics that properly includes effects of spatial fluctuations for the attractive Hubbard model. We consider both coherent and adiabatic limits. To model the coherent SC dynamics, we develop a real-space von Neumann equation based on the time-dependent Hartree-Fock-Bogoliubov theory. Applying our method to interaction quenches in the negative-UUU Hubbard model, we show that the relaxation of SC order at weak coupling is dominated by Landau-damping. At strong coupling, we find a two-stage relaxation of the pairing field: a collapse of the synchronized oscillation of Cooper pairs due to spatial inhomogeneity, followed by a slow relaxation to a quasistationary state. SC dynamics in the adiabatic limit is described by a quantum Landau-Lifshitz equation with Ginzburg-Landau relaxation. Numerical simulations of the pump-probe process show that long-time recovery of the pairing field is dominated by defect dynamics. Our results demonstrate the important role of spatial fluctuations in both limits.
We present a Gutzwiller von Neumann dynamics (GvND) method for simulating equilibrium and nonequilibrium phenomena in strongly correlated electron systems. Our approach is a real-space formulation of the time-dependent Gutzwiller approximation method. Applying the GvND method to simulate interaction quenches in the Peierls-Hubbard model, we demonstrate the amplification of initial inhomogeneities, which in turn results in the collapse of quench-induced synchronized oscillation. Moreover, we find a dynamical phase transition separating two quasi-stationary regimes with rather distinct spatial distributions of physical quantities after the collapsed oscillation. In particular, in the strong-coupling regime, the system exhibits a dynamic phase separation in the quasi-stationary state. Our results thus underscore the importance of spatial fluctuations in the nonequilibrium dynamics of strongly correlated systems.
We investigate the quantum dynamics of a conventional superconductor, described by the attrac- tive Hubbard model, driven far from equilibrium by quenching the Hubbard attraction. We evolve the full density matrix of the system in time using the Heisenberg equation of motion, tracking the spatiotemporal dynamics of the pairing and density fields. We use very weakly inhomogeneous initial conditions to probe the stability of the homogeneous dynamics to spatial perturbations. At initial times after the quench the pairing and density fields remain homogeneous in space while fluctuating in time. However, depending on the quench strength the spatially uniform dynamics can become unstable beyond a threshold time, with emergence of spatial textures. The onset time for the instability depends on the strength of the quench and the initial inhomogeneity. At long times the system evolves towards a state with significantly reduced spatial heterogeneity. The time dependence of the global superconducting order parameter suggests two broad categories of re- sponse to a quench. One of these involves strongly inhomogeneous pairing and bond currents in the nonequilibrium superconductor.
We investigate the dynamics of superconducting fluctuations in the attractive three-dimensional Hubbard model after a quench from the disordered phase to the ordered regime. While the long-time evolution is well understood in terms of dissipative time-dependent Ginzburg-Landau models with unstable potentials, early times are more demanding due to the inseparable dynamics of the pairing fluctuations and the electronic quasiparticles. Our simulation using the time-dependent fluctuation exchange approximation treats both degrees of freedom on the same footing and reveals a nonthermal electronic regime causing a nonmonotonous growth of the fluctuations. This feature is not directly captured by the Ginzburg-Landau theory but nevertheless remains observable beyond the thermalization time of the electrons. We further explore how the growth of the order parameter fluctuations leads to an opening of a pseudogap in the electronic spectrum and identify Andreev reflections as the dominant mechanism behind the gap opening.
Pump-probe experiments have turned out as a powerful tool in order to study the dynamics of competing orders in a large variety of materials. The corresponding analysis of the data often relies on standard linear-response theory generalized to nonequilibrium situations. Here we examine the validity of such an approach for the charge and pairing response of systems with charge-density wave and (or) superconducting (SC) order. Our investigations are based on the attractive Hubbard model which we study within the time-dependent Hartree-Fock approximation. In particular, we calculate the quench and pump-probe dynamics for SC and charge order parameters in order to analyze the frequency spectra and the coupling of the probe field to the specific excitations. Our calculations reveal that the ``linear-response assumption'' is justified for small to moderate nonequilibrium situations (i.e., pump pulses) in the case of a purely charge-ordered ground state. However, the pump-probe dynamics on top of a superconducting ground state is determined by phase and amplitude modes which get coupled far from the equilibrium state indicating the failure of the linear-response assumption.
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.
We study the crossover between the sudden quench limit and the adiabatic dynamics of superconducting states in the attractive Hubbard model. We focus on the dynamics induced by the change of the attractive interaction during a finite ramp time. The ramp time is varied in order to track the evolution of the dynamical phase diagram from the sudden quench to the equilibrium limit. Two different dynamical regimes are realized for quenches towards weak and strong coupling interactions. At weak coupling the dynamics depends only on the energy injected into the system, whereas a dynamics retaining memory of the initial state takes place at strong coupling. We show that this is related to a sharp transition between a weak and a strong coupling quench dynamical regime, which defines the boundaries beyond which a dynamics independent from the initial state is recovered. Comparing the dynamics in the superconducting and nonsuperconducting phases, we argue that this is due to the lack of an adiabatic connection to the equilibrium ground state for nonequilibrium superconducting states in the strong coupling quench regime.