{"id":"8811697b-eabb-40b8-b050-de5773665790","arxiv_id":"1908.10077","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"HEDIS, a GENIE extension, implements NLO neutrino-nucleon cross sections with PYTHIA6 hadronization, extending neutrino-telescope event simulation to EeV energies.","lead":"The KM3NeT collaboration presents HEDIS, a new extension of the GENIE neutrino event generator that simulates high-energy neutrino interactions up to 10^10 GeV using next-to-leading-order cross sections. It could give neutrino telescopes like KM3NeT and IceCube a more reliable way to simulate the rare, energetic events they are built to detect.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The heavy-quark final-state claim rests on an unvalidated LO flavor-assignment approximation: NLO gluon-initiated channels are not represented in Eq. 2.3, so simulated bottom/top yields may be biased.","rationale":"The reader's weakest-assumption identification is exactly the load-bearing concern: the flavor assignment of Eq. 2.3 uses LO structure functions after NLO kinematics are sampled. This is not merely a cosmetic approximation; it changes the physics of the final state. At NLO the concept of a single 'outgoing quark' is not well-defined, and gluon-initiated contributions are an important part of the cross section at low x, the region relevant at EeV energies. Because the paper's abstract explicitly claims that heavy-quark production, including bottom and top, is 'correctly simulated,' and because Section 2.2 identifies leptonic decays of heavy flavor as a notable observable, the unvalidated LO flavor assignment is directly load-bearing for the central claim. The 1% total-cross-section agreement with CSMS and BRG is genuine independent support for the inclusive rate calculation, but it does not constrain flavor composition. The paper itself flags the approximation as pragmatic and offers no event-level validation; this is an internal admission that the needed support is missing, not an artifact of the reading. A concrete heavy-flavor yield comparison against the same NLO formalism would settle the issue. Since the reader's CONDITIONAL verdict already captures this gap, the recommendation is UNCHANGED: the concern reinforces the conditional status rather than moving it to accept or reject.","tokens_in":7527,"tokens_out":3207,"duration_ms":38025,"concrete_test":"Generate 10^6 CC neutrino-Oxygen events with HEDIS-BGR at E_nu = 10^7 GeV and count the number of primary b- and t-quark events (or the rate of leptons from heavy-flavor decays). Compare these event-level yields with the heavy-quark production cross sections predicted by the same NLO formalism, e.g., the b- and t-quark component of the FONLL structure functions used in the BRG calculation. If the generated heavy-quark yield deviates by more than ~10% from the NLO prediction, the LO flavor-assignment scheme biases the central claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim that heavy-quark production, including bottom and top, is 'correctly simulated' depends on Section 2.2's flavor-assignment scheme. After x and Q2 are sampled from the NLO double-differential cross section, Eq. 2.3 assigns the outgoing quark flavor using LO structure functions only. At NLO, the inclusive cross section receives substantial contributions from gluon-initiated subprocesses (e.g., g -> c cbar followed by c scattering), where there is no single 'struck quark'; these channels are entirely absent from the LO-based probability in Eq. 2.3. The scheme forces every NLO event into a LO q -> q' assignment, so heavy quarks produced through NLO gluon splitting are mis-modeled as arising from sea quarks. In addition, the LO 'slow rescaling' used for the flavor pick may not reproduce the NLO mass-scheme threshold behavior for top production, especially near threshold. The paper presents no event-level validation of heavy-quark rates: the 1% agreement in Fig. 2 validates only the inclusive total cross section, which is insensitive to the flavor composition of the final state. Section 2.2 itself acknowledges the method is a 'first pragmatic step' and notes that heavy-flavor decays can produce additional leptons, making this directly relevant to neutrino-telescope observables. Thus the claim that heavy-quark production is correctly simulated is unsupported by the evidence shown, and the missing validation is a load-bearing gap.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":7811,"tokens_out":4905,"duration_ms":51599,"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":[{"comment":"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":"Section 2.2, Eq. (2.3)"},{"comment":"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.","section":"Section 2.2, slow rescaling"}],"minor_comments":[{"comment":"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.","section":"Throughout"},{"comment":"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.","section":"Eq. (2.3)"},{"comment":"The caption contains a typo, 'muon-netrino'; it should read 'muon-neutrino'.","section":"Fig. 4 caption"},{"comment":"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":"Section 2.3"},{"comment":"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.","section":"Section 2.1, Fig. 2"},{"comment":"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.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The core cross-section benchmarking is sound and useful, but the manuscript's strongest advertised feature — NLO final states with correct heavy-quark content — is presently unsupported. I would be willing to accept after the authors either provide event-level validation of bottom and top production rates against a trusted NLO framework, or explicitly soften the abstract and Section 2.2 to state that only the kinematics and inclusive rates are NLO while flavor assignment remains a leading-order approximation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is the first GENIE extension that samples NLO double-differential DIS cross sections for neutrino telescopes, and the total-rate validation against CSMS and BRG is real. The heavy-quark final-state claim in the abstract is stronger than what Section 2.2 actually delivers.\n\nWhat is new: NLO structure functions from APFEL, cross-checked with QCDNUM; NLO PDFs through LHAPDF; PYTHIA6 hadronization inside GENIE; three configurations (CSMS, BGR, LOI) with total cross sections matching published calculations at the 1% level. That is a concrete step forward for KM3NeT and IceCube-Gen2, which currently rely on LO generators. The LO flavor-picking scheme (Eq. 2.3) is described openly, including slow rescaling; the paper does not hide that it is a first pragmatic step.\n\nThe soft spots are real but localized. The abstract's \"correctly simulated\" heavy-quark production is not supported by the evidence shown. After x and Q2 are sampled from the NLO cross section, flavor is assigned using LO structure functions only; gluon-initiated NLO channels have no single struck quark and are absent from Eq. 2.3. The 1% agreement validates inclusive total cross sections, which are insensitive to flavor composition. There is no event-level check of bottom or top yields, or of lepton rates from heavy-flavor decays. That gap is load-bearing for the abstract's claim, although Section 2.2 partially anticipates it with its \"first pragmatic step\" caveat. Secondary issues: no code link or commit hash, no uncertainty bands on the validation plots, and the choice to switch off the top mass to match CSMS is a tuning choice rather than a prediction. These are minor for a proceedings paper but would matter in a full technical write-up.\n\nBottom line: the central rate calculation is credible and useful. The event-level heavy-quark simulation is an announced approximation, not a validated result. A serious referee should ask for either softening the abstract or adding a validation of heavy-flavor yields, e.g. against a NLO flavor decomposition or a tuned Monte Carlo. I would bring this to a reading group concerned with neutrino-telescope simulation and would cite it as the GENIE/NLO bridge for high-energy event generation. It deserves peer review; it is a tool paper with one clear overclaim that revision can fix.","headline":"Useful NLO event-generator bridge with solid rate validation; the heavy-quark final-state claim outruns the evidence.","tokens_in":8396,"tokens_out":2700,"would_cite":true,"duration_ms":26645,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"High-energy neutrino events can now be simulated at NLO to EeV energies.","keywords":["neutrino deep inelastic scattering","NLO cross sections","event generator","heavy quark production","neutrino telescope simulation","parton density functions","structure functions","EeV neutrinos"],"falsifier":"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.","tokens_in":7280,"feed_emoji":"🔭","tokens_out":11289,"duration_ms":101295,"temperature":0.7,"pith_summary":"The paper presents a package that lets the standard few-GeV neutrino event generator produce deep-inelastic neutrino–nucleus scattering events at neutrino energies up to $10^{10}$ GeV. It replaces the generator's leading-order cross-section calculation with next-to-leading-order structure functions, computed from modern parton density functions and stored in precomputed lookup tables, and it supplies a pragmatic scheme for assigning the flavour of the struck and outgoing quark. The result is a generator whose total charged- and neutral-current cross sections reproduce two published NLO calculations and the LO calculation used in a deep-sea telescope's letter of intent at the 1% level, while also producing heavy-flavour final states including bottom and top quarks. The paper claims this removes the large theoretical uncertainties attached to LO generators and extends reliable simulation into the energy range relevant for high-energy neutrino observatories.","feed_headline":"Neutrino-event simulation extends to EeV energies with NLO accuracy","feed_subtitle":"A generator extension reproduces published NLO cross sections at 1 percent and adds heavy-quark final states.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Published NLO calculation used as a benchmark and reproduced by the HEDIS configuration.","marker":"[11]"},{"why":"Second published NLO calculation used as a benchmark; its mass scheme treats heavy-quark production differently.","marker":"[12]"},{"why":"Hadronization model that turns the struck-quark system into final-state particles.","marker":"[13]"},{"why":"The host generator to which HEDIS is added; it previously had a stated validity limit near 1 TeV.","marker":"[14]"},{"why":"PDF interface used to read and interpolate the parton densities into the lookup tables.","marker":"[16]"},{"why":"Numerical DGLAP-evolution code used to compute the NLO structure functions.","marker":"[18]"},{"why":"Alternative numerical DGLAP-evolution code tested to check the structure-function results.","marker":"[19]"},{"why":"PDF set used in the benchmark configuration that reproduces the first NLO calculation.","marker":"[23]"},{"why":"PDF set used in the benchmark configuration that reproduces the second NLO calculation.","marker":"[24]"},{"why":"Slow-rescaling prescription that suppresses heavy-quark production in the flavour-assignment step.","marker":"[29]"}],"fun_headline_variants":["Neutrino simulation hits NLO up to 10^10 GeV","GENIE gets NLO accuracy for high-energy neutrinos","Heavy-quark final states now in NLO neutrino events","NLO neutrino cross sections validated to 1% at EeV","EeV neutrino simulation with NLO and top quarks"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Neutrino simulation hits NLO up to 10^10 GeV","GENIE gets NLO accuracy for high-energy neutrinos","Heavy-quark final states now in NLO neutrino events","NLO neutrino cross sections validated to 1% at EeV","EeV neutrino simulation with NLO and top quarks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000247,"raw_usage":{"total_tokens":1554,"prompt_tokens":967,"completion_tokens":587,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":500}},"tokens_in":583,"tokens_out":587,"duration_ms":6327,"temperature":1.0,"reasoning_tokens":500,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:53:22.421383+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"ANIS: High Energy Neutrino Generator for Neutrino Telescopes","cited_arxiv_id":"astro-ph/0406439","evidence_quote":"The host generator to which HEDIS is added; it previously had a stated validity limit near 1 TeV."},{"cited_title":"Proton Structure from HERA to LHC","cited_arxiv_id":"1012.1438","evidence_quote":"PDF set used in the benchmark configuration that reproduces the first NLO calculation."},{"cited_title":"Precision determination of the small-$x$ gluon from charm production at LHCb","cited_arxiv_id":"1610.09373","evidence_quote":"PDF set used in the benchmark configuration that reproduces the second NLO calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Slow-rescaling prescription that suppresses heavy-quark production in the flavour-assignment step."}],"review_version":1}