Gravitational Waves and Cosmological Observables from First-Order Phase Transitions: Thermal Corrections at Low Temperature

1. Problem Statement

Cosmological first-order phase transitions (FOPTs) in the early Universe are potential sources of stochastic gravitational waves (GWs) and can alter cosmological observables (e.g., via dark radiation). The characteristics of these signals depend critically on the thermal effective potential, $V(\phi, T)$.

While high-temperature thermal corrections (arising from light degrees of freedom where $m(\phi) \ll T$) are well-understood and typically modeled via high-$T$ expansions, there is a gap in understanding low-temperature thermal corrections. These arise from heavy degrees of freedom where the field-dependent mass $m(\phi)$ satisfies:

• $m(\phi) \ll T_N$ (nucleation temperature) in the false vacuum (symmetry-preserving phase, $\phi \approx 0$).
• $m(\phi) \gg T_N$ in the true vacuum (symmetry-breaking phase, $\phi = \phi_+$).

In the true vacuum, these heavy fields are Boltzmann-suppressed, but in the false vacuum, they contribute significantly to the potential. Calculating the exact form of these corrections requires knowing the specific heavy degrees of freedom integrated out, which is often unknown in Beyond the Standard Model (BSM) scenarios. The authors aim to determine how these unknown low-$T$ corrections generally impact FOPT parameters and the resulting GW signal without needing a specific microscopic model.

2. Methodology

The authors employ a combination of analytic approximation and numerical validation:

• Parameterization Strategy:
  Instead of calculating specific loop corrections for unknown heavy fields, they model the net effect of low-$T$ corrections as a shift in the effective potential at the false vacuum ($\phi=0$). They introduce a single dimensionless parameter, $f$, such that the shift is:
  $$\Delta V_{\text{low-}T} \approx f \left(\frac{\Lambda}{v}\right)^4 T^4$$
  where $\Lambda$ is the energy scale and $v$ is the vacuum expectation value. They argue that the precise behavior of this correction between the false and true vacua is less important than its magnitude at the false vacuum, provided the thin-wall approximation holds.

• Analytic Analysis:
  Using the thin-wall approximation for the bubble nucleation action $S_3(T)$, they derive analytic expressions for the shifts in key thermal parameters caused by a non-zero $f$:

  • Nucleation Temperature ($T_N$): Determined by solving $S_3(T_N)/T_N = C(T_N)$.
  • Transition Rate ($\beta/H$): Related to the time derivative of the action.
  • Latent Heat ($\xi$): Related to the difference in energy density and entropy between vacua.
    They derive linear relationships showing how $\delta T_N$, $\delta(\beta/H)$, and $\delta\xi$ scale with $f$.

• Numerical Validation:
  To test the analytic approximations, they construct a specific quartic scalar potential model. They implement the low-$T$ correction using a "window function" $w(\tilde{\phi})$ that smoothly turns off the correction as $\phi$ approaches the true vacuum (mimicking the Boltzmann suppression). They use the CosmoTransitions package to numerically calculate the bounce solution, $S_3(T)$, and the resulting thermal parameters for various values of $f$ and compare them to the analytic predictions.

• GW Signal Calculation:
  They map the shifts in thermal parameters to the peak frequency ($f_{sw}$) and peak amplitude ($h^2\Omega_{sw}$) of the gravitational wave signal generated by sound waves in the plasma.

3. Key Contributions

1. Universal Parameterization: The paper establishes that the complex effects of heavy degrees of freedom (which are light in the false vacuum but heavy in the true vacuum) can be effectively captured by a single parameter $f$. This allows for a "dictionary" between unknown UV physics and observable GW signatures.
2. Analytic Shift Formulas: The authors derive explicit analytic formulas for the fractional shifts in $T_N$, $\beta/H$, and $\xi$ as functions of $f$.
   • $T_N$ increases linearly with $f$.
   • $\beta/H$ decreases linearly with $f$.
   • $\xi$ decreases linearly with $f$.
3. Insensitivity to Barrier Details: They demonstrate that the shifts in thermal parameters depend primarily on the magnitude of the potential shift at the false vacuum, not on the detailed shape of the correction within the barrier between vacua.
4. Cosmological Implications: They highlight the impact on the effective number of relativistic species ($N_{\text{eff}}$), showing that even a single heavy degree of freedom can significantly alter $N_{\text{eff}}$ in dark sector transitions, potentially violating observational constraints if not accounted for.

4. Key Results

• Thermal Parameter Shifts:
  • Nucleation Temperature ($T_N$): Increases with $f$. The shift is generally small ($\delta T_N / T_N \ll 1$) but positive.
  • Transition Rate ($\beta/H$): Decreases significantly with $f$. Since $\beta/H$ sets the duration of the transition, a lower value implies a longer-lasting transition.
  • Latent Heat ($\xi$): Decreases with $f$. Counter-intuitively, while the potential energy difference ($\Delta V$) increases, the latent heat (free energy difference) decreases because the entropy term dominates the shift.
• Gravitational Wave Signal:
  • Frequency ($f_{sw}$): Since $f_{sw} \propto \beta/H$, and $\beta/H$ decreases as $f$ increases, the GW frequency shifts to lower values.
  • Amplitude ($h^2\Omega_{sw}$): The amplitude scales roughly as $\xi^2 (\beta/H)^{-1}$. Since $\xi$ decreases and $\beta/H$ decreases, the net effect depends on the competition. However, the authors find that the dependence on $\xi$ usually dominates, leading to a decrease in GW amplitude as $f$ increases.
• Numerical Agreement: The numerical results from CosmoTransitions for the quartic potential model align remarkably well with the analytic predictions derived without assuming the specific window function shape. This validates the assumption that the barrier details are sub-dominant.
• Magnitude of Effect: For a single degree of freedom, $f \sim \mathcal{O}(0.1)$. This can lead to shifts in $\xi$ of up to 10% and significant changes in the GW signal, which are non-negligible for precision cosmology.

5. Significance

• Robustness of GW Predictions: The study provides a robust framework for interpreting GW signals from FOPTs. It suggests that uncertainties regarding heavy, unknown particles in the dark sector can be systematically parameterized, rather than ignored.
• Model Independence: By showing that the results hold for a general class of potentials (under the thin-wall approximation), the paper offers a tool applicable to a wide range of BSM theories without requiring a full recalculation of the thermal potential for every specific model.
• Observational Constraints: The findings imply that future GW observatories (like LISA, DECIGO, or PTA) could constrain the spectrum of heavy particles in the early Universe. Conversely, if a GW signal is detected, the low-$T$ corrections must be accounted for to correctly infer the underlying particle physics parameters.
• Dark Radiation Constraints: The paper underscores that low-$T$ corrections can drastically alter the latent heat released, thereby changing the constraints on dark sector transitions from $N_{\text{eff}}$ measurements. This is crucial for models where the dark sector temperature is comparable to the Standard Model temperature.

In summary, the paper bridges the gap between unknown high-energy physics (heavy degrees of freedom) and low-energy cosmological observables, providing a simplified yet accurate parameterization to predict how these "invisible" corrections modify the gravitational wave signatures of the early Universe.