Lower Limits on the Strengths of Gamma Ray Lines from WIMP Dark Matter Annihilation

We study the spectra of gamma ray signals that arise from dark matter annihilation in the universe. We focus on the large class of theories where the photon spectrum includes both continuum spectrum of gamma rays that arise from annihilation into Standard Model states at tree level, as well as monochromatic gamma rays arising from annihilation directly into two photons at the one loop level. In this class of theories we obtain lower bounds on the ratio of the strength of the gamma ray line relative to the gamma ray continuum as a function of the dark matter mass and spin. These limits are obtained from the unitarity relation between the tree level amplitude of the primary annihilation channel and the imaginary part of the loop level amplitude for annihilation directly into photons, with the primary decay products running in the loop. These results are exact in the limit that dark matter annihilation is exclusively to a single Standard Model species, occurs through the lowest partial wave and respects CP. Away from this limit the bounds are approximate. Our conclusions agree with the known results in the literature in the case of the Minimal Supersymmetric Standard Model (MSSM). We use the Fermi-LAT observations to translate these limits into upper bounds on the dark matter annihilation cross section into any specific Standard Model state.

💡 Research Summary

The paper investigates the gamma‑ray signatures produced by annihilation of weakly interacting massive particle (WIMP) dark matter, focusing on the relationship between the monochromatic line signal (γγ) generated at one‑loop order and the broad continuum spectrum arising from tree‑level annihilation into Standard Model (SM) final states. By exploiting the unitarity condition of the S‑matrix, S†S = 1, and assuming CP invariance together with dominance of the s‑wave (L = 0) partial wave, the authors derive a model‑independent lower bound on the ratio of the line cross‑section σ(χχ → γγ) to the continuum cross‑section σ(χχ → X X̄), where X denotes a single SM species.

The derivation proceeds by writing S = 1 + iT and using −i(T−T†)=T†T. Taking matrix elements between the initial dark‑matter pair |ii⟩ and the two‑photon final state |ff⟩, and invoking time‑reversal (or CP) symmetry, one finds that the imaginary part of the loop amplitude for χχ → γγ is directly related to the product of the tree‑level amplitude for χχ → X X̄ and the one‑loop amplitude for X X̄ → γγ. Squaring this relation yields

 4 |Im⟨ff|T|ii⟩|² = |⟨ff|T†|X⟩|² |⟨X|T|ii⟩|².

Since the left‑hand side is proportional to σ(χχ → γγ) and the factor |⟨X|T|ii⟩|² to σ(χχ → X X̄), the ratio of the two cross‑sections is bounded from below by the purely SM quantity |⟨ff|T|X⟩|²/4, i.e. the absolute square of the imaginary part of the X X̄ → γγ loop amplitude. This quantity can be computed exactly within the SM, independent of the details of the dark‑matter coupling to X.

To make the bound concrete, the authors introduce a fictitious boson Φ with mass m_Φ = 2 m_χ that decays at tree level exclusively to X X̄ and at one loop to γγ. The ratio Γ_im(Φ → γγ)/Γ(Φ → X X̄) is identical to the expression above, establishing

 σ(χχ → γγ)/σ(χχ → X X̄) ≥ Γ_im(Φ → γγ)/Γ(Φ → X X̄).

The right‑hand side depends only on SM parameters (charges, masses, spins) and can be evaluated for any chosen final state X.

The paper then analyses how the bound applies to different dark‑matter spins and CP properties. For scalar dark matter (real or complex) the initial state is CP‑even with J = 0; annihilation to light fermions is helicity‑suppressed, so the dominant channels are the heaviest kinematically accessible fermion or W⁺W⁻. In the fermionic case the final state must have L = 1, S = 1, allowing the bound. For Majorana fermion dark matter the initial state is CP‑odd, J = 0, and the allowed final state is L = 0, S = 0, again yielding a bound for fermion or W⁺W⁻ channels. Dirac fermion dark matter can annihilate with J = 0 or J = 1; the J = 1 channel is forbidden to γγ by the Landau‑Yang theorem, so the weakest bound comes from the J = 0 component. For real vector dark matter the initial state can have J = 0 or J = 2; the J = 0 case mirrors the scalar analysis, while J = 2 involves many (L, S) combinations, preventing a universal bound except in special kinematic limits (e.g., near the W‑mass threshold).

Having established the theoretical framework, the authors apply it to real data. Using the Fermi‑LAT measurements of both line searches and diffuse gamma‑ray fluxes, they compute Γ_im/Γ for several SM final states (b b̄, t t̄, W⁺W⁻, etc.) and translate the line‑search limits into upper limits on σ(χχ → X X̄). The resulting constraints are consistent with, and in many cases reproduce, the known limits obtained in specific models such as the Minimal Supersymmetric Standard Model (MSSM), Universal Extra Dimensions (UED), and the Littlest Higgs model with T‑parity.

In summary, the work provides a robust, model‑independent lower bound on the strength of gamma‑ray lines relative to the continuum for a wide class of WIMP scenarios. By grounding the bound in SM loop calculations and unitarity, it sidesteps uncertainties associated with unknown new‑physics particles in the loop. The methodology can be directly applied to any dark‑matter model that annihilates dominantly into a single SM final state via s‑wave, CP‑conserving processes, offering a powerful tool for interpreting current and future gamma‑ray observations.