---
title: 'Kibble-Zurek Freeze-Out: Scaling & Defects'
url: https://www.emergentmind.com/topics/kibble-zurek-freeze-out
type: topic
---

# Kibble-Zurek Freeze-Out: Scaling & Defects

The Kibble-Zurek freeze-out is a foundational concept in non-equilibrium statistical mechanics that describes how universal scaling laws emerge for defect formation and correlation lengths when a system is driven through a continuous (second-order) phase transition at a finite rate. This process is intimately linked to critical slowing down: as the critical point is approached, both the relaxation time and spatial correlation length diverge, making it impossible for the system to maintain adiabatic evolution. The breakdown of adiabaticity produces a “frozen” regime—the freeze-out interval—where dynamics are effectively arrested, and fluctuations become imprinted. This sets the scale for symmetry-breaking domains and topological defects. The freeze-out scenario is quantitatively specified by universal formulas relating defect densities and frozen-in correlation lengths to the driving protocol and underlying critical exponents. Recent research has revealed that these scaling relations universally fail in the limit of rapid quenches, where the freeze-out time and defect density are determined by the final quench depth rather than the quench rate, establishing a sharp dynamical regime change.

## 1. Foundations: Critical Slowing Down and Freeze-Out

Near a continuous phase transition, the system’s relaxation time $\tau(\epsilon)$ and correlation length $\xi(\epsilon)$ diverge algebraically as a function of the reduced distance to criticality, $\epsilon$:
\[
\xi(\epsilon) = \xi_0 |\epsilon|^{-\nu}\,, \qquad \tau(\epsilon) = \tau_0 |\epsilon|^{-\nu z}
\]
with $\nu$ the (static) correlation-length exponent and $z$ the dynamical exponent. When a control parameter $\lambda(t)$ is swept linearly across its critical value $\lambda_c$ (e.g.\ $\lambda(t) = \lambda_c (1 - t/\tau_Q)$, so $\epsilon(t) = t/\tau_Q$), the system’s intrinsic relaxation time eventually exceeds the time remaining to the transition, $\tau[\epsilon(\hat t)] \sim \hat t$, defining the freeze-out time $\hat t$. Solving,
\[
\tau_0 \left( \frac{\hat t}{\tau_Q} \right)^{-\nu z} \sim \hat t \implies
\hat t \sim \left(\tau_0 \tau_Q^{\nu z}\right)^{1/(1+\nu z)}
\]
The “frozen” correlation length at this instant is
\[
\hat \xi = \xi_0 \left| \frac{\hat t}{\tau_Q} \right|^{-\nu} \sim \xi_0 \left( \frac{\tau_Q}{\tau_0} \right)^{\nu/(1+\nu z)}
\]
These relations are universal and apply to both classical and quantum phase transitions, provided the system is homogeneous and the quench is sufficiently slow [2204.13529].

## 2. Defect Density and Scaling Laws in the Freeze-Out Regime

The physical consequence of freeze-out is that domains of size $\sim \hat \xi$ choose their broken-symmetry configuration independently, and defects appear at boundaries between uncorrelated regions. In $d$ dimensions, the defect density then scales as
\[
n_d \sim \hat \xi^{-d} \sim \tau_Q^{-d \nu/(1+\nu z)}
\]
This regime and its power-laws have been confirmed across a wide swath of experimental platforms: ultracold gases crossing BEC transitions (3D, 2D, and annular geometries) [1611.01145], colloidal and spin systems [1303.6821, 1503.08698], cold atomic self-organization [1610.00917], superconducting and quantum electronic models [1505.07362, 2203.12438, 2412.01436], extended Gross–Pitaevskii/BEC models [2411.18395], and Ginzburg–Landau simulations for symmetry-breaking in superfluids and chiral condensates [2412.03544, 2511.22583].

In inhomogeneous systems, the theory generalizes by introducing a local effective quench time $\tau_Q(\vec r)$ and local freeze-out scales [2010.04374, 1301.4770]. The basic exponents and scaling combinations $\nu/(1+\nu z),\ z\nu/(1+\nu z)$ remain universal, but the spatial variation of the critical parameter induces a position-dependent freeze-out boundary (“front”) between nonadiabatic and adiabatic regions.

## 3. Fast-Quench Regime: Breakdown of Universal Kibble-Zurek Scaling

Standard Kibble-Zurek scaling applies only for quench times $\tau_Q$ much larger than a critical value $\tau_Q^*$. If the quench ends at a finite value $\epsilon_f$ before the relaxation time has reached the time remaining, the system never re-enters the adiabatic regime and remains frozen until the final parameter value is reached. In this limit, the freeze-out time and correlation length become independent of $\tau_Q$ and depend only on $\epsilon_f$:
\[
t_{\rm fo} \sim \tau_0 |\epsilon_f|^{-\nu z},\qquad
\xi_{\rm fo} \sim \xi_0 |\epsilon_f|^{-\nu}
\]
The defect density “saturates” at a universal plateau controlled by the quench depth [2204.13529]:
\[
n_{\rm plateau} \sim \xi_{\rm fo}^{-d} \propto |\epsilon_f|^{d\nu}
\]
The crossover between Kibble-Zurek and plateau regimes occurs at $\tau_Q^* \propto |\epsilon_f|^{-(1+\nu z)}$. This universal breakdown and crossover behavior have been established both analytically and numerically for classical $\phi^4$ and stochastic Gross–Pitaevskii models, quantum Ising chains, and in experiment with trapped ions and superconducting qubits [2204.13529, 2110.07969, 2506.06841].

## 4. Dynamical Regimes and Crossover Structure

A summary of dynamical regimes is provided in the table below:

| Quench Regime           | Scaling (Freeze-Out Time)                        | Defect Density                   | Dependence         |
|-------------------------|--------------------------------------------------|-----------------------------------|--------------------|
| Slow ($\tau_Q \gg \tau_Q^*$) | $\hat t \sim (\tau_0 \tau_Q^{\nu z})^{1/(1+\nu z)}$        | $n_d \sim \tau_Q^{-d\nu/(1+\nu z)}$     | Universal KZM $\tau_Q$-power law |
| Fast ($\tau_Q \ll \tau_Q^*$) | $t_{fo} \sim \tau_0 |\epsilon_f|^{-\nu z}$                 | $n_{plateau} \sim |\epsilon_f|^{d\nu}$  | Universal plateau (no $\tau_Q$ dependence) |
| Crossover                | Smooth interpolation between above regimes         | $(1 + \tau_Q/K(\epsilon_f))^{-d\nu/(1+\nu z)}$ | Universal function of $\tau_Q/\tau_Q^*$ |

Experimental data confirm the universal plateau and critical-scaling crossovers, including in holographic AdS/CFT realizations of rapid quenches [2110.07969], and trapped-ion and qubit simulators [2506.06841].

## 5. Analytical and Numerical Models: Verification and Extensions

The standard freeze-out scenario has been analytically solved in paradigmatic models. In the classical $\phi^4$ model (mean-field $\nu=1/2, z=2$), classical and quantum Ising chains ($\nu=z=1$), and the stochastic Gross–Pitaevskii equation, the full sequence of KZ scaling for slow quenches and saturation for fast quenches is observed [2204.13529]. Analytical solutions of the overdamped time-dependent Ginzburg–Landau ODE provide direct access to the time-resolved evolution of the order parameter, exhibiting the sharp inflection at $\pm\hat t$ that marks freeze-out [2511.22583]. These solutions further predict experimental signatures such as the super-exponential growth of the order parameter, temporal overlap correlation curves, and domain size scaling, all of which are observable in real-time measurements.

Finite-rate Kibble–Zurek scaling and the breakdown of adiabaticity have also been directly mapped onto the Landau–Zener two-level transition (LZT) framework, demonstrating the formal equivalence between KZM and LZT for two-level systems [1505.07362]. This connection has allowed the freeze-out boundary and associated scaling behaviors to be probed with exceptional temporal resolution in superconducting qubits.

## 6. Generalizations, Universality, and Limitations

The universality of Kibble-Zurek freeze-out has been confirmed over a vast range of systems, symmetry classes, and even out-of-equilibrium initial states [1503.02762]. The introduction of finite-time scaling (FTS) theory allows the freeze-out process to be extended to far-from-equilibrium starting points, showing that universal behavior (defined by underlying equilibrium exponents) persists even in the absence of an initial adiabatic regime.

Limitations of the standard freeze-out picture have been identified in the context of slow annealing [1001.0693]. If the system is given time to coarsen during and after the freeze-out, the actual density of defects observed at late times can fall below the “frozen-in” density predicted by KZM, a result that is quantitatively captured by scaling forms combining adiabatic, critical coarsening, and regular coarsening regimes. Additionally, modifications are necessary in lower-dimensional systems governed by non-power-law (e.g., exponential, KTHNY-type) divergences, and in cases (such as strongly anisotropic or inhomogeneous systems) where spatial dependence alters the freeze-out boundaries [2010.04374, 1303.6821, 1301.4770, 2310.18216].

## 7. Applications and Significance

Kibble-Zurek freeze-out scaling, both in its standard and breakdown regimes, has been instrumental in interpreting experiments across condensed matter, ultracold atomic gases, quantum electronic platforms, and synthetic quantum simulators. Its universality provides a means to extract and verify critical exponents, classify universality classes, and understand the microscopic mechanisms underpinning the generation of non-equilibrium topological defects. The freeze-out paradigm also provides operational protocols for optimizing dynamic control (e.g., minimizing irreversible work [2412.01436]), and for engineering non-equilibrium states with tunable correlation properties.

Recent experimental and analytical confirmation of the universal breakdown—the emergence of a $\tau_Q$-independent plateau in defect density for rapid quenches—has set firm boundaries on dynamical universality, clarifying the interplay between rate, range, and universality of non-equilibrium scaling across phase transitions [2204.13529, 2110.07969, 2506.06841].

Source: https://www.emergentmind.com/topics/kibble-zurek-freeze-out