🔭 astrophysics

Primordial black holes induced stochastic axion-photon oscillations in primordial magnetic field

This paper investigates the stochastic oscillations between axion-like particles emitted by ultra-light primordial black holes and photons within primordial magnetic fields, analyzing their probability distributions and potential impacts on the cosmic microwave, X-ray, and gamma-ray backgrounds.

Original authors: Hai-Jun Li

Published 2026-01-27

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Hai-Jun Li

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: Invisible Ghosts and Tiny Black Holes

Imagine the very early universe as a chaotic, high-energy party. In this paper, the author, Hai-Jun Li, explores a specific scenario involving two mysterious characters: Primordial Black Holes (PBHs) and Axion-Like Particles (ALPs).

Think of Primordial Black Holes as tiny, invisible specks of darkness that formed right after the Big Bang. They are so small (ranging from the weight of a mountain to the weight of a grain of sand) that they have already "evaporated" and disappeared by now. Because they vanished before the universe cooled down enough to form atoms, we can't see them directly.

However, as these tiny black holes evaporated, they didn't just disappear; they spat out energy. One of the things they might have spit out are ALPs. Think of ALPs as "ghost particles." They are incredibly light, almost massless, and they are very shy. They don't like to interact with normal matter, but they have a secret handshake with photons (particles of light).

The Main Event: The Cosmic Dance

The core of this paper is about what happens when these "ghost" ALPs travel through the universe and encounter a Primordial Magnetic Field (PMF).

To understand this, imagine the universe is filled with a giant, invisible ocean of magnetic force.

The ALPs are like dancers who can only move in one way.
The Photons are dancers who move in a different way.
The Magnetic Field is the music.

When the ALP dancers hear the magnetic music, they can magically transform into Photon dancers, and vice versa. This is called oscillation. It's like a dancer switching costumes mid-performance.

The paper asks: If these ALPs were created by the evaporating tiny black holes, how much of them would turn into light (photons) as they travel through this cosmic magnetic ocean?

The Two Scenarios: A Smooth Road vs. A Bumpy Path

The author investigates two different ways the "magnetic music" (the magnetic field) could be arranged in the universe:

The Smooth Road (Homogeneous Field): Imagine the magnetic field is like a perfectly flat, calm lake. The water is still, and the direction is the same everywhere. In this scenario, the ALPs turn into light in a very predictable, rhythmic pattern.
The Bumpy Path (Stochastic Field): This is the scenario the author focuses on more. Imagine the magnetic field is like a choppy sea with waves going in random directions. The strength and direction of the magnetic "wind" change randomly as you move through space. In this case, the ALPs' transformation into light becomes a stochastic (random) process. It's like trying to dance on a ship in a storm; the outcome is less predictable and depends on the specific "waves" you hit.

The Rules of the Game (The Limits)

The author doesn't just make up numbers; they use the strictest rules available from the Planck satellite data (specifically a 2019 study). These rules say: "The magnetic field in the early universe cannot be stronger than a certain tiny amount."

If the field is smooth, it must be weaker than a specific limit (about 47 pG).
If the field is choppy/random, it must be even weaker (about 8.9 pG).

The author runs computer simulations using these strict limits to see how often ALPs turn into photons under these conditions.

What Did They Find?

The paper presents a lot of charts and numbers, but the main takeaways are:

Energy Matters: The chance of an ALP turning into a photon depends heavily on how much energy the particle has. At very high energies (like those found in cosmic rays), the transformation becomes much less likely.
Strength Matters: A stronger magnetic field makes the transformation happen more often. However, because the universe's magnetic field limits are so strict, the transformation probability is generally quite low.
Randomness Matters: In the "choppy sea" (stochastic) scenario, the results vary wildly depending on the specific random arrangement of the magnetic fields. The author shows that while some paths might allow for a high chance of transformation, others might block it almost entirely.

Why Does This Matter? (According to the Paper)

The paper concludes that even though these tiny black holes are gone, the "ghost particles" (ALPs) they created might still be around. If these ALPs turned into photons while traveling through the early universe's magnetic fields, they could leave a faint fingerprint on the light we see today.

The author suggests that this process could affect:

The Cosmic Microwave Background (CMB): The afterglow of the Big Bang.
The Cosmic X-ray Background: The faint glow of X-rays filling the sky.
The Extragalactic Gamma-ray Background: High-energy light from deep space.

Essentially, the paper argues that if we look closely at the background light of the universe, we might be able to see the "echo" of these tiny black holes and their ghostly particles, provided we can detect the subtle changes caused by this cosmic dance in the magnetic fields.

Summary

The paper is a theoretical investigation. It calculates how likely it is for "ghost particles" (ALPs), created by the death of tiny ancient black holes, to turn into light as they travel through the random, magnetic turbulence of the early universe. It uses strict observational limits to show that while this happens, the probability is generally low and depends heavily on the energy of the particles and the specific "randomness" of the magnetic fields they encounter.

Technical Summary: Primordial Black Hole Induced Stochastic Axion-Photon Oscillations in Primordial Magnetic Field

Problem Statement
Primordial Black Holes (PBHs) formed in the very early Universe due to large density fluctuations can serve as non-thermal sources of Axion-Like Particles (ALPs). While PBHs with masses $M_{\text{PBH}} \gtrsim 10^{15}$ g may constitute dark matter, lighter PBHs ( $10 \text{ g} \lesssim M_{\text{PBH}} \lesssim 10^9 \text{ g}$ ) would have completely evaporated via Hawking radiation before the onset of Big Bang Nucleosynthesis (BBN). Consequently, these ultra-light PBHs evade direct constraints but could leave observable imprints through the ALPs they emit. The central problem addressed is the propagation and conversion of these PBH-induced ALPs into photons as they traverse the Primordial Magnetic Field (PMF). Specifically, the paper investigates the stochastic nature of these oscillations, moving beyond simplified homogeneous magnetic field models to account for the non-homogeneous, stochastic background of the cosmic plasma.

Methodology
The study employs a multi-stage theoretical framework:

PBH Formation and Evaporation: The authors model PBH formation in the radiation-dominated epoch, establishing a mass range of $10 \text{ g} \lesssim M_{\text{PBH}} \lesssim 10^9 \text{ g}$ based on evaporation timescales relative to BBN and inflation constraints. They calculate the emission spectra of spin-0 particles (ALPs) from both Schwarzschild and Kerr black holes using the public code BlackHawk. This includes deriving instantaneous spectra and mass loss rates, accounting for greybody factors and the effects of black hole spin ( $a_*$ ) on the emission.
ALP-Photon Interaction: The interaction is modeled under the minimal scenario where ALPs interact solely with photons via the coupling term $\mathcal{L} \supset -\frac{1}{4}g_{a\gamma}a F_{\mu\nu}\tilde{F}^{\mu\nu}$ . The propagation of the ALP-photon system is described by a mixing matrix formalism, incorporating plasma effects ( $\Delta_{pl}$ ), QED vacuum polarization ( $\Delta_{QED}$ ), and CMB photon dispersion ( $\Delta_{CMB}$ ).
Magnetic Field Modeling: Two PMF scenarios are analyzed:
- Homogeneous Field: A simplified model with a constant field strength, used for comparison.
- Stochastic Field: The primary focus, modeled as a stochastic background with non-helical and helical components. The magnetic field is simulated as a domain-like structure where the field strength is constant within domains but the direction varies randomly between domains.
Constraints: The analysis utilizes the latest stringent limits on PMF derived from CMB radiation (Jedamzik & Saveliev 2019, JS 2019), specifically the scale-invariant limit ( $B < 47$ pG for $n=0$ ) and the violet Batchelor spectrum limit ( $B < 8.9$ pG for $n=5$ ). These are compared against older Planck 2015 limits ( $B \sim 5$ nG).
Numerical Simulation: The oscillation probability is calculated using a density matrix approach and a transport matrix formalism. For the stochastic case, the authors simulate 1000 random magnetic field configurations over 1 Mpc (divided into 100 domains) to derive probability distributions and confidence intervals ( $1\sigma$ and $2\sigma$ ).

Key Contributions and Results

Spectral Characterization: The paper provides instantaneous emission spectra for ALPs from Schwarzschild and Kerr PBHs across the mass range $10^5$ g to $10^9$ g. It demonstrates that for Kerr black holes, increasing spin shifts the high-energy spectral distribution to the right due to the competition between the greybody factor and the Boltzmann factor.
Oscillation Probability Distributions:
- Homogeneous Scenario: Numerical results show that the ALP-photon oscillation probability ( $P_{a\gamma}$ ) in the $(m_a, g_{a\gamma})$ parameter space is highly sensitive to energy and magnetic field strength. The probability is significantly suppressed at energies above $\sim 10^4$ GeV. A comparison between the JS 2019 limits ($47$ pG) and Planck 2015 limits ($5$ nG) reveals substantial differences in the predicted probability distributions.
- Stochastic Scenario: In the stochastic magnetic field model (using the $B < 8.9$ pG limit), the oscillation probability exhibits distinct behaviors depending on ALP parameters. The coupling constant $g_{a\gamma}$ dramatically enhances the probability, while the ALP mass $m_a$ primarily affects the probability in low-energy regions, with the probability stabilizing at high energies. The study presents $1\sigma$ and $2\sigma$ contours derived from 1000 random configurations, illustrating the statistical spread of oscillation probabilities.
Redshift and Absorption Effects: The paper outlines the necessary modifications to the mixing matrix to account for cosmological redshift and photon absorption by hydrogen and helium during recombination, noting that these effects would alter the present-day photon flux, though detailed numerical results for these specific corrections are not the primary output of this work.

Significance and Claims
The authors claim that their work provides a detailed investigation into the stochastic oscillations between ALPs and photons induced by evaporating ultra-light PBHs within the context of the most stringent current PMF limits. By moving from homogeneous to stochastic magnetic field models, the paper offers a more realistic assessment of the ALP-photon conversion probability in the early Universe.

The paper posits that these PBH-induced stochastic oscillations could have observable effects on several cosmological backgrounds, specifically the Cosmic Microwave Background (CMB), the Cosmic X-ray Background (CXB), and the Extragalactic Gamma-ray Background (EGB). The authors suggest that these effects could potentially be detected as signals in next-generation experiments such as ADMX, IAXO, ALPS-II, LiteBIRD, CMB-S4, THESEUS, eXTP, LHAASO, and CTA. However, the paper explicitly states that it focuses on the calculation of the oscillation probabilities and does not investigate the specific phenomenological consequences or observational signatures in these backgrounds in detail. The work serves as a foundational step for future studies linking PBH evaporation, ALP physics, and primordial magnetism.