Supercooled First-Order Phase Transition

• Supercooled FOPT is a thermodynamic process where systems remain in a metastable phase well below equilibrium until nucleation triggers a rapid phase change.
• The transition is governed by nucleation dynamics, with critical parameters like ΔT, α, and β/H dictating the timing, bubble growth, and phase conversion rates.
• Applications span condensed matter, cosmology, and quantum field theory, leading to observable effects such as gravitational waves, primordial black holes, and altered axion dynamics.

A supercooled first-order phase transition (FOPT) is a thermodynamic process where a system remains trapped in a metastable phase well below the equilibrium transition temperature or critical point, with the actual phase conversion—mediated by nucleation and growth of bubbles or domains of the stable phase—initiating only after considerable undercooling or supersaturation. Such transitions are fundamental to a broad range of condensed matter, cosmological, and quantum field-theoretic settings, with recent developments clarifying how supercooling qualitatively affects dynamical, statistical, and observable features of the phase transition.

1. Phenomenology and Classification

A first-order phase transition is characterized by discontinuous changes in order parameters and the release or absorption of latent heat. In the context of supercooling, the system’s temperature, pressure, or chemical potential is driven far into the metastable regime before nucleation becomes probable. The degree of supercooling is quantified by the difference ΔT = T_c − T_n, where T_c is the temperature at which the free energies of the two phases are degenerate and T_n is the nucleation (or percolation) temperature at which the rate of critical bubble formation per Hubble volume per Hubble time (for cosmological settings) or per characteristic system volume crosses unity (Wang et al., 2020).

The strength of the phase transition can be classified by the α parameter—the ratio of latent heat (vacuum energy released) to the background energy density at the transition: $$\alpha = \frac{\Delta \rho_{\rm vac}}{\rho_{\rm rad}(T_{*})}$$ with thresholds such as slight supercooling (α ≲ 0.1), mild (0.1 ≲ α ≲ 0.5), strong (0.5 ≲ α ≲ 1), and ultra-supercooling (α ≫ 1) (Wang et al., 2020). Cosmological and laboratory systems display similar dynamical features despite system-specific order parameters and microscopic mechanisms.

2. Nucleation, Growth, and Dynamics of Supercooled First-Order Transitions

Nucleation theory governs the onset of phase conversion in supercooled FOPTs. The nucleation rate per unit volume and time takes the form

$$\Gamma(T) \approx T^4 \exp\left[-\frac{S_3(T)}{T}\right]$$

where $S_3(T)$ is the O(3)-symmetric bounce action for the critical bubble, obtained by minimizing the three-dimensional Euclidean action using the appropriate finite-temperature effective potential (Levi et al., 2022, Wang et al., 2020). For cosmological transitions, percolation and completion are possible only if bubble nucleation outpaces cosmic expansion; in the limit of deep supercooling, the Universe can become dominated by vacuum (latent) energy, transiently expanding quasi-de Sitter–like until bubbles nucleate in sufficient density to percolate (Lewicki et al., 19 Nov 2025, Athron et al., 2023).

The mean bubble nucleation timing can be characterized by the nucleation temperature ($T_n$) and a rapidity parameter β, defined as

$$\beta / H \equiv T_* \left. \frac{d}{dT} \left( \frac{S_3}{T} \right) \right|_{T_*}$$

where H is the Hubble parameter. Lower β/H corresponds to more prolonged (slower) transitions with larger bubbles; the distinction between nucleation and percolation temperature ($T_p$) becomes pivotal in ultra-supercooled cases ($\alpha \gtrsim 1$), as the bubbles can only overlap and complete the transition at much later times and lower temperatures compared to the onset of nucleation (Wang et al., 2020).

3. Realizations: Water, Electron Systems, Helium-3, and Cosmology

Supercooled Transitions in Water and Amorphous Ices

Supercooled water provides a canonical example of a FOPT in a metastable regime. Two-state models incorporating LDL-HDL coexistence and ferroelectric ordering capture both the LL critical point and singularities in observables such as heat capacity and dielectric constant near the λ-point (~233 K), with transitions modeled as (Fedichev et al., 2012): $$G(c, s; T, P) = c G_{LDL}(s; T, P) + (1-c) G_{HDL}(T, P) + Uc(1-c) + RT[c\ln c + (1-c)\ln(1-c)]$$ where c is the fraction of LDL, and only LDL undergoes a Landau-type ferroelectric transition. The sharp dielectric anomaly and heat capacity peak observed in supercooled water are reproduced by this mean-field picture, aligning with simulation and experiment (Giovambattista et al., 2021).

Experiments and simulations on ultrafast heating and decompression of amorphous ice directly probe the non-equilibrium remnant of an equilibrium LLPT: rapid decompression of high-density amorphous ice (HDA) heated above its glass transition leads to a sharp, nearly discontinuous HDL→LDL density drop that maps the underlying coexistence and spinodal boundaries, even in strongly out-of-equilibrium conditions (Giovambattista et al., 2021).

Electronic Supercooled Liquid and Glassy Arrest

In certain electron systems (spinless electrons at half-filling on a triangular lattice), the first-order stripe-ordering transition can be supercooled in the absence of disorder. Extended dynamical mean-field theory (EDMFT) yields a first-order transition at a temperature T_s with a coexistence window terminated by spinodal points. The supercooled electronic liquid at T < T_s is metastable, exhibits a deep pseudogap in the local DOS, and displays kinetic arrest (glassy behavior), mirroring experimental observations in ultra-pure θ-(BEDT-TTF) compounds (Rademaker et al., 2015).

Superfluid ³He A→B Transition

The A→B transition in superfluid ³He (complex 3×3 order parameter) is the leading laboratory analogue for supercooled FOPTs relevant to the early universe. Experimental work shows that the A phase can be cooled significantly below the equilibrium A–B boundary, with the degree of supercooling and the locus of the "catastrophe line" sensitive to thermal history and elimination of B-phase seeds (heterogeneous nucleation centers). Homogeneous nucleation rates remain negligible for experimentally accessible timescales (Tian et al., 2022).

4. Supercooled FOPTs in Quantum Field Theory: Cosmological and Beyond-SM Scenarios

Supercooled first-order transitions naturally arise in both weakly coupled (e.g., Coleman–Weinberg) and strongly coupled (Randall–Sundrum, conformal field theories) settings (Levi et al., 2022, 0706.3388). The finite-temperature potential for a generic scalar order parameter φ can be written as

$$V(\phi, T) = \frac{1}{2} m^2(T) \phi^2 - \frac{1}{3} \delta(T) \phi^3 + \frac{1}{4} \lambda(T) \phi^4$$

with temperature-dependent parameters governed by gauge, Yukawa, or higher-dimensional interactions. Supercooling occurs when nucleation is delayed by a high barrier or a nearly flat potential, with completion possible only within a "supercooling window" in parameter space—bounded below by the requirement that bubbles do nucleate and above by the condition that the transition does not complete in the radiation-dominated era (Levi et al., 2022).

Extensions of the Standard Model with approximate classical conformal symmetry often exhibit extreme supercooling: the phase transition may not complete until QCD-scale temperatures, with six quark flavors remaining massless. The QCD chiral phase transition becomes first order, triggering EW symmetry breaking through the top-Higgs Yukawa coupling and manifesting a strongly supercooled, dynamical sequence of FOPTs with distinctive cosmological signatures (Bodeker, 2021, Iso et al., 2017, Guan et al., 2024).

5. Supercooling in Cosmological Phase Transitions: Gravitational Waves, PBHs, and Axion Physics

Gravitational Wave Production

Supercooled FOPTs generically generate stochastic gravitational wave (GW) backgrounds, with amplitudes and spectral peaks highly sensitive to α, β/H, nucleation and percolation dynamics, and the cosmic expansion history. In the strong/ultra-supercooled regime (β/H ≲ 50–100, α ≫ 1), cosmic expansion modifies standard GW observables: frequency scaling becomes sublinear in β, amplitudes are suppressed compared to Minkowski estimates, and a pronounced causality-induced f³ tail appears below the bubble collision horizon scale (Lewicki et al., 19 Nov 2025, Wang et al., 2020).

Notably, supercooled transitions shift the GW peak towards higher frequencies and lower amplitudes than naïvely expected, constraining their detectability in planned interferometers (LISA, DECIGO, BBO) and disfavoring scenarios attempting to explain nHz GW backgrounds in pulsar timing arrays via MeV–GeV scale supercooled FOPTs (Athron et al., 2023).

Formation of Primordial Black Holes

During supercooled FOPTs, stochastic bubble nucleation results in rare, large regions of unconverted false vacuum. The large pressure imbalance between expanding bubble walls and vacuum-dominated patches can cause such regions to collapse into primordial black holes (PBHs) once they re-enter the horizon, with the energy stored in relativistically moving domain walls driving collapse. The PBH mass spectrum is extremely narrow, set by the timing of collapse (T ~ T_PBH), and may account for observed microlensing events in the Earth-mass window (Gonçalves et al., 2024, Flores et al., 2024). The same underlying transition also produces a GW signature and sets constraints for cosmological relics.

QCD Axion Dynamics and Domain Walls

Supercooled FOPTs that reheat the universe above Λ_QCD can alter the cosmological evolution of axion fields. The abrupt change in energy density and temperature during bubble percolation induces kinetic misalignment of the axion and, by spatial inhomogeneities in reheating, can form closed or infinite QCD axion domain walls whose eventual collapse affects the relic density and may seed further PBH formation (Lyu et al., 24 Jun 2025).

6. Nucleation Kinetics and Holographic Approaches

Holographic models of FOPTs using Einstein–scalar systems in AdS backgrounds have allowed fully nonlinear simulation of supercooled transitions. Dynamical studies reveal threshold behavior (critical nuclei) with universal three-stage time evolution: rapid approach to the critical solution, logarithmic dwelling, and exponential departure, with a scaling law τ ∝ −ln|p − p*|. Both local seed perturbations and collision-induced phase separation exhibit the same universality class and critical exponents, with barrier heights and transition rates directly computable in the holographic model (Chen et al., 2022, Chen et al., 2022). Kinetic arrest and glass-like behavior emerge naturally as the system lingers near or below the nucleation barrier.

7. Theoretical and Experimental Constraints, Open Problems

The supercooled FOPT scenario is subject to several key theoretical and experimental constraints:

• Transitions must complete (percolation) within the available time; otherwise, the universe cannot exit from a vacuum-dominated phase, leading to a "graceful-exit" problem (Athron et al., 2023).
• Successful baryogenesis, dark matter genesis, and compatibility with nucleosynthesis impose bounds on the transition strength, reheating temperature, and the details of critical fluctuation spectra (0706.3388).
• Macroscopic signatures include stochastic GW spectra in the mHz–100 Hz band, mass spectrum and abundance of PBHs, and nonthermal shifts in relic axion or dark matter populations. Current and future observational programs (LISA, PTA, microlensing, CMB) provide the relevant experimental probes (Gonçalves et al., 2024, Lewicki et al., 19 Nov 2025, Lyu et al., 24 Jun 2025).

The dependence of completion and observable signatures on model deformations, vacuum structure, and coupling constants is now quantifiable both in weakly-coupled field theories and strongly-coupled holographic duals (Levi et al., 2022). Laboratory analogues such as supercooled ³He and organic charge-ordered compounds further inform the complex interplay between nucleation, disorder, geometry, and kinetic arrest in FOPTs (Tian et al., 2022, Rademaker et al., 2015).

References Table: Representative Systems Exhibiting Supercooled First-Order Phase Transitions

Physical System	Mechanism / Order Parameter	Key Features / Observables
Supercooled water	LDL-HDL fraction, ferroelectric order	λ-point singularity, Cp, ε(T), two-state mean-field
Electron liquids (θ-MM')	Stripe order parameter (charge density)	Metastable pseudogap, glassy arrest, EDMFT
Superfluid ³He A→B	3×3 order parameter, domain-wall pinning	Strong path, seed elimination, reproducible supercooling
Cosmological FOPTs	Scalar vacuum expectation value (φ), bubbles	Gravitational wave spectrum, PBH, baryogenesis
Holographic QFT (AdS/CFT)	Bulk scalar φ, dual order parameter	Nonlinear dynamical nucleation, critical scaling

In sum, supercooled first-order phase transitions are a unifying, quantitatively tractable concept spanning condensed matter, statistical mechanics, and cosmology, with ramifications for gravitational-wave astronomy, dark matter, baryogenesis, and the dynamical emergence of new phases under extreme nonequilibrium conditions.