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

=== DRAFT ===

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

Problem Statement
Despite over half a century of research, the production mechanism of heavy quarkonium remains partially mysterious. While the Non-Relativistic QCD (NRQCD) factorization formalism provides a rigorous framework separating perturbative short-distance coefficients (SDCs) from non-perturbative long-distance matrix elements (LDMEs), the universality of LDMEs across different production processes is still under intense scrutiny. Previous fits to single prompt $J/\psi$ production data have yielded LDME sets that generally agree with world data, with notable exceptions in specific kinematic regions (e.g., low/medium $p_T$, inelastic photoproduction). Furthermore, previous analyses often relied on the assumption of Heavy-Quark Spin Symmetry (HQSS) to relate $J/\psi$ and $\eta_c$ LDMEs. A significant gap exists in the theoretical description of prompt $J/\psi$ pair hadroproduction at large transverse momenta ($p_T^{\psi\psi}$). Leading-order (LO) and Next-to-Leading Order (NLO) calculations in the Color Singlet (CS) channel have been shown to undershoot experimental data from CMS and ATLAS by an order of magnitude, a discrepancy that cannot be resolved by Double Parton Scattering (DPS) alone.

Methodology
The authors perform the first complete $O(\alpha_s^5)$ analysis for the prompt hadroproduction of $J/\psi$ pairs with large transverse momenta within the NRQCD factorization framework.

• Scope of Channels: The calculation includes all possible Fock state combinations $c\bar{c}(m) + c\bar{c}(n)$ where $m, n \in \{ {}^3S_1^{[1,8]}, {}^1S_0^{[8]}, {}^3P_J^{[1,8]} \}$. This encompasses 28 possible pairings, reduced to 21 active channels after suppressing the ${}^3P_0^{[1]}$ channel due to the small branching fraction of $\chi_{c0} \to J/\psi\gamma$.
• Computational Framework:
  • SDC Calculation: Using QGRAF for diagram generation and FORM for Dirac and SU(3) color algebra, the authors computed the SDCs. Special attention was paid to double P-wave channels, which required sequential differentiation of scattering amplitudes with respect to relative momenta, resulting in complex analytical expressions.
  • Numerical Integration: A quadruple-precision C++ program utilizing the Cuba library was developed to achieve 1% numerical precision.
  • Power Counting: The channels were classified into Leading Power (LP, $N=4$), Next-to-Leading Power (NLP, $N=6$), and Next-to-Next-to-Leading Power (NNLP, $N=8$) based on their scaling behavior with $p_T^{\psi\psi}$ and the perturbative QCD power counting rules.
• Experimental Data: The analysis targets high-$p_T^{\psi\psi}$ bins from CMS ($\sqrt{s}=7$ TeV) and ATLAS ($\sqrt{s}=8$ TeV) experiments, specifically focusing on regions where $p_T^{\psi\psi} > 18$ GeV (CMS) and $> 20$ GeV (ATLAS), where Color Octet (CO) mechanisms are expected to dominate.
• Fit Strategy:
  1. Fit the linear combination $M_{0,1.8}^{J/\psi} = \langle O_{J/\psi}(^1S_0^{[8]}) \rangle + 1.8 \frac{1}{m_c^2} \langle O_{J/\psi}(^3P_0^{[8]}) \rangle$ to the 14 high-$p_T^{\psi\psi}$ bins where the $3S_1^{[1]} + (1S_0^{[8]}, 3P_J^{[8]})$ channels dominate.
  2. Combine this result with two other linear combinations ($M_{0,3.9}^{J/\psi}$ and $M_{1,-0.56}^{J/\psi}$) previously extracted from single prompt $J/\psi$ hadroproduction fits.
  3. Solve for the individual LDMEs: $\langle O_{J/\psi}(^1S_0^{[8]}) \rangle$, $\langle O_{J/\psi}(^3S_1^{[8]}) \rangle$, and $\langle O_{J/\psi}(^3P_0^{[8]}) \rangle$.
  4. Verify stability by including all 21 channels in the fit.

Key Contributions and Results

• Complete $O(\alpha_s^5)$ Analysis: The paper provides the first complete calculation of the $O(\alpha_s^5)$ contributions for all relevant Fock state combinations in double $J/\psi$ production, moving beyond the limited CS-only or LO studies.
• Constraint on LDMEs: The analysis reveals that CMS and ATLAS data at large $p_T^{\psi\psi}$ specifically constrain the linear combination $M_{0,1.8}^{J/\psi}$.
• New LDME Set: By combining the new constraint with existing single-production data, the authors derive a new set of CO LDMEs (Table II). This set yields a $\chi^2/\text{dof} \approx 0.40$ when applied to the full dataset, indicating excellent agreement.
• Agreement with Data:
  • The $O(\alpha_s^5)$ predictions successfully describe the CMS and ATLAS data in the large $p_T^{\psi\psi}$ region, resolving the order-of-magnitude deficit seen in previous NLO CS calculations.
  • The predictions also agree well with data in the lowest bins of the $J/\psi$ pair invariant mass ($m_{\psi\psi}$), extending beyond the strict kinematics of the leading-order calculation.
• Testing HQSS: Crucially, the new LDME set was derived without assuming Heavy-Quark Spin Symmetry (HQSS). When this set is used to predict the prompt $\eta_c$ yield (assuming HQSS relations between $J/\psi$ and $\eta_c$ LDMEs), the results agree well with LHCb data ($\chi^2/11 \approx 1.38-1.52$).

Significance
The paper claims a first genuine experimental verification of Heavy-Quark Spin Symmetry in heavy-quarkonium physics. Unlike previous works (e.g., Ref. [24]) that assumed HQSS to relate $J/\psi$ and $\eta_c$ LDMEs to achieve a global fit, this work derives the $J/\psi$ LDMEs purely from $J/\psi$ data (both single and pair production). The subsequent agreement of these data-driven LDMEs with $\eta_c$ yield data serves as a robust, independent test of HQSS. Additionally, the work establishes that the Color Octet mechanism is essential for describing high-$p_T$ double quarkonium production and provides a more reliable basis for extracting DPS parameters ($\sigma_{\text{eff}}$) and understanding backgrounds for exotic tetraquark searches.