Constraints on Primordial Black Holes from Galactic Diffuse Synchrotron Emissions

This paper demonstrates that Galactic diffuse synchrotron emissions, particularly when accounting for diffusive re-acceleration of Hawking-radiated electrons and positrons, provide stringent constraints on the abundance of primordial black holes with masses above $10^{15}$ g that are significantly stronger than those derived from current cosmic-ray electron and positron data.

The Gist

The Big Picture: Hunting Ghosts with Radio Waves

Imagine the universe is filled with invisible "ghosts" called Dark Matter. We know they are there because they have gravity, but we can't see them. One popular theory suggests these ghosts might be Primordial Black Holes (PBHs). These aren't the black holes formed by dying stars; they are tiny, ancient black holes that formed in the very first split-second after the Big Bang.

If these tiny black holes exist, they shouldn't be silent. According to a famous theory by Stephen Hawking, they should be slowly "evaporating" by shooting out particles, like a very slow, cold leak. The authors of this paper wanted to see if we could find these leaks by listening to the radio static of our galaxy.

The Problem: The Particles are Too "Lazy"

Here is the catch:

1. The Leak: Tiny black holes (around the mass of a mountain) shoot out electrons and positrons (anti-electrons).
2. The Energy: When they first shoot out, these particles are "lazy." They have very low energy (around 10 MeV).
3. The Limitation: If you just look at the radio waves these lazy particles make, they are too weak to be seen against the background noise of the galaxy. It's like trying to hear a whisper in a hurricane.

The Solution: The "Cosmic Treadmill"

The paper's main discovery relies on a specific theory about how particles move through the galaxy, called Diffusive Re-acceleration.

Think of the galaxy's magnetic field as a giant, chaotic ocean.

• The Old View: Particles just drift through this ocean, slowly losing energy.
• The New View (This Paper): The ocean isn't calm; it's turbulent. The magnetic fields are constantly jiggling and moving (like waves crashing). When the lazy electrons from the black holes hit these moving magnetic waves, they get a boost.

The authors used data from the AMS-02 (a particle detector on the International Space Station) and Voyager-1 (a probe at the edge of our solar system) to prove that this "turbulent ocean" theory is correct. They found that the magnetic waves are strong enough to act like a cosmic treadmill, kicking the lazy electrons up to much higher speeds (around 100 MeV).

The Result: Turning Up the Volume

Once these electrons get that speed boost from the magnetic treadmill, they start spinning around the galaxy's magnetic fields much faster. When charged particles spin fast, they emit synchrotron radiation—which is basically a bright radio signal.

• Before the boost: The radio signal was a whisper.
• After the boost: The radio signal is a shout.

The authors used radio telescopes that listen to frequencies between 22 MHz and 1.4 GHz (low-frequency radio waves) to look for this shout.

The "No-Go" Zone: Setting the Limits

The researchers didn't find a specific "smoking gun" signal that says, "Here is a black hole!" Instead, they did something even more powerful: They set a limit.

They calculated: "If there were this many black holes, the radio shout would be so loud that it would drown out the natural background noise of the galaxy."

Since the radio telescopes don't hear a shout that loud, the authors can say: "There cannot be more than X amount of these black holes."

Key Findings:

• Better than before: Their new limits are much stricter (stronger) than previous attempts. For example, for black holes heavier than a certain mass, their new rules are 10 times stricter than what we knew from looking at just the Voyager-1 data.
• The Sweet Spot: This method works best for black holes that are heavy enough to still exist today but light enough to still be evaporating particles.
• The "Conservative" Approach: The authors were very careful. They didn't try to subtract the background noise perfectly. They assumed all the radio noise we see could be from black holes. Even with this super-cautious approach, they still managed to rule out huge amounts of black holes.

The Takeaway

This paper is like a detective saying: "We know the suspect (the black hole) leaves a specific kind of footprint (radio waves). We checked the crime scene (the galaxy's radio sky) with a very sensitive microphone. We didn't hear the suspect's footsteps loud enough to match our theory. Therefore, the suspect cannot be hiding in the crowd in the numbers we thought."

By proving that the galaxy's magnetic fields act like a treadmill that speeds up these particles, the authors turned a faint, unhearable whisper into a loud radio signal, allowing us to set much tighter rules on how many of these ancient black holes can exist in our universe.

Technical Summary

1. Problem Statement

The nature of Dark Matter (DM) remains unknown. Primordial Black Holes (PBHs) are a leading candidate, particularly in the mass range $M_{\text{PBH}} \gtrsim 10^{15}$ g, where they are stable against Hawking evaporation but light enough to emit significant Standard Model particles.

• The Challenge: PBHs in this mass range emit electrons and positrons (all-electrons) with initial energies of $\mathcal{O}(10 \text{ MeV})$ via Hawking radiation. Direct detection of these low-energy particles is difficult due to solar modulation and the dominance of astrophysical backgrounds.
• The Gap: Previous constraints on PBH abundance ($f_{\text{PBH}}$) from cosmic-ray (CR) data (e.g., Voyager-1, AMS-02) were limited by the inability to detect low-energy particles or by the assumption of standard CR propagation models without significant re-acceleration.
• The Opportunity: If CR propagation involves diffusive re-acceleration driven by a significant Alfvén velocity ($V_a \sim \mathcal{O}(10)$ km/s), the initial $\sim 10$ MeV electrons from PBHs can be boosted to $\sim 100$ MeV. These higher-energy electrons emit synchrotron radiation in the Galactic Magnetic Field (GMF) at radio frequencies ($\sim 20$ MHz to 1.4 GHz), offering a new indirect detection channel.

2. Methodology

A. PBH Evaporation and Source Modeling

• Hawking Radiation: The authors use the BlackHawk code to calculate the energy spectra of primary and secondary all-electrons emitted by PBHs.
• Mass Functions: Two scenarios are considered:
  1. Monochromatic: All PBHs have a single mass $M_c$.
  2. Log-normal: A distribution centered at $M_c$ with width $\sigma$.
• DM Profiles: Both NFW (cuspy) and Burkert (cored) DM density profiles are tested to account for uncertainties in the Galactic center density.

B. Cosmic Ray Propagation

• Propagation Models: The authors employ the Galprop code to solve the diffusion equation for CRs in the Galaxy. They fit parameters to AMS-02 and Voyager-1 Boron-to-Carbon (B/C) flux ratio data to establish benchmark models.
• Key Models:
  • Diffusive Re-acceleration Models (DRz4-10): Four models with different diffusion halo half-heights ($z_h = 4, 6, 8, 10$ kpc) and a free Alfvén velocity parameter ($V_a$).
  • Diffusion Break Model (DBz4): A model without re-acceleration ($V_a=0$) but with a low-energy break in the diffusion coefficient, serving as a control.
  • GH Model: A model combining Galprop and Helmod (for solar modulation) with parameters fitted in previous literature.
• Fit Results: The fits confirm that $V_a \sim 20$ km/s is favored in the re-acceleration models. This velocity is sufficient to boost evaporated PBH electrons from $\sim 10$ MeV to $\sim 100$ MeV during propagation.

C. Synchrotron Emission Calculation

• GMF Models: Two realistic Galactic Magnetic Field models are used: SUNE and JF12. These include regular, striated random, and isotropic random field components.
• Radiative Processes: The calculation includes synchrotron emission, free-free absorption (significant below $\sim 100$ MHz), and Faraday rotation.
• Observational Data: Constraints are derived by comparing predicted synchrotron intensities against radio continuum survey data from 22 MHz to 1.4 GHz (e.g., 22 MHz, 45 MHz, 408 MHz, 1420 MHz).
• Conservative Approach: To ensure robustness, the authors do not subtract astrophysical backgrounds. They assume the entire observed radio intensity could potentially be attributed to PBHs, setting limits where the predicted PBH signal does not exceed the observed intensity plus $2\sigma$ errors.

3. Key Contributions

1. Validation of Re-acceleration Impact: The paper quantitatively demonstrates that diffusive re-acceleration ($V_a \sim 20$ km/s) significantly enhances the flux of PBH-induced electrons at $\sim 100$ MeV, making them detectable via low-frequency radio synchrotron emission.
2. New Constraint Channel: It establishes Galactic diffuse synchrotron emission as a powerful probe for PBHs in the mass range $10^{15} \text{ g} \lesssim M_{\text{PBH}} \lesssim 10^{17}$ g, a range previously less constrained by radio data.
3. Benchmark Model Fitting: The authors provide a set of rigorously fitted CR propagation models (DRz4-10, DBz4, GH) that reproduce B/C data, explicitly quantifying the impact of the diffusion halo height ($z_h$) and Alfvén velocity on PBH constraints.
4. Robustness Analysis: The study systematically tests the sensitivity of constraints to:
   • Different DM profiles (NFW vs. Burkert).
   • Different GMF models (SUNE vs. JF12).
   • Different PBH mass functions (Monochromatic vs. Log-normal).
   • Magnetic field strength uncertainties.

4. Results

• Stringent Constraints: In the diffusive re-acceleration scenario, the constraints on the PBH fraction $f_{\text{PBH}}$ are extremely tight.
  • For $M_{\text{PBH}} \gtrsim 1 \times 10^{16}$ g, the constraints are more than one order of magnitude stronger than those derived from Voyager-1 all-electron data.
  • For $M_{\text{PBH}} \gtrsim 2 \times 10^{16}$ g, the constraints are stronger than previous limits derived from AMS-02 positron data.
• Mass Dependence:
  • Monochromatic Mass Function: The strongest constraints occur around $M_c \sim 10^{16}$ g. For example, using the GH model, $f_{\text{PBH}} < 1.25 \times 10^{-4}$ for $M_c = 1.9 \times 10^{16}$ g.
  • Log-normal Mass Function: Constraints are even stronger for broader distributions ($\sigma \ge 1$) at higher masses because the extended mass function retains a significant contribution from lighter PBHs, which emit more high-energy electrons.
• Model Dependence:
  • Re-acceleration vs. No Re-acceleration: The diffusion break model (DBz4, $V_a=0$) yields much weaker constraints (only constraining $M_c < 10^{15}$ g). This highlights that the existence of diffusive re-acceleration is crucial for the strength of these limits.
  • Halo Height ($z_h$): Constraints strengthen by approximately one order of magnitude as $z_h$ increases from 4 kpc to 10 kpc, as a larger halo contains more PBHs contributing to the signal.
• Frequency Sensitivity: The strongest constraints generally come from lower frequencies (e.g., 22 MHz for SUNE, 45–150 MHz for JF12) because the PBH synchrotron spectrum is softer than the background, and low-frequency observations probe the $\sim 100$ MeV electron population most effectively.

5. Significance

• Competitive Limits: This work provides some of the most stringent limits to date on PBHs in the $10^{15}-10^{17}$ g mass range, surpassing previous direct CR measurements and competing with CMB and gamma-ray constraints.
• Testing CR Physics: The results are critically dependent on the assumption of diffusive re-acceleration. If future observations (e.g., by X-ray missions like e-ASTROGAM) disprove the existence of significant re-acceleration, these specific constraints would weaken. Conversely, the strength of these radio constraints serves as a consistency check for the re-acceleration hypothesis.
• Future Prospects: The authors note that current limits are primarily constrained by the astrophysical radio background and its spatial variance. Future radio facilities like SKA-low and improved LOFAR surveys, with higher angular resolution and better background subtraction capabilities, could improve these constraints by another order of magnitude. Additionally, extending observations below the atmospheric cutoff (10 MHz) would probe even lower electron energies, potentially constraining heavier PBHs.

In summary, this paper leverages the synergy between updated CR propagation models (favoring re-acceleration) and low-frequency radio surveys to place the most stringent constraints yet on the abundance of primordial black holes in the sub-solar mass regime.