Nucleosynthesis and CMB bounds on photophilic ALPs: a fresh look

This paper provides a model-independent re-evaluation of cosmological constraints on photophilic axion-like particles with masses below 10 GeV and lifetimes under $10^4$ seconds by incorporating rare hadronic decays and reheating temperature dependencies, revealing extended bounds and potential parameter space to resolve tensions in $N_{\rm eff}$ and deuterium abundance.

The Gist

The Big Picture: Searching for Invisible Ghosts

Imagine the early universe as a giant, chaotic party that happened just after the Big Bang. For a long time, scientists have been looking for "Axion-Like Particles" (ALPs). Think of these ALPs as invisible ghosts that might have been hanging out at this party.

These ghosts have a special trick: they love to talk to light (photons). If they exist, they might have been created in huge numbers during the early universe, lived for a while, and then vanished, turning back into light.

This paper is a fresh investigation into what happens if these ghosts show up at the party. The authors, Miguel Escudero Abenza, Clara Garcia-Perez, and Maksym Ovchynnikov, are asking: If these ghosts were there, how would they have messed up the "birth certificates" of the universe?

The Two Main Clues: The Baby Photos

To figure out if these ghosts were there, the scientists look at two "baby photos" of the universe:

1. Big Bang Nucleosynthesis (BBN): This is the moment when the universe was a hot soup, and the first atomic nuclei (like Helium and Deuterium) were cooked up. It's like the universe's first kitchen.
2. The Cosmic Microwave Background (CMB): This is the "afterglow" of the Big Bang, a faint light that fills the universe today. It's like a snapshot of the universe when it was a toddler.

The scientists are checking these photos to see if the "ghosts" left any fingerprints.

The New Twist: The Rare "Side-Effect"

In the past, scientists mostly assumed these ghosts only turned into pure light (two photons) when they died. It was like a ghost popping and turning into a flash of light.

However, this paper says: "Wait a minute! Sometimes, these ghosts don't just turn into light. They can also turn into a tiny, unstable particle called a meson (specifically pions)."

Think of it this way:

• Old View: The ghost pops and turns into a harmless flashbulb.
• New View: The ghost usually turns into a flashbulb, but rarely (maybe 1 in 10,000 times), it turns into a tiny, angry firecracker (the meson).

Why does this matter? Because while the flashbulb just adds light, the firecracker can actually change the chemistry of the soup.

The Kitchen Chaos: How the Firecrackers Change the Recipe

The universe's "kitchen" (BBN) was very sensitive. It needed a perfect balance of ingredients to make the right amount of Helium and Deuterium.

1. The Neutron-Proton Swap: In the early universe, protons and neutrons were constantly swapping places.
2. The Firecracker Effect: When those rare meson "firecrackers" exploded, they interacted with the protons and neutrons. They acted like a chef's hand reaching into the soup, forcing more neutrons to turn into protons (or vice versa).
3. The Result: This changed the recipe. Instead of the standard amount of Helium, the universe might have ended up with too much or too little.

The authors found that even though these "firecracker" decays are rare, they are so powerful that they can ruin the recipe much faster than previously thought. This means we have to rule out a much larger area of "ghost territory" than before.

The "Reheating" Temperature: How Hot Was the Party?

The paper also looks at how hot the universe was when the party started (the "reheating temperature").

• High Temperature: If the party was super hot, the ghosts were everywhere, and the constraints on them are very strict.
• Low Temperature: If the party was cooler, fewer ghosts were made.

The authors discovered something interesting: Even if the party was only moderately hot, those rare "firecracker" decays still leave a mark. They found a specific "island" of ghost properties that was previously thought to be safe (allowed) but is now forbidden because of these rare decays.

The Silver Lining: Fixing a Small Glitch

While the main goal was to find where these ghosts cannot exist, the authors also found a small "sweet spot."

There are currently two tiny puzzles in cosmology:

1. The measured number of neutrino types is slightly lower than expected.
2. The measured amount of Deuterium is slightly higher than some theories predict.

The paper suggests that if these ALPs exist with very specific properties (a specific mass and lifetime), they could simultaneously fix both puzzles. It's like finding a single key that unlocks two different doors. While this isn't a proven fact, it's a fascinating possibility that the authors highlight.

The Toolkit: A Public Recipe Book

Finally, the authors didn't just do the math; they built a digital kitchen simulator (a computer code called BBNEasyALP). They made this code publicly available on GitHub.

This means any other scientist can download their "recipe book," plug in their own theories about different types of particles, and see if those particles would have ruined the universe's baby photos.

Summary

• The Subject: Invisible particles (ALPs) that love light.
• The Discovery: Even rare decays into "firecracker" particles (mesons) have a huge impact on the early universe's chemistry.
• The Result: We have to banish these particles from a much larger area of the universe's history than we thought before.
• The Bonus: There is a tiny, specific region where these particles might actually help solve small mysteries in our current data.
• The Gift: The authors shared their computer code so others can test their own ideas against these new, stricter rules.

Technical Summary

Technical Summary: Nucleosynthesis and CMB bounds on photophilic ALPs: a fresh look

Problem Statement
This work revisits the cosmological constraints on Axion-Like Particles (ALPs) that couple predominantly to photons, specifically focusing on the mass range $m_a \lesssim 10$ GeV and lifetimes $\tau_a \lesssim 10^4$ s. While previous studies (e.g., Refs. [23, 26]) have scanned this parameter space, they often relied on older observational data, assumed infinite reheating temperatures ($T_{\text{reh}}$), or neglected specific decay channels. Crucially, prior analyses frequently treated ALPs as decaying exclusively into photons, ignoring rare but cosmologically significant hadronic decay modes into light mesons (pions, kaons) that occur via off-shell photons when $m_a > 2m_{\pi^\pm}$. This paper aims to fill these gaps by incorporating precision cosmological data, exploring the dependency on unknown reheating temperatures, and rigorously accounting for the impact of hadronic decays on Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB).

Methodology
The authors developed a comprehensive pipeline to model the evolution of the early Universe in the presence of photophilic ALPs:

1. ALP Production: The evolution of the ALP abundance ($Y_a = n_a/s$) is calculated by integrating the Boltzmann equation from the reheating temperature ($T_{\text{reh}}$) down to $T = 20$ MeV. Production mechanisms include Primakoff scattering ($\gamma + f^\pm \to a + f^\pm$) and photon fusion ($\gamma\gamma \to a$). The production rate is calculated using Maxwell-Boltzmann approximations for Standard Model fermions, validated against asymptotic results.
2. Decay Widths and Branching Ratios: While the dominant decay is $a \to \gamma\gamma$, the authors calculate sub-dominant hadronic decay widths ($a \to \gamma + \text{hadrons}$) using a data-driven approach. They relate the partial widths to the experimentally measured $R$-ratio ($\sigma_{e^+e^- \to \text{hadrons}}/\sigma_{e^+e^- \to \mu^+\mu^-}$) via the "Dalitz trick," effectively summing over exclusive hadronic final states (e.g., $\pi^+\pi^-$, $K^+K^-$).
3. Thermodynamics and Expansion: The coupled evolution of the electromagnetic plasma, neutrino bath, and scale factor is solved using an integrated Boltzmann approach. This accounts for the ALP energy density contribution to the expansion rate ($H$) and the modification of the photon-to-neutrino temperature ratio ($T_\gamma/T_\nu$) due to ALP decays. The code tracks the ALP momentum distribution using Gauss-Laguerre quadrature to handle cases where ALPs remain in partial equilibrium at MeV temperatures.
4. BBN Network: A custom Mathematica code, BBNEasyALP, was developed to simulate light element synthesis. This code incorporates:
   • Modified time-temperature relations and scale factors.
   • Modified weak interaction rates ($p \leftrightarrow n$) based on the altered neutrino distribution.
   • Meson-driven processes: Crucially, the code includes the injection of metastable mesons ($\pi^\pm, K^\pm, K_L$) from ALP decays. These mesons induce rapid $p \leftrightarrow n$ conversions (e.g., $\pi^- + p \to n + \dots$) and nuclear dissociation (e.g., $\pi^- + {}^4\text{He} \to t + n$), significantly altering the neutron-to-proton ratio and light element yields.
5. Observational Constraints: The results are compared against:
   • Primordial Helium-4 ($Y_P = 0.2458 \pm 0.0013$) and Deuterium ($10^5 \text{D/H}|_P = 2.547 \pm 0.029$) abundances.
   • Effective number of neutrino species ($N_{\text{eff}} = 2.81 \pm 0.12$) derived from Planck, ACT, and SPT data.

Key Contributions and Results

• Impact of Rare Hadronic Decays: The study demonstrates that even tiny hadronic branching ratios ($\text{Br}(a \to \text{hadrons}) \sim 10^{-5}$) can drastically strengthen cosmological bounds. For $m_a \gtrsim 400$ MeV, the injection of charged pions drives the neutron-to-proton ratio above Standard Big-Bang Nucleosynthesis (SBBN) values, increasing the primordial helium yield ($Y_P$). This leads to excluded regions in parameter space that were previously considered viable, particularly for moderate reheating temperatures ($T_{\text{reh}} \sim 1$ TeV).
• Reheating Temperature Dependency: The authors map constraints across a wide range of $T_{\text{reh}}$ (from 5 MeV to $10^{15}$ GeV). They identify a new "island" of excluded parameter space for $10 \text{ GeV} \lesssim T_{\text{reh}} \lesssim 1 \text{ TeV}$ and $m_a \in [0.5, 10]$ GeV, driven primarily by $Y_P$ bounds enhanced by meson interactions.
• Model-Independent Limits: The results are presented in a model-independent framework (mass vs. lifetime and mass vs. coupling), allowing them to be recast for other photophilic or electrophilic relics.
• Tension Resolution: The paper identifies a specific region of parameter space ($m_a \sim 500$ MeV, $\tau_a \sim 300\text{--}5000$ s) where ALPs could simultaneously alleviate two minor cosmological tensions:
  1. The CMB measurement of $N_{\text{eff}}$ being slightly lower than the Standard Model prediction ($N_{\text{eff}} \approx 3.04$).
  2. The observed deuterium abundance being higher than some theoretical SBBN predictions (specifically those using rates from Refs. [57, 71]).
     In this region, ALP decays reduce $N_{\text{eff}}$ to $2.81 \pm 0.12$ and increase the deuterium yield by 2–6%, bringing it into agreement with observations.

Significance and Claims
The authors claim that this work provides a "fresh look" by simultaneously addressing three limitations of previous studies: the use of outdated observational data, the assumption of infinite reheating temperatures, and the neglect of rare hadronic decay channels. By including the meson-driven nuclear reaction network, they establish that cosmological bounds on photophilic ALPs are significantly stronger and more complex than previously thought, particularly for masses above the pion threshold.

The paper emphasizes that while the identified parameter space for resolving cosmological tensions is not statistically definitive, it highlights a viable region where ALPs could play a role in current cosmological anomalies. The authors provide their BBNEasyALP code publicly to facilitate further model testing and generalizations to longer lifetimes or different coupling patterns. They maintain a modest tone regarding the tension resolution, noting that the hints are not statistically significant and could disappear with improved nuclear cross-section data or theoretical predictions.