Pith. sign in

REVIEW 4 major objections 6 minor 39 references

A confirmed 20 GeV high-latitude Galactic halo excess matches sub-TeV WIMP annihilation once systematics are folded in.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-10 05:29 UTC pith:42HK6QGU

load-bearing objection Independent pixel-level confirmation of Totani's 20 GeV high-latitude residual is solid; the WIMP reading is carefully framed as consistency and remains load-bearing on diffuse-model choice. the 4 major comments →

arxiv 2607.08552 v1 pith:42HK6QGU submitted 2026-07-09 astro-ph.HE hep-ex

The 20 GeV Galactic Halo Excess: Pixel-Level Confirmation and Consistency with Sub-TeV WIMP Annihilation

Trinity Rosebud Stenhouse , Chamkaur Ghag , Frank F. Deppisch This is my paper
classification astro-ph.HE hep-ex PACS 95.35.+d98.70.Rz95.85.Pw
keywords Fermi-LATGalactic halo excess20 GeV residualWIMP annihilationNFW profiledwarf spheroidal limitsSommerfeld enhancementdark matter indirect detection
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper independently recovers a spherically symmetric, high-latitude gamma-ray residual in Fermi-LAT data that peaks near 20 GeV, first reported in a cell-aggregated analysis. Extending the fit to the native 0.125-degree pixel maps with energy-dependent PSF folding and bright-source masking leaves the same spectral feature, with a normalisation about 20 percent higher. The residual is centrally concentrated, spectrally and spatially distinct from the few-GeV Galactic-centre excess, and strongly disfavours an isotropic extragalactic origin. Prompt s-wave annihilation spectra fit the feature with best-fit masses of roughly 0.55 TeV (W+W−) and 0.72 TeV (b b-bar) and a cross section near 10−24 cm3 s−1. Face-value tension with dwarf-spheroidal limits is a factor of 4–5, but foreground-modelling and J-factor uncertainties widen the window to 1.6–9.3, leaving the s-wave interpretation viable; the only velocity structure that also matches relic density is low-velocity-enhanced annihilation supplying an approximate 45-fold boost from the thermal value.

Core claim

Independent cellwise and pixel-level template fits of fifteen-plus years of Fermi-LAT data both recover a high-latitude, NFW-like residual peaking near 20 GeV. The feature is morphologically preferred over isotropic emission, survives disk-included versus disk-excluded tests that isolate it from the inner-Galaxy excess, and is consistent with prompt s-wave WIMP annihilation at sub-TeV masses once the full systematic budget on foregrounds and J-factors is included.

What carries the argument

The global pixel-level Poisson likelihood on native 0.125° maps, with energy-dependent PSF forward-folding of the halo and diffuse templates plus bright-source masking; this retains spatial information suppressed by 10° cell aggregation while controlling point-source leakage, and supplies the per-bin halo spectrum that is then fit by PPPC annihilation (or decay) yields and compared to velocity-dependent effective J-factors.

Load-bearing premise

The residual left after subtracting the conventional gas, inverse-Compton, bubble, Loop I and point-source templates is not simply residual mismodelling of the Galactic interstellar emission, especially the choice between GALPROP and the Fermi GIEM diffuse model that alone can halve the recovered flux.

What would settle it

A future diffuse-model suite or independent high-latitude analysis that removes the 20 GeV residual entirely while preserving the fit quality of the conventional templates, or a simultaneous antiproton spectrum that cannot be reconciled even with a light-mediator cascade.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • If the residual is dark matter, the preferred particle is a sub-TeV WIMP whose present-day annihilation rate is enhanced by a light mediator relative to freeze-out.
  • Standard thermal s-wave annihilation remains allowed at the lower edge of the systematic band; pure p-wave or pure decay are ruled out by relic density or the isotropic gamma-ray background.
  • The high-latitude morphology supplies a cleaner target than the Galactic centre for future gamma-ray and multi-messenger tests of the same mass window.
  • Fully eliminating the residual dwarf tension would require a resonance whose enhancement peaks near Galactic-halo velocities and declines for colder dwarf systems.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Because the absolute spectral fit remains poor (χ²/ν ~ 8) from coherent diffuse residuals, the mass is better constrained by peak location than the overall normalisation is by amplitude.
  • A light-mediator realisation that softens the antiproton spectrum could simultaneously ease the AMS-02 tension cited in the paper and still produce the observed gamma-ray yield.
  • The same pixel-level pipeline could be re-run on future Pass-8 or Pass-9 data releases with updated diffuse models to test whether the 20 GeV peak migrates or vanishes.
  • If confirmed, the required ~45× boost at freeze-out versus today would tighten model-building around sub-GeV mediators whose saturation velocity sits between dwarf and halo scales.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. The manuscript independently reproduces Totani’s cell-aggregated Fermi-LAT analysis of a ~20 GeV high-latitude Galactic residual and extends it to a pixel-level Poisson likelihood on native 0.125° maps, with energy-dependent PSF forward-folding and bright-source masking. Both pipelines recover a similar halo spectrum (pixel-level normalisation ~20% higher) for NFW emissivity scalings ρ^p (p=1,2,2.5). The residual is argued to be a high-latitude, centrally concentrated feature distinct from the few-GeV Galactic Centre excess, disfavouring isotropic/extragalactic emission. Prompt s-wave annihilation fits give m_χ ≃ 0.55 TeV (W⁺W⁻) and 0.72 TeV (b b-bar) with ⟨σv⟩ ≃ 1×10⁻²⁴ cm³ s⁻¹, in ~4–5× tension with dSph limits; foreground and J-factor systematics are said to widen the tension to R ≃ 1.6–9.3. Velocity-suppressed annihilation and decay are disfavoured (relic density and IGRB, respectively); low-velocity-enhanced annihilation (resonant Sommerfeld/Breit–Wigner) is presented as the only structure reconciling halo rate, dSph limits, and thermal relic density, at the cost of a light-mediator dark sector.

Significance. If the residual is genuine and not residual interstellar-emission mismodelling, this is a high-latitude, spectrally hard counterpart to the long-debated Galactic Centre excess, extracted where source confusion and disk systematics are weaker. An independent pipeline confirmation with pixel-level likelihood, PSF folding, bright-source masks, north–south split, injection-null test, and a documented systematics suite is valuable, as is the public analysis code. The careful separation of the empirical residual from the model-dependent DM claim, and the explicit mapping of s-wave, p-wave, decay, and Sommerfeld/Breit–Wigner options against dSph limits and relic density, are strengths. The result would motivate targeted light-mediator model building and multi-messenger checks (antiprotons, future gamma-ray instruments) if the diffuse-model dependence can be controlled.

major comments (4)
  1. Sec. IIID and Fig. 4: the GALPROP↔GIEM choice roughly halves the recovered halo flux and moves the best-fit mass over ~0.2–0.9 TeV, driving R from ~5 to ~2. This is correctly identified as the leading systematic, but the central claim that the residual is a distinct, DM-compatible component remains load-bearing on it. Please quantify, under the GIEM-only (and other science-band) variants, whether the ring-profile preference for centrally concentrated NFW over isotropic (Fig. 3) and the per-bin ΔlnL preference for ρ² (Fig. 2) survive at comparable strength, and report the corresponding best-fit m_χ and ⟨σv⟩ with the same J_pole convention. Without that, the lower edge of the R ≃ 1.6–9.3 window is hard to interpret as leaving s-wave viable rather than as residual IEM freedom.
  2. Secs. IIIB and IVA: the absolute spectral fit is poor (χ²/ν ~ 8), attributed to spatially coherent diffuse residuals rather than Poisson noise, and the ring-profile fits are also imperfect (best χ²/ν ≃ 5). The mass is said to be well localised by the peak position, but the paper should show explicitly how the χ² landscape in m_χ behaves when the per-bin halo errors are inflated to absorb the coherent residual (or when the science-band envelope of Fig. 4a is used as a systematic covariance). Otherwise the quoted best-fit masses (0.55/0.72 TeV) and the claim of spectral consistency with PPPC prompt yields overstate the constraining power of the data.
  3. Sec. IVB, Eq. (IV.2) and the R ≃ 1.6–9.3 window: the combination of foreground range (×0.4–1.1), halo J (~48% from ρ_⊙), and dwarf J (~54%) needs a transparent recipe. Are the three terms treated as independent and multiplied/added in quadrature on a linear or log scale? Does the lower edge R ≃ 1.6 correspond to a single coherent choice (e.g. GIEM + high ρ_⊙ + high dwarf J) or to an envelope of extremes? A short table or appendix giving R for each science-band variant at fixed J, then after J rescaling, would make the “approaches the dwarf limit” statement falsifiable.
  4. Sec. IVC and the concluding claim that low-velocity-enhanced annihilation is “the only viable velocity structure”: p-wave and decay are cleanly excluded by relic density and IGRB respectively, which is useful. The required ~45× boost is within published Sommerfeld/Breit–Wigner benchmarks, but fully removing the dwarf tension (R ≤ 1) is stated to need a fine-tuned above-threshold resonance peaking at halo velocity. Please quantify that tuning (mediator mass window, α_χ, saturation velocity relative to σ_v,halo and dwarf moments) and confront the AMS-02 antiproton bound cited in the discussion more directly for the light-mediator case, rather than only noting that softened spectra “have been argued” to evade it. As written, the particle-physics resolution risks reading as more unique and less costly than the calculation supports.
minor comments (6)
  1. Fig. 1 and Fig. 5: the low-energy sign change of the halo amplitude is visible but not discussed in the text beyond a brief mention; a sentence on whether negative amplitudes are physical (over-subtraction) or prior artefacts would help non-specialist readers.
  2. Sec. IIC: the decision not to forward-fold energy dispersion is justified by bin width vs resolution; a one-line estimate of the residual bias on the peak energy (or a reference to a Pass 8 IRF study) would close the point cleanly.
  3. Eq. (IV.1) and the J_pole value: both GeV² cm⁻⁵ and M_⊙² kpc⁻⁵ units are given; ensure the conversion factor is stated once so that cross-checks against Totani and standard NFW conventions are immediate.
  4. Fig. 6 caption is dense; separating the definition of R_s, R_p, R_d, R_S into a short table or inline list would improve readability.
  5. Data availability: the GitHub/Zenodo links are a clear strength; please confirm that the exact energy-bin edges, mask definitions, and GALPROP configuration tags used for the baseline are version-pinned in the archive.
  6. Typographical: “Fermi–LAT” hyphenation and “b ¯b” spacing are inconsistent in a few places (abstract vs body); unify.

Circularity Check

0 steps flagged

No significant circularity: free per-bin template amplitudes extract an empirical spectrum that is then compared to external PPPC spectra and published dSph limits.

full rationale

The derivation chain is data-driven and non-circular. Fermi-LAT counts are fit with free normalizations heta_mk for a fixed family of spatial templates (GALPROP/GIEM gas+ICS, isotropic, 4FGL point sources, Loop I, structured bubbles, NFW ho^p) independently per energy bin (Eqs. II.1–II.5, Secs. IIB–IIC). The resulting halo spectrum is an output of that fit, not an input. Prompt annihilation/decay spectra from the external PPPC4DMID library are then least-squares fit with a single free normalization to that extracted spectrum (Eq. IV.1, Sec. IVA), yielding best-fit m_\chi and angle; the same spectrum is compared to independent published dSph limits (Hoof et al., McDaniel et al.) and velocity-dependent effective J-factors (Boddy et al.). Relic-density and Sommerfeld/Breit–Wigner benchmarks are likewise taken from the external literature. No equation reduces a claimed prediction to a fitted input by construction, no uniqueness theorem is imported from overlapping authors, and the sole citation of Totani is for independent replication of a prior empirical claim rather than a load-bearing definitional premise. Systematics (GALPROP vs GIEM, J-factor uncertainties) are propagated openly and widen rather than close the tension window. The analysis is therefore self-contained against external data and models.

Axiom & Free-Parameter Ledger

5 free parameters · 6 axioms · 1 invented entities

The empirical claim rests on standard Fermi analysis assumptions (Pass 8 IRFs, Poisson likelihood, linear template model) plus the completeness of the chosen diffuse/bubble/Loop I/point-source template set. The DM claim further rests on NFW morphology, PPPC prompt yields, published dSph limits and effective J-factors, a Jeans-modelled halo velocity field, and the thermal-relic normalisation. No new particle is invented by the authors; resonant Sommerfeld/Breit–Wigner is an existing model class invoked to reconcile rates. Free parameters are the per-bin template normalisations and the fitted (m_χ, ⟨σv⟩) or lifetime.

free parameters (5)
  • Per-bin template normalisations θ_mk (halo, gas, ICS, isotropic, bubbles, Loop I, point sources)
    All sky-model amplitudes are free (or free with sign constraints) in each of 13 energy bins; the halo spectrum is entirely determined by these fits.
  • Best-fit WIMP mass m_χ
    Scanned and selected by χ² minimum against the extracted halo spectrum; central values 0.55 TeV (W⁺W⁻) and 0.72 TeV (b b-bar).
  • Annihilation cross section ⟨σv⟩ (or decay lifetime τ)
    Single normalisation fit to the halo pole flux at each mass; best-fit ~1×10⁻²⁴ cm³ s⁻¹ (annihilation) or τ~5–7×10²⁶ s (decay).
  • NFW pole J-factor / local density ρ_⊙
    J_pole fixed at 4.0×10²¹ GeV² cm⁻⁵ with ρ_⊙=0.42±0.10 GeV cm⁻³ driving the ~48% halo J uncertainty that enters the tension budget.
  • Emissivity index p ∈ {1, 2, 2.5}
    Discrete morphology choices for the halo template; p=2 is fiducial for annihilation, p=1 for decay, p=2.5 as a cuspier stress test.
axioms (6)
  • domain assumption The true high-latitude emission is adequately spanned by the linear combination of GALPROP (or GIEM) gas+ICS, isotropic, 4FGL point sources, Loop I, structured Fermi bubbles, and an NFW-like halo template.
    Sec. IIA–IIC; residual structure not captured by this set is what is attributed to the halo. The GALPROP↔GIEM swap halves the signal, showing the assumption is load-bearing.
  • domain assumption Dark-matter density follows a smooth NFW profile with the stated pole normalisation; substructure is neglected in the halo J-factor.
    Sec. IIA, IVA; used for both morphology ranking and ⟨σv⟩ conversion.
  • domain assumption Prompt PPPC4DMID photon yields correctly describe the gamma-ray spectrum from WIMP annihilation/decay to SM final states at the fitted masses.
    Sec. IVA, IVD; secondary inverse-Compton and other non-prompt contributions are not refit.
  • domain assumption Published dSph annihilation limits (Hoof et al. 2020; McDaniel et al. 2024) and velocity-dependent effective J-factors (Boddy et al.) apply as external benchmarks, with ~54% stacked J uncertainty.
    Sec. IVB–IVC; define the tension ratio R.
  • domain assumption Thermal freeze-out with ⟨σv⟩_th ≃ 2.2×10⁻²⁶ cm³ s⁻¹ sets the relic density for a single-species WIMP; velocity moments follow a Jeans model of the NFW halo.
    Sec. IVC; used to exclude thermal p-/d-wave and to compute the ~45× boost required for low-velocity enhancement.
  • standard math Poisson likelihood on counts (cellwise or pixelwise) with flat priors on template amplitudes is the correct statistical model; energy dispersion can be neglected relative to bin width.
    Sec. IIB–IIC; standard for Fermi template analyses.
invented entities (1)
  • Resonant low-velocity-enhanced annihilator (Sommerfeld or Breit–Wigner with light mediator) no independent evidence
    purpose: To supply the ~45× boost from thermal freeze-out to the present-day halo rate while remaining compatible with dSph limits and relic density.
    Not invented by this paper; it is an existing model class (Feng et al., Ibe et al., Guo & Wu, and recent arXiv follow-ups). The paper invokes it as the only velocity structure that works, and notes that fully closing R≤1 needs a fine-tuned above-threshold resonance. independent_evidence is false within this paper because no specific mediator mass or coupling is predicted and tested externally here.

pith-pipeline@v1.1.0-grok45 · 24589 in / 4598 out tokens · 45073 ms · 2026-07-10T05:29:17.055202+00:00 · methodology

0 comments
read the original abstract

A recent analysis of 15 years of Fermi-LAT data reported a spherically symmetric, halo-like component of the Galactic diffuse emission that peaks near 20GeV. We independently reproduce this cell-aggregated analysis, then extend it to a pixel-level likelihood on the native $0.125^\circ$ maps, adding energy-dependent point-spread-function forward folding and masking bright sources. Both methods replicate the 20GeV halo spectrum, with the pixel-level normalisation ${\sim}20\%$ above the cellwise fit across NFW emissivity scalings $\rho^p$, $p \in 1,2,2.5$. This 20GeV halo is a high-latitude feature, distinct from the inner-Galaxy excess, and consistent with sub-TeV dark matter (WIMP) annihilation. It is centrally concentrated, strongly disfavouring extragalactic emission. Fitting prompt $s$-wave annihilation spectra, best-fit masses are $m_\chi \simeq 0.55$TeV ($W^+W^-$) and $0.72$TeV ($b\bar{b}$) with $\langle\sigma v\rangle \simeq 1\times10^{-24}~\mathrm{cm^3\,s^{-1}}$, in $\sim\!4$-$5\times$ tension with dwarf spheroidal galaxy limits. However, accounting for foreground modelling and $J$-factor systematic uncertainties widens the tension window to $R\simeq1.6$-$9.3$, leaving the $s$-wave interpretation viable. To close the tension, we consider alternative particle dark matter models. $p$-wave annihilation misses relic abundance constraints by $\sim\!7$ orders of magnitude. A decay interpretation evades dwarf limits but is disfavoured by the isotropic gamma-ray background. The only viable velocity structure consistent with dwarf limits, present-day halo rates, and relic density is low-velocity-enhanced annihilation (resonant Sommerfeld or Breit-Wigner). This supplies the required $\approx\!45\times$ boost from a thermal relic. Fully resolving the dwarf tension requires a fine-tuned resonance peaking at the halo velocity and falling for colder systems.

Figures

Figures reproduced from arXiv: 2607.08552 by Chamkaur Ghag, Frank F. Deppisch, Trinity Rosebud Stenhouse.

Figure 1
Figure 1. Figure 1: FIG. 1. Comparison of the halo spectra recovered in the [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Per-energy-bin log-likelihood improvement from [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Model-independent angular profile of the 20 GeV [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Systematic stability of the [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Spectral fits of prompt PPPC annihilation models [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Velocity-dependent dSph tension comparison for the global [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Spectral fits of prompt PPPC decay models to the [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages · 30 internal anchors

  1. [1]

    Deterministic Global Fit: Initialisation and Systematics As a fast diagnostic stage, we run a deterministic pixel- level maximum-likelihood fit using the L-BFGS-B al- gorithm [20, 21]. In this optimiser only, the template columns are internally rescaled for numerical condition- ing, φmk =θ mkSmk, S mk = P j Pm,jk P j djk ,(II.6) with the amplitudes transf...

    work page

  2. [2]

    Production Pixelwise MCMC The production pixelwise fit is sampled with an affine- invariant ensemble MCMC. The fitted amplitudes are as- signed flat bounds: standard emission templates are con- strained to be non-negative, while the NFW halo coef- ficient and the negativeFermi-Bubble residual template are allowed to take either sign, as in Sec. IIB. The c...

    work page

  3. [3]

    Uncertainties, Residuals, and the Halo Likelihood Ratio For each template amplitude we report the maximum- likelihoodsampleandthecentralposteriorquantilesfrom the MCMC chain. These intervals marginalise over the other fitted template amplitudes within the same energy bin and replace the local Fisher/Laplace errors used in the deterministic diagnostic fit....

    work page doi:10.5281/zenodo.21280724

  4. [4]

    Particle Dark Matter: Evidence, Candidates and Constraints

    G. Bertone, D. Hooper, and J. Silk, Particle dark matter: Evidence, candidates and constraints, Phys. Rep.405, 279 (2005), arXiv:hep-ph/0404175 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2005

  5. [5]

    Precise Relic WIMP Abundance and its Impact on Searches for Dark Matter Annihilation

    G. Steigman, B. Dasgupta, and J. F. Beacom, Pre- cise relic WIMP abundance and its impact on searches for dark matter annihilation, Phys. Rev. D86, 023506 (2012), arXiv:1204.3622 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  6. [6]

    J.L.Feng,DarkMatterCandidatesfromParticlePhysics and Methods of Detection, Annu. Rev. Astron. Astro- phys.48, 495 (2010), arXiv:1003.0904 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  7. [7]

    Possible Evidence For Dark Matter Annihilation In The Inner Milky Way From The Fermi Gamma Ray Space Telescope

    L. Goodenough and D. Hooper, Possible evidence for dark matter annihilation in the inner Milky Way from the Fermi Gamma Ray Space Telescope, arXiv e-prints (2009), arXiv:0910.2998 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2009

  8. [8]

    Background model systematics for the Fermi GeV excess

    F. Calore, I. Cholis, and C. Weniger, Background model systematics for the Fermi GeV excess, JCAP03, 038 (2015), arXiv:1409.0042 [astro-ph.CO]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  9. [9]

    The Fermi Galactic Center GeV Excess and Implications for Dark Matter

    M. Ackermannet al.(Fermi-LAT Collaboration), The Fermi Galactic Center GeV excess and implications for dark matter, Astrophys. J.840, 43 (2017), arXiv:1704.03910 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  10. [10]

    Strong support for the millisecond pulsar origin of the Galactic center GeV excess

    R. Bartels, S. Krishnamurthy, and C. Weniger, Strong support for the millisecond pulsar origin of the Galactic center GeV excess, Phys. Rev. Lett.116, 051102 (2016), arXiv:1506.05104 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2016

  11. [11]

    Totani, 20 GeV halo-like excess of the Galactic diffuse emission and implications for dark matter annihilation, JCAP11, 080 (2025), arXiv:2507.07209 [astro-ph.HE]

    T. Totani, 20 GeV halo-like excess of the Galactic diffuse emission and implications for dark matter annihilation, JCAP11, 080 (2025), arXiv:2507.07209 [astro-ph.HE]

    work page arXiv 2025

  12. [12]

    W. B. Atwoodet al.(Fermi-LAT Collaboration), The Large Area Telescope on the Fermi Gamma-ray Space Telescope mission, Astrophys. J.697, 1071 (2009), arXiv:0902.1089

    work page internal anchor Pith review Pith/arXiv arXiv 2009

  13. [13]

    M. Su, T. R. Slatyer, and D. P. Finkbeiner, Giant gamma-ray bubbles from Fermi-LAT: Active galactic nu- cleus activity or bipolar galactic wind?, Astrophys. J. 724, 1044 (2010), arXiv:1005.5480

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  14. [14]

    J. F. Navarro, C. S. Frenk, and S. D. M. White, The structure of cold dark matter halos, Astrophys. J.462, 563 (1996), arXiv:astro-ph/9508025

    work page internal anchor Pith review Pith/arXiv arXiv 1996

  15. [15]

    Searching for Dark Matter Annihilation from Milky Way Dwarf Spheroidal Galaxies with Six Years of Fermi-LAT Data

    M. Ackermannet al.(Fermi-LAT Collaboration), Search- ing for dark matter annihilation from Milky Way dwarf spheroidal galaxies with six years of Fermi Large Area Telescope data, Phys. Rev. Lett.115, 231301 (2015), arXiv:1503.02641 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  16. [16]

    S. Hoof, A. Geringer-Sameth, and R. Trotta, A global analysis of dark matter signals from 27 dwarf spheroidal galaxies using 11 years of Fermi-LAT observations, JCAP 02, 012 (2020), arXiv:1812.06986 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2020

  17. [17]

    Legacy Analysis of Dark Matter Annihilation from the Milky Way Dwarf Spheroidal Galaxies with 14 Years of Fermi-LAT Data

    A. McDaniel, M. Ajello, C. M. Karwin, M. Di Mauro, A. Drlica-Wagner, and M. A. Sánchez-Conde, Legacy analysis of dark matter annihilation from the Milky Way dwarf spheroidal galaxies with 14 years of Fermi-LAT data, Phys. Rev. D109, 063024 (2024), arXiv:2311.04982 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2024

  18. [18]

    Pass 8: Toward the Full Realization of the Fermi-LAT Scientific Potential

    W. Atwoodet al.(Fermi-LAT Collaboration), Pass 8: Toward the full realization of the Fermi-LAT scientific potential, arXiv e-prints (2013), arXiv:1303.3514 [astro- ph.IM]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  19. [19]

    Fermi-LAT Observations of the Diffuse Gamma-Ray Emission: Implications for Cosmic Rays and the Interstellar Medium

    M. Ackermann, M. Ajello, W. B. Atwood, L. Baldini, J. Ballet, G. Barbiellini,et al.(Fermi-LAT Collabora- tion), Fermi-LAT observations of the diffuse gamma-ray emission: Implications for cosmic rays and the interstel- larmedium,Astrophys.J.750,3(2012),arXiv:1202.4039 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  20. [20]

    Fermi Large Area Telescope Fourth Source Catalog

    S. Abdollahiet al.(Fermi-LAT Collaboration), Fermi Large Area Telescope Fourth Source Catalog, Astro- phys. J. Suppl.247, 33 (2020), arXiv:1902.10045 [astro- ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2020

  21. [21]

    Fermi Large Area Telescope Fourth Source Catalog Data Release 4 (4FGL-DR4)

    J. Ballet, P. Bruel, T. H. Burnett, and B. Lott (Fermi- LAT Collaboration), Fermi Large Area Telescope Fourth Source Catalog Data Release 4 (4FGL-DR4), arXiv e- prints (2023), arXiv:2307.12546 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2023

  22. [22]

    emcee: The MCMC Hammer

    D. Foreman-Mackey, D. W. Hogg, D. Lang, and J. Good- man,emcee: The MCMC Hammer, Publ. Astron. Soc. Pac.125, 306 (2013), arXiv:1202.3665 [astro-ph.IM]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  23. [23]

    R. H. Byrd, P. Lu, J. Nocedal, and C. Zhu, A Limited MemoryAlgorithmforBoundConstrainedOptimization, SIAM J. Sci. Comput.16, 1190 (1995)

    work page 1995

  24. [24]

    C. Zhu, R. H. Byrd, P. Lu, and J. Nocedal, Algorithm 778: L-BFGS-B: Fortran Subroutines for Large-Scale Bound-Constrained Optimization, ACM Trans. Math. Softw.23, 550 (1997)

    work page 1997

  25. [25]

    A. A. Abdoet al.(Fermi-LAT Collaboration), The second Fermi large area telescope catalog of gamma- ray pulsars, Astrophys. J. Suppl.208, 17 (2013), arXiv:1305.4385 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  26. [26]

    PPPC 4 DM ID: A Poor Particle Physicist Cookbook for Dark Matter Indirect Detection

    M.Cirelli, G.Corcella, A.Hektor, G.Hütsi, M.Kadastik, P. Panci, M. Raidal, F. Sala, and A. Strumia, PPPC 4 DM ID: A poor particle physicist cookbook for dark mat- ter indirect detection, JCAP03, 051 (2011), [Erratum: JCAP 10, E01 (2012)], arXiv:1012.4515 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  27. [27]

    K. K. Boddy, J. Kumar, A. B. Pace, J. Runburg, and L. E. Strigari, Effective J-factors for Milky Way dwarf spheroidal galaxies with velocity-dependent annihilation, Phys. Rev. D102, 023029 (2020), arXiv:1909.13197 [astro-ph.CO]

    work page internal anchor Pith review Pith/arXiv arXiv 2020

  28. [28]

    K. K. Boddy, J. Kumar, L. E. Strigari, and M.- Y. Wang, Sommerfeld-enhanced J-factors for dwarf spheroidal galaxies, Phys. Rev. D95, 123008 (2017), arXiv:1702.00408 [astro-ph.CO]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  29. [29]

    J.L.Feng, M.Kaplinghat,andH.-B.Yu,SommerfeldEn- hancements for Thermal Relic Dark Matter, Phys. Rev. D82, 083525 (2010), arXiv:1005.4678 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  30. [30]

    M. Ibe, H. Murayama, and T. T. Yanagida, Breit-Wigner Enhancement of Dark Matter Annihilation, Phys. Rev. D 79, 095009 (2009), arXiv:0812.0072 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2009

  31. [31]

    Enhancement of Dark Matter Annihilation via Breit-Wigner Resonance

    W.-L. Guo and Y.-L. Wu, Enhancement of Dark Matter Annihilation via Breit-Wigner Resonance, Phys. Rev. D 79, 055012 (2009), arXiv:0901.1450 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2009

  32. [32]

    Murayama, Resonant annihilation of WIMP dark matter for halo gamma ray signal, arXiv e-prints (2025), arXiv:2512.01404 [hep-ph]

    H. Murayama, Resonant annihilation of WIMP dark matter for halo gamma ray signal, arXiv e-prints (2025), arXiv:2512.01404 [hep-ph]

    work page arXiv 2025

  33. [33]

    Y. Jho, J. Park, M. G. Park, and S. C. Park, Aspects of Sommerfeld enhancement in the light of halo gamma-ray excess, arXiv e-prints (2025), arXiv:2512.24662 [hep-ph]

    work page arXiv 2025

  34. [34]

    Electroweak Doublet Dark Matter for a Galactic Halo Gamma-Ray Excess

    Y. Nomura and T. Totani, Electroweak doublet dark matter for a Galactic halo gamma-ray excess, arXiv e- prints (2026), arXiv:2604.05016 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  35. [35]

    A Theory of Dark Matter

    N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer, and N. Weiner, A Theory of Dark Matter, Phys. Rev. D79, 015014 (2009), arXiv:0810.0713 [hep-ph]. 14

    work page internal anchor Pith review Pith/arXiv arXiv 2009

  36. [36]

    Constraints on decaying dark matter from the extragalactic gamma-ray background

    S. Ando and K. Ishiwata, Constraints on decaying dark matter from the extragalactic gamma-ray background, JCAP05, 024 (2015), arXiv:1502.02007 [astro-ph.CO]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  37. [37]

    Gamma-ray Constraints on Decaying Dark Matter and Implications for IceCube

    T. Cohen, K. Murase, N. L. Rodd, B. R. Safdi, and Y. Soreq,γ-ray Constraints on Decaying Dark Matter and Implications for IceCube, Phys. Rev. Lett.119, 021102 (2017), arXiv:1612.05638 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2017

  38. [38]

    Constraints on Decaying Dark Matter from the Isotropic Gamma-Ray Background

    C. Blanco and D. Hooper, Constraints on Decaying Dark Matter from the Isotropic Gamma-Ray Background, JCAP03, 019 (2019), arXiv:1811.05988 [astro-ph.HE]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  39. [39]

    Wang and K.-K

    X. Wang and K.-K. Duan, Constraining the dark mat- ter origin of the halo-like 20 GeVγ-ray excess with the AMS-02 antiproton data, JCAP05, 087 (2026), arXiv:2512.12176 [astro-ph.HE]

    work page arXiv 2026