Primordial Black Holes Evaporating before Big Bang Nucleosynthesis

This study establishes a transparent framework to constrain the initial mass fraction of primordial black holes that evaporate before primordial nucleosynthesis, demonstrates that observable effects require masses greater than $10^9$ g, with maximum sensitivity at $2\times10^9$ g, and provides public code to facilitate future updates.

The Gist

Imagine the early universe as a giant, boiling pot of soup. For a long time, scientists believed this soup was perfectly smooth and predictable. Yet this article poses a fascinating question: What if there were tiny, invisible "black hole spots" floating in this soup, evaporating just before the soup began to cook into the first stars and elements?

Here is a simple summary of what the authors, Quan-feng Wu and Xun-Jie Xu, discovered.

1. The Tiny Time-Traveling Black Holes

Normally, we think of black holes as massive monsters formed from dying stars. However, the authors discuss Primordial Black Holes (PBHs). Imagine these as microscopic spots of darkness that formed instantly at the very beginning of time, like bubbles bursting in a soda can before the soda is even poured.

These spots are unstable. Thanks to a rule discovered by Stephen Hawking, they slowly lose energy and shrink until they vanish completely. The smaller they are, the faster they disappear.

2. The "Big Bang Cooking" (BBN)

About one second after the Big Bang, the universe cooled just enough for a "cooking process" called Big Bang Nucleosynthesis (BBN) to begin. This is the moment when the universe's basic ingredients—protons and neutrons—started fusing to form the first light elements like helium and deuterium (heavy hydrogen).

Imagine this as the moment the baker pushes the dough into the oven. If you then change the temperature or add a strange ingredient, the bread comes out completely different.

3. The "Evaporation Party"

The authors investigated what happens when these tiny black holes evaporate just before this cooking process begins. When a black hole evaporates, it does not simply vanish quietly; it throws a massive party, shooting out a burst of high-energy particles (such as protons, neutrons, and pions).

Imagine a firework exploding in the middle of a quiet kitchen.

• The Chaos: These particles collide with the "dough" (the protons and neutrons).
• The Swap: Some of these particles act like a mischievous chef, swapping protons for neutrons or vice versa.
• The Result: This changes the recipe. If you have too many neutrons, you end up with too much helium. If you have too few, you get too little.

4. The Big Discovery: A "Goldilocks" Threshold

Previous studies suggested that even tiny black holes (lighter than a mountain) could mess up the recipe. Yet this article says: "Not so fast."

The authors found a strict threshold.

• Too small (under 10⁹ grams): If the black holes are lighter than a large asteroid (about one billion grams), they evaporate too early. They disappear before the "cooking" begins. Their energy is diluted and washed away by the hot soup, leaving no trace on the final bread. It is like throwing a pebble many kilometers upstream into a river; by the time the water reaches the city, the wave has vanished.
• Just right (over 10⁹ grams): Only black holes heavier than this threshold survive long enough to throw their "party" just before cooking begins. Only then can they actually alter the amount of helium in the universe.

5. The "Sweet Spot" of Sensitivity

The authors discovered something interesting about the size of these black holes.

• If they are just heavy enough to play a role, they have a small effect.
• If they become slightly heavier (around 2 billion grams), their ability to spoil the recipe reaches its peak. Here, they are most sensitive.
• If they become even heavier, the effect diminishes again.

It is like tuning a radio: there is a specific frequency where the signal is loudest. The authors found that the "loudest" signal for these black holes occurs at a mass of about 2 billion grams.

6. The Conclusion: The Net is Drawn Tighter

The article concludes that for these black holes to leave a trace in today's universe, they must be heavier than previously assumed.

• The Constraint: The authors calculated that if these black holes existed within the mass range they studied, they must be incredibly rare. Their initial abundance (how many there were relative to the total energy of the universe) had to be less than one part in a hundred quadrillion ($10^{-17}$ to $10^{-19}$).
• The Tool: To ensure their mathematics is solid and verifiable by anyone, they made their computer code publicly available. This allows other scientists to update the recipe with new data in the future.

In summary:
This article is a detailed examination of the universe's "cookbook." It tells us that tiny black holes can only influence the cooking of the first elements if they are heavy enough to survive until the very last second before the oven is turned on. If they are too light, they vanish too early to play a role. The authors have drawn a new, sharper line on the map of possibilities, ruling out many scenarios that previous studies still assumed were open for debate.

Technical Summary

Technical Summary: Primordial Black Holes Evaporating Before Primordial Nucleosynthesis

Problem Statement
Primordial Black Holes (PBHs), formed from density fluctuations in the early universe, offer a unique window into pre-structure formation physics. While PBHs with masses below $\sim 10^{15}$ g have completely evaporated by the present day, those with masses below $\sim 10^9$ g evaporate before primordial nucleosynthesis (BBN). Although these lighter PBHs do not survive to influence later cosmological processes, their Hawking radiation releases particles and entropy that could alter the thermal and particle history of the universe during the BBN era. Previous studies have provided constraints on such PBHs, yet the literature exhibits significant discrepancies regarding mass thresholds and the strength of these constraints. In particular, there is no consensus on whether PBHs with masses below $10^9$ g can leave observable imprints on the abundances of light elements, necessitating a transparent, reproducible re-examination of the calculation steps.

Methodology
The authors present a detailed analytical and numerical framework to evaluate the impact of PBHs evaporating before BBN. The methodology includes a step-by-step "anatomy" of the calculation, where numerical outputs at each stage are cross-checked against analytical estimates:

1. Background Evolution: The authors parameterize the initial PBH abundance via the mass fraction $\beta$. They calculate the evolution of the PBH energy density relative to the standard radiation background and determine the impact on the Hubble expansion rate ($H$) and the effective number of neutrino species ($N_{\text{eff}}$).
2. Hadronization of Hawking Radiation: The study focuses on the hadronic sector of Hawking radiation, as mesons drive the most significant non-standard effects. The authors use the BlackHawk package (v2.3) in combination with analytical estimates to determine meson fluxes ($\pi^\pm, K^\pm, K^0_L$). They address a specific limitation in BlackHawk regarding particle energy bins ($10^5$ GeV) by supplementing it with analytical scaling laws derived from PYTHIA simulations (v8).
3. Meson-Driven $n \leftrightarrow p$ Conversion: The core mechanism analyzed is the interaction of mesons emitted by PBHs with the thermal bath. The authors calculate thermally averaged cross-sections for meson-induced neutron-proton conversions (e.g., $\pi^+ + n \to p + \pi^0$) and include Sommerfeld enhancement factors. These rates are added to the standard weak interaction rates in the Boltzmann equation governing the neutron-to-proton ratio ($X_n$).
4. Nucleon Annihilation: The study accounts for the annihilation of antinucleons produced by PBHs with thermal nucleons. While $p\bar{p}$ annihilation has a negligible net effect, $n\bar{p}$ annihilation effectively drives a conversion of neutrons to protons, altering the evolution of $X_n$.
5. Numerical Implementation: The authors solve the coupled Boltzmann equations for $X_n$ evolution, incorporating the modified conversion rates. They calculate the primordial Helium-4 mass fraction ($Y_P$) both via a simplified approximation ($Y_P \approx 2X_n(T_{\text{nuc}})$) and through a full nuclear reaction network (via BBN-simple). The code is publicly available to ensure reproducibility.

Main Contributions

• Transparent Framework: The work provides a fully reproducible codebase and detailed analytical derivations for every step of the calculation, from Hawking radiation emission to the final $Y_P$ constraint.
• Refined Hadronization: The authors correct and calibrate hadronization fluxes by combining BlackHawk with PYTHIA-based analytical scalings, thereby addressing energy bin issues in existing tools.
• Systematic Resolution of Discrepancies: By rigorously examining each physical step, the authors identify the sources of disagreement with previous literature, particularly regarding the treatment of PBH evaporation scales and hadronic emissivity.

Results

• Mass Threshold: The analysis finds that PBHs must have an initial mass $m_{\text{BH},i} > 10^9$ g to produce observable effects on BBN. PBHs lighter than this threshold evaporate at temperatures $T \gtrsim 1.8$ MeV, where their effects are washed out by rapid thermal interactions prior to neutron freeze-out. This contradicts earlier studies (e.g., Kohri et al. 2000, Keith et al. 2020) that suggested constraints below $10^9$ g.
• Sensitivity Profile: For masses just above $10^9$ g, BBN sensitivity to PBHs increases rapidly with mass, peaks at $m_{\text{BH}} \approx 2 \times 10^9$ g, and then declines.
• Constraints on $\beta$: For the mass range $[10^9, 10^{10}]$ g, current measurements of $Y_P$ (specifically $Y_P = 0.245 \pm 0.003$) set an upper limit on the initial mass fraction $\beta$ in the range of $10^{-17}$ to $10^{-19}$.
• Comparison with Literature: The derived constraints are generally weaker than those in Kohri et al. (2000) and Keith et al. (2020), but stronger than those in Carr et al. (2020) for most of the mass range, except for a specific interval ($m_{\text{BH}} \in [3.2, 7.4] \times 10^9$ g), where the Carr et al. limit is tighter. The most restrictive limit in this work occurs at $m_{\text{BH}} \approx 2 \times 10^9$ g, compared to $6 \times 10^9$ g for Carr et al.
• Background Effects: The study confirms that modifications to the background expansion rate ($H$) and $N_{\text{eff}}$ by PBHs are subdominant to the direct hadronic effects on nucleon abundances for the mass range considered.

Significance
The work asserts that its primary significance lies in providing a robust, transparent, and reproducible framework for BBN constraints on PBHs evaporating before BBN. By resolving discrepancies in the literature and establishing a clear mass threshold of $10^9$ g, the work clarifies the parameter space in which PBHs can be probed via light element abundances. The authors emphasize that their results are derived from a careful disentangling of physical effects (background evolution vs. hadronic interactions) and that the public availability of their code enables future updates with improved nuclear reaction rates and cosmological data. The study explicitly notes that its analysis is restricted to PBHs evaporating before nucleosynthesis; constraints for heavier PBHs involving photo- and hadro-dissociation of light elements are reserved for future work.