REVIEW 2 major objections 6 minor 30 references
High-energy neutrino event simulation at NLO in Genie for KM3NeT and other observatories
T0 review · 2 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read High-energy neutrino events can now be simulated at NLO to EeV energies.
desk verdict Useful NLO event-generator bridge with solid rate validation; the heavy-quark final-state claim outruns the evidence. 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 load-bearing machinery is a lookup table of next-to-leading-order structure functions. For each interaction channel HEDIS precomputes $F_1$, $F_2$, and $F_3$ with a numerical DGLAP-evolution code at NLO for a given PDF set and stores them as functions of $\log_{10}x$ and $\log_{10}Q^2$, turning the expensive NLO calculation into fast interpolation during event generation. A stored maximal differential cross section on the same grid feeds an acceptance–rejection sampler over $\log_{10}x$ and $\log_{10}y$. After the kinematics are fixed, the flavour of the outgoing quark is picked using LO structure functions via Eq. (2.3), with slow rescaling suppressing heavy-quark production near threshold, and a hadronization model turns the resulting quark system into final-state particles. This separation is what lets an NLO cross section drive the rate and kinematics while retaining a simple LO rule for flavour.
What would settle it
Generate HEDIS event samples at $E_\nu = 10^5$ and $10^7$ GeV and compare the charm-, bottom-, and top-production fractions, or the rate of events with leptonic heavy-flavour decays, against an event-level NLO calculation that determines the outgoing flavour from NLO matrix elements rather than from LO structure functions using the same PDFs. A difference larger than the PDF uncertainties would falsify the paper's flavour-assignment approximation.
Extended reading notes
Core claim
On the paper's own terms, the central claim is that NLO event generation for neutrino deep inelastic scattering can be built by separating the cross-section computation from the flavour assignment. HEDIS computes the NLO double-differential cross section from structure functions obtained with numerical DGLAP evolution for a chosen PDF set, stores them as functions of $\log_{10}x$ and $\log_{10}Q^2$, and samples $x$ and $y$ by acceptance–rejection over the full phase space. The outgoing quark flavour is then drawn using LO structure functions with slow rescaling, so charm, bottom, and top production are included with kinematic suppression near threshold. The paper validates the approach by showing that the total charged- and neutral-current cross sections on an isoscalar target reproduce two published NLO calculations and the LO calculation used in the telescope's letter of intent within about 1%, and it presents event-level kinematics from neutrino–oxygen interactions at $10^3$ and $10^7$ GeV.
Load-bearing premise
The load-bearing premise is that assigning the struck and outgoing quark flavour with leading-order structure functions, after the NLO kinematics are chosen, does not meaningfully change the simulated heavy-quark content or the secondary leptons from heavy-flavour decays; if NLO flavour composition differs sharply from LO, especially through gluon-initiated processes at low $x$, those final-state yields would be biased.
Editorial extensions
If this is right
- The generator's usable energy range extends from roughly 1 TeV to $10^{10}$ GeV, so neutrino-telescope analyses can simulate events across the full range of detected and expected fluxes.
- Total charged- and neutral-current cross sections match two independent published NLO calculations and the LO calculation used in the deep-sea telescope's letter of intent at the 1% level over the tested range.
- NLO PDFs are used consistently with NLO cross sections, removing the known inconsistency in generators that pair LO cross sections with NLO parton densities.
- The flavour-assignment scheme includes charm, bottom, and top in the final state, so secondary leptons from heavy-flavour decays appear in simulated event samples.
- Precomputed structure-function tables make NLO sampling fast enough for routine event-by-event generation.
Reading between the lines
- The paper validates total cross sections, not differential event-level distributions; the most uncertain region is very low $x$, where PDFs are extrapolated beyond their fitted range, and this is where EeV cosmogenic-neutrino rates live.
- The same precomputed-table machinery could incorporate nuclear PDFs, using the procedure the paper cites, to add nuclear shadowing corrections that the current free-nucleon treatment omits.
- A direct comparison with an NLO-matched parton shower would isolate how much the LO flavour-assignment approximation biases heavy-flavoured final states; this is the most exposed part of the generator.
- Charm production is abundant in the EeV range, so using HEDIS to model the prompt-muon background from charm decays could sharpen astrophysical-neutrino searches in current and future telescopes.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper reports HEDIS, an extension of the GENIE neutrino event generator that provides deep-inelastic-scattering (DIS) event simulation up to EeV energies for neutrino telescopes such as KM3NeT. HEDIS computes NLO structure functions with APFEL (cross-checked with QCDNUM), uses LHAPDF grids with NLO or NNLO PDFs, samples event kinematics from the NLO double-differential cross section, and interfaces PYTHIA6 for hadronization. The total CC and NC cross sections are compared with the published CSMS, BRG, and KM3NeT LoI calculations; the authors report agreement at the 1% level. For the hadronic final state, the flavor of the outgoing quark is assigned using LO structure functions after the NLO kinematics are generated, with a slow-rescaling prescription for heavy quarks. The paper claims that this scheme correctly simulates heavy-quark production, including bottom and top.
Significance. If fully validated, HEDIS would fill a genuine need: existing neutrino-generator DIS implementations are leading-order and valid only to about 1 TeV, while modern PDFs are extracted at NLO/NNLO. The inclusive cross-section engine is benchmarked against two independent published calculations (CSMS and BRG) using independent DGLAP evolution tools (APFEL and QCDNUM) and LHAPDF, which is a solid basis for the total-rate calculation. The paper is also useful as a description of a tool that will be available in GENIE v4. However, the final-state claim of the abstract — that heavy-quark production, including bottom and top, is 'correctly simulated' — is not established by the evidence presented, because the only quantitative validation concerns inclusive cross sections, which are insensitive to the flavor composition of the final state.
major comments (2)
- [Section 2.2, Eq. (2.3)] The abstract's claim that heavy-quark production, including bottom and top, is 'correctly simulated' is not supported by the evidence in the paper. After x and Q2 are sampled from the NLO cross section, the outgoing quark flavor is assigned by Eq. (2.3) using leading-order structure functions. This LO probability contains no explicit NLO gluon-initiated subprocesses (e.g., gluon splitting to ccbar followed by charm scattering), so such events are forced into a LO quark-to-quark assignment and the simulated bottom and top yields may be biased. The 1% agreement in Fig. 2 validates only the inclusive total cross section, which is insensitive to the flavor composition of the final state. Since Section 2.2 itself notes that heavy-flavor decays can produce additional leptons, the paper needs either an event-level validation of heavy-quark rates or a restriction of the claim to the inclusive cross section.
- [Section 2.2, slow rescaling] The slow-rescaling prescription used for the LO flavor pick may not reproduce the threshold behavior of the NLO massive scheme used for the cross-section calculation, particularly for top production. Figure 4 shows that the top content in HEDIS-BGR differs strongly from HEDIS-CSMS, but no comparison is made to any validated NLO event-level prediction. Without such a comparison, the statement that top production is correctly simulated remains an assumption rather than a demonstrated property of the generator.
minor comments (6)
- [Throughout] The labels 'BRG' and 'BGR' are used inconsistently (e.g., 'BGR' in Fig. 2 but 'BRG' in the text); please standardize the abbreviation for the Bertone-Gauld-Rojo calculation.
- [Eq. (2.3)] The notation in Eq. (2.3) is not fully defined: the summation index j, the CKM factor, and the distinction between sigma_LO_f and sigma_LO_tot should be written explicitly so that the probability is unambiguous.
- [Fig. 4 caption] The caption contains a typo, 'muon-netrino'; it should read 'muon-neutrino'.
- [Section 2.3] The text says that x and y are sampled, but the actual sampling is over log10 x and log10 y; this should be stated consistently.
- [Section 2.1, Fig. 2] The claim of agreement at the 1% level would benefit from a precise definition of the quoted difference and from a discussion of PDF uncertainties, since the authors identify PDFs as the dominant uncertainty source.
- [References] The arXiv identifier attached to Ref. [14] appears to be inconsistent with the journal citation; please verify it, and similarly check the completeness of the other references.
Circularity Check
No significant circularity: the NLO cross-section benchmark is checked against independent published calculations and the LO flavour-picking is an acknowledged approximation.
full rationale
The central cross-section result is self-contained: HEDIS computes NLO structure functions with APFEL/QCDNUM and LHAPDF PDFs, and the total cross sections are compared with CSMS, BRG, and KM3NeT-LOI, which are independent published calculations using different codes and mass schemes. The 1% agreement is an implementation benchmark, not a construction: HEDIS-CSMS reproduces CSMS by choosing the same PDFs and mass scheme, including switching off the top mass above the 5-flavour threshold, but that choice is disclosed and is a matching condition, not a fitted parameter. The generator is also validated against BRG (FONLL) and the LOI, so no single input is predicted from itself. The flavour-assignment scheme in Eq. (2.3) explicitly uses LO structure functions as a "first pragmatic step"; this is an acknowledged approximation, and while it leaves the heavy-quark-final-state claim not fully validated, it is not circular because HEDIS does not derive the NLO flavour composition from that equation. The only self-citation, KM3NeT LoI, is one of three benchmarks and is not load-bearing; CSMS and BRG are external, so no self-citation chain supports the central claim. No fitted parameter is renamed as a prediction.
Assumptions & free parameters
free parameters (1)
- Top-quark mass treatment in ZM-VFNS scheme =
On/off switch: top mass off above 5-flavour threshold
assumptions (6)
- standard math NLO QCD factorization and DGLAP evolution describe neutrino-nucleon DIS structure functions.
- domain assumption PDFs are valid outside their fitted x and Q2 ranges through extrapolation.
- ad hoc to paper Outgoing quark flavor can be sampled from LO structure functions at NLO-sampled kinematics.
- domain assumption PYTHIA6 hadronization is adequate for the hadronic final state up to 10^10 GeV.
- domain assumption Nuclear effects (shadowing) are negligible for the simulated targets.
- domain assumption The chosen Q2 lower limits (Table 1) correctly separate the perturbative and non-perturbative regimes.
Cite this review
Pith. "Pith review of High-energy neutrino event simulation at NLO in Genie for KM3NeT and other observatories." pith.science (2026). https://pith.science/paper/SK6EOPCM
@misc{pith2026190810077,
author = {Pith},
title = {Pith review of: High-energy neutrino event simulation at NLO in Genie for KM3NeT and other observatories},
year = {2026},
howpublished = {\url{https://pith.science/paper/SK6EOPCM}},
note = {Machine review of arXiv:1908.10077}
}
abstract
High-energy neutrino astronomy is a key pillar of multi-messenger astronomy and has the potential to also advance fundamental neutrino physics. An accurate simulation of the neutrino interactions is key in any analysis of neutrino telescope data. The currently available generator codes of neutrino interactions use leading order expressions for the differential cross sections which govern rate and kinematics of the events. These calculations are affected by well-known large theoretical uncertainties, hence the need of beyond leading order formalism. Also, a consistent use of modern Parton Density Functions, which are derived using NLO (or NNLO) frameworks, is required. For this reason and others, the GENIE event generator of neutrino interactions, which is a standard tool in the few-GeV energy regime, has so far been stated to be valid up to only about 1 TeV. The work reported here consists of a high-energy (up to 10$^{10}$ GeV) extension, which will be available in GENIEv4. It interfaces PYTHIA6 for the hadronization and employs a pragmatic method to assign the flavour of the struck and outgoing quarks, so that the heavy quark production, including bottom and top, is also correctly simulated.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
IceCube Collaboration, Phys. Rev. Lett. 113 (2014) 101101 [1405.5303]
arXiv 2014
-
[2]
ANTARES Collaboration, Astrophys. J. Lett 853 (2018) 1 [1711.07212]
arXiv 2018
-
[3]
IceCube Collaboration, Science 361 (2018) 147 [1807.08794]
arXiv 2018
-
[4]
ANTARES Collaboration, JHEP 06 (2019) 113 [1812.08650]
work page Pith review arXiv 2019
-
[5]
IceCube Collaboration, Nature 551 (2017) 596 [1711.08119]
arXiv 2017
-
[6]
KM3NeT Collaboration, J. Phys. G43 (2016) 084001 [1601.07459]
arXiv 2016
-
[7]
IceCube-Gen2 Collaboration, PoS FRAPWS2016 (2017) 004 [1412.5106]
arXiv 2017
-
[8]
KM3NeT/ANTARES Collaborations, Workshop on MC simulation for neutrino telescopes (2018)
work page 2018
Show all 30 references
-
[9]
Ingelman et al., Comput
G. Ingelman et al., Comput. Phys. Commun. 101 (1997) 108 [hep-ph/9605286]
1997 arXiv
-
[11]
Cooper-Sarkar et al., JHEP 08 (2011) 042 [1106.3723]
A. Cooper-Sarkar et al., JHEP 08 (2011) 042 [1106.3723]
2011 arXiv
-
[12]
Bertone et al., JHEP 01 (2018) 217 [1808.02034]
V . Bertone et al., JHEP 01 (2018) 217 [1808.02034]
2018 arXiv
-
[13]
Sjostrand et al., JHEP 05 (2006) 2006 [hep-ph/0603175]
T. Sjostrand et al., JHEP 05 (2006) 2006 [hep-ph/0603175]
2006 arXiv
-
[14]
Andreopoulos et al., Nucl
C. Andreopoulos et al., Nucl. Instrum. Meth. A614 (2010) 1 [astro-ph/0406439]
2010 arXiv
-
[15]
Altarelli et al., Nucl
G. Altarelli et al., Nucl. Phys. B126 (1977) 298
1977
- [16]
-
[17]
V ogt et al., Nucl
A. V ogt et al., Nucl. Phys. Proc. Suppl. 160 (2006) 44 [hep-ph/0608307]
2006 arXiv
-
[18]
Bertone et al., Comput
V . Bertone et al., Comput. Phys. Commun. 185 (2014) 1647 [1310.1394]
2014 arXiv
- [19]
- [20]
-
[21]
Collins et al., Nucl
J. Collins et al., Nucl. Phys. B278 (1986) 934
1986
-
[22]
R. D. Ball et al., Phys. Lett. B754 (2016) 49 [1510.00009]
2016 arXiv
- [23]
- [24]
-
[25]
Pumplin et al., JHEP 07 (2002) 2002 [hep-ph/0201195]
J. Pumplin et al., JHEP 07 (2002) 2002 [hep-ph/0201195]
2002 arXiv
- [26]
- [27]
- [28]
-
[29]
R. M. Barnett, Phys. Rev. Lett. 36 (1976) 1163
1976
-
[30]
Sjostrand et al., Comput
T. Sjostrand et al., Comput. Phys. Commun. 178 (2008) 852 [0710.3820]
2008 arXiv
-
[31]
Gauld, Nikhef 2019-011 (2019) [1905.03792]
R. Gauld, Nikhef 2019-011 (2019) [1905.03792]. 7
2019 arXiv
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.