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

Fate of dynamical phases of a BCS superconductor beyond the dissipationless regime

H. P. Ojeda Collado, Gonzalo Usaj, J. Lorenzana, C. A. Balseiro

Published May 15, 2019• 21 citations
arXivHTMLPDF

Abstract

The BCS model of an isolated superconductor initially prepared in a nonequilibrium state predicts the existence of interesting dynamical phenomena in the time-dependent order parameter as decaying oscillations, persistent oscillations, and overdamped dynamics. To make contact with real systems remains an open challenge as one needs to introduce dissipation due to the environment in a self-consistent computation. Here, we reach this goal with the use of the Keldysh formalism to treat the effect of a thermal bath. We show that, contrary to the dissipationless case, all dynamical phases reach the equilibrium order parameter in a characteristic time that depends on the coupling with the bath. Remarkably, as time evolves, the overdamped phase shows a fast crossover where the superconducting order parameter recovers to reach a state with a well-developed long range order that tends towards equilibrium with the damped Higgs mode oscillations. Our results provide a benchmark for the description of the dynamics of real out-of-equilibrium superconductors relevant for quantum technological applications.

DOI
10.1103/physrevb.99.174509

Related papers

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

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

Nonequilibrium relaxation in neutral BCS superconductors: Ginzburg-Landau approach with Landau damping in real time

S. M. Alamoudi, D. Boyanovsky, Shang-Yung Wang

We present a field-theoretical method to obtain consistently the equations of motion for small amplitude fluctuations of the order parameter directly in real time for a homogeneous, neutral BCS superconductor. This method allows to study the nonequilibrium relaxation of the order parameter as an initial value problem. We obtain the Ward identities and the effective actions for small phase the amplitude fluctuations to one-loop order. Focusing on the long wavelength, low-frequency limit near the critical point, we obtain the time-dependent Ginzburg-Landau effective action to one-loop order, which is nonlocal as a consequence of Landau damping. The nonequilibrium relaxation of the phase and amplitude fluctuations is studied directly in real time. The long-wavelength phase fluctuation (Bogoliubov-Anderson-Goldstone mode) is overdamped by Landau damping and the relaxation time scale diverges at the critical point, revealing critical slowing down.

arXiv:cond-mat/0204239•Physical review. B, Condensed matter•Nov 5, 2002•4 citations

Novel Non–Equilibrium Dynamics in Superconductors: Induced Superconductivity and Higgs Modes

Nikolaj Bittner

This thesis reports on theoretical investigations of the ultrafast nonequilibrium dynamics in superconductors. In particular, two nonequilibrium scenarios are modeled by combining several analytical and numerical methods: transient superconducting phase transition and collective excitations of a superconducting condensate. In the first part of the thesis the possibility to induce singlet superconducting correlations far from equilibrium is explored by applying time-dependent Lanczos algorithm to the one-dimensional extended Hubbard model. Here, the system is driven into a nonequilibrium state by (i) interaction quench and (ii) quench by pulse. In both cases the buildup of the superconducting correlations is studied in view of several correlation functions. Moreover, a transient Meissner effect, which reflects fingerprints of the induced superconductivity, is observed in calculations of the important for experiments time-dependent optical conductivity. In addition, the differences between nonequilibrium dynamics of the extended Hubbard model initiated by interaction quench and by pulse quench are discussed in detail. The second part is devoted to the excitation of collective modes in a superconductor out of equilibrium. First of all, focus is given to the properties of collective modes close to equilibrium conditions, where a linearized response can be expected. Within the framework of the (Nambu-) matrix kinetic theory in the clean limit a comprehensive gauge-invariant description of the collective modes associated with phase fluctuations of the order parameter is analytically formulated for the generalized case of two- band superconducting systems. Particular emphasis is on the application to the recently discovered noncentrosymmetric superconductors. A detailed analysis of the collective modes includes specification of the gauge mode (common to all superconductors) and the Leggett mode (collective fluctuations of the interband phase difference), and investigation of their subtle interplay in view of validity of the charge conservation law and the participation in the Anderson-Higgs mechanism. Further, the nonequilibrium excitation of the superconducting gaps for the case of two-band superconductors with singlet pairing correlations is investigated in a view of optical pump-probe experiments. For this purpose a semi-numerical method based on the density- matrix theory is formulated. In this case additional amplitude (Higgs) modes were identified, which allow to directly detect the properties of the superconducting condensates as a function of time. Moreover, a new interplay between the phase (Leggett) and amplitude (Higgs) modes is studied in detail. In particular it was found that the nonadiabatic excitation process leads to strong coupling between both modes. Finally, based on calculations of the time-dependent optical conductivity the prediction is made that the coupled Higgs and Leggett modes can be clearly visible in the pump-probe absorption spectra.

•Jan 1, 2017•0 citations

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

Signatures of nonadiabatic BCS state dynamics in pump-probe conductivity

H. Krull, Dirk Manske, Götz Silvester Uhrig +1 more

We theoretically study the pump-probe response of nonequilibrium BCS superconductors coupled to optical phonons. For ultrashort pump pulses a nonadiabatic regime emerges, which is characterized by oscillations of the superconducting order parameter as well as by the generation of coherent phonons. Using the density-matrix formalism, we compute the pump-probe response in the nonadiabatic regime of the coupled Bogoliubov quasiparticle-phonon system and determine the signatures of the order parameter and the phonon oscillations in the pump-probe conductivity. We find that the nonadiabatic dynamics of the BCS superconductor reflects itself in oscillations of the pump-probe response as functions of delay time $\ensuremath{\delta}t$ between pump and probe pulses. We argue that from the analysis of this oscillatory behavior both frequency and decay time of the algebraically decaying order-parameter oscillations can be inferred. Similarly, the coherent phonons are evidenced in the pump-probe conductivity by oscillations with the frequency of the phonons. Remarkably, we find that the oscillatory response in the pump-probe conductivity is resonantly enhanced when the frequency of the order-parameter oscillations is tuned to the phonon energy.

arXiv:1309.7318•Physical Review B•Jul 25, 2014•69 citations
The BCS model of an isolated superconductor initially prepared in a nonequilibrium state predicts the existence of interesting dynamical phenomena in the time-dependent order parameter as decaying oscillations, persistent oscillations, and overdamped dynamics. To make contact with real systems remains an open challenge as one needs to introduce dissipation due to the environment in a self-consistent computation. Here, we reach this goal with the use of the Keldysh formalism to treat the effect of a thermal bath. We show that, contrary to the dissipationless case, all dynamical phases reach the equilibrium order parameter in a characteristic time that depends on the coupling with the bath. Remarkably, as time evolves, the overdamped phase shows a fast crossover where the superconducting order parameter recovers to reach a state with a well-developed long range order that tends towards equilibrium with the damped Higgs mode oscillations. Our results provide a benchmark for the description of the dynamics of real out-of-equilibrium superconductors relevant for quantum technological applications.
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.
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 present a field-theoretical method to obtain consistently the equations of motion for small amplitude fluctuations of the order parameter directly in real time for a homogeneous, neutral BCS superconductor. This method allows to study the nonequilibrium relaxation of the order parameter as an initial value problem. We obtain the Ward identities and the effective actions for small phase the amplitude fluctuations to one-loop order. Focusing on the long wavelength, low-frequency limit near the critical point, we obtain the time-dependent Ginzburg-Landau effective action to one-loop order, which is nonlocal as a consequence of Landau damping. The nonequilibrium relaxation of the phase and amplitude fluctuations is studied directly in real time. The long-wavelength phase fluctuation (Bogoliubov-Anderson-Goldstone mode) is overdamped by Landau damping and the relaxation time scale diverges at the critical point, revealing critical slowing down.
This thesis reports on theoretical investigations of the ultrafast nonequilibrium dynamics in superconductors. In particular, two nonequilibrium scenarios are modeled by combining several analytical and numerical methods: transient superconducting phase transition and collective excitations of a superconducting condensate. In the first part of the thesis the possibility to induce singlet superconducting correlations far from equilibrium is explored by applying time-dependent Lanczos algorithm to the one-dimensional extended Hubbard model. Here, the system is driven into a nonequilibrium state by (i) interaction quench and (ii) quench by pulse. In both cases the buildup of the superconducting correlations is studied in view of several correlation functions. Moreover, a transient Meissner effect, which reflects fingerprints of the induced superconductivity, is observed in calculations of the important for experiments time-dependent optical conductivity. In addition, the differences between nonequilibrium dynamics of the extended Hubbard model initiated by interaction quench and by pulse quench are discussed in detail. The second part is devoted to the excitation of collective modes in a superconductor out of equilibrium. First of all, focus is given to the properties of collective modes close to equilibrium conditions, where a linearized response can be expected. Within the framework of the (Nambu-) matrix kinetic theory in the clean limit a comprehensive gauge-invariant description of the collective modes associated with phase fluctuations of the order parameter is analytically formulated for the generalized case of two- band superconducting systems. Particular emphasis is on the application to the recently discovered noncentrosymmetric superconductors. A detailed analysis of the collective modes includes specification of the gauge mode (common to all superconductors) and the Leggett mode (collective fluctuations of the interband phase difference), and investigation of their subtle interplay in view of validity of the charge conservation law and the participation in the Anderson-Higgs mechanism. Further, the nonequilibrium excitation of the superconducting gaps for the case of two-band superconductors with singlet pairing correlations is investigated in a view of optical pump-probe experiments. For this purpose a semi-numerical method based on the density- matrix theory is formulated. In this case additional amplitude (Higgs) modes were identified, which allow to directly detect the properties of the superconducting condensates as a function of time. Moreover, a new interplay between the phase (Leggett) and amplitude (Higgs) modes is studied in detail. In particular it was found that the nonadiabatic excitation process leads to strong coupling between both modes. Finally, based on calculations of the time-dependent optical conductivity the prediction is made that the coupled Higgs and Leggett modes can be clearly visible in the pump-probe absorption spectra.
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 theoretically study the pump-probe response of nonequilibrium BCS superconductors coupled to optical phonons. For ultrashort pump pulses a nonadiabatic regime emerges, which is characterized by oscillations of the superconducting order parameter as well as by the generation of coherent phonons. Using the density-matrix formalism, we compute the pump-probe response in the nonadiabatic regime of the coupled Bogoliubov quasiparticle-phonon system and determine the signatures of the order parameter and the phonon oscillations in the pump-probe conductivity. We find that the nonadiabatic dynamics of the BCS superconductor reflects itself in oscillations of the pump-probe response as functions of delay time $\ensuremath{\delta}t$ between pump and probe pulses. We argue that from the analysis of this oscillatory behavior both frequency and decay time of the algebraically decaying order-parameter oscillations can be inferred. Similarly, the coherent phonons are evidenced in the pump-probe conductivity by oscillations with the frequency of the phonons. Remarkably, we find that the oscillatory response in the pump-probe conductivity is resonantly enhanced when the frequency of the order-parameter oscillations is tuned to the phonon energy.