{"id":"67379e2a-80c7-4415-96ed-f416d8b9466e","arxiv_id":"2508.17475","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A simple two-parameter bi-effect formula reproduces muon deflections in simulated air showers by combining geomagnetic Lorentz deflection with atmospheric absorption.","lead":"This paper proposes a two-parameter formula that describes how the Earth's magnetic field and atmospheric absorption jointly shift muon landing positions in cosmic-ray air showers. It fits the formula to thousands of CORSIKA simulations and argues it can guide detector placement at cosmic-ray observatories.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Independence Hypothesis is untested: additive superposition may hide B-dependent coupling; a B-on/B-off CORSIKA test would settle it.","rationale":"The reader identified the Independence Hypothesis as the weakest assumption; I agree that this is the most load-bearing unverified element. The paper's azimuthal-pattern validation is a genuine and nontrivial test of the geometric form of Eq. (3), and the reported 11%/8% errors indicate that the model can fit the simulated mean deflections. However, that test does not isolate the superposition assumption: the two free parameters can absorb a B-dependent coupling, so the good fit does not establish that the model is predictive outside the simulated conditions. The proposed CORSIKA test with B switched off and B scaled would directly check whether the magnetic contribution is additive and proportional to B, which is required for the model to be generic across sites and for the 'simple applications' in the paper. If the test passes, the conditional acceptance is justified; if it fails, the model would need to be revised or restricted to the calibration site. The paper would also benefit from clearly stating whether the 1728-showers error set was used in fitting Table II, but the independence test is the decisive missing check.","tokens_in":7855,"tokens_out":15989,"duration_ms":174692,"concrete_test":"Run CORSIKA for a fixed set of showers (e.g., proton, E = 10^16 eV, theta = 70 deg, FIXHEI = 20 km, all 72 azimuths) at the LHAASO site in three configurations: (a) true geomagnetic field, (b) B = 0, and (c) B scaled by a factor of 2 with the same direction. Fit catm from configuration (b) using the atmospheric-only model. Then compute delta(a) = Delta r(a) - catm a and delta(c) = Delta r(c) - catm a. If the Independence Hypothesis and B-linearity hold, delta(a) and delta(c)/2 should be equal, for every azimuth, to the same magnetic term cmag [b.(ep x B)] b of Eq. (3), with a single cmag matching the direct fit from (a). If the residuals differ in magnitude or direction, or cmag changes between (a) and (c), the additive superposition is invalid.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central formula (Eq. 3) relies on the Independence Hypothesis (Sec. II): Delta r = Delta r_mag + Delta r_atm. The paper validates Eq. (3) by fitting cmag and catm to simulated mean deflections and showing that the azimuthal curve is reproduced. But a two-parameter fit to a smooth curve cannot distinguish true linear superposition from an effective parametrization in which cmag and catm absorb coupling between the Lorentz-force displacement and atmospheric absorption. When a muon is magnetically deflected, it samples a different atmospheric depth, so the number- and energy-weighted means of the surviving muon population are not obtained by simply adding the B=0 displacement to an independent magnetic displacement. The validation uses only the local real magnetic field, so the additivity of the two effects and the claimed independence of cmag from B are never directly tested. This matters because the model is then used predictively (detector layout, mass composition) and claimed to be generic across sites. If superposition fails at large zenith angles or high-field sites, the fitted cmag and catm would be effective parameters that vary with B and atmospheric conditions in ways not captured by Eqs. (5) and (7). A direct test of the superposition is therefore the most load-bearing missing check.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a 'bi-effect model' for the ground-plane lateral deflection of muons in extensive air showers. The model combines two effects linearly: the geomagnetic Lorentz-force deflection and an atmospheric-absorption-induced displacement along the shower asymmetry direction, leading to Eq. (3) for the mean deflection vector as a function of the primary direction, the local magnetic field, and two scalar parameters cmag and catm. The parameters are fitted to CORSIKA simulations for proton and iron primaries, energies between 10^14 and 1.024×10^17 eV, zenith angles 50°–80°, 72 azimuth angles, and the magnetic-field geometries of LHAASO, SKA, and the Pierre Auger Observatory. The authors propose empirical analytic forms cmag=A τ^B (cosθ)^C and catm=A τ^B tanθ exp(Cτ/tanθ), where τ is the atmospheric depth, and tabulate fitted coefficients. They report that the model reproduces the simulated azimuthal double-circle structure and quote 68% relative-error intervals of 11% for the number deflection and 8% for the energy deflection, and discuss applications to detector layout and mass-composition reconstruction. The central claim is that this simple formula can describe the overall deflection of muons and accurately fit the deflection in simulated air showers, thereby validating the two underlying hypotheses.","tokens_in":8169,"tokens_out":6336,"duration_ms":63390,"significance":"If the central claim holds, the paper offers a compact analytic description of a nontrivial geometric effect, with practical value for detector design and shower reconstruction at current and future observatories. The work has notable strengths: it uses a large and systematic simulation library (about 11,000 showers, three sites, 72 azimuthal bins), tests azimuthal independence over a full circle, checks the energy dependence over three decades, and provides tabulated coefficients that allow immediate use of the model. The geometric part of the model, namely the projection of the Lorentz deflection onto the ground plane with a single azimuth-independent parameter, is physically well motivated and is supported by the observed double-circle pattern. However, the empirical validation is currently in-sample for the fitted parameters, and the linear superposition of magnetic and atmospheric effects is assumed rather than directly tested. The significance of the paper therefore hinges on the additional holdout and field-switch tests requested below.","major_comments":[{"comment":"The quoted 68% errors of 11% and 8% are computed with ε = |Δr_model − Δr_simulation| / |Δr_simulation| using simulations and fitted parameters from the same analysis chain. As written, the text does not state whether the 'several hundreds showers for testing' mentioned in Sec. II are actually held out from the fits that determine cmag and catm and from the empirical coefficient fits of Eqs. (5) and (7). If they are not, the quoted errors are in-sample residuals and cannot support the predictive claim in the Abstract. Please specify the train/test split explicitly and report the relative-error distribution separately for the fitted and held-out samples; if no true holdout was used, rerun the validation on an independent set, for example a random subset of azimuth angles or a different set of FIXHEI values.","section":"Sec. III, Model uncertainties; Fig. 6"},{"comment":"The Independence Hypothesis states Δr = Δr_mag + Δr_atm. Because a magnetically deflected muon traverses a different atmospheric depth, the number- and energy-weighted means of the detected muon population are not guaranteed to be additive in this way. The validation included in the paper uses only the real geomagnetic fields of three sites, so the possible B-dependence of cmag and catm through coupling is never probed. The central predictive claim, that the model is generic across sites and zenith angles, therefore requires a direct test: rerun a representative set of showers with the magnetic field switched off, compute Δr(B=0), and compare Δr(B) − Δr(B=0) with the magnetic term in Eq. (3); also vary the field strength (for example B=0, B, and 2B) to verify that cmag remains independent of B and that catm is unchanged. If the superposition fails at large zenith angles or high-field sites, cmag and catm are effective parameters that vary with the magnetic field, and the application in Fig. 7 overstates the predictive reach.","section":"Sec. II, Model hypotheses; Eq. (3)"},{"comment":"The functional forms for cmag and catm are chosen to satisfy the boundary conditions in Eqs. (4) and (6), but boundary condition (iv), 'FIXHEI→112 km, |Δr|>0', is not checked against simulations, and the model's own text notes that the energy deflection of highly inclined showers initiated at very high altitudes deviates from the fitted form because of muon decay, with the discrepancy dismissed because the magnetic effect dominates. This breakdown regime is exactly the regime where the independence hypothesis is least safe, so it cannot be dismissed without a quantitative estimate. Please state the ranges of zenith angle and FIXHEI over which Eqs. (5) and (7) reproduce the fitted parameters to better than, say, 10%, and quantify the muon-decay effect, for example by comparing against a simulation run with muon decay switched off, rather than ignoring it.","section":"Sec. III, Fitting formula for deflection parameters; Eqs. (5), (7), (4), (6)"}],"minor_comments":[{"comment":"The figure cross-references are inconsistent: 'FIG. III' in Sec. III should be Fig. 3, and the relative-error distribution is described in the text as Fig. 5 but the caption identifies it as Fig. 6.","section":"Throughout"},{"comment":"The units in the Table II header are unclear: the column header '([m/μT])' is attached to all three coefficients A, B, and C, but these coefficients have different dimensions in Eqs. (5) and (7). Please state the units of each coefficient explicitly.","section":"Table II"},{"comment":"The sentence 'Since the definition of τ contains cos(θ)−1, when B > C, and B > 0, the boundary conditions are satisfied' is too terse; please spell out how each boundary condition in Eq. (4) follows from the parameter conditions.","section":"Sec. III, Fitting formula for deflection parameters"},{"comment":"References [13] and [21] are identical (Bae and Chatzidakis, PTEP 2022, 043F01); please remove the duplicate.","section":"References"},{"comment":"The notation for the magnetic deflection parameter is inconsistent: the text uses both 'cmag' and 'Cmag' in the validation paragraph; please unify to a single symbol.","section":"Sec. III, Model validation"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and contains a useful geometric formula, but I cannot support acceptance in the current form. The two load-bearing issues are the lack of a demonstrated train/test separation for the quoted errors and the absence of a direct test of the linear superposition of magnetic and atmospheric effects. Both are fixable with additional simulations and analysis, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline is that this is a plausible empirical parameterization, and the azimuthal double-circle test is the real gem: a single pair of parameters (cmag, catm) reproduces the deflection pattern across 72 azimuths. That is a genuine geometric check, not a trivial fit. The novelty is moderate — combining Lorentz force and slant-depth absorption in one formula is new, but the ingredients are standard and the output is an engineering tool, not new physics.\n\nWhat the paper does well: the projection in Eqs. (2)–(3) is clean, the boundary conditions for the fitting formulas are sensible, and Table II gives coefficients for three sites, proton/iron primaries, and a decent energy span. The reported 8–11% relative errors are honest, though they are in-sample.\n\nThe soft spot is exactly what the stress test flagged: the Independence Hypothesis in Sec. II is load-bearing and never directly tested. The authors validate only at each site's real magnetic field, so a two-parameter fit can absorb B-dependent coupling between magnetic displacement and atmospheric absorption without revealing it. A CORSIKA run with B set to zero, and another with, say, twice the field, would settle whether cmag is truly B-independent and whether the superposition holds. Without that, the predictive uses for detector layout and mass composition are conditional. The lack of released code or data makes it harder to check; and there is no baseline comparison against a simpler one-parameter model. The energy range stops at 1e17 eV, which is fine for the stated scope but not for UHECR regimes.\n\nMinor quibbles: the muon-decay discrepancy at high zenith and high FIXHEI is waved off, and the empirical formulas are matched to boundary conditions, which is acceptable but limits extrapolation.\n\nVerdict: worth a serious referee. The geometric form is likely right and the model can inform array design at LHAASO, SKA, and Auger. I would send it to peer review with an explicit request for a B-switch test and a data release; I would not use it predictively before that check, but I would cite it as a useful parameterization.","headline":"A useful two-parameter muon deflection model with a convincing azimuthal shape test, but the unresolved B-coupling question makes it a conditional tool rather than a proven predictive one.","tokens_in":8627,"tokens_out":2982,"would_cite":true,"duration_ms":32277,"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 two-term vector formula, split between geomagnetic bending and atmospheric absorption, predicts where muons land in inclined air showers.","keywords":["air showers","muons","geomagnetic deflection","atmospheric absorption","lateral distribution","detector layout","zenith angle dependence","cosmic ray composition"],"falsifier":"Run air-shower simulations at a site with a substantially stronger magnetic field, or at zenith angles beyond 80°, and compare the simulated ground-plane muon deflection with Eq. (3) using parameters fitted at current sites; if the residual grows with the size of the magnetic deflection, the linear superposition fails. A second check is to add a cross term proportional to $c_{\\mathrm{mag}} c_{\\mathrm{atm}}$ to Eq. (3) and see whether it improves the fit by more than the reported 8-11% uncertainties.","tokens_in":7694,"feed_emoji":"🔭","tokens_out":6607,"duration_ms":59125,"temperature":0.7,"pith_summary":"This paper tries to establish that the sideways spreading of muons reaching the ground from inclined air showers reduces to a two-piece vector formula. One piece comes from the Lorentz force of Earth's magnetic field bending the muon's path; the other comes from atmospheric absorption, which removes muons preferentially on one side of the shower axis. The two effects are assumed to act independently and add linearly. Compared against Monte Carlo air-shower simulations at three sites with different geomagnetic fields, the formula reproduces the simulated average deflections of muon number and energy with 68% relative-error intervals of 11% and 8%. If correct, this gives a simple predictive tool for placing detector stations and for reconstructing the direction and composition of the primary cosmic ray.","feed_headline":"Muon landings match a two-term formula to 11 percent","feed_subtitle":"One term is geomagnetic bending, the other atmospheric absorption; the sum reproduces simulated air showers.","key_machinery":"The central object is the vector decomposition in Eq. (3): every muon's ground-plane deflection is split into a magnetic term along the ground-plane north-south axis and an atmospheric term along the projection of the shower direction. The magnetic linearity hypothesis approximates Larmor motion as projectile motion, making the magnetic deflection proportional to $\\vec{e}_p \\times \\vec{B}$ with a single scalar $c_{\\mathrm{mag}}$; the independence hypothesis makes atmospheric absorption contribute a displacement along $-\\vec{a}$ with scalar $c_{\\mathrm{atm}}$. The pair $(c_{\\mathrm{mag}}, c_{\\mathrm{atm}})$ is constant over azimuth, so once they are known the whole deflection pattern for a site is a geometric projection, and the empirical formulas for $c_{\\mathrm{mag}}$ and $c_{\\mathrm{atm}}$ connect those parameters to the primary's zenith angle, first-interaction height, and atmospheric depth.","core_discovery":"The central claim is Eq. (3): the ground-plane deflection vector $\\Delta \\vec{r}$ equals $c_{\\mathrm{mag}} [\\vec{b}\\cdot(\\vec{e}_p \\times \\vec{B})]\\vec{b} + [c_{\\mathrm{atm}} - c_{\\mathrm{mag}} \\vec{c}\\cdot(\\vec{e}_p \\times \\vec{B})/(\\vec{e}_p\\cdot\\vec{n})]\\vec{a}$, where $\\vec{b},\\vec{a},\\vec{n}$ are the ground-plane and shower-plane axes, $\\vec{e}_p$ is the shower direction, $\\vec{B}$ is the geomagnetic field, and $c_{\\mathrm{mag}}$, $c_{\\mathrm{atm}}$ are two scalar parameters independent of the azimuth angle of the primary. The paper reports that one fixed pair $(c_{\\mathrm{mag}}, c_{\\mathrm{atm}})$ reproduces the simulated deflection pattern for all azimuths, including its two-lobed 'double-circle' structure, and that the empirical functions $c_{\\mathrm{mag}} = A\\tau^{B}(\\cos\\theta)^{C}$ and $c_{\\mathrm{atm}} = A\\tau^{B}\\tan\\theta \\exp(C\\tau/\\tan\\theta)$ capture the dependence on zenith angle, first-interaction height, and site. The stated accuracy is a 68% relative-error interval of 11% for the number-weighted deflection and 8% for the energy-weighted deflection.","pith_inferences":["Read as an extension: if the linear superposition holds, the same two-parameter vector form should describe muon deflections at any site, so the model's core physics content is really the scaling of $c_{\\mathrm{mag}}$ and $c_{\\mathrm{atm}}$ with atmospheric depth.","A natural stress test beyond the paper: push the model into the regime where muon decay is non-negligible, which the paper notes is where the energy-weighted fit degrades; separating decay losses from magnetic-atmospheric coupling would clarify whether the linear superposition or the muon physics is responsible.","An untested consequence of the geometry: the double-circle deflection pattern is essentially the projection of a fixed shower-plane vector onto the ground plane, which implies the same pattern should appear in electromagnetic particles if their absorption can be modelled, a topic the paper leaves for separate work."],"forward_implications":["Detector upgrade footprints can be computed directly from the model: the required extra area in each direction falls out of the deflection curve, without running new simulations.","Because $c_{\\mathrm{mag}}$ and $c_{\\mathrm{atm}}$ are azimuth-independent, the deflection pattern for any arrival direction at a given site is predicted by the same two parameters.","The fitted parameters differ between proton and iron primaries, so measured deflection patterns could be used to infer the mass composition of the primary cosmic ray.","The model explicitly covers zenith angles from 50° to 80°, a regime beyond the 60° limit of earlier universality descriptions of particle distributions."],"supporting_citations":[{"why":"Supplies the Monte Carlo air-shower generator whose simulated muon deflections are compared with the model.","marker":"[18]"},{"why":"Provides the standard atmosphere model and the atmospheric-depth calculation used to build the fitting functions for the two deflection parameters.","marker":"[20]"},{"why":"Establishes the roughly constant ionization energy-loss rate (2 MeV per gram per square centimeter) that underlies the simple atmospheric absorption term.","marker":"[19]"},{"why":"Provides the recent evidence that geomagnetic effects on inclined air showers are larger than previously expected, motivating the new model.","marker":"[16]"}],"fun_headline_variants":["Two-term formula nails muon deflection to 11%","11% accuracy: bi-effect muon model","One formula, two effects: muon showers decoded","Magnetic and atmospheric bending combined: 11% match","Bi-effect model: muon deflection, one formula"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that magnetic bending and atmospheric absorption act independently and add as vectors, even though a magnetically bent muon travels through a different amount of atmosphere than an unbent one.","fun_headline_variants_meta":{"raw":{"variants":["Two-term formula nails muon deflection to 11%","11% accuracy: bi-effect muon model","One formula, two effects: muon showers decoded","Magnetic and atmospheric bending combined: 11% match","Bi-effect model: muon deflection, one formula"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000737,"raw_usage":{"total_tokens":3302,"prompt_tokens":961,"completion_tokens":2341,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":577,"completion_tokens_details":{"reasoning_tokens":2264}},"tokens_in":577,"tokens_out":2341,"duration_ms":17069,"temperature":1.0,"reasoning_tokens":2264,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:03:31.363476+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run air-shower simulations at a site with a substantially stronger magnetic field, or at zenith angles beyond 80°, and compare the simulated ground-plane muon deflection with Eq. (3) using parameters fitted at current sites; if the residual grows with the size of the magnetic deflection, the linear superposition fails. A second check is to add a cross term proportional to $c_{\\mathrm{mag}} c_{\\mathrm{atm}}$ to Eq. (3) and see whether it improves the fit by more than the reported 8-11% uncertainties.","supporting_citations":[{"cited_title":"Chiche, C","cited_arxiv_id":null,"evidence_quote":"Supplies the Monte Carlo air-shower generator whose simulated muon deflections are compared with the model."},{"cited_title":"Heck et al","cited_arxiv_id":null,"evidence_quote":"Provides the standard atmosphere model and the atmospheric-depth calculation used to build the fitting functions for the two deflection parameters."},{"cited_title":"Fuehne and T","cited_arxiv_id":null,"evidence_quote":"Establishes the roughly constant ionization energy-loss rate (2 MeV per gram per square centimeter) that underlies the simple atmospheric absorption term."},{"cited_title":"Stadelmaier et al","cited_arxiv_id":null,"evidence_quote":"Provides the recent evidence that geomagnetic effects on inclined air showers are larger than previously expected, motivating the new model."}],"review_version":2}