Near a critical point, the equilibrium relaxation time of a system diverges and any change of control/thermodynamic parameters leads to nonequilibrium behavior. The Kibble-Zurek problem is to determine the dynamical evolution of the system parametrically close to its critical point when the change is parametrically slow. The nonequilibrium behavior in this limit is controlled entirely by the critical point and the details of the trajectory of the system in parameter space (the protocol) close to the critical point. Together, they define a universality class consisting of critical exponents, discussed in the seminal work by Kibble and Zurek, and scaling functions for physical quantities, which have not been discussed hitherto. In this article, we give an extended and pedagogical discussion of the universal content in the Kibble-Zurek problem. We formally define a scaling limit for physical quantities near classical and quantum transitions for different sets of protocols. We report computations of a few scaling functions in model Gaussian and large-
N problems and prove their universality with respect to protocol choice. We also introduce a protocol in which the critical point is approached asymptotically at late times with the system marginally out of equilibrium, wherein logarithmic violations to scaling and anomalous dimensions occur even in the simple Gaussian problem.
We study the out-of-equilibrium Kibble-Zurek (KZ) dynamics in quantum Ising chains in a transverse field, driven by a time-dependent longitudinal field
h(t)=t/ts (
ts is the time scale of the protocol), across their first-order quantum transitions (FOQTs) at
h=0. The KZ protocol starts at time
ti<0 from the negatively magnetized ground state for
hi=ti/ts<0. Then, the system evolves unitarily up to a time
tf>0, such that the magnetization of the state at time
tf is positive. In finite-size systems, the KZ dynamics develops out-of-equilibrium finite-size scaling (OFSS) behaviors. Their scaling variables depend either exponentially or with a power law on the size, depending on the boundary conditions (BC). The OFSS functions can be computed in effective models restricted to appropriate low-energy (magnetized and/or kink) states. The KZ scaling behavior drastically changes in the thermodynamic limit (TL), defined as the infinite-size limit keeping
t and
ts fixed, which appears substantially unrelated with the OFSS regime, because it involves higher-energy multi-kink states, which are irrelevant in the OFSS limit. The numerical analyses of the KZ dynamics in the TL show the emergence of a quantum spinodal-like scaling behavior at the FOQTs for all considered BC, which is independent of the BC. The longitudinal magnetization changes sign at
h(t)=h∗>0, where
h∗ decreases with increasing
ts, as
h∗∼1/lnts. Moreover, in the large-
ts limit, the time-dependence of the magnetization is described by a universal function of
Ω=t/τs, with
τs=ts/lnts.
A system gradually driven through a symmetry-breaking phase transition is subject to the Kibble-Zurek mechanism (KZM). As a consequence of the critical slowing down, its state cannot follow local equilibrium, and its evolution becomes nonadiabatic near the critical point. In the simplest approximation, that stage can be regarded as an ``impulse'' where the state of the system remains unchanged. It leads to the correct KZM scaling laws. However, such a ``freeze-out'' might suggest that the coherence length of the nascent order parameter remains unchanged as the critical region is traversed. By contrast, the original causality-based discussion emphasized the role of the sonic horizon: Domains of the broken symmetry phase can expand with a velocity limited by the speed of the relevant sound. This effect was demonstrated in the quantum Ising chain where the dynamical exponent
z=1 and quasiparticles excited by the transition have a fixed speed of sound. To elucidate the role of the sonic horizon, in this paper we study two systems with
z>1 where the speed of sound is no longer fixed, and the fastest excited quasiparticles set the size of the sonic horizon. Their effective speed decays with increasing transition time. In the extreme case, the dynamical exponent
z can diverge such as in the Griffiths region of the random Ising chain where the localization of excited quasiparticles freezes the growth of the correlation range when the critical region is traversed. Of particular interest is an example with
z<1, the long-range extended Ising chain, where there is no upper limit to the velocity of excited quasiparticles with small momenta. Initially, the power-law tail of the correlation function grows adiabatically, but in the nonadiabatic stage it lags behind the adiabatic evolution.
We study the effects of symmetry-breaking perturbations in the out-of-equilibrium quantum dynamics of many-body systems slowly driven by a time-dependent symmetry-preserving parameter, across the quantum critical regime associated with a continuous quantum transition (CQT). For this purpose, we analyze the out-of-equilibrium dynamics arising from generalized Kibble-Zurek (KZ) protocols, within a dynamic renormalization-group framework allowing for finite-size systems. We show that the time dependence of generic observables develops an out-of-equilibrium finite-size scaling (FSS) behavior, arising from the interplay between the timescale
ts of the parameter variations in the KZ protocol, the size
L of the system, and the strength
h of the symmetry-breaking perturbation, in the limit of large
ts and
L. Moreover, scaling arguments based on the first-order adiabatic approximation of slow variations in quantum systems allow us to characterize the approach to the adiabatic regimes for some limits of the model parameters (for example, when we take ${t}_{s}\ensuremath{\rightarrow}\ensuremath{\infty}$ before $L\ensuremath{\rightarrow}\ensuremath{\infty}$), predicting asymptotic power-law suppressions of the nonadiabatic behaviors in the adiabatic limits. This out-of-equilibrium FSS is supported by numerical analyses for the paradigmatic quantum Ising chain along generalized KZ protocols, with a time-dependent transverse field crossing its CQT, in the presence of a static longitudinal field breaking the
Z2 symmetry.
Kibble-Zurek mechanism is widely known to appear in the transverse-field quantum Ising chain in the thermodynamic limit at zero temperature, having notorious characteristics, like the divergence of its relaxation time. In this work, I present the same effect in a simple system, the driven underdamped Brownian motion. Using linear response theory, I show the appropriate limits where the Kibble-Zurek mechanism happens. The divergence of the relaxation time, the high-temperature condition in the initial thermal equilibrium state, a new Kibble-Zurek scaling at sudden processes, and the pausing effect in the optimal protocol are presented as consequences.
The transverse field in the quantum Ising chain is linearly ramped from the para- to the ferromagnetic phase across the quantum critical point at a rate characterized by a quench time ${\ensuremath{\tau}}_{Q}$. We calculate a connected kink-kink correlator in the final state at zero transverse field. The correlator is a sum of two terms: a negative (antibunching) Gaussian that depends on the Kibble-Zurek correlation length only and a positive term that depends on a second longer scale of length. The second length is made longer by dephasing of the state excited near the critical point during the following ramp across the ferromagnetic phase. This interpretation is corroborated by considering a linear ramp that is halted in the ferromagnetic phase for a finite waiting time and then continued at the same rate as before the halt. The extra time available for dephasing increases the second scale of length that asymptotically grows linearly with the waiting time. The dephasing also suppresses magnitude of the second term, making it negligible for waiting times much longer than ${\ensuremath{\tau}}_{Q}$. The same dephasing can be obtained with a smooth ramp that slows down in the ferromagnetic phase. Assuming sufficient dephasing, we also obtain higher order kink correlators and the ferromagnetic correlation function.