{"id":"19ed07ab-d667-43bc-95e9-4b9c26447323","arxiv_id":"2608.02718","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A beam-dynamics-aware simulation of a 10 TeV muon collider finds an intense forward neutrino beam with about 10^9 neutrino interactions per year in a 3.2 tonne detector and about two rock-produced muons per bunch crossing.","lead":"This paper simulates the neutrino beam made by decaying muons in a 10 TeV muon collider ring, finding about a billion neutrino interactions per year in a small downstream detector. It also predicts roughly two high-energy muons from the surrounding rock will cross the detector each bunch crossing, an important background for future measurements.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline rates and prism-washout claim rest on a single non-public, unvaried IMCC lattice input; a perturbation scan of the IR optics is needed to establish that the O(10^9) event rate and 1-m spot size are not lattice-specific.","rationale":"The paper's central claim is a detailed simulation result: realistic muon beam dynamics make the forward neutrino spot size set by beam divergence, giving O(10^9) interactions/yr in a 3.2 t detector at 5 km and about two rock-produced secondary muons per bunch crossing. The described Monte Carlo, the use of standard PDFs and cross sections, and the muon-propagation validation in Appendix C are appropriate support. The qualitative prism-washout conclusion follows robustly from the baseline optics, since the relevant divergence is 0.1-1 mrad, far above 1/gamma ~ 0.02 mrad. I found no internal inconsistency that would invalidate the physics conclusions. One text issue: Sec. II.A quotes gamma c tau_mu = 31.2 km, but for a 5 TeV muon the value is 31,200 km; the subsequent 1.92 decay lengths and 85% decay are consistent with the correct value, so this is a units typo rather than a physics error. The weakest assumption is the external, unversioned hybrid lattice, exactly as the reader identified. Because the paper does not quantify how the accepted flux varies with lattice parameters, the O(10^9) rate and the degree of prism washout carry an unquantified optics systematic. This is a condition on the headline numbers, not a demonstrated error; it warrants the CONDITIONAL verdict already issued, and my stress-test does not move that verdict.","tokens_in":31337,"tokens_out":23938,"duration_ms":233874,"concrete_test":"Rerun MINT with a public, versioned IMCC lattice and also with a perturbation scan: vary the straight-section beta* by factors of 0.5, 1, and 2 and the normalized emittance by +/-20%, keeping the benchmark detector unchanged. Compare the radial flux profile at L=5 km, the integrated event rate in the 1.3-m-radius fiducial face, and the energy-radius map of Fig. 7. If the accepted flux and event rate stay within about 30% of the hybrid-lattice result and the energy-radius map remains broad with no re-emergent prism structure, the concern is resolved. If the event rate shifts by more than about 50% or the prism effect reappears for any perturbation within the IMCC design range, the headline rates should be reported as a band over lattice variants rather than as single numbers.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing condition for the central claim is that the hybrid v0.6+v0.9 IR lattice (Ref. [15], private communication) is representative of the eventual 10 TeV MuC optics. This lattice is the sole input for the muon-beam divergence and the chicane/arc geometry that set the ~1-m neutrino spot size, the washing out of the energy-angle prism effect shown in Fig. 7, the 60% detector acceptance, and hence the O(10^9) yr^-1 event rate and the ~2 secondary muons per bunch crossing. The paper does not vary the lattice or release its parameters or a commit hash, so the sensitivity of these claims to the beta functions, emittance, and chicane geometry is unquantified. The claim is not internally inconsistent: the public MINT code, the standard cross-section treatment, and the nuPyProp validation in Appendix C provide independent support. Still, the forward flux profile is optics-dominated: a factor-of-two change in the straight-section divergence changes the spot area and the accepted fraction, and a substantially larger beta* could partially restore the energy-angle correlation that Fig. 7 uses to demonstrate washout. This is a genuine correctness risk, and it is the least secure load-bearing link in the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper computes the forward neutrino flux produced by muon decays in the straight section and nearby arcs of a 10 TeV muon collider, using a new Monte Carlo code MINT that combines a MAD-X/IMCC lattice description of the muon orbit and optics with standard muon-decay kinematics. The central results are: at a 5 km benchmark detector, the neutrino spot is set by the muon-beam divergence (0.1-1 mrad) rather than by 1/gamma, so the energy-angle \"prism\" effect is washed out; the benchmark detector with a 3.2 t fiducial volume accumulates about 10^9 neutrino-nucleus interactions per year; roughly two high-energy secondary muons per bunch crossing reach each detector from neutrino interactions in the upstream rock; secondary tau-neutrino and wrong-sign neutrino fluxes are small; and neutrino upscattering in rock or detector can give interesting event-yield targets for heavy neutral leptons with mixing or dipole couplings. The paper includes appendices on the neutrino angular distribution, the Courant-Snyder beam dynamics used in MINT, and validation of the muon energy-loss treatment with nuPyProp.","tokens_in":31577,"tokens_out":10376,"duration_ms":99541,"significance":"If the quoted rates hold, this is a useful quantitative basis for the forward-physics program of a 10 TeV muon collider. The paper has clear strengths: MINT is made public; the Standard Model calculation is parameter-free in the sense that no fitted quantity enters the flux or rate; the cross sections use CT18NNLO PDFs and standard CC/NC formulas; and the muon propagation is validated against nuPyProp in Appendix C. The identification of beam-divergence-dominated spot size and the washout of the prism effect is a concrete, falsifiable prediction that goes beyond earlier simplified flux estimates. The secondary-muon, wrong-sign neutrino, and tau-neutrino calculations are useful reference predictions, and the HNL event-yield maps are explicitly framed as targets for future detector and background studies. The main caveat is that the numerical results are anchored to a single non-public, unvaried IMCC lattice, and a few secondary claims rest on estimates that are less fully developed than the primary flux calculation.","major_comments":[{"comment":"The headline numbers - the ~1 m neutrino spot, the 60% detector acceptance, the O(10^9) interactions/year, the washout of the prism effect in Fig. 7, and the roughly two secondary muons per bunch crossing - all depend on a single hybrid v0.6+v0.9 IMCC lattice supplied as a private communication (Ref. [15]). The lattice file is not released with MINT and no variation over beta*, emittance, chicane bending angles, or straight-section length is presented. The constant-divergence curves in Fig. 5 are useful for understanding the role of divergence, but they are not an optics scan: they replace the full lattice by a single Gaussian width rather than perturbing the actual lattice parameters. Because the beam divergence and chicane/arc geometry set the angular profile, a factor-of-two change in the straight-section divergence would change both the spot area and the accepted flux, and a substantially larger beta* could partially restore the energy-angle correlation that Fig. 7 uses to demonstrate washout. I ask the authors to either release the lattice input (or a parameterized surrogate) together with MINT, or add a systematic scan over the relevant optics parameters and report how the event rate, acceptance, spot size, and secondary-muon rate vary. Without this, the numerical claims are tied to a single unverifiable external input.","section":"Sec. II.A; Figs. 3 and 7; Ref. [15]"},{"comment":"The paper advertises the secondary muons as \"highly polarized\" in the abstract and concludes |P_mu| > 0.999, but the supporting estimate is a single sentence: folding the Highland angle over the rock column gives a flux-averaged depolarization of ~2 x 10^-4. No formula for the spin precession, no treatment of energy-dependent depolarization during stochastic energy losses, and no uncertainty estimate are given, and the text itself concedes that a dedicated simulation would be needed. Since the polarization enters the abstract, the conclusions, and the positron spectrum in Fig. 12, this claim should either be backed by a quantitative spin-depolarization calculation or softened to \"expected to remain highly polarized, pending a dedicated simulation study.\"","section":"Sec. III.A; Fig. 12; Appendix C"}],"minor_comments":[{"comment":"Several labels render the mu-minus symbol as \"mu box\" (e.g., Figs. 6, 10, 13 and the \"mu box beam\" labels); please fix the glyph encoding throughout.","section":"Figures 6, 10, 13 and text"},{"comment":"The sentence \"The lattice accounts for about 1.5 km of the entire ring, which we assume to be total circumference of 10 km\" is ambiguous: clarify that MINT propagates the 1.5 km lattice segment and that the multi-turn storage and 85% decay fraction are folded in through the stated normalization, or describe the alternative procedure explicitly.","section":"Sec. II.A, p. 3"},{"comment":"The integration measure in Eq. (25) uses dPi and dL without definitions; please define the phase-space and propagation-length variables, or simplify the notation.","section":"Eq. (25)"},{"comment":"The two panels use different x-axis ranges and color scales; the caption notes this, but a common scale or an inset would help readers compare the prism effect before and after including beam dynamics.","section":"Fig. 7 caption"}],"recommendation":"major_revision","confidential_remarks":"The central physics is sound and the paper is a valuable contribution, but the quantitative claims rest on a single private lattice that is not released or varied. If the authors can make the lattice input available or provide a credible optics-sensitivity scan, the paper would be suitable for publication. The polarization claim in Sec. III.A is the second point I would insist on addressing before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper delivers what the title promises: a realistic forward neutrino flux for a 10 TeV muon collider, based on a public Monte Carlo (MINT) that includes the actual lattice optics, and the central numbers look credible. The genuinely new results are that beam divergence (0.1–1 mrad) rather than the 1/gamma decay cone sets the ~1 m spot at 5 km, which washes out the energy-angle prism effect; and that the rock upstream produces roughly two high-energy secondary muons per bunch crossing in the benchmark detector. The detector design is concrete and the fiducial-volume event rates (O(10^9)/yr) follow from standard cross sections and stated assumptions.\n\nThe paper earns credit for several things: MINT is public, the muon propagation in rock is validated against nuPyProp and a CSDA check, the PDF and cross-section choices are standard, the charm and tau secondary estimates are carefully built up from known processes, and the BSM sensitivity plots are honestly labeled as event-yield targets before backgrounds. The note about the independent simultaneous work by Burk et al. is appropriate and not a red flag.\n\nThe main soft spot is exactly where the stress-test note points: the entire forward flux profile depends on one IMCC hybrid v0.6+v0.9 lattice supplied by private communication, and the paper does not release its parameters or vary it. That said, the paper partially mitigates this by comparing against constant-divergence Gaussians (Fig. 5), and the rates change by only tens of percent for the plausible divergence variations shown. What is less tested is the prism-washout conclusion: a factor ~100 larger beta* would restore the energy-angle correlation, as the authors themselves note, and a quantitative scan of beta* and chicane geometry would make the claim much firmer. This is a moderate concern, not a fatal one; the physics of a high-intensity neutrino beam from decay-in-straight-section is robust.\n\nMinor issues: no systematic uncertainty budget for PDF, lattice, or rock-model variations; the astrophysical normalization of charm fragmentation parameters (epsilon_P, beta) is taken from older references and could be revisited; and the polarization estimate for secondary muons rests on a simplified Highland-style treatment, though the authors flag this as an estimate.\n\nWho is this for? Anyone planning forward neutrino physics at a muon collider, and designers of the interaction region who care about neutrino backgrounds and possible neutrino-dedicated straight sections. It deserves a serious referee and will be a standard reference, but the lattice reproducibility issue should be addressed in a revision—ideally with the TFS file or at least a versioned public lattice and a sensitivity scan. I would send it to peer review without hesitation.","headline":"First beam-dynamics-aware forward neutrino flux for a 10 TeV muon collider, with a robust headline rate and a real but curable soft spot in the single private lattice input.","tokens_in":32187,"tokens_out":1882,"would_cite":true,"duration_ms":22356,"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":"Realistic beam optics, not decay kinematics, set the size and shape of the forward neutrino beam at a 10 TeV muon collider, and this determines what a forward detector can measure.","keywords":["muon collider","forward neutrino flux","beam dynamics","neutrino interactions","secondary muons","heavy neutral leptons","tau neutrinos","neutrino scattering"],"falsifier":"Measure the transverse profile and energy-radius correlation of the forward neutrino beam at a 5 km detector in a real 10 TeV muon collider. If the neutrino spot turns out to be set by the intrinsic 1/gamma opening angle (about a few centimeters at 5 km) rather than by the 0.1-1 mrad muon-beam divergence (about a meter), and if the energy-angle prism correlation survives, then the claim that beam dynamics dominate the forward flux is wrong. A cheaper check: replace the hybrid lattice with a hypothetical straight-section lattice with 100 times larger beta function at the interaction point and recompute the radial flux; the paper's own argument predicts the prism effect would reappear.","tokens_in":31075,"feed_emoji":"⚛️","tokens_out":9414,"duration_ms":78916,"temperature":0.7,"pith_summary":"Muon decays in a 10 TeV muon collider ring create an intense, flavor-pure neutrino beam with energies up to several TeV. This paper argues that a realistic treatment of the muon beam dynamics, rather than simple decay kinematics, determines what a forward detector actually sees, and it quantifies the consequences. With a benchmark detector 5 km from the interaction point, the flux gives about $10^{9}$ neutrino-nucleus interactions per year in a 3.2 tonne fiducial volume, with a beam spot size of about a meter set by the muon beam divergence. The same calculation predicts a steady stream of highly polarized secondary muons from neutrino interactions in the rock, and it shows that the large neutrino exposure can probe heavy neutral leptons in parameter space that other experiments have not covered. The bottom line is that the forward neutrino program is a quantitative, real physics opportunity of the muon collider, but its geometry is controlled by accelerator optics.","feed_headline":"10 TeV muon collider yields 10^9 neutrino interactions yearly","feed_subtitle":"A forward detector at 5 km sees a meter-wide beam plus two polarized TeV muons per bunch crossing.","key_machinery":"The central engine is MINT, a Monte Carlo that places muon decays along the central orbit of the supplied collider lattice, samples transverse offsets and angles from Gaussian envelopes set by the local optics parameters (beta functions, emittance, dispersion), and propagates the decay neutrinos to detector planes. Beam divergence enters through the local angular width of the muon beam, so the neutrino profile is the convolution of the intrinsic 1/gamma decay cone with the muon angular spread; an accompanying semi-analytic convolution in the appendix gives the same picture. The same machinery, extended with deep-inelastic-scattering kinematics, charm production, muon propagation with energy loss and multiple Coulomb scattering, and heavy-neutral-lepton upscattering cross sections, produces all the secondary and new-physics rates.","core_discovery":"The paper establishes that the forward neutrino beam from muon decays in the 10 TeV muon collider ring is shaped more by the accelerator lattice than by the decay kinematics. Because the muon beam is strongly focused at the interaction point, its angular divergence reaches 0.1-1 mrad, far above the intrinsic 1/gamma ~ 0.02 mrad opening angle of the decay, so the neutrino spot at a 5 km detector is about a meter wide and the energy-angle 'prism' correlation is washed out. With a 3.2 tonne fiducial volume (a compact vertex tracker plus gaseous argon TPCs), the benchmark detector records roughly $10^{9}$ neutrino-nucleus interactions per year from about 7 x $10^{18}$ neutrinos crossing its face. Neutrino interactions in the upstream rock generate about two highly polarized (>0.999) TeV secondary muons crossing the detector per bunch crossing, a secondary tau-neutrino flux of roughly 0.2 charged-current events per year across both detectors, and wrong-sign neutrino rates at the $10^{-9}$ level relative to the primary beam. The same neutrino-on-target exposure gives event-yield targets for heavy neutral leptons produced by mixing or dipole-portal upscattering that can reach unexplored parameter space at masses up to tens of GeV.","pith_inferences":["Beyond the paper: if the prism effect is as thoroughly washed out as claimed, forward detectors will need tracking calorimetry or other energy estimators rather than radial position to reconstruct neutrino energy, a consequence the paper leaves implicit.","Beyond the paper: the same beam-optics sensitivity makes the neutrino spot a design lever; a dedicated straight section with much larger beta function could shrink the spot and restore the energy-angle correlation, at the cost of larger apertures, and the paper explicitly leaves that feasibility study to future work.","Beyond the paper: the highly polarized, roughly 1.2 TeV secondary muon beam, about 10^12 muons per year, is itself a physics resource, effectively a fixed-target polarized muon source, though the paper only notes the comparison and does not develop an experimental program around it.","Beyond the paper: since wrong-sign neutrino rates are at most 10^-9 of the primary rate, any observed excess of wrong-sign flavor would be a sharp new-physics indicator, a conclusion the paper does not draw explicitly."],"forward_implications":["A 5 km forward detector with a 3.2 tonne fiducial volume will record about 10^9 neutrino-nucleus interactions per year, with roughly 5 x 10^5 neutrino-electron scattering events, opening precision electroweak and QCD measurements.","The neutrino spot size is set by the 0.1-1 mrad muon beam divergence rather than by 1/gamma, so detector apertures must match a roughly meter-scale beam and the energy-angle prism correlation cannot be used to tag neutrino energies.","About two high-energy, highly polarized (P > 0.999) secondary muons cross each detector per bunch crossing, so vetoes must reject them event by event rather than bunch by bunch.","Secondary tau neutrinos from neutrino interactions in the rock give only about 0.2 charged-current events per year across both detectors at 5 km, making tau-appearance studies from rock interactions likely out of reach.","A ten-year exposure yields event-yield targets for heavy neutral leptons, reaching mixings as low as |U_muN|^2 ~ 2 x 10^-11 near m_N ~ 15 GeV and dipole couplings down to roughly 4 x 10^-10 GeV^-1 near m_N ~ 6 GeV, before detector backgrounds are applied."],"supporting_citations":[{"why":"Supplies the baseline beam parameters (bunch population, emittance, repetition rate, bunch length) used to normalize all rates.","marker":"[7]"},{"why":"Supplies the collider lattice design whose central orbit and focusing optics define the neutrino geometry; this is the load-bearing external input.","marker":"[15]"},{"why":"Public release of the MINT Monte Carlo that produces the flux, rate, and secondary-particle results.","marker":"[48]"},{"why":"Provides the lattice table format and optics engine that defines the beam line MINT reads in.","marker":"[49]"},{"why":"Supplies the parton distribution functions used for neutrino-nucleon and charm-production kinematics.","marker":"[54]"},{"why":"Provides the neutrino-nucleon cross-section tables used to set the absolute event rates.","marker":"[55]"},{"why":"Supplies muon energy-loss and survival tables in standard rock used to propagate secondary muons to the detector.","marker":"[57]"},{"why":"Supplies the multiple-Coulomb-scattering formula used to add angular deflections to propagating muons.","marker":"[58]"},{"why":"Provides the heavy-neutral-lepton upscattering cross sections used for the new-physics event-yield estimates.","marker":"[88]"},{"why":"Provides the resonant vector-meson production calculation used in the secondary tau-neutrino estimate.","marker":"[52]"}],"fun_headline_variants":["Billion neutrino hits per year from 10 TeV muon collider decays","Forward neutrino flux at muon collider: 10^9 interactions/year","MINT simulation reveals meter-wide neutrino beam at 10 TeV collider","Two polarized TeV muons per crossing from rock-induced showers"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the specific arrangement of magnets and straight sections used to model the beam is representative of the eventual 10 TeV muon collider; that arrangement was supplied privately and is not varied in the study, so a different final design could change the neutrino spot size, the loss of the energy-angle correlation, and every downstream rate.","fun_headline_variants_meta":{"raw":{"variants":["Billion neutrino hits per year from 10 TeV muon collider decays","Forward neutrino flux at muon collider: 10^9 interactions/year","MINT simulation reveals meter-wide neutrino beam at 10 TeV collider","Two polarized TeV muons per crossing from rock-induced showers"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000842,"raw_usage":{"total_tokens":3759,"prompt_tokens":1128,"completion_tokens":2631,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":744,"completion_tokens_details":{"reasoning_tokens":2550}},"tokens_in":744,"tokens_out":2631,"duration_ms":18486,"temperature":1.0,"reasoning_tokens":2550,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:00:28.996613+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the transverse profile and energy-radius correlation of the forward neutrino beam at a 5 km detector in a real 10 TeV muon collider. If the neutrino spot turns out to be set by the intrinsic 1/gamma opening angle (about a few centimeters at 5 km) rather than by the 0.1-1 mrad muon-beam divergence (about a meter), and if the energy-angle prism correlation survives, then the claim that beam dynamics dominate the forward flux is wrong. A cheaper check: replace the hybrid lattice with a hypothetical straight-section lattice with 100 times larger beta function at the interaction point and recompute the radial flux; the paper's own argument predicts the prism effect would reappear.","supporting_citations":[],"review_version":2}