{"id":"8efb7869-1d85-4f1d-afc5-39c86b50edfb","arxiv_id":"2507.12074","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":12,"one_line_summary":"A Milky Way simulation with anisotropic cosmic-ray diffusion shows that magnetic field geometry controls where gamma-ray emission is enhanced or suppressed, producing templates that could help constrain the diffusion anisotropy.","lead":"The authors simulate how cosmic rays spread through the Milky Way when diffusion is faster along magnetic field lines than across them. They compute the resulting gamma-ray sky and show that magnetic field geometry turns diffusion anisotropy into visible spatial changes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed inner/outer gamma-ray contrast is tied to the UF23 baseline GMF geometry, which is untested against alternative field models; this model dependence is the most load-bearing unverified premise.","rationale":"The reader's weakest_assumption identifies exactly the GMF model dependence, and I agree that this is the most load-bearing issue. The paper's own explanation of the spatial contrast in Fig. 4 is phrased entirely in terms of field-line orientation relative to the Galactic disc, so the qualitative conclusion inherits every geometric assumption of the UF23 baseline model. The paper flags that JF12 and UF23 differ 'significantly in shape' and that other UF23 variants exist, but does not test them, making the limitation explicit rather than hidden. The other weaknesses noted by the reader, such as the absence of uncertainty estimates, the local-data fitting procedure, the post-hoc heavy-nuclei factor, and missing point-source masking, are real but secondary: they affect the quantitative comparison to observed fluxes, not the qualitative spatial pattern that constitutes the paper's main claim. Because the reader already set CONDITIONAL on this basis, my independent stress test does not move the verdict; the condition should be that the GMF-geometry dependence is checked against at least one alternative model before the effect is presented as a Milky Way property. The proposed concrete test, rerunning with solenoidal JF12, is the minimal computation that would settle whether the inner/outer contrast is a robust consequence of anisotropic transport or an artifact of the chosen field model.","tokens_in":6550,"tokens_out":3199,"duration_ms":40478,"concrete_test":"Rerun the CRPropa transport with the solenoidal JF12 model (and, if feasible, an alternative UF23 variant) at epsilon = 10^-1, 10^-2, and 10^-3, applying the same source-fit procedure and HERMES gas maps. Then recompute the relative gamma-ray maps of Fig. 4 and the region SEDs of Fig. 5. If the inner-Galaxy dimming / outer-Galaxy brightening pattern and the strong global suppression at epsilon = 10^-3 persist, the model-dependence concern is resolved; if the pattern reverses, disappears, or relocates, the central claim is specific to the UF23 baseline field geometry and must be restated as such.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that diffusion anisotropy strongly reshapes the all-sky gamma-ray spatial distribution rests on the UF23 baseline Galactic magnetic field model. The mechanism offered in Section 4 depends on specific UF23 field-line orientations: in the inner Galaxy, field lines pointing out of the disc let parallel-dominated diffusion escape faster and dim the emission, while in the outer Galaxy, in-plane field lines confine CRs and brighten the emission. Anisotropic diffusion follows the local field direction, so any GMF model with different vertical structure, pitch-angle distribution, or halo transition could reduce, reverse, or relocate this contrast. The authors state in Section 2 that JF12 and UF23 'show significant differences in the shape' and that other UF23 variants exist, but all transport runs in Section 3 and all maps in Figs. 3-5 use only the UF23 baseline. The paper explicitly defers testing other GMF models to future work, yet the abstract and summary present the spatial effect as a general feature of anisotropic CR transport in the Milky Way. The quantitative comparison to observations (Fig. 5) also relies on a post-hoc factor-of-2 allowance for heavier nuclei and does not apply point-source masking, but those issues affect absolute normalization rather than the geometric pattern itself. Therefore, the single most load-bearing unverified premise is that the UF23 baseline geometry is representative enough for the qualitative inner/outer effect to be robust.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents CRPropa simulations of anisotropic diffusive transport of cosmic-ray protons in the Milky Way, using the UF23 baseline Galactic magnetic field model and three values of the anisotropy ratio epsilon = D_perp/D_parallel (10^-1, 10^-2, 10^-3). The injection spectrum is fitted to local CR proton observations, and the resulting 3D CR distribution is used together with the HERMES line-of-sight integrator to compute all-sky gamma-ray maps at 1 TeV and spectral energy distributions in four Galactic-plane regions. The authors report that decreasing epsilon dims the inner Galaxy and brightens the outer Galaxy, that epsilon = 10^-3 globally suppresses gamma-ray production, and that the SEDs are broadly compatible with TibetAS-gamma, ARGO, and LHAASO data.","tokens_in":6812,"tokens_out":5082,"duration_ms":58068,"significance":"If robust, the result would provide a useful demonstration that all-sky hadronic gamma-ray maps can, in principle, constrain both the diffusion anisotropy and the large-scale GMF geometry, complementing local CR spectral fits. The use of publicly available CRPropa and HERMES, the explicit scan over epsilon, and the clear side-by-side maps are strengths. However, the central claim currently rests on a single GMF model, and the simulated maps and SEDs are presented without uncertainty estimates, so the quantitative generality of the result is not yet established.","major_comments":[{"comment":"The central spatial effect (inner Galaxy dimmer, outer Galaxy brighter for lower epsilon) is explained in Section 4 by the orientation of field lines in the UF23 baseline model. The authors state in Section 2 that JF12 and UF23 \"show significant differences in the shape\" and explicitly defer tests of other GMF models to future work, but all transport runs in Section 3 and all maps in Figs. 3-5 use only the UF23 baseline. Since anisotropic diffusion follows the local field direction, a different GMF geometry could reduce, reverse, or relocate the effect. Please run at least one alternative GMF model (e.g., the solenoidal JF12 shown in Fig. 1) for epsilon = 10^-2 and 10^-3 and compare the ratio maps, or alternatively restrict the conclusions in the abstract and summary to the UF23 baseline model.","section":"Section 2 and Section 4, Fig. 4"},{"comment":"The simulated all-sky maps and SEDs are presented without statistical or systematic uncertainties. The differences in the relative maps in Fig. 4 are of order tens of percent, while the SDE simulation has finite particle statistics and the source injection parameters are refitted for each epsilon case; part of the contrast between models could therefore reflect shot noise or normalization differences rather than the physical mechanism. Please provide an estimate of the pixel-level statistical uncertainty in the ratio maps (for example, from bootstrap resampling or from the number of simulated particles per spatial cell), or at least state the statistical precision explicitly.","section":"Section 4, Figs. 3-5"},{"comment":"The comparison with observed SEDs is only qualitative: no point-source masking is applied to the simulated maps while the observations mask sources, and a factor-of-2 allowance for heavier nuclei is invoked post hoc to absorb the normalization deficit. The statement that the observed gamma-ray spectra are \"in agreement\" with the prediction therefore overstates the constraining power of the comparison. Please either apply a source mask to the simulation maps or explicitly state that the comparison is indicative only, with the normalization uncertainty dominated by unmodelled nuclei and unmasked sources.","section":"Section 4, Fig. 5 and SED comparison"}],"minor_comments":[{"comment":"The caption contains the typo \"pannel\" instead of \"panel\".","section":"Section 2, Fig. 1"},{"comment":"The reference to Jansson & Farrar (2012) is misspelled as \"Ransson, G. Farrar\"; it should be \"Jansson, R.\" and \"Farrar, G.\".","section":"References, Ref. [9]"},{"comment":"The horizontal axis label in Fig. 5 reads \"Energy [T eV]\"; it should read \"Energy [TeV]\".","section":"Section 4, Fig. 5 axis label"},{"comment":"The sentence \"All predictions are slightly below the measured fluxes, which can be expected as the model presented here, contains only protons\" has an ungrammatical comma after \"here\"; please rephrase for clarity.","section":"Section 4, paragraph after Eq. (3)"}],"recommendation":"major_revision","confidential_remarks":"This is a clearly written proceedings-style contribution, and the authors are transparent about the single-GMF limitation in Section 2, which supports a conditional rather than a negative assessment. For a full-length journal article, however, the absence of an alternative-GMF test and the lack of uncertainties on the simulated maps would be more serious concerns. The requested revision is feasible within the manuscript's scope: adding one or two UF23 variants or the solenoidal JF12 field for the key epsilon values, and providing a rough statistical uncertainty estimate, would substantially strengthen the claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the genuinely new piece is the all-sky application of anisotropic diffusion transport to hadronic gamma-ray emission with the UF23 baseline GMF, scanning epsilon = 0.1, 0.01, 0.001. Second, the central spatial claim — inner Galaxy dims and outer Galaxy brightens as diffusion becomes more field-aligned — is clearly demonstrated in the maps, and it is tied to UF23's field geometry; that geometry dependence is both the source of the signal and the main unverified premise.\n\nWhat works: the paper is a clean, honest extension of the CMZ study, built on public tools (CRPropa, HERMES) and a fit to local CR spectra. The SEDs show the geometry-driven contrast, and the authors do not overclaim a measurement of epsilon. The paper reads as a template paper, which is a legitimate contribution.\n\nSoft spots, in order. The UF23-baseline dependence is real and should be flagged clearly. JF12 and UF23 differ significantly in shape, and the authors say they will test other models in the future; until then, the abstract-level statement that anisotropic diffusion and magnetic field geometry 'strongly influence' the all-sky emission is true only conditional on UF23. The stress-test note is right: a GMF with different vertical or pitch structure could reduce or relocate the effect. That is a load-bearing caveat, not a fatal flaw, because the mechanism is physically sensible and the maps support it for this model.\n\nNext, no statistical or systematic uncertainties are given for the maps or SEDs, and no configs/outputs are provided for exact reproduction. For an ICRC proceedings that is acceptable, but it means the quantitative statements are illustrative rather than definitive. The injection spectrum is fitted to the same local CR data that Fig. 2 reproduces; that is standard calibration rather than circularity, but it should be stated as calibration, not validation. The factor-of-2 allowance for heavier nuclei is post-hoc, and the citation token [7?] is broken.\n\nBottom line: for readers working on Galactic CR transport or diffuse gamma-ray predictions, this is worth reading and citing as a template study. It deserves a serious referee if submitted as a journal article; the referee should push for UF23-alternative GMF runs and uncertainty estimates.","headline":"Useful template paper: the inner/outer gamma-ray contrast from anisotropic diffusion is real in UF23, but its robustness across GMF models is untested and should be the main referee target.","tokens_in":7437,"tokens_out":2176,"would_cite":true,"duration_ms":26929,"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 that the all-sky gamma-ray emission from hadronic cosmic-ray interactions shifts measurably with the anisotropy of cosmic-ray diffusion, with the inner Galaxy dimming and the outer Galaxy brightening as diffusion becomes…","keywords":["anisotropic cosmic-ray diffusion","Galactic magnetic field","diffuse gamma-ray emission","cosmic-ray transport","all-sky gamma-ray maps","hadronic interactions","diffusion anisotropy ratio","Galactic cosmic rays"],"falsifier":"Run the same transport and gamma-ray calculation with an alternative large-scale field model whose geometry differs substantially; if the inner/outer brightness contrast does not change sign or position, the magnetic-field-geometry mechanism proposed here is wrong. Alternatively, measure the 1 TeV inner-to-outer Galactic disc intensity ratio after masking resolved sources and compare with the $\\epsilon = 10^{-1}$ and $\\epsilon = 10^{-2}$ predictions: a ratio that matches neither, or that varies with longitude in a way the modeled field cannot produce, would falsify the claim.","tokens_in":6278,"feed_emoji":"🌌","tokens_out":10167,"duration_ms":100715,"temperature":0.7,"pith_summary":"Galactic cosmic rays should diffuse faster along the large-scale magnetic field than across it, and this paper asks whether that anisotropy leaves a measurable imprint on the all-sky gamma-ray emission produced when cosmic rays collide with interstellar gas. The authors solve the anisotropic diffusive transport of protons in a baseline model of the Galactic magnetic field, varying the anisotropy ratio $\\epsilon = D_\\perp/D_\\parallel$ from $10^{-1}$ (nearly isotropic) to $10^{-3}$ (strongly parallel), fit the cosmic-ray source spectrum to the observed local proton flux, and then integrate the resulting three-dimensional proton flux along lines of sight to produce 1 TeV gamma-ray maps and spectral energy distributions. They find that every tested anisotropy reproduces the local proton spectrum, but the gamma-ray sky changes: the inner Galaxy becomes dimmer and the outer Galaxy brighter as diffusion becomes more parallel-dominated, and at $\\epsilon = 10^{-3}$ the global emission is suppressed because cosmic rays are essentially confined to the spiral arms. A sympathetic reader would take the paper's point to be that gamma-ray morphology, not the local cosmic-ray spectrum, is the observable that can constrain diffusion anisotropy.","feed_headline":"Parallel cosmic-ray diffusion dims inner, brightens outer Galaxy","feed_subtitle":"Varying how tightly cosmic rays follow Galactic magnetic fields visibly reshapes the 1-TeV gamma-ray sky.","key_machinery":"The central object is the anisotropic-diffusion ratio $\\epsilon = D_\\perp/D_\\parallel$, the ratio of cosmic-ray diffusion perpendicular to and along the local Galactic magnetic field. The argument carries it through a 3+1-dimensional stochastic differential equation solution of cosmic-ray diffusion, with a parallel diffusion coefficient that follows a fitted broken power law in energy and a perpendicular coefficient equal to $\\epsilon$ times that value, and then through a line-of-sight integral of the resulting proton flux against the gas distribution with a hadronic production cross section. The work this machinery does is to convert the single parameter $\\epsilon$ into a spatial pattern: because field lines leave the disc in the inner Galaxy but lie in the plane in the outer Galaxy, smaller $\\epsilon$ dims the inner Galaxy, brightens the outer Galaxy, and at $10^{-3}$ confines cosmic rays so strongly to spiral arms that the integrated emission drops.","core_discovery":"On the paper's own terms, the central discovery is that the spatial distribution of diffuse hadronic gamma-ray emission is a sensitive probe of anisotropic cosmic-ray transport. Lowering $\\epsilon$ from $10^{-1}$ to $10^{-2}$ raises the predicted gamma-ray intensity in the outer Galactic disc and lowers it in the inner Galactic disc, and lowering it to $10^{-3}$ suppresses the all-sky emission as a whole. The mechanism is the geometry of the adopted Galactic magnetic field: in the inner Galaxy field lines point out of the disc, so parallel-dominated cosmic rays escape quickly and produce fewer gamma rays, while in the outer Galaxy field lines lie inside the disc, so confinement and gamma-ray production increase. Because the same source model fits the observed local cosmic-ray spectrum for all three values of $\\epsilon$, the paper concludes that the anisotropy is invisible in the local spectrum but visible on the sky.","pith_inferences":["Inference: because the sign of the predicted contrast is set by field-line geometry, replacing the adopted large-scale field model with an alternative geometry (the paper shows two viable models with visibly different fields) could weaken, reverse, or relocate the inner/outer asymmetry, so gamma-ray maps can also discriminate among field models.","Inference: the local cosmic-ray spectrum leaves $\\epsilon$ degenerate with the source distribution, but adding secondary-to-primary ratios such as boron-to-carbon would break that degeneracy; the paper does not compute them, but its transport solution could.","Inference: if the spiral-arm confinement at $\\epsilon = 10^{-3}$ is real, gamma-ray emissivity per unit gas should be much lower between arms than inside arms; a targeted comparison along arm and inter-arm sight lines would test this directly.","Inference: since leptonic emission is neglected, the hadronic component at GeV energies could be isolated by subtracting inverse-Compton and bremsstrahlung templates from all-sky maps; the residual should then show the predicted $\\epsilon$-dependent inner/outer asymmetry."],"forward_implications":["Matching the local cosmic-ray proton spectrum cannot constrain the diffusion anisotropy, while matching gamma-ray morphology can.","For smaller $\\epsilon$, the inner-to-outer Galactic gamma-ray intensity ratio at 1 TeV should fall at fixed gas column, giving a concrete, spatially resolved prediction.","At $\\epsilon = 10^{-3}$, the model predicts a global suppression of hadronic gamma-ray brightness relative to more isotropic transport, a signature that can be compared with all-sky intensity measurements.","The spectral energy distributions in inner and outer sky regions stay close to the observed band for all tested anisotropies, so the anisotropy should be sought in maps rather than spectral slopes.","The same hadronic interactions produce neutrinos, so the predicted anisotropy-dependent spatial pattern transfers to the neutrino sky and can be checked against Galactic-plane neutrino measurements."],"supporting_citations":[{"why":"Supplies the baseline UF23 Galactic magnetic field model whose field-line geometry drives the predicted gamma-ray contrast.","marker":"[11]"},{"why":"Supplies the solenoidal JF12 field model used to illustrate how strongly alternative field geometries differ.","marker":"[10]"},{"why":"Supplies the stochastic differential equation method used to obtain steady-state cosmic-ray distributions from burst injections.","marker":"[15]"},{"why":"Supplies the pulsar-based spatial source distribution used for cosmic-ray injection in the Galaxy.","marker":"[16]"},{"why":"Supplies the fitted broken power law for the parallel diffusion coefficient used in all transport runs.","marker":"[13]"},{"why":"Supplies the line-of-sight integration framework that converts the cosmic-ray distribution into all-sky gamma-ray maps.","marker":"[25]"},{"why":"Supplies the hadronic gamma-ray production cross section used in the line-of-sight integral.","marker":"[26]"},{"why":"Supplies part of the observed high-energy proton data used to fit the source spectrum.","marker":"[24]"},{"why":"Establishes the epsilon parameter and shows the compact source region favors nearly isotropic diffusion, the starting point this paper extends to the full sky.","marker":"[8]"}],"fun_headline_variants":["Anisotropic CR transport reshapes the gamma-ray sky","Cosmic-ray anisotropy flips inner-outer Galaxy gamma-ray brightness","Lower epsilon dims inner Galaxy, brightens outer in gamma rays","Tighter field-line alignment dims inner Galaxy, brightens outer"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the adopted baseline Galactic magnetic field model (labelled UF23) has the right large-scale geometry, specifically field lines that point out of the disc in the inner Galaxy and lie in the plane in the outer Galaxy; if the real field geometry differs, the predicted inner-dimmer/outer-brighter contrast could shrink, shift, or reverse.","fun_headline_variants_meta":{"raw":{"variants":["Anisotropic CR transport reshapes the gamma-ray sky","Cosmic-ray anisotropy flips inner-outer Galaxy gamma-ray brightness","Lower epsilon dims inner Galaxy, brightens outer in gamma rays","Tighter field-line alignment dims inner Galaxy, brightens outer"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000522,"raw_usage":{"total_tokens":2564,"prompt_tokens":1026,"completion_tokens":1538,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":642,"completion_tokens_details":{"reasoning_tokens":1465}},"tokens_in":642,"tokens_out":1538,"duration_ms":13691,"temperature":1.0,"reasoning_tokens":1465,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:55:51.525445+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same transport and gamma-ray calculation with an alternative large-scale field model whose geometry differs substantially; if the inner/outer brightness contrast does not change sign or position, the magnetic-field-geometry mechanism proposed here is wrong. Alternatively, measure the 1 TeV inner-to-outer Galactic disc intensity ratio after masking resolved sources and compare with the $\\epsilon = 10^{-1}$ and $\\epsilon = 10^{-2}$ predictions: a ratio that matches neither, or that varies with longitude in a way the modeled field cannot produce, would falsify the claim.","supporting_citations":[{"cited_title":"Unger, G","cited_arxiv_id":null,"evidence_quote":"Supplies the baseline UF23 Galactic magnetic field model whose field-line geometry drives the predicted gamma-ray contrast."},{"cited_title":"Kleimann, T","cited_arxiv_id":null,"evidence_quote":"Supplies the solenoidal JF12 field model used to illustrate how strongly alternative field geometries differ."},{"cited_title":"Merten, J","cited_arxiv_id":null,"evidence_quote":"Supplies the stochastic differential equation method used to obtain steady-state cosmic-ray distributions from burst injections."},{"cited_title":"Blasi, E","cited_arxiv_id":null,"evidence_quote":"Supplies the pulsar-based spatial source distribution used for cosmic-ray injection in the Galaxy."},{"cited_title":"Dörner PhD Thesis, Ruhr-Universität Bochum (2025)","cited_arxiv_id":null,"evidence_quote":"Supplies the fitted broken power law for the parallel diffusion coefficient used in all transport runs."},{"cited_title":"Dundovic, C","cited_arxiv_id":null,"evidence_quote":"Supplies the line-of-sight integration framework that converts the cosmic-ray distribution into all-sky gamma-ray maps."},{"cited_title":"Kachelrieß, S","cited_arxiv_id":null,"evidence_quote":"Supplies the hadronic gamma-ray production cross section used in the line-of-sight integral."},{"cited_title":"Dörner, J","cited_arxiv_id":null,"evidence_quote":"Establishes the epsilon parameter and shows the compact source region favors nearly isotropic diffusion, the starting point this paper extends to the full sky."}],"review_version":1}