{"id":"213d1583-221f-4526-a9c4-178544bde93c","arxiv_id":"2507.20654","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The gluon-gluon fusion process gg -> J/psi + c cbar is computed at NLO QCD, with large corrections, nearly unpolarized J/psi, and improved agreement with LHC data, while a normalization gap of about an order of magnitude persists.","lead":"This paper reports the first complete next-to-leading-order QCD calculation of the process where a J/psi meson is produced together with a charm-anticharm pair in proton-proton collisions. The finding: the higher-order corrections are large, change the predicted polarization of the J/psi, and bring the color-singlet theory closer to LHC data, though a gap of about a factor of ten remains.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Missing cutoff-independence validation for the two-cutoff real-emission slicing: Eq. (6)-(7) show no δ_s, δ_c dependence or independent check, so the reported NLO rates and λθ values are not yet numerically secured.","rationale":"The reader's weakest_assumption—numerical independence of the arbitrary cutoffs in the two-cutoff phase-space slicing method—is exactly the load-bearing concern I identify. The paper is a technically demanding NLO calculation with a nontrivial analytic IR-cancellation check in Eq. (7), and the authors do acknowledge a residual gap to data. However, the finite part of the real-emission contribution is where the method's cutoffs must cancel, and the manuscript provides no numerical evidence of that cancellation. The reported enhancement by several times and the nearly unpolarized polarization pattern are quantitative claims that cannot be validated from the printed equations alone. I also note the secondary issue that the 'first complete calculation' claim is not reconciled with Refs. [8] and [23], which the paper itself cites as earlier NLO studies of J/ψ + c cbar production; this is a real bibliographic concern, but it is secondary to the numerical reliability question because even a perfectly novel calculation would still need the cutoff-independence check. Accordingly, the reader's CONDITIONAL verdict remains appropriate: the paper should be accepted only if the authors supply cutoff-dependence tables or an independent integration check, and clarify the novelty statement relative to prior NLO work.","tokens_in":8064,"tokens_out":14052,"duration_ms":145799,"concrete_test":"Recompute the NLO dσ/dpT and λθ at a representative point (e.g., pT ≈ 50 GeV at √s = 13 TeV) with δ_s and δ_c varied over at least three orders of magnitude (e.g., 10^-3, 10^-4, 10^-5) while keeping all other settings fixed; require that the total cross section and each polarization component dσ11 and dσ00 change by less than about 1%. As a further check, reproduce the hard non-collinear integral with an independent adaptive Monte Carlo integrator and verify that the sum in Eq. (6) is cutoff-independent. If no plateau exists, the reported numbers are not trustworthy.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical claims—NLO enhancement by several times and nearly unpolarized J/ψ—rest on the finite remainder of the 2→4 real-emission contribution in Eq. (6), where the real corrections are separated with the two-cutoff phase-space slicing method into soft (S), hard-collinear (HC), and hard non-collinear regions. In this method, the analytically integrated soft and collinear pieces carry logarithms of the cutoffs δ_s and δ_c, while the hard non-collinear part is obtained by numerical integration over the 4-body phase space; only after summing all pieces does the δ dependence cancel in the limit δ→0. The paper's only stated check, Eq. (7), verifies cancellation of the IR poles between virtual and real amplitudes; it says nothing about the finite remainders, which are exactly where cutoff leakage or an incorrect 4-body integration would appear. No numerical values of δ_s and δ_c are given, and no demonstration is shown that σ, dσ/dpT, or λθ are stable when the cutoffs are varied. Without this, the substantial enhancement, the near-zero polarization, and the comparison with ATLAS/CMS data in Figs. III.1–III.3 could be numerical artifacts rather than physics. A separate bibliographic point is that the assertion of a 'first complete' NLO calculation is not reconciled with the earlier NLO literature cited as Refs. [8] and [23]; resolving that would also require comparing to those results.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter reports what it claims to be the first complete next-to-leading-order QCD calculation, at O(α_s^5), of the hadroproduction process gg → J/ψ + c cbar in the NRQCD color-singlet framework. The authors compute the J/ψ transverse-momentum distribution and the polarization parameter λθ at LO and NLO, compare their results with ATLAS data at 13 TeV and CMS data at 7 TeV, include ψ(2S) and χcJ feeddown estimates, and compare the associated-production channel with charm-quark fragmentation. The main findings are that the NLO corrections enhance the color-singlet yield substantially, that the J/ψ + c cbar channel is nearly unpolarized over a wide pT range, and that including this process improves agreement with data while a gap of about one order of magnitude remains.","tokens_in":8365,"tokens_out":6525,"duration_ms":70916,"significance":"If the numerical results are correct, this is a significant step forward for NRQCD phenomenology: it extends the NLO program for inclusive quarkonium hadroproduction to a 2→4 partonic final state, and it shows that this channel can alter the predicted polarization pattern, which is directly testable with existing LHC data. The paper has notable strengths: the analytic IR-cancellation check in Eq. (7), the use of modern amplitude techniques such as integrand reduction and NeatIBP, the use of external inputs for PDFs, fragmentation functions, and the CS LDME, and the explicit caveat that a residual factor-of-ten gap with data remains. However, the absence of any numerical validation of the two-cutoff phase-space slicing, the lack of numerical tables, and the vague treatment of the ψ(2S) feeddown fit mean that the quantitative claims are not yet fully secured. The central prediction is not tuned to the data it is compared with, and the polarization prediction is falsifiable, which is to the paper's credit.","major_comments":[{"comment":"The central numerical results in §III rest on the two-cutoff phase-space slicing separation in Eq. (6), but the manuscript gives no evidence that the finite remainders are independent of the slicing parameters. The analytic check in Eq. (7) verifies cancellation of the 1/ε IR poles between virtual and real corrections; it does not test for cutoff leakage or errors in the numerical 2→4 body phase-space integration, which affect exactly the finite parts that produce the reported K-factors and λθ values. Please report the values of δ_s and δ_c used, and show that σ, dσ/dpT, and λθ are stable over a range of cutoffs, or provide an independent subtraction or Monte Carlo validation. Without this, the claimed \"several times\" enhancement and the nearly unpolarized polarization could be numerical artifacts.","section":"§II, Eqs. (6)–(7)"},{"comment":"The paper states in §I that the full NLO QCD corrections to gg→J/ψ+c cbar \"have never been computed\" and in §IV that it presents \"the first calculation,\" yet §I cites Refs. [8] and [23] as earlier NLO studies of J/ψ+c cbar associated production. This apparent conflict must be resolved: please specify precisely what those earlier calculations did and did not include, and compare at least one common observable, such as the total NLO cross section or pT distribution, to demonstrate that the present calculation is new and consistent with the earlier results. Without this, the novelty claim and the assessment of the size of the NLO corrections cannot be evaluated.","section":"§I and §IV"},{"comment":"The ψ(2S) feeddown is said to be \"fitted up to 120 GeV using ATLAS data [35],\" but the fit function, the fitted parameters, and the propagation of the associated uncertainties into the \"direct J/ψ\" data are not described. Because the same ATLAS sample is then used to compare both direct J/ψ and prompt ψ(2S) in Fig. III.1, the comparison is partly circular for the feeddown component. Please specify the fitting procedure and its uncertainty, and ideally show that the main NLO conclusions are insensitive to reasonable variations of this feeddown subtraction.","section":"§III, Fig. III.1 and feeddown procedure"}],"minor_comments":[{"comment":"The polarization parameter is quoted in the HX frame in Fig. III.2, but Eq. (5) defines the longitudinal axis from the partonic momentum k1; the transformation to the experimental HX frame should be stated explicitly so that the comparison with CMS data is unambiguous.","section":"§III, Eq. (5) and Fig. III.2"},{"comment":"Several typographical errors should be corrected: \"was found plays\" in the abstract, \"caculable\" in §I, \"extracted form\" in §III, \"evquation\" in §III, and \"convolution production\" in the description of Eq. (11).","section":"Throughout"},{"comment":"No numerical values are given for the K-factor, for the integration errors on the plotted cross sections, or for the phase-space slicing cutoffs; since the figures are log-scale, the reader cannot verify the claimed \"several times\" enhancement or the size of the residual gap.","section":"§III, Figs. III.1–III.3"},{"comment":"The charm-quark fragmentation estimate mixes an NLO cross section with a LO fragmentation function plus an LO cross section with a NLO fragmentation function; a sentence explaining why this combination, rather than a fully NLO convolution, is appropriate and does not introduce double counting would be helpful.","section":"§III, Eq. (11)"},{"comment":"Reference [40] is cited only as an arXiv preprint; please update it to the published version, if one exists, and check that all cited branching ratios and feeddown fractions are the most recent values.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The two points I would most want addressed before publication are the missing cutoff-independence validation of the two-cutoff slicing and the unresolved relation of the \"first complete calculation\" claim to Refs. [8] and [23]. I would also encourage the editor to ask for an ancillary or supplementary file containing numerical tables of the cross sections and λθ, since the Letter currently does not permit independent reproduction of the central quantitative results."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a real calculation, not a toy. The authors push through a 2-to-4 body NLO computation for gg -> J/psi + c cbar in the color-singlet channel, with an analytic check that the IR poles cancel among the color-decomposed Born amplitudes. That is a serious technical step, and if the numbers hold, it gives the NRQCD global-fit community a new CS benchmark that could change how much room is left for color-octet matrix elements. The comparison to charm-quark fragmentation is also useful and undercuts the usual assumption that fragmentation reproduces the associated production process.\n\nThe soft spots are about numerical verification, not about the formal framework. The two-cutoff phase-space slicing method has arbitrary parameters delta_s and delta_c; the paper states the analytic cancellation of IR divergences but never shows that the finite remainders -- the NLO enhancement, the pT shapes, the polarization -- are independent of those cutoffs. Without that, the central numerical claims are not yet secured. A table of cross sections at two or three cutoff values, or a plot of dsigma/dpT versus delta, would settle it. There are also no numerical tables at all, just figures, which makes independent checking hard. I would also like a sentence in the introduction that explicitly says what is new relative to Refs. [8] and [23]; as written, the 'first complete' claim sits uneasily next to those earlier NLO papers and the reader shouldn't have to guess whether they computed related but different observables.\n\nThe feeddown treatment is acceptable: fitting the psi(2S) contribution to ATLAS data is a subtraction, not a tune of the J/psi central prediction, and the residual order-of-magnitude gap is honestly acknowledged. The internal consistency of Eq. (7) counts for something, but it is not a substitute for a Monte Carlo sanity check on the 4-body phase space.\n\nBottom line: this deserves a serious referee, roughly in its current form, with a request for supplemental numeric detail. If the cutoff independence holds, the paper will be a useful benchmark for the quarkonium community. I would not cite it myself until the numbers are independently checkable, but I would bring it to the reading group to see how the calculation was set up.","headline":"A technically substantial NLO calculation that likely fills a real gap in quarkonium hadroproduction, but the numerical results are not yet secured without cutoff-independence checks and a clearer novelty statement.","tokens_in":8912,"tokens_out":2493,"would_cite":false,"duration_ms":26449,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Bx","13.85.Ni","14.40.Pq"],"model":"deepseek-v4-flash","headline":"The first complete NLO QCD calculation of gluon fusion to J/psi plus a charm pair shows the NLO corrections multiply the color-singlet yield severalfold and make the polarization nearly unpolarized, matching LHC data.","keywords":["J/psi production","NRQCD factorization","next-to-leading-order QCD","color-singlet channel","quarkonium polarization","associated heavy-quark pair production","hadroproduction","QCD corrections"],"falsifier":"Recompute the NLO cross section and $\\lambda_\\theta$ with substantially different phase-space-slicing cutoffs, or with an independent subtraction scheme; if the results move beyond integration uncertainties, the reported enhancement and near-unpolarized polarization are artifacts of the slicing prescription.","tokens_in":7848,"feed_emoji":"⚛️","tokens_out":12089,"duration_ms":117592,"temperature":0.7,"pith_summary":"The paper addresses a standing puzzle in charmonium production: the color-singlet mechanism alone seemed unable to reproduce measured J/psi yields and polarization. It presents the first complete next-to-leading-order QCD calculation of the gluon-fusion process $gg\\to J/\\psi+c\\bar{c}$ within NRQCD factorization, at $\\mathcal{O}(\\alpha_s^5)$. The result is that the NLO corrections are large, enhancing the leading-order color-singlet yield by several times, so that this associated-production channel dominates over $J/\\psi+g$ among color-singlet mechanisms. The predicted $p_T$ spectrum matches the shape of direct-J/psi data and the polarization stays nearly unpolarized over a wide $p_T$ range, consistent with CMS measurements, while a normalization gap of roughly an order of magnitude remains. If the calculation is right, global fits of color-octet matrix elements need to include this channel, and the J/psi-plus-open-charm signature becomes a clean way to probe the production mechanism.","feed_headline":"J/psi + charm-pair channel explains flat LHC polarization","feed_subtitle":"The full next-to-leading-order cross section jumps severalfold and matches the observed shapes, narrowing a long-standing gap.","key_machinery":"The load-bearing object is the short-distance coefficient for the $2\\to4$ partonic process $gg\\to c\\bar{c}[{}^3S_1]+c\\bar{c}$, evaluated in NRQCD factorization with the color-singlet matrix element $\\langle \\mathcal{O}^{J/\\psi}({}^3S_1)\\rangle$. The calculation combines helicity amplitudes for the finite parts with a projection tensor that isolates the longitudinal polarization, and the finiteness of the result rests on an infrared-cancellation identity: the virtual and real divergences are matched through color-decomposed Born amplitudes and a $3\\times3$ kernel $\\Delta$, which lets the two-cutoff phase-space slicing of the real corrections cancel the one-loop poles. This matching turns a formally divergent expression into a finite cross section and a finite polarization parameter.","core_discovery":"The central discovery is that the $\\mathcal{O}(\\alpha_s^5)$ corrections to $gg\\to J/\\psi+c\\bar{c}$ are not a small refinement. In the color-singlet NRQCD channel they enhance the leading-order cross section by several times and flatten the $p_T$ distribution. The polarization parameter $\\lambda_\\theta$ comes out close to zero from low to high transverse momentum in the helicity frame, matching the mildly varying values observed by CMS, whereas the $J/\\psi+g$ channel predicts strong longitudinal polarization. The paper concludes that associated $J/\\psi+c\\bar{c}$ production is a dominant color-singlet source of prompt $J/\\psi$ and should be included when extracting color-octet long-distance matrix elements from global fits, even though a residual gap of about one order of magnitude between theory and data remains.","pith_inferences":["If the NLO enhancement is as large as reported, the same calculation should be repeated for $\\psi(2S)+c\\bar{c}$ and for bottomonium counterparts; those channels would receive comparable corrections and could be tested with existing data.","The near-zero polarization implies that a global fit that includes this process will shift the fitted color-octet matrix elements downward, which would alter predictions for other observables such as $J/\\psi$ production in photoproduction or in proton-nucleus collisions.","A direct experimental test would be to reconstruct $J/\\psi$ mesons accompanied by open-charm hadrons ($D$ or $\\bar{D}$ mesons) at the LHC; the predicted $p_T$ shape, normalization, and polarization of that subsample separate the associated-production mechanism from fragmentation."],"forward_implications":["The color-singlet contribution to inclusive $J/\\psi$ production is no longer negligible: at NLO, the $J/\\psi+c\\bar{c}$ process dominates over $J/\\psi+g$ in the color-singlet channel and falls off more slowly with $p_T$.","Global fits of color-octet long-distance matrix elements should include this process; without it, octet matrix elements will absorb yield that is actually color-singlet.","The near-unpolarized prediction for this channel brings color-singlet calculations into closer agreement with the CMS polarization measurements across the measured $p_T$ range.","At high $p_T$, charm-quark fragmentation eventually overtakes the associated-production process but with a different shape and normalization, so the two mechanisms can be distinguished experimentally.","A residual gap of roughly one order of magnitude remains after including this channel, so additional prompt-$J/\\psi$ production mechanisms are still required."],"supporting_citations":[{"why":"Supplies the NRQCD factorization that separates the short-distance coefficients computed here from the nonperturbative color-singlet matrix element.","marker":"[1]"},{"why":"Earlier leading-order study of J/psi+c cbar associated production whose finding of a sizable, nearly unpolarized gluon-fusion contribution motivates the full NLO calculation.","marker":"[7]"},{"why":"Earlier NLO work on associated charmonium production whose large corrections set the expectation the present calculation verifies at full O(alpha_s^5).","marker":"[8]"},{"why":"Previous calculation by the authors whose computational strategy, helicity-amplitude handling, and loop treatment the present work largely follows.","marker":"[25]"},{"why":"Provides the two-cutoff phase-space slicing method used to isolate and subtract the soft and collinear singularities of the real-emission contributions.","marker":"[31]"},{"why":"Gives the method used to extract the infrared divergences of the loop amplitudes and check their cancellation against the real corrections.","marker":"[32]"},{"why":"ATLAS prompt J/psi differential cross-section data used as the experimental benchmark and source for the psi(2S) feeddown subtraction.","marker":"[35]"},{"why":"CMS 13 TeV measurement of J/psi polarization that the calculated lambda_theta is compared against and found to match.","marker":"[21]"}],"fun_headline_variants":["NLO charm-pair channel flattens J/psi polarization","J/psi + charm pair: NLO boosts cross section, flattens polarization","Full NLO reduces gap for J/psi with charm pair","Charm-pair channel reshapes J/psi production at NLO","NLO corrections flatten J/psi polarization, shrink gap"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The final cross sections are reliable only if the artificial cutoffs used to separate soft and collinear real emissions leave the result unchanged; the paper states the infrared cancellation analytically but does not demonstrate numerical cutoff independence.","fun_headline_variants_meta":{"raw":{"variants":["NLO charm-pair channel flattens J/psi polarization","J/psi + charm pair: NLO boosts cross section, flattens polarization","Full NLO reduces gap for J/psi with charm pair","Charm-pair channel reshapes J/psi production at NLO","NLO corrections flatten J/psi polarization, shrink gap"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000763,"raw_usage":{"total_tokens":3371,"prompt_tokens":915,"completion_tokens":2456,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":531,"completion_tokens_details":{"reasoning_tokens":2363}},"tokens_in":531,"tokens_out":2456,"duration_ms":15300,"temperature":1.0,"reasoning_tokens":2363,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:41:34.133802+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the NLO cross section and $\\lambda_\\theta$ with substantially different phase-space-slicing cutoffs, or with an independent subtraction scheme; if the results move beyond integration uncertainties, the reported enhancement and near-unpolarized polarization are artifacts of the slicing prescription.","supporting_citations":[{"cited_title":"Artoisenet, J","cited_arxiv_id":null,"evidence_quote":"Earlier leading-order study of J/psi+c cbar associated production whose finding of a sizable, nearly unpolarized gluon-fusion contribution motivates the full NLO calculation."},{"cited_title":"NLO corrections to $J/\\psi+c+\\bar{c}$ photoproduction","cited_arxiv_id":"2405.05683","evidence_quote":"Previous calculation by the authors whose computational strategy, helicity-amplitude handling, and loop treatment the present work largely follows."}],"review_version":2}