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

Emergent Spatial Textures from Interaction Quenches in the Hubbard Model

Sankha Subhra Bakshi, Gia-Wei Chern

Published Jan 29, 2026• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.48550/arxiv.2601.22053

Related papers

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

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

Quantum quenches in the Hubbard model: Time-dependent mean-field theory and the role of quantum fluctuations

Marco Schirò, Michele Fabrizio

We study the nonequilibrium dynamics in the fermionic Hubbard model after a sudden change of the interaction strength. To this scope, we introduce a time-dependent variational approach in the spirit of the Gutzwiller ansatz. At the saddle-point approximation, we find at half filling a sharp transition between two different regimes of small and large coherent oscillations, separated by a critical line of quenches where the system is found to relax. Any finite doping washes out the transition, leaving aside just a sharp crossover. In order to investigate the role of quantum fluctuations, we map the model onto an auxiliary quantum Ising model in a transverse field coupled to free fermionic quasiparticles. Remarkably, the Gutzwiller approximation turns out to correspond to the mean-field decoupling of this model in the limit of infinite coordination lattices. The advantage is that we can go beyond mean field and include Gaussian fluctuations around the non-equilibrium mean-field dynamics. Unlike at equilibrium, we find that quantum fluctuations become massless and eventually unstable before the mean-field dynamical critical line, which suggests they could even alter qualitatively the mean-field scenario.

arXiv:1102.1658•Physical Review B•Apr 12, 2011•118 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

Time-Dependent Mean Field Theory for Quench Dynamics in Correlated Electron Systems

Marco Schirò, Michele Fabrizio

A simple and very flexible variational approach to the out-of-equilibrium quantum dynamics in strongly correlated electron systems is introduced through a time-dependent Gutzwiller wave function. As an application, we study the simple case of a sudden change of the interaction in the fermionic Hubbard model and find at the mean-field level an extremely rich behavior. In particular, a dynamical transition between small and large quantum quench regimes is found to occur at half-filling, in accordance with the analysis of Eckstein, Phys. Rev. Lett. 103, 056403 (2009)10.1103/PhysRevLett.103.056403, obtained by dynamical mean-field theory, that turns into a crossover at any finite doping.

arXiv:1005.0992•Physical Review Letters•Aug 10, 2010•208 citations

Extended dynamic Mott transition in the two-band Hubbard model out of equilibrium

Malte Behrmann, Michele Fabrizio, Frank Lechermann

We reformulate the time-dependent Gutzwiller approximation by M. Schir\'o and M. Fabrizio, [Phys. Rev. Lett. 105, 076401 (2010)] in the framework of slave-boson mean-field theory, which is used to investigate the dynamical Mott transition of the generic two-band Hubbard model at half filling upon an interaction quench. Interorbital fluctuations lead to notable changes with respect to the single-band case. The singular dynamical transition is replaced by a broad regime of long-lived fluctuations between metallic and insulating states, accompanied by intriguing precursor behavior. A mapping to a spin model proves helpful to analyze the different regions in terms of the evolution of an Ising-like order parameter. Contrary to the static case, singlet occupations remain vital in the Mott-insulating regime with finite Hund's exchange.

arXiv:1304.6013•Physical Review B•Jul 15, 2013•11 citations
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 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 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 study the nonequilibrium dynamics in the fermionic Hubbard model after a sudden change of the interaction strength. To this scope, we introduce a time-dependent variational approach in the spirit of the Gutzwiller ansatz. At the saddle-point approximation, we find at half filling a sharp transition between two different regimes of small and large coherent oscillations, separated by a critical line of quenches where the system is found to relax. Any finite doping washes out the transition, leaving aside just a sharp crossover. In order to investigate the role of quantum fluctuations, we map the model onto an auxiliary quantum Ising model in a transverse field coupled to free fermionic quasiparticles. Remarkably, the Gutzwiller approximation turns out to correspond to the mean-field decoupling of this model in the limit of infinite coordination lattices. The advantage is that we can go beyond mean field and include Gaussian fluctuations around the non-equilibrium mean-field dynamics. Unlike at equilibrium, we find that quantum fluctuations become massless and eventually unstable before the mean-field dynamical critical line, which suggests they could even alter qualitatively the mean-field scenario.
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.
A simple and very flexible variational approach to the out-of-equilibrium quantum dynamics in strongly correlated electron systems is introduced through a time-dependent Gutzwiller wave function. As an application, we study the simple case of a sudden change of the interaction in the fermionic Hubbard model and find at the mean-field level an extremely rich behavior. In particular, a dynamical transition between small and large quantum quench regimes is found to occur at half-filling, in accordance with the analysis of Eckstein, Phys. Rev. Lett. 103, 056403 (2009)10.1103/PhysRevLett.103.056403, obtained by dynamical mean-field theory, that turns into a crossover at any finite doping.
We reformulate the time-dependent Gutzwiller approximation by M. Schir\'o and M. Fabrizio, [Phys. Rev. Lett. 105, 076401 (2010)] in the framework of slave-boson mean-field theory, which is used to investigate the dynamical Mott transition of the generic two-band Hubbard model at half filling upon an interaction quench. Interorbital fluctuations lead to notable changes with respect to the single-band case. The singular dynamical transition is replaced by a broad regime of long-lived fluctuations between metallic and insulating states, accompanied by intriguing precursor behavior. A mapping to a spin model proves helpful to analyze the different regions in terms of the evolution of an Ising-like order parameter. Contrary to the static case, singlet occupations remain vital in the Mott-insulating regime with finite Hund's exchange.