Coherent Freeze-Out Mechanism

• Coherent Freeze-Out Mechanism is a non-equilibrium framework that tracks slow, collective modes influencing the freeze-out dynamics in both QCD and cosmological settings.
• It employs extended hydrodynamic and Boltzmann methods to match evolving field correlators to observable fluctuations, ensuring accurate predictions.
• The approach reveals how conservation laws and critical slowing down suppress fluctuation observables and modify dark matter relic abundances compared to standard models.

The coherent freeze-out mechanism describes a set of frameworks across diverse physical contexts in which the non-equilibrium evolution of long-wavelength, coherent, or collective modes substantially modifies the standard freeze-out process. In contrast to conventional treatments—where local equilibrium is assumed and fluctuations are neglected or treated perturbatively—coherent freeze-out procedures consistently track the out-of-equilibrium evolution and imprint of slow, coupled degrees of freedom through to observable final states. This is essential, for example, near the QCD critical point in heavy-ion collisions, and in cosmological scenarios where WIMPs interact with a light field such as an ALP. The mechanism requires a formalism that matches non-equilibrium statistical information from hydrodynamic or mean-field evolution to the kinetic and fluctuation observables measured after freeze-out. Recent developments have demonstrated that coherent freeze-out is necessary to correctly predict the magnitude and structure of fluctuation observables, as well as to accurately compute relic abundances where standard thermal freeze-out assumptions fail (Pradeep et al., 2021, Ferrante et al., 20 Nov 2025, Pradeep et al., 2022).

1. Key Features of the Coherent Freeze-Out Concept

The central tenet of the coherent freeze-out mechanism is that slow, critical, or coherent modes—which evolve on longer timescales than the microscopic relaxation—retain out-of-equilibrium imprints that persist across freeze-out and influence observable states. In the QCD context, this refers to critical entropy-per-baryon fluctuations in an expanding quark-gluon plasma; in cosmological applications, to an ALP field coherently coupled to a WIMP density.

A distinguishing feature of coherent freeze-out is the matching of the full non-equilibrium state (encoded in correlators or field configurations) onto a kinetic description at freeze-out. This match is accomplished via explicit prescriptions, such as introducing a critical σ-field with fluctuations determined by the pre-freeze-out hydrodynamics, or by evolving the effective potential and coupled Boltzmann equations in extended field–particle systems.

In all cases, baryon number or analogous global conservation laws induce suppression and memory effects that prevent naive equilibrium predictions from holding—most notably, suppressing the divergence of fluctuation measures as the correlation length grows.

2. Hydro+, Critical Fluctuations, and Coherent Freeze-Out at the QCD Critical Point

Hydro+ formalism extends conventional relativistic hydrodynamics by including deterministic evolution equations for correlators of slow, non-hydrodynamic modes. Specifically, the slow entropy-per-baryon fluctuation $\hat s(x) \equiv s(x)/n(x)$ is characterized via its Wigner-transformed two-point function $\phi_Q(x)$ (Pradeep et al., 2022, Pradeep et al., 2021).

The Hydro+ equations combine:

• Hydrodynamic conservation laws for energy–momentum and baryon number,
• Non-equilibrium mode evolution for $\phi_Q(x)$:

$$u^\mu(x)\partial_\mu \phi_Q(x) = -\Gamma(Q; \xi(x)) \left[ \phi_Q(x) - \bar\phi_Q(\xi(x)) \right],$$

where $\bar\phi_Q(\xi)$ is the local equilibrium Ornstein–Zernike form and $\Gamma(Q;\xi)$ is a critical slowing-down relaxation rate with $\Gamma(Q\to 0) \to 0$ due to baryon number conservation.

At freeze-out, these non-equilibrium fluctuations are mapped onto the correlator of a scalar σ-field, which is minimally coupled to hadronic species, shifting their masses and thereby inducing correlated fluctuations in the observed multiplicities. The generalized Cooper–Frye prescription incorporates the σ-induced corrections to the particle distribution function:

$$f_A(x,p) = f_A^0(p\cdot u(x)) + g_A \frac{\partial f_A^0}{\partial m_A} \sigma(x),$$

where $g_A$ is the coupling of the species $A$ to the σ-field.

Cumulants of hadron number are computed as surface integrals weighted by the σ-correlators, which are fixed by the pre-freeze-out $\phi_Q(x)$ and thermodynamic prefactors. Notably, the out-of-equilibrium suppression of long-wavelength fluctuations (i.e., that $\phi_Q(x) < \bar\phi_Q(\xi(x))$ for modes with $Q\xi \lesssim 1$) leads to a reduction of fluctuation observables—e.g., proton variance—relative to equilibrium expectations and diminishes sensitivity to the freeze-out temperature (Pradeep et al., 2021, Pradeep et al., 2022).

3. Boltzmann Mean-Field Coupling and Dark Matter Coherent Freeze-Out

In cosmological settings, coherent freeze-out emerges in coupled WIMP–ALP models where mean-field and forward-scattering effects mediate strong coupling between otherwise weakly interacting sectors (Ferrante et al., 20 Nov 2025). The model involves a WIMP χ (mass $m_\chi$) and a light ALP φ (mass $m_\phi$), with leading quadratic interaction:

$${\cal L} \supset \tfrac12 (\partial_\mu\phi)^2 - \tfrac12 m_\phi^2 \phi^2 + \tfrac12 \bar\chi i\!\!\not\!\partial \chi - m_\chi \bar\chi\chi + \frac{1}{2\Lambda} \bar\chi\chi\,\phi^2,$$

where $\Lambda$ is a large cutoff scale.

The thermal WIMP bath induces temperature-dependent mass shifts for both species:

$$m_\phi^2(T) = m_\phi^2 - \frac{\langle\bar\chi\chi\rangle_T}{\Lambda}, \quad m_\chi^{\rm eff}(\phi) = |m_\chi - \frac{\phi^2}{2\Lambda}|.$$

High-temperature symmetry breaking displaces the ALP field from the origin. Upon cooling, symmetry is restored via a first-order phase transition or crossover, depending on the coupling parameter $\kappa$. During the high-$T$ broken phase, the WIMP effective mass is reduced and the freeze-out is delayed to lower temperatures, with the WIMP staying semi-relativistic ($m_\chi^{\rm eff}/T \approx 1.33$ for $\kappa\ll1$). At the transition, the WIMP mass snaps back and standard freeze-out proceeds, but at a much lower $T$, greatly enhancing the relic yield compared to the standard thermal scenario.

Coherent freeze-out here directly results from the coupled, non-thermal evolution of the ALP and WIMP sectors, with non-equilibrium symmetry restoration dynamics crucial for determining the final relic abundance—permitting s- or p-wave WIMPs with much larger annihilation cross sections than conventionally allowed (Ferrante et al., 20 Nov 2025).

4. Matching Non-Equilibrium Fields to Observable Fluctuations

The coherent freeze-out prescription in both heavy-ion and cosmological contexts connects a non-equilibrium, coupled fluid or field evolution to observable quantities. This is achieved by:

• Identifying the relevant slow or collective mode (e.g., $\hat s$ near the QCD critical point, the homogeneous ALP field in the early Universe),
• Evolving its dynamics using extended hydrodynamics or coupled Boltzmann equations,
• Matching these dynamics to an effective field (σ or φ) whose fluctuations or occupation numbers directly modulate observables,
• Computing fluctuation observables (cumulants, relic density) as functionals of the evolved non-equilibrium correlators or field values.

For critical fluctuations, higher-order cumulants such as skewness and kurtosis can be computed by integrating the corresponding higher-point correlators of the effective field over the freeze-out (kinetic) hypersurface, employing the thermodynamic matching prescribed by universality and the equilibrium–non-equilibrium map (Pradeep et al., 2022, Pradeep et al., 2021).

5. Phenomenological Consequences and Suppression Effects

A key universal prediction of coherent freeze-out frameworks is the suppression of fluctuation observables—relative to equilibrium expectations—due to memory and conservation effects. In the QCD scenario, baryon number conservation results in the zero mode of critical fluctuations being frozen, despite the equilibrium correlator diverging as the correlation length grows. This leads to finite, reduced cumulants for proton number fluctuations and others, and robust predictions that are less sensitive to freeze-out temperature (Pradeep et al., 2022, Pradeep et al., 2021).

In the dark matter context, the delayed freeze-out dynamics and phase transitions can significantly modify the allowed parameter space for thermal relics, admitting larger cross sections and masses while still matching cosmological data. In regimes where the phase transition is a crossover, the coupled field dynamics can yield a relic ALP abundance that is independent of the initial displacement ("ALP miracle"), set only by the coupling and scale $\Lambda$ (Ferrante et al., 20 Nov 2025).

6. Practical Implementation and Comparison to Standard Freeze-Out

The coherent freeze-out procedure is designed for practical implementation on top of existing hydrodynamic or Boltzmann codes:

• In heavy-ion simulations, one solves Hydro+ or equivalent extended hydrodynamics, evolves the fluctuation spectrum, and computes cumulants via generalized Cooper–Frye integrals using the mapped σ-correlators.
• In cosmological settings, one solves coupled Boltzmann and mean-field evolution equations to extract the late-time relic densities and fluctuation spectrum.

Quantitative analysis demonstrates that the coherent freeze-out procedure captures non-equilibrium suppression and predicts observable quantities ready for comparison with experimental data, such as Beam Energy Scan (BES) measurements of higher-order cumulants in heavy-ion collisions (Pradeep et al., 2022, Pradeep et al., 2021), or indirect and direct DM detection constraints in the WIMP–ALP paradigm (Ferrante et al., 20 Nov 2025).

Table: Key Distinctions Between Standard and Coherent Freeze-Out

Aspect	Standard Freeze-Out	Coherent Freeze-Out
Fluctuations	Assumes equilibrium; neglected	Non-equilibrium, coupled, evolved
Conservation Effects	Ignored or approximate	Exactly encoded, enforce suppression
Observable Sensitivity	Strong dependence on $T_f$	Suppressed, robust to $T_f$ choices
Higher Cumulants	Not reliably computed	Derived from field correlators
DM Relic Cross Section	Narrowly fixed ("WIMP miracle")	Enlarged parameter space via delay

7. Interrelations and Broader Context

Coherent freeze-out mechanisms bridge theoretical approaches across systems marked by the interplay of slow collective modes and rapidly changing environments. They unify fluctuation observables in QCD matter with the dynamics of dark sector relic formation. The adoption of Hydro+ and mean-field extended Boltzmann formalisms, together with explicit thermodynamic matching at freeze-out, represents a substantial advance in the quantitative description of fluctuation and abundance observables in both laboratory and cosmological settings (Ferrante et al., 20 Nov 2025, Pradeep et al., 2022, Pradeep et al., 2021).