Kibble–Zurek Scaling in Phase Transitions

• Kibble–Zurek scaling is a framework that unifies non-equilibrium dynamics and defect formation in systems undergoing continuous phase transitions.
• It employs critical exponents and quench rates to predict power-law scaling of topological defect density and correlation lengths across diverse regimes.
• Experimental studies in ultracold gases, ion crystals, and superconductors validate its universal application and precise control of defects in phase transitions.

The Kibble–Zurek scaling describes universal non-equilibrium dynamics and defect formation in systems driven through continuous (second-order) phase transitions at finite rates. It arises from the interplay of critical slowing down near the transition and the finite speed of external driving, leading to characteristic power-law scaling laws for the density of topological defects, correlation lengths, and dynamical observables as functions of the quench rate. Kibble–Zurek scaling unifies phenomena across cosmology, condensed matter physics, ultracold gases, quantum criticality, and beyond, and is determined solely by equilibrium critical exponents and universality classes.

1. Fundamental Mechanism and Scaling Laws

In a continuous phase transition, both the equilibrium correlation length $\xi$ and relaxation time $\tau$ diverge algebraically as the system approaches the critical point: $$\xi(\epsilon) \sim \xi_0 |\epsilon|^{-\nu}, \qquad \tau(\epsilon) \sim \tau_0 |\epsilon|^{-\nu z}$$ where $\epsilon$ measures the distance from criticality, $\nu$ and $z$ are, respectively, the static and dynamical critical exponents. For a linear quench $\epsilon(t)=t/\tau_Q$ with quench time $\tau_Q$, the system cannot maintain equilibrium near criticality due to critical slowing down. The "freeze-out" time $\hat t$ and length $\hat\xi$ are determined by equating the instantaneous relaxation time to the time remaining to the transition: $\tau$0 After traversing the transition, domains of linear size $\tau$1 form, producing topological defects whose density in $\tau$2 spatial dimensions scales as: $\tau$3 These laws apply to both classical and quantum continuous transitions, provided the critical exponents are replaced by the appropriate universality class values (Lee et al., 2023, Natsuume et al., 2017, Silvi et al., 2015).

2. Universality, Experimental Verification, and Protocol Dependence

Kibble–Zurek scaling is universal: the scaling exponents and functions depend only on intrinsic critical exponents $\tau$4, spatial dimension $\tau$5, and protocol exponent $\tau$6 (for $\tau$7). Scaling functions for one- and two-point observables,

$\tau$8

are universal within the scaling limit and collapse for different protocols sharing the same $\tau$9 (Chandran et al., 2012). The theoretical framework is validated in quantum field theory (Ising, $$\xi(\epsilon) \sim \xi_0 |\epsilon|^{-\nu}, \qquad \tau(\epsilon) \sim \tau_0 |\epsilon|^{-\nu z}$$0), cold atomic gases, Josephson junctions, ion crystals, Rydberg atom arrays, and holographic superconductors (Lee et al., 2023, Hódsági et al., 2020, Su et al., 2013, Ulm et al., 2013, Zhang et al., 12 May 2025, Bu et al., 2019, Natsuume et al., 2017).

Experimental implementation involves controlling temperature, interaction, or external fields to enforce the prescribed quench, and direct detection of vortex lines, kinks, domain walls, or excitations. For the superfluid transition in a homogeneous Fermi gas, both temperature and interaction quenches yield the same Kibble–Zurek exponent $$\xi(\epsilon) \sim \xi_0 |\epsilon|^{-\nu}, \qquad \tau(\epsilon) \sim \tau_0 |\epsilon|^{-\nu z}$$1 consistent with 3D U(1) universality ($$\xi(\epsilon) \sim \xi_0 |\epsilon|^{-\nu}, \qquad \tau(\epsilon) \sim \tau_0 |\epsilon|^{-\nu z}$$2, $$\xi(\epsilon) \sim \xi_0 |\epsilon|^{-\nu}, \qquad \tau(\epsilon) \sim \tau_0 |\epsilon|^{-\nu z}$$3) (Lee et al., 2023). In ion crystals undergoing a symmetry-breaking zigzag transition, inhomogeneity and finite size alter the scaling exponent to $$\xi(\epsilon) \sim \xi_0 |\epsilon|^{-\nu}, \qquad \tau(\epsilon) \sim \tau_0 |\epsilon|^{-\nu z}$$4 as predicted for the inhomogeneous Kibble–Zurek mechanism (Ulm et al., 2013).

3. Extensions: Finite Size, Symmetry Breaking, and Exceptional Criticality

Kibble–Zurek scaling persists in finite size and near-critical crossovers, provided relevant variables are rescaled in accord with the ramp speed. In Rydberg atom arrays, precise scaling is observed if the system size $$\xi(\epsilon) \sim \xi_0 |\epsilon|^{-\nu}, \qquad \tau(\epsilon) \sim \tau_0 |\epsilon|^{-\nu z}$$5 and symmetry-breaking field $$\xi(\epsilon) \sim \xi_0 |\epsilon|^{-\nu}, \qquad \tau(\epsilon) \sim \tau_0 |\epsilon|^{-\nu z}$$6 obey $$\xi(\epsilon) \sim \xi_0 |\epsilon|^{-\nu}, \qquad \tau(\epsilon) \sim \tau_0 |\epsilon|^{-\nu z}$$7 and $$\xi(\epsilon) \sim \xi_0 |\epsilon|^{-\nu}, \qquad \tau(\epsilon) \sim \tau_0 |\epsilon|^{-\nu z}$$8 as the ramp speed $$\xi(\epsilon) \sim \xi_0 |\epsilon|^{-\nu}, \qquad \tau(\epsilon) \sim \tau_0 |\epsilon|^{-\nu z}$$9 is varied. This joint finite-size scaling ensures collapse of correlation functions and restoration of exponents even when the critical point is smeared (near-critical Kibble–Zurek scaling) (Zhang et al., 12 May 2025).

Tricritical points present multiple relevant operators; "tangential" Kibble–Zurek ramps can selectively probe subleading scaling exponents (e.g., $\epsilon$0) by following the subleading direction in parameter space, as realized in Rydberg ladders of Ising and Potts symmetry (Wang et al., 19 May 2025). At the Yang–Lee edge singularity, Kibble–Zurek scaling holds with exponents determined by the underlying finite-size (0+1)D or (1+1)D non-Hermitian universality class, despite the absence of topological-defect formation (Yin et al., 2016).

4. Breakdown, Crossover, and Limits of Universal Scaling

Kibble–Zurek scaling precisely describes defect production only in the slow-quench (adiabatic–impulse–adiabatic) regime. For finite-depth or fast quenches (i.e., quench rate above a critical threshold $\epsilon$1 set by the finite quench range $\epsilon$2), the freeze-out point may lie outside the critical region, and the standard scaling laws cease to apply. In this regime, the density of topological defects saturates at a value determined by the control range rather than the quench rate: $\epsilon$3 with a crossover at $\epsilon$4 (Rao et al., 7 Jun 2025, Zeng et al., 2022). This breakdown reflects the finite maximal relaxation time and is experimentally confirmed in trapped-ion simulations of the Landau–Zener and Rice–Mele models.

In first-order phase transitions, standard Kibble–Zurek scaling of defect density fails due to the presence of an intrinsic, co-moving symmetry-breaking field $\epsilon$5 that preselects one ordered phase, preventing universal defect formation. Nevertheless, full finite-time scaling for observables such as the order parameter and correlation length is preserved with exponents controlled by the off-equilibrium renormalization group near the spinodal (Zhong, 2 Jan 2025).

5. Quantum, Classical, and Crossover Regimes

The Kibble–Zurek mechanism operates in both classical and quantum critical regimes, but the governing exponents can switch depending on the quench rate and quantum-to-classical crossover scale $\epsilon$6. For slow quenches, scaling exponents of the quantum phase transition (e.g., $\epsilon$7 for 1D Ising) apply. For fast quenches, classical (mean-field) exponents (e.g., $\epsilon$8) are recovered, with the crossover point calculable from Ginzburg-criterion or RG arguments (Silvi et al., 2015). This duality clarifies observations in cold atom and ion dynamics where both regimes are manifest.

6. Dissipation, Open Systems, and Non-Markovian Effects

In open quantum systems, Kibble–Zurek scaling is affected by the nature of system-bath coupling. Markovian dissipation typically leads to anti-Kibble–Zurek (AKZ) behavior (faster quenches producing fewer defects), and generally degrades universality. At specific limits, such as a loss difference between sublattices in a Lindblad setting, universal "dissipative Kibble–Zurek" power-law scaling can be restored, while the defect density measured via total residual particles remains immune to the AKZ effect (Kou et al., 2024). In non-Markovian environments, such as an open quantum Rabi model coupled to an Ohmic bath, criticality can be induced by the environment itself (e.g., a BKT transition), and Kibble–Zurek scaling is preserved with power-law scaling of excitation energy evaluated at the freeze-out time (Pirozzi et al., 17 Mar 2026). When the system's thermal bath is quenched to a critical point, defect density scales as $\epsilon$9 for thermal critical points and $\nu$0 for quantum critical points, with lower defect production than in the isolated scenario (Bácsi et al., 2022).

7. Irreversible Entropy Production and Higher Cumulants

Beyond defect density, Kibble–Zurek scaling governs other non-equilibrium quantities, such as irreversible entropy production and higher-order statistics. For slow quenches,

$\nu$1

where $\nu$2 is the exponent of the equilibrium susceptibility with respect to the ramped parameter (e.g., specific heat or magnetic susceptibility exponents for temperature or field quenches) (Deffner, 2017). Higher moments and cumulants of the excess work, such as in the Ising field theory, follow similar scaling with exponents determined by the universality class and ramp protocol (Hódsági et al., 2020).

---

References:

• (Lee et al., 2023) Observation of universal Kibble-Zurek scaling in an atomic Fermi superfluid
• (Natsuume et al., 2017) Kibble-Zurek scaling in holography
• (Silvi et al., 2015) Crossover from the classical to the quantum Kibble-Zurek scaling
• (Zhong, 2 Jan 2025) Is there Kibble-Zurek scaling of topological defects in first-order phase transitions?
• (Hódsági et al., 2020) Kibble-Zurek mechanism in the Ising Field Theory
• (Rao et al., 7 Jun 2025) Verified Universal Breakdown of Kibble-Zurek Scaling in Fast Quenches
• (Zeng et al., 2022) Universal breakdown of Kibble-Zurek scaling in fast quenches across a phase transition
• (Su et al., 2013) Kibble-Zurek scaling and its breakdown for spontaneous generation of Josephson vortices in Bose-Einstein condensates
• (Zhang et al., 12 May 2025) Observation of Near-Critical Kibble-Zurek Scaling in Rydberg Atom Arrays
• (Wang et al., 19 May 2025) Tricritical Kibble-Zurek Scaling in Rydberg Atom Ladders
• (Chandran et al., 2012) The Kibble-Zurek Problem: Universality and the Scaling Limit
• (Deffner, 2017) Kibble-Zurek scaling of the irreversible entropy production
• (Yin et al., 2016) Kibble-Zurek scaling in the Yang-Lee edge singularity
• (Bu et al., 2019) Kibble-Zurek Scaling in a Holographic p-wave Superconductor
• (Garcia et al., 2024) The quantum Kibble-Zurek mechanism: the role of boundary conditions, endpoints and kink types
• (Pirozzi et al., 17 Mar 2026) Kibble-Zurek Mechanism in the Open Quantum Rabi Model
• (Nazé et al., 2022) Kibble-Zurek scaling from linear response theory
• (Kou et al., 2024) Kibble-Zurek scaling immune to anti-Kibble-Zurek behavior in driven open systems at the limit of loss difference
• (Ulm et al., 2013) Observation of the Kibble-Zurek scaling law for defect formation in ion crystals
• (Bácsi et al., 2022) Kibble-Zurek scaling due to environment temperature quench in the transverse field Ising model