{"id":"aa8ba8cf-722e-480c-9c6e-9fab3d19e076","arxiv_id":"2510.15634","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A cosmological simulation with on-the-fly cosmic-ray spectral tracking predicts diffuse gamma-ray emission from local galaxy clusters that lies 3-5 orders of magnitude below Fermi-LAT upper limits, implying a sensitivity of F_gamma < 10^-11 gamma s^-1 cm^-2 is needed to detect Coma.","lead":"Using a new simulation that tracks cosmic-ray protons in a realistic model of the local universe, this paper computes the faint gamma-ray glow that galaxy clusters and cosmic-web filaments should emit and finds it is 3-5 orders of magnitude below current Fermi-LAT limits. The result quantifies how little non-thermal pressure cosmic-ray protons exert in clusters and sets a concrete sensitivity target for detecting Coma.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'lower limit' label and the Coma sensitivity floor are not robust: the fixed injection momentum p_inj=0.1 overproduces gamma-ray emission for the steep-spectrum weak shocks that dominate cluster emission, so the true DSA-only flux could be far lower than quoted (or higher, if shock-finder/obliqu","rationale":"The reader's weakest assumption focused on the shock-finder resolution limit and obliquity bias, which cause underprediction. I agree those are serious systematics. However, the more load-bearing issue for the paper's own headline 'lower limit' is the fixed injection momentum p_inj=0.1, which the authors explicitly state leads to increased emission relative to physically motivated values. Because this bias is opposite in sign to the shock-finder/obliquity biases, the net direction of the systematic error is unknown, and the paper's claim to provide a lower limit is not logically supported. The qualitative result (fluxes far below Fermi-LAT limits) likely survives, but the quantitative Coma sensitivity floor and the 'lower limit' framing are conditional. This reinforces the reader's CONDITIONAL verdict rather than changing it, so I mark UNCHANGED. My disagreement with the reader is partial: they emphasize underprediction while I emphasize the mixed/unknown sign of the systematics, with p_inj as the clearest internal evidence against the lower-limit interpretation.","tokens_in":19257,"tokens_out":8732,"duration_ms":75036,"concrete_test":"Re-run the Coma cluster zoom-in at 8x resolution (as proposed in Sec. 6.2.1) with two variants: (a) fiducial fixed p_inj=0.1, and (b) per-shock p_inj from Eq. 15. Compare the resulting 0.5-200 GeV fluxes. If the flux changes by more than a factor of 10, the quoted F<10^-11 sensitivity floor and the 'lower limit' label are not robust. Alternatively, in post-processing, quantify the fraction of Coma's gamma-ray luminosity contributed by weak shocks (M_s≲4); if this fraction is >50%, the p_inj sensitivity shown in Fig. 7 propagates directly into the headline number.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's central claim, 'This provides a lower limit for diffuse emission from CR protons accelerated at structure formation shocks,' is undermined by the paper's own parameter-choice analysis. In Sec. 6.2.1, Fig. 7 shows that for the weak shocks (M_s≈2-3, spectral slope q≈4.5-5) which dominate the cluster-center CR population (Fig. 3), the gamma-ray emissivity at fixed CR energy density varies by 2-3 orders of magnitude as p_inj is varied over the physically expected range (~0.01-1). The authors state that their fixed p_inj=0.1 'leads to an increased emission compared to typical values of p_inj when computed directly from Eq. 15.' Thus the simulated fluxes are likely overestimates of pure DSA emission, contradicting the 'lower limit' interpretation. Conversely, the shock-finder resolution limits and obliquity bias (Sec. 6.1) cause underprediction. With biases of opposite sign and unknown magnitude, the net direction is uncertain. The specific Coma sensitivity floor, F_gamma<10^-11 γ s^-1 cm^-2, therefore is not a reliable detection threshold: the true required sensitivity could be substantially stricter (if p_inj effects dominate) or looser (if shock-finder corrections dominate). The qualitative conclusion that fluxes lie below current Fermi-LAT limits is probably robust, but the quantitative 'lower limit' and the exact sensitivity number are not supported by the evidence presented.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents predictions for diffuse gamma-ray emission from the local cosmic web and galaxy clusters, based on the first cosmological MHD simulation with an on-the-fly spectral cosmic-ray (CR) proton model. The simulation uses constrained initial conditions matching the local universe and follows CR injection at shocks (Ryu et al. 2019; Pais et al. 2018), adiabatic energy changes, and advection, with gamma rays computed from pion decay. The authors find CR protons accelerated at structure-formation shocks, producing diffuse gamma-ray fluxes a few orders of magnitude below current Fermi-LAT upper limits, and they state that a sensitivity of F_gamma < 10^-11 gamma s^-1 cm^-2 would be required to detect diffuse emission in Coma. The paper concludes that this provides a lower limit for diffuse emission from CR protons accelerated at structure-formation shocks.","tokens_in":19484,"tokens_out":5277,"duration_ms":41813,"significance":"If the quantitative predictions hold, this work would establish that shock-accelerated CR protons in the local universe radiate far below the current observable floor, and that the cluster-averaged CR-to-thermal pressure ratio is ~1e-4 to 1e-3, well below the ~1% Fermi-LAT bound. The paper's strengths include the first on-the-fly spectral CR treatment in a large cosmological MHD simulation, a gamma-ray emissivity calculation independently benchmarked against the minot package in App. B (agreement at the tens-of-percent level), and a direct comparison to external Fermi-LAT upper limits without fitting to the target observable. However, as discussed below, the central 'lower limit' claim and the specific Coma sensitivity threshold are not supported by the paper's own systematic-uncertainty analysis, so the significance of the exact numbers is lower than the abstract suggests.","major_comments":[{"comment":"The abstract's central claim that the results provide 'a lower limit for diffuse emission from CR protons accelerated at structure formation shocks' is contradicted by the paper's own parameter study. For weak shocks (q≈4.5–5) that dominate the cluster-center CR population (Fig. 3), the gamma-ray emissivity at fixed CR energy density varies by 2–3 orders of magnitude as p_inj is varied over the physically expected range (~0.01–1). The authors state that their fixed p_inj=0.1 'leads to an increased emission compared to typical values of p_inj when computed directly from Eq. 15.' Thus the simulated fluxes are likely overestimates, not lower limits, of pure DSA emission. The opposite bias from shock-finder resolution (Sec. 6.1) leaves the net direction uncertain. The 'lower limit' language should be revised or replaced with a statement about model-dependent prediction.","section":"Sec. 6.2.1, Fig. 7, Abstract"},{"comment":"The predicted Coma detection floor, F_gamma < 10^-11 gamma s^-1 cm^-2, is presented as a quantitative result, but the paper's own estimates of systematic uncertainties are large and of opposite sign. The shock finder misses roughly 1–2 orders of magnitude of dissipated energy at low Mach number (Sec. 6.1), and the obliquity bias suppresses proton acceleration by a factor 5–10. Conversely, the fixed injection momentum p_inj=0.1 overproduces emission for weak shocks by up to 2–3 orders of magnitude (Sec. 6.2.1). Without propagating these biases into a conservative range, the specific sensitivity number is not reliable. The authors should either derive a lower limit using the lower envelope of their parameter study or explicitly label the Coma sensitivity as an order-of-magnitude estimate with large systematic uncertainty.","section":"Sec. 6.1, Eq. (14); Sec. 4, Fig. 2"},{"comment":"The pressure ratio X_cr is reported as a lower limit in Sec. 4 (due to the ultra-relativistic approximation), while the gamma-ray emission is claimed as a lower limit in the abstract. This is internally inconsistent: the same numerical treatment that underestimates CR pressure also affects the low-momentum part of the spectrum, but the gamma-ray emissivity for steep weak-shock spectra is dominated by the choice of injection momentum, which the authors find overestimates emission relative to Eq. (15). The direction of the bias should be stated explicitly for each observable, and the 'lower limit' attribution for the gamma-ray flux should be corrected.","section":"Sec. 4 and Sec. 2.3"}],"minor_comments":[{"comment":"The flux formula seems to use Ω both as a solid angle and as a volume; please clarify the notation and define the integration domain.","section":"Eq. (7)"},{"comment":"Typo: 'scattering offthe gas' should be 'scattering off the gas'.","section":"Abstract"},{"comment":"The y-axis label 'F [0.5−200] GeV [ cm−2 s−1]' should include 'γ' before 'cm^-2' to indicate photon units.","section":"Fig. 2"},{"comment":"The piecewise definition of A_max(E_p) uses a threshold 'E_th^p' that is not explicitly defined; please define the threshold kinetic energy.","section":"App. A, Eq. (A.5)"},{"comment":"The reference for Stecker (1971) appears incomplete/odd: '1, Vol. 249, Cosmic gamma rays'. Please format it correctly.","section":"References"},{"comment":"The statement that the closed lower boundary 'mimics low-momentum cooling on adiabatic compression' is unclear; please explain how a closed boundary represents cooling.","section":"Sec. 2.3"},{"comment":"The caption does not describe what the dashed lines represent; please add a sentence explaining the contours.","section":"Fig. 7"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is technically solid and the simulation is a valuable resource, but the headline claim of a 'lower limit' is overstated relative to the paper's own error budget. The authors need to either revise the central claim to a model-dependent prediction with a quantified range, or remove the 'lower limit' language entirely. The comparison to Fermi-LAT data is appropriate and the App. B benchmark is a strong element. I recommend major revision rather than rejection because the core qualitative result (fluxes below current limits) is probably robust, but the quantitative load-bearing numbers must be reframed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, what to know: this is a genuine technical first — on-the-fly spectral CR protons in a constrained cosmological MHD simulation — and the qualitative conclusion that DSA-only cluster gamma-ray emission lies well below current Fermi-LAT limits is probably right. The specific \"lower limit\" claim and the Coma sensitivity number are not supported by the paper's own parameter study.\n\nWhat the paper does well: the Crescendo solver is a real advance. Getting actual CR spectra from shock injection, adiabatic changes, and advection in a 500 Mpc constrained box, then computing gamma-ray emissivity directly from those spectra, is new. The emissivity is cross-checked against the minot package in Appendix B with reasonable agreement. The comparison against Fermi-LAT upper limits and against earlier post-processing work (Pinzke, Vazza, Ha) is fair, and the qualitative agreement with those papers is credible. The extensive self-citations are to the simulation series and the CR model papers; they are appropriate.\n\nWhere it gets soft: the paper undercuts its own headline. The shock finder is resolution-limited, and the authors estimate an under-prediction of dissipated energy by 1-2 orders of magnitude at low Mach numbers; the obliquity efficiency model further suppresses proton acceleration by a factor 5-10 because of numerically broadened shocks. Those effects push the predicted flux down. But the fixed injection momentum p_inj=0.1 pushes the other way. The authors themselves state in Sec. 6.2.1 that this choice \"leads to an increased emission compared to typical values of p_inj\" from Eq. 15, and Fig. 7 shows a 2-3 order of magnitude increase in emissivity for the weak shocks (q~4.5-5) that dominate cluster-center CR populations. So the net bias is unknown — possibly up, possibly down, with orders of magnitude at stake. That makes the abstract's \"This provides a lower limit\" statement wrong, or at least not established. The Coma sensitivity threshold F<10^-11 ph/s/cm^2 should be read as model-dependent, not a robust detection requirement.\n\nA secondary issue: no code or data products are released. For a 30-million-CPU-hour simulation that's understandable, but it does mean the quantitative claims are a re-implementation exercise, and the propagated systematic uncertainty is not quantified anywhere in the paper.\n\nWho should read this: anyone working on nonthermal ICM physics, CR acceleration at structure formation shocks, or gamma-ray observability of clusters. It is a pathfinder with a real capability and an honest discussion of numerics, but the headline numbers need a caveat pass.\n\nMy recommendation: send it to peer review. The technical advance and the emissivity cross-check deserve referee time, and the main qualitative result is probably robust. A good referee should ask for a quantitative uncertainty budget or a softening of the \"lower limit\" language before publication.","headline":"First on-the-fly spectral CR-proton cosmological simulation: qualitative conclusion likely right, but the 'lower limit' and Coma sensitivity claims are not supported by the paper's own parameter study.","tokens_in":20204,"tokens_out":3920,"would_cite":true,"duration_ms":33094,"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":"A new simulation finds that shock-accelerated protons in local galaxy clusters produce gamma-rays far below what Fermi-LAT can detect.","keywords":["cosmic rays","gamma-ray emission","galaxy clusters","cosmological simulations","magnetohydrodynamics","cosmic web","pion decay","Fermi-LAT"],"falsifier":"Re-run the same constrained local-volume simulation at 8x resolution (as the authors plan for zoom-ins) and compare the predicted Coma gamma-ray flux; if it rises by more than a factor of ~100, the stated lower limit is not robust. Alternatively, a gamma-ray telescope reaching a sensitivity of 10^-12 gamma s^-1 cm^-2 that detects diffuse emission from Coma would directly falsify the paper's prediction that this emission lies below 10^-11.","tokens_in":18959,"feed_emoji":"🌌","tokens_out":6471,"duration_ms":51721,"temperature":0.7,"pith_summary":"This paper presents the first cosmological magnetohydrodynamic simulation that follows cosmic-ray proton spectra on the fly and uses it to predict diffuse gamma-ray emission from the local universe's galaxy clusters and cosmic web. The central finding is that protons accelerated at structure-formation shocks produce a steady glow via pion-decay, but the predicted flux lies a few orders of magnitude below current Fermi-LAT upper limits. For the Coma cluster, the model says detection would require a sensitivity below 10^-11 gamma s^-1 cm^-2, which the authors interpret as a lower limit for diffuse emission from shock-accelerated protons. If correct, cluster cosmic-ray pressures are only about 10^-4 to 10^-3 of thermal pressure, far below the roughly 1% bound, meaning non-thermal pressure support is negligible.","feed_headline":"Cluster gamma-ray glow predicted far below detection limits","feed_subtitle":"First spectral cosmic-ray simulation says Coma's diffuse glow needs sensitivity below 10^-11 photons/cm2/s.","key_machinery":"The central mechanism is a spectral cosmic-ray model coupled to the magnetohydrodynamic simulation: a Fokker-Planck solver that represents the proton distribution as piecewise power laws in momentum, injects them at shocks according to diffusive shock acceleration with an obliquity-dependent efficiency, follows adiabatic compression/expansion and advection with the gas, and closes the low-momentum boundary. Gamma-ray emission is computed directly from these spectra using a parametrized pion-decay cross section, so the prediction does not rely on post-processing assumptions beyond the CR transport model itself.","core_discovery":"Using a constrained simulation of a 500 h^-1 Mpc volume that reproduces the local large-scale structure, the authors evolve cosmic-ray protons with an on-the-fly Fokker-Planck solver, injecting them at diffusive shocks with a Mach-number- and obliquity-dependent efficiency and then advecting them with the gas. They find proton acceleration at essentially all structure-formation and accretion shocks around clusters and filaments, producing diffuse gamma-ray halos that extend to and beyond the virial radius. The cluster-averaged cosmic-ray-to-thermal pressure ratio comes out 2-3 orders of magnitude below the Fermi-LAT bound, and the corresponding gamma-ray flux and luminosity are 3-5 orders be","pith_inferences":["If the shock finder's resolution gap is closed by higher-resolution zoom-ins (as the authors propose), the missing low-Mach-number shocks could raise the predicted flux by 1-2 orders of magnitude, potentially narrowing but not eliminating the gap to Fermi-LAT limits.","Recent 3D particle-in-cell results indicating proton acceleration at quasi-perpendicular shocks could add another factor of 5-10 to the emission, still likely leaving it below current bounds but changing the 'lower limit' status.","The closed low-momentum boundary and advection-only transport mean any additional source of low-energy protons—such as turbulent re-acceleration or a supra-thermal injection—would raise emission; a direct test would be running the same volume with an open boundary or streaming transport and comparing Coma's predicted flux.","If a future telescope detects cluster diffuse emission at the level predicted here, it would indicate that shock acceleration in clusters is more efficient than the conservative model assumes, or that non-shock CR sources contribute."],"forward_implications":["If the prediction holds, cosmic-ray protons contribute negligibly to cluster pressure, so hydrostatic mass estimates need no significant correction for non-thermal support at this level.","Detecting diffuse gamma-ray emission from even the brightest cluster, Coma, will require a next-generation telescope with sensitivity below 10^-11 gamma s^-1 cm^-2; current instruments cannot reach it.","Emission from cosmic-web filaments and cluster outskirts is 4-5 orders of magnitude below cluster centers, making any near-term detection in those regions essentially impossible.","The simulated gamma-ray spectra are universal above 1 GeV with slope about -1, flatter than earlier analytic expectations, giving a concrete spectral prediction to test if emission is ever observed.","The result provides a lower-limit background for diffuse cluster emission, useful for searches for dark-matter annihilation signals in clusters."],"fun_headline_variants":["Simulated cosmic-ray glow in clusters stays hidden from Fermi","Cluster gamma-ray glow predicted 4 orders below Fermi limits","First CR spectral simulation sets gamma-ray detection benchmark","Local universe gamma-ray glow needs 10^-11 sensitivity","Cosmic-ray protons in clusters: gamma-ray flux below detection"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The prediction depends on the simulation's shock finder capturing the shocks that accelerate cosmic rays, and the authors estimate it misses a factor of 1-2 orders of magnitude of dissipated energy at low Mach numbers and suppresses proton acceleration by another factor of 5-10 due to numerically biased shock obliquities.","fun_headline_variants_meta":{"raw":{"variants":["Simulated cosmic-ray glow in clusters stays hidden from Fermi","Cluster gamma-ray glow predicted 4 orders below Fermi limits","First CR spectral simulation sets gamma-ray detection benchmark","Local universe gamma-ray glow needs 10^-11 sensitivity","Cosmic-ray protons in clusters: gamma-ray flux below detection"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000195,"raw_usage":{"total_tokens":1249,"prompt_tokens":852,"completion_tokens":397,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":596,"completion_tokens_details":{"reasoning_tokens":317}},"tokens_in":596,"tokens_out":397,"duration_ms":3752,"temperature":1.0,"reasoning_tokens":317,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T09:21:30.244412+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the same constrained local-volume simulation at 8x resolution (as the authors plan for zoom-ins) and compare the predicted Coma gamma-ray flux; if it rises by more than a factor of ~100, the stated lower limit is not robust. Alternatively, a gamma-ray telescope reaching a sensitivity of 10^-12 gamma s^-1 cm^-2 that detects diffuse emission from Coma would directly falsify the paper's prediction that this emission lies below 10^-11.","supporting_citations":[],"review_version":1}