The impact of gamma-ray propagation effects on indirect dark matter searches

This paper demonstrates that neglecting the regeneration of high-energy gamma rays via inverse Compton scattering by secondary electrons and positrons during propagation can lead to significant underestimations of the observed flux from distant dark matter annihilation, thereby necessitating a more detailed treatment of these effects to derive accurate indirect detection limits.

The Gist

The Big Picture: Hunting Ghosts with Flashlights

Imagine the universe is filled with invisible "ghosts" called Dark Matter. Scientists have never seen them directly, but they know they are there because of how they tug on stars and galaxies. One popular theory is that these ghosts are actually heavy particles called WIMPs (Weakly Interacting Massive Particles).

The big idea is: if two WIMPs crash into each other, they might annihilate (disappear) and explode into a shower of other particles, including gamma rays (super-high-energy light). If we can detect these gamma rays coming from deep space, we might finally catch a glimpse of Dark Matter.

The Problem: The "Fog" of the Universe

For a long time, scientists have been playing a game of "catch the flashlight beam." They look at distant galaxies or galaxy clusters, expecting to see a specific pattern of gamma rays.

However, the universe isn't empty; it's filled with a faint, invisible "fog" made of background light (like the afterglow of the Big Bang and starlight from other galaxies).

The Old Way of Thinking:
Scientists used to think that when a high-energy gamma ray travels through this fog, it hits a background photon and simply dies. It turns into an electron and a positron (matter and antimatter), and the gamma ray is gone forever. They calculated the signal based on this "attenuation" (loss of light).

The New Discovery:
This paper says, "Wait a minute! We forgot what happens after the gamma ray dies."

The Analogy: The Broken Vase and the Bouncing Ball

Let's use an analogy to understand the physics:

1. The Original Signal: Imagine a WIMP annihilation is like a giant vase being smashed. It shatters into thousands of sharp shards (high-energy gamma rays) flying toward Earth.
2. The Old View (Attenuation): As these shards fly through the "fog" of the universe, they hit invisible walls (background light) and shatter further, turning into dust (electrons). The old models assumed the dust just fell to the ground and disappeared. The signal got weaker and weaker.
3. The New View (Inverse Compton Scattering): The authors realized that when the vase shatters into dust (electrons), those dust particles don't just sit there. They are moving fast! As they zoom through the fog, they kick the invisible fog particles.
   • Think of it like a fast-moving billiard ball (the electron) hitting a stationary ping-pong ball (a background photon).
   • The kick sends the ping-pong ball flying at high speed, turning it into a new high-energy gamma ray.

The Result: The original gamma ray died, but it created a new gamma ray. The energy didn't disappear; it just got repackaged.

What the Paper Found

The authors used powerful computer simulations (like a video game engine for the universe) to track these particles from the source to Earth. They found three major things:

1. The Signal Changes Shape: Because the "repackaging" process turns high-energy gamma rays into slightly lower-energy ones, the spectrum (the "color" of the light) changes. The high-energy end gets dimmer, but the lower-energy end gets much brighter.
2. Distance Matters: The farther away the source is, the more times this "kick" happens. For very heavy WIMPs (the "super-heavy" ghosts) and very distant sources (like the Perseus galaxy cluster), the difference is massive.
   • The Math: In some cases, the amount of light we see at lower energies can be 1,000 times (3 orders of magnitude) higher than scientists previously thought because they forgot about the "kicking" electrons.
3. The "Leptonic" vs. "Hadronic" Difference:
   • If the WIMPs break apart into heavy particles (like bottom quarks), the effect is moderate.
   • If they break apart into lighter particles (like tau leptons), the "kicking" is much more violent, creating a huge surge of new gamma rays.

Why This Changes Everything

This isn't just a small correction; it changes the rules of the game for hunting Dark Matter.

• The "Exclusion Limits": Scientists have been drawing lines on a graph saying, "If Dark Matter existed with these properties, we would have seen it by now. Since we didn't, it doesn't exist."
• The Shift: Because the new model predicts more gamma rays at certain energies (due to the regeneration effect), the "lines" on the graph need to move.
  • Some areas that scientists thought were "ruled out" (proven empty) might actually still be valid hunting grounds.
  • Conversely, if we do see a signal, we might have been misinterpreting it because we didn't account for this regeneration.

The Takeaway

Think of this paper as realizing that when you shout across a canyon, the sound doesn't just fade away. The echo bounces off the walls, hits the ground, and creates new sounds that reach your ear.

If you are trying to figure out how loud the original shout was, you have to account for all those echoes. If you ignore them, you'll get the wrong answer.

In short: The authors are telling the scientific community, "We've been ignoring the echoes (secondary gamma rays) for too long. If we want to find Dark Matter, we need to listen to the whole symphony, not just the first note." This is crucial for upcoming telescopes like the CTAO and Fermi-LAT, which are about to start listening very closely.

Technical Summary

1. Problem Statement

Indirect dark matter (DM) searches rely on detecting gamma rays produced by the annihilation of Weakly Interacting Massive Particles (WIMPs) in high-density astrophysical regions. Current theoretical models used to interpret observational data from instruments like Fermi-LAT and Imaging Atmospheric Cherenkov Telescopes (IACTs) typically account for the attenuation of high-energy gamma rays due to pair production ($\gamma + \gamma_{bg} \to e^+ + e^-$) on the Extragalactic Background Light (EBL) and Cosmic Microwave Background (CMB).

However, these standard models often neglect the fate of the resulting secondary electrons and positrons. Specifically, they fail to account for Inverse Compton Scattering (ICS), where these secondary leptons upscatter background photons (CMB/EBL) back into high-energy gamma rays. This omission ignores the regeneration of the gamma-ray flux via electromagnetic cascades, potentially leading to significant inaccuracies in predicted spectra and derived exclusion limits, particularly for distant sources and heavy WIMP masses.

2. Methodology

The authors developed a comprehensive framework combining two distinct simulation tools to model both the source emission and the propagation effects:

• Source Generation (CosmiXs):

  • Used to compute the intrinsic "prompt" gamma-ray spectra ($dN_\gamma/dE_\gamma$) resulting from WIMP annihilation.
  • Simulated annihilation channels: $b\bar{b}$ (hadronic) and $\tau^+\tau^-$ (leptonic).
  • WIMP masses ($m_\chi$) ranged from 1 TeV to 100 TeV.
  • The code incorporates Quantum Electrodynamics (QED) and Quantum Chromodynamics (QCD) effects to generate accurate secondary particle spectra.

• Propagation Modeling (CRPropa 3):

  • Used to simulate the propagation of gamma rays and secondary particles through the intergalactic medium (IGM).
  • Physics included: Pair production on CMB, EBL, and Cosmic Radio Background (CRB); Inverse Compton Scattering (ICS) of secondary $e^\pm$ on the same backgrounds.
  • Physics neglected: Higher-order interactions (double/triplet pair production) were excluded as their cross-sections are negligible in the TeV energy range. Magnetic field deflections were also neglected (1D propagation mode).
  • Simulation parameters:
    • Distances: 0.5 Mpc, 5 Mpc, 100 Mpc, and 600 Mpc.
    • Specific Case Study: The Perseus Galaxy Cluster ($z \approx 0.017$, $\sim 70$ Mpc).
    • Initial energies: Up to 120 TeV.
    • Minimum energy threshold: 1 GeV.

• Analysis Strategy:

  • The intrinsic spectra from CosmiXs were used to reweight CRPropa propagation simulations.
  • The authors compared fluxes calculated with ICS ($F_{ICS}$) versus without ICS ($F_{nICS}$).
  • They evaluated the impact on the velocity-averaged annihilation cross-section ($\langle \sigma v \rangle$) limits, specifically referencing the projected sensitivity of the Cherenkov Telescope Array Observatory (CTAO).

3. Key Contributions

• Integrated Framework: Established a consistent pipeline linking intrinsic DM annihilation spectra (CosmiXs) with full electromagnetic cascade propagation (CRPropa), explicitly including the regeneration of gamma rays via ICS.
• Quantification of Regeneration: Demonstrated that neglecting ICS leads to a systematic underestimation of the gamma-ray flux at lower energies (GeV–TeV) for distant sources, as the energy lost by primary gamma rays is partially recovered by secondary photons.
• Channel Dependence: Highlighted that the effect is significantly more pronounced for leptonic channels ($\tau^+\tau^-$) compared to hadronic channels ($b\bar{b}$), due to the higher production rate of high-energy electrons/positrons in leptonic annihilation.

4. Key Results

A. Spectral Modification

• Flux Redistribution: Including ICS results in a "softening" of the spectrum. High-energy photons are attenuated, but their energy is transferred to secondary electrons, which then upscatter background photons to produce a regenerated flux at lower energies.
• Magnitude of Effect:
  • For distant sources (e.g., 600 Mpc) and heavy WIMPs ($m_\chi \sim 100$ TeV), the regenerated flux can exceed the non-ICS flux by three orders of magnitude ($10^3$) in the lower-energy regime (below 1 TeV).
  • For nearby sources (e.g., 0.5 Mpc) or lighter WIMPs (GeV scale), the effect is negligible, and the observed spectrum closely matches the intrinsic source spectrum.
• Channel Comparison: The $\tau^+\tau^-$ channel shows much larger deviations (flux ratios up to $10^3$) compared to $b\bar{b}$ (deviations typically up to one order of magnitude) because leptonic channels produce more energetic $e^\pm$ pairs capable of efficient ICS.

B. Impact on Exclusion Limits

• Perseus Cluster Case: When applying these propagation effects to the Perseus cluster:
  • The derived limits on $\langle \sigma v \rangle$ are modified.
  • For the $\tau^+\tau^-$ channel and high masses ($m_\chi \sim 100$ TeV), the inclusion of ICS shifts the exclusion limits by a factor of $\sim 0.4$ (i.e., the constraints become less stringent).
  • This implies that previous limits derived without ICS may be overly conservative for specific mass ranges and channels.
• Energy Window Dependence:
  • Limits derived from Fermi-LAT (sensitive to $\sim 1-100$ GeV) are highly sensitive to these propagation effects for high-mass WIMPs.
  • Limits from IACTs (sensitive to $\gtrsim 1$ TeV) show minimal dependence on ICS for the Perseus cluster, as the regenerated flux peaks at lower energies.

5. Significance and Implications

• Re-evaluation of DM Constraints: The study suggests that current indirect detection limits, particularly those relying on distant targets or heavy WIMP candidates, may need revision. Ignoring ICS can lead to systematic biases, potentially ruling out WIMP parameter space that is actually viable.
• Necessity for Future Observatories: For next-generation instruments like CTAO, which will probe the multi-TeV to PeV range with high sensitivity, accurate modeling of electromagnetic cascades is not just a refinement but a necessity for reliable signal prediction.
• Superheavy DM: The findings are critical for searches involving superheavy DM candidates ($m_\chi > 100$ TeV), where electromagnetic cascades dominate the observable signal.
• Methodological Shift: The paper advocates for a shift in standard analysis pipelines to include full propagation treatments (pair production + ICS) rather than simple attenuation models, ensuring that the "fate" of secondary particles is accounted for in the interpretation of gamma-ray data.

In conclusion, the authors demonstrate that the regeneration of gamma rays via inverse Compton scattering is a dominant effect for distant, heavy WIMP sources, fundamentally altering the expected spectral shape and weakening the derived exclusion limits for specific annihilation channels.