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

Basic Information

• Paper ID: 2510.11791
• Title: BBN Constraints on the Hadronic Annihilation of sub-GeV Dark Matter
• Authors: Afif Omar (University of Victoria & TRIUMF), Adam Ritz (University of Victoria)
• Classification: hep-ph (High Energy Physics - Phenomenology)
• Publication Date: October 14, 2025
• Paper Link: https://arxiv.org/abs/2510.11791

Abstract

This paper investigates the effects of residual annihilation of sub-GeV mass thermal relic dark matter candidates during Big Bang Nucleosynthesis (BBN). The focus is on candidates with p-wave annihilation channels, demonstrating hadronic injection of pions and kaons after freeze-out and their interactions with protons and neutrons before the deuterium bottleneck, providing annihilation sensitivity beyond Cosmic Microwave Background (CMB) and galactic indirect detection.

Research Background and Motivation

Problem Background

1. Challenges of sub-GeV Dark Matter: Traditional thermal relic dark matter models focus on the electroweak scale, but dark matter in the sub-GeV mass range faces severe constraints, particularly from CMB observations, requiring thermal relic models with mass below ~5 GeV to suppress annihilation at low velocities.
2. Limitations of Existing Constraints:
   • CMB constraints exclude thermal relic models below ~5 GeV (unless low-velocity annihilation is suppressed)
   • Indirect detection provides weak constraints in the sub-GeV region due to the "MeV gap"
   • The potential of BBN as a detection tool remains underdeveloped
3. Research Motivation:
   • sub-GeV dark matter models typically require dominant p-wave annihilation channels to avoid CMB constraints
   • BBN may be the region with the highest characteristic temperature, providing optimal sensitivity to p-wave annihilation
   • Need to reassess BBN's constraint capability on sub-GeV dark matter

Core Contributions

1. Established a new BBN constraint framework: Developed a BBN constraint method based on hadronic injection for p-wave annihilation dark matter in the sub-GeV mass range (100-1000 MeV).
2. Discovered key effects of hadronic injection: Demonstrated that the long lifetimes of charged pions and kaons enable them to produce significant effects during the BBN phase before the deuterium bottleneck.
3. Sensitivity beyond existing detection methods: Showed that BBN constraints exceed CMB and galactic indirect detection sensitivity in certain parameter regions.
4. Provided benchmark model analysis: Conducted detailed analysis of benchmark models with scalar dark matter and dark photon mediators, demonstrating complementarity with ground-based experiments.

Methodology Details

Task Definition

Investigate the effects of charged mesons (π±, K±) produced through p-wave annihilation of sub-GeV thermal relic dark matter during BBN on the abundances of light elements, particularly corrections to deuterium (D) and helium-4 (⁴He) abundances.

Theoretical Framework

1. p-wave Annihilation Model

For p-wave dark matter annihilation, the annihilation cross section is:

⟨σv⟩_{p-wave} = (3/2) × (T/m_χ) × b

where b characterizes the magnitude of the p-wave annihilation rate, T is temperature, and m_χ is the dark matter mass.

2. Boltzmann Equation System

The Boltzmann equation for meson evolution takes the form:

dn_π/dt + 3Hn_π = DM injection - meson decay + meson/nucleon interactions + ...

Specifically including three main terms:

• Dark Matter Injection: Γ_χ × X²_χ/2, where X_χ is the dark matter number density ratio relative to baryons
• Meson Decay: X_π±/τ_π±, τ_π± = 2.6×10⁻⁸ s
• Charge Exchange Reactions: Modifying the proton-neutron ratio through interactions with nucleons

3. Key Interaction Processes

Pion Charge Exchange Reactions:

π⁻ + p → n + γ,     (σv) ≃ 0.57 mb
π⁻ + p → n + π⁰,    (σv) ≃ 0.88 mb  
π⁺ + n → p + π⁰,    (σv) ≃ 1.7 mb

Kaon Charge Exchange Reactions:

K⁻ + p → n + X:     (σv) ≃ 32 mb
K⁻ + n → p + X:     (σv) ≃ 13 mb

Numerical Calculation Methods

BBN Network Solution

• Employed an adapted BBN numerical framework including neutrons, protons, deuterium, tritium, ³He, and ⁴He
• Used standard BBN assumptions: radiation-dominated ΛCDM cosmology, three neutrino species, instantaneous decoupling
• Key parameters: baryon-to-photon ratio η_b ≃ 6.1×10⁻¹⁰, neutron lifetime τ_n ≃ 878.4 s

Coupled Evolution System

Coupled the meson Boltzmann equations with the BBN reaction network:

Ẋ_π⁻ = Γ_χπ × X²_χ/2 - X_π⁻/τ_π± - Ẋ_pπ⁻
Ẋ_π⁺ = Γ_χπ × X²_χ/2 - X_π⁺/τ_π± - Ẋ_nπ⁺

Experimental Setup

Parameter Space

• Dark Matter Mass Range: m_π ≤ m_χ ≤ 10 GeV (limited by kinematic thresholds)
• Annihilation Strength Parameter: b varying from 10⁻²⁵ to 10⁻¹⁵ cm³/s
• Branching Ratios: Considered both BR = 1 and BR = 0.1 cases

Observational Constraints

Used precise light element abundance observational data:

• Deuterium Abundance: D/H × 10⁻⁵ = 2.547 ± 0.025
• Helium-4 Mass Fraction: Y_p = 0.245 ± 0.003

Benchmark Model Parameters

For scalar dark matter models with dark photon mediators:

• Dark coupling constant: α_D = 0.5
• Mass ratio: m_A'/m_χ = 3
• Kinetic mixing parameter: ε (parameter to be constrained)

Experimental Results

Main Constraint Results

1. Universal p-wave Annihilation Constraints

• Pion Channel: At m_χ ~ 400 MeV, 95% CL constraint is b < 10⁻²³ cm³/s
• Kaon Channel: Due to larger charge exchange cross sections, constraints are stronger, approximately b < 5×10⁻²⁴ cm³/s
• Constraint strength varies with mass, strongest at intermediate masses (~500 MeV)

2. Comparison with Other Detection Methods

• Exceeds CMB Constraints: In specific parameter regions, BBN constraints are more stringent than CMB
• Exceeds Indirect Detection: Due to the MeV gap, BBN provides better sensitivity in the sub-GeV region
• Complementary to Cosmic Ray Constraints: Comparable to Voyager 1 and AMS-02 constraints in certain mass ranges

3. Benchmark Model Constraints

For dark photon mediator models at m_χ ~ 400 MeV:

• 90% CL constraint: ε < 10⁻⁷
• Forms complementary constraints with ground-based experiments (BaBar, NA64, CRESST-II)
• Partial parameter space approaches the thermal relic line

Physical Mechanism Analysis

Timescales of Meson Evolution

From Figure 3, the typical evolution process shows:

1. t ~ 10⁻³ s: Mesons decouple from dark matter chemical equilibrium
2. t ~ 0.1 s: Charge exchange interactions dominate
3. t > 0.1 s: Meson decay becomes the dominant process

Dynamic Changes in Nucleon Ratios

• Coulomb enhancement of π⁻ favors proton → neutron conversion
• Increases available neutron abundance before the bottleneck
• Subsequently rebalances through weak interactions, affecting final light element abundances

Related Work

Historical Development of BBN Constraints

1. Long-lived Particle Decay: BBN as a probe of new physics has been widely studied 5-12
2. Annihilation Process Constraints: Relatively understudied, typically considered less sensitive than other observations 13-16
3. MeV Dark Matter: Recent BBN constraints on MeV-scale dark matter estimated at 3-10 MeV 17-19

sub-GeV Dark Matter Research

• Lee-Weinberg bound requires new force mediators
• p-wave suppression or mass threshold mechanisms avoid CMB constraints
• Multiple mediator models including dark photons and dark scalars

Conclusions and Discussion

Main Conclusions

1. BBN Provides Competitive Constraints: For sub-GeV p-wave annihilation dark matter, BBN constraints exceed CMB and indirect detection in certain parameter regions.
2. Importance of Hadronic Injection: The long lifetimes and strong interactions of charged mesons produce significant effects during the early BBN phase.
3. Complementary Value: BBN constraints form good complementarity with ground-based experiments, particularly in dark photon mediator models.

Limitations

1. Observational Precision Limits: Current observational precision of D and ⁴He abundances is insufficient to fully constrain benchmark thermal relic models.
2. Theoretical Simplifications: Employed stopped meson approximation, neglecting certain secondary processes.
3. Model Dependence: Results depend on specific annihilation channels and branching ratio assumptions.

Future Directions

1. Enhanced Observational Precision: More precise light element abundance measurements will strengthen constraint capabilities.
2. Extended Models: Consider other models such as pseudo-Dirac dark matter and inelastic dark matter.
3. New Mediators: Explore other mediator mechanisms beyond dark photons.

In-Depth Evaluation

Strengths

1. Innovative Approach: First systematic connection between p-wave annihilation sub-GeV dark matter and BBN hadronic injection.
2. Rigorous Calculations: Established a complete Boltzmann equation system, correctly handling coupled evolution of multiple systems.
3. Practical Value: Provides new constraint tools for sub-GeV dark matter models, filling detection gaps.
4. Clear Physical Picture: Effectively explains the physical mechanisms of meson injection affecting BBN.

Weaknesses

1. Limited Parameter Space: Primarily focuses on specific mass ranges, with insufficient coverage of lower mass regions.
2. Branching Ratio Treatment: Relatively simplified treatment of branching ratios for complex annihilation channels.
3. Systematic Uncertainties: Insufficient discussion of BBN theoretical uncertainties.

Impact

1. Academic Value: Provides new theoretical tools for dark matter phenomenology research.
2. Experimental Guidance: Offers theoretical predictions for future precision BBN observations and dark matter search experiments.
3. Methodological Contribution: The established computational framework can be generalized to other new physics scenarios.

Applicable Scenarios

• Theoretical constraints on sub-GeV thermal relic dark matter models
• Phenomenological research on dark sector models
• Cross-disciplinary research between precision cosmology and particle physics
• Development of detection strategies for new physics in the early universe

References

The paper cites abundant relevant literature, including key works in BBN theoretical development 5-12, important research on dark matter constraints 13-19, and theoretical foundations of sub-GeV dark matter models 40-48, reflecting deep understanding of the field's development trajectory.