First-Order Phase Transitions

• First-order phase transitions are defined by discontinuous jumps in first derivatives (e.g., entropy, volume) and involve latent heat release.
• They are characterized by nucleation processes where droplets of the stable phase grow and coalesce, affecting remnant dynamics and scaling behavior.
• These transitions occur in various systems from lattice models and fluids to quantum and cosmological phenomena, each with distinct experimental signatures.

A first-order phase transition (FOPT) is a thermodynamic transformation between distinct phases of matter characterized by the coexistence of two (or more) phases and a discontinuity in at least one first derivative of the thermodynamic potential—typically in entropy, volume, or order parameter. Such transitions involve latent heat and nucleation processes, and encompass phenomena ranging from the boiling of water, magnetic and structural transitions in solids, and fluid–solid transitions, to cosmological events and quantum many-body systems.

1. Thermodynamic Definition and Classification

A FOPT is defined by a nonanalyticity in the equilibrium free energy such that a first derivative (e.g., entropy $S=-\partial G/\partial T$, volume $V=\partial G/\partial P$, or magnetization $M=-\partial G/\partial H$) jumps discontinuously at the transition. The coexistence point is determined by equality of the relevant thermodynamic potentials between phases, such as $G_1(P_c,T_c)=G_2(P_c,T_c)$ for two phases with respective Gibbs free energies. The corresponding entropy and volume jump obey

$$\Delta S = S_2 - S_1 = -\left(\frac{\partial G}{\partial T}\right)_P^{2}+\left(\frac{\partial G}{\partial T}\right)_P^{1}\,, \quad \Delta V = V_2 - V_1 = \left(\frac{\partial G}{\partial P}\right)_T^{2}-\left(\frac{\partial G}{\partial P}\right)_T^{1}$$

and are related by the Clausius–Clapeyron equation

$\frac{dP_c}{dT} = \frac{\Delta S}{\Delta V}$

(Henry et al., 2017). FOPTs commonly display latent heat $L = T_c\,\Delta S$.

Importantly, the manifestation of discontinuity (jump versus continuous evolution) depends on the number of constrained extensive variables relative to the number of coexisting phases; only for $n<k$ does a strict discontinuity and a well-defined latent heat emerge in the global variables, where $n$ is the number of fixed extensive variables and $k$ the number of coexisting phases (Hempel, 18 Feb 2026).

2. Nucleation, Growth, and Remnant Dynamics

During a FOPT, the transition advances by nucleation of droplets ("bubbles") of the thermodynamically stable phase within the metastable phase, typically described by a nucleation rate per unit volume per unit time

$V=\partial G/\partial P$0

where $V=\partial G/\partial P$1 is a bounce action determined from instanton or critical droplet solutions (Lu et al., 2022). The stable phase bubbles expand at a wall velocity $V=\partial G/\partial P$2 and eventually percolate, fragmenting the metastable phase into disconnected shrinking pockets, or "old-phase remnants." These remnants collapse under the pressure of the stable phase, with their size and number distributions determined by integrating the time-reversed nucleation process. The analytic formalism for their statistics—including explicit formulas for remnant nucleation rates and size distributions—has been developed for general cosmological and statistical FOPTs and is essential for modeling phenomena such as primordial black hole or soliton production (Lu et al., 2022).

The entire domain-growth, coalescence, and wall dynamics can be highly nontrivial, especially when strong metastabilities and energy barriers lead to long-lived, meta-stable or critical patterns in driven or frustrated systems, as observed in numerical holographic simulations and nano-resolved experiments (Bellantuono et al., 2019, Choi et al., 2014). Remnant statistics and their coalescence history directly impact the outcome of cosmological phase transitions and condensed matter pattern formation.

3. Universal Scaling, Renormalization-Group Theory, and Kinetics

While FOPTs are traditionally associated only with nonuniversal nucleation/barrier physics, research in the Ising model, Langevin-class field theories, and related microscopic models has demonstrated "universal scaling" in metastable and dynamically driven FOPTs. The core theoretical structure is a coarse-grained Landau–Ginzburg theory expanded about the mean-field spinodal, leading to a cubic (φ³) effective free energy at the instability, with associated "spinodal" exponents (e.g., $V=\partial G/\partial P$3, $V=\partial G/\partial P$4, $V=\partial G/\partial P$5 for 2D Ising) (Zhang et al., 31 May 2025, Zhong, 2024).

Key universal scaling forms are:

• Static finite-size scaling: $V=\partial G/\partial P$6
• Finite-time scaling (FTS) under a field ramp $V=\partial G/\partial P$7: $V=\partial G/\partial P$8, where $V=\partial G/\partial P$9

Empirically, such scaling and data collapse have been observed in microscopic simulations over broad parameter regimes, unifying FOPTs with continuous transitions under a broader renormalization-group umbrella and yielding universal scaling functions similar to those at continuous critical points (Zhang et al., 31 May 2025, Zhong, 2024, Zhong, 2018, Zhong, 2017). Classical nucleation-and-growth kinetics remain valid as corrections away from the scaling regime.

4. Models, Mechanisms, and Realizations

FOPTs occur across a broad class of systems and models:

a. Lattice Systems and Fluids: Standard paradigms include the Potts model (strong first-order transitions in $M=-\partial G/\partial H$0 in 2D), Bell–Lavis triangles, associating lattice gases, and athermal hard-core lattice gases (Fiore et al., 2012). Semi-analytic approaches leveraging partition function decomposition into phase branches allow precise determination of coexistence properties, order parameter discontinuities, and locations of phase boundaries, with small-system Monte Carlo input (Fiore et al., 2012).

b. Liquid–Liquid and Crystalline Transitions: FOPTs occur not only in solid–liquid or solid–solid transitions, but also between liquid phases. In sulfur, a sharp first-order liquid–liquid transition driven primarily by density is observed, marked by a 7% density jump, abrupt disappearance of molecular Raman features, and vanishing of mid-range XRD correlations—directly confirming a change in the first derivative of the Gibbs free energy (Henry et al., 2017). Similar mechanisms and diagnostics apply in other single-component liquids.

c. Continuum Systems: FOPTs in the continuum, such as liquid–gas criticality in continuum Gibbs point processes, can be rigorously established using Pirogov–Sinai–Zahradník theory, provided the interactions satisfy a locality "saturation" condition and a Peierls-type bound; examples include the Quermass model and diluted pairwise interactions (Dereudre et al., 11 Feb 2026).

d. Ionization/Dissociation-Driven Transitions: A distinct class of FOPTs is driven by chemical composition changes, such as discontinuous jumps in ionization or molecular dissociation (as in the plasma phase transition in dense hydrogen). These transitions have unique thermodynamic and topological features: S-shaped three-valued $M=-\partial G/\partial H$1 isotherms, isolated metastable segments, negative coexistence-curve slopes, and small latent heats—distinct from van der Waals–type transitions (Norman et al., 1 Sep 2025).

e. Quantum Phase Transitions: In quantum systems, FOPTs manifest as discontinuities in the first derivative of the ground-state energy per particle as a function of a control parameter (e.g., a coupling constant). These can be framed as "condensation in state space": as the parameter is tuned, the ground state collapses from spreading over the full Hilbert space to a vanishingly small subspace, generating a jump in $M=-\partial G/\partial H$2 (Ostilli et al., 2017).

5. Experimental and Cosmological Implications

FOPTs are ubiquitous in condensed matter, quantum simulators, and cosmology:

• Nanoscopic Dynamics and Wetting: Direct imaging of the FeRh magneto-structural transition reveals that strong interfacial frustration and wetting layers can slow or freeze