---
title: Gelmini-Gleiser Bubble Kinetics
url: https://www.emergentmind.com/topics/gelmini-gleiser-subcritical-bubble-kinetics
type: topic
---

# Gelmini-Gleiser Bubble Kinetics

Searching arXiv for the specified paper and closely related work on subcritical bubble kinetics and weak first-order phase transitions.
Gelmini-Gleiser subcritical bubble kinetics is a kinetic framework for describing the production, evolution, and disappearance of thermally generated broken-phase fluctuations with radius below the critical radius during a weakly first-order thermal phase transition. In the formulation analyzed in "Subcritical bubble prehistory in weak first-order phase transition" [2605.24891], the framework is used to test the self-consistency of the standard assumption that critical bubbles nucleate on a homogeneous symmetric-phase background. The central result is that, in sufficiently weak transitions, subcritical bubbles can already occupy a percent-level volume fraction by the standard nucleation temperature $T_n$, so the background is better regarded as apparently mixed rather than homogeneous [2605.24891].

## 1. Physical setting and definition of subcritical bubbles

In a weakly first-order thermal phase transition, the finite-temperature potential is taken to be
$$
V(\phi,T)=D\,(T^2-T_0^2)\phi^2-E\,T\,\phi^3+\frac{\lambda}{4}\phi^4.
$$
For temperatures in the interval $T_1>T>T_c$, this potential develops two minima, a symmetric minimum at $\phi_s=0$ and a broken minimum at $\phi_b>0$, separated by a barrier. Here $T_1$ is the temperature at which the nonzero extrema first appear, and $T_c$ is defined by the degeneracy condition $V(0,T_c)=V(\phi_c,T_c)$ [2605.24891].

Within this setting, subcritical bubbles are compact thermal fluctuations of the broken phase generated inside the symmetric phase whose radius $R$ remains below the critical radius $R_c(T)$ defined by the usual $O(3)$ bounce. Because they do not exceed $R_c$, these configurations collapse back to $\phi_s$ rather than grow into expanding true-vacuum domains. They are therefore transient objects, continuously created by thermal noise and destroyed by shrinkage, reverse fluctuations, or thermal agitation [2605.24891].

This distinction is essential. A subcritical bubble is not a critical bubble that has begun successful nucleation; it is a fluctuation that remains below the growth threshold. The relevance of Gelmini-Gleiser kinetics is precisely that a large population of such transient objects can modify the assumed background before standard critical-bubble nucleation becomes efficient.

## 2. Birth-death kinetics in an FRW background

Let $n(R,t)\,dR$ denote the comoving number density of subcritical bubbles with physical radius in the interval $[R,R+dR]$ at cosmic time $t$. In an FRW universe with Hubble rate $H(t)$, the Gelmini-Gleiser evolution equation is
$$
\frac{\partial n}{\partial t}+3H\,n
=-\frac{\partial}{\partial R}\left[\left(\frac{dR}{dt}\right)n\right]
+p_s^{\rm sub}(T)\,\Gamma_{s\to b}^{\rm sub}(R,T)
-\Gamma_{b\to s}^{\rm sub}(R,T)\,n
-\Gamma_{\rm TN}^{\rm sub}(R,T)\,n.
$$
The shrinkage law is
$$
\frac{dR}{dt}=-v,
$$
with $v\approx 1$ in units $c=1$. The factor $p_s^{\rm sub}(T)=e^{-f_b(T)}$ is the fraction of volume still in the symmetric phase. The source term $\Gamma_{s\to b}^{\rm sub}$ creates broken-phase fluctuations of radius $R$ inside the symmetric background, while $\Gamma_{b\to s}^{\rm sub}$ and $\Gamma_{\rm TN}^{\rm sub}$ erase them through reverse fluctuations and thermal noise, respectively [2605.24891].

Using a Gaussian profile
$$
\phi_b^{\rm sub}(r)=\phi_b(T)e^{-r^2/R^2},
$$
the free energy of a broken-phase subcritical fluctuation is
$$
F_b(R,T)=A_\nabla(T)\,R+B_b(T)\,R^3,
$$
with
$$
A_\nabla(T)=\frac{3\pi^{3/2}}{4\sqrt{2}}\,\phi_b^2(T),
$$
and
$$
B_b(T)=\pi^{3/2}\left[\frac{D(T^2-T_0^2)\phi_b^2}{2^{3/2}}-\frac{E\,T\,\phi_b^3}{3^{3/2}}+\frac{\lambda\,\phi_b^4}{4^{5/2}}\right].
$$
The corresponding source rate is
$$
\Gamma_{s\to b}^{\rm sub}(R,T)
=w_b(T)\,A_{\rm sc}\,T^4\,e^{-F_b(R,T)/T}\,
\Theta(R-\xi_s)\Theta(R_c-R),
$$
where
$$
\xi_s=[V''(0,T)]^{-1/2},
$$
$$
w_b(T)=\left[\frac{\phi_b-\phi_{\rm top}}{\phi_b}\right]^2\Theta(T_1-T),
$$
and
$$
R_c=\sqrt{\frac{-A_\nabla}{3B_b}}
\quad \text{if } B_b<0,
$$
with $R_c\to\infty$ otherwise.

For the reverse process, the profile
$$
\phi_s^{\rm sub}(r)=\phi_b[1-e^{-r^2/R^2}]
$$
gives
$$
F_s(R,T)=A_\nabla(T)\,R+B_s(T)\,R^3,
$$
and
$$
\Gamma_{b\to s}^{\rm sub}(R,T)
=A_{\rm sc}\,T^4\,e^{-F_s(R,T)/T}\,
\Theta(R-\xi_b)\Theta(R_{c,s}-R),
$$
with
$$
\xi_b=[V''(\phi_b,T)]^{-1/2},
$$
and
$$
R_{c,s}=\sqrt{\frac{-A_\nabla}{3B_s}}
\quad \text{when } B_s<0.
$$
The thermal-noise erasure rate is
$$
\Gamma_{\rm TN}^{\rm sub}= \frac{a\,T}{4\pi R^3/3}.
$$
This structure makes the kinetics explicitly a competition between fluctuation production, collapse, reverse conversion, and Hubble dilution [2605.24891].

## 3. Subcritical volume fraction and the single-bin approximation

The geometric broken-phase volume fraction associated with the radius distribution is
$$
f_b(T)=\int_0^\infty dR\,\frac{4\pi R^3}{3}\,n(R,t).
$$
To incorporate random overlap when $f_b$ is not asymptotically small, the Poisson expression
$$
p_b^{\rm sub}(T)=1-e^{-f_b(T)}
$$
is used for the actual broken-phase volume fraction, so that
$$
p_s^{\rm sub}(T)=e^{-f_b(T)}
$$
is the symmetric-phase volume fraction [2605.24891].

A simplified treatment replaces the full distribution by a single representative bin at $R=\xi_s$. Writing
$$
\mu(t)=\int dR\,n(R,t)\approx n(\xi_s,t)\,\Delta R,
$$
and neglecting boundary flux, one obtains
$$
\frac{d\mu}{dt}=S_{\rm tot}(T)-K(\xi_s,T)\,\mu,
$$
with
$$
S_{\rm tot}\simeq \Gamma_{s\to b}^{\rm sub}(\xi_s,T)\,\Delta R,
$$
and
$$
K(\xi_s,T)=3H+a\,T+\left(\frac{4\pi \xi_s^3}{3}\right)\left[\Gamma_{b\to s}^{\rm sub}(\xi_s,T)+\Gamma_{\rm TN}^{\rm sub}(\xi_s,T)\right].
$$
In quasi-steady state, $d\mu/dt\approx 0$, so
$$
\mu\simeq \frac{S_{\rm tot}}{K}.
$$
The resulting fast estimate for the subcritical volume fraction is
$$
\hat f_\xi(T)\simeq \left(\frac{4\pi \xi_s^3}{3}\right)\frac{S_{\rm tot}(\xi_s,T)}{K(\xi_s,T)}.
$$

The single-bin approximation is not presented as a replacement for the full kinetic equation in all regimes. Rather, it is a compact estimator designed to identify parameter regions in which subcritical occupancy is large enough to threaten the homogeneous-background approximation. A plausible implication is that the approximation is most useful as a scan-level diagnostic, with the full $n(R,t)$ evolution reserved for validation near the boundary.

## 4. Relation to standard critical-bubble nucleation

The standard critical-bubble rate is expressed through the three-dimensional Euclidean action
$$
S_3(T)=4\pi\int dr\,r^2\left[\frac{1}{2}(\partial_r\phi)^2+V(\phi,T)-V(0,T)\right],
$$
with nucleation rate
$$
\Gamma_c(T)=A_c\,e^{-S_3/T},
$$
where $A_c\approx T^4$ [2605.24891].

The expected number of critical bubbles per Hubble volume is
$$
N_c(T)=\int_T^{T_c}\frac{dT'}{T'}\,\frac{\Gamma_c(T')}{H^4(T')},
$$
and the nucleation temperature $T_n$ is defined by
$$
N_c(T_n)=1.
$$
The fraction of space converted by critical bubbles is
$$
p_{\rm crit}(T)=1-\exp[-I_c(T)],
$$
where
$$
I_c(T)=\int_T^{T_c}\frac{dT'}{T'}\,\Gamma_c(T')\left(\frac{a(T')}{a(T)}\right)^3\frac{4\pi}{3}R(T,T')^3.
$$

The conceptual comparison is then straightforward. Standard analyses assume that before $T_n$ the system remains effectively in a homogeneous symmetric vacuum, with critical bubbles nucleating on top of that background. Gelmini-Gleiser kinetics tests that premise by evolving the subcritical population from a starting temperature $T_{\rm start}>T_c$ down to $T_n$ and evaluating $p_b^{\rm sub}(T_n)$. The homogeneous-background approximation is classified as valid when
$$
p_b^{\rm sub}(T_n)\lesssim 10^{-3},
$$
with small subcritical corrections for
$$
10^{-3}\lesssim p_b^{\rm sub}(T_n)\lesssim 10^{-2},
$$
and breakdown of the homogeneous-background assumption when
$$
p_b^{\rm sub}(T_n)\gtrsim 10^{-2}.
$$
In the latter case, the system is categorized as a mixed background rather than an ordinary homogeneous bounce point [2605.24891].

A frequent simplification in phase-transition phenomenology is to treat all pre-nucleation fluctuations as negligible if they are subcritical. The analysis here shows that this is not generically self-consistent in weak transitions: subcriticality prevents indefinite growth, but not macroscopic occupancy.

## 5. Numerical criterion for mixed-background candidates

A central quantitative result is a simple diagnostic based on the fast estimator $\hat f_\xi(T_n)$. At the nucleation temperature,
$$
\hat f_\xi(T_n)\simeq \left(\frac{4\pi \xi_s^3(T_n)}{3}\right)
\frac{\Gamma_{s\to b}^{\rm sub}(\xi_s,T_n)\,\Delta R}
{3H_n+a\,T_n+\left(\frac{4\pi \xi_s^3}{3}\right)\left(\Gamma_{b\to s}^{\rm sub}+\Gamma_{\rm TN}^{\rm sub}\right)}.
$$
Across a scan of weak transitions in which $(D,\lambda,r_c)$ are varied and $T_n$ is determined, the comparison between the fast estimate and the full kinetic result for points with $f_b<10^{-2}$ yields the fit
$$
\log_{10} f_b(T_n)\simeq 0.20+1.13\,\log_{10}\hat f_\xi(T_n).
$$
Imposing $f_b=10^{-2}$ gives
$$
1.13\,\log_{10}\hat f_\xi(T_n)+0.20=-2,
$$
and therefore
$$
\log_{10}\hat f_\xi(T_n)\approx -1.95.
$$
This leads to the practical rule
$$
\log_{10}\hat f_\xi(T_n)\gtrsim -1.95
\quad \Rightarrow \quad
p_b^{\rm sub}(T_n)\gtrsim 1\%.
$$
Points above this boundary are flagged as mixed-background candidates rather than ordinary homogeneous bounce points [2605.24891].

The significance of this result is operational. It replaces a full kinetic evolution with a compact threshold test that can be embedded in parameter scans. This suggests a two-stage workflow: fast classification through $\hat f_\xi(T_n)$, followed by direct solution of the kinetic equation in marginal or phenomenologically important regions.

## 6. Parameter dependence, interpretation, and scope

The parameter scan identifies the regime in which sizable subcritical volume fractions arise at $T_n$. These occur when the two phases are nearly degenerate, the barrier is low, the free-energy difference between phases is moderate, and the transition is weak [2605.24891]. More specifically, large subcritical fractions are found when:

- **Small free-energy splitting**: $\Delta V_{ft}=V(0,T_n)-V(\phi_b,T_n)$ is small relative to $T_n^4$, so broken patches have a low volume-energy cost.
- **Low barrier height**: $\Delta V_{bar}=V(\phi_{\rm top},T_n)-V(0,T_n)$ is low, reducing Boltzmann suppression.
- **Small order parameter**: $\phi_n/T_n$ is small, which reduces the gradient term because $A_\nabla\propto \phi_b^2$.
- **Weak transition strength**: $\alpha_n=\Delta\rho/\rho_{\rm rad}\lesssim 10^{-3}$, so cooling is slow and the system remains near $T_c$ for longer.

Conversely, stronger transitions with large $\alpha_n$, large $\phi_n/T_n$, and a large barrier suppress the subcritical population. The resulting interpretation is narrowly targeted but consequential: the issue is not whether critical-bubble theory fails in general, but whether the assumed prehistory of the background remains homogeneous in weak first-order transitions. In that sense, Gelmini-Gleiser kinetics functions as a consistency test for the usual bounce-based treatment rather than a rejection of it.

The scope is similarly specific. The analysis is formulated for the finite-temperature quartic potential given above, employs Gaussian subcritical profiles, and uses both the full birth-death equation and the single-bin estimate. Within that setup, the conclusion is that one must evaluate $\hat f_\xi(T_n)$, or solve the full $n(R,t)$ kinetics, in model scans whenever weak transitions are present, because some points conventionally treated as standard nucleation events should instead be regarded as mixed-background configurations [2605.24891].

Source: https://www.emergentmind.com/topics/gelmini-gleiser-subcritical-bubble-kinetics