Gelmini-Gleiser Bubble Kinetics

• The paper demonstrates that subcritical bubbles—transient broken-phase fluctuations—can occupy a percent-level volume at the nucleation temperature, questioning homogeneous background assumptions.
• It employs a kinetic framework based on birth-death dynamics with Gaussian profiles to model bubble evolution, shrinkage, and competing processes under cosmic expansion.
• The analysis introduces a practical single-bin diagnostic to flag mixed-background regimes, offering a compact test for weak first-order phase transitions.

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" (Chen et al., 24 May 2026), 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 (Chen et al., 24 May 2026).

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)$ (Chen et al., 24 May 2026).

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 $$V(\phi,T)=D\,(T^2-T_0^2)\phi^2-E\,T\,\phi^3+\frac{\lambda}{4}\phi^4.$$0 bounce. Because they do not exceed $$V(\phi,T)=D\,(T^2-T_0^2)\phi^2-E\,T\,\phi^3+\frac{\lambda}{4}\phi^4.$$1, these configurations collapse back to $$V(\phi,T)=D\,(T^2-T_0^2)\phi^2-E\,T\,\phi^3+\frac{\lambda}{4}\phi^4.$$2 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 (Chen et al., 24 May 2026).

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 $$V(\phi,T)=D\,(T^2-T_0^2)\phi^2-E\,T\,\phi^3+\frac{\lambda}{4}\phi^4.$$3 denote the comoving number density of subcritical bubbles with physical radius in the interval $$V(\phi,T)=D\,(T^2-T_0^2)\phi^2-E\,T\,\phi^3+\frac{\lambda}{4}\phi^4.$$4 at cosmic time $$V(\phi,T)=D\,(T^2-T_0^2)\phi^2-E\,T\,\phi^3+\frac{\lambda}{4}\phi^4.$$5. In an FRW universe with Hubble rate $$V(\phi,T)=D\,(T^2-T_0^2)\phi^2-E\,T\,\phi^3+\frac{\lambda}{4}\phi^4.$$6, the Gelmini-Gleiser evolution equation is

$$V(\phi,T)=D\,(T^2-T_0^2)\phi^2-E\,T\,\phi^3+\frac{\lambda}{4}\phi^4.$$7

The shrinkage law is

$$V(\phi,T)=D\,(T^2-T_0^2)\phi^2-E\,T\,\phi^3+\frac{\lambda}{4}\phi^4.$$8

with $$V(\phi,T)=D\,(T^2-T_0^2)\phi^2-E\,T\,\phi^3+\frac{\lambda}{4}\phi^4.$$9 in units $T_1>T>T_c$0. The factor $T_1>T>T_c$1 is the fraction of volume still in the symmetric phase. The source term $T_1>T>T_c$2 creates broken-phase fluctuations of radius $T_1>T>T_c$3 inside the symmetric background, while $T_1>T>T_c$4 and $T_1>T>T_c$5 erase them through reverse fluctuations and thermal noise, respectively (Chen et al., 24 May 2026).

Using a Gaussian profile

$T_1>T>T_c$6

the free energy of a broken-phase subcritical fluctuation is

$T_1>T>T_c$7

with

$T_1>T>T_c$8

and

$T_1>T>T_c$9

The corresponding source rate is

$\phi_s=0$0

where

$\phi_s=0$1

$\phi_s=0$2

and

$\phi_s=0$3

with $\phi_s=0$4 otherwise.

For the reverse process, the profile

$\phi_s=0$5

gives

$\phi_s=0$6

and

$\phi_s=0$7

with

$\phi_s=0$8

and

$\phi_s=0$9

The thermal-noise erasure rate is

$\phi_b>0$0

This structure makes the kinetics explicitly a competition between fluctuation production, collapse, reverse conversion, and Hubble dilution (Chen et al., 24 May 2026).

3. Subcritical volume fraction and the single-bin approximation

The geometric broken-phase volume fraction associated with the radius distribution is

$\phi_b>0$1

To incorporate random overlap when $\phi_b>0$2 is not asymptotically small, the Poisson expression

$\phi_b>0$3

is used for the actual broken-phase volume fraction, so that

$\phi_b>0$4

is the symmetric-phase volume fraction (Chen et al., 24 May 2026).

A simplified treatment replaces the full distribution by a single representative bin at $\phi_b>0$5. Writing

$\phi_b>0$6

and neglecting boundary flux, one obtains

$\phi_b>0$7

with

$\phi_b>0$8

and

$\phi_b>0$9

In quasi-steady state, $T_1$0, so

$T_1$1

The resulting fast estimate for the subcritical volume fraction is

$T_1$2

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 $T_1$3 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

$T_1$4

with nucleation rate

$T_1$5

where $T_1$6 (Chen et al., 24 May 2026).

The expected number of critical bubbles per Hubble volume is

$T_1$7

and the nucleation temperature $T_1$8 is defined by

$T_1$9

The fraction of space converted by critical bubbles is

$T_c$0

where

$T_c$1

The conceptual comparison is then straightforward. Standard analyses assume that before $T_c$2 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_c$3 down to $T_c$4 and evaluating $T_c$5. The homogeneous-background approximation is classified as valid when

$T_c$6

with small subcritical corrections for

$T_c$7

and breakdown of the homogeneous-background assumption when

$T_c$8

In the latter case, the system is categorized as a mixed background rather than an ordinary homogeneous bounce point (Chen et al., 24 May 2026).

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 $T_c$9. At the nucleation temperature,

$V(0,T_c)=V(\phi_c,T_c)$0

Across a scan of weak transitions in which $V(0,T_c)=V(\phi_c,T_c)$1 are varied and $V(0,T_c)=V(\phi_c,T_c)$2 is determined, the comparison between the fast estimate and the full kinetic result for points with $V(0,T_c)=V(\phi_c,T_c)$3 yields the fit

$V(0,T_c)=V(\phi_c,T_c)$4

Imposing $V(0,T_c)=V(\phi_c,T_c)$5 gives

$V(0,T_c)=V(\phi_c,T_c)$6

and therefore

$V(0,T_c)=V(\phi_c,T_c)$7

This leads to the practical rule

$V(0,T_c)=V(\phi_c,T_c)$8

Points above this boundary are flagged as mixed-background candidates rather than ordinary homogeneous bounce points (Chen et al., 24 May 2026).

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 $V(0,T_c)=V(\phi_c,T_c)$9, 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 $R$0. 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 (Chen et al., 24 May 2026). More specifically, large subcritical fractions are found when:

• Small free-energy splitting: $R$1 is small relative to $R$2, so broken patches have a low volume-energy cost.
• Low barrier height: $R$3 is low, reducing Boltzmann suppression.
• Small order parameter: $R$4 is small, which reduces the gradient term because $R$5.
• Weak transition strength: $R$6, so cooling is slow and the system remains near $R$7 for longer.

Conversely, stronger transitions with large $R$8, large $R$9, 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 $R_c(T)$0, or solve the full $R_c(T)$1 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 (Chen et al., 24 May 2026).