{"id":"3e91fee0-7165-4dde-bd20-dd950c9cafaf","arxiv_id":"1908.04150","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"A Monte Carlo and KASCADE data study uses local age and segmented slope parameters from electron and muon lateral densities to infer a light-to-heavy cosmic ray composition change across the knee.","lead":"This paper simulates cosmic ray air showers and then uses two new shape parameters, the local age and the segmented slope, to compare simulated proton and iron showers with data from the KASCADE experiment. The authors conclude that the cosmic ray mass composition changes gradually from light to heavy across the knee, but the supporting equations and the strength of the evidence contain serious inconsistencies.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The local-age definition in Eq. (1) with α1=4.5, α2=2 does not return the NKG age s; it returns s + 2.5/(1+2x), so every reported LAP value and its radial trend is biased.","rationale":"The reader's weakest assumption is exactly the same algebraic issue, and direct algebra from Eq. (3) confirms it: the printed exponents are swapped. This is the most load-bearing concern because every numerical result in Sec. V is expressed in terms of slocal and βss, and the radial variation of slocal (the scaling claim) is partly a product of the 2.5/(1+2x) term. A constant offset would not change the p-versus-Fe separation, but the offset varies with r and with Rm, so the shapes and minima are distorted. Thus the CONDITIONAL verdict is appropriate: the authors must correct Eq. (1) (or the exponent assignment), confirm the implementation, and demonstrate that Figs. 4-6 and the composition conclusion survive. I do not see an additional concern that would change the verdict to reject; the KASCADE data are publicly available and the simulation description is sufficient for the requested re-check.","tokens_in":14171,"tokens_out":7775,"duration_ms":78657,"concrete_test":"Take a simulated NKG density profile with known age s (or any CORSIKA electron event with a fitted s), compute alocal from Eq. (1) exactly as printed at r=44 m; if it returns s+1.26 rather than s, the formula error is confirmed. Then redo the analysis of Figs. 4-6 with α1=2, α2=4.5, and check whether the KASCADE points still lie between the p and Fe bands and whether the gradual light-to-heavy trend persists; if the ordering changes, the composition claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (3) defines the NKG profile as f(r)=C x^(s-2)(1+x)^(s-4.5) with x=r/Rm. Substituting this into Eq. (1) with the paper's stated α1=4.5, α2=2 gives Fij Xij^4.5 Yij^2 = Xij^(s+2.5) Yij^(s-2.5); in the continuous limit this yields alocal(r)=s+2.5/(1+2x), not s. The choice that cancels s in both factors is α1=2, α2=4.5; the manuscript has the two exponents swapped. At r=44 m for electron data (Rm=89 m), the offset is about 1.26; for muon data with Rm=320/420 m it is about 2.0. Because the offset is radial-distance dependent, the entire shape of the slocal(r) curves, the location/depth of the minimum, the Rm dependence in Fig. 3, and the quantitative 'age' scales in Figs. 4-6 are affected. If the implementation follows the printed formula, the paper's observable is not the lateral shower age, and the composition comparison is performed with a biased variable; if the implementation used the swapped exponents, the text misstates the method. Either way the central claim needs correction and re-checking. The paper's own admission in Sec. V.A that the muon-LAP trend is 'not obvious' makes the broad conclusion in Sec. VI.4 overreach until the corrected observable is computed.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces two local shape observables for extensive air showers: the local age parameter (LAP), computed from the NKG lateral density function, and the segmented slope parameter (SSP), computed from a Greisen-like muon lateral density function. Using CORSIKA proton and iron simulations at KASCADE energies, the authors study radial profiles of these parameters and define the mean minimum LAP and the mean maximum SSP as functions of electron shower size or truncated muon size. They then compute the same observables from public KASCADE (KCDC) data and interpret the comparison as evidence for a gradual light-to-heavy change of cosmic-ray mass composition across the knee. The paper also reports a first study of radial LAP/SSP profiles for muon lateral densities.","tokens_in":14506,"tokens_out":8252,"duration_ms":76589,"significance":"If the proposed shape parameters were cleanly defined, the approach could offer a useful composition-sensitive observable that uses only lateral density data, and the comparison with public KASCADE data would be a valuable external benchmark. The paper has concrete strengths: it uses a dedicated CORSIKA simulation sample, checks two high-energy hadronic interaction models, and is the first to report radial LAP/SSP profiles for muon densities. However, the central algebraic definition of the LAP is inconsistent with the NKG form as stated, so the method as written does not return the lateral shower age. Because the LAP is the main observable behind Figures 1, 3, 4, and 5 and the primary composition conclusion, this is a load-bearing technical error that must be corrected and the analysis redone before the claims can be evaluated.","major_comments":[{"comment":"The stated exponents α1=4.5 and α2=2 are inconsistent with the NKG lateral distribution in Eq. (3). Substituting f(r)=C x^(s-2)(1+x)^(s-4.5) into Eq. (1) gives ln(Fij Xij^{α1} Yij^{α2}) = (s-2+α1) ln Xij + (s-4.5+α2) ln Yij. For this to equal s (ln Xij + ln Yij), one needs α1=2 and α2=4.5, not the values printed in the text. With the printed values, the continuous limit in Eq. (2) yields slocal(r)=s+2.5/(1+2r/Rm), a radial-dependent offset. Quantitatively, at r=44 m with Rm=89 m the offset is about 1.26, and with Rm=320 or 420 m it is about 2.0. Thus every reported LAP value, the radial position and depth of the minimum, the Rm dependence shown in Figure 3, and the quantitative age scales in Figures 4 and 5 are biased. If the implementation actually used the swapped exponents, the text misstates the method; if it followed the text, the reported observables are not lateral shower ages. The paper must correct the definition, recompute all LAP-based quantities, and re-examine the composition conclusions.","section":"II, Eqs. (1)-(2)"},{"comment":"Section V.A first states that the KASCADE muon data 'also indicates a heavier domination with increasing muon size/energy across the knee,' but the immediately following sentence says that the idea of a gradual transition from lighter to heavier composition is 'somewhat at least not obvious in experimental results.' These statements are mutually contradictory. Given that the muon-LAP trend is acknowledged not to be obvious, the strong conclusion in Section VI.4 that Figures 4, 5, and 6 indicate a gradual light-to-heavy change is an overstatement, independent of the algebraic issue in Eq. (2). The authors should either provide a quantitative measure of the trend or temper the conclusion.","section":"V.A and VI.4"},{"comment":"The composition inference is only a qualitative bracketing of the KASCADE data between proton and iron simulation curves. No mixture fraction, goodness-of-fit statistic, or compatibility test between the data and the p/Fe templates is presented. Since the central claim is a 'gradual change' in composition, the paper should quantify how well the data points are described by mixtures of the two templates, with uncertainties, rather than relying on visual interpolation.","section":"V, Figs. 4-6"}],"minor_comments":[{"comment":"The text contains the typo 'fist time' in the sentence introducing the first muon LAP/SSP study; it should be 'first time.'","section":"IV"},{"comment":"The caption refers to 'Fig. c and Fig. d'; these should be labeled '(c)' and '(d)' consistently with the other panels.","section":"Fig. 1 caption"},{"comment":"The caption contains the typo 'QGSJeT'; it should be 'QGSJet' for consistency with the text.","section":"Fig. 5 caption"},{"comment":"The phrase 'each point' and the notation alocal(i,j) for distinct radii is a little confusing; the paper would benefit from an explicit statement that Eq. (2) is obtained in the limit ri -> rj and how the discrete experimental estimates are assigned to radial positions.","section":"II"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper does something genuinely new—applying the local age parameter to muon lateral distributions and introducing a segmented slope parameter, then testing both on public KASCADE data. That is worth a look. But there is a load-bearing algebraic problem in the definition of the local age that the authors need to fix before the results can be trusted.\n\nEq. (1) defines the local age from Fij and Xij^α1 Yij^α2. For a pure NKG profile, f(r)=C x^(s-2)(1+x)^(s-4.5), substituting the stated α1=4.5, α2=2 gives alocal = s + 2.5/(1+2x), not s. The correct exponents for NKG are α1=2, α2=4.5; the paper has them swapped. The error shifts the local age by a distance-dependent offset—about 1.3 at 44 m for electrons (Rm=89 m), about 2 for muons with Rm=320/420 m. That changes the minimum location, the shape of the slocal(r) curves, and the composition comparison. If the code used the printed formula, the KASCADE data points in Figs. 4–6 are plotted with a biased variable. If the code used the corrected exponents, the text misstates the method. Either way, the central methodological claim must be re-checked.\n\nOn the positive side: the simulation sample is described in enough detail to reproduce (CORSIKA, EPOS 1.99/QGSJet/EPOS-LHC, energy cuts), the KCDC data handling is explicit, and the paper is honest that the composition conclusion is not new—it agrees with KASCADE's Ne–Nµ result. The extension to muon LAP and SSP is a legitimate idea and the first time this is done with KASCADE muon data. Also, Sec. V.A concedes the experimental muon-LAP trend is 'not obvious', which is candid.\n\nSoft spots beyond the algebra: the composition inference is qualitative—no mixture fit, just template comparison of proton vs iron curves. The statistical errors are reported as ±0.03–0.05 for electrons and ±0.10–0.15 for muons, which is fine, but no systematic uncertainties are discussed. The self-citation pattern is heavy, but the cited earlier work is directly relevant.\n\nMy take: this deserves a serious referee, because the underlying idea is publishable and the flaws are fixable. But I would not accept it in current form. The authors need to correct or clarify Eq. (1), recompute the LAP/SSP with the right exponents, and confirm the radial trends and the composition inference survive. If they can show that, the paper would be a solid, incremental contribution to the air-shower composition toolbox.","headline":"A useful extension of local-age methods to muon LDFs, but the exponent swap in Eq. (1) biases the central observable and must be fixed before the results can be trusted.","tokens_in":15046,"tokens_out":3985,"would_cite":false,"duration_ms":35550,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["96.50.sd","95.75.z","96.50.S"],"model":"deepseek-v4-flash","headline":"The paper claims that two pointwise shape parameters of electron and muon lateral density profiles—local age and segmented slope—show cosmic-ray composition changing gradually from proton- to iron-dominated across the knee.","keywords":["cosmic-ray composition","extensive air showers","lateral density distribution","local age parameter","segmented slope parameter","KASCADE","knee of cosmic-ray spectrum","NKG function"],"falsifier":"Simulate showers with a known single primary mass and a known NKG age, then apply the paper's exact equations (1) and (3) to the simulated densities; if the output local age differs from the input age by the radius-dependent term $2.5/(1+2x)$, the KASCADE trends in figures 4–6 cannot be interpreted as measuring true shower age, and the composition conclusion would need to be re-derived with the corrected formula.","tokens_in":13931,"feed_emoji":"🌌","tokens_out":8887,"duration_ms":86585,"temperature":0.7,"pith_summary":"This paper argues that the mass composition of cosmic rays around the knee of the energy spectrum can be read off from the radial shapes of the lateral density profiles of air-shower electrons and muons, without needing a full reconstruction of each shower. Using proton- and iron-initiated Monte Carlo showers together with the KASCADE experiment's public electron and muon density data, the authors propose two pointwise shape estimators—the local age parameter for the electron (NKG) profile and the segmented slope parameter for the muon (Greisen) profile—and show that their mean extrema (minimum local age near 44 m, maximum segmented slope near 71 m) separate light and heavy primaries. Comparing simulation to data, they report a gradual change from a proton-dominated to an iron-dominated composition as shower size and muon size increase across the knee. If this is right, the lateral shapes of particle densities provide a simple, model-light handle on cosmic-ray composition in an energy region where the transition from galactic to extragalactic sources is still debated.","feed_headline":"Shower shapes trace cosmic rays turning heavier at the knee","feed_subtitle":"Local-age and slope data from KASCADE support a gradual proton-to-iron shift across the knee.","key_machinery":"The central objects are the local age parameter (LAP) and segmented slope parameter (SSP), pointwise shape estimators computed from density ratios at two core distances. For any lateral distribution $f(r)$, equation (1) defines $$a_{\\rm local}=\\frac{\\ln(F_{ij}X_{ij}^{\\alpha_1}Y_{ij}^{\\alpha_2})}{\\ln(X_{ij}Y_{ij})}$$ with $F_{ij}=\\rho(r_i)/\\rho(r_j)$, $X_{ij}=r_i/r_j$, and $Y_{ij}=(r_i/R+1)/(r_j/R+1)$. For the NKG electron form $f(r)=C(r/R)^{s-2}(1+r/R)^{s-4.5}$ the authors take $\\alpha_1=4.5$, $\\alpha_2=2$ and call the result the LAP; for the Greisen-type muon form $f_\\mu(r)=A(r/R_G)^{-\\beta}(1+r/R_G)^{-2.5}$ they take $\\alpha_1=0$, $\\alpha_2=2.5$ and call it the SSP. These estimators are meant to convert local density shapes into the age or slope of the cascade at a given radius, and the analysis uses their extrema, the minimum LAP near 44 m and the maximum SSP near 71 m, as mass-sensitive observables.","core_discovery":"The central claim is that the local age parameter and the segmented slope parameter, estimated directly from density ratios at neighboring core distances, are mass-sensitive observables and that their averages over the KASCADE data track a composition transition around the knee. The paper's Figure 4 shows that the mean minimum local ages of electrons and muons from observed events fall between the proton and iron simulation bands and move from the proton side toward the iron side with increasing shower size and muon size; Figures 5 and 6 make the same point for muon local age and for the maximum segmented slope. The authors interpret this as evidence that the cosmic-ray composition gradually shifts from lighter (proton-dominated) to heavier (iron-dominated) nuclei as primary energy increases through the knee, consistent with the earlier KASCADE result based on the electron-size versus muon-size correlation. They also report that the radial variation of the local age and segmented slope is approximately independent of shower size, which they call scaling behavior.","pith_inferences":["A direct algebra test of equation (1) with the paper's own NKG exponents gives $s_{\\rm local}=s+2.5/(1+2x)$ for a pure NKG profile rather than the true age $s$, so the 'minimum at 44 m' may partly reflect this offset; a corrected estimator would need to be tested before the composition conclusion is taken quantitatively.","A stronger test would apply the same LAP and SSP recipe to Monte Carlo showers with known mixed compositions (helium, CNO, and silicon as well as proton and iron) and check whether the mean extrema still separate mass groups with the same ranking.","The 44 m and 71 m feature locations may depend on detector spacing and on the chosen Moliere and Greisen radii; repeating the analysis at a denser array would show whether the features are a property of the shower or of the sampling geometry."],"forward_implications":["The mean minimum local age near 44 m and the mean maximum segmented slope near 71 m can serve as mass discriminators for experiments that measure lateral densities but not the depth of shower maximum.","Because the radial LAP and SSP curves are nearly size-independent, composition might be inferred from a single shape extremum per event, reducing the need for precise primary-energy reconstruction.","Extending the local-age analysis from electrons to muons gives a second, independent composition probe that is less affected by atmospheric attenuation.","If the reported composition trend is correct, the knee of the cosmic-ray spectrum is not a sudden switch but a gradual light-to-heavy transition, in line with rigidity-dependent cutoffs of galactic sources."],"supporting_citations":[{"why":"Supplies the CORSIKA Monte Carlo generator used to produce proton and iron showers for the simulated LAP and SSP templates.","marker":"[4]"},{"why":"Provides the public KASCADE electron and muon density data that the paper compares against the simulations.","marker":"[5]"},{"why":"Introduced the local-age and scaling concept for electron lateral densities that this paper extends to muons.","marker":"[2]"},{"why":"Establishes the local shape parameter framework and the NKG-based formulation used in equation (1).","marker":"[3]"},{"why":"Gives the NKG structure function that underlies the LAP definition for electrons.","marker":"[14]"},{"why":"Supplies the lateral energy correction functions used to convert KASCADE detector energy deposits into densities.","marker":"[25]"},{"why":"Provides the earlier KASCADE Ne-Nmu composition result that the paper's LAP and SSP trend is cross-checked against.","marker":"[27]"}],"fun_headline_variants":["Local shower age reveals cosmic-ray mass shift at knee","New method uses shower slopes to weigh cosmic rays","KASCADE data link cosmic-ray composition change to knee","Shower age and slope pin down cosmic-ray mix at knee","From protons to iron: cosmic-ray mix shifts at knee"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing step is equation (1)'s claim that $s_{\\rm local}$ recovered from a two-point density ratio equals the true shower age for the standard electron profile; direct algebra with the paper's chosen exponents gives $s_{\\rm local}=s+2.5/(1+2x)$ instead, so the reported 'minimum at 44 m' and the composition trend built on it may be measuring this offset.","fun_headline_variants_meta":{"raw":{"variants":["Local shower age reveals cosmic-ray mass shift at knee","New method uses shower slopes to weigh cosmic rays","KASCADE data link cosmic-ray composition change to knee","Shower age and slope pin down cosmic-ray mix at knee","From protons to iron: cosmic-ray mix shifts at knee"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000179,"raw_usage":{"total_tokens":1261,"prompt_tokens":868,"completion_tokens":393,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":484,"completion_tokens_details":{"reasoning_tokens":315}},"tokens_in":484,"tokens_out":393,"duration_ms":4289,"temperature":1.0,"reasoning_tokens":315,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:50:40.546942+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate showers with a known single primary mass and a known NKG age, then apply the paper's exact equations (1) and (3) to the simulated densities; if the output local age differs from the input age by the radius-dependent term $2.5/(1+2x)$, the KASCADE trends in figures 4–6 cannot be interpreted as measuring true shower age, and the composition conclusion would need to be re-derived with the corrected formula.","supporting_citations":[{"cited_title":"The characteristic feature of SSP versus the core distance curves satisfy a complete opposite configuration compared to the radial variation of the LAP","cited_arxiv_id":null,"evidence_quote":"Supplies the CORSIKA Monte Carlo generator used to produce proton and iron showers for the simulated LAP and SSP templates."},{"cited_title":"Haungs et al., European Phys.J","cited_arxiv_id":null,"evidence_quote":"Provides the public KASCADE electron and muon density data that the paper compares against the simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduced the local-age and scaling concept for electron lateral densities that this paper extends to muons."},{"cited_title":"The frequency distributions of the minimum values of slocal(r) for two different Moli´ ere radii possess different mean val- ues","cited_arxiv_id":null,"evidence_quote":"Establishes the local shape parameter framework and the NKG-based formulation used in equation (1)."},{"cited_title":"Nishimura, and K","cited_arxiv_id":null,"evidence_quote":"Gives the NKG structure function that underlies the LAP definition for electrons."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the lateral energy correction functions used to convert KASCADE detector energy deposits into densities."},{"cited_title":"Bleicher, et al., J","cited_arxiv_id":null,"evidence_quote":"Provides the earlier KASCADE Ne-Nmu composition result that the paper's LAP and SSP trend is cross-checked against."}],"review_version":1}