Prethermal Dynamics Research Articles

Systematic large flavor fTWA approach to interaction quenches in the Hubbard model

Abstract We study the nonequilibrium dynamics after an interaction quench in the two-dimensional Hubbard model using the recently introduced fermionic truncated Wigner approximation (fTWA). To assess the range of validity of the method in a systematic way, we consider the SU(N) Hubbard model with the fermion degeneracy N as a natural semiclassical expansion parameter. Using both a numerical and a perturbative analytical approach we show that fTWA is exact at least up to and including the prethermalization dynamics. We discuss the limitations of the method beyond this regime.

Slowdown of thermalization and the emergence of prethermal dynamics in disordered optical lattices

By utilizing the theoretical tools of optical thermodynamics, we investigate thermalization in nonlinear disordered lattices beyond the diffusion regime. Even under extreme levels of disorder, we analytically predict the expected thermal value of entropy and the associated Rayleigh-Jeans distribution, once thermalization ensues. In this context, we reveal a crossover point of disorder beyond which we observe a pronounced slowdown of thermalization in a regime characterized by a logarithmic scaling of equilibration times. By employing the same analytical and numerical tools, we investigate the scaling of thermalization times with respect to both the number of sites and optical temperatures. By exploring the physics of the underlying nonequilibrium response, we unveil a multitude of emerging thermal phenomena, including the development of prethermal Rayleigh-Jeans states in the presence of an optical heat bath and the deceleration of self-heating in Floquet photonic lattices. Published by the American Physical Society 2024

Screening the Coulomb interaction leads to a prethermal regime in two-dimensional bad conductors

The absence of thermalization in certain isolated many-body systems is of great fundamental interest. Many-body localization (MBL) is a widely studied mechanism for thermalization to fail in strongly disordered quantum systems, but it is still not understood precisely how the range of interactions affects the dynamical behavior and the existence of MBL, especially in dimensions D > 1. By investigating nonequilibrium dynamics in strongly disordered D = 2 electron systems with power-law interactions ∝ 1/rα and poor coupling to a thermal bath, here we observe MBL-like, prethermal dynamics for α = 3. In contrast, for α = 1, the system thermalizes, although the dynamics is glassy. Our results provide important insights for theory, especially since we obtained them on systems that are much closer to the thermodynamic limit than synthetic quantum systems employed in previous studies of MBL. Thus, our work is a key step towards further studies of ergodicity breaking and quantum entanglement in real materials.

Critical slowing down in sudden quench dynamics

We reveal a prethermal dynamical regime upon suddenly quenching to the vicinity of a quantum phase transition in the time evolution of one-dimensional spin chains. The prethermal regime is analytically found to be self-similar and its duration is governed by the ground-state energy gap. Based on analytical insights and numerical evidence, we show that this dynamical regime universally exists independently of the location of the probe site, the presence of weak interactions, or the initial state. The resulting prethermal dynamics leads to an out-of-equilibrium scaling function of the order parameter in the vicinity of the transition. Our theory suggests that sudden quench dynamics, besides probing quantum phase transitions, may give rise to a universal critical slowing down near the critical point.

Prethermalization in one-dimensional quantum many-body systems with confinement

Unconventional nonequilibrium phases with restricted correlation spreading and slow entanglement growth have been proposed to emerge in systems with confined excitations, calling their thermalization dynamics into question. Here, we show that in confined systems the thermalization dynamics after a quantum quench instead exhibits multiple stages with well separated time scales. As an example, we consider the confined Ising spin chain, in which domain walls in the ordered phase form bound states reminiscent of mesons. The system first relaxes towards a prethermal state, described by a Gibbs ensemble with conserved meson number. The prethermal state arises from rare events in which mesons are created in close vicinity, leading to an avalanche of scattering events. Only at much later times a true thermal equilibrium is achieved in which the meson number conservation is violated by a mechanism akin to the Schwinger effect. The discussed prethermalization dynamics is directly relevant to generic one-dimensional, many-body systems with confined excitations.

Realization of a discrete time crystal on 57 qubits of a quantum computer.

Unconventional dynamical phases that violate ergodicity have been a subject of extensive research in recent years. A periodically driven system is naively expected to lose all memory of its initial state due to thermalization, yet this can be avoided in the presence of many-body localization. A discrete time crystal represents a driven system whose local observables spontaneously break time translation symmetry and retain memory of the initial state indefinitely. Here, we report the observation of a discrete time crystal on a chain consisting of 57 superconducting qubits on a state-of-the-art quantum computer. We probe random initial states and compare the cases of vanishing and finite disorder to distinguish many-body localization from prethermal dynamics. We further report results on the dynamical phase transition between the discrete time crystal and a thermal regime, which is observed via critical fluctuations in the system’s subharmonic frequency response and a substantial speedup of spin depolarization.

Realization of a discrete time crystal on 57 qubits of a quantum computer

Novel dynamical phases that violate ergodicity have been a subject of extensive research in recent years. A periodically driven system is naively expected to lose all memory of its initial state due to thermalization, yet this can be avoided in the presence of many-body localization. A discrete time crystal represents a driven system whose local observables spontaneously break time translation symmetry and retain memory of the initial state indefinitely. Here we report the observation of a discrete time crystal on a chain consisting of 57 superconducting qubits on a state--of--the--art quantum computer. We probe random initial states and compare the cases of vanishing and finite disorder to distinguish many-body localization from pre-thermal dynamics. We further report results on the dynamical phase transition between the discrete time crystal and a thermal regime, which is observed via critical fluctuations in the system's sub-harmonic frequency response and a significant speed-up of spin depolarisation.

Prethermalization, thermalization, and Fermi's golden rule in quantum many-body systems

We study the prethermalization and thermalization dynamics of local observables in weakly perturbed nonintegrable systems, with Hamiltonians of the form $\hat{H}_0+g\hat{V}$, where $\hat{H}_0$ is nonintegrable and $g\hat{V}$ is a perturbation. We explore the dynamics of far from equilibrium initial states in the thermodynamic limit using a numerical linked cluster expansion (NLCE), and in finite systems with periodic boundaries using exact diagonalization. We argue that generic observables exhibit a two-step relaxation process, with a fast prethermal dynamics followed by a slow thermalizing one, only if the perturbation breaks a conserved quantity of $\hat{H}_0$ and if the value of the conserved quantity in the initial state is $\mathcal{O}(1)$ different from the one after thermalization. We show that the slow thermalizing dynamics is characterized by a rate $\propto g^2$, which can be accurately determined using a Fermi golden rule (FGR) equation. We also show that during such a slow dynamics, observables can be described using projected diagonal and Gibbs ensembles, and we contrast their accuracy.

Universal Prethermal Dynamics in Heisenberg Ferromagnets.

We study the far from equilibrium prethermal dynamics of magnons in Heisenberg ferromagnets. We show that such systems exhibit universal self-similar scaling in momentum and time of the quasiparticle distribution function, with the scaling exponents independent of microscopic details or initial conditions. We argue that the SU(2) symmetry of the Hamiltonian, which leads to a strong momentum-dependent magnon-magnon scattering amplitude, gives rise to qualitatively distinct prethermal dynamics from that recently observed in Bose gases. We compute the scaling exponents using the Boltzmann kinetic equation and incoherent initial conditions that can be realized with microwave pumping of magnons. We also compare our numerical results with analytic estimates of the scaling exponents and demonstrate the robustness of the scaling to variations in the initial conditions. Our predictions can be tested in quench experiments of spin systems in optical lattices and pump-probe experiments in ferromagnetic insulators such as yttrium iron garnet.

Prethermalization without Temperature

While a clean driven system generically absorbs energy until it reaches `infinite temperature', it may do so very slowly exhibiting what is known as a prethermal regime. Here, we show that the emergence of an additional approximately conserved quantity in a periodically driven (Floquet) system can give rise to an analogous long-lived regime. This can allow for non-trivial dynamics, even from initial states that are at a high or infinite temperature with respect to an effective Hamiltonian governing the prethermal dynamics. We present concrete settings with such a prethermal regime, one with a period-doubled (time-crystalline) response. We also present a direct diagnostic to distinguish this prethermal phenomenon from its infinitely long-lived many-body localised cousin. We apply these insights to a model of the recent NMR experiments by Rovny et al., [Phys. Rev. Lett. 120, 180603 (2018)] which, intriguingly, detected signatures of a Floquet time crystal in a clean three-dimensional material. We show that a mild but subtle variation of their driving protocol can increase the lifetime of the time-crystalline signal by orders of magnitude.

Suppression and control of prethermalization in multicomponent Fermi gases following a quantum quench

We investigate the mechanisms of control and suppression of pre-thermalization in $N$-component alkaline earth gases. To this end, we compute the short-time dynamics of the instantaneous momentum distribution and the relative population for different initial conditions after an interaction quench, accounting for the 11 peffect of initial interactions. We find that switching on an interaction that breaks the SU$(N)$ symmetry of the initial Hamiltonian, thus allowing for the occurrence of spin-changing collisions, does not necessarily lead to a suppression of pre-thermalization. However, the suppression will be most effective in the presence of SU$(N)$-breaking interactions provided the number of components $N \ge 4$ and the initial state contains a population imbalance that breaks the SU$(N)$ symmetry. We also find the conditions on the imbalance initial state that allow for a pre-thermal state to be stabilized for a certain time. Our study highlights the important role played by the initial state in the pre-thermalization dynamics of multicomponent Fermi gases. It also demonstrates that alkaline-earth Fermi gases provide an interesting playground for the study and control of pre-thermalization.

Universal Prethermal Dynamics in Gross-Neveu-Yukawa Criticality.

We study the prethermal dynamics of the Gross-Neveu-Yukawa quantum field theory, suddenly quenched in the vicinity of a critical point. We find that the universal prethermal dynamics is controlled by two fixed points depending on the size of the quench. Besides the usual equilibrium chiral Ising fixed point for a shallow quench, a dynamical chiral Ising fixed point is identified for a deep quench. Intriguingly, the latter is a nonthermal fixed point without any equilibrium counterpart due to the participation of gapless fermionic fields. We also find that in the scaling regime controlled by the equilibrium fixed point, the initial slip exponent is rendered negative if there are enough flavors of fermions, thus providing a unique signature of fermionic prethermal dynamics. We then explore the temporal crossover between the universal scaling regimes governed by the two universality classes. Possible experimental realizations are also discussed.

Universal prethermal dynamics of Bose gases quenched to unitarity.

Understanding strongly correlated phases of matter, such as the quark-gluon plasma and neutron stars, and in particular the dynamics of such systems, for example, following a Hamiltonian quench(a sudden change in some Hamiltonian parameter, such as the strength of interparticle interactions) is a fundamental challenge in modern physics. Ultracold atomic gases are excellent quantum simulators for these problems, owing to their tunable interparticle interactions and experimentally resolvable intrinsic timescales. In particular, they provide access to the unitary regime, in which the interactions are as strong as allowed by quantum mechanics. This regime has been extensively studied in Fermi gases1,2. The less-exploredunitary Bose gases3-11 offer possibilities12 such as universal physics controlled solely by the gas density13,14 and new forms of superfluidity15-17. Here, through momentum- and time-resolved studies, we explore degenerate and thermal homogeneous Bose gases quenched to unitarity. In degenerate samples, we observe universal post-quench dynamics in agreement with the emergence of a prethermal state18-24 with a universal non-zero condensed fraction22,24. In thermal gases, the dynamic and thermodynamic properties generally depend on the gas density and the temperature, but we find that they can still be expressed in terms of universal dimensionless functions. Surprisingly, we find that the total quench-induced correlation energy is independent of the gas temperature. These measurements provide quantitative benchmarks and challenges for the theory of unitary Bose gases.

Dynamical Crossovers in Prethermal Critical States.

We study the prethermal dynamics of an interacting quantum field theory with an N-component order parameter and O(N) symmetry, suddenly quenched in the vicinity of a dynamical critical point. Depending on the initial conditions, the evolution of the order parameter, and of the response and correlation functions, can exhibit a temporal crossover between universal dynamical scaling regimes governed, respectively, by a quantum and a classical prethermal fixed point, as well as a crossover from a Gaussian to a non-Gaussian prethermal dynamical scaling. Together with a recent experiment, this suggests that quenches may be used in order to explore the rich variety of dynamical critical points occurring in the nonequilibrium dynamics of a quantum many-body system. We illustrate this fact by using a combination of renormalization group techniques and a nonperturbative large-N limit.

Thermalization in the initial stage of heavy ion collisions

The high density non-abelian matter produced in heavy ion collisions is extremely anisotropic. Prethermal dynamics for the anisotropic and weakly coupled matter is discussed. Thermalization is realized with the effective kinetic theory in the leading order accuracy of the weakly coupled expansion. With the initial condition from color glass condensate, hydrodynamization time for the LHC energies is realized to be about 1 fm/c, while the thermalization happens much later than the hydrodynamization.

Initialization of hydrodynamics in relativistic heavy ion collisions with an energy-momentum transport model

A key ingredient of hydrodynamical modeling of relativistic heavy ion collisions is thermal initial conditions, an input that is the consequence of a pre-thermal dynamics which is not completely understood yet. In the paper we employ a recently developed energy-momentum transport model of the pre-thermal stage to study influence of the alternative initial states in nucleus-nucleus collisions on flow and energy density distributions of the matter at the starting time of hydrodynamics. In particular, the dependence of the results on isotropic and anisotropic initial states is analyzed. It is found that at the thermalization time the transverse flow is larger and the maximal energy density is higher for the longitudinally squeezed initial momentum distributions. The results are also sensitive to the relaxation time parameter, equation of state at the thermalization time, and transverse profile of initial energy density distribution: Gaussian approximation, Glauber Monte Carlo profiles, etc. Also, test results ensure that the numerical code based on the energy-momentum transport model is capable of providing both averaged and fluctuating initial conditions for the hydrodynamic simulations of relativistic nuclear collisions.

Quantum quench and prethermalization dynamics in a two-dimensional fermi gas with long-range interactions.

We study the effect of suddenly turning on a long-range interaction in a spinless Fermi gas in two dimensions. The short- to intermediate-time dynamics is described using the method of bosonization of the Fermi surface. The space-time dependence of the nonequilibrium fermion density matrix as well as the evolution after the quench of the discontinuity at the Fermi momentum of the momentum distribution are computed. We find that the asymptotic state predicted by bosonization is consistent with the existence of a prethermalization plateau, which is also predicted by a perturbative approach in terms of the fermionic degrees of freedom. The bosonized representation, however, explicitly allows for the construction of the generalized Gibbs ensemble describing the prethermalized state.

Pre-Thermalization Dynamics: Initial Conditions for QGP at the LHC and RHIC from Perturbative QCD

I discuss how the initial conditions for QGP-production in ultrarelativistic heavy ion collisions at the LHC and RHIC can be computed from perturbative QCD.