{"id":"4a619916-15d2-4003-9e39-56e051bf1bee","arxiv_id":"2411.16384","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A single set of NRQCD color-octet matrix elements fitted to LHC J/psi and eta_c data predicts a wide range of quarkonium observables, including Upsilon(3S) production and low-momentum photoproduction, with residual failures at low pT hadroproduction and high z.","lead":"This paper fits three nonperturbative NRQCD parameters to LHC J/psi and eta_c data, then uses them to predict many other quarkonium production processes. The same parameters describe J/psi polarization, very high-momentum production, and unexpectedly low-momentum photoproduction, but fail at low hadroproduction momentum and high inelasticity.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"High-pT extrapolation to 360 GeV rests on an indirect inference that log(m_c^2/p_T^2) resummation is negligible for the fitted LDMEs; the resummed calculation should be repeated with the authors' own LDMEs.","rationale":"The central claim of broad LDME universality is supported by several independent predictions, and the paper is transparent about its failure regions. I focus on the high-pT resummation issue because the ATLAS 360 GeV prediction is the most extreme extrapolation and the only place where the authors explicitly rely on an indirect argument rather than a calculation. A failure there would remove one of the headline successes, while the 30% theory-error issue mainly affects the statistical interpretation and the low-pT photoproduction success has large experimental uncertainties. The proposed check is directly feasible because the resummation framework exists in Ref. [53] and the authors already have fixed-order codes for the same observables. Since the reader already identified this as the weakest assumption and the concern does not change the overall conditional assessment, I recommend no change to the verdict.","tokens_in":12241,"tokens_out":8304,"duration_ms":86051,"concrete_test":"Take the resummed high-pT formalism from Ref. [53] (or equivalent) and recompute the pT-differential J/ψ cross section for ATLAS kinematics (sqrt(s) = 13 TeV, |y| < 0.75) using each of the three LDME sets from Table I. Compare the resummed curves to the fixed-order NLO curves and to the ATLAS data. Quantify the shift relative to the yellow band for pT > 120 GeV; if the shift is smaller than the band and the central line remains within data, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is in the section 'Prediction of J/ψ hadroproduction at very high pT': the paper claims that the fitted LDMEs describe ATLAS data up to pT = 360 GeV while using fixed-order NLO without resummation of log(m_c^2/p_T^2) terms. The justification is an inference from Fig. 3 of Ref. [53] that the effect is negligible because the Chao et al. LDMEs [16] 'do not significantly differ' from the authors' default LDMEs. This inference is not self-evidently valid: the size of the resummation correction depends on the full Fock-state decomposition, and the fitted LDMEs vary strongly across the three scale choices in Table I (e.g. ⟨O(3S1[8])⟩ ranges from 0.59 to 1.38 × 10^-2 GeV^3 and ⟨O(3P0[8])⟩/m_c^2 from 0.70 to 3.27 × 10^-2 GeV^3), so a statement about another LDME set does not directly bound the correction for the authors' set. Moreover, agreement of a resummed curve with data does not by itself show that the fixed-order curve is unchanged; it could be that resummation shifts the prediction and the conclusion would require re-evaluation. If the resummation shift exceeds the yellow uncertainty band at pT > 120 GeV, the headline claim of description up to 360 GeV would lose its support.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper performs a combined NLO NRQCD fit of three J/psi color-octet LDMEs to 42 data points from CMS J/psi and LHCb eta_c hadroproduction, using three scale choices (mT/2, mT, 2mT) and a global 30% theory error. The fitted LDMEs are then used to predict J/psi polarization, ATLAS high-pT J/psi production up to 360 GeV, low-pT LHCb J/psi production, Upsilon(3S) production via pNRQCD-based relations, J/psi+Z production, gamma gamma scattering at LEP, and HERA photoproduction in five z bins. The central claim is that a single set of CO LDMEs can describe a wide range of quarkonium observables, with failures only at low/medium pT hadroproduction and at large photoproduction inelasticity z; the paper interprets these patterns as evidence for NRQCD factorization and LDME universality, while also identifying regions where factorization appears to break down.","tokens_in":121,"tokens_out":1896,"duration_ms":81348,"significance":"If the central claim holds, this would be a strong demonstration of LDME universality across hadroproduction, photoproduction, two-photon scattering, and bottomonium production. The paper's fit-and-predict procedure, which propagates scale variations through the fit and into predictions via an envelope over three scale choices, is a methodological improvement over fits that ignore scale dependence or treat it as a purely global uncertainty. The successful Upsilon(3S) prediction using pNRQCD relations is a particularly nontrivial cross-check. The paper also provides explicit falsifiable statements: low-pT hadroproduction and z>0.6 photoproduction are not described, and J/psi+Z production shows a possible discrepancy. These positive and negative results together make the paper valuable for the quarkonium-production community, provided the high-pT extrapolation and the treatment of theory uncertainty are adequately justified.","major_comments":[{"comment":"The claim that the fitted LDMEs describe ATLAS data up to pT = 360 GeV rests on fixed-order NLO without resummation of log(m_c^2/pT^2) terms, justified only by an inference from Fig. 3 of Ref. [53] that the effect is negligible because the Chao et al. LDMEs 'do not significantly differ' from the authors' default LDMEs. This inference is not demonstrated quantitatively: Table I shows that the fitted LDMEs vary strongly across the three scale choices (e.g., <O(3S1[8])> ranges from 0.59 to 1.38 x 10^-2 GeV^3 and <O(3P0[8])>/m_c^2 from 0.70 to 3.27 x 10^-2 GeV^3), so the size of resummation corrections, which depends on the Fock-state decomposition, cannot be reliably inferred from another LDME set. The authors should either compute the resummed result with their own LDMEs or provide a quantitative bound (e.g., compare the resummed and fixed-order predictions for the Chao et al. set and show the difference is small relative to the yellow band for the authors' LDME range). Without this, the headline high-pT success is a load-bearing unresolved assumption.","section":"Prediction of J/psi hadroproduction at very high pT (p. 3, Fig. 2c)"},{"comment":"The global 30% theory error is chosen by hand and applied uniformly, yet the resulting chi^2/d.o.f. values (0.21-0.34) are far below 1. This indicates that the data are fitted with substantially more freedom than the quoted errors suggest, and the yellow bands in the prediction plots are correspondingly wide. The central claim of 'good description' is therefore weakened: a chi^2/d.o.f. of 0.2-0.3 with an ad hoc 30% error does not by itself demonstrate that the model is correct, only that it is not excluded. The authors should justify the 30% value (e.g., by comparing relative sizes of relativistic corrections estimated from lattice or potential-model calculations) and discuss how the conclusions change if the theory error is reduced or eliminated; the current presentation makes the quality of the fit appear better than the input assumptions warrant.","section":"Fit and predict procedure, beginning of 'Introduction and overview' and 'Fit results' (p. 2, Table I)"},{"comment":"The paper states that the two highest pT bins lie about two experimental standard deviations below data and argues that underestimated DPS is 'unlikely, albeit not impossible,' based on the pocket formula. However, the text then acknowledges that the pocket formula itself is 'subject to debate.' This is an internal tension: if the DPS estimation is uncertain, the discrepancy with the SPS-subtracted data cannot be crisply interpreted, and the statement 'J/psi + Z production remains intricate to interpret' is appropriate but should be presented as a limitation of the prediction test, not as a confirmed failure or success. The authors should clarify whether the conclusion on J/psi+Z is robust under reasonable variations of sigma_eff and the DPS procedure.","section":"Prediction of J/psi + Z production (Fig. 2f)"},{"comment":"The paper explains the failure at low pT by the sign change of the 3PJ[8] SDCs below pT ~ 7 GeV, which turns the cancellation between 3S1[8] and 3PJ[8] into an amplification. This is a plausible mechanism, but the figure appears to show that the theoretical prediction is orders of magnitude above the data at pT = 1-3 GeV. The authors should state explicitly whether this failure is a violation of NRQCD factorization or an expected artifact of the fixed-order calculation in the endpoint region pT^2 << 4m_c^2, and whether resummation or shape-function effects (as cited in Refs. [55-57]) would be expected to cure the discrepancy. As written, the paper leaves the interpretation ambiguous.","section":"Prediction of low-pT hadroproduction (Fig. 2d)"}],"minor_comments":[{"comment":"The notation <O(3PJ[8])> in Fig. 2(d)-(f) instead of <O(3P0[8])> is confusing because the text consistently uses the J-summed <O(3P0[8])>/m_c^2 as the fit parameter; please define the plotted quantity.","section":"Throughout"},{"comment":"The sentence 'The latter approach thus appears to be the most promising fit strategy, and is the one we adopt here' is slightly misleading because the authors' approach differs from Refs. [16,17] in the treatment of scale variations; this difference should be stated in the same paragraph.","section":"Introduction and overview, p. 2"},{"comment":"Reference [53] is cited as the source of the resummed curve, but the text does not specify whether Ref. [53] uses the same CS LDME and feeddown treatment as the present paper; a brief note on compatibility would help.","section":"Prediction of J/psi hadroproduction at very high pT (p. 3)"},{"comment":"The LHCb polarization data extend down to pT = 2 GeV, while the fit is based on pT >= 7 GeV (eta_c) and pT >= 10 GeV (CMS J/psi); the paper should state whether the low-pT polarization points were included in the fit or are predictions, since the text is ambiguous.","section":"Prediction of J/psi polarization (p. 3, Fig. 2b)"},{"comment":"The abstract mentions 'the summary reveals an interesting pattern,' but the main text's Discussion section is short; a short table listing each observable, the data set, and the verdict (described/not described) would make the pattern more accessible.","section":"Summary, p. 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is competent and addresses an important question, but the high-pT prediction to 360 GeV is currently supported by an indirect inference about resummation that the authors should verify with their own LDMEs. The global 30% theory error also makes the very low chi^2/d.o.f. less informative than it appears. These are fixable within the scope of a revision, so I recommend major revision rather than rejection. The paper would be strengthened by a direct resummation calculation or a quantitative check of the size of the log(m_c^2/pT^2) corrections for the authors' LDME sets."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should read this paper if you care about the NRQCD LDME universality question. The authors fit three J/psi color-octet LDMEs to 42 data points from CMS J/psi and LHCb eta_c production, then use the same LDMEs to predict a wide menu of observables: J/psi polarization, ATLAS hadroproduction up to pT = 360 GeV, Upsilon(3S) production via pNRQCD relations, J/psi+Z, and low-pT gamma-gamma and gamma-p production. The headline result is that one set of LDMEs describes most of these, while failing in the regions the authors explicitly flag: low/medium-pT hadroproduction and z > 0.6 photoproduction. That honesty about failures is a real strength and makes the paper more trustworthy than the average fit paper.\n\nWhat is actually new: the combined CMS J/psi + LHCb eta_c fit with three separate scale choices and an envelope of predictions is a genuine methodological step, even if the general strategy is close to Refs. [16,17]. The gamma-gamma and gamma-p predictions down to pT = 1 GeV and the Upsilon(3S) prediction from the pNRQCD relations are new tests. The small chi2/dof is less impressive than it looks because a global 30% theory error is added to the fit by hand, which inflates the fit quality. The qualitative agreement across many independent observables is the central evidence, and it does not depend much on that error choice.\n\nThe softest spot is the high-pT tail. To claim description up to 360 GeV at fixed order, the authors rely on an indirect inference from Ref. [53] that resummation of log(m_c^2/pT^2) terms is negligible because the Chao et al. LDMEs used there are similar to theirs. The stress-test note has this right: Table I shows their fitted LDMEs vary by factors of two across scale choices, and the resummation correction depends on the full Fock-state decomposition, so a statement about another LDME set is not a direct bound on their set. This is not fatal - the universality claim is supported by many other panels - but it is load-bearing for the specific 'up to 360 GeV' statement and should be fixed by actually running the resummed calculation with their own LDMEs.\n\nMinor points: no code or data release; PDF uncertainties are ignored, though the scale envelope likely dominates; and the Upsilon(3S) prediction uses self-cited pNRQCD relations, but the cited derivation is formal and the prediction is a legitimate test. The J/psi+Z discrepancy is discussed honestly, including the DPS caveats.\n\nBottom line: this deserves a serious referee. Send it to review, with a request for the resummed check and a sensitivity study of the 30% theory error. The paper may well be an important piece of evidence for LDME universality, but I would not cite the 360 GeV claim as settled until the resummation question is closed.","headline":"A broad, honest NRQCD universality study whose main claim holds up, but the 360 GeV extrapolation rests on an indirect resummation argument that should be pinned down.","tokens_in":13193,"tokens_out":2910,"would_cite":true,"duration_ms":25511,"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":"One set of nonrelativistic-QCD matrix elements, fitted to 42 data points, predicts quarkonium production across hadron, photon, and photon-photon collisions.","keywords":["nonrelativistic QCD factorization","color-octet LDMEs","quarkonium production","J/psi","Upsilon production","eta_c","LDME universality","next-to-leading order"],"falsifier":"Compute the $p_T$-differential $J/\\psi$ cross section at $\\sqrt{s}=13$ TeV up to $p_T=360$ GeV using the LDMEs fitted in this paper, now including resummation of the $\\log(m_c^2/p_T^2)$ terms; the claim that the fit describes the ATLAS data fails if the resummed prediction at $p_T>100$ GeV moves outside the paper's quoted uncertainty envelope by more than the experimental precision.","tokens_in":11946,"feed_emoji":"⚛️","tokens_out":21028,"duration_ms":149727,"temperature":0.7,"pith_summary":"The paper tests whether the long-distance matrix elements (LDMEs) of nonrelativistic QCD (NRQCD) factorization are universal: the same numbers describing how a heavy quark-antiquark pair becomes a quarkonium, regardless of how the pair is created. It fits three color-octet LDMEs to 42 data points of $J/\\psi$ and $\\eta_c$ production at the LHC, then uses those fixed numbers to predict a wide range of other quarkonium observables. The predictions work: $J/\\psi$ hadroproduction up to $p_T = 360$ GeV, $\\Upsilon(3S)$ production via potential-NRQCD relations, and $J/\\psi$ production in $\\gamma\\gamma$ and $\\gamma p$ collisions down to $p_T = 1$ GeV for inelasticity $z < 0.6$. The notable failures are low-$p_T$ hadroproduction and photoproduction at large $z$, which the authors trace to kinematic endpoint regions where fixed-order perturbation theory is expected to lose convergence. If the universality claim holds, a single set of LDMEs becomes a predictive tool for quarkonium yields at future colliders, and the residual discrepancies define where higher-order or resummed calculations are needed.","feed_headline":"One 42-point fit predicts quarkonium yields across collision types","feed_subtitle":"The same matrix elements that fit 42 data points also predict Upsilon(3S) yields and low-pT photoproduction.","key_machinery":"The central object is the set of three $J/\\psi$ color-octet LDMEs, $\\langle \\mathcal{O}^{J/\\psi}({}^3S_1^{[8]})\\rangle$, $\\langle \\mathcal{O}^{J/\\psi}({}^1S_0^{[8]})\\rangle$ and $\\langle \\mathcal{O}^{J/\\psi}({}^3P_J^{[8]})\\rangle/m_c^2$, the nonperturbative parameters of NRQCD factorization that encode how a heavy quark-antiquark pair in a color-octet state evolves into a physical quarkonium. The machinery carrying the argument is the fit-and-predict procedure: three independent NLO least-squares fits are performed with renormalization and factorization scales set to $m_T/2$, $m_T$, and $2m_T$, the predicted error bands from each fit are superposed, and their envelope is the quoted theory uncertainty. The LHCb $\\eta_c$ data directly constrain the ${}^1S_0^{[8]}$ matrix element, breaking the degeneracy that plagues fits using only $J/\\psi$ yields. Heavy-quark spin symmetry relates the $J/\\psi$ and $\\eta_c$ LDMEs, and pNRQCD relations transfer the fitted charmonium LDMEs to $\\Upsilon(nS)$ production. The cancellation between the large positive ${}^3S_1^{[8]}$ and large negative ${}^3P_J^{[8]}$ contributions—not fine-tuning, because NLO mixing makes only their sum physical—is what keeps the high-$p_T$ $J/\\psi$ prediction stable.","core_discovery":"The authors claim that a 'fit-and-predict' procedure using a single set of $J/\\psi$ color-octet long-distance matrix elements (LDMEs)—determined by a next-to-leading-order (NLO) fit to 42 points of CMS $J/\\psi$ and LHCb $\\eta_c$ data—can quantitatively describe quarkonium production across very different collision environments. The key new results are that the same LDMEs reproduce ATLAS $J/\\psi$ hadroproduction up to $p_T = 360$ GeV, predict $\\Upsilon(3S)$ production through potential-NRQCD (pNRQCD) derived relations, and describe $J/\\psi$ production in $\\gamma\\gamma$ and $\\gamma p$ collisions down to $p_T = 1$ GeV, provided the inelasticity $z$ is below about 0.6. The paper treats scale variation systematically by performing three separate fits at scales $m_T/2$, $m_T$, and $2m_T$ (with $m_T = \\sqrt{p_T^2 + 4m_c^2}$) and taking the envelope of the resulting error bands as the theory uncertainty. On this basis the authors argue that LDME universality, a core conjecture of NRQCD factorization, survives these nontrivial tests, while the remaining discrepancies (low-$p_T$ hadroproduction, $z>0.6$ photoproduction, and possibly the highest-$p_T$ $J/\\psi+Z$ bins) coincide with kinematic regions where fixed-order calculations are expected to lose convergence.","pith_inferences":["If the observed universality extends to double-quarkonium channels, the same LDMEs should describe high-pT J/psi plus J/psi production at the LHC; this is a testable corollary the paper does not pursue.","The success of fixed-order NLO in z < 0.6 photoproduction down to pT = 1 GeV suggests the color-octet mechanism dominates even where the transverse momentum is below the charm mass, a regime where one might have expected the factorization to break; a dedicated high-statistics measurement at the EIC could confirm or refute this boundary.","The paper's justification for neglecting resummation at very high pT rests on a comparison with a different LDME set; a direct resummed calculation with the paper's own LDMEs would either close this gap or reveal that the high-pT agreement is partly accidental.","The sharpness of the z ~ 0.6 boundary in photoproduction, independent of pT, suggests the missing physics is governed by inelasticity rather than transverse momentum, pointing to universal endpoint behavior that could be modeled with shape functions or soft-gluon resummation."],"forward_implications":["The same fitted LDMEs can be used to predict quarkonium yields at the EIC and the high-luminosity LHC, including J/psi photoproduction in the z < 0.6 region down to pT = 1 GeV.","The successful Upsilon(3S) prediction supports the pNRQCD relations linking charmonium and bottomonium matrix elements, so Upsilon data can be added as constraints in future global LDME fits.","The fit-and-predict scale-envelope procedure, by correlating scale variations with LDME uncertainties, offers a template for other NRQCD analyses that currently ignore scale variation or treat it as a global theory error.","The residual discrepancies in low-pT hadroproduction and high-z photoproduction identify specific kinematic windows where resummed or shape-function-improved calculations must be tested, and where new measurements would be most discriminating."],"supporting_citations":[{"why":"Defines NRQCD factorization and the long-distance matrix elements under test.","marker":"[1]"},{"why":"Supplies the pNRQCD relations (3.47)-(3.48) that turn fitted charmonium LDMEs into Upsilon(nS) predictions.","marker":"[4]"},{"why":"Introduced the fit strategy of including LHCb eta_c data to constrain the 1S0[8] LDME, adopted here.","marker":"[16]"},{"why":"Provides LHCb eta_c production data at 7 and 8 TeV used as fit input.","marker":"[23]"},{"why":"Provides LHCb eta_c production data at 13 TeV used as fit input.","marker":"[24]"},{"why":"Provides CMS J/psi hadroproduction data used as the primary hadroproduction fit input.","marker":"[25]"},{"why":"ATLAS J/psi data up to pT=360 GeV that the fitted LDMEs are shown to predict.","marker":"[26]"},{"why":"Resummed calculation cited as evidence that log(m_c^2/p_T^2) resummation is negligible for the high-pT prediction.","marker":"[53]"}],"fun_headline_variants":["One fit predicts quarkonium yields from pp to gamma-p","42-point NRQCD fit reproduces quarkonium across processes","Universal quarkonium matrix elements pass NLO test","Single set of matrix elements predicts quarkonium from hadrons to photons","Quarkonium universality: one fit to rule many collisions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that fixed-order next-to-leading-order calculations, with the three-scale envelope and a flat 30% theory error, adequately capture the true theoretical uncertainty; in particular, the paper infers rather than computes that resummation of $\\log(m_c^2/p_T^2)$ terms is negligible for its own LDMEs at very high $p_T$, citing a resummed result for a different LDME set.","fun_headline_variants_meta":{"raw":{"variants":["One fit predicts quarkonium yields from pp to gamma-p","42-point NRQCD fit reproduces quarkonium across processes","Universal quarkonium matrix elements pass NLO test","Single set of matrix elements predicts quarkonium from hadrons to photons","Quarkonium universality: one fit to rule many collisions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000213,"raw_usage":{"total_tokens":1815,"prompt_tokens":1102,"completion_tokens":713,"prompt_tokens_details":{"cached_tokens":1024},"prompt_cache_hit_tokens":1024,"prompt_cache_miss_tokens":78,"completion_tokens_details":{"reasoning_tokens":630}},"tokens_in":78,"tokens_out":713,"duration_ms":280443,"temperature":1.0,"reasoning_tokens":630,"cache_read_input_tokens":1024,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:10:23.240021+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the $p_T$-differential $J/\\psi$ cross section at $\\sqrt{s}=13$ TeV up to $p_T=360$ GeV using the LDMEs fitted in this paper, now including resummation of the $\\log(m_c^2/p_T^2)$ terms; the claim that the fit describes the ATLAS data fails if the resummed prediction at $p_T>100$ GeV moves outside the paper's quoted uncertainty envelope by more than the experimental precision.","supporting_citations":[{"cited_title":"Measurement of $J/\\psi$ production in association with a $W^\\pm$ boson with $pp$ data at 8 TeV","cited_arxiv_id":"1909.13626","evidence_quote":"Provides CMS J/psi hadroproduction data used as the primary hadroproduction fit input."}],"review_version":1}