{"id":"b3a8ba5e-6563-4c9e-9e19-476b283aa1e4","arxiv_id":"1908.07695","paper_version":2,"verdict":"REJECT","confidence":"HIGH","novelty_score":2.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"A prime-number lookup table, with energy treated as the prime index and the prime value scaled by fitted constants, is proposed as an estimator of neutrino-nucleon and proton-proton total cross sections.","lead":"This paper proposes that neutrino-nucleon and proton-proton collision cross sections can be estimated by matching particle energies to prime numbers and multiplying by fixed scaling factors. A physicist might read it because it claims a simple number-theoretic formula matches IceCube neutrino data where the Standard Model appears to under-predict events.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"PeV predictions are not out-of-sample: the normalization constants are fitted to the same data, so the IceCube agreement is weak evidence for the prime-index claim.","rationale":"The reader's rejection is justified. The central claim is that prime-index parametrization provides an accurate estimate over many decades, but the quantitative content reduces to a fitted constant multiplying an E ln E curve. The paper itself concedes that the estimates are empirically driven and that the normalization factors are a first-step approximation. No formal verification or independent validation is offered. The key weakness is the absence of any out-of-sample test: the constants are calibrated on data that overlap the validation range, and the high-energy cascade normalization is explicitly chosen to match the IceCube data being compared. The muon-neutrino comparison, while using the 0.48 constant, is a low-statistics event count where the astrophysical flux and cross-section are degenerate; scaling the SM expectation by a cross-section ratio while holding the flux fixed is not a clean test. A simple calibration-split test would settle whether the constants extrapolate. Because no such test exists in the paper, the REJECT verdict remains appropriate.","tokens_in":5808,"tokens_out":12894,"duration_ms":123817,"concrete_test":"Calibrate c_nu and c_antinu only on data below 100 GeV; use those constants to predict the cross-section data above 100 GeV and the IceCube >0.5 PeV event count (with the same flux model as the paper). Report the validation residuals and whether the predicted IceCube count remains near 11. A robust parametrization must reproduce the high-energy data within the quoted uncertainties; if the constants shift materially or the event count moves outside the error band, the central extrapolation claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The prime-index model contains no physical parameters; the constants 0.70, 0.26, 0.48, and 0.69 are free normalization factors adjusted to data. Since p_i ~ i ln i, the model is functionally sigma(E) ~ c E ln E; the only load-bearing content is that c stays constant over six decades. The paper provides no out-of-sample test for this. The 0.69 constant used for the cascade data is explicitly called a 'first step approximation', meaning the Fig. 3 agreement is circular. For the muon-neutrino Table 3, the expected count is obtained by scaling the SM rate by the cross-section ratio. But the SM rate is computed with the same IceCube observation constraining the flux; fixing the flux and inflating the cross-section double-counts the events used to set the flux. Even taken at face value, 11 +/- 3.3 vs 9 +/- 3 is a roughly 0.6-sigma fluctuation, not 'strong evidence'. The +/- 3.3 is Poisson counting only; no uncertainty on the fitted constants is propagated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes an empirical mapping between prime numbers and their indices and neutrino-nucleon cross sections: taking the prime's index as the neutrino energy in MeV, the cross section in units of 10^-42 cm^2 is 0.70 times the prime for neutrinos and 0.26 times the prime for antineutrinos. The author claims this reproduces charged-current cross-section data from MeV to PeV energies. For IceCube events above 0.5 PeV, the paper reports 11 ± 3.3 expected events from the prime-index method versus 9 ± 3 observed events, compared with 1 ± 1 from the Standard Model. A normalization of 0.48, the average of 0.70 and 0.26, is used for muon neutrinos, and 0.69 is used for cascade data. The paper further claims that a twin-prime companion ratio, multiplied by 0.25, reproduces total pp cross sections and 'explains' the (ln s)^2 growth first proposed by Heisenberg. The text itself states that the estimates are empirically driven and are not intended as replacements for physics-based approaches.","tokens_in":6035,"tokens_out":7425,"duration_ms":160950,"significance":"If the prime-index mapping were a genuine out-of-sample predictor across six decades, it would be a striking empirical curiosity. As it stands, however, the five normalization constants (0.70, 0.26, 0.48, 0.69, and 0.25) are all fixed to the data they are later said to reproduce, so the agreement shown in Figs. 1-4 is not independent evidence. No uncertainties are assigned to these constants, no goodness-of-fit statistic is reported, and the IceCube comparison has overlapping error bars and depends on flux assumptions. The pp section builds on the Hardy-Littlewood twin-prime conjecture rather than an established theorem. The paper is honest about its empirical character and calls the cascade normalization a 'first step approximation,' but these caveats point to the central weakness: the claims reduce to curve fitting with free parameters. The only concrete forward-looking statement is a qualitative appeal to a future 10 km^3 detector; no quantitative prediction or decision threshold is given.","major_comments":[{"comment":"The claimed six-decade agreement is by construction. The model output is sigma = 0.70 * p_n for neutrinos and 0.26 * p_n for antineutrinos, with n equal to the energy in MeV and p_n the n-th prime; the constants 0.70 and 0.26 are chosen so that the prime sequence lines up with the data in Fig. 1, and no derivation fixes them. Since p_n ~ n ln n, the functional form is sigma(E) = c E ln E, and with one free constant c it is unsurprising that a smooth curve can be matched to a monotone data set. The paper reports no chi-square, no residuals, and no uncertainty on c, so the statement in Sec. 3 that the method provides a 'quick and accurate estimate ... over many decades of energy scales' is quantitatively unsupported.","section":"Sec. 3, Eq. (6), Fig. 1"},{"comment":"The IceCube event comparison is not an out-of-sample test. The normalization 0.48 used for the muon-neutrino rate is the average of the two lower-energy fitted constants, 0.70 and 0.26, so it inherits the earlier fits. The expected rate is obtained by scaling the SM rate by the cross-section ratio; to the extent that the SM rate is normalized using the same IceCube event sample to determine the astrophysical flux, the same events are used both to fix the flux and to test the cross-section enhancement. Even taken at face value, 11 ± 3.3 versus 9 ± 3 is a sub-one-sigma difference, not 'strong evidence' as claimed in Sec. 4, and the stated ±3.3 is only Poisson counting with no propagation of the uncertainty in the fitted normalization constants.","section":"Sec. 4, Table 3"},{"comment":"The cascade-data comparison is admitted to be circular. The text states that the 0.69 normalization factor 'should be only considered as a first step approximation' until more IceCube data are available, meaning the constant is adjusted to make the curve match Fig. 3. Agreement obtained by fixing the overall normalization in this way carries no evidential weight, and the relationship between the factor 0.7 in Eq. (7) and the 0.69 used for Fig. 3 is not explained.","section":"Sec. 4, Eq. (7), Fig. 3"},{"comment":"The pp 'explanation' is not an explanation. Equation (10) is the Hardy-Littlewood twin-prime asymptotic, which is a conjecture, not a theorem; the Brun bound cited immediately before it gives only an upper bound of the form C N/(ln N)^2, not an asymptotic equality. The (ln s)^2 dependence is therefore assumed rather than derived, and the TPC ratio is multiplied by 0.25 'to normalize it to the experimental data' (Fig. 4 caption). Thus the level of agreement in Fig. 4 is set by a fitted constant, and the paper does not explain Heisenberg's parametrization.","section":"Sec. 5, Eq. (10), Fig. 4"}],"minor_comments":[{"comment":"The abstract contains a typo ('prim e numbers'), and the manuscript uses a nonstandard apostrophe in 'Gauss′s'.","section":"Abstract"},{"comment":"The statement that Eq. (6) introduces 'a 22% error at low energies and about 7% error at high energies' refers to the accuracy of the prime-index approximation itself, not to a comparison with measured cross sections; the two types of error should not be conflated.","section":"Sec. 3"},{"comment":"Table 2 would be clearer if it stated which neutrino species (nu_mu and antinu_mu) and which target (isoscalar nucleon) are assumed, and whether radiative corrections are included in the SM cross sections.","section":"Sec. 4, Table 2"},{"comment":"In Eq. (10), the phrase 'x is a pair of twin primes' is imprecise: x is the upper limit of the counting function pi_2(x), not a pair of primes.","section":"Sec. 5, Eq. (10)"},{"comment":"The term 'Twin Prime Companion (TPC)' is used without a formal definition; the text says only that it is the composite sandwiched between a pair of twin primes, so a reader cannot reproduce the ratio plotted in Fig. 4.","section":"Sec. 5"},{"comment":"Figs. 2 and 3 are reproductions of published IceCube figures with curves superimposed; without the underlying numerical data or a description of how the published curves were digitized, readers cannot independently verify the comparisons.","section":"Figs. 2 and 3"},{"comment":"The text refers to a '10−km3 upgrade' to IceCube; the notation should be '10 km^3'.","section":"Sec. 3"}],"recommendation":"reject","confidential_remarks":"This manuscript is an empirical curve-fitting exercise. The central claims rest on normalization constants that are adjusted to the data being 'predicted'; no uncertainties are propagated; the IceCube comparison is not out-of-sample; and the pp explanation relies on a conjecture. I do not see a route to fixing these issues within the scope of the paper, which is why I recommend rejection. A substantially revised version would need to provide a genuine out-of-sample prediction with quantified uncertainties and a statistical comparison to the data."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a curve fit wearing a prime-number costume. The method is simple: take the prime-index relation p_i ~ i ln i, call i the neutrino energy in MeV, and multiply by a constant fitted to measured cross sections. That is the whole content. The claimed IceCube agreement is not evidence. The prediction of 11 ± 3.3 events versus the observed 9 ± 3 is a sub-1-sigma fluctuation, and the predicted rate is obtained by scaling the SM rate by the cross-section ratio without self-consistently adjusting the flux — a flux that was inferred from the same events under SM cross sections. That double-counts the events. The pp section is even weaker: it leans on the Hardy-Littlewood twin-prime conjecture, a conjecture rather than a theorem, and the factor 0.25 is explicitly chosen to match the data.\n\nWhat is new here is the label, not the math. The prime number theorem gives p_i ~ i ln i, so the model is sigma(E) ∝ E ln E with a fitted normalization. The paper is honest about being empirically driven and describes its normalization factors as a 'first step approximation.' The figures are readable, and the extrapolation to six decades is clearly shown. Credit for clarity, but not for physical content.\n\nThe central flaw is that a fitted curve is presented as a predictive method. No uncertainties are propagated from the fitted constants, no out-of-sample test is performed, and the high-energy 'prediction' is simply the same function with the same constants. The IceCube comparison does nothing to validate the model because the muon-neutrino constant 0.48 is the average of the two low-energy constants, and the statistical agreement is poor. The cascade normalization 0.69 is also acknowledged as a first-order approximation.\n\nThe pp section is a category error: the twin-prime conjecture, even if true, describes the distribution of twin primes; using it to 'explain' the (ln s)^2 behavior of pp cross sections through a fitted normalization is numerology, not explanation. The self-citations for the TPC generation are to an unpublished conference note and a self-published book, which don't add credibility.\n\nWho is this for? Possibly someone collecting examples of how arbitrary math can be made to fit data. It is not a paper that advances neutrino or hadron physics. I would not cite it, and I would not send it to peer review without expecting a quick rejection. A desk rejection is the right call.","headline":"A curve fit dressed as a prime-number discovery; the IceCube agreement is statistically weak and the pp 'explanation' leans on an unproven conjecture.","tokens_in":6566,"tokens_out":5726,"would_cite":false,"duration_ms":78801,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["13.15.+g","25.30.Pt","02.10.De"],"model":"deepseek-v4-flash","headline":"Prime-index lookup tables reproduce neutrino-nucleon cross sections over six decades, and a twin-prime ratio yields the pp rise.","keywords":["prime-index parametrization","neutrino-nucleon cross section","charged-current interactions","IceCube astrophysical neutrinos","twin primes","pp total cross section","(ln s)^2 growth","cross-section estimation"],"falsifier":"Measure or extract the neutrino-nucleon charged-current cross section at a single energy near 1 PeV; the method expects roughly ten times the Standard Model value at that energy, so a measurement near the Standard Model prediction would refute the high-energy extrapolation. In the same spirit, an updated IceCube sample with several times the current exposure that finds events above 0.5 PeV consistent with the Standard Model estimate rather than the $11\\pm3.3$ prediction would settle the claim.","tokens_in":5586,"feed_emoji":"⚛️","tokens_out":8796,"duration_ms":75620,"temperature":0.7,"pith_summary":"This paper proposes that the total charged-current neutrino-nucleon cross section at a given energy can be read from a two-column table of primes: set the neutrino energy in MeV to be the index of a prime, multiply that prime by 0.70 for neutrinos and by 0.26 for antineutrinos, and the result is the cross section in units of $10^{-42}\\,\\mathrm{cm}^2$. The author shows that this prescription follows accelerator data from MeV to PeV energies and claims it provides a quick estimate of interaction rates in neutrino experiments. At IceCube energies above 0.5 PeV, the method predicts $11\\pm3.3$ astrophysical muon-neutrino events against the observed $9\\pm3$, while the Standard Model estimate is $1\\pm1$. A companion twin-prime construction, in which the ratio of a twin-prime companion to its index, times 0.25, gives the $pp$ cross section in mb, is offered as an explanation of the $(\\ln s)^2$ growth of high-energy $pp$ cross sections. The paper is explicit that these estimates are empirical and not a replacement for physics-based calculations.","feed_headline":"Prime table reproduces neutrino cross sections from MeV to PeV","feed_subtitle":"A prime-index lookup with one fitted constant also predicts 11 IceCube events above 0.5 PeV, where the Standard Model predicts 1.","key_machinery":"The central object is the prime-index relation: the $n$th prime $p(n)$ as a function of its rank $n$, used as a lookup table in which $n$ is taken to be the particle energy in MeV (or GeV for $pp$) and $p(n)$ times a fitted coefficient gives the cross section. The paper also uses the twin-prime companion (TPC), the composite number sandwiched between a twin-prime pair, whose ratio to its index, times 0.25, gives the $pp$ cross section; it is tied to the twin-prime asymptotic $\\pi_2(x)\\sim 2C_2\\,x/(\\ln x)^2$, which supplies the $(\\ln x)^2$ factor. The normalization constants 0.70, 0.26, 0.48, 0.69, and 0.25 are what convert prime magnitudes into physical cross-section units, and they are fixed by matching available data.","core_discovery":"The central claim is that the sequence of primes, read through their positional indices, is a ready-made parametrization of neutrino-nucleon cross sections: with energy in MeV as the index and the prime itself as the unnormalized cross section, charged-current neutrino and antineutrino cross sections are obtained simply by multiplying by 0.70 and 0.26 respectively, in units of $10^{-42}\\,\\mathrm{cm}^2$. The paper asserts that this prime-index method reproduces measured neutrino and antineutrino total cross sections over six decades of energy and, applied to the published IceCube muon-neutrino data, predicts $11\\pm3.3$ events above 0.5 PeV compared with $9\\pm3$ observed and $1\\pm1$ from the Standard Model. It further claims that the total $pp$ cross section's $(\\ln s)^2$ rise is explained by the twin-prime companion ratio, with the index as energy in GeV and the ratio times 0.25 matching data, thereby tying the logarithmic-squared form to the twin-prime asymptotic.","pith_inferences":["Because the index is assigned to the energy only after choosing MeV as the unit, the same data could be fit with a different normalization if energies were expressed in GeV; the method's reach therefore depends on that unit choice plus the fitted constants, not on a derived dynamical law.","If a larger IceCube sample settles near the Standard Model rate, the constants would have to become energy-dependent, which would leave the low- and mid-energy interpolation intact but remove the high-energy extrapolation.","The same construction could be tested on neutral-current neutrino data: equation (7) fixes the CC/(CC+NC) ratio near 0.7, so a precise measurement of that ratio at PeV energies would be a separate check of the method.","The $pp$ explanation inherits the status of the twin-prime asymptotic, which is a conjecture rather than a theorem; a different true growth rate for twin-prime counts would decouple the $(\\ln s)^2$ form from the number-theoretic ratio."],"forward_implications":["Any neutrino experiment with a known energy spectrum can obtain cross-section estimates from a precomputed prime table without integrating structure functions.","For IceCube, the method favors a neutrino-nucleon cross section at PeV energies well above the Standard Model value, so a larger exposure should see an event rate above 0.5 PeV closer to 11 than to 1.","At ultra-high energies near $10^{12}$ GeV, the prime-index extrapolation gives cross sections roughly a million times larger than Standard Model predictions, implying strong Earth-absorption effects that would be absent in the Standard Model.","The $pp$ analysis asserts that the total $pp$ cross section's logarithmic-squared rise follows from the density of twin primes, connecting a particle-physics regularity to a number-theoretic asymptotic."],"supporting_citations":[{"why":"Provides the low-energy (MeV) neutrino cross-section approximation that the prime-index prescription must match in the lowest decade.","marker":"[5]"},{"why":"Supplies a second low-energy analytic approximation used as a comparison anchor in the few-hundred-MeV region.","marker":"[6]"},{"why":"Gives the measured linear energy dependence in the GeV-TeV range that fixes the overall slope the prime-index method reproduces.","marker":"[7]"},{"why":"States a standard high-energy approximation that the prime-index extrapolation is compared against at ultra-high energies.","marker":"[8]"},{"why":"Provides the digitized experimental neutrino and antineutrino cross-section data plotted in figure 1, the main data anchor for the method.","marker":"[9]"},{"why":"Reports the IceCube astrophysical muon-neutrino rates used to test the prime-index prediction at PeV energies.","marker":"[10]"},{"why":"Gives the Standard Model cross-section calculations used in table 2 as the baseline the prime-index values are compared with at high energy.","marker":"[11]"},{"why":"Introduces the $(\\ln s)^2$ growth of $pp$ cross sections that the twin-prime companion ratio is claimed to explain.","marker":"[13]"},{"why":"Supplies the twin-prime asymptotic $\\pi_2(x)\\sim 2C_2\\,x/(\\ln x)^2$ that produces the logarithmic-squared factor in the $pp$ argument.","marker":"[15]"},{"why":"Defines the twin-prime companion (TPC) object used to build the $pp$ cross-section curve.","marker":"[16, 17]"}],"fun_headline_variants":["Prime-index parametrization fits neutrino cross sections from MeV to PeV","Twin primes trace the (ln s)^2 rise in pp cross sections","One constant plus primes predicts neutrino cross sections over six decades","Primes as index: neutrino cross sections from MeV to PeV with single factor","Prime-number table reproduces neutrino and pp cross sections"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the fitted normalization constants (0.70, 0.26, 0.48, 0.69, 0.25) remain valid outside the energy ranges where they were set, so the high-energy extrapolations and the IceCube event count inherit them; the $pp$ side also leans on the twin-prime asymptotic, which is a conjecture rather than a theorem.","fun_headline_variants_meta":{"raw":{"variants":["Prime-index parametrization fits neutrino cross sections from MeV to PeV","Twin primes trace the (ln s)^2 rise in pp cross sections","One constant plus primes predicts neutrino cross sections over six decades","Primes as index: neutrino cross sections from MeV to PeV with single factor","Prime-number table reproduces neutrino and pp cross sections"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000704,"raw_usage":{"total_tokens":3157,"prompt_tokens":912,"completion_tokens":2245,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":528,"completion_tokens_details":{"reasoning_tokens":2153}},"tokens_in":528,"tokens_out":2245,"duration_ms":14650,"temperature":1.0,"reasoning_tokens":2153,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:59:46.401700+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure or extract the neutrino-nucleon charged-current cross section at a single energy near 1 PeV; the method expects roughly ten times the Standard Model value at that energy, so a measurement near the Standard Model prediction would refute the high-energy extrapolation. In the same spirit, an updated IceCube sample with several times the current exposure that finds events above 0.5 PeV consistent with the Standard Model estimate rather than the $11\\pm3.3$ prediction would settle the claim.","supporting_citations":[{"cited_title":"The rate for these muon ne utrinos have been compared with the SM calculations of Cooper-Sarkar, et al., [ 11]","cited_arxiv_id":null,"evidence_quote":"Provides the low-energy (MeV) neutrino cross-section approximation that the prime-index prescription must match in the lowest decade."},{"cited_title":"69 for the data in figures 2 and 3, respectively","cited_arxiv_id":null,"evidence_quote":"Supplies a second low-energy analytic approximation used as a comparison anchor in the few-hundred-MeV region."},{"cited_title":"(9) The suggestion, originally by Heisenberg[ 13], proposed a universal logarithmic increase in the pp cross sections of the form (ln s)2","cited_arxiv_id":null,"evidence_quote":"Gives the measured linear energy dependence in the GeV-TeV range that fixes the overall slope the prime-index method reproduces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"States a standard high-energy approximation that the prime-index extrapolation is compared against at ultra-high energies."},{"cited_title":"This represents the neutrino-nucleon cross section in ab","cited_arxiv_id":null,"evidence_quote":"Provides the digitized experimental neutrino and antineutrino cross-section data plotted in figure 1, the main data anchor for the method."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the IceCube astrophysical muon-neutrino rates used to test the prime-index prediction at PeV energies."},{"cited_title":"Neutrino nucleus scattering","cited_arxiv_id":"nucl-th/9901027","evidence_quote":"Gives the Standard Model cross-section calculations used in table 2 as the baseline the prime-index values are compared with at high energy."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the $(\\ln s)^2$ growth of $pp$ cross sections that the twin-prime companion ratio is claimed to explain."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the twin-prime asymptotic $\\pi_2(x)\\sim 2C_2\\,x/(\\ln x)^2$ that produces the logarithmic-squared factor in the $pp$ argument."}],"review_version":1}