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:2609.13426

Spatial instability of a homogeneous superconductor after an interaction quench

Medha Roy Choudhury, Sankha Subhra Bakshi, Pinaki Majumdar, Amit Ghosal

Published Sep 11, 2026• 0 citations
arXivHTMLPDF

Abstract

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.

Related papers

Emergent Spatial Textures from Interaction Quenches in the Hubbard Model

Sankha Subhra Bakshi, Gia-Wei Chern

Interaction quenches in strongly correlated electron systems provide a powerful route to probe nonequilibrium many-body dynamics. For the Hubbard model, nonequilibrium dynamical mean-field theory has revealed coherent post-quench oscillations, dynamical crossovers, and long-lived transient regimes. However, these studies are largely restricted to spatially homogeneous dynamics and therefore neglect the role of spatial structure formation during ultrafast evolution. Here we investigate interaction quenches in the half-filled Hubbard model using a real-space time-dependent Gutzwiller framework. We show that homogeneous nonequilibrium dynamics is generically unstable: even arbitrarily weak spatial fluctuations grow dynamically and drive the system toward intrinsically inhomogeneous states. Depending on the interaction strength, the post-quench evolution exhibits spatial differentiation, nucleation, and slow coarsening of Mott-like domains. Our results establish spatial self-organization as a generic feature of far-from-equilibrium correlated matter and reveal a fundamental limitation of spatially homogeneous nonequilibrium theories.

arXiv:2601.22053•Jan 29, 2026•0 citations

Nonequilibrium dynamics of superconductivity in the attractive Hubbard model

Gia-Wei Chern, Kipton M. Barros

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.

arXiv:1803.04118•Physical review. B./Physical review. B•Jan 30, 2019•11 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

Negative superfluid density and spatial instabilities in driven superconductors

Andrey Grankin, Victor M. Galitski, Vadim Oganesyan

We consider excitation of Higgs modes via the modulation of the BCS coupling within the Migdal-Eliashberg-Keldysh theory of time-dependent superconductivity. Despite the presence of phonons, which break integrability, we observe Higgs amplitude oscillations reminiscent of the integrable case. The dynamics of quasiparticles follows from the effective Bogolyubov-de Gennes equations, which represent a Floquet problem for the Bogoliubov quasiparticles. We find that when the Floquet-Bogoliubov bands overlap, the homogeneous solution formally leads to a negative superfluid density, which is no longer proportional to the amplitude of the order parameter. This result indicates an instability, which we explore using spatially-resolved BdG equations. Spontaneous appearance of spatial inhomogeneities in the order parameter is observed and they first occur when the superfluid density becomes unphysical. We conclude that the homogeneous solution to time-dependent superconductivity is generally unstable and breaks up into a complicated spatial landscape via an avalanche of topological excitations.

arXiv:2501.08216•Physical review. B./Physical review. B•Feb 25, 2026•0 citations

Comment on "Nonequilibrium dynamics of superconductivity in the attractive Hubbard model"

Maxim Dzero, Emil A. Yuzbashyan, B. L. Altshuler

In a recent preprint [arXiv:1803.04118v2] Chern and Barros report numerical simulations of the mean-field interaction quench dynamics, $U_i\to U_f$, of the attractive Hubbard model that confirm our earlier prediction [Europhys. Lett. 85, 20004 (2008), arXiv:0805.2798] of spontaneous eruption of spatial inhomogeneities in the post-quench state with periodically oscillating superconducting order. Chern and Barros attribute this instability with respect to spatial fluctuations to the large magnitude of the final Hubbard coupling $U_f$. We point out that this interpretation is inaccurate and discuss further work necessary to numerically verify the mechanism of the instability and the nature of the steady state.

arXiv:1806.03474•Jun 9, 2018•0 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
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.
Interaction quenches in strongly correlated electron systems provide a powerful route to probe nonequilibrium many-body dynamics. For the Hubbard model, nonequilibrium dynamical mean-field theory has revealed coherent post-quench oscillations, dynamical crossovers, and long-lived transient regimes. However, these studies are largely restricted to spatially homogeneous dynamics and therefore neglect the role of spatial structure formation during ultrafast evolution. Here we investigate interaction quenches in the half-filled Hubbard model using a real-space time-dependent Gutzwiller framework. We show that homogeneous nonequilibrium dynamics is generically unstable: even arbitrarily weak spatial fluctuations grow dynamically and drive the system toward intrinsically inhomogeneous states. Depending on the interaction strength, the post-quench evolution exhibits spatial differentiation, nucleation, and slow coarsening of Mott-like domains. Our results establish spatial self-organization as a generic feature of far-from-equilibrium correlated matter and reveal a fundamental limitation of spatially homogeneous nonequilibrium theories.
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 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.
We consider excitation of Higgs modes via the modulation of the BCS coupling within the Migdal-Eliashberg-Keldysh theory of time-dependent superconductivity. Despite the presence of phonons, which break integrability, we observe Higgs amplitude oscillations reminiscent of the integrable case. The dynamics of quasiparticles follows from the effective Bogolyubov-de Gennes equations, which represent a Floquet problem for the Bogoliubov quasiparticles. We find that when the Floquet-Bogoliubov bands overlap, the homogeneous solution formally leads to a negative superfluid density, which is no longer proportional to the amplitude of the order parameter. This result indicates an instability, which we explore using spatially-resolved BdG equations. Spontaneous appearance of spatial inhomogeneities in the order parameter is observed and they first occur when the superfluid density becomes unphysical. We conclude that the homogeneous solution to time-dependent superconductivity is generally unstable and breaks up into a complicated spatial landscape via an avalanche of topological excitations.
In a recent preprint [arXiv:1803.04118v2] Chern and Barros report numerical simulations of the mean-field interaction quench dynamics, Ui→UfU_i\to U_fUi​→Uf​, of the attractive Hubbard model that confirm our earlier prediction [Europhys. Lett. 85, 20004 (2008), arXiv:0805.2798] of spontaneous eruption of spatial inhomogeneities in the post-quench state with periodically oscillating superconducting order. Chern and Barros attribute this instability with respect to spatial fluctuations to the large magnitude of the final Hubbard coupling UfU_fUf​. We point out that this interpretation is inaccurate and discuss further work necessary to numerically verify the mechanism of the instability and the nature of the steady state.
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.