Refining Galactic primordial black hole evaporation constraints

This paper presents a comprehensive analysis of cosmic ray signals from asteroid-mass primordial black holes, incorporating full transport effects and new multi-wavelength data to derive leading constraints on their contribution to dark matter while testing various spin and mass distribution assumptions.

The Gist

Imagine the universe is a vast, dark ocean. For decades, scientists have been trying to figure out what makes up the "dark" part of this ocean (Dark Matter). They've looked for tiny, invisible fish (particles), but haven't caught any yet. So, they are starting to wonder: Could the dark matter actually be made of ghostly, ancient black holes that formed right after the Big Bang? These are called Primordial Black Holes (PBHs).

This paper is like a team of detectives (Pedro, Jordan, and Shyam) going back to the crime scene to re-examine the evidence, but this time with much better magnifying glasses and a more realistic understanding of how things move in the galaxy.

Here is the story of their investigation, broken down into simple concepts:

1. The Suspects: Tiny, Dying Black Holes

Most black holes are huge monsters that swallow stars. But these suspects are different. They are "asteroid-mass" black holes—tiny, invisible specks, but incredibly dense.

• The Hawking "Leak": According to physics, black holes aren't truly black; they slowly leak energy and particles, like a hot cup of coffee cooling down. This is called Hawking Radiation.
• The Clue: As these tiny black holes get smaller, they get hotter and leak faster. When they are the size of a mountain or an asteroid, they are leaking so much energy that they shoot out high-speed electrons and positrons (anti-electrons) into space.

2. The Chase: How Particles Travel

In the past, scientists tried to catch these leaking particles by looking straight at the black holes. But the galaxy is messy. It's full of magnetic fields, gas clouds, and solar winds that act like a bouncy castle or a pinball machine.

• The Old Way: Previous studies assumed the particles traveled in straight lines or followed simple rules.
• The New Way: This team used a super-computer simulation (like a high-tech weather forecast for space) to track exactly how these particles bounce, drift, and get re-energized as they travel through the Milky Way. They realized that the "bounciness" of the galaxy changes the clues significantly.

3. The Three Detectives (The Evidence)

The team looked for the "footprints" of these black holes in three different places:

• Detective Voyager 1 (The Local Scout):
  Voyager 1 is a probe that has flown past our Sun's protective bubble (the heliosphere) into deep space. It can catch low-energy particles that the Sun usually blocks.

  • The Finding: The team checked if the particles Voyager 1 saw matched the "leakage" from our asteroid-mass black holes. They found that if there were too many black holes, Voyager 1 would have seen a flood of particles it didn't see. This sets a limit on how many black holes can exist.

• Detective Xmm-Newton (The X-Ray Eye):
  When the high-speed particles from the black holes hit the "fog" of the galaxy (starlight and gas), they scatter and glow in X-rays.

  • The Finding: The team looked at X-ray maps of the center of our galaxy. They realized that if there were too many black holes, the galaxy would be glowing much brighter in X-rays than it actually is.
  • The Twist: The paper includes an Erratum (a correction). They realized they initially overestimated how much X-ray light they should see because of a math error regarding the size of the telescope's view. After fixing this, the X-ray evidence is less strict than they first thought, but still important.

• Detective Integral (The 511 keV Line):
  When a positron (anti-electron) meets a normal electron, they annihilate each other and explode in a flash of light with a very specific color: 511 keV.

  • The Finding: There is a massive glow of this specific light coming from the center of our galaxy. The team asked: "Could this glow be caused by our black holes?"
  • The Result: They found that if the black holes were too common, the glow would be much brighter and shaped differently than what we see. This turned out to be their strongest evidence, ruling out a huge chunk of the "asteroid" black hole population.

4. The Verdict

The team didn't just say "No black holes." They said, "We can't have that many."

• The "Goldilocks" Zone: They found that while black holes could make up some of the dark matter, they cannot make up all of it in the mass range they studied (between the mass of a small asteroid and a large mountain).
• Spin Matters: They also tested if the black holes were spinning like tops (Kerr black holes) or sitting still (Schwarzschild). Spinning black holes leak more energy, which makes the "no-go" zones even stricter.
• The Uncertainty: They emphasized that our understanding of how particles bounce around the galaxy (the "pinball machine") is still a bit fuzzy. If the galaxy is "bouncier" than we think, the limits change. But even with these uncertainties, the evidence is strong.

The Big Picture

Think of this paper as refining a map of a treasure hunt.

• Before: "The treasure (Dark Matter) might be hidden in this whole forest."
• Now: "We've checked the trees, the rivers, and the wind. We know the treasure isn't hidden in this specific part of the forest (asteroid-mass black holes). It might be elsewhere, or it might be a different kind of treasure entirely."

By using better math, real data from space probes, and correcting their own mistakes, the authors have tightened the noose around the idea that tiny, ancient black holes are the main ingredient of Dark Matter. They haven't ruled it out completely, but they've made it much harder for that theory to survive.

Technical Summary

1. Problem Statement

Primordial Black Holes (PBHs) remain a compelling candidate for Dark Matter (DM), particularly in the "asteroid-mass" range ($10^{14} - 10^{17}$ g). While gravitational lensing constrains heavier PBHs, lighter PBHs are difficult to detect via lensing due to finite-source effects. However, PBHs in this mass range are expected to emit Hawking radiation, producing high-energy particles (electrons, positrons, photons, neutrinos) that can be detected indirectly.

Previous studies have faced limitations:

• Propagation Simplifications: Many analyses used semi-analytical codes or simplified propagation models that did not fully account for cosmic ray (CR) reacceleration and diffusion effects within the Milky Way.
• Data Scope: Constraints often relied on single probes (e.g., only the 511 keV line or only X-rays) without a unified treatment of the full CR transport chain.
• Systematic Uncertainties: The impact of varying PBH spin distributions (Schwarzschild vs. Kerr), mass distributions (monochromatic vs. log-normal), and propagation parameters (Alfvén velocity, halo height) was not comprehensively quantified.
• Data Re-evaluation: Recent re-analyses of XMM-Newton data and 511 keV positronium formation calculations have highlighted potential overestimations in previous flux predictions.

2. Methodology

The authors perform a comprehensive, numerical analysis of PBH evaporation signals, integrating state-of-the-art astrophysical ingredients.

• Hawking Evaporation Modeling:

  • Used the BlackHawk v2.2 code to compute primary particle emission spectra.
  • Incorporated Hazma for secondary processes (hadronization, decay of $\mu^\pm, \pi^\pm$, and final state radiation) valid for sub-GeV energies ($M_{PBH} \gtrsim 10^{14.5}$ g).
  • Investigated two extreme spin benchmarks: Schwarzschild ($a_* = 0$) and near-extremal Kerr ($a_* = 0.9999$).
  • Modeled mass distributions as both monochromatic and log-normal (varying standard deviation $\sigma$).

• Cosmic Ray Propagation:

  • Solved the diffusion-advection-convection-loss equation using the DRAGON2 numerical code.
  • Explicitly modeled reacceleration (momentum diffusion), which is crucial for low-energy $e^\pm$ produced by PBHs, potentially boosting them to GeV energies.
  • Tested multiple propagation scenarios to quantify uncertainties:
    • Benchmark: $H = 8$ kpc, $V_A = 13.4$ km/s.
    • Realistic: Variations within $3\sigma$ of benchmark fits.
    • General/Conservative: Extreme ranges ($H \in [3, 16]$ kpc, $V_A \in [0, 40]$ km/s).

• Observational Probes:

  1. Local $e^\pm$ Flux: Compared predicted spectra against Voyager 1 data (which is outside the heliosphere and unaffected by solar modulation).
  2. Diffuse X-ray Emission: Analyzed XMM-Newton data in the 2.5–8 keV band. Calculated inverse Compton scattering of Galactic ambient photons by PBH-produced $e^\pm$.
     • Correction (Erratum): The authors corrected a geometric factor error in the solid angle calculation for XMM-Newton Ring 3, significantly reducing the predicted flux and weakening the constraints.
  3. 511 keV Line Profile: Compared the longitudinal profile of the 511 keV emission (from positron annihilation) with INTEGRAL/SPI data.
     • Correction (Erratum): Updated the calculation to only count thermalized positrons, excluding high-energy positrons that annihilate in-flight, leading to more conservative limits.

3. Key Contributions

• Full Numerical Transport: First study to apply a fully numerical CR transport code (DRAGON2) to PBH evaporation signals, capturing the full impact of reacceleration on the $e^\pm$ spectrum.
• Multi-Probe Synthesis: Simultaneously constrains PBH DM fraction ($f_{PBH}$) using Voyager 1, XMM-Newton, and INTEGRAL/SPI within a single, consistent framework.
• Parameter Space Exploration: Systematically quantified how spin ($a_*$), mass distribution width ($\sigma$), and propagation uncertainties affect the derived limits.
• Data Corrections: Provided an erratum correcting geometric normalization errors in XMM-Newton analysis and thermalization assumptions in 511 keV calculations, ensuring the robustness of the final constraints.

4. Results

• Voyager 1 Constraints:
  • The local $e^\pm$ flux is highly sensitive to reacceleration. Without reacceleration, constraints on high-mass PBHs ($>10^{16}$ g) vanish because the injected particles are too low energy to reach the detector.
  • With standard reacceleration, Voyager 1 provides competitive constraints for $M_{PBH} \lesssim 10^{16}$ g.
• XMM-Newton Constraints (Corrected):
  • After correcting the solid angle error, the X-ray constraints weakened substantially (by orders of magnitude) compared to the initial draft.
  • However, they remain relevant for specific mass ranges and are sensitive to the inner Galaxy DM profile.
• 511 keV Constraints (Corrected):
  • The corrected 511 keV limits are the most stringent for PBH masses between $10^{16}$ g and $2 \times 10^{17}$ g.
  • The constraints are robust against changes in the DM density profile (NFW vs. Moore vs. cNFW) because the high-longitude data points (which drive the limits) are less sensitive to the central cusp.
• Impact of PBH Properties:
  • Spin: Near-extremal Kerr black holes ($a_* \approx 1$) produce harder spectra and higher fluxes, leading to tighter constraints than Schwarzschild black holes.
  • Mass Distribution: A log-normal distribution with $\sigma > 0$ introduces a population of lower-mass PBHs. Since lower-mass PBHs evaporate more efficiently, this significantly strengthens the constraints compared to a monochromatic distribution.
• Uncertainties:
  • Propagation uncertainties (specifically $V_A$ and $H$) can shift the derived $f_{PBH}$ limits by up to an order of magnitude.
  • The 511 keV constraints are the most robust against propagation model variations.

5. Significance and Conclusion

This work establishes the current leading constraints on the fraction of Dark Matter that can be composed of asteroid-mass Primordial Black Holes.

• Exclusion Limits: The study concludes that PBHs can constitute a significant fraction of DM only in the mass gap between $10^{18} - 10^{21}$ g (for monochromatic distributions) or $\sim 10^{19} - 10^{21}$ g (for log-normal distributions). The region below $10^{18}$ g is tightly constrained, effectively ruling out PBHs as the sole component of DM in this range.
• Methodological Benchmark: The paper sets a new standard for indirect detection studies by demonstrating the critical importance of full numerical propagation (especially reacceleration) and rigorous handling of systematic errors in observational data.
• Future Outlook: The authors highlight that future data from eROSITA (Galactic diffuse X-rays) and improved 21-cm observations could further tighten these bounds, particularly in the $10^{16} - 10^{17}$ g range.

In summary, by refining the transport models and correcting data analysis errors, the authors provide a conservative yet robust exclusion of PBHs as the dominant Dark Matter candidate in the $10^{14} - 10^{18}$ g mass range.