Pith. sign in

REVIEW 4 major objections 6 minor 78 references

Probing Heavy-Quark Spin Symmetry in Double $J/\psi$ Hadroproduction

T0 review · 4 major / 6 minor · reviewed 2026-08-05 · deepseek-v4-flash

Pith's one-line read This paper claims that LHC data on pairs of J/psi mesons provide the first experimental verification of heavy-quark spin symmetry in quarkonium, without putting that symmetry into the fit.

desk verdict First complete O(alpha_s^5) NRQCD calculation for double J/psi hadroproduction, with a clever LDME-combination argument and a novel HQSS test; the central result is solid but the unpropagated spread in the ratio r0 is a real weakness. read the letter →

arxiv 2608.03510 v1 pith:5DUEZV5M submitted 2026-08-04 hep-ph hep-ex

classification hep-phhep-ex
keywords heavyquarkoniumproductionNRQCDfactorizationcolor-octetmatrixelementsJ/psipairhadroproductionheavy-quarkspinsymmetryeta_cLHCphenomenologylong-distance
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper attempts to establish that double J/psi hadroproduction at the LHC can resolve a long-standing ambiguity in how quark-antiquark pairs become J/psi mesons, and can test a symmetry that theorists have often assumed. It reports the first complete calculation of prompt J/psi pair production at order alpha_s^5 in nonrelativistic QCD, including all 28 pairings of quark-antiquark Fock states. The key finding is that high-transverse-momentum CMS and ATLAS data do not fix the two dominant color-octet matrix elements separately, but only their linear combination M_{0,1.8}; combining this with two combinations already fixed by single-J/psi production yields a new set of matrix elements. That set agrees with world data on J/psi yield and polarization and, when spin symmetry is applied, with LHCb data on eta_c production. The paper reads this agreement as evidence that heavy-quark spin symmetry is real rather than an input.

What carries the argument

The load-bearing object is the linear combination M_{0,1.8}^{J/psi} = <O^{J/psi}(1S_0^[8])> + (1.8/m_c^2)<O^{J/psi}(3P_0^[8])>, defined in Eq. (3). It exists because the short-distance coefficient ratio R_{3S1[1],3PJ[8]}/R_{3S1[1],1S0[8]} stays close to 1.80 (observed range 1.77-1.86) over the 14 CMS/ATLAS bins at large pair pT. This lets two nonperturbative parameters collapse into the single combination the data can actually constrain; the companion combinations M_{0,3.9} and M_{1,-0.56} from single-J/psi production then break the remaining degeneracy. The full O(alpha_s^5) calculation of all 28 Fock-state pairings is the computational carrier that makes the comparison quantitative.

What would settle it

A reanalysis that lets <O(1S0^[8])> and <O(3P0^[8])> float separately in the CMS/ATLAS bins, or a measurement of d sigma/d pT^{psi psi} at lower pair pT where the 2->2 and 2->3 contributions differ, would show whether the 1.80 collapse is legitimate; if the extracted combination changes bin-by-bin or substantially improves the fit when split into two free parameters, the constancy assumption fails.

Watch

Extended reading notes

Core claim

The paper performs the first complete O(alpha_s^5) NRQCD computation of prompt J/psi pair hadroproduction, keeping every Fock-state combination of the two charm-anticharm pairs. Comparing with CMS and ATLAS data at large pair transverse momentum, it finds that the short-distance coefficients for the 3P_0^[8] and 1S_0^[8] color-octet channels have an almost constant ratio of about 1.80 across the 14 fitted bins. As a result, the data constrain only the combination <O^{J/psi}(1S_0^[8])> + (1.8/m_c^2)<O^{J/psi}(3P_0^[8])>, which the fit fixes at (3.49 +/- 0.59) x 10^-2 GeV^3. Combined with two LDME combinations from single-J/psi hadroproduction, this pins down the three color-octet matrix eleme

Load-bearing premise

The extraction treats the short-distance coefficient ratio of the 3P0^[8] to 1S0^[8] channels as fixed at 1.80 across all 14 fitted bins; if that ratio actually shifts with kinematics, scale choice, or higher orders, the merged matrix element and the spin-symmetry test would be biased.

Editorial extensions

If this is right

  • Large-pT double-J/psi production at the LHC is now described at order alpha_s^5 with all Fock states, closing the order-of-magnitude deficit of earlier color-singlet-only predictions through the color-octet mechanism.
  • The CMS and ATLAS data can be fitted with the color-singlet channel plus the 3S1^[8], 1S0^[8], and 3P0^[8] channels; adding all 21 retained channels changes the fit score by only about 5%.
  • The extracted LDME set reproduces the world prompt J/psi yield and polarization data and, under heavy-quark spin symmetry, the LHCb prompt eta_c data within uncertainties.
  • The lowest invariant-mass bins, which lie beyond leading-order kinematics, are also described, showing that the O(alpha_s^5) calculation covers real-radiation phase space.
  • Alternative LDME sets from the literature produce worse agreement on the same fit range, indicating that the pair-production data discriminate among existing extractions.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the constancy of the 1.80 ratio survives scale and kinematic variations, the same collapse of 1S0^[8] and 3P0^[8] matrix elements could be applied to other charmonium pair observables, such as J/psi + psi' or J/psi + chi_c, to sharpen LDME determinations.
  • A dedicated measurement at lower pT^{psi psi}, or with finer bins, could test the constancy assumption directly and might separate <O(1S0^[8])> from <O(3P0^[8])>, giving a more stringent check of the spin-symmetry relation than the current comparison.
  • Because the fit fixes only the linear combination, the agreement with eta_c data depends on applying HQSS; a future precise eta_c measurement that disagrees could force a reanalysis even if the pair-production fit remains good.
  • With the single-parton-scattering contribution computed at this order, the same framework could provide an independent extraction of the double-parton-scattering effective cross section after subtracting the SPS background.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. The paper presents the first complete O(α_s^5) NRQCD analysis of prompt J/ψ pair hadroproduction at large p_T^{ψψ}, including all allowed Fock-state pairings of the two c-cbar pairs. The authors compute the relevant short-distance coefficients with QGRAF/FORM/Cuba, classify the channels by power counting, and observe that CMS and ATLAS data constrain the linear combination M_{0,1.8}^{J/ψ}=<O^{J/ψ}(^1S_0^{[8]})> + (1.8/m_c^2)<O^{J/ψ}(^3P_0^{[8]})>. Combining this with two combinations extracted from single-J/ψ hadroproduction, they obtain a new CO LDME set. They then compare NLO NRQCD predictions for prompt η_c production with LHCb data, using HQSS to relate J/ψ and η_c LDMEs, and report good agreement (χ²/11 = 1.52 or 1.38 depending on the set). They also note that their O(α_s^5) predictions reproduce the lowest m^{ψψ} bins in CMS and ATLAS data, which are kinematically inaccessible at O(α_s^4). The central claim is that this constitutes the first genuine, data-driven verification of HQSS in heavy quarkonium production, because the J/ψ-pair fit itself does not assume HQSS.

Significance. If the extraction is robust, this is a valuable step. The computation covers a much larger set of Fock-state channels than previous double-J/ψ studies, and the observation that high-p_T^{ψψ} data pin down only a specific LDME combination is a useful and nontrivial finding. The paper also makes genuine predictions, particularly in the low-m^{ψψ} bins, and it provides an independent route to HQSS testing that does not build HQSS into the J/ψ fit. The agreement with the world J/ψ data and with LHCb η_c data is encouraging. However, the strength of the HQSS claim is currently limited by an unpropagated systematic uncertainty in the coefficient r0=1.8 and by a post-hoc choice of fit range. These issues are local and fixable, but they are load-bearing for the advertised conclusion.

major comments (4)
  1. [Section 'quantitative verification', Eq. (3)] The extraction of M_{0,1.8}^{J/ψ} rests on the near-constancy of the ratio R_{3S1[1],3PJ[8]}/R_{3S1[1],1S0[8]}, quoted as 1.77–1.86 with average 1.80 across the 14 fit bins. The spread is not propagated into the fit. This matters because the subsequent separation of <O(^1S0[8])> and <O(^3P0[8])> is a subtraction of two nearly equal combinations: with the central values, <O(^1S0[8])> = (3.9 M_{0,1.8} - 1.8 M_{0,3.9})/(3.9-1.8) ≈ 0.14×10^{-2} GeV^3. A change of r0 by ±0.06 shifts <O(^1S0[8])> by about ∓0.11×10^{-2} GeV^3, which is comparable to the difference between this fit and Ref. [24] that the paper itself identifies as responsible for the η_c χ² difference (1.52 vs 1.38). Please propagate the observed 1.77–1.86 range as a systematic uncertainty and show that the HQSS conclusion is stable under this variation.
  2. [Section 'quantitative verification', fit-range choice] The two ATLAS-I bins with p_T^{ψψ}>40 GeV are excluded because 'the data systematically undershoot the theoretical evaluations,' and the fit is then quoted with χ²/dof = 0.38. Later, including all bins gives χ²/dof = 0.40, and the excluded points are said to be 'very well described' except for two bins with small gaps. This is a post-hoc selection of the fit range based on the same data that are being fit. A principled criterion is needed, for example an a priori upper p_T^{ψψ} cutoff motivated by the validity of the power-counting or NRQCD factorization, or a systematic error obtained by repeating the fit with several cut choices (e.g., p_T>18, >20, >25, >30 GeV). Without this, the central value and uncertainty of M_{0,1.8}, and hence of the extracted LDME set, are not fully justified.
  3. [Section 'quantitative verification', η_c comparison] The claim of a 'first genuine verification of HQSS' is based on χ²/11 = 1.52 and 1.38 for the two LDME sets. These two values are very close, and the paper states that the small difference is due to <O(^1S0[8])>. This is not presented as a statistical test of HQSS, but as an affirmation of agreement. Please quantify the significance of the test: e.g., give the χ² for a hypothetical fit without HQSS, or the χ² obtained when the alternative LDME sets discussed in the paper are used for the η_c comparison. As written, the evidence for 'verification' is suggestive but not quantitatively calibrated.
  4. [Eq. (2) and Fig. 2] The paper's central quantitative results—the scaling ratios R_{m,n}, the ratio 1.77–1.86, and the fitted M_{0,1.8}—all rely on the multi-gigabyte numerical computation performed with QGRAF/FORM/Cuba, but no numerical tables or public code are provided. Since the SDC calculation is not independently checkable from the manuscript, the authors should provide at least the values of R_{m,n} for the 14 fit bins (and ideally a benchmark table for the CMS, ATLAS-I, ATLAS-II setups) in an appendix or supplementary material. This would also allow readers to reproduce the constancy of the crucial ratio.
minor comments (6)
  1. [Eq. (2)] The notation 'mdm+dn−6c' is unclear; it presumably means m_c^{d_m+d_n-6}. Please clarify.
  2. [Author affiliations] Typo: 'Mathemati cs' in the Beijing University affiliation.
  3. [References] Reference [26] appears twice in the bibliography list, with the second entry being a duplicate of the LHCb η_c paper.
  4. [Figure 3] The caption says 'dσ/p_T^{ψψ}' but the plots likely show dσ/dp_T^{ψψ}; please make the derivative explicit and specify the bin widths and units consistently.
  5. [Text near Eq. (3)] The combinations M_{0,3.9}^{J/ψ} and M_{1,-0.56}^{J/ψ} are taken from Refs. [11,17], but the definitions of r0=3.9 and r1=-0.56 are not restated. Please define them or refer explicitly to the equations in the original papers.
  6. [Fit statistics] The χ²/dof values are quoted as 0.38, 0.40, 1.31. Please give the number of data points, the number of fitted parameters, and the absolute χ² values so the reader can judge the goodness of fit.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the double-J/psi fit, the single-J/psi constraints, and the eta_c comparison are independent inputs used in a genuine prediction chain.

full rationale

The paper's derivation chain is self-contained and does not reduce to its inputs by construction. The double-J/psi SDCs are computed in this paper (with the CS channel taken from the authors' previous parameter-free calculation, Ref. [56]), and the near-constancy of the ratio R_{3S1[1],3P0[8]}/R_{3S1[1],1S0[8]} is a numerical observation from these SDCs, not an assumed output. The fitted quantity M_{0,1.8} is defined using this observed ratio and then constrained by CMS and ATLAS data; the individual LDMEs are obtained by combining this with two independent single-J/psi constraints from Refs. [11,17], which are external to the present paper. The eta_c comparison is a genuine prediction: LHCb eta_c data are not used in any step of the fit, and the HQSS test therefore has independent empirical content. The authors' reliance on their own earlier SDC calculations (Refs. [42,43,56]) is normal building-block citation of parameter-free, externally falsifiable results, not circular. The reader's concern about the unpropagated spread of the ratio r0 = 1.77–1.86 is a robustness/correctness issue, not a circularity, because the ratio is an input derived from the calculation and its uncertainty does not make the fitted combination logically equivalent to the data. No equation is defined in terms of the result it purports to derive, and no fitted parameter is relabeled as a prediction.

Assumptions & free parameters 6 free parameters · 6 assumptions · 0 invented entities

The central claim rests on standard NRQCD factorization, external LDME constraints from single-J/psi fits, and a few modeling choices (fixed r=1.80, psi(2S) feed-down factor, DPS handling, and the post-hoc exclusion of two ATLAS-I bins). No new beyond-SM entities are introduced.

free parameters (6)
  • M_{0,1.8}^{J/psi} = <O(^1S_0^{[8]})> + 1.8/m_c^2 <O(^3P_0^{[8]})> = (3.49 +/- 0.59) x 10^-2 GeV^3
    Fitted to 14 CMS and ATLAS bins with p_T^psi psi > 18/20 GeV; the only quantity directly constrained by the double-J/psi data.
  • M_{0,3.9}^{J/psi} = (7.4 +/- 1.9) x 10^-2 GeV^3
    External constraint from single prompt J/psi hadroproduction fits [11,17]; used as input to separate the LDMEs.
  • M_{1,-0.56}^{J/psi} = (0.05 +/- 0.02) x 10^-2 GeV^3
    External constraint from single prompt J/psi hadroproduction fits [11,17]; used as input.
  • r = coefficient in M_{0,1.8} = 1.80 (average of computed 1.77-1.86)
    Derived from the computed SDC ratios; fixed without propagating its spread, a possible source of bias.
  • Charm quark mass m_c = 1.5 GeV
    Standard input; enters kinematics and the LDME combinations; not varied.
  • psi(2S) feed-down factor = 7/6 per J/psi
    Approximation from known branching fractions, applied to each J/psi to account for psi(2S) feed-down.
assumptions (6)
  • domain assumption NRQCD factorization applies to double J/psi hadroproduction with universal LDMEs (Eq. (1)).
    The calculation and fit are entirely within NRQCD factorization; known factorization-breaking for double P-wave pairings [61] is mentioned but not resolved for the included channels.
  • domain assumption Heavy-quark spin symmetry (HQSS) relations connect J/psi and eta_c LDMEs.
    Used to convert the fitted J/psi LDME set into prompt eta_c yield predictions for comparison with LHCb data; the agreement is presented as a test of this symmetry.
  • domain assumption CTEQ6L1 LO PDFs and alpha_s with n_f=4, Lambda_QCD^(4)=215 MeV.
    Standard fixed inputs; only the scale xi is varied for uncertainty.
  • domain assumption DPS contamination is negligible (CMS) or already subtracted (ATLAS) in the selected large-p_T^psi psi bins.
    Relies on Ref [51] for CMS and on the ATLAS DPS extraction; if wrong, the fitted LDMEs are biased.
  • ad hoc to paper The two ATLAS-I bins with p_T^psi psi > 40 GeV are excluded from the fit because data undershoot theory.
    Post-hoc selection; changes chi^2/dof from 1.31 to 0.38 and affects M_{0,1.8}; no a priori physical reason for exclusion.
  • domain assumption Higher-order v^2 and alpha_s corrections beyond those included do not alter the fitted combination.
    The fit uses LO-in-v SDCs for the CO channels and NLO for the CS channel; missing v^2 corrections could shift the ratio r and the combination.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Probing Heavy-Quark Spin Symmetry in Double $J/\psi$ Hadroproduction." pith.science (2026). https://pith.science/paper/5DUEZV5M

@misc{pith2026260803510,
  author       = {Pith},
  title        = {Pith review of: Probing Heavy-Quark Spin Symmetry in Double $J/\psi$ Hadroproduction},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5DUEZV5M}},
  note         = {Machine review of arXiv:2608.03510}
}
abstract

We perform the first complete $\mathcal{O}(\alpha_s^5)$ analysis of the prompt hadroproduction of $J/\psi$ pairs with large transverse momenta $p_T^{\psi\psi}$ in the nonrelativistic-QCD factorization framework, including all possible Fock state combinations $c\bar{c}(m)+c\bar{c}(n)$, with $m,n={}^3S_1^{[1,8]},{}^1S_0^{[8]},{}^3P_J^{[1,8]}$. We observe that CMS and ATLAS data constrain a specific linear combination of the long-distance matrix elements (LDMEs) $\langle\mathcal{O}^{J/\psi}({}^1S_0^{[8]})\rangle$ and $\langle\mathcal{O}^{J/\psi}({}^3P_0^{[8]})\rangle$. In conjunction with two other combinations fixed by single prompt $J/\psi$ hadroproduction, we gain a new LDME set, which turns out to be largely compatible with the world data of prompt $J/\psi$ yield and polarization and to probe heavy-quark spin symmetry, by which agreement is established with LHCb data of prompt $\eta_c$ yield. Our $\mathcal{O}(\alpha_s^5)$ predictions also nicely agree with CMS and ATLAS data in the lowest bins of $J/\psi$ pair invariant mass $m^{\psi\psi}$, beyond leading-order kinematics.

Figures

Figures reproduced from arXiv: 2608.03510 by the authors.

Figure 1
Figure 1. FIG. 1: Typical Feynman diagrams for [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Ratios [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: CMS [39] (left), ATLAS-I (center), and ATLAS-II [40] [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

78 extracted references · 29 canonical work pages

  1. [8]

    Butenschoen and B

    M. Butenschoen and B. A. Kniehl, Next-to-Leading- Order Tests of Non-Relativistic-QCD Factorization with J/ψ Yield and Polarization, Mod. Phys. Lett. A 28, 1350027 (2013) [arXiv:1212.2037 [hep-ph]]

  2. [24]

    Butensch¨ on and B

    M. Butensch¨ on and B. A. Kniehl, Complete Next-to- Leading-Order Corrections to J/ψ Photoproduction in Nonrelativistic Quantum Chromodynamics, Phys. Rev. Lett. 104, 072001 (2010) [arXiv:0909.2798 [hep-ph]]

  3. [1]

    In single prompt J/ψ production, the feed-down contributions from χcJ andψ(2S) mesons are about 30% and 10%, respectively [51]

    = 0.107 GeV 5 corresponding to the evaluation of the wave function at the origin in the Buchm¨ uller-Tye poten- tial [68], and ⟨OχcJ (3S[8] 1 )⟩/(2J + 1) = 2.2 × 10−3 GeV5 [69]. In single prompt J/ψ production, the feed-down contributions from χcJ andψ(2S) mesons are about 30% and 10%, respectively [51]. This implies that the ψ(2S) feed-down contribution ...

  4. [2]

    W. E. Caswell and G. P. Lepage, Effective lagrangians for bound state problems in QED, QCD, and other field theories, Phys. Lett. B 167, 437–442 (1986)

  5. [3]

    and noticing that Br( χc1,2 → J/ψγ ) are numerically of O(v2), we can then roughly estimate the relative im- portance of each channel in the large- pψψ T limit; see Ta- ble I. In our numerical analysis, we take the on-shell mass to be mc = 1.5 GeV, adopt the LO CTEQ6L1 PDF set [67], and evaluate α(nf ) s (µr) with nf = 4 and asymptotic scale parameter Λ (...

  6. [4]

    data: (i) include only the 3S[1] 1 + (3S[1] 1 , 1S[8] 0 ) chan- nels and fit MJ/ψ 0,1.8 to the 14 bins with pψψ T > 18 GeV; (ii) combine this with MJ/ψ 0,3.9 andMJ/ψ 1,−0.56 from a fit to single prompt J/ψ hadroproduction [11, 17], including LP and NLP channels [71, 72], and so pin down ⟨OJ/ψ (1S[8] 0 )⟩, ⟨OJ/ψ (3S[8] 1 )⟩, and ⟨OJ/ψ (3P [8] 0 )⟩ separately...

  7. [5]

    Brambilla et al

    N. Brambilla et al. (Quarkonium Working Group), Heavy Quarkonium Physics, CERN Yellow Report No. CERN- 2005-005 [arXiv:hep-ph/0412158 [hep-ph]]

  8. [6]

    and ATLAS [40] data at large pψψ T provides a first genuine verification of HQSS in heavy quarkonium. Com- paring NLO predictions for prompt ηc yield [21, 24] with LHCb data at √ s = 7 , 8, 13 TeV [26, 27], we obtain χ2/11 = {1.52, 1.38} for the two LDME sets in Table II, the small difference being due to ⟨OJ/ψ (1S[8] 0 )⟩. In summary, we performed the first ...

Show all 78 references
  1. [7]

    G. T. Bodwin, E. Braaten, and G. P. Lepage, Rigorous QCD analysis of inclusive annihilation and production of heavy quarkonium, Phys. Rev. D 51, 1125–1171 (1995); 55, 5853(E) (1997) [arXiv:hep-ph/9407339 [hep-ph]]

  2. [9]

    G. P. Lepage, L. Magnea, C. Nakhleh, U. Magnea, and K. Hornbostel, Improved nonrelativistic QCD for heavy-quark physics, Phys. Rev. D 46, 4052–4067 (1992) [arXiv:hep-lat/9205007 [hep-lat]]

  3. [10]

    Butenschoen and B

    M. Butenschoen and B. A. Kniehl, Constraints on Nonrelativistic-QCD Long-Distance Matrix Elements fromJ/ψ PlusW/Z Production at the LHC, Phys. Rev. Lett. 130, 041901 (2023) [arXiv:2207.09366 [hep-ph]]

  4. [11]

    Y.-Q. Ma, K. Wang, and K.-T. Chao, J/ψ (ψ′) Pro- duction at the Tevatron and LHC at O(α4 sv4) in Non- relativistic QCD, Phys. Rev. Lett. 106, 042002 (2011) [arXiv:1009.3655 [hep-ph]]

  5. [12]

    Lansberg, J/ψ, ψ′ and Υ Production at Hadron Colliders: A Review, Int

    J.-P. Lansberg, J/ψ, ψ′ and Υ Production at Hadron Colliders: A Review, Int. J. Mod. Phys. A 21, 3857–3916 (2006) [arXiv:hep-ph/0602091 [hep-ph]]

  6. [13]

    Brambilla et al

    N. Brambilla et al. , Heavy Quarkonium: progress, puz- zles, and opportunities, Eur. Phys. J. C 71, 1534 (2011) [arXiv:1010.5827 [hep-ph]]

  7. [14]

    Chao, Y.-Q

    K.-T. Chao, Y.-Q. Ma, H.-S. Shao, K. Wang, and Y.- J. Zhang, J/ψ Polarization at Hadron Colliders in Non- relativistic QCD, Phys. Rev. Lett. 108, 242004 (2012) [arXiv:1201.2675 [hep-ph]]

  8. [15]

    Brambilla et al

    N. Brambilla et al. , QCD and strongly coupled gauge theories: challenges and perspectives, Eur. Phys. J. C 74, 2981 (2014) [arXiv:1404.3723 [hep-ph]]

  9. [16]

    Lansberg, New observables in inclusive produc- tion of quarkonia, Phys

    J.-P. Lansberg, New observables in inclusive produc- tion of quarkonia, Phys. Rept. 889, 1–106 (2020) [arXiv:1903.09185 [hep-ph]]

  10. [17]

    H.-S. Shao, H. Han, Y.-Q. Ma, C. Meng, Y.-J. Zhang, and K.-T. Chao, Yields and polarizations of prompt J/ψ and ψ(2S) production in hadronic collisions, JHEP 05, 103 (2015) [arXiv:1411.3300 [hep-ph]]

  11. [18]

    Butensch¨ on and B

    M. Butensch¨ on and B. A. Kniehl, Reconciling J/ψ Pro- duction at HERA, RHIC, Tevatron, and LHC with Non- relativistic QCD Factorization at Next-to-Leading Order, Phys. Rev. Lett. 106, 022003 (2011) [arXiv:1009.5662 [hep-ph]]

  12. [19]

    Butenschoen and B

    M. Butenschoen and B. A. Kniehl, World data of J/ψ production consolidate nonrelatvistic QCD factorization at next-to-leading order, Phys. Rev. D 84, 051501(R) (2011) [arXiv:1105.0820 [hep-ph]]

  13. [20]

    Butenschoen and B

    M. Butenschoen and B. A. Kniehl, J/ψ Polarization at the Tevatron and the LHC: Nonrelativistic-QCD Factor- ization at the Crossroads, Phys. Rev. Lett. 108, 172002 (2012) [arXiv:1201.1872 [hep-ph]]

  14. [21]

    Gong, L.-P

    B. Gong, L.-P. Wan, J.-X. Wang, and H.-F. Zhang, Po- larization for Prompt J/ψ and ψ(2s) Production at the Tevatron and LHC, Phys. Rev. Lett. 110, 042002 (2013) [arXiv:1205.6682 [hep-ph]]

  15. [22]

    G. T. Bodwin, H. S. Chung, U-R. Kim, and J. Lee, Frag- mentation Contributions to J/ψ Production at the Teva- tron and the LHC, Phys. Rev. Lett. 113, 022001 (2014) 6 [arXiv:1403.3612 [hep-ph]]

  16. [23]

    Butenschoen and B

    M. Butenschoen and B. A. Kniehl, Dipole subtraction vs. phase space slicing in NLO NRQCD heavy-quarkonium production calculations, Nucl. Phys. B 957, 115056 (2020) [arXiv:2003.01014 [hep-ph]]

  17. [25]

    Butenschoen and B

    M. Butenschoen and B. A. Kniehl, Probing Nonrelativis- tic QCD Factorization in Polarized J/ψ Photoproduction at Next-to-Leading Order, Phys. Rev. Lett. 107, 232001 (2011) [arXiv:1109.1476 [hep-ph]]

  18. [26]

    Aaij et al

    R. Aaij et al. (LHCb Collaboration), Measurement of the ηc(1S) production cross-section in proton–proton colli- sions via the decay ηc(1S) →p¯p, Eur. Phys. J. C 75, 311 (2015) [arXiv:1409.3612 [hep-ex]]

  19. [27]

    Butenschoen, Z.-G

    M. Butenschoen, Z.-G. He, and B. A. Kniehl, ηc Pro- duction at the LHC Challenges Nonrelativistic QCD Factorization, Phys. Rev. Lett. 114, 092004 (2015) [arXiv:1411.5287 [hep-ph]]

  20. [28]

    Butenschoen and B

    M. Butenschoen and B. A. Kniehl, Dipole subtraction at next-to-leading order in nonrelativistic-QCD factoriz a- tion, Nucl. Phys. B 950, 114843 (2020) [arXiv:1909.03698 [hep-ph]]

  21. [29]

    Zhang, Z

    H.-F. Zhang, Z. Sun, W.-L. Sang, and R. Li, Impact of ηc Hadroproduction Data on Charmonium Production and Polarization within the Nonrelativistic QCD Framework, Phys. Rev. Lett. 114, 092006 (2015) [arXiv:1412.0508 [hep-ph]]

  22. [30]

    Brambilla, M

    N. Brambilla, M. Butenschoen, and X.-P. Wang, How well does nonrelativistic QCD factorization work at next- to-leading order?, Phys. Rev. D 112, L011902 (2025) [arXiv:2411.16384 [hep-ph]]

  23. [31]

    Khachatryan et al

    V. Khachatryan et al. (CMS Collaboration), Measure- ment of J/ψ and ψ(2S) Prompt Double-Differential Cross Sections in pp Collisions at √ s = 7 TeV, Phys. Rev. Lett. 114 (2015), 191802 [arXiv:1502.04155 [hep- ex]]. [26]

  24. [32]

    Li, Y.-J

    R. Li, Y.-J. Zhang. and K.-T. Chao, Pair production of heavy quarkonium and B(∗) c mesons at hadron colliders, Phys. Rev. D 80, 014020 (2009) [arXiv:0903.2250 [hep- ph]]

  25. [33]

    Aaij et al

    R. Aaij et al. (LHCb Collaboration), Measurement of the ηc(1S) production cross-section in pp collisions at √ s = 13 TeV, Eur. Phys. J. C 80, 191 (2020) [arXiv:1911.03326 [hep-ex]]

  26. [34]

    Han, Y.-Q

    H. Han, Y.-Q. Ma, C. Meng, H.-S. Shao, and K.- T. Chao, ηc Production at LHC and Indications on the Understanding of J/ψ Production, Phys. Rev. Lett. 114, 092005 (2015) [arXiv:1411.7350 [hep-ph]]

  27. [35]

    Aaij et al

    R. Aaij et al. (LHCb Collaboration), Observation of structure in the J/ψ-pair mass spectrum, Sci. Bull. 65, 1983–1993 (2020) [arXiv:2006.16957 [hep-ex]]

  28. [36]

    Barger, S

    V. Barger, S. Fleming, and R. J. N. Phillips, Double gluon fragmentation to J/ψ pairs at the Tevatron, Phys. Lett. B 371, 111–116 (1996) [arXiv:hep-ph/9510457 [hep-ph]]

  29. [37]

    Qiao, J/ψ pair production at the Fermilab Teva- tron, Phys

    C.-F. Qiao, J/ψ pair production at the Fermilab Teva- tron, Phys. Rev. D 66, 057504 (2002) [arXiv:hep- ph/0206093 [hep-ph]]

  30. [38]

    Aaij et al

    R. Aaij et al. (LHCb Collaboration), Measurement of J/ψ-pair production in pp collisions at √s = 13 TeV and study of gluon transverse-momentum dependent PDFs, JHEP 03, 088 (2024) [arXiv:2311.14085 [hep-ex]]

  31. [39]

    and ATLAS [40] data at large pψψ T may be gained by consideringMJ/ψ 0,1.8. The most frequently used alternative sets [13, 15, 16] have MJ/ψ 0,1.8 = {2.3±0.37, 8.0±1.0, 10.8± 2.3} × 10−2 GeV3 and thus yield central predictions that correspondingly undershoot or overshoot the da...

  32. [40]

    cases, however, the NLO predictions turned out to undershoot the total cross sections by about one order of magnitude, despite sizable QCD correction factors of about 3, and to poorly describe the various cross section distributions [56]. The NLO predictions fall particularly ...

  33. [41]

    P. Ko, J. Lee, and C. Yu, Inclusive double-quarkonium production at the Large Hadron Collider, JHEP 01, 070 (2011) [arXiv:1007.3095 [hep-ph]]

  34. [42]

    C. H. Kom, A. Kulesza, and W. J. Stirling, Pair Pro- duction of J/ψ as a Probe of Double Parton Scat- tering at LHCb, Phys. Rev. Lett. 107, 082002 (2011) [arXiv:1105.4186 [hep-ph]]

  35. [43]

    Aaij et al

    R. Aaij et al. (LHCb Collaboration), Observation of J/ψ- pair production in pp collisions at √ s = 7 TeV, Phys. Lett. B 707, 52–59 (2012) [arXiv:1109.0963 [hep-ex]]

  36. [44]

    Aaij et al

    R. Aaij et al. (LHCb Collaboration), Measurement of the J/ψ pair production cross-section in pp collisions at√s = 13 TeV, JHEP 06, 047 (2017); 10, 068(E) (2017) [arXiv:1612.07451 [hep-ex]]

  37. [45]

    Khachatryan et al

    V. Khachatryan et al. (CMS Collaboration), Measure- ment of prompt J/ψ pair production in pp collisions at√ s = 7 Tev, JHEP 09, 094 (2014) [arXiv:1406.0484 [hep- ex]]

  38. [46]

    Aaboud et al

    M. Aaboud et al. (ATLAS Collaboration), Measurement of the prompt J/ψ pair production cross-section in pp collisions at √s = 8 TeV with the ATLAS detector, Eur. Phys. J. C 77, 76 (2017) [arXiv:1612.02950 [hep-ex]]

  39. [47]

    Acharya et al

    S. Acharya et al. (ALICE Collaboration), Measurement of inclusive J/ψ pair production cross section in pp colli- sions at √s = 13 TeV, Phys. Rev. C 108, 045203 (2023) [arXiv:2303.13431 [hep-ex]]

  40. [48]

    He and B

    Z.-G. He and B. A. Kniehl, Complete Nonrelativistic- QCD Prediction for Prompt Double J/ψ Hadropro- duction, Phys. Rev. Lett. 115, 022002 (2015) [arXiv:1609.02786 [hep-ph]]

  41. [49]

    Z.-G. He, B. A. Kniehl, M. A. Nefedov, and V. A. Saleev, Double Prompt J/ψ Hadroproduction in the Parton Reggeization Approach with High-Energy Resummation, Phys. Rev. Lett. 123, 162002 (2019) [arXiv:1906.08979 [hep-ph]]

  42. [50]

    B. A. Kniehl, M. A. Nefedov, and V. A. Saleev, Prompt- photon plus jet associated photoproduction at HERA in the parton Reggeization approach, Phys. Rev. D 89, 114016 (2014) [arXiv:1404.3513 [hep-ph]]

  43. [51]

    A. V. Karpishkov, M. A. Nefedov, and V. A. Saleev, B ¯B angular correlations at the LHC in parton Reggeiza- tion approach merged with higher-order matrix elements, Phys. Rev. D 96, 096019 (2017) [arXiv:1707.04068 [hep- ph]]

  44. [52]

    E. A. Kuraev, L. N. Lipatov, and V. S. Fadin, Multi- regge processes in the Yang–Mills theory, Zh. Eksp. Teor. Fiz. 71, 840–855 (1976) [Sov. Phys. JETP 44, 443–451 (1976)]

  45. [53]

    Ya. Ya. Balitski ˘i and L. N. Lipatov, The Pomeranchuk singularity in quantum chromodynamics, Yad. Fiz. 28, 1597–1611 (1978) [Sov. J. Nucl. Phys. 28, 822–829 (1978)]

  46. [54]

    S. P. Baranov, A. M. Snigirev, N. P. Zotov, A. Szczurek, 7 and W. Sch¨ afer, Interparticle correlations in the produc- tion of J/ψ pairs in proton-proton collisions, Phys. Rev. D 87, 034035 (2013) [arXiv:1210.1806 [hep-ph]]

  47. [55]

    Lansberg and H

    J.-P. Lansberg and H. S. Shao, Production of J/ψ +ηc versus J/ψ + J/ψ at the LHC: Importance of Real α5 s Corrections, Phys. Rev. Lett. 111, 122001 (2013) [arXiv:1308.0474 [hep-ph]]

  48. [56]

    L.-P. Sun, H. Han, and K.-T. Chao, Impact of J/ψ pair production at the LHC and predictions in nonrelativistic QCD, Phys. Rev. D 94, 074033 (2016) [arXiv:1404.4042 [hep-ph]]

  49. [57]

    Lansberg and H.-S

    J.-P. Lansberg and H.-S. Shao, J/ψ-pair production at large momenta: Indications for double parton scatterings and large α5 s contributions, Phys. Lett. B 751, 479–486 (2015) [arXiv:1410.8822 [hep-ph]]

  50. [58]

    Lansberg and H

    J.-P. Lansberg and H. S. Shao, Double-quarkonium production at a fixed-target experiment at the LHC (AFTER@LHC), Nucl. Phys. B 900, 273–294 (2015) [arXiv:1504.06531 [hep-ph]]

  51. [59]

    Lansberg, H.-S

    J.-P. Lansberg, H.-S. Shao, N. Yamanaka, and Y.- J. Zhang, Prompt J/ψ-pair production at the LHC: im- pact of loop-induced contributions and of the colour- octet mechanism,” Eur. Phys. J. C 79, 1006 (2019) [arXiv:1906.10049 [hep-ph]]

  52. [60]

    Sun, J/ψ pair hadroproduction at next-to-leading order in nonrelativistic-QCD at CMS, Chin

    L.-P. Sun, J/ψ pair hadroproduction at next-to-leading order in nonrelativistic-QCD at CMS, Chin. Phys. C 47, 093105 (2023) [arXiv:2307.02809 [hep-ph]]

  53. [61]

    Sun, J/ψ pair hadroproduction at next-to-leading order in nonrelativistic-QCD at ATLAS, Chin

    L.-P. Sun, J/ψ pair hadroproduction at next-to-leading order in nonrelativistic-QCD at ATLAS, Chin. Phys. C 48, 033102 (2024)

  54. [62]

    He, X.-B

    Z.-G. He, X.-B. Jin, and B. A. Kniehl, Next-to- leading-order QCD corrections to double prompt J/ψ hadroproduction, Phys. Rev. D 111, 094040 (2025) [arXiv:2505.04357 [hep-ph]]

  55. [63]

    Li, G.-Z

    Y.-J. Li, G.-Z. Xu, K.-Y. Liu, and Y.-J. Zhang, Rela- tivistic correction to J/ψ and Υ pair production, JHEP 07, 051 (2013) [arXiv:1303.1383 [hep-ph]]

  56. [64]

    He, X.-B

    Z.-G. He, X.-B. Jin, and B. A. Kniehl, Relativistic cor- rections to prompt double charmonium hadroproduc- tion near threshold, Phys. Rev. D 109, 094013 (2024) [arXiv:2402.07773 [hep-ph]]

  57. [65]

    Z.-G. He, B. A. Kniehl, M. A. Nefedov, and V. A. Saleev, Double promptJ/ψ production at hadron colliders, Mod. Phys. Lett. A 36, 2130018 (2021)

  58. [66]

    V. M. Abazov et al. (D0 Collaboration), Observation and studies of double J/ψ production at the Tevatron, Phys. Rev. D 90, 111101 (2014) [arXiv:1406.2380 [hep-ex]]

  59. [67]

    Z.-G. He, B. A. Kniehl, and X.-P. Wang, Breakdown of Nonrelativistic QCD Factorization in Processes Involv- ing Two Quarkonia and its Cure, Phys. Rev. Lett. 121, 172001 (2018) [arXiv:1809.07993 [hep-ph]]

  60. [68]

    Nogueira, Automatic Feynman Graph Generation, J

    P. Nogueira, Automatic Feynman Graph Generation, J. Comput. Phys. 105, 279–289 (1993)

  61. [69]

    J. A. M. Vermaseren, New features of FORM, [arXiv:math-ph/0010025 [math-ph]]

  62. [70]

    Hahn, CUBA–a library for multidimensional numer- ical integration, Comput

    T. Hahn, CUBA–a library for multidimensional numer- ical integration, Comput. Phys. Commun. 168, 78–95 (2005) [arXiv:hep-ph/0404043 [hep-ph]]

  63. [71]

    Navas et al

    S. Navas et al. (Particle Data Group), Review of Particle Physics, Phys. Rev. D 110, 030001 (2024)

  64. [72]

    Kang, J.-W

    Z.-B. Kang, J.-W. Qiu, and G. Sterman, Heavy Quarko- nium Production and Polarization, Phys. Rev. Lett. 108, 102002 (2012) [arXiv:1109.1520 [hep-ph]]

  65. [73]

    Pumplin, D

    J. Pumplin, D. R. Stump, J. Huston, H.-L. Lai, P. Nadol- sky, and W.-K. Tung, New generation of parton distribu- tions with uncertainties from global QCD analysis, JHEP 07, 012 (2002) [arXiv:hep-ph/0201195 [hep-ph]]

  66. [74]

    E. J. Eichten and C. Quigg, Quarkonium wave func- tions at the origin, Phys. Rev. D 52, 1726–1728 (1995) [arXiv:hep-ph/9503356 [hep-ph]]

  67. [75]

    Y.-Q. Ma, K. Wang, and K.-T. Chao, QCD radiative corrections to χcJ production at hadron colliders, Phys. Rev. D 83, 111503(R) (2011) [arXiv:1002.3987 [hep-ph]]

  68. [76]

    Aaij et al

    R. Aaij et al. (LHCb Collaboration), Measurement of associated J/ψ-ψ(2S) production cross-section in pp collisions at √ s = 13 TeV, JHEP 05, 259 (2024) [arXiv:2311.15921 [hep-ex]]

  69. [77]

    Kang, Y.-Q

    Z.-B. Kang, Y.-Q. Ma, J.-W. Qiu, and G. Sterman, Heavy quarkonium production at collider energies: Fac- torization and evolution, Phys. Rev. D 90, 034006 (2014) [arXiv:1401.0923 [hep-ph]]

  70. [78]

    Ma, J.-W

    Y.-Q. Ma, J.-W. Qiu, G. Sterman, and H. Zhang, Fac- torized Power Expansion for High- pT Heavy Quarko- nium Production, Phys. Rev. Lett. 113, 142002 (2014) [arXiv:1407.0383 [hep-ph]]

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.