Phenomenology of bubble size distributions in a first-order phase transition

Basic Information

• Paper ID: 2407.15419
• Title: Phenomenology of bubble size distributions in a first-order phase transition
• Authors: Danny Marfatia, Po-Yan Tseng, Yu-Min Yeh
• Classification: hep-ph astro-ph.CO hep-ex
• Publication Date: July 2024 (arXiv:2407.15419v2, revised September 23, 2025)
• Paper Link: https://arxiv.org/abs/2407.15419

Abstract

This paper investigates the phenomenological effects of bubble size distributions in cosmological first-order phase transitions (FOPT). In contrast to the commonly assumed monochromatic distribution, the radius distributions of vacuum and false vacuum bubbles should be extended. This leads to extended mass distributions for Fermi balls (FBs) and primordial black holes (PBHs) produced by FOPT in the dark matter sector. The study demonstrates that extended distributions shift the gravitational wave spectrum peak toward lower frequencies and broaden it below the peak frequency; microlensing observations can probe PBH distributions with mean masses below $10^{-10}M_\odot$; the extragalactic gamma-ray spectrum from PBH Hawking evaporation may exhibit power-law spectral breaks in the 5-10 MeV range.

Research Background and Motivation

Core Problem

Traditional studies typically assume that vacuum bubbles in FOPT have fixed sizes, resulting in monochromatic mass distributions for the produced FBs and PBHs. However, this assumption neglects the inherent distribution characteristics of bubble sizes.

Significance

1. Theoretical Completeness: Bubble size distribution is an intrinsic feature of FOPT; the monochromatic assumption is overly simplified
2. Observational Predictions: Extended mass distributions significantly alter the characteristics of various observational signals
3. Parameter Space: Extended distributions may expand the observable FOPT parameter range

Existing Limitations

• Previous phenomenological studies were primarily based on monochromatic mass distribution assumptions
• Lack of systematic investigation of how bubble size distributions affect multi-messenger observational signals
• Gravitational wave spectrum calculations insufficiently account for bubble size distribution uncertainties

Core Contributions

1. Established a theoretical framework connecting false vacuum bubble radius distributions to FB/PBH mass distributions
2. Systematically studied the effects of extended distributions on gravitational wave spectra, revealing spectral peak frequency shifts and asymmetric broadening
3. Analyzed the impact of extended mass distributions on microlensing observations, expanding the detectable parameter space
4. Predicted characteristic spectral breaks in extragalactic gamma-ray spectra from PBH Hawking evaporation
5. Provided seven benchmark points for concrete phenomenological analysis

Methodology Details

Task Definition

Investigate how bubble size distributions in FOPT affect the mass distributions of FBs/PBHs, and analyze the impact of these extended distributions on gravitational wave, microlensing, and gamma-ray observational signals.

Theoretical Framework

1. Bubble Radius Distribution

At the percolation temperature $T_\star$, the radius distributions of vacuum (TV) and false vacuum (FV) bubbles are respectively:

False Vacuum Bubble Distribution: $\frac{dn_{FV}}{dR_f}(t_\star) \simeq \frac{\beta^4}{192v_w^4} x^4 e^{-x}(1-e^{-x})$

Vacuum Bubble Distribution: $\frac{dn_{TV}}{dR}(t_\star) \simeq \frac{\beta^4}{8\pi v_w^4} \frac{1}{x} e^{-1/x}$

where $x \equiv \exp(\beta R_{(f)}/v_w) > 1$, $\beta^{-1}$ is the phase transition duration, and $v_w$ is the bubble wall velocity.

2. Mass Distribution Mapping

The relationship between FB/PBH mass and false vacuum bubble radius is $M_{FB/PBH} \propto R_f^3$, yielding:

$\frac{dn_{FB/PBH}}{dM_{FB/PBH}} = \frac{R_f}{3M_{FB/PBH}} \frac{dn_{FV}}{dR_f}(t_\star)$

Normalized mass distribution function: $f(M_{FB/PBH}) = \frac{\frac{4\pi R_f^3}{3} \frac{dn_{FB/PBH}}{dM_{FB/PBH}}}{\int dM_{FB/PBH} \frac{4\pi R_f^3}{3} \frac{dn_{FB/PBH}}{dM_{FB/PBH}}}$

Phenomenological Analysis Methods

1. Gravitational Wave Spectrum Calculation

Employing a double-broken power-law (DBPL) template, convolving with extended bubble distributions:

$\Omega_{GW}(f) = \frac{\int dR \frac{dn_{TV}(t_\star)}{dR} \Omega_{GW}^{DBPL}(f,R)}{\int dR \frac{dn_{TV}(t_\star)}{dR}}$

2. Microlensing Event Rate

Number of microlensing events for extended mass distributions:

$N_{events} = N_S T_{obs} \int dM_{FB/PBH} \frac{dn_{FB/PBH}}{dM_{FB/PBH}} \int_{t_{min}}^{t_{max}} dt_E \int dR_S \int dx \frac{d^2\Gamma(M_{FB/PBH})}{dx dt_E} \frac{dn}{dR_S}$

3. Extragalactic Gamma-Ray Spectrum

Extragalactic photon flux considering PBH mass evolution:

$\frac{d\Phi}{dE} = \int_{t_{CMB}}^{\min(t_{eva},t_0)} dt \, c[1+z(t)]n_{PBH}|_0 \int dM_{PBH} \frac{1}{n_{PBH}|_\star} \frac{dn_{PBH}|_\star}{dM_{PBH}} \frac{d^2N_\gamma(M_{PBH})}{d\tilde{E}dt}\bigg|_{\tilde{E}=[1+z(t)]E}$

Experimental Setup

Benchmark Point Parameters

The paper establishes seven benchmark points (BP-1 through BP-7) covering different FOPT parameter combinations:

• Energy scale: $B^{1/4} \sim 0.04-20$ MeV
• Phase transition temperature: $T_\star \sim 0.06-6$ MeV
• Dark matter asymmetry parameter: $\eta_\chi \sim 10^{-18} - 10^{-9}$

BP-1 through BP-5 produce PBHs, while BP-6 and BP-7 produce stable FBs.

Observational Constraints

1. Gravitational Waves: Detector sensitivities for LISA, NANOGrav, SKA, etc.
2. Microlensing: Subaru-HSC observational data (8.7×10^7 stars, 7-hour observation)
3. Gamma Rays: COMPTEL and Fermi-LAT limits, expected sensitivities of AMEGO-X/e-ASTROGAM

Evaluation Metrics

• 95% C.L. exclusion regions and sensitivity predictions
• 3σ signal detection capability
• Consistency with existing observational data

Experimental Results

Main Results

1. Gravitational Wave Spectrum Characteristics

• Frequency Shift: Extended distributions shift the GW spectrum peak frequency toward lower frequencies
• Spectral Broadening: Significant asymmetric broadening appears below the peak frequency
• Amplitude Changes: Extended distributions show slightly higher amplitude than monochromatic distributions

Mean bubble radius relationship: $\langle R_f \rangle_\star = 1.32v_w/\beta > R_{f\star} = 6^{1/4}v_w/\beta$

2. Microlensing Observational Capability

Fermi Ball Microlensing:

• Sensitive mass range: Most sensitive around $\sim 10^{-9} M_\odot$
• Energy scale: $1$ keV $\lesssim B^{1/4} \lesssim 10$ MeV
• Extended distributions show better sensitivity to larger mean masses compared to monochromatic distributions

PBH Microlensing:

• Monochromatic distribution limits: $10^{-10} \lesssim M_{PBH}/M_\odot \lesssim 10^{-4}$
• Extended distributions can probe PBH distributions with mean masses $< 10^{-10}M_\odot$

3. Extragalactic Gamma-Ray Signals

• Spectral Break Features: BP-1 and BP-2 show moderate spectral breaks in the 5-10 MeV range
• Detection Capability: All benchmark points produce 3σ signals with AMEGO-X/e-ASTROGAM
• Mass Correlation: $0.1$ MeV $\lesssim B^{1/4} \lesssim 10$ MeV corresponds to $10^{13}$ g $\lesssim \langle M_{PBH} \rangle \lesssim 10^{20}$ g

Extended Distribution vs. Monochromatic Distribution Comparison

Related Work

Main Research Directions

1. FOPT Theory: Bubble nucleation, percolation theory, phase transition dynamics
2. PBH Formation Mechanisms: Density fluctuation collapse, phase transition-induced formation
3. Multi-Messenger Astrophysics: Gravitational wave, microlensing, and gamma-ray observations

Innovations in This Work

• First systematic investigation of phenomenological effects of bubble size distributions
• Established complete theoretical chain from bubble distributions to observational signals
• Provided observational criteria to distinguish extended distributions from monochromatic distributions

Conclusions and Discussion

Main Conclusions

1. Bubble size distribution is an indispensable factor in FOPT phenomenology, significantly affecting various observational signals
2. Gravitational wave spectrum calculations must account for bubble distribution uncertainties, introducing additional systematic errors
3. Extended mass distributions expand the observable FOPT parameter space, particularly in microlensing observations
4. Spectral break features in extragalactic gamma-ray spectra can serve as unique signals for identifying extended distributions

Limitations

1. Gravitational Wave Spectrum Template Limitations: DBPL template is based on numerical simulations of simultaneous nucleation, differing from exponential nucleation
2. Model Dependence: FB cooling rates and collapse criteria exhibit strong model dependence
3. Approximation Treatments: Tetrahedral approximation, geometric optics approximation, etc., introduce systematic errors

Future Directions

1. Improved Numerical Simulations: Develop large-scale hydrodynamic simulations accounting for exponential nucleation
2. Observational Verification: Utilize future gravitational wave, microlensing, and gamma-ray observations to verify theoretical predictions
3. Parameter Estimation: Develop methods to infer FOPT parameters and bubble distributions from observational data

In-Depth Evaluation

Strengths

1. Theoretical Completeness: Systematically addresses the complete physical chain from bubble distributions to observational signals
2. Phenomenological Richness: Encompasses multiple observational methods including gravitational waves, microlensing, and gamma rays
3. Quantitative Analysis: Provides specific benchmark points and numerical predictions for easy comparison with observations
4. Methodological Innovation: First to incorporate bubble size distribution effects into FOPT phenomenological studies

Weaknesses

1. Template Limitations: Gravitational wave spectrum calculations depend on DBPL template with restricted applicability range
2. Parameter Selection: Benchmark point selection may be incomplete, insufficiently covering parameter space
3. Uncertainty Quantification: Quantitative analysis of theoretical uncertainties is insufficient

Impact

• Academic Value: Provides new theoretical framework for FOPT phenomenological research
• Practical Value: Offers theoretical guidance for future observational experiment data analysis
• Field Advancement: May promote interdisciplinary research between phase transition cosmology and multi-messenger astrophysics

Applicable Scenarios

1. Theoretical Research: Dark matter model construction, early universe phase transition studies
2. Observational Analysis: Gravitational wave detector data analysis, microlensing surveys, gamma-ray astronomy
3. Experiment Design: Scientific objective setting and observational strategy optimization for future detectors

References

The paper cites 15 core references, encompassing:

• FOPT theoretical foundations 10,11,12
• Numerical simulation studies 13,14,16
• Observational constraints 17,25,26
• Phenomenological analysis 1,2,3,6,22

These references form the solid theoretical foundation and observational basis for this research.