⚛️ phenomenology

Big Bang Nucleosynthesis as a Probe of First-Order Phase Transitions

This paper demonstrates that precise measurements of deuterium abundances from Big Bang Nucleosynthesis constrain the parameter space of supercooled first-order phase transitions, specifically disfavoring scenarios compatible with the NANOGrav gravitational wave signal unless the baryon asymmetry is generated after or during the transition rather than before.

Original authors: Masanori Tanaka

Published 2026-09-24

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Masanori Tanaka

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the earliest moments of the universe, a fraction of a second after the Big Bang, the cosmos was a seething, ultra-hot soup of fundamental particles. As this primordial fireball expanded and cooled, it underwent a dramatic transformation known as a phase transition, much like water freezing into ice. In the standard story of our universe, this change happened smoothly and uniformly everywhere at once. However, some theories suggest that under certain conditions, this transition could have been violent and uneven, occurring in a "first-order" fashion where bubbles of the new, stable state of matter nucleated and expanded within the old, supercooled state. If such a chaotic transition occurred, it would have left a distinct fingerprint on the universe's thermal history, potentially altering the way light and matter interacted in the first few minutes of existence.

This early era is crucial because it set the stage for Big Bang Nucleosynthesis, the process that forged the first atomic nuclei. During this brief window, protons and neutrons fused to create the lightest elements: hydrogen, helium, and a tiny but vital amount of deuterium. The amount of deuterium produced is exquisitely sensitive to the density of matter in the universe at that time. Today, astronomers can measure the abundance of deuterium in ancient gas clouds with incredible precision, providing a strict test for any theory about the universe's infancy. If the conditions during nucleosynthesis were not uniform, the resulting chemical recipe would look different from what we observe.

A recent study by Masanori Tanaka explores how a specific type of violent phase transition could have disrupted this chemical balance. The research focuses on a scenario where the universe supercooled, meaning it dropped below the temperature required for the phase transition to occur before the change actually happened. When the transition finally triggered, it did so stochastically, with bubbles of the new phase appearing at random times and locations. This randomness meant that some regions of space reheated—gained heat from the energy released by the transition—earlier than others. Because the universe was expanding, these regions of different temperatures would have different densities of photons, the particles of light.

If the imbalance between matter and antimatter, known as the baryon asymmetry, had already been established before this chaotic transition, the uneven heating would have distorted the ratio of matter to light in different regions. In simpler terms, the "recipe" for the universe would have varied from place to place. The researchers calculated how these variations would evolve as the universe continued to expand and cool. They found that while the random motion of protons and neutrons tends to smooth out these differences over time, the smoothing is not always complete by the time nucleosynthesis begins.

The study reveals that if such a transition occurred at temperatures between 0.01 and 100 billion electron volts, the resulting unevenness in the matter-to-light ratio would likely have produced a deuterium abundance that contradicts current astronomical observations. Specifically, the data disfavors scenarios where the transition released a large amount of hidden energy and took a long time to complete. The researchers quantified this by simulating the power spectrum of the transition times and the subsequent diffusion of particles, showing that the constraints are tightest for transitions that are both energetic and slow.

This finding has a significant implication for a separate line of inquiry involving gravitational waves. Recent observations by the NANOGrav collaboration have detected a background hum of gravitational waves, which some theorists have tried to explain as the echo of these very same violent phase transitions. Tanaka's work suggests that if the baryon asymmetry existed before the transition, the specific parameters required to match the gravitational wave signal are likely ruled out by the precise measurements of deuterium. The universe simply would not have the right chemical composition today if those two events happened in that specific order.

However, the study does not close the door on these theories entirely. It offers a way out: if the imbalance between matter and antimatter was generated after the transition, or if it was created in a way that perfectly tracked the local release of energy during the transition, the dangerous unevenness would not have formed. In this alternative scenario, the "recipe" remains uniform, and the deuterium measurements remain consistent with the gravitational wave signals. Thus, the timing of how matter came to dominate over antimatter in the early universe may be constrained by combining the new gravitational wave data with the ancient chemical fingerprints left in the light elements. The study concludes that while the simplest version of these violent transitions is incompatible with our current understanding of the early universe's chemistry, more complex timing for the origin of matter could still be viable.

Technical Summary: Big Bang Nucleosynthesis as a Probe of First-Order Phase Transitions

Problem Statement
Standard cosmology relies on Big Bang Nucleosynthesis (BBN) to explain primordial light element abundances, providing a critical probe of the early universe's thermal history. Recent work has highlighted that precision measurements of these abundances can constrain baryogenesis scenarios that generate large spatial inhomogeneities in the baryon-to-photon ratio ( $\eta$ ). While diffusion typically smears such inhomogeneities, Ref. [5] demonstrated that inhomogeneities generated at temperatures below a few TeV can survive until the BBN epoch if they span the horizon scale.

This paper addresses a specific mechanism: supercooled cosmological first-order phase transitions (FOPTs). In these scenarios, the stochastic timing of bubble nucleation creates significant temperature inhomogeneities via the reheating process. If a baryon asymmetry exists prior to the FOPT and is conserved during the transition, these temperature fluctuations are imprinted onto the baryon-to-entropy ratio ( $Y_B$ ) or $\eta$ . Previous analyses, such as Ref. [6], constrained FOPT parameters based on temperature fluctuations but neglected the specific dependence of reheating on the induced $Y_B$ fluctuations. Furthermore, the interpretation of the NANOGrav gravitational wave (GW) signal via FOPTs requires understanding whether such models are compatible with BBN constraints on baryon inhomogeneity.

Methodology
The author extends the formalism of Ref. [6] by incorporating the reheating dependence of induced $Y_B$ fluctuations and combining the transition-time power spectrum from finite-bubble statistics with proton and neutron diffusion dynamics at the BBN epoch.

FOPT Dynamics: The study assumes a stochastic bubble nucleation rate $\Gamma(t) = \Gamma_0 e^{\beta(t-t_n)}$ , where $\beta$ is the inverse duration of the transition and $t_n$ is the onset time. The fraction of the false vacuum is calculated using finite-bubble statistics.
Reheating and Temperature Inhomogeneity: The author models the radiation energy density evolution in regions where the transition completes at different times ( $t_c$ ). They define a reference radiation energy for the average transition time ( $\bar{t}$ ) and compare it to regions undergoing early or late transitions. The vacuum energy released ( $\rho_V$ ) is converted into radiation, creating a temperature contrast ( $T_{dev}$ vs. $T_{ref}$ ).
Baryon Asymmetry Fluctuations: Assuming the baryon asymmetry is generated before the FOPT and the comoving baryon number is conserved, the author derives the fluctuation in the baryon-to-entropy ratio ( $Y_B$ ). They establish that $Y_B$ fluctuations are driven by the ratio of local entropy densities, which depends on the local reheating temperature.
Power Spectrum and Diffusion: To account for spatial correlations, the author utilizes the transition-time power spectrum $P_{\delta t}(k)$ derived from finite-bubble statistics (Ref. [29]). They construct an effective parametrization for the baryon-asymmetry power spectrum, $\Delta^2(k)$ , which scales with the total fluctuation amplitude $\epsilon_0^2$ (dependent on latent heat $\alpha_n$ , $\beta/H_n$ , and $T_n$ ) and the spectral shape of the transition time.
BBN Constraints: The root-mean-square (RMS) fluctuation at the BBN epoch ( $\epsilon_{BBN}$ ) is calculated by integrating the power spectrum over wavenumbers, weighted by a diffusion damping factor $|T_p(k)|^2$ . The author considers two cutoff scales: an optimistic case ( $k_{cut} = \infty$ ) and a conservative case ( $k_{cut} = k_{th}$ , the inverse mean bubble separation).

Key Contributions

Reheating Dependence: The paper quantifies the impact of reheating on temperature and baryon-asymmetry fluctuations, showing that the fluctuation amplitude $\epsilon_0$ depends explicitly on the latent heat parameter $\alpha_n$ , a dependency missing in previous analyses (Ref. [6]).
Spectral Shape Integration: By incorporating the transition-time power spectrum from finite-bubble statistics, the author provides a more rigorous treatment of the scale dependence of baryon inhomogeneities, moving beyond simple variance estimates.
Updated Formalism: The author clarifies the relationship between the baryon-to-photon ratio ( $\eta$ ) and the baryon-to-entropy ratio ( $Y_B$ ), noting that while they differ immediately after the FOPT due to $g_*$ variations, they converge by the BBN epoch, allowing $Y_B$ to serve as the fiducial matching variable.

Results

BBN Constraints on FOPT Parameters: For nucleation temperatures in the range $10^{-2} \text{ GeV} \lesssim T_n \lesssim 10^2 \text{ GeV}$ , precise measurements of deuterium abundances disfavor parameter regions characterized by a sizable latent heat ( $\alpha_n$ ) and a long transition duration (small $\beta/H_n$ ) at the $1\sigma$ level.
Dependence on $T_n$ : The bound weakens at higher nucleation temperatures as diffusion becomes more efficient.
No Constraint for Small $\alpha_n$ : No significant constraint is obtained for $\alpha_n < 0.1$ .
NANOGrav Compatibility: The study finds that the FOPT parameter regions compatible with the NANOGrav gravitational wave signal are disfavored by light element abundance measurements, provided the power spectrum from Ref. [29] remains valid down to sub-bubble scales.

Significance and Claims
The paper concludes that if the NANOGrav signal is explained by FOPTs, the baryon asymmetry cannot simply predate the transition in the standard manner assumed in the analysis. To evade the BBN bounds, the author proposes two alternative scenarios:

Post-Reheating Baryogenesis: The baryon asymmetry is generated after the reheating process is complete.
Transition-Period Baryogenesis: The asymmetry is generated during the FOPT, provided the local baryon production tracks the local entropy injection, thereby suppressing the resulting $Y_B$ fluctuation.

Thus, the paper claims that combining gravitational wave observations with light element abundance measurements offers a method to constrain the timing of baryogenesis. The results highlight the critical importance of the reheating effect when deriving BBN-based bounds on first-order phase transitions.

More like this