{"id":"0b2f00ab-b2de-491c-9daf-69c86d29314c","arxiv_id":"2501.14267","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":13,"one_line_summary":"A nearby cosmic-ray source would make the gamma-ray spectral indices of molecular clouds evolve differently with energy, providing a new observational test.","lead":"A nearby cosmic-ray source proposed to explain the TeV bump in local spectra would make the gamma-ray spectral indices of nearby molecular clouds depend on each cloud's distance to the source. The paper predicts these differences are observable with LHAASO and future observatories, providing a direct test of the nearby-source hypothesis.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The proposed GMC test is unique only if the ambient Galactic CR sea is uniform over the ~1 kpc volume spanned by the ten clouds; the paper states this assumption without quantification, so spatial background fluctuations could mimic or mask the nearby-source signature.","rationale":"The reader's weakest assumption was that GMCs are passive targets with a uniform background CR sea plus a single source. I agree, and I sharpen the concern to the 'universal background' component, which is the basis for the test's uniqueness. The paper's qualitative logic is internally consistent: a single young source produces a distance-dependent CR enhancement that shifts the gamma-ray index minimum to lower energies for nearby clouds. The ad hoc error rescaling and missing uncertainty bands are real weaknesses, but they weaken the fit's credibility rather than the uniqueness of the proposed observable. The free-penetration assumption is standard and likely valid at the >10 GeV energies relevant here, so I do not treat it as the primary risk. The extended-source sensitivity issue is acknowledged by the authors and affects detectability rather than the physical interpretation. Thus the uniform-background assumption is the most load-bearing: if it fails, the proposed test loses its discriminatory power. A dedicated propagation simulation can settle whether ambient CR gradients over 1 kpc are negligible compared to the predicted nearby-source signal. Because this is a concrete, testable assumption and the paper otherwise presents a clear model prediction, the CONDITIONAL verdict remains appropriate; my analysis does not move the verdict.","tokens_in":11946,"tokens_out":16323,"duration_ms":168433,"concrete_test":"Run a 3D CR propagation code (GALPROP or DRAGON2) with a canonical supernova source distribution and no extra nearby source, using diffusion parameters consistent with the B/C ratio, and compute proton spectra at the positions of the ten GMCs in Table II. Convert these spectra to gamma-ray spectral-index curves using Eq. (5). If the spread in the minimum index energy among the ten clouds is less than the factor-of-ten shift predicted for Perseus/Taurus (all minima above ~1 TeV), the uniform-background assumption is safe; if the spread is comparable to or exceeds that shift, the proposed test is degenerate and cannot uniquely confirm the nearby-source hypothesis.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central discriminator of the paper is the energy of the gamma-ray spectral-index minimum: below 100 GeV for Perseus/Taurus and above 1 TeV for the other clouds (Fig. 4, lower-right). This mapping rests on the assumption, stated in Sec. III near Eq. (5), that the CR flux at each GMC is the sum of a single burst-like nearby source and a universal background sea that is identical for all clouds. If the ambient CR sea is not uniform over the 100-800 pc distances spanned by the selected GMCs, then a distant cloud could develop its own low-energy index minimum (e.g., from an older local source, a gradient across the Local Bubble, or spiral-arm structure), while a cloud near the fitted source could be diluted by an unrelated background gradient. In either case, the observed pattern of index minima would no longer map uniquely onto the single nearby-source hypothesis. The paper provides no quantitative estimate of background CR spectral fluctuations on these spatial scales, so the claimed one-to-one discrimination between 'nearby source' and 'widespread bump' is not fully secured. This is the most load-bearing concern because the entire test is an interpretation of spatial non-uniformity: if the background itself is non-uniform, the test's conclusion is ambiguous.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes using gamma-ray observations of nearby giant molecular clouds (GMCs) as a spatial probe of a hypothesized nearby cosmic-ray (CR) source. The authors first fit a burst-like, point-source diffusion model to the local proton/helium spectra and dipole anisotropy data (Section II), obtaining best-fit parameters for fixed source distances rs = 100, 250, and 400 pc (Table I). They then calculate the CR flux at ten GMCs within ~1 kpc as the sum of a universal background sea and the source contribution (Section III, Eq. 5), and use a pion-decay model to predict each GMC's gamma-ray spectrum (Section IV). The central claim is that the energy dependence of the gamma-ray spectral index distinguishes GMCs close to the source (Perseus, Taurus), whose index minimum falls below 100 GeV, from distant GMCs, whose minimum lies above 1 TeV (Fig. 4); if the TeV bump is a widespread Galactic phenomenon, all GMCs should show a uniform index evolution. The authors argue this provides a direct, falsifiable test of the nearby-source scenario with LHAASO, CTA, SWGO, and future space detectors.","tokens_in":12278,"tokens_out":3364,"duration_ms":34715,"significance":"If the proposed test is robust, it would convert the indirect, local evidence for a nearby CR source into a spatially resolved and falsifiable prediction, a genuinely valuable step. The paper's specific new element is the use of the energy dependence of the gamma-ray spectral index, which is insensitive to the uncertain A-factor normalization of each GMC and therefore isolates the CR spectral shape at each cloud. The model reproduces the local CR and anisotropy data with chi2/dof = 259/270, and the qualitative separation between nearby and distant GMCs is clearly illustrated. The prediction is concrete and testable with forthcoming instruments. However, the test's validity rests on assumptions that are not fully quantified, in particular the uniformity of the background CR sea over the ~1 kpc volume and the treatment of extended-source sensitivities; these are the main weaknesses. The paper is clearly written and the derivation follows standard diffusion and gamma-ray production formalism, which makes the presented predictions easy to scrutinize and reproduce.","major_comments":[{"comment":"The paper's central discriminator is the energy of the gamma-ray spectral-index minimum for different GMCs (below 100 GeV for Perseus/Taurus, above 1 TeV for the others). This mapping assumes that the background CR sea is identical at all GMC locations, as stated near Eq. (5): \"the background CR flux remains consistent across all GMCs.\" The paper provides no quantitative estimate of how much the background CR spectrum can vary over the 100-800 pc scales spanned by the selected clouds. Spatial fluctuations in the ambient CR sea (e.g., gradients across the Local Bubble, older local sources, or spiral-arm structure) could produce a low-energy index minimum in a distant cloud or mask the minimum in a nearby cloud, breaking the claimed one-to-one correspondence between the observed index pattern and the single nearby-source hypothesis. Please provide an estimate or an upper limit on background spectral fluctuations on these scales, or incorporate a spatially varying background model, before the test can be considered decisive.","section":"Section III, Eq. (5) and Section IV, Fig. 4"},{"comment":"The detectability statements compare the predicted GMC flux with point-source sensitivity curves from Ref. [41], but the text itself notes that nearby GMCs have significant extensions (~1 degree) and that the sensitivity is reduced by a factor sqrt(1 + (theta/sigma_PSF(E))^2). This correction is not applied to any of the sensitivity curves or to the claim that \"many of these chosen GMCs are likely to be detectable.\" For clouds such as Taurus and Orion A, which are emphasized as the best LHAASO targets, please quantify the expected reduction in sensitivity at the relevant energies and state whether the proposed spectral-index measurement remains feasible after this correction.","section":"Section IV, top panel of Fig. 4"},{"comment":"The predicted gamma-ray fluxes and spectral-index curves in Fig. 4 are shown as single lines without any uncertainty bands. The best-fit model parameters in Table I have statistical uncertainties (e.g., D0 of 2.73 +/- 0.20 x 10^26 cm^2 s^-1, ts of 7.2 +/- 0.5 x 10^5 yr), and the text mentions that the GMC A factors have ~30% uncertainties. The key claim that Perseus and Taurus show an earlier hardening with a minimum below 100 GeV, while all other GMCs reach their minimum above 1 TeV, should be accompanied by a propagation of these uncertainties. Without such bands, it is unclear whether the predicted separation between the two groups is statistically significant, or whether the index-minimum energies could overlap within 1 sigma.","section":"Section IV, Fig. 4 and Table I"}],"minor_comments":[{"comment":"In the sentence about LHAASO, \"future high-energy high-energy gamma-ray detectors\" contains a duplicated \"high-energy\". The same paragraph also ends with \"leading to diverse observed gamma-ray.\", which is missing a noun such as \"spectra\".","section":"Section I"},{"comment":"The sentence \"We maintain rs fixed and and explore three cases\" contains a duplicated \"and\".","section":"Section II C"},{"comment":"The phrase \"This distinct is evident\" should read \"This distinction is evident\".","section":"Section IV"},{"comment":"The manual rescaling of DAMPE helium and p+He energies (delta = 1.037 and 1.029) and the inflation of anisotropy error bars to 35% and 25 degrees are ad hoc and are not tested for their influence on the best-fit parameters. Since the fitted source parameters directly determine the GMC predictions, a short robustness check (e.g., varying delta or the error rescaling factors within reasonable ranges) would strengthen confidence in the results.","section":"Section II A and Appendix A"},{"comment":"Reference [21] is listed as \"arxiv eprint (2023)\"; please provide the journal or arXiv identifier in the standard format used by other references.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is a well-motivated and clearly presented proposal for a falsifiable test of the nearby CR source hypothesis. The main concern is not circularity (the gamma-ray predictions are genuinely external to the fit), but rather the unquantified assumption of a spatially uniform background CR sea, which is load-bearing for the spectral-index mapping. The detectability claims also need to be corrected for extended-source losses. These are fixable with additional modeling and uncertainty propagation, so I recommend major revision rather than rejection. Please also ask the authors to be clear that their conclusions about the background being a 'widespread Galactic phenomenon' assume that the background CR sea is identical across all ten GMCs; this point should be stated more carefully in the abstract and conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid model-based prediction paper. The genuinely new piece is the gamma-ray spectral-index diagnostic: if the TeV bump comes from a nearby burst-like source, GMCs close to the source should show the minimum of their gamma-ray spectral index at lower energies than distant clouds; if the bump is a global Galactic feature, all clouds should look the same. That contrast is clean and worth testing.\n\nWhat the paper does well: it fits a simple burst-like source to CR proton/helium spectra plus dipole anisotropy (chi2/dof 259/270), then propagates the best-fit CR flux to ten local GMCs. The proposed observable, the energy dependence of the gamma-ray spectral index, is insensitive to the poorly known GMC mass/distance factor A, which is a smart choice. The predictions are genuinely external — the fit does not use gamma-ray data — so the test is not circular.\n\nSoft spots, in order of importance. First, the test rests on the assumption that the background CR sea is the same at every cloud. The paper states this but does not quantify how large spatial fluctuations in the background would have to be to mimic or erase the signature. Given that the clouds span roughly 100-800 pc, this is a real caveat, though not a disqualifying one; it defines the regime where the test is clean. Second, the detectability estimates compare predicted fluxes to point-source sensitivities, while the authors themselves note the clouds are extended (about one degree) and suffer a sensitivity reduction; the curves in Fig. 4 may be optimistic. Third, there are no uncertainty bands on the predicted spectral index curves, so it is hard to see how significant the predicted differences are given parameter errors. Fourth, the anisotropy error-bar rescaling (35% amplitude, 25 degrees phase) is ad hoc, though the appendix explains a defensible rationale.\n\nNone of these sink the central argument. The paper is a predictive framework, not a measurement claim. It deserves a serious referee and, if the caveats are addressed (quantifying background non-uniformity, folding in extended-source losses, propagating errors), it will be a useful reference for LHAASO/CTA/SWGO analyses.\n\nRecommendation: send it to peer review. I would engage with it.","headline":"A clean, genuinely predictive gamma-ray test of the nearby-source explanation of the TeV CR bump, with caveats about background uniformity and detectability that a good referee can push on.","tokens_in":12794,"tokens_out":3031,"would_cite":true,"duration_ms":27975,"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":"A single nearby cosmic-ray source can be confirmed or ruled out by the energy-dependent gamma-ray spectra of nearby molecular clouds: close clouds should harden early, distant clouds late, while a Galaxy-wide bump would make all clouds…","keywords":["cosmic-ray spectral bump","nearby cosmic-ray source","giant molecular clouds","gamma-ray spectral index","dipole anisotropy","cosmic-ray diffusion","LHAASO","TeV bump"],"falsifier":"Measure the gamma-ray spectral index as a function of energy for the ten selected clouds, especially Taurus and Orion A, with LHAASO or CTA. If Taurus and Orion A both reach their hardest index at nearly the same energy, or both only above 1 TeV, the predicted distance ordering fails; if all ten clouds show indistinguishable index-versus-energy curves, the nearby-source explanation for the TeV bump and dipole anisotropy is contradicted.","tokens_in":11734,"feed_emoji":"🔭","tokens_out":8763,"duration_ms":73032,"temperature":0.7,"pith_summary":"This paper argues that the observed \"TeV bump\" in cosmic-ray proton and helium spectra, together with the amplitude and phase evolution of the cosmic-ray dipole anisotropy, can be explained by a single nearby cosmic-ray source, and that nearby giant molecular clouds (GMCs) can serve as natural detectors for that source through their gamma-ray emission. The central claim is that the energy dependence of each cloud's gamma-ray spectral index depends on the cloud's distance from the source: the closest clouds, Perseus and Taurus, should show an early hardening with a minimum index below 100 GeV, while more distant clouds should not reach their minimum until above 1 TeV. If the bump is instead a Galaxy-wide phenomenon, all clouds should show the same index-versus-energy behavior. A sympathetic reader would care because this converts an indirect, local measurement into a spatially resolved prediction that LHAASO and upcoming gamma-ray observatories can test.","feed_headline":"Nearby cosmic-ray source leaves a fingerprint in molecular clouds","feed_subtitle":"Close clouds like Perseus and Taurus should harden earlier; identical curves would rule the local source out.","key_machinery":"The central object is the energy-dependent gamma-ray spectral index, $\\Gamma(E_\\gamma) = -d\\log F_\\gamma / d\\log E_\\gamma$, computed from the predicted gamma-ray flux of each cloud. The flux is built from the proton intensity at the cloud location, which is the sum of a universal background sea and the contribution of a point-like burst source propagated with the spherically symmetric diffusion solution, then convolved with the proton-proton gamma-ray production cross section. The index curve does the argument's work because it is independent of the cloud mass-to-distance factor $A$ that sets the absolute flux level, while still carrying the distance-to-source information through the diffusion suppression factor $\\exp(-r_s^2/(4Dt))$: closer clouds feel the source's contribution at lower energies, and that shift is what should appear as an earlier minimum in the index.","core_discovery":"In the paper's own framing, the discovery is a predicted observable signature. A burst-like nearby source with a cutoff near tens of TeV, fitted to the AMS-02, DAMPE, and GRAPES-3 proton and helium spectra and to the dipole anisotropy amplitude and phase, produces gamma-ray spectra from ten GMCs within 1 kpc whose spectral-index curves are ordered by the clouds' distances to the source. For a source distance around 250 pc, Perseus and Taurus receive enough flux from the source to develop a pronounced bump and an early hardening, reaching their hardest spectral index below 100 GeV, whereas clouds such as Orion A, Cepheus, and Mon R2, being farther from the source, keep their index minimum above 1 TeV. Because the absolute gamma-ray flux of each cloud depends on the uncertain factor $A = M/d^2$, the paper identifies the energy dependence of the spectral index, not the flux normalization, as the feature that carries the information about the nearby source. The same calculation shows that for a source at 100 pc only Taurus is strongly affected, while for 400 pc the pattern is similar to 250 pc, so which clouds deviate can also help determine the source distance.","pith_inferences":["The same spectral-index technique could be applied to gamma rays from the diffuse interstellar medium around the source direction, giving a continuous two-dimensional map of the source's diffusion footprint instead of a handful of cloud samples.","If future observations find no distance ordering, the near-source hypothesis would not immediately die: magnetic shielding inside clouds or a non-uniform background cosmic-ray sea within 1 kpc could dilute the predicted signal, so cloud transport physics would need to be checked before concluding the source is absent.","The offset between close and distant clouds' minimum-index energies depends on the diffusion coefficient and source age, so precise measurements could yield an independent local measurement of the diffusion coefficient.","Applying the same framework to electrons and positrons from the same source would predict cloud-dependent gamma-ray or synchrotron emission that could be cross-checked with multi-wavelength observations."],"forward_implications":["LHAASO should detect a measurable difference between Taurus and Orion A: Taurus's gamma-ray spectrum should show a pronounced bump and an early hardening, while Orion A's spectrum should track the background until much higher energies.","CTA's energy resolution and full-sky coverage should allow the minimum-index energy to be measured for several clouds; a clean split between close clouds hardening below 100 GeV and distant clouds hardening above 1 TeV would support the nearby-source scenario.","If all selected clouds show the same index-versus-energy behavior, the TeV bump would be a Galaxy-wide feature and the nearby-source explanation of the local cosmic-ray data would be ruled out.","The pattern of which clouds deviate, with only Taurus affected for a 100 pc source and Perseus and Taurus affected for 250 or 400 pc, can break degeneracies in determining the source distance.","The predicted fluxes place many of the ten GMCs within reach of current and planned instruments, so the test can be carried out within realistic exposure times."],"supporting_citations":[{"why":"DAMPE proton spectrum: supplies the observed softening near 10 TV that defines the TeV bump.","marker":"[8]"},{"why":"DAMPE helium spectrum: provides the corresponding helium measurement and spectral softening used in the fit.","marker":"[9]"},{"why":"Ahlers and Mertsch compilation: supplies the dipole anisotropy amplitude and phase data and the smooth background anisotropy parameterization.","marker":"[13]"},{"why":"Previous nearby-source models: establish that a nearby source can naturally explain the TeV bump and anisotropy co-evolution, motivating this paper's source prescription.","marker":"[15–17]"},{"why":"Aharonian et al.: source of the GMC A factors and nuclear enhancement factor, and the earlier uniform-cosmic-ray-flux approach this work extends.","marker":"[19]"},{"why":"Albert et al.: provides additional GMC A factors and gamma-ray observation context for clouds within 1 kpc.","marker":"[20]"},{"why":"Lv et al.: previous fit of the nearby source to cosmic-ray spectra, including the DAMPE helium energy-scale rescaling and the proton-to-helium injection ratio assumed here.","marker":"[23]"},{"why":"Zucker et al.: Gaia-based distances and positions of the clouds used to compute each cloud's distance from the fitted source direction.","marker":"[38]"},{"why":"Kafexhiu et al.: differential proton-proton to gamma-ray cross section used in Eq. (5) to convert cosmic-ray flux into gamma-ray flux.","marker":"[40]"},{"why":"Peron and Aharonian: instrument sensitivity curves for LHAASO, SWGO, CTA, and space-borne detectors used to judge detectability.","marker":"[41]"}],"fun_headline_variants":["Cosmic-ray source nearby? Check cloud spectral index curves","Molecular clouds reveal if cosmic-ray bump is local or Galactic","Gamma-ray cloud spectra could fingerprint a nearby cosmic-ray source","Taurus and Perseus hold the key to a cosmic-ray source test","Spectral index in clouds: a decisive probe for a nearby CR source"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that each giant molecular cloud is a passive target: cosmic rays penetrate freely, with no magnetic shielding or internal gradients, and the proton flux at a cloud is exactly the universal background plus the diffusion-delayed contribution of the single nearby source, with no other local sources interfering.","fun_headline_variants_meta":{"raw":{"variants":["Cosmic-ray source nearby? Check cloud spectral index curves","Molecular clouds reveal if cosmic-ray bump is local or Galactic","Gamma-ray cloud spectra could fingerprint a nearby cosmic-ray source","Taurus and Perseus hold the key to a cosmic-ray source test","Spectral index in clouds: a decisive probe for a nearby CR source"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000218,"raw_usage":{"total_tokens":1467,"prompt_tokens":997,"completion_tokens":470,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":613,"completion_tokens_details":{"reasoning_tokens":382}},"tokens_in":613,"tokens_out":470,"duration_ms":4216,"temperature":1.0,"reasoning_tokens":382,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T15:14:52.274159+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the gamma-ray spectral index as a function of energy for the ten selected clouds, especially Taurus and Orion A, with LHAASO or CTA. If Taurus and Orion A both reach their hardest index at nearly the same energy, or both only above 1 TeV, the predicted distance ordering fails; if all ten clouds show indistinguishable index-versus-energy curves, the nearby-source explanation for the TeV bump and dipole anisotropy is contradicted.","supporting_citations":[{"cited_title":"Implications on the origin of cosmic rays in light of 10 TV spectral softenings","cited_arxiv_id":"1909.12857","evidence_quote":"Albert et al.: provides additional GMC A factors and gamma-ray observation context for clouds within 1 kpc."},{"cited_title":"Probing the \"Sea\" of Galactic Cosmic Rays with Fermi-LAT","cited_arxiv_id":"1811.12118","evidence_quote":"Lv et al.: previous fit of the nearby source to cosmic-ray spectra, including the DAMPE helium energy-scale rescaling and the proton-to-helium injection ratio assumed here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Zucker et al.: Gaia-based distances and positions of the clouds used to compute each cloud's distance from the fitted source direction."}],"review_version":1}