Pith. sign in

REVIEW 5 major objections 4 minor 17 references

All-Sky Cosmic-Ray Anisotropy Update at Multiple Energies

T0 review · 5 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash

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

Pith's one-line read The paper claims the first all-sky cosmic-ray anisotropy study above 10 TeV, built by combining HAWC and IceCube maps, and reports an energy-dependent anisotropy with a rapid phase transition between 40 and 76 TV.

desk verdict A useful first all-sky HAWC+IceCube combination above 10 TeV, but the rigidity-dependence claim rests on a matching procedure that needs a direct post-matching compatibility check. read the letter →

arxiv 2507.07070 v1 pith:UF6QA542 submitted 2025-07-09 astro-ph.HE

Juan Carlos Díaz-Vélez , Riya Yogesh Kore , Paolo Desiati , Ferris Wolf (for the HAWC , IceCube Collaborations) This is my paper · ORCID
classification astro-ph.HE
keywords cosmic-rayanisotropyall-skymapangularpowerspectrumHAWCIceCuberigiditydependenceenergy-dependentTeVcosmicrays
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Cosmic rays do not arrive uniformly from all directions, but measuring the pattern is hard because every observatory sees only part of the sky, and partial coverage distorts the spherical-harmonic components of the anisotropy. This paper combines eight years of HAWC data with twelve years of IceCube data to build sky maps covering 93 percent of the sky across eleven bins from 0.6 to 280.5 TV. The authors claim that this is the first all-sky anisotropy study with primary energies above 10 TeV, and that the combined maps and their angular power spectra remove most of the bias caused by partial sky coverage. The central result is an energy-dependent anisotropy whose large-scale structure changes rapidly between 40 and 76 TV, matching earlier partial-sky measurements.

What carries the argument

The relative-intensity map $δI_j = (N_j - \langle N_j\rangle)/\langle N_j\rangle$ per HEALPix pixel, reconstructed from combined observatories with the maximum-likelihood method of Ahlers et al.; the matching of HAWC energy bins to IceCube cuts by Kolmogorov–Smirnov comparison of rigidity distributions under the GSF composition model; and the pseudo-angular power spectrum $C_\ell = \frac{1}{2\ell+1}\sum_m |a_{\ell m}|^2$, whose low-$\ell$ modes carry the large-scale anisotropy. The combination step is what removes the partial-sky bias: with nearly full coverage the $a_{\ell m}$ correlations are suppressed, so the dipole and quadrupole terms can be trusted.

What would settle it

Recompute the combined maps and the phase-transition location using the H3a composition model instead of GSF; if the rapid change between 40 and 76 TV shifts, smears, or disappears, the transition is an artifact of the composition assumption. A second check: take the 126 TeV overlapping region where the two detectors are statistically incompatible and see whether the incompatibility vanishes when bins are re-matched by rigidity.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that a rigidity-matched combination of HAWC and IceCube data yields a nearly full-sky view of the cosmic-ray arrival-direction distribution from 0.6 to 280.5 TV, and that this view shows a rapid phase transition in the large-scale anisotropy between 40 and 76 TV. Because the two detectors have complementary fields of view and the combined map covers 93 percent of the sky, the authors argue that the angular power spectrum is largely free of the mode-correlation bias that afflicts partial-sky analyses. The energy dependence is confirmed by HAWC alone up to about 0.5 PeV, and the agreement between the two instruments in their overlapping field of view improves when bins are matched by rigidity using the GSF composition model rather than by energy.

Load-bearing premise

The result depends on the assumption that cosmic-ray arrival directions organize by rigidity rather than energy, and that the GSF composition model converts HAWC energy cuts into IceCube rigidity cuts with enough accuracy for the matching.

Editorial extensions

If this is right

  • If the claims hold, future analyses can treat the combined HAWC–IceCube maps as a near-full-sky reference for cosmic-ray anisotropy from roughly TeV to PeV energies.
  • The 40–76 TV phase transition becomes a fixed feature that theories of cosmic-ray transport in the local interstellar medium will have to reproduce.
  • Rigidity, not energy, should be used as the organizing variable when comparing anisotropy measurements from different experiments.
  • The angular power spectrum from full-sky maps can be compared directly with partial-sky results to quantify and remove the bias.
  • IceCube's anisotropy structures evolving faster with energy than HAWC's is expected if both are ordered by rigidity; the combined maps make this effect explicit.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the rigidity hypothesis is right, the dominant systematic uncertainty in any future all-sky anisotropy measurement shifts from statistics to the assumed cosmic-ray composition model; choosing a different model should move the phase-transition energy.
  • The same maximum-likelihood combination could be applied to observatories with complementary fields of view at lower and higher energies, extending this analysis beyond the 0.6–280.5 TV window the paper covers.
  • A direct test of the phase transition is to split the combined data by season or by detector and look for a stable 40–76 TV transition in each subset, since a transition caused by the matching procedure would not appear consistently.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

5 major / 4 minor

Summary. This proceedings paper reports a preliminary all-sky cosmic-ray anisotropy analysis combining 8 years of HAWC data with 12 years of IceCube data, using a likelihood-based combination method. Eleven energy/rigidity bins span 0.6 TV to 280.5 TV, with IceCube energy cuts adjusted to match HAWC rigidity distributions under the GSF composition model. The paper claims to be the first all-sky anisotropy study above 10 TeV and reports a rapid phase transition in the anisotropy pattern between 40 and 76 TV.

Significance. If the rigidity-matching procedure is valid, this is a valuable step: it extends all-sky anisotropy studies to higher energies, uses a longer HAWC exposure than earlier work, and makes a concrete cross-check of energy-dependent versus rigidity-dependent anisotropy. The underlying combination method (Ahlers et al. 2016) is well established, and the paper makes the new IceCube energy cuts explicit in Table 1. However, the central high-energy conclusion rests on the rigidity-matching assumption, and the current manuscript does not provide a post-matching validation of that assumption, so the significance cannot yet be fully assessed.

major comments (5)
  1. [Sec. 5, Fig. 2, Table 1] The rigidity-matching procedure is circular as a test of the claim that angular distributions depend on rigidity. The IceCube energy cuts are selected by minimizing the KS distance between GSF-based rigidity distributions, so the rigidity distributions agree by construction. The manuscript reports no post-matching comparison of the measured anisotropy in the overlapping field of view, such as a chi-square or KS value between the IceCube and HAWC one-dimensional projections after applying the new cuts. Without this check, the statement in Sec. 6 that 'the best-matching energy bins are consistent with this hypothesis' is not supported, and the 40-76 TV phase transition could be an artifact of the matching procedure. Please report the post-matching compatibility statistics for all seven combined bins, or show the overlay projections after matching.
  2. [Sec. 5.1, Fig. 5] The paper claims that the combined maps and angular power spectra 'largely eliminate biases that result from partial sky coverage,' but it shows only the pseudo-angular power spectrum, with no deconvolution or window-function correction. Differences between the HAWC-only (black) and combined (red) spectra can therefore reflect the different sky masks rather than the intrinsic anisotropy. Please either apply a deconvolution or mode-coupling correction, or explicitly label the result as the pseudo-APS and discuss how the window function affects the comparison.
  3. [Sec. 5, Fig. 1] There is an internal inconsistency between the text and the Fig. 1 caption on whether the 10 TeV overlap is compatible. The text states 'At 10 TeV energy, the distributions are qualitatively different but are statistically compatible,' while the caption states that 'at lower and higher energies' the distributions 'are not statistically compatible.' Please clarify, and report the actual test statistic values. If the 10 TeV point is incompatible, this matters directly for the boundary between HAWC-only and combined maps at 4.1 TV in Table 1.
  4. [Sec. 5, Fig. 3] The 'rapid phase transition between 40 TV and 76 TV' is based on visual inspection of the maps. To make this claim quantitatively testable, please provide a measure such as the dipole phase and amplitude evolution across bins, or a change-point statistic with uncertainties. Without this, the claim cannot be rigorously compared with previous partial-sky measurements.
  5. [Sec. 6, Fig. 5] Only statistical uncertainties are shown, and the sensitivity of the rigidity-matched bins to the composition model (e.g., H3a versus GSF, which are both shown in Fig. 2) is not propagated into the combined maps or the phase-transition claim. The paper acknowledges the need for composition-model systematics in Sec. 6, but that is precisely the input that currently prevents a quantitative assessment of the central result.
minor comments (4)
  1. [Table 1 caption] The caption says 'GST composition model' but the text and Fig. 2 refer to the GSF model; please make the notation consistent.
  2. [Fig. 5] The y-axis label 'c' should be written as C_ell for clarity, and the caption should state more explicitly that the first four panels are HAWC-only while the last seven combine IceCube and HAWC.
  3. [Sec. 2] The sentence 'The IceCube dataset is described in detail in [1], though energy cuts have been adjusted' could state explicitly which IceCube analysis (event selection, declination range) is used, since the rigidity matching changes the effective energy range.
  4. [Abstract and Sec. 6] The phrase '93% coverage of the sky' is not further justified; the maps are described as covering 70N to 90S. Please state how the 93% is computed, or soften the wording to 'approximately 93% of the sky.'

Circularity Check

1 steps flagged · score 5.0 of 10

The rigidity-dependence justification for combining IceCube and HAWC is partly circular: the IceCube cuts are optimized so that rigidity distributions match under the GSF model, the promised chi-squared anisotropy check is never reported, and 'the best-matching energy bins are consistent with this hypothesis' then cites the fit itself as evidence; the HAWC-only maps and the cross-check against…

  1. fitted input called prediction [Section 5 (Combined Anisotropy with IceCube), supported by Section 6 (Conclusions)]
    "In order to find the optimal IceCube energy cuts to match the rigidity for each HAWC bin, we perform a scan and use a Kolmogorov-Smirnov test to compare the rigidity distributions assuming the GSF [16] composition model, as shown in Fig 2. ... Monte Carlo studies suggest that the angular distributions depend on rigidity rather than energy. The best-matching energy bins are consistent with this hypothesis."

    The IceCube energy cuts are fit by a KS scan whose objective is agreement between IceCube and HAWC rigidity distributions under the GSF composition model, so the selected bins agree in rigidity by construction. The paper's in-paper support for rigidity dependence ('The best-matching energy bins are consistent with this hypothesis') then appeals to that same selection. The promised anisotropy-level test ('using a chi2 test and assess whether we achieve better agreement') is never quantified, and Fig. 1 shows that before matching, the 126 TeV distributions are not statistically compatible. Thus the combined maps above 10 TeV and the 40-76 TV phase transition rest on a hypothesis whose in-paper validation reduces, as written, to the fitting procedure.

full rationale

The paper's core data products (the 11 HAWC maps and the combined relative-intensity and significance maps) are measured quantities that do not reduce to any fit, and the rapid phase transition is explicitly cross-checked against previous independent measurements ([1]), so the analysis is not globally circular. The circular element is confined to the rigidity-matching step that justifies combining IceCube and HAWC above 10 TeV. The IceCube energy cuts are selected by minimizing the KS distance between rigidity distributions computed assuming the GSF composition model; the matching bins therefore agree in rigidity by construction. The paper announces a chi-squared comparison of the actual one-dimensional anisotropy distributions but never reports its outcome, and instead concludes that 'the best-matching energy bins are consistent with this hypothesis.' As written, the in-paper evidence for rigidity dependence is the success of a selection procedure that already assumed rigidity dependence, which is the fitted-input-called-prediction pattern. The external Monte Carlo studies invoked in the conclusions could supply independent support, but they are not cited or detailed here, so they cannot break the circularity within this manuscript. Because the lower-energy maps (bins 0-3, HAWC only) and the external consistency of the phase transition do not depend on the fit, the circularity is partial, supporting a score of 5 rather than 6-8. The paper's honest acknowledgement that composition-model uncertainties require further investigation further limits the severity, but the missing chi-squared result leaves the high-energy combination's validity resting on the fitted matching.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The central claim rests on external estimators, a composition model, and an assumed rigidity dependence, but no new particles, forces, or entities are introduced. The main maps are data products, so the circularity burden comes from the rigidity-matching optimization rather than from a derivation.

free parameters (4)
  • Atmospheric pressure coefficient beta = -0.0086 hPa^-1
    Determined experimentally from HAWC local measurements (Sec. 2.1) and used in the event weighting formula Eq. 2.
  • Spectral index gamma(E) for Compton-Getting correction = from HAWC all-particle spectrum [11]
    Used in Eqs. 1 and 2 to model the solar dipole correction; the value is taken from an external fit to HAWC data.
  • IceCube energy cuts for rigidity-matched bins = 10.0, 14.8, 30.2, 53.4, 310.5, 725.3, 1716 TeV (Table 1)
    Selected in Sec. 5 by scanning and minimizing a Kolmogorov-Smirnov test between IceCube and HAWC rigidity distributions. These cuts define the combined high-energy maps.
  • Smoothing radii and significance thresholds = 10, 20, 30 degrees depending on energy bin
    Adjusted by hand for higher energy bins 'to compensate for decreasing statistics' (Sec. 5, Figs. 3-4), affecting the apparent significance of structures.
assumptions (4)
  • domain assumption The arrival direction distribution of cosmic rays is rigidity-dependent rather than energy-dependent across the matched bins.
    Sec. 5 states 'If we assume... we can compare...' and uses this to pair IceCube and HAWC energy bins; the rigidity-dependence conclusion draws strength from this assumption.
  • domain assumption The GSF (and H3a) composition models correctly describe the mean charge of cosmic rays as a function of energy for both detectors.
    Used in Fig. 2 to convert energy to rigidity and to select matching bins; errors in composition propagate directly into the rigidity-matched maps.
  • domain assumption The Ahlers et al. likelihood reconstruction method recovers the true anisotropy without detector simulations and without systematic bias in the overlap region.
    Sec. 3 invokes [13]; the combined maps rely on this method's unbiasedness.
  • standard math The standard Compton-Getting and atmospheric pressure correction model accounts for solar-dipole and weather effects.
    Eqs. 1-2 implement this standard correction, with beta measured and gamma(E) taken from an external HAWC spectrum fit.

how reviews work

0 comments
Cite this review

Pith. "Pith review of All-Sky Cosmic-Ray Anisotropy Update at Multiple Energies." pith.science (2026). https://pith.science/paper/UF6QA542

@misc{pith2026250707070,
  author       = {Pith},
  title        = {Pith review of: All-Sky Cosmic-Ray Anisotropy Update at Multiple Energies},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/UF6QA542}},
  note         = {Machine review of arXiv:2507.07070}
}
abstract

We present preliminary results on an updated full-sky analysis of the cosmic-ray arrival direction distribution with data collected by the High-Altitude Water Cherenkov (HAWC) Observatory and IceCube Neutrino Observatory with complementary field of views covering a large fraction of the sky. This study extends the energy range to higher energies. The HAWC Observatory, located at 19$^{\circ}$N has analyzed 8 years of cosmic-ray data over an energy range between 3.0 TeV and 1.0 PeV and confirms an energy-dependent anisotropy in the arrival direction distribution of cosmic rays seen by other experiments. Combined with recently published results from IceCube with 12 years of data, the combined sky maps with 93\% coverage of the sky -- between 70$^{\circ}$N and 90$^{\circ}$S -- and the corresponding angular power spectra largely eliminate biases that result from partial sky coverage.

Figures

Figures reproduced from arXiv: 2507.07070 by the authors.

Figure 1
Figure 1. One-dimensional projection of relative intensity in the overlapping field of view of IceCube and HAWC. At 48 TeV energies, there is good statistical agreement between both experiments (top right). However at lower and higher energies (top left and center bottom), the two distributions are not statistically compatible. 5. Combined Anisotropy with IceCube The overlapping FoV between the two observatories serves as a c… view at source ↗
Figure 2
Figure 2. Left: mean logarithm of the particle charge Z for cosmic rays detected by IceCube (blue) and HAWC (red) assuming the Gaisser H3a [15] (dashed) and GSF [16] (solid) composition models. IceCube data is dominated by protons and light elements at lower energies than HAWC, but becomes heavier above about 200 TeV. Right: Assuming a rigidity–dependent angular distribution of cosmic rays, we find the most compatible energy … view at source ↗
Figure 3
Figure 3. Relative intensity for 11 HAWC energy bins [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: Li–Ma significance for 11 HAWC energy bins [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 5
Figure 5. Figure 5: Angular power spectra for 11 HAWC sky maps binned in energy (black circles) and 7 combined all￾sky IceCube+HAWC maps using rigidity-driven pairs of energy bins (red squares). Error bars represent statistical uncertainties. The shaded regions correspond to the isotropic…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

17 extracted references · 16 canonical work pages

  1. [1]

    Abbasi, M

    R. Abbasi, M. Ackermann, J. Adams,et al.,Astrophys. J.981 no. 2, (Mar, 2025) 182

    work page 2025

  2. [2]

    Sommers,Astropart

    P. Sommers,Astropart. Phys.14no. 4, (2001) 271 – 286. 7 All-Sky Cosmic-Ray Anisotropy Update at Multiple Energies 0 5 10 15 20 25 30 35 40 10 11 10 10 10 9 10 8 10 7 10 6 c PRELIMINARY 0.6 TV HAWC 180 180 45 20 10 5 0 5 10 15 20 25 30 35 40 10 11 10 10 10 9 10 8 10 7 10 6 c PRELIMINARY 1.1 TV HAWC 180 180 45 20 10 5 0 5 10 15 20 25 30 35 40 10 11 10 10 10...

    work page 2001

  3. [3]

    Ahlers,ApJL886 no

    M. Ahlers,ApJL886 no. 1, (Nov, 2019) L18

    work page 2019

  4. [4]

    IceCubeCollaboration, M. G. Aartsenet al.,JINST12 no. 03, (2017) P03012

    work page 2017

  5. [5]

    Abeysekara, A

    A. Abeysekara, A. Albert, R. Alfaro,et al.,NIMPA1052 (2023) 168253

    work page 2023

  6. [6]

    Whitehornet al.,Comp

    N. Whitehornet al.,Comp. Phys. Com.184 no. 9, (2013) 2214 – 2220

    work page 2013

  7. [7]

    M. G. Aartsenet al., Astrophys. J.826 (Aug., 2016) 220

    work page 2016

  8. [8]

    Verpoest, D

    S. Verpoest, D. Soldin, and P. Desiati,Astropart. Phys.161 (2024) 102985

    work page 2024

Show all 17 references
  1. [9]

    A. H. Compton and I. A. Getting,Phys. Rev.47 (Jun, 1935) 817–821

    1935

  2. [10]

    PAL: A Positional Astronomy Library,

    T. Jenness and D. S. Berry, “PAL: A Positional Astronomy Library,” inAstro. Data Analysis Software and Systems XXII, D. N. Friedel, ed., vol. 475 ofASP Conf. Series, p. 307. 2013

    2013

  3. [11]

    J. A. Morales-Soto and J. C. Arteaga-Velázquez,PoSICRC2023(2023) 364

    2023

  4. [12]

    K. M. Górskiet al., Astrophys. J.622 (2005) 759–771

    2005

  5. [13]

    Ahlerset al., Astrophys

    M. Ahlerset al., Astrophys. J.823 no. 1, (2016) 10

    2016

  6. [14]

    T. P. Li and Y. Q. Ma,Astrophys. J.272 (1983) 317–324

    1983

  7. [15]

    T. K. Gaisser,Astropart. Phys.35 (2012) 801–806

    2012

  8. [16]

    Dembinski, R

    H. Dembinski, R. Engel, A. Fedynitch,et al., PoSICRC2017(2017) 533

    2017

  9. [17]

    A. U. Abeysekaraet al., Astrophys. J.865 no. 1, (2018) 57. 8 All-Sky Cosmic-Ray Anisotropy Update at Multiple Energies Full Author List: HAWC Collaboration R.Alfaro 1,C.Alvarez 2,A.Andrés 3,E.Anita-Rangel 3,M.Araya 4,J.C.Arteaga-Velázquez 5,D.AvilaRojas 3,H.A.AyalaSolares 6, R...

    2018

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.