REVIEW 3 major objections 5 minor 22 references
A bi-effect model of muon deflections in air showers
T0 review · 3 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read A two-term vector formula, split between geomagnetic bending and atmospheric absorption, predicts where muons land in inclined air showers.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [Sec. III, Model uncertainties; Fig. 6] 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.
- [Sec. II, Model hypotheses; Eq. (3)] 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.
- [Sec. III, Fitting formula for deflection parameters; Eqs. (5), (7), (4), (6)] 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.
minor comments (5)
- [Throughout] 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.
- [Table II] 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.
- [Sec. III, Fitting formula for deflection parameters] 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.
- [References] References [13] and [21] are identical (Bae and Chatzidakis, PTEP 2022, 043F01); please remove the duplicate.
- [Sec. III, Model validation] 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.
Circularity Check
Partial circularity: cmag and catm are fitted to the same CORSIKA data used to report the 11%/8% errors, and the Independence Hypothesis is embedded in the fitted form; the azimuthal double-circle test provides independent content.
-
fitted input called prediction
[Sec. III, 'Fitting formula for deflection parameters' and 'Model uncertainties' (Eq. (3))]
"By fitting the particle deflection data using these parameters, cmag and catm are expected to be expressed as functions of the parameters of the primary particles. ... Using the parameters of the Pierre Auger observatory site, by comparing theoretical and simulated deflections, the distribution of the relative error for ∆⃗ rN ∆⃗ rE are presented by Fig. 5, with 68% confidence intervals for these errors are 11% and 8%, respectively — confirming the robustness of this model."
The two parameters defining the model, cmag and catm, are obtained by fitting the CORSIKA mean deflections, and the same simulated deflections are then used to compute the reported relative error ε = |∆⃗ r_model − ∆⃗ r_simulation| / |∆⃗ r_simulation|. The 68% intervals (11% and 8%) are therefore in-sample residuals of a two-parameter fit, not out-of-sample predictions. Calling these residuals 'confirming the robustness of this model' reduces the headline accuracy claim to a goodness-of-fit statement. The azimuthal double-circle test, which uses one (cmag, catm) pair across 72 azimuth angles, is an independent check of the model's angular form, so the circularity is partial rather than total.
-
self definitional
[Sec. II 'Model hypotheses' and Sec. III 'Model validation']
"We assume that the magnetic deflection (∆⃗ rmag) caused by the Lorentz force and the atmospheric absorption-induced displacement (∆⃗ ratm) are independent. This relationship is expressed as: ∆⃗ r= ∆⃗ rmag + ∆⃗ ratm. ... using a single pair of cmag, catm, the model effectively reproduce the double-circle structure of the magnetic deflection, further validating the hypotheses and the azimuth angle independence."
The model's total deflection in Eq. (3) is constructed from exactly this Independence Hypothesis, with the atmospheric contribution represented by a fitted vector catm a and the magnetic contribution by a fitted scalar cmag. Fitting this two-parameter additive form to simulated mean deflections cannot validate additivity: the simulations always use the full real local magnetic field and never switch off B to isolate ∆⃗ ratm or vary B to isolate ∆⃗ rmag. A two-parameter additive form can absorb some coupling between the two effects while still matching the smooth azimuthal curve. Thus the claim that the data are 'further validating the hypotheses' is the independence assumption re-encountered through the fitted model, not an independent test of that assumption.
full rationale
The paper's central formula is a geometric projection of two assumed effects, and its validation rests on fits to the same CORSIKA showers used to evaluate the model's errors. This makes the reported 11% (number) and 8% (energy) relative-error intervals in-sample residuals rather than predictive tests. The azimuthal double-circle structure is nonetheless a genuine check of the functional form: one parameter pair reproduces the angle dependence across 72 azimuths, which a poorly chosen model form would not do. Similarly, the empirical formulas for cmag and catm in Eqs. (5) and (7) are openly presented as fitting functions satisfying boundary conditions; they are not disguised predictions. The self-citation [16] is motivational and not load-bearing for the derivation. Overall, the model form has independent geometric content, but the quantitative validation and the claimed confirmation of the Independence Hypothesis are partly circular because the fitted parameters and the fitted model embody the hypotheses being 'validated'. This warrants a moderate score of 5, not a higher score, because the azimuthal test and the explicit fit-based empirical parameterizations provide independent substance.
Assumptions & free parameters
free parameters (4)
- cmag =
Site- and species-dependent; e.g., AUGER proton Delta r_N Amag = 3.463 ± 0.010 m/uT
- catm =
Site- and species-dependent; e.g., AUGER proton Delta r_N Aatm = -130.4 ± 1.0 m
- Empirical coefficients A, B, C for cmag =
Table II, e.g., Amag = 3.463, Bmag = 3.024, Cmag = 1.381 for AUGER proton Delta r_N
- Empirical coefficients A, B, C for catm =
Table II, e.g., Aatm = -130.4, Batm = 1.768, Catm = -1.013 for AUGER proton Delta r_N
assumptions (4)
- domain assumption Independence Hypothesis: Delta r = Delta rmag + Delta ratm; magnetic deflection and atmospheric absorption displacement are independent and add linearly.
- domain assumption Magnetic Linearity Hypothesis: Delta l_mag = cmag (ep cross B), treating Larmor motion as projectile motion with constant Lorentz-force direction.
- domain assumption Linsley's standard atmosphere model provides the density profile for computing atmospheric depth tau via GDASTOOL.
- ad hoc to paper The functional forms cmag = A tau^B (cos theta)^C and catm = A tau^B tan theta exp(C tau / tan theta) are chosen to satisfy boundary conditions and to fit simulations.
Cite this review
Pith. "Pith review of A bi-effect model of muon deflections in air showers." pith.science (2026). https://pith.science/paper/WZZHZ2UM
@misc{pith2026250817475,
author = {Pith},
title = {Pith review of: A bi-effect model of muon deflections in air showers},
year = {2026},
howpublished = {\url{https://pith.science/paper/WZZHZ2UM}},
note = {Machine review of arXiv:2508.17475}
}
abstract
Recent progress has shown that the geomagnetic field exerts a more significant impact than expected on the behavior of charged secondary particles in inclined air showers. In this study, we for the first time combine it with atmospheric effects to construct a bi-effect model, aiming to investigate the lateral distribution of particles on the ground plane. Despite the complex physical interactions during the development of air showers, a simple formula can describe the overall deflection of $\mu^{\pm}$ and accurately fit the deflection in simulated air showers, thereby validating the hypotheses about these effects in this study. Furthermore, we have obtained the relationship between model parameters and primary particle information for different experimental sites. This new model is highly successful and is promising to provide new insights for improving detector layout design and air shower reconstruction.
Figures
Reference graph
Works this paper leans on
-
[1]
This relationship is expressed as: ∆⃗ r= ∆⃗ rmag + ∆⃗ ratm
Independence Hypothesis: We assume that the magnetic deflection (∆⃗ rmag) caused by the Lorentz force and the atmospheric absorption-induced displacement (∆⃗ ratm) are independent. This relationship is expressed as: ∆⃗ r= ∆⃗ rmag + ∆⃗ ratm
-
[2]
Magnetic Linearity Hypothesis: Charged particles undergo Larmor motion in the geomagnetic field. When the change in the direction of particle motion is mini- mal, we approximate the direction of the Lorentz force as constant; under this condition, Larmor motion is approx- imated as projectile motion. The deflection of muons in the shower plane ∆⃗l is prop...
- [3]
-
[4]
R. Alves Batista et al. , Front. Astron. Space Sci. 6, 23 (2019). 7
work page 2019
- [5]
- [6]
-
[7]
T. Abu-Zayyad et al. (Telescope Array), Astrophys. J. Lett. 768, L1 (2013)
work page 2013
- [8]
Show all 22 references
-
[9]
Aab et al
A. Aab et al. (Pierre Auger), Phys. Rev. Lett. 126, 152002 (2021)
2021
-
[10]
Pierog, PoS ICRC2017, 1100 (2018)
T. Pierog, PoS ICRC2017, 1100 (2018)
2018
- [11]
-
[12]
S. P. Knurenko and I. S. Petrov, Phys. Rev. D 102, 023036 (2020)
2020
-
[13]
Mayotte, in ICRC2025, (2025) arXiv:2507.10292 [astro-ph.HE]
E. Mayotte, in ICRC2025, (2025) arXiv:2507.10292 [astro-ph.HE]
2025 arXiv
-
[14]
Aab et al
A. Aab et al. (Pierre Auger), JINST 16 (07), P07016
-
[16]
Stadelmaier et al
M. Stadelmaier et al. , Phys. Rev. D 110, 023030 (2024)
2024
-
[17]
Stadelmaier, in ICRC2025, (2025) arXiv:2507.13209 [astro-ph.HE]
M. Stadelmaier, in ICRC2025, (2025) arXiv:2507.13209 [astro-ph.HE]
2025 arXiv
-
[18]
Chiche, C
S. Chiche, C. Zhang, F. Schl¨ uter,et al., Phys. Rev. Lett. 132, 231001 (2024)
2024
-
[19]
Fuehne and T
D. Fuehne and T. Heibges (JEM-EUSO), PoS ICRC2023, 363 (2023)
2023
-
[20]
Heck et al
D. Heck et al. , CORSIKA: A Monte Carlo code to simu- late extensive air showers, FZKA-6019 (1998)
1998
-
[21]
Barazandeh et al
C. Barazandeh et al. , J. Phys. Conf. Ser. 770, 012050 (2016)
2016
-
[22]
Mitra et al., Astroparticle Physics 123, 102470 (2020)
P. Mitra et al., Astroparticle Physics 123, 102470 (2020)
2020
-
[23]
Bae and S
J. Bae and S. Chatzidakis, PTEP 2022, 043F01 (2022)
2022
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.