REVIEW 3 major objections 3 minor 11 references
CosmiXs: Improved spectra for dark matter indirect detection
T0 review · 3 major / 3 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read Dark matter annihilation spectra can be computed to about 10% accuracy in the energy ranges that indirect-detection telescopes use, by tracking spin through the full shower, adding off-shell WW/ZZ final states, and retuning hadronization…
desk verdict Useful proceedings summary of the CosmiXs spectra work, but the central '10% precision' claim is extrapolated beyond what this text actually shows. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the helicity-dependent antenna shower, a parton-shower algorithm in which every branching carries information about the helicity of the emitting and emitted particles, so polarization is propagated through the entire radiation history. Around this core, the method uses matrix elements that preserve spin information, four-body treatments for WW and ZZ final states below the on-shell threshold, and a re-tuned version of the Lund string fragmentation model whose parameters are fitted to data taken at the Z-boson pole. The shower is what adds the previously missing electroweak radiation, including all trilinear boson interactions and soft-coherence effects, while the fragmentation tune fixes the non-perturbative step that converts colored partons into the stable hadrons whose decays feed the cosmic-messenger fluxes.
What would settle it
Compare the event generator, with the paper's re-tuned hadronization parameters, against measured inclusive hadron momentum spectra from electron-positron annihilation at center-of-mass energies between 10 and 90 GeV; a systematic disagreement larger than 10% in the yields of pions, kaons, protons or photons at these lower energies would show that the Z-pole tune does not transfer to dark matter masses near 5 GeV.
Extended reading notes
Core claim
The central claim is that the source spectra of stable particles from dark matter annihilation are not as uncertain as the spread between existing public tools suggests, provided three physical effects are treated correctly. The paper accomplishes this by generating annihilation matrix elements with complete spin information, evolving them with a helicity-dependent antenna shower so that electroweak radiation is tracked through the full radiation history rather than added as an afterthought, computing the WW and ZZ channels as full four-body processes down to DM masses of 5 GeV, and re-tuning the Lund string fragmentation parameters to electron-positron annihilations measured at the Z-boson pole. The result is a set of spectra for gamma rays, positrons, antiprotons and neutrinos, for DM masses between 5 and 100 GeV, released publicly. The paper states that this yields a precision below or around 10% in the energy regions important for dark matter indirect detection experiments, and that deviations with respect to previous sets reach tens of percent in some kinematic regions.
Load-bearing premise
The claim of below-or-around 10% precision rests on the assumption that the hadronization parameters, fixed once from electron-positron collisions at the Z-pole, describe fragmentation equally well for every dark matter annihilation channel and for energies down to about 10 GeV; if that transfer fails, the error budget for the lowest dark matter masses in the set near 5 GeV would be larger than claimed.
Editorial extensions
If this is right
- For dark matter masses between 5 and 100 GeV, the new spectra reduce the theory uncertainty on gamma-ray, positron, antiproton and neutrino yields to about 10% in the energy windows where searches are most sensitive.
- Annihilation into $W^+W^-$ and $ZZ$ is described as a full four-body process, so predictions remain valid below the on-shell gauge-boson threshold down to $m_\chi=5$ GeV.
- The helicity-aware treatment changes the low-energy ($x\lesssim 10^{-2}$) yields of photons and positrons relative to standard spectra, with differences reaching factors of 2-3 in some channels, which directly affects the interpretation of any putative signal.
- Two new annihilation channels ($HZ$ and $\gamma Z$) are added, and the loop-induced channels $gg$, $\gamma\gamma$ and $\gamma Z$ are computed with full one-loop matrix elements rather than effective couplings.
- Because the spectra are public and the same results can be rescaled to dark matter decay, the improved predictions can be adopted immediately by current and planned experiments.
Reading between the lines
- Editorial inference: the same helicity-dependent showering pipeline is directly transferable to collider processes that produce polarized $W$ and $Z$ bosons from heavy new-physics resonances, where the default treatment averages over spin and may bias kinematic distributions.
- Editorial inference: the claimed 10% precision is conditional on a single Z-pole tune; independent low-energy fragmentation data at center-of-mass energies far below the Z pole would provide a decisive cross-check of the 5-10 GeV DM mass range.
- Editorial inference: because the method removes the on-shell approximation for the weak-boson channels, it suggests the applicable mass window could be extended beyond 100 GeV, where the matching between the shower and analytic high-mass fragmentation would become the next limiting uncertainty.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings contribution describes CosmiXs, a publicly released set of spectra for stable particles produced in dark matter (DM) annihilation, covering gamma rays, positrons, antiprotons, neutrinos, and antinuclei for DM masses between 5 and 100 GeV. The pipeline combines MadDM matrix elements with Pythia 8 and the Vincia helicity-dependent antenna shower, includes off-shell WW/ZZ four-body decays, adds two new annihilation channels (HZ and gamma Z), and introduces a new tuning of the Lund fragmentation parameters to LEP Z-pole data. The central claim is that the resulting spectra achieve a precision below or around 10% in the energy regions critical for DM indirect detection experiments.
Significance. If the 10% precision claim holds, CosmiXs would be a valuable public resource for analyses with Fermi-LAT, AMS-02, CTA, and HAWC, improving on PPPC and QCDUnc in its treatment of helicity-dependent electroweak radiation and off-shell gauge-boson production. The paper's strengths are concrete: it ships a public GitHub repository, uses an established toolchain (MadDM + Pythia 8 + Vincia), and compares directly with PPPC, HDMSpectra, and QCDUnc. However, the quantitative precision claim is the central new statement, and it is not established by the evidence included in this proceedings.
major comments (3)
- [Section 4] The concluding statement that 'Our results lead to a precision below or around 10%' is a quantitative precision claim that is not substantiated in this manuscript. The only quantitative evidence is Fig. 1, which shows ratio comparisons without uncertainty bands; the left panel (mχ = 100 GeV, e+e− to gamma rays) does not probe the low-x hadronization region where the 10-50% hadronization uncertainties quoted in Section 2 are largest, and the right panel (mχ = 10 TeV, W+W− to positrons) lies outside the declared 5-100 GeV mass range. Since Section 3 defers the detailed derivation to ref. [10], the proceedings should either include the validation from [10] (for example, uncertainty bands or a closure test at low center-of-mass energy) or restate the 10% figure as a result of [10] rather than as a new conclusion of this paper.
- [Section 2] The new tuning of the Lund fragmentation parameters is described only qualitatively: no fitted parameter values, goodness-of-fit, or validation against data below the Z pole are reported. Given that the paper itself quotes hadronization uncertainties of 10% to about 50% depending on the kinematical region, DM mass, and annihilation channel (refs. 7,8), the transfer of a Z-pole tune to annihilations with mχ = 5 GeV (effective center-of-mass energy around 10-20 GeV) is not self-evident and is a central assumption behind the 10% precision claim. Please add at least a brief quantitative summary of the tune's performance and a low-energy validation, or explicitly point to the section of [10] where this is provided.
- [Figure 1] The right panel of Fig. 1 uses mχ = 10 TeV, which conflicts with the stated mass range of the provided spectra (5-100 GeV) and therefore does not support the precision claimed in that range. Please replace it with an in-range example or clearly label it as an illustration of high-mass behavior intended only for comparison with HDMSpectra.
minor comments (3)
- [Abstract] The phrase 'publicly distributed this GitHub repository' is missing the preposition 'in'.
- [Figure 1 caption] Please specify whether the ratio curves are computed with respect to the VINCIA spectra using the same retuned fragmentation parameters, and state the statistical precision of the Monte Carlo samples so that the reader can judge the significance of the deviations.
- [Sections 1 and 3] The companion paper [10] is cited only in Section 3; citing it earlier in the Introduction would help readers locate the full derivation and validation of the pipeline.
Circularity Check
Minor self-citation behind the 10% precision claim; the spectral derivation itself is not circular.
-
self citation load bearing
[Section 4 (Conclusions), see also Section 2 and refs. [7,8,10]]
"We also improved on the tuning of the Pythia 8 parameters for the Lund fragmentation function using LEP data at the Z-pole. Our results lead to a precision below or around 10% in the energy regions critical for DM indirect detection experiments and provide vital input for LHAASO, HAWC, CTA and AMS experiments."
The paper asserts, without in-text derivation, that its spectra have 'a precision below or around 10%'. The only support offered is the immediately preceding statement of a self-authored Lund tune and a pointer to companion paper [10] by the same six authors; Section 2's uncertainty range (10-50%) is cited to refs. [7,8], which share authors Jueid and Ruiz de Austri. The 10% precision claim therefore rests on a chain of self-citations estimating the same quantity, rather than on a test against data not used in the fit. It is minor because the spectral pipeline itself is not fit to DM observables and uses external tools (MadDM, Pythia, Vincia).
full rationale
The core spectral calculation is not circular: the Lund fragmentation parameters are fitted to public LEP Z-pole data, not to dark-matter indirect-detection observables, and the spectra are produced by independent external event generators. The paper's central quantitative selling point, the 'precision below or around 10%' claim, is however not demonstrated in this proceedings; it is deferred to the self-authored companion paper [10] and to QCDUnc uncertainty estimates [7,8] whose author lists overlap with the present work. Additionally, the positron comparison shown in Fig. 1 is for m_chi = 10 TeV, outside the stated 5-100 GeV range, and the 100 GeV gamma-ray example does not exercise the low-energy hadronization regime where the Z-pole tune is least tested. These are evidentiary gaps and self-citation concerns rather than a fit-renamed-as-prediction or an equation-level circular reduction; the score of 2 reflects one minor self-citation issue in an otherwise self-contained derivation.
Assumptions & free parameters
free parameters (1)
- Lund string fragmentation parameters (four parameters of the left-right symmetric fragmentation function) =
not specified in this proceedings; tuned to LEP Z-pole data
assumptions (3)
- domain assumption Universality of the Lund fragmentation function: parameters fitted at the Z pole transfer to DM annihilation across all channels and at lower energies.
- domain assumption Vincia's helicity-dependent antenna shower correctly resums electroweak radiation and matches the MadGraph/MadDM matrix elements.
- domain assumption Factorization of hard process, parton shower, and hadronization holds at DM energies down to 5 GeV.
Cite this review
Pith. "Pith review of CosmiXs: Improved spectra for dark matter indirect detection." pith.science (2026). https://pith.science/paper/U3MQILRY
@misc{pith2026250113281,
author = {Pith},
title = {Pith review of: CosmiXs: Improved spectra for dark matter indirect detection},
year = {2026},
howpublished = {\url{https://pith.science/paper/U3MQILRY}},
note = {Machine review of arXiv:2501.13281}
}
abstract
The spectra of stable particles produced from dark matter (DM) are one of the most important ingredients to calculate the fluxes for DM indirect detection experiments. At energies above a few GeV, most of the particles are produced following a complex sequence of phenomena including resonance decays, QED and QCD final-state radiation, radiation of weak gauge bosons, hadronization and hadron decays. In this contribution, we discuss improvements on the calculation of the energy spectra at the source using state-of-the-art tools that include effects that were not taken previously into account. We include helicity information of the particles produced in the annihilation channels, which leads to proper inclusion of electroweak radiation during the entire showering history. These effects are taken into account using the Vincia, which is based on the helicity-dependent antenna shower formalism. Off-shell contributions are also taken into account for annihilation channels into $WW$ and $ZZ$ through the four-body processes into fermions and for DM masses below the gauge boson mass. We also revisit the tune of the Lund fragmentation function parameters in Pythia using LEP data at the $Z$-boson pole. The spectra of cosmic messengers are provided for DM masses between 5 GeV and 100 GeV and are publicly distributed in this \href{https://github.com/ajueid/CosmiXs.git}{GitHub repository}.
Reference graph
Works this paper leans on
- [10]
-
[1]
J. M. Gaskins, Contemp. Phys. 57 (2016) no.4, 496-525 doi:10.1080/00107514.2016.1175160 [arXiv:1604.00014 [astro-ph.HE]]
arXiv 2016
-
[2]
R. K. Leane, [arXiv:2006.00513 [hep-ph]]
arXiv 2006
-
[3]
M. Cirelli, G. Corcella, A. Hektor, G. Hutsi, M. Kadastik, P. Panci, M. Raidal, F. Sala and A. Strumia, JCAP 03 (2011), 051 [erratum: JCAP 10 (2012), E01] doi:10.1088/1475- 7516/2012/10/E01 [arXiv:1012.4515 [hep-ph]]
arXiv 2011
-
[4]
C. W. Bauer, N. L. Rodd and B. R. Webber, JHEP 06 (2021), 121 doi:10.1007/JHEP06(2021)121 [arXiv:2007.15001 [hep-ph]]
arXiv 2021
-
[5]
N. Fischer, S. Prestel, M. Ritzmann and P. Skands, Eur. Phys. J. C 76 (2016) no.11, 589 doi:10.1140/epjc/s10052-016-4429-6 [arXiv:1605.06142 [hep-ph]]
arXiv 2016
-
[6]
S. Amoroso, S. Caron, A. Jueid, R. Ruiz de Austri and P. Skands, JCAP 05 (2019), 007 doi:10.1088/1475-7516/2019/05/007 [arXiv:1812.07424 [hep-ph]]
arXiv 2019
- [7]
Show all 11 references
-
[8]
Jueid, J
A. Jueid, J. Kip, R. R. de Austri and P. Skands, JHEP 02 (2024), 119 doi:10.1007/JHEP02(2024)119 [arXiv:2303.11363 [hep-ph]]
2024 arXiv
-
[9]
Ambrogi, C
F. Ambrogi, C. Arina, M. Backovic, J. Heisig, F. Maltoni, L. Mantani, O. Mattelaer and G. Mohlabeng, Phys. Dark Univ. 24 (2019), 100249 doi:10.1016/j.dark.2018.11.009 [arXiv:1804.00044 [hep-ph]]
2019 arXiv
-
[11]
Di Mauro, N
M. Di Mauro, N. Fornengo, A. Jueid, R. R. de Austri and F. Bellini, [arXiv:2411.04815 [astro-ph.HE]]
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.