{"id":"c179b48e-4532-45c1-a4f8-3008ab243580","arxiv_id":"2507.07070","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"HAWC and IceCube data combine to produce all-sky cosmic-ray anisotropy maps from 0.6 to 280 TV, extending the energy range and supporting a rigidity-based interpretation of the changing anisotropy.","lead":"This paper combines eight years of HAWC data with twelve years of IceCube data to map the arrival directions of cosmic rays over 93 percent of the sky. It reports that the anisotropy pattern changes with energy and likely depends on magnetic rigidity rather than on energy alone.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The combined maps above 10 TeV require the rigidity-matching assumption to be validated, but the paper never reports whether overlap compatibility improves after matching; if the IceCube energy cuts are tuned to absorb composition-model differences, the 40–76 TV phase transition could be an…","rationale":"The reader's weakest assumption identifies exactly the load-bearing point: the comparability of IceCube and HAWC data above 10 TeV rests on rigidity-dependent anisotropy and on the GSF composition model used to choose IceCube cuts. My read of the manuscript confirms this and sharpens it. The matching procedure is not merely a calibration detail; it is the step that defines which IceCube events enter each combined map, and the paper reports no direct evidence that the matched bins are compatible in the overlap region. The stated incompatibility at 126 TeV makes this concern concrete. A composition-model cross-check and an overlap chi2 before/after matching would settle whether the rigidity-matching creates a real physical match or an artifact. Because the paper is explicitly preliminary and the reader already conditions acceptance on such checks, my stress-test does not move the verdict.","tokens_in":14044,"tokens_out":4847,"duration_ms":61675,"concrete_test":"Rerun the KS-based IceCube energy-cut scan and the subsequent map combination twice: once with the H3a composition model and once with an extreme composition assumption (pure protons and pure iron), keeping all HAWC bins fixed. If the median IceCube energies in Table 1 shift by more than roughly 0.1 in log10(E/GeV), or if the combined dipole and quadrupole coefficients change by more than the statistical errors, the rigidity-matched result is not robust. Additionally, report the overlap-region chi2 between IceCube and HAWC relative-intensity projections before and after the matched cuts; if chi2 does not drop significantly, the rigidity-dependence assumption is not supported by the data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's \"first all-sky anisotropy study above 10 TeV\" and its phase-transition claim depend on IceCube and HAWC being combined at the same physical rigidity. Section 5 states the only direct compatibility test: in the overlapping FoV, the relative-intensity projections at 10 and 126 TeV are \"not statistically compatible\" (Fig. 1). The response is to assume rigidity dependence, slide IceCube's energy cuts to match the rigidity distributions of HAWC bins under the GSF composition model (Fig. 2, right), and then combine. This is circular in effect: the IceCube energy cut is selected precisely to make the rigidity distributions agree, and no post-matching chi2 or KS value is reported to demonstrate that the measured anisotropy maps now agree. The key unvalidated link is therefore not the GSF model by itself but the claim that the selected cuts make the two experiments observe the same physical anisotropy. If the true angular distribution is partly energy-dependent, or if GSF is wrong for the relevant mass groups, the IceCube events entering a given combined bin are at the wrong rigidity, and the combined maps and the reported rapid phase transition between 40 and 76 TV are artifacts of the matching procedure rather than measurements.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":14252,"tokens_out":5601,"duration_ms":60264,"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":[{"comment":"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.","section":"Sec. 5, Fig. 2, Table 1"},{"comment":"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.","section":"Sec. 5.1, Fig. 5"},{"comment":"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.","section":"Sec. 5, Fig. 1"},{"comment":"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.","section":"Sec. 5, Fig. 3"},{"comment":"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.","section":"Sec. 6, Fig. 5"}],"minor_comments":[{"comment":"The caption says 'GST composition model' but the text and Fig. 2 refer to the GSF model; please make the notation consistent.","section":"Table 1 caption"},{"comment":"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.","section":"Fig. 5"},{"comment":"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.","section":"Sec. 2"},{"comment":"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.'","section":"Abstract and Sec. 6"}],"recommendation":"major_revision","confidential_remarks":"This is a preliminary ICRC proceedings contribution, and the authors are appropriately cautious in calling the results preliminary. The main issue is that the central high-energy claim depends on a rigidity-matching procedure that is not validated against the actual measured anisotropy in the overlap region. If the post-matching compatibility check fails, the phase-transition claim would collapse, so I recommend requiring that check and a deconvolved or explicitly pseudo-APS treatment before acceptance. The 'first all-sky above 10 TeV' claim is also strong for a proceedings paper and should be carefully worded."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Have you seen this? It's the first all-sky anisotropy combination from HAWC and IceCube above 10 TeV, using 8 years of HAWC and 12 years of IceCube. The HAWC-only maps extend to 1 PeV and confirm the known energy-dependent anisotropy; the genuinely new piece is seven rigidity-matched combined maps and their pseudo-angular power spectra up to ~280 TV.\n\nThe authors are transparent about the method: they assume the angular distribution is rigidity-dependent, then slide IceCube energy cuts to match HAWC rigidity distributions under the GSF composition model, using a KS test to pick the best bins. The phase transition near 40-76 TV is consistent with earlier partial-sky work, so that's not a new discovery.\n\nThe soft spots are real but manageable for a preliminary proceedings. The overlap region at 126 TeV is statistically incompatible before matching; the text says 10 TeV is compatible, so that's a minor discrepancy in the reader's summary. After matching, they never show a post-matching chi2 or KS value on the actual anisotropy maps. That's the key missing number. Second, only statistical uncertainties are shown; no systematic band from composition model choice, pressure corrections, or reconstruction. Third, the pseudo-APS is not deconvolved, so the claim that 93% coverage 'largely eliminates biases' is stronger than what a pseudo-APS can demonstrate. And the procedure is circular in effect: you choose cuts to make rigidity distributions agree, then interpret agreement as evidence for rigidity dependence. Plausible, but this paper doesn't close the loop.\n\nThe writing is clear and the authors explicitly flag the need for systematic studies, including composition-model uncertainties. For an ICRC proceedings, this is a reasonable contribution. For a journal, it would need the post-matching compatibility test, a systematic band, and a deconvolved APS. I'd send it to a referee if it crossed my desk, because the all-sky combination above 10 TeV is important enough to merit careful scrutiny, but I'd expect heavy revision. I'd bring it to reading group to discuss the circularity, and I'd cite the HAWC maps if I needed an updated anisotropy measurement.","headline":"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.","tokens_in":14880,"tokens_out":4805,"would_cite":true,"duration_ms":46951,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["cosmic-ray anisotropy","all-sky map","angular power spectrum","HAWC","IceCube","rigidity dependence","energy-dependent anisotropy","TeV cosmic rays"],"falsifier":"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.","tokens_in":13799,"feed_emoji":"🌌","tokens_out":4966,"duration_ms":48354,"temperature":0.7,"pith_summary":"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.","feed_headline":"Cosmic-ray sky map, 93% full, shows a shift at 40–76 TV","feed_subtitle":"Merging HAWC and IceCube data reveals how the arrival-direction pattern reorganizes across eleven energy bins.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the method for combining multiple observatories with overlapping exposure, which the paper applies to HAWC and IceCube.","marker":"[3]"},{"why":"Provides the maximum-likelihood reconstruction that estimates the isotropic expectation and the combined relative-intensity maps without relying on detector simulations.","marker":"[13]"},{"why":"Supplies the published IceCube 12-year anisotropy data that is combined with HAWC and the previous partial-sky measurement whose energy dependence is extended.","marker":"[1]"},{"why":"Establishes that partial sky coverage creates large correlations among the spherical-harmonic coefficients, motivating the full-sky combination.","marker":"[2]"},{"why":"Provides the GSF composition model used to convert HAWC energy cuts to IceCube rigidity cuts for the bin-matching scan.","marker":"[16]"},{"why":"Supplies the generalized Li–Ma method used to compute the statistical significance of anisotropy features in the reconstructed maps.","marker":"[14]"},{"why":"Provides the H3a composition model used as a comparison in the rigidity-matching figure, helping to show the sensitivity of the matching to the composition choice.","marker":"[15]"}],"fun_headline_variants":["Full-sky cosmic-ray view reveals anisotropy shift at 40–76 TV","Combined HAWC-IceCube maps expose cosmic-ray pattern change","93% sky coverage recalibrates cosmic-ray anisotropy picture","Cosmic-ray anisotropy flips at 40–76 TV in full-sky map","HAWC + IceCube join to map cosmic-ray sky with 93% coverage"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Full-sky cosmic-ray view reveals anisotropy shift at 40–76 TV","Combined HAWC-IceCube maps expose cosmic-ray pattern change","93% sky coverage recalibrates cosmic-ray anisotropy picture","Cosmic-ray anisotropy flips at 40–76 TV in full-sky map","HAWC + IceCube join to map cosmic-ray sky with 93% coverage"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00022,"raw_usage":{"total_tokens":1408,"prompt_tokens":870,"completion_tokens":538,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":486,"completion_tokens_details":{"reasoning_tokens":440}},"tokens_in":486,"tokens_out":538,"duration_ms":5245,"temperature":1.0,"reasoning_tokens":440,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:47:37.536368+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Ahlers,ApJL886 no","cited_arxiv_id":null,"evidence_quote":"Supplies the method for combining multiple observatories with overlapping exposure, which the paper applies to HAWC and IceCube."},{"cited_title":"Ahlerset al., Astrophys","cited_arxiv_id":null,"evidence_quote":"Provides the maximum-likelihood reconstruction that estimates the isotropic expectation and the combined relative-intensity maps without relying on detector simulations."},{"cited_title":"Abbasi, M","cited_arxiv_id":null,"evidence_quote":"Supplies the published IceCube 12-year anisotropy data that is combined with HAWC and the previous partial-sky measurement whose energy dependence is extended."},{"cited_title":"Sommers,Astropart","cited_arxiv_id":null,"evidence_quote":"Establishes that partial sky coverage creates large correlations among the spherical-harmonic coefficients, motivating the full-sky combination."},{"cited_title":"Dembinski, R","cited_arxiv_id":null,"evidence_quote":"Provides the GSF composition model used to convert HAWC energy cuts to IceCube rigidity cuts for the bin-matching scan."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the generalized Li–Ma method used to compute the statistical significance of anisotropy features in the reconstructed maps."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the H3a composition model used as a comparison in the rigidity-matching figure, helping to show the sensitivity of the matching to the composition choice."}],"review_version":1}