{"id":"44a186ee-d0a3-4a1a-886c-d767a7c92a9b","arxiv_id":"2507.10823","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Diffuse gamma-ray predictions from locally measured cosmic rays exceed LHAASO observations above 100 TeV, suggesting the PeV cosmic-ray knee originates from a local bubble roughly 0.7 to 1.5 kpc across.","lead":"Using precise new cosmic-ray measurements from LHAASO, the authors predict how much gamma-ray and neutrino light the Milky Way's gas should emit, assuming the cosmic-ray spectrum seen near Earth is the same everywhere. The prediction overshoots LHAASO's measured diffuse gamma-ray flux above 100 TeV, which they explain by a nearby 'PeV cosmic-ray bubble' whose size they estimate.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The >30 TeV excess is a robust spectral-shape mismatch, but the inferred local bubble radii are underconstrained because a spatially varying CR spectrum is not tested; a longitude-profile check is needed.","rationale":"The reader's weakest-assumption statement focuses on gas column density and the 10 GeV pion-decay normalization. However, because the model spectra are normalized at 10 GeV to the Fermi-LAT data, the gas column density cancels in the ratio of predicted flux at high energy to that at 10 GeV, so the gas-distribution uncertainty does not directly inflate the >30 TeV excess; the excess is a spectral-shape effect. The 10 GeV normalization uncertainty is bounded by the paper at 10% and is too small to erase a factor-of-several excess. The more consequential weakness is the interpretation of the excess as a local bubble: the region-averaged spectra cannot discriminate a bubble from a spatially varying CR spectrum, and the paper does not confront the bubble's predicted longitude structure with LHAASO data. This leaves the bubble radii underconstrained, but it does not invalidate the central discrepancy, so the reader's CONDITIONAL verdict remains the appropriate outcome. The proposed longitude-profile test would settle whether the bubble interpretation is genuinely required or merely one of several viable explanations.","tokens_in":13661,"tokens_out":22499,"duration_ms":272081,"concrete_test":"Compute the >100 TeV diffuse gamma-ray longitude profile predicted by the best-fit bubble model (Model A3) and by a model with an assumed smooth radial gradient of the CR spectral slope but no local bubble, using the same gas maps and the same 10 GeV Fermi-LAT normalization. Compare both profiles to the LHAASO longitude-binned diffuse flux from Ref. [54]. If the bubble model's localized bump is present in the data and the smooth-gradient model is clearly disfavored, the bubble inference is supported; if the smooth-gradient model fits at least as well, the bubble radius claim is not uniquely constrained and should be treated as illustrative rather than a firm measurement.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The claimed discrepancy in the inner Galaxy above 30 TeV and in the outer Galaxy above 100 TeV is primarily a spectral-shape effect: the diffuse spectra are normalized at 10 GeV to the Fermi-LAT diffuse flux, so the gas column density cancels in the ratio of predicted flux at high energy to that at 10 GeV. The reader's identified gas-distribution uncertainty is therefore not the most limiting factor for the existence of the excess. The genuinely load-bearing weak point is the next step: interpreting the excess as a local PeV CR bubble with radii 1.5 +/- 0.3 kpc (inner) and 0.7 +/- 0.3 kpc (outer). The two broad region-averaged spectra used in the fit cannot distinguish a local spherical bubble from a spatially varying CR spectrum, for example a position-dependent knee or escape that softens the high-energy proton spectrum in the inner and outer disk. The bubble radius is a one-zone toy fit to the very data it is meant to explain, and the paper does not test the alternative of a smooth CR spectral gradient across the Galaxy. The bubble model predicts a distinctive longitude profile of >100 TeV gamma-ray emission, with a localized enhancement around the Sun, but the paper does not compare this prediction to the LHAASO spatial data. Without such a test, the quantitative bubble extension is not uniquely determined and should be regarded as illustrative rather than measured.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs two independent models of the Galactic cosmic-ray (CR) spectrum up to the knee (a three-Peters-cycle fit, Model A, and a DRAGON-2 diffusion-break fit, Model B), both calibrated to AMS-02, DAMPE, CALET, LHAASO, IceTop and Auger data. Assuming the local CR spectrum is universal throughout the Milky Way, the authors compute diffuse gamma-ray and neutrino emission from the inner (15°<l<125°) and outer (125°<l<235°) Galactic plane, normalizing the gamma-ray prediction to Fermi-LAT at 10 GeV with a 0.85 pion-decay fraction. Comparing to LHAASO diffuse gamma-ray data, they find the predicted pion-decay flux exceeds the measured diffuse flux above ~30 TeV (inner) and ~100 TeV (outer). They interpret this excess as evidence for a \"local PeV CR bubble\" and fit its radius, obtaining 1.5±0.3 kpc toward the inner Galaxy and 0.7±0.3 kpc toward the outer Galaxy; they also compare neutrino predictions to rescaled IceCube template fluxes.","tokens_in":13961,"tokens_out":5533,"duration_ms":65862,"significance":"If the claimed discrepancy is real, it challenges the standard assumption of a universal CR spectrum in the Galaxy and suggests that the knee is produced by a local, recent CR injection. The use of two independent CR models (Peters cycles and DRAGON-2) to predict the diffuse emission is a strength, and the consistency of the two predictions supports the robustness of the spectral-shape mismatch. The predicted gamma-ray and neutrino fluxes are falsifiable and the comparison to IceCube is a useful cross-check. However, the inference of bubble radii is not on the same footing as the discrepancy itself: it relies on a one-zone toy model with no spatial test.","major_comments":[{"comment":"The bubble-radii claim (R = 1.5 ± 0.3 kpc inner, 0.7 ± 0.3 kpc outer) is based only on the region-averaged spectra of the inner and outer Galaxy. The model predicts a characteristic longitude profile of >100 TeV emission that is localized around the Sun, but the paper does not compare this prediction with LHAASO's spatial data. A position-dependent CR spectrum (e.g., a radial gradient of the knee energy or of the escape rate) could produce the same two region-averaged spectra with different \"bubble\" sizes. Without a longitude-profile test, the quoted radii are not uniquely determined; they should be presented as illustrative, not as measured extensions.","section":"Section V, Fig. 3"},{"comment":"The uncertainties on the bubble radii (0.3 kpc) are not derived from a statistical fit or from propagation of the input uncertainties (CR fit parameters, gas models, 10 GeV normalization, gamma-ray absorption). The text states only that the bubble radius was \"adjusted\" until the model \"fits\" the LHAASO data. Please specify the fitting procedure, the chi-square or likelihood used, and how the 0.3 kpc uncertainty was obtained; alternatively, state explicitly that these are systematic ranges from varying model assumptions.","section":"Section V"},{"comment":"The decomposition of the local CR spectrum into a steady-state component (Peters cycles 1 and 2) and a local knee component (Peters cycle 3) is an assumption of the toy model, not a result. The same LHAASO proton spectrum that defines the third cycle is then used as the spectrum inside the bubble, and the bubble radius is fit to the gamma-ray data. This leaves a degeneracy between the amplitude of the local component and its spatial extent. A minimal test would be to fix the bubble spectrum from an independent CR propagation model (e.g., a recent source injection) and fit only the geometry, or to include an alternative spatially varying CR model and compare fits.","section":"Sections III and IV"}],"minor_comments":[{"comment":"The caption reads \"155◦ < l <125◦\", which appears to be a typo; it should be \"15◦ < l <125◦\" consistent with the text in Section III.","section":"Section III, Fig. 3 caption"},{"comment":"There is a typo \"tot he\" in the sentence \"the CRs with the spectrum corresponding tot he 3rd Peters cycle\"; it should be \"to the\".","section":"Section V"},{"comment":"The word \"intergrals\" in \"line-of-sight intergrals\" should be \"integrals\".","section":"Section III"},{"comment":"The 0.85 pion-decay fraction at 10 GeV and its quoted 10% uncertainty are not propagated to the predicted fluxes above 100 TeV; please add a systematic band or state explicitly that the excess is insensitive to this normalization.","section":"Section III"},{"comment":"The 10% upward shift of Auger energies is described, but its effect on the fitted CR spectra (and hence on the predicted gamma-ray flux) is not quantified; a short discussion would help.","section":"Appendix A"},{"comment":"The figure captions and legends are dense; consider adding a clear legend directly in the figure distinguishing the diffuse (masked) LHAASO data, total flux, Model A1/A2/A3, Model B, and IceCube templates.","section":"Figures 3 and 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal. The central discrepancy between the CR-model predictions and LHAASO diffuse gamma-ray data is likely robust and worth publishing, but the bubble-radius inference is not supported by the present spatial analysis. The abstract and conclusion overstate the certainty of the radii. I recommend major revision to either add a longitude-profile test or substantially soften the claim about inferred extensions. No concerns about novelty or citation practice; the reference to the authors' previous work [32] is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the headline: the paper's main discrepancy—that extrapolating the locally measured CR spectrum to the whole Galaxy overpredicts LHAASO's diffuse gamma-ray flux above ~100 TeV—holds up. This is a real result, a useful confirmation and sharpening of the earlier excess seen by this group's Ref. [32]. The new LHAASO proton spectrum, plus the cross-check with two very different CR models (Peters cycles and DRAGON-2), makes the spectral-shape mismatch fairly convincing.\n\nWhere the paper is genuinely new is the attempt to turn the excess into a geometry: a 'local PeV CR bubble' with radii 1.5 ± 0.3 kpc inward and 0.7 ± 0.3 kpc outward. I agree with the stress-test note that this is the load-bearing weak point. The two region-averaged spectra cannot distinguish a local spherical bubble from a smooth spatial gradient in the CR spectrum—say, a position-dependent knee or rigidity-dependent escape that softens the high-energy proton population in the inner and outer disk. The bubble radius is a one-zone parameter fitted to the same data it explains, and the paper does not test the predicted longitude profile against LHAASO's spatial information. So treat the radii as illustrative, as the authors themselves say when they call it a toy model, but the abstract and conclusion still present them as inferred quantities.\n\nThe reader's concern about gas distributions is real but secondary. Because the predicted spectra are normalized to Fermi-LAT at 10 GeV, the gas column density largely cancels in the shape of the excess. The more important systematics are the 10 GeV pion-decay fraction and gamma-ray absorption, both bracketed at a rough 10% level. The excess is large enough that these don't kill the main claim, but they are not propagated into the bubble radii, whose systematic width is wider than the quoted statistics.\n\nOn balance, this is a serious paper. The contradiction between a universal CR spectrum and the LHAASO diffuse flux above 100 TeV is an important clue for Galactic CR physics, and the paper makes a reproducible case for it. The bubble interpretation is not established, but the paper is honest about its limitations and points to the specific future test—longitude profiles of >100 TeV flux. I'd send it to a good referee and would cite it for the confirmation of the excess, not for the bubble geometry. Recommendation: engage, but push for an explicit longitude-profile check or a clear statement that alternative spatial CR gradients are not tested.","headline":"The universal-CR-spectrum excess above 100 TeV is robust and worth citing; the inferred bubble radii are illustrative, not measured.","tokens_in":14552,"tokens_out":3849,"would_cite":true,"duration_ms":39234,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Gamma-ray data reveal a local PeV cosmic-ray bubble around the Solar System.","keywords":["cosmic-ray knee","diffuse gamma-ray emission","PeV cosmic-ray bubble","LHAASO","galactic cosmic rays","neutrino emission from the Milky Way","Peters cycles"],"falsifier":"Measure the >100 TeV diffuse gamma-ray flux in a grid of longitude bins along the inner Galactic plane and compare it with the gas-column-weighted universal-cosmic-ray prediction. If any distant, gas-rich direction shows the full universal-prediction flux, the Solar System cannot be inside a small local bubble; if the flux drops below the prediction everywhere above about 100 TeV, the local-bubble interpretation is supported.","tokens_in":13421,"feed_emoji":"🌌","tokens_out":6658,"duration_ms":82195,"temperature":0.7,"pith_summary":"The paper tries to show that the cosmic-ray spectrum measured locally with LHAASO, if assumed to be the same everywhere in the Milky Way, predicts more diffuse gamma-ray emission from the Galactic plane above 30 to 100 TeV than LHAASO actually observes. It argues that the cleanest explanation is that the highest-energy part of the local cosmic-ray spectrum is not universal: the Solar System sits inside a recent 'PeV bubble' of freshly injected cosmic rays, about 1.5 kpc in extent toward the inner Galaxy and 0.7 kpc toward the outer Galaxy. If this is right, the cosmic-ray knee is not a permanent, Galaxy-wide feature but a transient local event, and the spatial pattern of very-high-energy gamma-ray emission becomes a map of recent cosmic-ray injections.","feed_headline":"Milky Way gamma rays reveal a local PeV cosmic-ray bubble","feed_subtitle":"A universal-knee model overshoots LHAASO's >100 TeV diffuse flux; a recent nearby source fits the deficit.","key_machinery":"The central object is the local cosmic-ray spectrum decomposed into three 'Peters cycles', rigidity-scaled source populations whose spectral shape is identical for all nuclear groups, with the third, highest-energy cycle accounting for the knee region seen by LHAASO. The argument works by computing diffuse gamma-ray and neutrino yields from these spectra through pion-production calculations and gamma-ray absorption corrections, comparing the results with Fermi-LAT and LHAASO diffuse measurements, and then re-fitting the third cycle as emission from a finite spherical bubble around the Solar System. The bubble radius is the parameter that reconciles the model with the observed deficit.","core_discovery":"The paper establishes that the pion-decay component of diffuse gamma-ray emission computed from a cosmic-ray spectrum fitted to LHAASO's precise knee-region measurements exceeds the measured LHAASO diffuse flux in the inner Galaxy above about 30 TeV and in the outer Galaxy above about 100 TeV. Removing the highest-energy rigidity component of the fitted cosmic-ray spectrum from the emission calculation removes the excess, which the authors interpret as evidence that these particles occupy a limited volume around the Solar System rather than the whole Galactic disk. Attributing that highest-energy component to a local PeV cosmic-ray bubble and fitting its size to the LHAASO data gives an extension of 1.5 ± 0.3 kpc toward the inner Galaxy and 0.7 ± 0.3 kpc toward the outer Galaxy.","pith_inferences":["Extension: the same decomposition could be applied to a longitude-resolved map of 100 TeV to 1 PeV gamma-ray emission, turning each high-energy spectral component that appears locally but not globally into a tracer of a recent cosmic-ray injection event.","Extension: the local-bubble interpretation predicts a small arrival-direction anisotropy in local cosmic rays at and above the knee, roughly aligned with the inferred bubble direction, which current and future observatories could search for.","Extension: the inner-Galaxy radius of 1.5 ± 0.3 kpc may be better interpreted as a lower bound if the injected particles have not had time to fill the whole bubble; combining the radius with a diffusion model could yield an age estimate for the recent injection event."],"forward_implications":["If the local-bubble scenario is right, the cosmic-ray knee is not the same feature everywhere in the Galaxy; other regions have their own knee features set by recent local injections.","LHAASO's deeper exposure should reveal strong variation of the >100 TeV diffuse gamma-ray flux from one line of sight to another, tracking where fresh PeV cosmic rays currently reside.","Because gamma-rays and neutrinos share the same pion-production parentage, the same excess should appear as a suppression of the predicted high-energy neutrino flux from the same sky regions, testable with future neutrino telescopes.","The fitted bubble sizes of 1.5 kpc inward and 0.7 kpc outward locate the Solar System asymmetrically inside this recent injection, pointing roughly toward the inner Galaxy and giving a target for identifying the source."],"supporting_citations":[{"why":"Supplies the LHAASO proton spectrum that fixes the shape of the highest-energy cycle in the cosmic-ray fit.","marker":"[12]"},{"why":"Provides the LHAASO all-particle spectrum and mean logarithmic mass used to constrain the composition fit across the knee.","marker":"[11]"},{"why":"Supplies the LHAASO diffuse gamma-ray flux measurements from the inner and outer Galactic plane that the universal-spectrum model overshoots.","marker":"[54]"},{"why":"Earlier version of the same modeling chain that first reported the excess without the LHAASO proton data.","marker":"[32]"},{"why":"Provides the atomic and molecular hydrogen column densities and cosmic-ray density approximations used to convert spectra into line-of-sight emission.","marker":"[47]"},{"why":"Fermi-LAT GeV diffuse emission analysis used to set the overall normalization of the predicted flux at 10 GeV.","marker":"[19]"},{"why":"Determines the 0.85 pion-decay fraction of the low-latitude diffuse flux used for that normalization.","marker":"[20]"},{"why":"IceCube's Galactic neutrino detection whose template flux estimates are rescaled for comparison with the model predictions.","marker":"[25]"},{"why":"Provides the pion-production yields used to compute both the gamma-ray and neutrino spectra from cosmic-ray interactions.","marker":"[50]"}],"fun_headline_variants":["Local PeV bubble explains LHAASO's gamma-ray excess","Knee-region cosmic rays point to a nearby source","A cosmic-ray bubble around Earth fits LHAASO data","Milky Way's gamma glow reveals a local PeV bubble"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole excess calculation leans on a model of how much hydrogen gas lies along each line of sight, and the uncertainties in that gas model are not propagated; if the gas is overestimated, the predicted glow is too bright and the inferred bubble could be an artifact.","fun_headline_variants_meta":{"raw":{"variants":["Local PeV bubble explains LHAASO's gamma-ray excess","Knee-region cosmic rays point to a nearby source","A cosmic-ray bubble around Earth fits LHAASO data","Milky Way's gamma glow reveals a local PeV bubble"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000285,"raw_usage":{"total_tokens":1674,"prompt_tokens":933,"completion_tokens":741,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":549,"completion_tokens_details":{"reasoning_tokens":670}},"tokens_in":549,"tokens_out":741,"duration_ms":8672,"temperature":1.0,"reasoning_tokens":670,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:24:21.711822+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the >100 TeV diffuse gamma-ray flux in a grid of longitude bins along the inner Galactic plane and compare it with the gas-column-weighted universal-cosmic-ray prediction. If any distant, gas-rich direction shows the full universal-prediction flux, the Solar System cannot be inside a small local bubble; if the flux drops below the prediction everywhere above about 100 TeV, the local-bubble interpretation is supported.","supporting_citations":[{"cited_title":"Alemanno, Q","cited_arxiv_id":null,"evidence_quote":"Provides the LHAASO all-particle spectrum and mean logarithmic mass used to constrain the composition fit across the knee."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the atomic and molecular hydrogen column densities and cosmic-ray density approximations used to convert spectra into line-of-sight emission."}],"review_version":1}