{"id":"d93b29d2-5b66-4ffa-a4ba-311354e864e0","arxiv_id":"2607.10403","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Upward-going muons from lunar regolith, simulated with FLUKA, yield altitude- and void-sensitive fluxes that enable shallow subsurface tomography with minute-scale exposures.","lead":"Simulations show that cosmic-ray interactions in lunar regolith produce a usable flux of upward-going muons that escape into vacuum and can be detected on the surface or in orbit. The flux depends strongly on effective altitude and shallow density anomalies, so it could map cavities and water deposits non-invasively for future lunar missions.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.5","headline":"The minute-scale cavity/water claims rest on an idealized overburden model that may not survive realistic showering and lateral geometry.","rationale":"The Reader correctly isolates the flat-terrain / equivalent-altitude approximation as the weakest link and already assigns CONDITIONAL. The present stress-test simply sharpens that concern: the quantitative minute-scale numbers that constitute the paper’s strongest claim are generated precisely under the approximation the Reader flags. Because the paper is transparent about the idealization and the underlying physics (displaced meson decay) is sound, no stronger verdict change is warranted; the concrete finite-cavity re-simulation would settle whether the exposure times survive or must be revised upward. Photon validation against Fermi-LAT still anchors the absolute scale, so the conceptual channel remains interesting even if the practical timescales lengthen.","tokens_in":11397,"tokens_out":603,"duration_ms":9679,"concrete_test":"Re-run the FLUKA setup of §II with a finite cylindrical cavity (radius 5–20 m, top at 0.2–0.5 m depth, height 10 m) embedded in the MP2007 regolith sphere; keep the same 10^6 primaries and solar-modulation settings. Compare the surface muon spectrum (E_k > 100 MeV) above the cavity center versus the pure layered-slab result of Fig. 4. If the cavity-to-nominal flux ratio drops below ~2 (or the water contrast falls below ~1.2), the <1 min / ~10 min exposure claims no longer hold for realistic geometries.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The strongest claim (abstract + §III.D) is that cavity-induced flux variations are observable in <1 min and water signals after ~10 min with a 1 m^{2} ideal detector. That claim is obtained by treating shallow anomalies as a thin (0.1–0.5 m) slab of nominal regolith over a 10 m void/water region, or equivalently as a pure change in detector altitude (Appendix B). Both approximations neglect (i) the additional hadronic showering that occurs inside any real overburden and (ii) the finite lateral extent of a cavity, which changes the solid angle of escaping mesons and the path-length distribution of the resulting muons. Because the surface muon flux is itself produced by mesons that have already escaped the last interaction length, even modest extra interactions or geometric dilution can erase the order-of-magnitude enhancement shown in Fig. 4. The paper itself flags the idealization (“highly idealized and not directly applicable in practice,” App. B), yet still quotes the minute-scale numbers as the headline sensitivity. If those numbers shrink by more than a factor of a few once realistic overburden and finite-size effects are restored, the practical claim collapses while the qualitative altitude dependence remains intact.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper proposes lunar subsurface tomography using upward-going secondary muons generated by galactic cosmic-ray interactions in the regolith. Because the Moon has no atmosphere, charged mesons that escape the surface can decay in vacuum, yielding a detectable muon flux at the surface or in near-lunar orbit. Using FLUKA (DPMJET/RQMD) with Badhwar–O’Neill spectra, the MP2007 regolith model, and solar-modulation potentials of 465 MV and 1440 MV, the authors compute energy spectra, angular distributions, and integrated fluxes for 10^6 primaries. Photon fluxes are validated against Fermi-LAT data (Appendix A). They show strong altitude dependence under a flat-terrain assumption and present idealized case studies of shallow cavities and water deposits, claiming cavity-induced flux variations are observable in less than a minute and water signals after about ten minutes with a 1 m^{2} ideal detector. Potential applications to Chang’e-7/8 and related missions are discussed.","tokens_in":11736,"tokens_out":1249,"duration_ms":25449,"significance":"If the altitude and void sensitivity survive more realistic geometry, the work would introduce a genuinely new, non-invasive probe of shallow lunar density structure and resources that complements existing X-ray and gamma-ray methods. The use of standard, publicly documented tools (FLUKA, DPMJET/RQMD, Badhwar–O’Neill, MP2007), the independent Fermi-LAT photon cross-check, and the explicit MC statistical uncertainties constitute clear methodological strengths. The qualitative recognition that displaced meson decays produce an altitude-dependent upward muon flux is sound and timely for upcoming polar and resource-utilization missions. The quantitative minute-scale claims, however, rest on highly idealized overburden and flat-terrain approximations that the authors themselves flag; the practical impact therefore remains prospective rather than demonstrated.","major_comments":[{"comment":"Abstract and §III.D (Figs. 4 and 6, Appendix B): The headline statements that cavity-induced flux variations are observable in less than a minute and water signals after ~10 min with a 1 m^{2} detector are obtained by modeling anomalies as a thin (0.1–0.5 m) slab of nominal regolith overlying a 10 m void/water region, or equivalently as a pure change in detector altitude. Both approximations neglect additional hadronic showering inside any real overburden and the finite lateral extent of a cavity (which alters solid angle and path-length distributions). The paper itself labels the scenario “highly idealized and not directly applicable in practice,” yet still quotes the minute-scale numbers as the central sensitivity result. Either more realistic 3-D overburden/finite-size simulations must be supplied, or the abstract and §III.D claims must be rewritten as qualitative upper bounds under i","section":"Abstract, §III.D, Appendix B"},{"comment":"§II and §III.C: Detector response, acceptance, and particle identification are neglected under the assumption of “high detection and particle-identification efficiencies” and a “nearly background-free environment.” The TOF timing resolutions and dual-readout arguments given for µ/π separation are order-of-magnitude estimates only; no efficiency curves, mass/power constraints, or residual background rates under lunar conditions are provided. Because the tomography feasibility argument relies on clean muon samples, at least a simplified detector-response study or a clear statement that the quoted fluxes are ideal upper limits is required.","section":"§II, §III.C"},{"comment":"Appendix B (microscopic vs. macroscopic tomography): The illustrations of sensitivity assume a flat surface and negligible cosmic-ray showering in the overburden. For macroscopic (orbital) tomography the text notes that deeper voids eventually compensate the muon-decay loss by additional production, but no quantitative depth range or contrast-versus-depth curve is given. Without this, the claimed reach of orbital muon tomography remains unquantified and cannot support mission-planning statements.","section":"Appendix B"}],"minor_comments":[{"comment":"Several reference titles contain obvious OCR/transcription errors (“cosmic yay muon,” “cosmic-yay spectrometer,” “Fermisatel-lite”). These should be corrected before publication.","section":"References"},{"comment":"The E^{2}-weighted spectra are defined with E^{2} = E_low E_high; a brief sentence clarifying that this is the conventional geometric-mean weighting used by Fermi-LAT would help non-specialist readers.","section":"§III"},{"comment":"Figure 1 caption and the schematic itself would benefit from an explicit indication of the detector locations (surface vs. altitude) that are later quantified.","section":"Fig. 1"},{"comment":"Earth shielding and surface topography are stated to be neglected “in this conceptual design.” A short quantitative estimate of the solid-angle fraction occulted by Earth (or a reference to existing lunar cosmic-ray maps) would strengthen the claim that the omission is harmless for the present purpose.","section":"§II"}],"recommendation":"major_revision","confidential_remarks":"The core physical idea (altitude-dependent upward muon flux from displaced meson decays) is novel and correctly simulated with standard tools; the photon validation is a genuine strength. The over-selling of minute-scale sensitivities under acknowledged idealizations is the main obstacle to acceptance. Once the abstract and §III.D are brought into line with the caveats already present in Appendix B, and a minimal detector-response discussion is added, the paper would be a solid contribution to lunar instrumentation concepts. Scope is appropriate for an astro-ph.EP or instrumentation-oriented journal."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The real contribution here is the first quantitative FLUKA prediction of the escaping upward-going muon flux on the airless Moon, its strong altitude dependence from displaced pion/kaon decays, and the recognition that this channel can serve as a passive density probe. Prior lunar-muon papers they cite are mostly about underground production or positioning; this one flips the geometry and shows the surface-to-orbit spectra with MC errors for 10^6 primaries.\n\nThey do the basics right. Photon spectra in Appendix A match Fermi-LAT well enough to anchor the MP2007 regolith model and the Badhwar–O’Neill normalization. Solar-modulation potentials are taken from standard values, not fitted. Energy and angular distributions (Figs. 2, 3, 6, 7) are clean, and the physical origin—pions escaping the last interaction length then decaying in vacuum—is transparent. Dual-readout and TOF PID remarks are practical for a future instrument paper.\n\nThe soft spot is exactly the one the stress-test flags, and the paper itself labels it “highly idealized.” The <1 min cavity and ~10 min water claims (abstract, §III.D, Fig. 4) rest on treating a shallow anomaly as a thin overburden slab or pure altitude shift, ignoring extra hadronic showering in real overburden and finite lateral size. Once those are restored the order-of-magnitude enhancement will shrink; how much is unknown. That does not kill the qualitative altitude/void sensitivity, but it means the headline timescales are optimistic upper bounds, not mission-ready numbers. No public decks either, so reproducibility is limited to re-running FLUKA yourself.\n\nThis is for people thinking about Chang’e-7/8 or ILRS payloads and for radiation-environment modelers. It is not a finished tomography method, but it is a clean conceptual MC study that deserves a serious referee. I would send it out; the idealizations are stated clearly enough that referees can demand the next-step geometry without desk-rejecting the idea.","headline":"Solid first MC map of the upward lunar muon channel; the minute-scale cavity/water numbers are idealized and should be read as order-of-magnitude only.","tokens_in":12320,"tokens_out":521,"would_cite":true,"duration_ms":7281,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Upward-going muons from the lunar regolith can map shallow subsurface voids and water on timescales of minutes with a one-square-meter detector.","keywords":["lunar muons","muon tomography","subsurface voids","regolith","cosmic-ray secondaries","lunar exploration","water prospecting"],"falsifier":"Place a muon telescope of known area and efficiency on a flat lunar surface or in low orbit above a region whose subsurface density is independently mapped (or above a deliberately excavated cavity) and measure whether the observed flux change matches the predicted altitude or overburden dependence within the claimed statistical precision.","tokens_in":12279,"feed_emoji":"🌑","tokens_out":607,"duration_ms":9252,"temperature":0.7,"pith_summary":"The Moon has no atmosphere, so cosmic-ray mesons that escape the dense regolith can decay in the vacuum above the surface and produce a usable upward-going muon flux. Monte Carlo simulations show that this flux changes strongly with detector altitude and with the presence of shallow density anomalies. Because the change is large, a simple one-square-meter detector on the surface or in low orbit can register cavity-induced flux variations in under a minute and weaker water-induced signals in roughly ten minutes. The same muons therefore offer a non-invasive way to prospect for underground cavities and water resources without drilling, using instruments that could fly on near-term lunar missions.","feed_headline":"Lunar muons map shallow voids in under a minute","feed_subtitle":"A 1 m^{2} detector can spot cavities fast and water after about ten minutes of counting.","key_machinery":"Altitude-dependent lunar muon flux generated by displaced decays of charged pions (and kaons) that escape the regolith into vacuum; the flux rise from surface to ~1 km and the subsequent fall at higher altitudes encode the vertical path length through vacuum or voids, turning effective height into a tomographic observable.","core_discovery":"Secondary cosmic-ray muons produced by the decay of charged mesons that escape the lunar regolith form a measurable upward-going flux whose intensity is highly sensitive to detector altitude and to shallow subsurface density anomalies. Under idealized flat-terrain conditions this altitude dependence can be exploited for both surface-based microscopic tomography of meter-scale voids and orbital macroscopic tomography of larger near-surface cavities, with cavity signals appearing in less than a minute and water signals after about ten minutes of exposure with a 1 m^{2} ideal detector.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Upward lunar muons map voids in under a minute","Regolith-escape muons probe Moon cavities fast","Altitude-sensitive muons image lunar voids","Moon regolith muons detect cavities and water","Upward-going muons tomography of lunar voids"],"cache_read_input_tokens":128,"weakest_assumption_plain":"Shallow voids or water layers can be treated as simple changes in effective detector height or as thin uniform overburden above a deeper anomaly, without realistic topography or full showering in the overburden.","fun_headline_variants_meta":{"raw":{"variants":["Upward lunar muons map voids in under a minute","Regolith-escape muons probe Moon cavities fast","Altitude-sensitive muons image lunar voids","Moon regolith muons detect cavities and water","Upward-going muons tomography of lunar voids"]},"model":"grok-4.5","effort":"low","cost_usd":0.005346,"raw_usage":{"total_tokens":1465,"prompt_tokens":765,"num_sources_used":0,"completion_tokens":76,"cost_in_usd_ticks":53460000,"prompt_tokens_details":{"text_tokens":765,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":624,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":765,"tokens_out":76,"duration_ms":6065,"temperature":1.0,"reasoning_tokens":624,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T12:00:15.655867+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Place a muon telescope of known area and efficiency on a flat lunar surface or in low orbit above a region whose subsurface density is independently mapped (or above a deliberately excavated cavity) and measure whether the observed flux change matches the predicted altitude or overburden dependence within the claimed statistical precision.","supporting_citations":[],"review_version":1}