{"id":"b9ba3b57-e383-4c78-928d-f9b1ed0a754a","arxiv_id":"1909.01406","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"IceCube's 7-year, half-trillion-event muon sample shows a correlation coefficient of 0.75 between muon rate and effective temperature, along with a newly visible nonlinear hysteresis.","lead":"IceCube has measured how the rate of atmospheric muons arriving at the South Pole varies with stratospheric temperature over seven years. The extremely high event rate reveals a subtle nonlinearity in the rate-temperature relationship and sharp short-term jumps.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed hysteresis in muon rate vs effective temperature lacks a quantitative significance test; without one, the 'first observation' of nonlinearity is not established.","rationale":"The reader's verdict is CONDITIONAL, citing model dependence of the Teff weighting and missing systematics. I agree that model dependence is a real concern, but I identify a more direct, load-bearing gap: the paper's central novel claim, the observed hysteresis, is never subjected to a statistical significance test. This is not just a missing systematic; it is missing evidence for the existence of the effect itself. If the apparent loop in Fig. 3(a) is consistent with a linear relation plus correlated noise, the 'first observation' claim collapses regardless of hadronic model choices. The concrete tests I propose would settle this directly. Since the reader's verdict already conditions acceptance on addressing quantitative gaps, my concern does not change the verdict category: it remains CONDITIONAL. I therefore set verdict_should_be to UNCHANGED, with partial agreement because the focus shifts from model dependence to the missing significance test, though both point to the need for additional analysis before the claim can be accepted.","tokens_in":5617,"tokens_out":8355,"duration_ms":86055,"concrete_test":"Using the 2012 daily IceCube muon rates and daily Teff values (as in Fig. 3(a)), perform a least-squares fit of the linear relation delta_R/<R> = alpha_T * delta_Teff/<Teff>. Compute the residuals and test whether they are significantly correlated with day-of-year after removing the seasonal Teff trend (e.g., Spearman rank correlation). Alternatively, define the hysteresis loop area by fitting an ellipse to the daily points in the (delta_Teff/<Teff>, delta_R/<R>) plane and estimate its uncertainty via bootstrap resampling of days. As a second check, bin the data by Teff and compare the mean rate in overlapping Teff bins between austral spring (September-November) and autumn (March-May) using a two-sample t-test or permutation test.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central novelty is the hysteresis in Fig. 3: the daily muon rate versus effective temperature for 2012 is presented as showing a nonlinear, loop-like relation, and the abstract/summary claims this is observed 'for the first time.' Yet the paper never quantifies this nonlinearity. No error bars are shown on the daily points in Fig. 3(a), the linear-fit slope alpha_T = 0.75 is quoted without an uncertainty, and no metric for the hysteresis (e.g., loop area, spring-versus-autumn rate difference at equal Teff) is defined or tested against the null hypothesis of a purely linear relation. The authors note that statistical fluctuations are at the 1e-4 level, which implies the daily points are very precise, but precision alone does not establish that the apparent loop is significant; correlated atmospheric or detector effects could produce a similar pattern. The model calculation in Fig. 3(b) shows a 'qualitatively similar, though slightly smaller' hysteresis, but this is not a quantitative comparison, and since both measured and calculated rates are plotted against the same model-derived Teff, it does not independently validate the effect. In addition, Sec. 4 reports a factor-of-two offset in absolute rate and a 2% discrepancy in seasonal amplitude, indicating the model is not fully validated; this reinforces that the claimed nonlinearity needs direct statistical support.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a seven-year (IC86, 2011-2018) measurement of the seasonal variation of the downgoing muon rate in IceCube, using roughly half a trillion muon events. The measured daily rate is correlated with an effective atmospheric temperature computed from AIRS temperature profiles and a muon production spectrum based on Sibyll 2.3c and the H3a nucleon flux model. The central claims are: (1) a linear correlation coefficient alpha_T = 0.75 for 2012, consistent with expectations for approximately TeV muons; (2) the first observation of a nonlinear, hysteresis-like relation between muon rate and effective temperature, most pronounced during the austral spring; and (3) a model calculation using Eqs. (3.1)-(3.5) that reproduces the short-term rate features and exhibits qualitatively similar but smaller hysteresis. The paper also discusses implications for atmospheric neutrino seasonal variations and the kaon-to-pion ratio.","tokens_in":5878,"tokens_out":1686,"duration_ms":17769,"significance":"If the central claim holds, this is the first high-precision demonstration that the seasonal response of atmospheric muons is not a single linear function of effective temperature, which would be an important constraint on atmospheric cascade models and on the interpretation of seasonal variations of atmospheric neutrinos. The paper's strengths include a very large data set (2.15 kHz trigger rate, half a trillion events) that reduces statistical errors to the 10^-4 level, and a transparent, physically motivated calculation of the effective temperature from first-principles production spectra. The analysis is also honest about known shortcomings: it explicitly states the absolute calculated rate is off by a factor of two and that the calculated seasonal amplitude is about 2% larger than observed. However, the headline novelty, the claimed hysteresis, is not quantified or tested for statistical significance, and the measured alpha_T is quoted without an uncertainty. Because the effective-temperature definition itself depends on the production model, the nonlinearity claim needs a direct statistical test against a linear null hypothesis before it can be considered established.","major_comments":[{"comment":"The claim that the hysteresis is observed 'for the first time' is the paper's central novelty, yet the paper provides no quantitative significance test of the nonlinearity. Fig. 3(a) shows daily points without error bars, the slope is quoted as alpha_T = 0.75 with no uncertainty, and no metric such as loop area or spring-versus-autumn rate difference at equal Teff is defined or compared with the null hypothesis of a purely linear relation. Since the statistical fluctuations are stated to be at the 10^-4 level, the authors should be able to provide a straightforward chi-square or likelihood test; without it, the apparent loop could be due to correlated atmospheric or detector effects.","section":"Section 4, Fig. 3"},{"comment":"The claimed hysteresis is model-dependent because both the measured and calculated rates are plotted against the same model-derived Teff from Eq. (1.2), weighted by the Sibyll 2.3c/H3a production spectrum. If the production spectrum is wrong, the Teff weights are wrong and the loop could be an artifact of the weighting rather than a physical property of the atmospheric response. The authors state the model comparison is 'qualitatively similar, though slightly smaller,' but this is not a quantitative comparison. The factor-of-two absolute rate offset and the 2% seasonal-amplitude discrepancy reported in Section 4 further indicate that the model is not fully validated, which strengthens the need for a model-robustness test of the nonlinearity, e.g., by recomputing Teff with alternative hadronic interaction models or by checking whether the loop survives when using a purely empirical temperature weighting.","section":"Section 4, Fig. 3(b) and Eqs. (1.2), (3.1)-(3.5)"},{"comment":"The correlation coefficient alpha_T is defined by a straight-line fit in Eq. (1.1), but the same section claims the relation is nonlinear. The paper should clarify whether alpha_T is a linear fit performed over the full-year data (including the hysteresis loop) or over a restricted period, and should report the fit uncertainty and the residuals from the linear fit. Without this, the reader cannot assess whether alpha_T = 0.75 is consistent with the model prediction alpha_T ~ 0.84, especially given the stated 2% discrepancy in seasonal amplitude.","section":"Eq. (1.1) and Section 4, Fig. 3(a)"}],"minor_comments":[{"comment":"The text states 'For T = 220° K, επ = 115 GeV and εK = 857 GeV,' but the equation uses a temperature-dependent epsilon; it would be helpful to state explicitly that these values are computed at that reference temperature and to give the numerical values of the other constants used.","section":"Section 1, Eq. (1.3)"},{"comment":"The caption 'Comparison of measured muon rate with calculated rate for 2012' does not state that the calculated rate is normalized to the observed rate; this is mentioned in the text but should be in the caption to avoid misinterpretation.","section":"Section 4, Fig. 2 caption"},{"comment":"The sentence 'The important region for downward muons is near the vertical cosθ≥ 0.5' contains a minor grammar issue; it should read 'is at cosθ ≥ 0.5' or 'is the region cosθ ≥ 0.5.'","section":"Section 5, last paragraph"},{"comment":"Reference [9] is formatted as 'PoS(ICRC2017)301 (2018). [35,301(2017)]'; the journal reference appears incomplete. Please provide the full arXiv identifier or published DOI for the Sibyll 2.3c paper.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"This is a conference proceedings paper (ICRC 2019), and the scope is appropriate for that venue. The central measurement and basic correlation are sound, but the novelty claim of hysteresis needs a quantitative statistical test. The paper should be encouraged to add such a test, or to soften the 'first observation' claim if a test cannot be provided in the proceedings format. I do not see grounds for rejection, but the current form has a load-bearing unsupported claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe thing to know: this is a solid IceCube measurement of seasonal muon-rate variation, with half a trillion events over seven years, and the authors are candid about where their model fails. But the advertised novelty—the \"first observation\" of a hysteresis in rate vs effective temperature—is not actually established, because the paper never quantifies the loop. That is the difference between a useful progress report and a finished result.\n\nThe genuinely good parts: the daily rate data, the effective-temperature weighting built on Sibyll 2.3c and H3a, the measured alpha_T = 0.75 for 2012, the reproduction of the sharp October rate jump (5.4% measured vs 5.9% calculated), and the angular comparison of muon vs neutrino correlation as a handle on the kaon-to-pion ratio. The paper also sits honestly in the literature, citing MINOS and earlier IceCube work. The factor-of-two absolute normalization offset is admitted up front, and the authors say the seasonal amplitude is overpredicted by about 2%; that is the kind of limitation you want to see stated.\n\nThe soft spots are real but specific. First, the hysteresis figure has no error bars, no loop-area metric, and no test against a linear null. The claim that statistical fluctuations are at the 10^-4 level does not by itself prove the loop is significant; daily rates are correlated with weather, and systematic effects could mimic the pattern. Second, Fig. 3(b) is calculated with the same model-derived Teff, so it does not independently validate the nonlinearity. Third, there is no systematic uncertainty budget, so alpha_T = 0.75 floats without a sigma. These are all fixable, but they are not cosmetic; they are exactly what a referee would need before the \"first time\" claim can stand.\n\nOverall, the central seasonal correlation is robust; the nonlinearity, as presented, is not. This paper is for atmospheric-muon specialists and people using IceCube neutrino seasonal corrections. It deserves peer review—serious review, not a desk reject—but with the expectation that the hysteresis claim gets a quantitative treatment and the systematics get written down.\n\nRecommendation: send it to review, and tell the authors their strongest result is the long-baseline rate–temperature correlation, not the loop.","headline":"A credible long-baseline measurement of seasonal muon-rate variation with honest model caveats, but the advertised 'first observation' of hysteresis is not quantified and therefore not established.","tokens_in":6398,"tokens_out":1826,"would_cite":true,"duration_ms":18778,"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":"IceCube data reveal, for the first time, that the seasonal muon rate is nonlinear in atmospheric temperature, forming a hysteresis loop with measured correlation coefficient 0.75 in 2012.","keywords":["seasonal variation","atmospheric muons","effective temperature","hysteresis","IceCube","pion decay","kaon decay","primary cosmic-ray spectrum"],"falsifier":"Recompute the 2012 daily effective temperatures with a substantially different hadronic interaction model, or with pion and kaon production moments varied within their uncertainties, and re-plot $\\delta R/\\langle R\\rangle$ against $\\delta T_\\mathrm{eff}/\\langle T_\\mathrm{eff}\\rangle$. If the hysteresis loop disappears, changes orientation, or changes size by more than the statistical precision of the IceCube rate, the claimed nonlinearity is an artifact of the production-spectrum weighting; if the loop persists in roughly the same shape, the atmospheric temperature profile is the cause.","tokens_in":5435,"feed_emoji":"🌡️","tokens_out":9753,"duration_ms":83208,"temperature":0.7,"pith_summary":"Seven years of full IceCube data, more than half a trillion muon events at an average trigger rate of 2.15 kHz, show that the seasonal variation of the atmospheric muon rate is not a single straight-line response to atmospheric temperature. The paper reports the first observation of a nonlinear, hysteresis-like relation between the measured muon rate and a model-weighted effective temperature, with a measured correlation coefficient $\\alpha_T = 0.75$ for 2012. A calculation using the same effective-temperature weighting reproduces the hysteresis qualitatively, with $\\alpha_T \\simeq 0.84$, and also matches sudden short-term rate jumps such as the 5.4% increase in five days in October 2012. The result matters because it indicates that the atmospheric response to temperature carries information about the vertical temperature profile, not just a single slope, and it connects directly to seasonal variations of atmospheric neutrinos in the same energy range.","feed_headline":"Seasonal muon rate loops around temperature, not a straight line","feed_subtitle":"Seven years and half a trillion muon events reveal a nonlinear rate-temperature relation, with slope 0.75 in 2012.","key_machinery":"The load-bearing object is the model-weighted effective temperature $T_\\mathrm{eff}$ of Eq. (1.2), which weights each atmospheric temperature $T(X)$ by the muon production spectrum $P_\\mu(E_\\mu,\\theta,X)$ along the muon trajectory and by the detector effective area. Temperature enters the physics only through the critical energies $\\varepsilon_\\pi \\simeq 115$ GeV and $\\varepsilon_K \\simeq 857$ GeV at $T=220$ K, Eq. (1.3), which control the competition between decay and re-interaction of charged pions and kaons. The production spectrum combines a low-energy pion and kaon decay form with a high-energy form, using the Sibyll 2.3c hadronic interaction model and the H3a nucleon flux model, and the correlation coefficient $\\alpha_T$ is defined as the slope of $\\delta R/\\langle R\\rangle$ versus $\\delta T_\\mathrm{eff}/\\langle T_\\mathrm{eff}\\rangle$. That weighting is what turns a one-dimensional temperature into a prediction for the muon rate, and it is the quantity whose nonlinearity the data reveal.","core_discovery":"The central discovery is that the fractional variation of the IceCube muon rate, plotted against the fractional variation of the model-weighted effective temperature for 2012, does not scatter around a straight line. It forms a loop, or hysteresis, with a fitted correlation coefficient $\\alpha_T = 0.75$, close to the value expected for the roughly TeV muons that dominate the InIce-SMT8 trigger. A forward calculation using the analytic production spectrum and the same daily temperature profiles reproduces a similar but slightly smaller loop, with $\\alpha_T \\simeq 0.84$, and the paper attributes the loop to the austral spring: the upper atmosphere warms quickly while deeper air remains cold, so the same $T_\\mathrm{eff}$ can correspond to different rates. The paper also notes that the calculated absolute rate is a factor of two higher than observed, which it attributes to the normalization of the primary nucleon spectrum, and that the calculated annual amplitude is about 2% larger than measured.","pith_inferences":["If the hysteresis is truly atmospheric, replacing the Sibyll 2.3c production moments with a different hadronic model should move the curve in Fig. 3 slightly but should not erase the loop; erasure would identify the nonlinearity as an artifact of the weighting scheme rather than a property of the atmosphere.","Because the austral-spring profile drives the loop, the same analysis applied to other years should show a loop of similar shape only in years whose October temperature profiles share the same upper-stratosphere warming; a year with a different profile would be a natural test of the mechanism.","The factor-two normalization offset suggests that absolute atmospheric-muon rate measurements at IceCube could serve as a check on primary cosmic-ray flux models, provided the detector effective area and the contribution of coincident or multiple muons are known to comparable precision."],"forward_implications":["The seasonal response of atmospheric muons cannot be summarized by a single $\\alpha_T$; the rate versus effective temperature relation has a loop structure that encodes the vertical temperature profile.","Short-term rate jumps of several percent on day-to-week time scales are reproduced by the same effective-temperature calculation, so the production-spectrum formalism can be used to predict daily muon rates from satellite temperature profiles.","The factor-two discrepancy between calculated and observed absolute rate, together with the 0.75 versus 0.84 difference in $\\alpha_T$, points to the normalization of the primary cosmic-ray nucleon spectrum as the main missing ingredient, not the temperature weighting.","Because the same formalism applies to atmospheric muon neutrinos, with kaon decay becoming the dominant neutrino source above roughly 100 GeV, simultaneous seasonal measurements of muons and neutrinos should be sensitive to the kaon-to-pion ratio."],"supporting_citations":[{"why":"Defines the standard correlation-coefficient fit of rate versus effective temperature and provides the MINOS far-detector context for alpha_T.","marker":"[2]"},{"why":"Earlier IceCube analysis of seasonal muon variations during detector construction, the baseline the full-detector IC86 measurement extends.","marker":"[6]"},{"why":"Supplies the daily South Pole atmospheric temperature profiles used to compute the effective temperature and the correlation.","marker":"[7]"},{"why":"Provides the analytic muon production spectrum formalism for low- and high-energy pion and kaon decay used in the weighting.","marker":"[8]"},{"why":"Supplies the Sibyll 2.3c hadronic interaction model used for the pion and kaon production moments in the production spectrum.","marker":"[9]"},{"why":"Supplies the H3a model for the nucleon flux that sets the primary spectrum and the normalization of the calculated rate.","marker":"[10]"}],"fun_headline_variants":["Muon rate vs temperature traces a loop in IceCube","IceCube records seasonal muon-temperature hysteresis","Half a trillion muons show rate-temperature loop","Seasonal muon rates loop with effective temperature","IceCube's muon rate seasonal loop: same temp, different rates"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result assumes that the particle-physics model of how cosmic-ray nuclei produce pions and kaons in the atmosphere — implemented here with the Sibyll 2.3c interaction model and the H3a nucleon flux — is accurate enough that the effective-temperature weighting reflects the real atmosphere rather than imposing the observed curvature on the rate-temperature relation.","fun_headline_variants_meta":{"raw":{"variants":["Muon rate vs temperature traces a loop in IceCube","IceCube records seasonal muon-temperature hysteresis","Half a trillion muons show rate-temperature loop","Seasonal muon rates loop with effective temperature","IceCube's muon rate seasonal loop: same temp, different rates"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000231,"raw_usage":{"total_tokens":1430,"prompt_tokens":835,"completion_tokens":595,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":451,"completion_tokens_details":{"reasoning_tokens":518}},"tokens_in":451,"tokens_out":595,"duration_ms":5931,"temperature":1.0,"reasoning_tokens":518,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:17:54.396441+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the 2012 daily effective temperatures with a substantially different hadronic interaction model, or with pion and kaon production moments varied within their uncertainties, and re-plot $\\delta R/\\langle R\\rangle$ against $\\delta T_\\mathrm{eff}/\\langle T_\\mathrm{eff}\\rangle$. If the hysteresis loop disappears, changes orientation, or changes size by more than the statistical precision of the IceCube rate, the claimed nonlinearity is an artifact of the production-spectrum weighting; if the loop persists in roughly the same shape, the atmospheric temperature profile is the cause.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the standard correlation-coefficient fit of rate versus effective temperature and provides the MINOS far-detector context for alpha_T."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier IceCube analysis of seasonal muon variations during detector construction, the baseline the full-detector IC86 measurement extends."},{"cited_title":"Desiati, T","cited_arxiv_id":null,"evidence_quote":"Supplies the daily South Pole atmospheric temperature profiles used to compute the effective temperature and the correlation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the analytic muon production spectrum formalism for low- and high-energy pion and kaon decay used in the weighting."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Sibyll 2.3c hadronic interaction model used for the pion and kaon production moments in the production spectrum."},{"cited_title":"Riehn, H","cited_arxiv_id":null,"evidence_quote":"Supplies the H3a model for the nucleon flux that sets the primary spectrum and the normalization of the calculated rate."}],"review_version":1}