BBN Constraints on the Hadronic Annihilation of sub-GeV Dark Matter

Problem Statement
Thermal relic dark matter (DM) models in the sub-GeV mass range face stringent constraints from Cosmic Microwave Background (CMB) observations, which typically rule out masses below $\sim 5$ GeV unless the annihilation rate is suppressed at low velocities (e.g., via p-wave channels). While indirect detection in the Milky Way and satellite galaxies provides constraints, the sub-GeV regime is subject to weaker limits due to the "MeV gap" in observational sensitivity. Consequently, Big Bang Nucleosynthesis (BBN) represents a potentially critical, yet under-explored, regime for testing these models. The specific problem addressed is whether residual hadronic annihilation from sub-GeV thermal relics (specifically those with dominant p-wave channels) occurring after freeze-out but before the deuterium bottleneck can significantly alter light element abundances ($D$ and $^4$He) through the injection of charged mesons.

Methodology
The authors investigate the impact of residual DM annihilation into charged pions ($\pi^\pm$) and kaons ($K^\pm$) on the BBN nuclear reaction network. The analysis employs the following steps:

1. Framework Adaptation: The study adapts the numerical framework from Ref. [23], solving a coupled network of Boltzmann equations for neutrons, protons, deuterium, tritium, $^3$He, and $^4$He.
2. Injection Modeling: The authors model DM freeze-out and subsequent residual annihilation ($\chi\chi \to \pi^+\pi^-$ or $K^+K^-$) assuming a fixed relic density ($\Omega_{DM}h^2 \simeq 0.12$). For p-wave annihilation, the rate is parameterized as $\langle \sigma v \rangle \propto T/m_\chi$.
3. Meson Evolution: The evolution of meson number densities is tracked via Boltzmann equations that include:
   • Residual injection from DM annihilation.
   • Meson decay (accounting for lifetimes of $\sim 10^{-8}$ s).
   • Charge conversion interactions with nucleons (e.g., $\pi^- + p \to n + \gamma$), which are enhanced by Coulomb factors.
   • Secondary production (e.g., kaon decay into pions).
4. Assumptions: The analysis assumes standard BBN (SBBN) conditions, including a $\Lambda$CDM expansion history, instantaneous neutrino decoupling (approximated by $N_\nu = 3.045$), and that mesons are efficiently stopped by Coulomb interactions before interacting with nuclei, thereby neglecting high-energy hadro-dissociation or significant final-state radiation effects.
5. Benchmarking: The results are compared against observational constraints on deuterium abundance ($D/H$) and helium mass fraction ($Y_p$). The authors also apply these constraints to a specific benchmark model: scalar DM annihilating via a dark photon mediator.

Key Contributions and Results

• Mechanism of Impact: The paper demonstrates that charged pions and kaons injected after DM freeze-out (but before the deuterium bottleneck, $T \sim \text{MeV}$) induce charge conversion reactions ($\pi^- + p \to n$ and $\pi^+ + n \to p$). Due to Coulomb enhancement, these reactions bias the neutron-to-proton ratio, increasing the neutron fraction relative to SBBN predictions before the bottleneck.
• Sensitivity Comparison: The authors find that BBN provides sensitivity to p-wave annihilation that surpasses current CMB constraints and indirect detection limits in the galaxy for the sub-GeV mass range ($\sim 100\text{--}1000$ MeV).
• Quantitative Limits:
  • For generic p-wave annihilation, the study derives 95% CL limits on the interaction strength parameter $b$ (related to $\langle \sigma v \rangle$) for DM masses between the kinematic thresholds of pions ($m_\pi$) and kaons ($m_K$) up to 10 GeV.
  • Kaons are found to have a stronger impact than pions due to larger charge-exchange cross-sections, though they are only kinematically accessible for $m_\chi \ge m_K$.
  • In the benchmark scalar DM model with dark photon mediation, the BBN constraints on the kinetic mixing parameter $\epsilon$ are shown to be competitive with, and complementary to, terrestrial probes like BaBar, NA64, and CRESST-II.
• Kinematic Thresholds: The analysis confirms that for $m_\chi \gtrsim m_\pi$, freeze-out occurs well before the deuterium bottleneck, allowing the injected mesons to influence the pre-bottleneck era significantly.

Significance and Claims
The paper claims that BBN serves as a robust and complementary tool for constraining sub-GeV thermal relic dark matter, particularly for models requiring p-wave suppression to evade CMB bounds. While the current observational precision on light element abundances is not yet sufficient to fully exclude all benchmark thermal relic models with p-wave annihilation, the sensitivity is superior to other indirect channels. The authors emphasize that this approach offers a unique window into the sub-GeV regime, where the combination of late freeze-out and the long lifetimes of charged mesons allows for significant modifications to the neutron-proton conversion rates. The study suggests that future improvements in the precision of light element abundance measurements could further enhance these constraints, potentially ruling out or confirming specific classes of sub-GeV dark matter candidates. Additionally, the authors note that BBN can constrain the mass splitting in pseudo-Dirac dark matter models if s-wave annihilation remains sizable.