Primordial Black Hole Formation in Dust-Radiation Bouncing Cosmologies

This paper establishes a unified framework for studying primordial black hole formation in dust-radiation bouncing cosmologies and demonstrates that while the curvature perturbation threshold is extremely low, the resulting abundance of primordial black holes is vanishingly small due to the significant suppression caused by radiation pressure and two-fluid collapse conditions.

The Gist

The Big Idea: Can a "Bouncing" Universe Make Black Holes?

Imagine the history of our universe not as a balloon that started with a tiny bang and has been inflating ever since, but as a giant rubber ball. In this "Bouncing Cosmology" theory, the universe was once shrinking (contracting), hit a tiny, dense point (the bounce), and then started expanding again.

The authors of this paper asked a specific question: During that shrinking phase, before the bounce, could the universe have squeezed itself so hard that it created tiny, ancient black holes called Primordial Black Holes (PBHs)?

These black holes are interesting because they could be the "Dark Matter" that holds galaxies together, or they could be the seeds for the supermassive black holes we see in the centers of galaxies today.

The Cast of Characters

To understand the paper, we need to meet the two main ingredients of the early universe:

1. Dust: Think of this as cold, slow-moving matter (like stars and gas clouds). It doesn't push back; it just falls in.
2. Radiation: Think of this as hot light and energy. It moves fast and pushes back hard (like steam in a pressure cooker).

In previous studies, scientists looked at a universe made only of Dust. They found that if the universe shrinks, the dust clumps together easily and forms black holes. But our real universe has both Dust and Radiation. This paper asks: What happens when you mix them?

The Story: A Tug-of-War

The authors built a mathematical model to simulate what happens when the universe shrinks with both Dust and Radiation inside. They used three main tools to tell the story:

1. The "Jeans Length" (The Size Limit)

Imagine you are trying to build a sandcastle on a beach. If the sand is too loose (like radiation), the waves wash it away before you can build a tower. You need a certain amount of sand packed tightly enough to hold its shape.
In physics, this is the Jeans Length. It's the minimum size a clump of matter needs to be to overcome the "pushing" force of radiation and collapse into a black hole.

• The Finding: The authors calculated this limit for a mix of dust and radiation. They found that radiation acts like a strong wind, making it much harder for the dust to clump together. You need a much bigger "sandcastle" to survive the wind.

2. The "Three-Zone Model" (The Collapse Test)

Once a clump gets big enough, does it actually collapse? The authors used a "Three-Zone Model" to test this. Imagine a bubble of dense matter in the middle of the universe.

• The Test: They asked, "How fast does the 'pressure' (sound waves) travel through this bubble compared to how fast the bubble is collapsing?"
• The Analogy: If you have a balloon inflating, and you poke a hole in it, the air rushes out. If the air rushes out faster than the balloon can shrink, the balloon survives. But if the balloon shrinks faster than the air can escape, it pops (collapses into a black hole).
• The Twist: They tested two different rules for when the "pop" happens. One rule is lenient; the other is strict. Both rules agreed on the outcome.

3. The "Power Spectrum" (The Noise Level)

For a black hole to form, the universe needs to be "noisy" enough. Imagine the universe is a calm lake. To make a splash (a black hole), you need a big wave.
The authors calculated how big the waves (fluctuations) were in this shrinking universe.

• The Result: The waves were incredibly small. The "noise" was almost non-existent.

The Verdict: A Resounding "No"

After running the numbers for black holes ranging from the size of a mountain to the size of a galaxy, the authors found a disappointing (but scientifically important) result:

Primordial Black Holes are almost impossible to form in this specific type of bouncing universe.

Here is why, using a simple metaphor:

• The Scenario: Imagine trying to build a sandcastle (a black hole) on a beach during a hurricane (the radiation pressure).
• The Problem: Even though the universe was shrinking (which usually helps build things), the "wind" from the radiation was so strong that it blew away the sand before it could pile up.
• The Numbers: The math showed that the chance of a black hole forming was so tiny it was effectively zero. It's like trying to win the lottery every single second for a billion years and never winning.

Why Does This Matter?

You might think, "So what? We just found out black holes don't form this way."

Actually, this is a huge discovery for two reasons:

1. It rules out a theory: If we want to explain Dark Matter using these ancient black holes, this specific "Dust-Radiation Bouncing" model doesn't work. The universe is too "smooth" and the radiation is too strong.
2. It sets a new bar: The authors showed that for these black holes to exist, the universe would need a "supercharger"—some other unknown mechanism to make the waves (fluctuations) much bigger than they naturally are.

Summary in One Sentence

The authors simulated a universe shrinking with both matter and light, and found that the light pushes back so hard that it prevents the matter from clumping together, meaning no significant number of ancient black holes can form this way.

Technical Summary

1. Problem Statement

Primordial Black Holes (PBHs) are a compelling candidate for Dark Matter and a probe of early-universe physics. While previous studies (e.g., Ref. [46]) investigated PBH formation in pure-dust bouncing cosmologies, those models predicted negligible PBH abundance for masses surviving to the present day, as the mass spectrum was dominated by low-mass PBHs that evaporate via Hawking radiation.

However, realistic early-universe scenarios involve a mixture of dust and radiation. The presence of radiation introduces significant physical complexities:

• Radiation pressure resists gravitational collapse, potentially suppressing PBH formation.
• The coupling between dust and radiation perturbations prevents the use of exact analytic solutions (like Lemaître–Tolman–Bondi) used in pure-dust models.
• It remains unclear whether the enhanced perturbation growth in bouncing cosmologies can overcome radiation pressure to produce a viable PBH abundance across a broad mass spectrum (from $10^{-17} M_\odot$ to $10^{13} M_\odot$).

2. Methodology

The authors developed a unified framework to study PBH formation during the classical contracting phase of a dust-radiation bouncing cosmology. The results are designed to be insensitive to the specific quantum details of the bounce itself. The methodology consists of three main components:

A. Background and Perturbation Evolution

• Model: A quantum bouncing cosmology with two fluids: dust (equation of state $w \approx 9.44 \times 10^{-22}$) and radiation ($w=1/3$).
• Perturbation Spectrum: The curvature power spectrum ($P_\zeta$) was computed using a stable semi-analytical (adiabatic) method. This approach avoids numerical instabilities caused by the extremely small dust equation of state.
• Approximation: The entropy modes were neglected because the curvature mode dominates during contraction, allowing the system to be treated as a single effective curvature perturbation.

B. Jeans Stability Analysis

• Dynamical System Approach: The authors generalized the dynamical-system analysis (Refs. [55–57]) to a coupled two-fluid system to derive the Jeans wavenumber ($k_J$).
• Validity Regime: They established that gravitational instability (collapse) is only possible when the Jeans length is smaller than the Hubble radius. They identified a transition scale factor ($\bar{a}_{Tr}$) where $k_J$ exceeds the Hubble scale. PBH formation is only analyzed for scale factors $\bar{a} > \bar{a}_{Tr}$.

C. Collapse Modeling (Three-Zone Model)

• Generalization: The single-fluid three-zone model (Harada et al.) was extended to a two-fluid (dust+radiation) system. The overdense region is modeled as a Closed Robertson-Walker (cRW) patch embedded in a flat background.
• Collapse Criteria: Two criteria were implemented to determine the critical threshold for collapse, accounting for sound waves that can disperse the overdensity:
  1. Criterion 1: The sound propagation time from the center to the surface must exceed the time of maximum expansion of the cRW patch.
  2. Criterion 2: The sound propagation time must exceed the time of apparent horizon formation (a stricter condition).
• Matching: The parameters of the cRW patch ($\eta_s, A_{eq}$) were matched to the quantum bouncing background at the time of PBH formation (when the mode crosses the Jeans scale).

D. Mass Fraction Calculation

• The Press-Schechter formalism was used to calculate the mass fraction ($\beta$) of PBHs.
• $\beta$ depends on the ratio of the critical curvature threshold ($\zeta_c$) to the root-mean-square of the curvature power spectrum ($\sqrt{P_\zeta}$) at the Jeans scale.

3. Key Contributions

1. Unified Two-Fluid Framework: Successfully extended PBH formation analysis from pure-dust to realistic dust-radiation bouncing cosmologies, incorporating the coupled dynamics of relativistic and non-relativistic fluids.
2. Stable Adiabatic Calculation: Developed a robust semi-analytical method to compute the curvature power spectrum in the presence of extreme parameter values (tiny $w$), ensuring numerical stability.
3. Generalized Jeans Scale: Derived the Jeans wavenumber for a coupled dust-radiation system using dynamical systems analysis, identifying the specific regime where sub-Hubble collapse is physically possible.
4. Refined Collapse Criteria: Applied and generalized two distinct collapse criteria (expansion-time vs. horizon-time) to the two-fluid case, demonstrating that radiation pressure significantly alters the threshold conditions.

4. Key Results

• Critical Threshold ($\zeta_c$): The critical curvature perturbation required for collapse is extremely small, $\zeta_c \sim 10^{-21}$, and is nearly mass-independent for masses between $10^{-14} M_\odot$ and $10^{13} M_\odot$.
• Power Spectrum Suppression: Despite the small threshold, the actual amplitude of the curvature power spectrum ($P_\zeta$) at the relevant formation times is vanishingly small ($10^{-152}$ to $10^{-82}$).
• Negligible Abundance: The mass fraction $\beta$ is effectively zero across all benchmark masses ($10^{-17} M_\odot$ to $10^{13} M_\odot$).
  • Reasoning: The ratio $\zeta_c / \sqrt{P_\zeta}$ is enormous ($10^{20}$ to $10^{54}$), making the probability of a perturbation exceeding the threshold statistically zero.
• Comparison with Pure-Dust: In contrast to previous pure-dust models where low-mass PBHs could form abundantly, the inclusion of radiation pressure and the two-fluid collapse conditions significantly suppresses PBH production.
• Criterion Comparison: Both Criterion 1 and Criterion 2 yield consistent results. Criterion 2 (horizon-based) allows for slightly easier formation in contracting backgrounds compared to expanding ones, but the difference is negligible in the final abundance calculation.

5. Significance and Conclusion

• Suppression Mechanism: The study demonstrates that in dust-radiation bouncing cosmologies, radiation pressure and the two-fluid collapse dynamics act as a strong suppression mechanism against PBH formation. The curvature perturbations amplified during contraction are insufficient to overcome these barriers.
• Implications for Dark Matter: Within the standard assumptions of this model (adiabatic perturbations, Gaussian statistics, and no additional amplification mechanisms), dust-radiation bouncing cosmologies cannot produce a significant abundance of PBHs to account for Dark Matter.
• Future Directions: To generate observable PBHs in such scenarios, additional mechanisms are required to further amplify the curvature power spectrum (e.g., features in the potential, non-Gaussianity, or modified gravity effects) to bridge the gap between the tiny spectrum amplitude and the critical threshold.

In summary, while bouncing cosmologies offer a mechanism to enhance perturbations, the specific inclusion of radiation in a dust-radiation mix renders the standard gravitational collapse channel ineffective for producing Primordial Black Holes without further model enhancements.