Pith. sign in

REVIEW 6 minor 40 references

Published dark-matter annihilation limits can be refreshed for revised J-factor priors with closed-form penalty factors, without access to the experimental likelihood.

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-30 18:24 UTC pith:POW7AJIQ

load-bearing objection Clean, usable closed-form J-factor penalties plus a public-data multi-target stacker; assumptions are stated and the validations hold.

arxiv 2607.26876 v1 pith:POW7AJIQ submitted 2026-07-29 hep-ph astro-ph.COastro-ph.HEastro-ph.IM

Revising Indirect Dark Matter Constraints with Updated Astrophysical J-Factor Priors

Giacomo D'Amico This is my paper
classification hep-ph astro-ph.COastro-ph.HEastro-ph.IM
keywords indirect dark matter detectionJ-factorastrophysical uncertaintiesprofile likelihoodlimit reinterpretationdwarf spheroidal galaxiesgamma-ray constraints
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.

Gamma-ray searches for annihilating dark matter publish upper limits that are locked to whatever astrophysical J-factor prior was assumed at the time. When better stellar kinematics revise that J-factor, the published numbers risk going stale unless the whole experimental analysis is redone. This paper shows that, under a quadratic approximation of the likelihood and a null-signal hypothesis, those limits can be updated with a simple rescaling: multiply by the ratio of old to new central J values and by the ratio of analytical penalty factors that encode the old and new J uncertainties. Closed forms are given for both Gaussian and log-normal J priors, and the same ingredients let multiple targets be combined numerically from public single-target limits alone. The point is to keep existing constraints scientifically usable as astrophysical knowledge improves, without proprietary data or instrument response functions.

Core claim

Under a quadratic expansion of the profile likelihood in the composite signal strength z = J x and the null best-fit assumption ẑ = 0, a published upper limit on the annihilation cross section updates as ⟨σv⟩_new = ⟨σv⟩_old × (J_old/J_new) × P(σ_new)/P(σ_old). The penalty P is 1/√(1−λ r²) for a Gaussian prior on J (with r = σ_J/J_0) and √ε e^ε / (√λ σ) for a log-normal prior on ln J, with ε fixed by the confidence-level threshold λ. The same curvature coefficients reconstructed from single-target limits allow multi-target combination by a one-dimensional numerical root find.

What carries the argument

The J-factor penalty P(σ): a closed-form multiplier, derived by profiling a quadratic likelihood plus an external Gaussian or log-normal prior on J, that converts a fixed-J upper limit into the limit that includes astrophysical uncertainty (and vice versa when priors are revised).

Load-bearing premise

The full experimental likelihood must be well described by a quadratic expansion in overall signal strength around zero signal, with the dark-matter spatial shape held fixed so that new astrophysics only rescales the normalization.

What would settle it

Take a published analysis whose full likelihood is public; re-profile it with a new J prior, then apply the paper’s analytical update to the original published limit and check whether the two updated limits agree within the few-percent bias the paper’s own successive-reinterpretation Monte Carlo reports.

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

If this is right

  • Existing single-target gamma-ray limits can be rewritten for any future J-factor catalog without reopening instrument pipelines.
  • Combined multi-dwarf limits can be rebuilt from public single-target curves and quoted J uncertainties alone.
  • The same penalty structure applies, with only a kinematic swap, to lower limits on dark-matter lifetime for decaying dark matter.
  • Successive J revisions can be applied directly to the original published limit, keeping cumulative bias at the percent level in realistic uncertainty ranges.

Where Pith is reading between the lines

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

  • Any search whose signal rate factors as (parameter of interest) × (externally constrained normalization) can reuse the same penalty algebra—not only dark-matter J-factors.
  • Public limit repositories could ship an automatic ‘re-J’ layer that ingests new kinematic posteriors and emits refreshed curves on the fly.
  • If upcoming surveys drive relative J uncertainties below the regime where the Gaussian validity cut fails, the log-normal form becomes the default practical choice for long-term reinterpretation.

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

0 major / 6 minor

Summary. The paper derives a practical prescription for updating published upper limits on the dark-matter annihilation cross section (or decay lifetime) when revised astrophysical J-factor priors become available, without access to the full experimental likelihood. Under a quadratic approximation of the likelihood in the composite signal strength z = J x and the null-signal hypothesis ẑ = 0, closed-form penalty factors P are obtained for Gaussian (Eq. 29) and log-normal (Eq. 38) J priors. Published limits are then rescaled as ⟨σv⟩_new = ⟨σv⟩_old × (J_old/J_new) × P(σ_new)/P(σ_old). The same framework is extended to multi-target combinations by reconstructing per-target curvature coefficients A_i from single-target limits and solving the combined profile-likelihood condition numerically. The analytics are validated with ON/OFF toy Monte Carlo (including successive-reinterpretation stability) and by reproducing published single-target and MAGIC combined limits from public information only.

Significance. J-factor systematics dominate many indirect-detection limits, and forthcoming kinematic and photometric surveys will revise current priors, risking obsolescence of published constraints. The work supplies a lightweight, statistically motivated reinterpretation tool that complements full experimental reanalyses and the author’s prior particle-physics recasting framework. Strengths include transparent profile-likelihood derivations with correct σ → 0 limits, quantitative MC validation (Figs. 1–2), independent reproduction of external published limits (Ref. [37] and MAGIC), a usable multi-target numerical procedure, and publicly released code. Within the stated regime (quadratic likelihood, fixed spatial morphology, null best fit), the result is immediately useful to the community.

minor comments (6)
  1. [Abstract / throughout] Abstract and several places in the text: minor grammar (“and demonstrate their accuracy”); also “Therefor,” “for conveience,” and “Upda ted” in a section heading. A light copy-edit pass would clean these.
  2. [Sec. IV] Sec. IV, around Eqs. (26)–(28): the update rule is written before the Gaussian UL is fully derived and cites Eq. (28) in a way that can confuse a first reading. A short forward pointer or reordering of one paragraph would help.
  3. [Sec. VI] Sec. VI / Fig. 3: the ~2σ tension for Ursa Minor is noted but not diagnosed. A sentence on whether it is consistent with a non-quadratic likelihood feature, a different J-prior implementation, or simply statistical fluctuation would strengthen the validation narrative.
  4. [Sec. VII A] Sec. VII A, Eq. (57): emphasize more explicitly in the main text (not only in the surrounding paragraph) that A_i must be reconstructed from the original J priors used in the published single-target limits, not from the updated priors. This is easy to misapply in practice.
  5. [Fig. 1 / Sec. V] Fig. 1 caption and right panel: “uncertainty on log_10 J” should be stated consistently with the analytic σ (natural log) used in Eq. (38), or the conversion should be made explicit so readers do not mix bases when applying P(σ).
  6. [Sec. VIII] Decay case (Sec. VIII): the claim P_decay = 1/P is brief. One line showing why the lower-limit form follows from the same profiled likelihood would make the extension fully self-contained.

Circularity Check

0 steps flagged

No significant circularity: penalties are profiled from a stated quadratic likelihood, and validations use external published limits as independent benchmarks.

full rationale

The central update rule and the Gaussian/log-normal penalty factors P follow from profiling an explicitly assumed quadratic likelihood in z = J x about the null ẑ = 0 (Secs. III–IV, Eqs. 9–38). That is a closed-form derivation under stated approximations, not a fit to the validation targets and not a renaming of an empirical pattern. Toy ON/OFF Monte Carlo (Fig. 1), successive-reinterpretation stability (Fig. 2), reproduction of Ref. [37] fixed-J → profiled-J limits, and the MAGIC multi-dSph combination (Fig. 4) all treat external published limits and J tables as independent checks of the approximation. The self-citation to the author’s prior recasting paper [28] is explicitly complementary (particle-physics model changes) and is not used to force or uniqueness-justify the J-penalty formulas; A_i coefficients for multi-target combination are reconstructed from the published single-target limits under the same approximation (Eq. 57), which is inversion of the derived single-target result rather than circular prediction. No step reduces the claimed result to its inputs by construction.

Axiom & Free-Parameter Ledger

2 free parameters · 7 axioms · 0 invented entities

The central update rule rests on standard profile-likelihood practice plus a small set of domain modeling choices (quadratic likelihood in z=Jx, null best-fit, fixed spatial template, Gaussian or log-normal J priors, independent targets). No new physical entities are postulated. Free parameters are only the conventional CL threshold λ and the externally supplied J prior parameters; nothing is fitted to force the penalties.

free parameters (2)
  • λ (likelihood-ratio threshold) = 2.71 (95% CL one-sided)
    Sets the confidence level of the upper limit (λ=2.71 for one-sided 95% CL). Standard choice, not fitted to data in this work, but the numerical value of P depends on it.
  • J0, σ_J (or σ on ln J) per target
    External astrophysical prior parameters taken from the literature or from the original experimental paper; they are inputs to the update, not fitted here.
axioms (7)
  • domain assumption The binned experimental log-likelihood is locally quadratic in the composite signal strength z = J x about its minimum (Sec. III, Eqs. 6–9).
    Enables closed-form profiling; standard asymptotic approximation but not exact for low-count or highly non-Gaussian analyses.
  • domain assumption No significant excess: best-fit signal normalization ẑ = 0 (null hypothesis).
    Used throughout Secs. IV–VII to simplify penalties; paper notes limits become approximate if a real excess is present.
  • domain assumption Revised astrophysics changes only the overall J normalization, not the spatial signal morphology (footnote in Sec. III).
    Required so that a single nuisance J (or y=ln J) captures the update; explicitly out of scope otherwise.
  • domain assumption J-factor uncertainty is modeled as Gaussian in J or Gaussian in ln J.
    Standard choices in the indirect-detection literature; penalties are derived specifically for these two forms.
  • domain assumption Multiple targets contribute independent likelihood factors (Sec. VII).
    Allows the combined profile likelihood to factorize target-by-target before the numerical root find.
  • standard math Profile-likelihood ratio with threshold λ yields the quoted frequentist upper limit (Wilks/asymptotic).
    Standard statistical machinery (cited via Rolke et al.); not re-proved.
  • standard math For multi-target ε_i(x), the principal branch of the Lambert-W function gives the unique physical root (Eq. 53).
    Standard special-function identity for the profiled log-normal equation.

pith-pipeline@v1.2.0-daily-grok45 · 20350 in / 3507 out tokens · 55112 ms · 2026-07-30T18:24:12.397811+00:00 · methodology

0 comments
read the original abstract

Indirect searches for particle dark matter with gamma-ray experiments have produced a large number of constraints on the annihilation cross section (or decay lifetime) over a wide range of dark matter masses. These constraints depend critically on the assumed astrophysical $J$ factor and its uncertainty, which encodes the dark matter distribution in the target and represents the dominant source of systematic uncertainty. As improved observational data and dynamical modeling are expected to revise current $J$-factor determinations, many published limits risk becoming obsolete unless the full experimental analyses are repeated. In this work we present a general and statistically consistent framework for updating published dark matter limits when revised $J$-factor estimates become available, without requiring access to the full experimental likelihood. We derive an analytical expression that quantifies the impact of astrophysical uncertainties on dark matter limits, treating both Gaussian and log-normal priors on the $J$ factor. The formalism is validated through toy Monte Carlo simulations, including dedicated studies of its numerical stability under successive reinterpretations, and demonstrate their accuracy by reproducing published limits. Lastly, we further show that the formalism naturally extends to the combination of multiple targets through a simple numerical procedure, allowing limits to be combined and updated using only publicly available information. The proposed method is intended as a complementary reinterpretation tool for situations in which a complete experimental reanalysis is impractical, offering a practical means to preserve and extend the scientific relevance of published dark matter constraints across present and future experiments.

Figures

Figures reproduced from arXiv: 2607.26876 by Giacomo D'Amico.

Figure 1
Figure 1. Figure 1: FIG. 1. Comparison between the analytical [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Mean relative difference between upper limits ob [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Upper limits on the dark matter annihilation cross section [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Upper limits at 95% CL on the dark matter annihi [PITH_FULL_IMAGE:figures/full_fig_p010_4.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

40 extracted references · 4 linked inside Pith

  1. [1]

    Up to an irrelevant additive constant, this prior contributes to the log-likelihood as −2 lnπ(J) = (J−J 0)2 σ2 J ,(11) which when substituted in Eq

    Gaussian prior We assume a Gaussian prior for the astrophysicalJ- factor with meanJ 0 and standard deviationσ J . Up to an irrelevant additive constant, this prior contributes to the log-likelihood as −2 lnπ(J) = (J−J 0)2 σ2 J ,(11) which when substituted in Eq. (10) gives −2 lnL(x, J)≃A(xJ−ˆz) 2 +B(J−J 0)2,(12) where we recall thatz≡J xand for convenienc...

  2. [2]

    Equivalently, one assumes that y≡lnJ(16) is Gaussian distributed

    Log-normal prior In many indirect-detection analyses, the astrophysical J-factor is modeled with alog-normaluncertainty [34]. Equivalently, one assumes that y≡lnJ(16) is Gaussian distributed. So that, up to an irrelevant nor- malization constant, −2 lnπ(y) = (y−y 0)2 σ2 .(17) As done in for Eq. (12), the log-likelihood is, up to an additive constant, −2 l...

    arXiv

  3. [3]

    VII A, we tested the method against published combined upper limits from the MAGIC Col- laboration

    Validation with published combined limits As a validation of the multi-target numerical procedure described in Sec. VII A, we tested the method against published combined upper limits from the MAGIC Col- laboration. In particular, we considered the limits re- ported in Fig. 4 of Ref. [3] for theτ +τ − annihilation channel, which are obtained by combining ...

    2023

  4. [4]

    M. Doro, M. A. S´ anchez-Conde, and M. H¨ utten, arXiv preprint arXiv:2111.01198 (2021)

    Pith/arXiv arXiv 2021

  5. [5]

    Ackermannet al.(Fermi-LAT), Phys

    M. Ackermannet al.(Fermi-LAT), Phys. Rev. Lett.115, 231301 (2015), arXiv:1503.02641 [astro-ph.HE]

    Pith/arXiv arXiv 2015

  6. [6]

    V. A. Acciariet al.(MAGIC), Phys. Dark Univ.35, 100912 (2022), arXiv:2111.15009 [astro-ph.HE]

    arXiv 2022

  7. [7]

    Abramowski, F

    A. Abramowski, F. Acero, F. Aharonian, A. Akhper- janian, G. Anton, A. Barnacka, U. Barres de Almeida, A. Bazer-Bachi, Y. Becherini, J. Becker,et al., Physical Review Letters106, 161301 (2011)

    2011

  8. [8]

    Archambault, A

    S. Archambault, A. Archer, W. Benbow, R. Bird, E. Bourbeau, T. Brantseg, M. Buchovecky, J. Buckley, V. Bugaev, K. Byrum,et al., Physical Review D95, 082001 (2017)

    2017

  9. [9]

    Caoet al.(LHAASO), Phys

    Z. Caoet al.(LHAASO), Phys. Rev. Lett.133, 061001 (2024), arXiv:2406.08698 [astro-ph.HE]

    arXiv 2024

  10. [10]

    Battaglia and C

    G. Battaglia and C. Nipoti, Nature Astronomy6, 659 (2022)

    2022

  11. [11]

    M. G. Walker, M. Mateo, E. W. Olszewski, O. Y. Gnedin, X. Wang, B. Sen, and M. Woodroofe, The Astrophysical Journal667, L53 (2007)

    2007

  12. [12]

    G. D. Martinez, Monthly Notices of the Royal Astronom- ical Society451, 2524 (2015)

    2015

  13. [13]

    L. E. Strigari, Reports on Progress in Physics81, 056901 (2018)

    2018

  14. [14]

    S´ anchez, D

    E. S´ anchez, D. Collaboration,et al., Nuclear and particle physics proceedings273, 302 (2016)

    2016

  15. [15]

    H. e. Flewelling, E. Magnier, K. Chambers, J. Heasley, C. Holmberg, M. Huber, W. Sweeney, C. Waters, A. Calamida, S. Casertano,et al., The Astrophysical Journal Supplement Series251, 7 (2020)

    2020

  16. [16]

    Aihara, N

    H. Aihara, N. Arimoto, R. Armstrong, S. Arnouts, N. A. Bahcall, S. Bickerton, J. Bosch, K. Bundy, P. L. Capak, J. H. Chan,et al., Publications of the Astronomical So- ciety of Japan70, S4 (2018)

    2018

  17. [17]

    Mellier, A

    Y. Mellier, A. Abdurro´ uf, J. A. Barroso, A. Ach´ ucarro, J. Adamek, R. Adam, G. Addison, N. Aghanim, M. Aguena, V. Ajani,et al., Astronomy & Astrophysics (2024)

    2024

  18. [18]

    C. C. van Damme, T. Prod’Homme, K. Isaak, T. R¨ uhl, and M. Sirianni, inSpace Telescopes and Instrumentation 2024: Optical, Infrared, and Millimeter Wave, Vol. 13092 (SPIE, 2024) pp. 136–144

    2024

  19. [19]

    D. E. S. Collaboration:, T. Abbott, F. Abdalla, J. Aleksi´ c, S. Allam, A. Amara, D. Bacon, E. Balbinot, M. Banerji, K. Bechtol,et al., Monthly Notices of the Royal Astronomical Society460, 1270 (2016)

    2016

  20. [20]

    Aihara, Y

    H. Aihara, Y. AlSayyad, M. Ando, R. Armstrong, J. Bosch, E. Egami, H. Furusawa, J. Furusawa, S. Ha- rasawa, Y. Harikane,et al., Publications of the Astro- nomical Society of Japan74, 247 (2022)

    2022

  21. [21]

    D. G. York, J. Adelman, J. E. Anderson, Jr, S. F. Ander- son, J. Annis, N. A. Bahcall, J. Bakken, R. Barkhouser, S. Bastian, E. Berman,et al., The Astronomical Journal 120, 1579 (2000)

    2000

  22. [22]

    A. Dey, D. J. Schlegel, D. Lang, R. Blum, K. Burleigh, X. Fan, J. R. Findlay, D. Finkbeiner, D. Herrera, S. Juneau,et al., The Astronomical Journal157, 168 (2019)

    2019

  23. [23]

    S. Jin, S. C. Trager, G. B. Dalton, J. A. L. Aguerri, J. Drew, J. Falc´ on-Barroso, B. T. G¨ ansicke, V. Hill, A. Iovino, M. M. Pieri,et al., Monthly Notices of the Royal Astronomical Society530, 2688 (2024)

    2024

  24. [24]

    X. Fan, G. Knapp, M. A. Strauss, J. E. Gunn, R. H. Lupton, ˇZ. Ivezi´ c, C. M. Rockosi, B. Yanny, S. Kent, D. P. Schneider,et al., The Astronomical Journal119, 928 (2000)

    2000

  25. [25]

    Prusti, J

    T. Prusti, J. De Bruijne, A. G. Brown, A. Vallenari, C. Babusiaux, C. Bailer-Jones, U. Bastian, M. Biermann, D. W. Evans, L. Eyer,et al., Astronomy & astrophysics 595, A1 (2016)

    2016

  26. [26]

    A. G. Brown, A. Vallenari, T. Prusti, J. H. De Bruijne, C. Babusiaux, C. A. Bailer-Jones, M. Biermann, D. W. Evans, L. Eyer, F. Jansen,et al., Astronomy & astro- physics616, A1 (2018)

    2018

  27. [27]

    Collaborationet al., Astronomy & Astrophysics674, A1 (2023)

    G. Collaborationet al., Astronomy & Astrophysics674, A1 (2023)

    2023

  28. [28]

    J. I. Read, M. Walker, and P. Steger, Monthly Notices of the Royal Astronomical Society484, 1401 (2019)

    2019

  29. [29]

    Hayashi, M

    K. Hayashi, M. Chiba, and T. Ishiyama, The Astrophys- ical Journal904, 45 (2020)

    2020

  30. [30]

    Errani, J

    R. Errani, J. F. Navarro, J. Pe˜ narrubia, B. Famaey, and R. Ibata, Monthly Notices of the Royal Astronomical So- ciety519, 384 (2023)

    2023

  31. [31]

    D’Amico, M

    G. D’Amico, M. Doro, and M. De Caria, Physics of the Dark Universe , 102154 (2025)

    2025

  32. [32]

    Rico, Galaxies8, 25 (2020)

    J. Rico, Galaxies8, 25 (2020)

    2020

  33. [33]

    Lefranc, E

    V. Lefranc, E. Moulin, P. Panci, F. Sala, and J. Silk, Journal of Cosmology and Astroparticle Physics2016 (09), 043

  34. [34]

    D’Amico, Universe8, 90 (2022), arXiv:2202.04590 [astro-ph.IM]

    G. D’Amico, Universe8, 90 (2022), arXiv:2202.04590 [astro-ph.IM]

    Pith/arXiv arXiv 2022

  35. [35]

    Conrad, Astroparticle Physics62, 165 (2015)

    J. Conrad, Astroparticle Physics62, 165 (2015)

    2015

  36. [36]

    Wald, Transactions of the American Mathematical so- ciety54, 426 (1943)

    A. Wald, Transactions of the American Mathematical so- ciety54, 426 (1943)

    1943

  37. [37]

    Cherenkov Telescope Array Observatory, Prospects for dark matter observations in dwarf spheroidal galaxies with the Cherenkov Telescope Array Observatory (2025), arXiv:2508.19120 [astro-ph.HE]

    arXiv 2025

  38. [38]

    W. A. Rolke, A. M. Lopez, and J. Conrad, Nucl. Instrum. Meth. A551, 493 (2005), arXiv:physics/0403059

    Pith/arXiv arXiv 2005

  39. [39]

    A. B. Pace and L. E. Strigari, Monthly Notices of the Royal Astronomical Society482, 3480 (2019)

    2019

  40. [40]

    K. Abe, S. Abe, J. Abhir, A. Abhishek, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella,et al., Monthly Notices of the Royal Astronomical Society544, 2946 (2025)

    2025