{"id":"77c4dd23-d1c1-4908-aba3-129fb9188d49","arxiv_id":"2608.09608","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A leading-order MadGraph study maps the expected heavy-quark and double-charm yields and nuclear modification ratios for early EIC beam configurations.","lead":"This paper computes leading-order predictions for charm and beauty production in electron-proton and electron-nucleus collisions at the future Electron-Ion Collider, including expected event yields for early running. It is a feasibility baseline that shows which kinematic regions the early EIC programme could access for heavy-quark and double-charm measurements.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's own arithmetic contradicts the like-sign D0D0 measurability claim: ~5 expected events versus the stated 20-event threshold, undermining the reach methodology that supports the headline conclusions.","rationale":"The paper provides a useful leading-order parton-level baseline with transparent kinematic cuts, scale and PDF uncertainties, and a reproducible MadGraph setup. The central claim has two pillars: the cross-section baseline and the statistical-reach/observability statements. The same-sign D0D0 conclusion is not a peripheral remark; the abstract and Section 5 present it as a new capability beyond HERA. Yet the paper's own Eq. (11) and Fig. 6 yield ~5 events, a factor of ~4 below the 20-event threshold that every observability line in the paper is based on. This is not a matter of external consensus or detector assumptions; it is internal arithmetic. If the threshold is taken seriously, the same-sign claim fails; if the threshold is not taken seriously, the reach conclusions are not quantitatively supported. Either way, the published claims need correction. The reader's verdict of CONDITIONAL is appropriate because the flaw is localized and correctable; the inclusive charm/beauty baseline and R_eA trends are not invalidated by this issue.","tokens_in":18555,"tokens_out":27230,"duration_ms":257179,"concrete_test":"Recompute the same-sign D0D0 yield from the published inputs: take σ = 1.02 pb (central value, Fig. 6), L = 5.85 fb^-1, ε = 0.03, and N = σ · L · ε^2 with unit conversion 1 pb·fb^-1 = 1000 events. If the result is below 20 (it is ~5.4), then either the observability threshold must be lowered or the like-sign measurability claim should be retracted; also scan the other yield estimates (e.g., eA reach lines) to confirm no other claim falls below the same threshold.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 5 and the Abstract claim that like-sign D0D0 pairs are 'potentially measurable as early as the first year of operation,' a capability highlighted as beyond HERA. However, Section 4.1 reports the same-sign D0D0 + D0barD0bar cross section as σ = 1.02^{+1.12}_{-0.46} pb (Fig. 6), and the paper's own yield formula N = σ·L·ε^2 with L = 5.85 fb^-1 and ε = 3% gives N ≈ 1.02 pb × 5.85 fb^-1 × 1000 × (0.03)^2 ≈ 5.4 events. This is below the observability threshold N_ev = 20 defined in Eq. (11) and used to draw every observability line in Section 4. The claim is therefore not supported by the stated assumptions. Because the same 20-event criterion is the quantitative basis for the 'substantial improvement in statistical reach' conclusions, the inconsistency also puts pressure on the central claim, not just on this side channel.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents leading-order predictions, made with MadGraph5 aMC@NLO, for inclusive charm and beauty pair electroproduction in ep and eA collisions at representative early-EIC energies, together with estimates of double-charm D0-pair production in ep. Differential distributions in Q^2, heavy-quark pT and rapidity are shown with scale and PDF uncertainties, and event-yield estimates are derived from projected luminosities, assumed fragmentation fractions and detection efficiencies. For eA collisions, the nuclear modification ratio R_eA is computed with EPPS21 nPDFs. The stated goal is to provide a parton-level baseline for early EIC heavy-flavour measurements and to demonstrate an improved statistical reach relative to HERA, including an explicit claim that like-sign D0D0 pairs may be measurable in the first year of EIC operation.","tokens_in":18831,"tokens_out":8059,"duration_ms":74658,"significance":"If taken as a baseline projection rather than a precision prediction, the study is useful for early EIC feasibility discussions. Its strengths are that the calculations use a standard public code with clearly stated PDFs, masses, scales, acceptance cuts, and explicit scale and PDF uncertainty bands, and that the setup is reproducible via the NLOAccess online platform. The paper does not fit any parameters to data, so there is no internal circularity in the cross-section calculation. The main limitations are the LO accuracy of the DIS matrix elements, the use of efficiencies inherited from ECCE and CMS rather than ePIC, and the fact that the observed R_eA pattern (shadowing, anti-shadowing, EMC-like suppression) is inherited from the EPPS21 input nPDFs rather than being an independent finding. These limitations reduce the strength of the paper's 'demonstrate substantial improvement' language but do not invalidate the baseline purpose.","major_comments":[{"comment":"The claim that like-sign D0D0 pairs are 'potentially measurable as early as the first year of operation' is contradicted by the paper's own arithmetic. The same-sign cross section is stated as σ = 1.02^{+1.12}_{-0.46} pb (Fig. 6), and the yield formula N = σ_DD · L · ε^2 with L = 5.85 fb^-1 and ε = 3% gives N ≈ 1.02 × 5.85 × 1000 × (0.03)^2 ≈ 5.4 events, i.e. below the N_ev = 20 observability threshold defined in Eq. (11). Even the upper scale-variation value, 2.14 pb, gives only about 11 events, still below the paper's own threshold. If instead σ_DD is intended not to include the fragmentation fractions, the yield is even smaller. The same issue affects the time line: according to Table 1, the ep 10×130 GeV configuration with L = 5.85 fb^-1 is a Year 2 scenario, not Year 1. The Section 5 statement, and the corresponding 'beyond HERA' capability claim, should be removed or replaced by a statement consistent with the stated threshold and luminosity schedule.","section":"Section 4.1, Fig. 6; Section 5"},{"comment":"The entire quantitative reach argument rests on the ad hoc choice N_ev = 20, which is justified only by the statement that it 'roughly corresponds to the statistical accuracy of the typical last bins of HERA measurements'. This is not substantiated with a HERA comparison, and the same-sign D0D0 example shows that the criterion is not applied consistently to all claimed measurability statements. The authors should either (i) derive N_ev from a specified target statistical uncertainty or significance, including a background model, or (ii) clearly label all observability lines as illustrative and avoid claiming that channels below the line are measurable. As written, the 'substantial improvement in statistical reach relative to HERA' conclusion in Section 5 is not quantitatively supported by the presented criterion.","section":"Section 3, Eq. (11); all observability lines in Section 4"},{"comment":"The paper states in several places that R_eA 'exhibits' shadowing, anti-shadowing and EMC-like suppression. Since R_eA is computed by dividing cross sections obtained with EPPS21 and CT18ANLO, this pattern is necessarily inherited from the EPPS21 nuclear gluon distribution; it is a prediction of the input nPDFs, not an observable finding. This is acceptable for a baseline study, but the wording should be changed (e.g. 'predicts', 'is consistent with EPPS21 expectations') so that readers are not led to think the calculation itself provides evidence for these nuclear effects. The current wording in Section 5 ('the nuclear modification factor R_eA exhibits a rich structure') overstates the novelty of the result.","section":"Section 4.2, Figs. 9 and 10"}],"minor_comments":[{"comment":"The phrase 'demonstrate the substantial improvement in statistical precision over HERA' is too strong for a LO parton-level study with assumed efficiencies from other detectors; 'suggest' or 'project' would be more appropriate.","section":"Abstract and Section 5"},{"comment":"The labels on the observability lines (e.g. '20 ev / 25 GeV^2 / 5.85 fb^-1') do not display the fragmentation fraction and efficiency values used for each flavour; including them would make the lines reproducible from the stated inputs.","section":"Section 3, Eq. (11)"},{"comment":"The dagger footnote states that the beauty upper range is 120–175 GeV^2, but the table row is labelled '120–200'; the bin label and the footnote should be made consistent.","section":"Table 2"},{"comment":"The legend uses σcc and σcbar cbar without kinematical subscript bars; using σ(c c) and σ(cbar cbar) with the same rapidity cuts as in the text would improve readability.","section":"Fig. 6"},{"comment":"The text says 'charm cross sections (in nb) are higher than those for beauty (in pb)' across the entire kinematic domain; this is true, but the statement would be clearer if the units were not compared directly without the numerical factors.","section":"Section 4.1, Table 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is a reasonable baseline projection and should be published after the measurability claim for same-sign D0D0 pairs is corrected and the observability criterion is made self-consistent. The contradiction between Fig. 6 and Section 5 is a load-bearing error because it directly supports one of the headline conclusions; it can be fixed within the scope of the manuscript by revising the claim and the threshold presentation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Useful paper for anyone planning EIC heavy-flavour measurements. It provides a clean, reproducible LO collinear-factorisation baseline: MadGraph with CT18ANLO and EPPS21, stated masses, scales, ePIC-like cuts, and scale/PDF uncertainties on every distribution. The new content is the three nuclear targets (Cu, Ag, Au), the early-science luminosities, and the double-charm-pair channel. I'd trust the central cross sections and the R_eA trends as a planning reference.\n\nThe soft spot is the same-sign D0D0 claim. The paper's own numbers give about 5 events for L=5.85 fb^-1 and ε=3%, below its N_ev=20 observability threshold, so the abstract and conclusion statement that like-sign pairs are 'potentially measurable as early as the first year' is not supported. Also, 5.85 fb^-1 is the Year 2-3 ep luminosity, not Year 1. The reach methodology itself works fine for single charm and beauty; the problem is specific to this side channel and needs a correction or removal.\n\nOther issues are minor. The comparison to HERA precision is asserted rather than quantified. The conclusion says charm cross sections are 'several hundreds of pb' while Table 2 lists them in nb. And the use of NLO PDF sets at LO is a scheme choice that deserves a comment. None of these affect the central baseline.\n\nThe parton-level baseline and the R_eA predictions for Cu/Ag/Au hold up. Send it to peer review, but ask for a revision that fixes the like-sign D0D0 event-yield claim and cleans up the prose errors. This is a citable reference for EIC heavy-flavor studies.","headline":"Useful LO baseline for EIC heavy-flavor planning, but the like-sign D0D0 'first-year' claim is contradicted by the paper's own arithmetic and should be corrected.","tokens_in":19332,"tokens_out":6386,"would_cite":true,"duration_ms":49408,"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":"Early EIC runs will map charm and beauty beyond HERA","keywords":["heavy-quark electroproduction","deep inelastic scattering","Electron-Ion Collider","nuclear parton distribution functions","gluon shadowing","EMC effect","double-charm production","leading-order QCD"],"falsifier":"A concrete check is to apply the paper's own numbers to the same-sign D0D0 channel: with $\\sigma \\approx 1.02$ pb, $L = 5.85$ fb$^{-1}$, and $\\epsilon = 3\\%$, the yield is about 5 events, below the 20-event threshold used elsewhere, so reaching 20 events would require $\\epsilon \\approx 5.8\\%$; early EIC runs that record fewer than 20 same-sign pairs, or a direct efficiency measurement below that level, would falsify the 'first-year like-sign' claim while leaving the other reach predictions testable against the differential spectra.","tokens_in":18336,"feed_emoji":"⚛️","tokens_out":13505,"duration_ms":108220,"temperature":0.7,"pith_summary":"This paper predicts the rates and kinematic reach of inclusive charm- and beauty-pair production in electron–proton and electron–nucleus collisions at the beam energies and luminosities planned for the first years of the Electron-Ion Collider. It uses leading-order QCD with modern proton and nuclear parton distribution functions to map the accessible ranges in photon virtuality, heavy-quark transverse momentum, and rapidity, and converts the cross sections into expected event yields. The goal is to supply a parton-level baseline that future higher-order calculations and the EIC's own data can be tested against, and to quantify how much the EIC will improve on earlier HERA measurements. For nuclear targets, the paper also computes the nuclear modification ratio and identifies the rapidity regions where gluon shadowing, anti-shadowing, and EMC-like suppression should be visible in early data.","feed_headline":"Early EIC runs will map charm and beauty beyond HERA","feed_subtitle":"A parton-level baseline predicts where gluon shadowing and anti-shadowing show up in nuclei","key_machinery":"The machinery is collinear-factorisation QCD at leading order, with the heavy quarks treated as massive final-state particles: a three-flavour scheme for charm and a four-flavour scheme for beauty, so the dominant production channel is virtual-photon–gluon fusion, $\\gamma^* g \\to Q\\bar{Q}$. Cross sections are generated with an automated leading-order matrix-element tool, using a common central scale $\\mu_0 = H_T/2$ and a nine-point scale variation to estimate the theoretical uncertainty. The paper's observability criterion is the event-count formula $N_{\\rm ev} = (d\\sigma/dX) \\times f \\times \\epsilon \\times \\Delta X \\times L$, with a threshold of 20 events per bin; the detection efficiencies (3% for charm, 5% for beauty) and fragmentation fractions are taken from other experiments or detector studies because detailed EIC performance estimates are not yet available. The nuclear modification ratio $R^{eA}$, computed with nuclear PDFs relative to a proton baseline at the same collision energy, is the key ratio that turns a measured eA/ep cross-section ratio into a statement about the nuclear gluon distribution.","core_discovery":"The paper's central claim is that early EIC running will make heavy-flavour electroproduction a high-statistics observable over a much wider kinematic range than HERA achieved, and that the same measurements in electron–nucleus collisions will expose the nuclear gluon distribution's main features. In ep collisions at $\\sqrt{s_{ep}}=72$ GeV, charm production is predicted to remain observable up to $Q^2 \\approx 400$ GeV$^2$ and across the full rapidity range, while beauty production reaches $Q^2 \\approx 175$ GeV$^2$; in e–Au collisions at $\\sqrt{s_{eN}}=63$ GeV, the corresponding reach is $Q^2 \\approx 224$ GeV$^2$ for charm and $Q^2 \\approx 75$ GeV$^2$ for beauty. The computed nuclear modification ratio $R^{eA} = (1/A)\\,\\sigma_{eA}/\\sigma_{ep}$, evaluated differentially in heavy-quark rapidity, shows a clear pattern: suppression at backward rapidity from gluon shadowing, an anti-shadowing rise of order 7–10% near central and forward rapidity, and an EMC-like downturn at the most forward rapidity. Because beauty is heavier, it probes larger gluon momentum fractions and mainly samples the anti-shadowing and EMC regions rather than the small-$x$ shadowing regime, giving the two flavours complementary sensitivity. The paper presents these results as a parton-level baseline, with kinematic reach and statistical precision quantified through observability lines defined by a 20-event threshold at the projected early-science luminosities.","pith_inferences":["Applying the paper's own event-count formula to its same-sign $D^0D^0$ cross section (about 1.02 pb), luminosity (5.85 fb$^{-1}$), and 3% efficiency gives roughly 5 events, below the paper's 20-event threshold; this suggests the stated 'first-year' like-sign claim would require a higher efficiency or a lower threshold than assumed.","The parton-level predictions could be turned into a sharper test by recomputing $R^{eA}$ with alternative nuclear PDF sets, which the paper notes are equally applicable; the spread among sets would show where early EIC data best discriminate nuclear gluon models.","Once real EIC detector efficiencies are known, the observability lines can be redrawn and the predicted differential cross sections compared directly with data, turning the baseline into a quantitative test of leading-order collinear factorisation in a multi-scale regime."],"forward_implications":["Charm electroproduction in ep at 72 GeV will be observable across essentially the full planned rapidity range and up to $Q^2 \\approx 400$ GeV$^2$, extending the $x_{\\rm Bj}$ reach well beyond what HERA's statistics permitted.","Beauty electroproduction, although suppressed by the larger quark mass, remains observable up to $Q^2 \\approx 175$ GeV$^2$ in ep and up to $y_b \\approx 2.5$ in e–Au, with narrower scale uncertainties than charm.","The $R^{eA}$ ratio for charm is predicted to show gluon shadowing at backward rapidity, anti-shadowing at central/forward rapidity, and an EMC-like downturn at the most forward rapidity for copper, silver, and gold targets.","Because beauty probes larger gluon momentum fractions at the same rapidity, its $R^{eA}$ stays near unity or mildly suppressed, complementing charm's sensitivity to the small-$x$ shadowing region.","Double-charm-pair electroproduction offers a new observable: opposite-sign $D^0\\bar{D}^0$ pairs at the nanobarn level and same-sign $D^0D^0$ pairs at the picobarn level, a channel that was beyond HERA's reach."],"supporting_citations":[{"why":"Supplies the proton PDF set used for the ep and eA cross-section baselines.","marker":"[15]"},{"why":"Supplies the nuclear PDF set used for all eA cross sections and the $R^{eA}$ ratios.","marker":"[27]"},{"why":"Provides the automated leading-order matrix-element generator used to compute the cross sections.","marker":"[64]"},{"why":"Defines the EIC science requirements and detector concepts that set the kinematic and luminosity context.","marker":"[25]"},{"why":"Provides the early-science beam-energy and luminosity matrix used for the event-yield estimates.","marker":"[65]"},{"why":"Gives the HERA combined heavy-flavour measurements that serve as the statistical benchmark the EIC is expected to surpass.","marker":"[5]"},{"why":"Supplies the $c \\to D^0$ fragmentation fraction used for the D-meson yield projections.","marker":"[70]"},{"why":"Motivates the 3% charm detection efficiency used in the observability lines.","marker":"[71]"},{"why":"Supplies the $b \\to B^0$ fragmentation fraction used for the beauty yield projections.","marker":"[72]"},{"why":"Motivates the electron and heavy-quark pseudorapidity acceptance cuts applied to the parton-level predictions.","marker":"[69]"}],"fun_headline_variants":["EIC early data to probe charm and beauty further than HERA","Early EIC runs to chart gluon shadows with heavy quarks","Charm and beauty at EIC: a wider window than HERA","EIC heavy flavors to expose nuclear gluon effects","Early EIC: heavy-quark reach beyond HERA's limits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The event-yield and observability claims rest on assumed D0 and B0 detection efficiencies (3% and 5%) and fragmentation fractions borrowed from other experiments, since detailed EIC detector performance is not yet available.","fun_headline_variants_meta":{"raw":{"variants":["EIC early data to probe charm and beauty further than HERA","Early EIC runs to chart gluon shadows with heavy quarks","Charm and beauty at EIC: a wider window than HERA","EIC heavy flavors to expose nuclear gluon effects","Early EIC: heavy-quark reach beyond HERA's limits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000657,"raw_usage":{"total_tokens":3101,"prompt_tokens":1136,"completion_tokens":1965,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":752,"completion_tokens_details":{"reasoning_tokens":1876}},"tokens_in":752,"tokens_out":1965,"duration_ms":12338,"temperature":1.0,"reasoning_tokens":1876,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:09:34.425907+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete check is to apply the paper's own numbers to the same-sign D0D0 channel: with $\\sigma \\approx 1.02$ pb, $L = 5.85$ fb$^{-1}$, and $\\epsilon = 3\\%$, the yield is about 5 events, below the 20-event threshold used elsewhere, so reaching 20 events would require $\\epsilon \\approx 5.8\\%$; early EIC runs that record fewer than 20 same-sign pairs, or a direct efficiency measurement below that level, would falsify the 'first-year like-sign' claim while leaving the other reach predictions testable against the differential spectra.","supporting_citations":[{"cited_title":"2025 EIC-France Workshop: Physics Highlights and Perspectives,","cited_arxiv_id":null,"evidence_quote":"Provides the early-science beam-energy and luminosity matrix used for the event-yield estimates."},{"cited_title":"Open Heavy Flavor Studies for the ECCE Detector at the Electron Ion Collider,","cited_arxiv_id":null,"evidence_quote":"Motivates the 3% charm detection efficiency used in the observability lines."},{"cited_title":"Recent open heavy flavor studies for the Electron-Ion Collider","cited_arxiv_id":"2501.18044","evidence_quote":"Motivates the electron and heavy-quark pseudorapidity acceptance cuts applied to the parton-level predictions."}],"review_version":1}