NRQCD colour-octet expansion vs. LHC quarkonium production: signs of a hierarchy puzzle?

The observation of unpolarized quarkonium production in high energy pp collisions, at mid rapidity, implies a significant violation of the non-relativistic QCD (NRQCD) velocity scaling rules. A precise experimental confirmation of this picture could definitely rule out the current formulation of the factorization expansion. This conclusion relies on current perturbative determinations of the short-distance kinematic factors and may be reverted if improved calculations would modify, in a very specific way, their transverse momentum dependences. That solution would result, however, in a full degeneracy in the presently assumed basis of 2S+1LJ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{2S+1}L_J$$\end{document} Fock states. Therefore, whatever the outcome, improved polarization measurements will challenge and improve our fundamental understanding of quarkonium production.

where M, p T , and y are the mass, transverse momentum, and rapidity of the Q Q.
Existing NLO SDC calculations can quantitatively describe the observed "universal unpolarized" scenario of LHC mid-rapidity data, where all states are seemingly produced with identical p T /M distributions and zero polarizations [24]. From a conceptual point of view, however, the structure of the NRQCD factorization expansion, where three kinematically very different octet terms ( 1 S [8] 0 , 3 S [8] 1 , and 3 P [8] J ) dominate the production of S-wave states and two singlet terms ( 3 P [1] 1 and 3 P [1] 2 ) are additionally necessary (summed to the 3 S [8] 1 term) to describe χ c1,b1 and χ c2,b2 production, does not seem to naturally match the exceptionally simple patterns seen in the experimental data. In fact, the theory accomplishes its formal success thanks to unexpected, precise cancellations.
While the seemingly superfluous complexity of the theory formulation is particularly apparent in the comparison with the remarkably simple experimental patterns, a certain degree of theory redundancy can be seen independently of any experimental observation. For S-wave quarkonium production, which is the focus of this paper, the dominating terms of the NRQCD expansion are the 1 S [8] 0 , 3    1 singlet term being negligible because of its small SDC). Throughout this paper we use the next-to-leading order (NLO) calculations of Refs. [17,18] for the SDCs as functions of p T , in pp collisions at √ s = 7 TeV and mid rapidity. Figure 1 shows that the shape of the 1 S [8] 0 SDC is indistinguishable from a linear combination of the two other terms, 3 S [8] 1 + κ 3 P [8] J , with κ between 1.8 and 1.85. While in principle one would expect to fit experimental data with a superposition of three independent terms, leading to a definite determination of the three corresponding LDMEs, it turns out that the presently-available NLO SDCs are not independent kinematic templates and even very precise measurements of p Tdifferential cross sections will only be able to determine, in the best case, two parameters (as already noted in Ref. [15]). In other words, the NLO description of the p T dependence of quarkonium production has three theory parameters, but only two degrees of freedom. This observation suggests the existence of a spurious element of complexity in the assumed base of subprocesses.
According to the present NLO knowledge of the SDCs, the polarization is a better discriminating observable. Figure 2 compares the calculated dilepton-decay polar anisotropy parameter, λ ϑ , as a function of p T /M, for the 1 S [8] 0 , 3 S [8] 1 , and 3 S [8] 1 + κ 3 P [8] J terms. Contrary to the case of the production yields, no value of κ exists for which λ ϑ ( 3 S [8] 1 + κ 3 P [8] J ) ≈ λ ϑ ( 1 S [8] 0 ) as a function of p T /M, and a momentumdependent polarization measurement should be able to dis- entangle the three components. However, for p T /M 10, the 3 S [8] 1 + κ 3 P [8] J combination with κ in the 1.8-1.85 range becomes unpolarized, just as the 1 S [8] 0 term. Therefore, at high p T , the polarization measurements are also unable to resolve more than two degrees of freedom in the space of the contributing LDMEs, at least for quarkonia produced at mid-rapidity in high-energy pp collisions.
It is worth noting that for κ 1.85 the 3 S J combination assumes unphysical behaviours, represented in Figs. 1 and 2 by the dotted lines: the p T /M distribution shows seemingly anomalous changes of curvature at high p T /M while the polarization parameter starts decreasing towards unphysical asymptotic values (λ ϑ < −1). Intriguingly, the "degeneracy condition" κ = 1.8-1.85 happens just before the border of the physical domain of positivity.
We will now see what the experimental measurements can add to the picture. Figures 3 and 4 show the mid-rapidity LHC measurements of cross sections and polarizations of different charmonium and bottomonium states.
As already mentioned, no significant differences are observed in the p T /M dependences of yields and polarizations for the different states, despite their varying feeddown contributions (ranging between 0 and 40%) from heavier quarkonia. This universal behaviour has been discussed in detail in Ref. [24], where it is shown that the χ c feeddown contributions do not alter significantly, within the current experimental precision, the kinematic patterns of the directly-produced J/ψ mesons. It is, therefore, reasonable to compare the inclusive data (including feed-down) to the Mid-rapidity quarkonium cross sections measured in pp collisions at √ s = 7 TeV by ATLAS (red markers) [3][4][5] and CMS (blue markers) [6,7], with normalizations arbitrarily adjusted to the J/ψ points to illustrate the universality of the p T /M dependence. The data are compared to the shapes of the 1 S [8] 0 octet (green) and of the 3 NRQCD curves computed for direct production (excluding feed-down). For the curves themselves we use the NLO SDC calculations [17,18] performed for a Q Q mass of 3 GeV, assuming the validity of the p T /M scaling observed in cross section data as a general rule for individual processes. This approximation is justified by dimensional analysis reasonings applied to a single production mechanism [25], when √ s is large with respect to the momentum and mass of the observed state.
The figure shows that the scenario κ = 1.8-1.85 is not just an hypothetical idea to be considered among many others but a scenario in rather good agreement with the presently existing quarkonium cross section and polarization measurements. For p T /M 10, the cross section measurements are compatible with the p T /M distribution of the 1 S [8] 0 term, only a small correction coming from other terms. The observed lack of polarization further supports the idea that the 1 S [8] 0 term is close to reproducing the data with no need of further contributions. For high p T /M, the NRQCD expansion appears, therefore, as fully degenerate: all the data can be described by a single one of the basis terms, the 1 S [8] 0 octet, or, equivalently, by the combination 3 S [8] 1 + κ 3 P [8] J , with κ = 1.8-1.85.
As previously mentioned, the degeneracy is broken for p T /M < 10, where the NLO SDCs lead to different polarization predictions. However, the measured polarizations remain independent of p T /M down to p T /M ∼ 1. This means that the comparison between the measured data and the NLO SDCs clearly excludes a significant contribution of the 3 S [8] 1 + κ 3 P [8] J term, which is expected to be either strongly p T dependent or strongly polarized, as shown in Fig. 2. Even a relatively small 3 S [8] 1 + κ 3 P [8] J contribution to the quarkonium production yields, of around 10%, is already excluded by the presently existing polarization data, independently of the value of κ. In other words, the mid-rapidity quarkonium polarization measurements made at the LHC indicate a strong hierarchy of cross section contributions, clearly enhancing the 1 S [8] 0 octet and suppressing the individual 3 S [8] 1 and 3 P [8] J terms, not only their partially cancelling combination. This observation constitutes a significant violation of the v-scaling rules, according to which the three LDMEs should have the same order of magnitude.
Such a strong and unexpected constraint relies on the accuracy of the present NLO calculations. Uncertainties from the factorization scale and from the quark masses, mainly affecting the SDC normalizations, have a negligible impact on our considerations, which are entirely based on shape differences. Since the shapes of the NLO 1 S [8] 0 and 3 S [8] 1 SDCs are similar to the leading order (LO) versions, it is reasonable to expect that future computations of higher-order contributions will not reveal significant changes, as indicated by the already known partial corrections due to fragmentation contributions [19,20]. The 3 P [8] J term, instead, shows drastic changes from LO to NLO, including a change of sign, and the fragmentation corrections are large [19,20]. One may argue, hence, that higher-order corrections could change the p T /M dependence of the 3 S [8] 1 + κ 3 P [8] J polarization, especially towards low p T /M. To reach a situation where the polarization measurements would no longer exclude this contribution, so that the v 2 hierarchy predicted by NRQCD ( 1 S [8] 0 ∼ 3 S [8] 1 ∼ 3 P [8] J ) would be recovered, those future higherorder corrections should necessarily lead to a flatter dependence of λ ϑ on p T /M. The presently available quarkonium polarization measurements are already sufficiently precise to suggest that those future computations need to reach a suspiciously high level of fine tuning, but it is clear that new polarization measurements, of significantly improved precision, are needed to provide the ultimate constraint on the NRQCD hierarchy puzzle.
In particular, high-precision measurements confirming beyond doubt the scenario of vanishing and p T -independent polarization would leave zero margin for a 3 S [8] 1 + κ 3 P [8] J contribution other than unpolarized and p T -independent. The crucial importance of such measurements can be appreciated by considering the two possible conclusions: (a) if the cal- J term remains significantly different from unpolarized and shows some level of p T dependence (as is the case of the NLO computations), then we will be forced to conclude that its contribution is zero and quarkonium production is driven entirely by the 1 S [8] 0 channel; (b) if improved (higher order) calculations lead to λ ϑ ( 3 S [8] 1 + κ 3 P [8] J ) ∼ 0, the validity of the NRQCD v-scaling rules will be formally rescued, but the three octet terms would collapse to a single one, given that they would then have identical observable kinematic properties. This outcome would clearly point to the existence of a more natural and fundamental formulation of the factorization expansion, reflecting less degrees of freedom, valid at least in the conditions of mid-rapidity quarkonium production in high-energy pp collisions.
The best way to accurately probe the lowp T /M region is to use ϒ(nS) measurements, ideally with well-resolved 1S, 2S and 3S states to indirectly examine possible effects from the feed-down decay contributions of the P-wave states. Accessing the highp T /M kinematical domain (easier to study with J/ψ and ψ(2S) data, given the lower masses and higher production cross sections) will be crucial to probe if quarkonium production is completely dominated by a single unpolarized production mechanism (the 1 S [8] 0 octet term, say) or if there are two (or more) polarized terms that seemingly cancel each other in the existing data. High-precision measurements revealing no polarization, independently of p T , could no longer be explained with coincidental cancellations between cross section terms: if the present hints of cancellations result from a conspiracy between the low experimental precision and the limited kinematical domain covered by data, a p T -dependent residual polarization should become visible at sufficiently high p T . It should be easy to vastly reduce the statistical uncertainties of the presently available mid-rapidity LHC quarkonium polarization measurements, which are exclusively based on the 7 TeV data collected by CMS in 2011, corresponding to an integrated luminosity of only 5 fb −1 , much less than the 140 fb −1 collected at 13 TeV between 2016 and 2018. Polarization measurements reported by the LHCb collaboration [12,13,26] will also be very valuable to precisely investigate the low p T region, once the p T /M scaling studies [24] will have been extended to the forward rapidity range covered by that experiment. The systematic uncertainties, likely to become the dominating ones in a large range of the measured p T spectrum, can be reduced by performing the measurements in at least two (orthogonal) polarization frames and by also reporting measurements of the λ ϕ and λ ϑϕ azimuthal anisotropy parameters [21]. Measuring theλ frame-invariant parameter [27] in several polarization frames will also help uncovering potential biases.
It is important to note that the experimental conditions of the LHC data considered in this paper may represent our best chance of exploring the hierarchy and possible limit degeneracies in the NRQCD expansion. Indeed, given the very high collision energies, we can access mid-rapidity data up to relatively high p T /M values with a good statistical precision, so that we can study a phase space window where quarkonium production is dominated by 2 → 2 processes of the kind gg → Q Q + g. In this configuration the Q Q pre-resonance can "freely" assume all J values up to, in principle, J = 3 (with S = 0 or S = 1); there is, in principle, no a priori topological constraint removing some of these terms, nor imposing a specific J z (i.e., a polarization). Instead, data obtained in experimental conditions dominated by specific individual channels would necessary be affected by additional constraints, as in the case of singlet-driven production via virtual photon in e + e − collisions, or in kinematic domains where 2 → 1 production becomes important and the Q Q directly inherits the angular momentum state of the system of colliding partons. In short, it could well be that only data collected in optimal conditions, where no special angular momentum constraints are directly imposed on the produced Q Q, allow us to probe the existence of an "unbroken degeneracy" at a fundamental level.
In summary, we have shown the existence of a kinematic domain in high-energy pp collisions, p T /M 10 at midrapidity, where the description of quarkonium production using the current formulation of a v 2 expansion in 2S+1 L J Fock-states, together with NLO SDC calculations, becomes degenerate. Indeed, the observed lack of polarization and universal p T /M dependence of the yields can be described, indifferently, by one of the expansion terms alone, the 1 S [8] 0 octet, or by a combination of the other two, the 3 S [8] 1 + κ 3 P [8] J linear combination (and, obviously, by a mixture of the two cases, combined with unobservable proportions). This degeneracy is broken for p T /M < 10 by the polarization prediction, because the 3 S [8] 1 + κ 3 P [8] J term acquires a strong p T /M dependence, differentiating it from the constant and unpolarized 1 S [8] 0 octet. However, such change of regime is not at all observed in the existing measurements, which seamlessly prolong their unpolarized behaviour down to the lowest data point, at p T /M ∼ 2. This implies a strong violation of the v-scaling rules and an unexpected dominance of the 1 S [8] 0 term. Future high-precision quarkonium polarization measurements have the potential to rule out the current formulation of the factorization expansion beyond any possible recovery, if they confirm the unpolarized scenario suggested by the currently available data. Fine-tuned higher-order corrections could provide SDC calculations that would resuscitate the 3 S [8] 1 + κ 3 P [8] J term, but at the expense of turning it, necessarily, into another unpolarized term, thereby exposing a full degeneracy of the currently postulated Fock-space expansion. Such an outcome, rather than solving the riddle, would raise an even more compelling question: does a more essential and natural description of the experimental observations exist? The search for explanations of why and how the degeneracy appears in the considered phase space (midrapidity and/or high p T /M, in high-energy pp collisions), will certainly bring advances in the fundamental understanding of quarkonium production.
We would like to thank Hua-Sheng Shao, who kindly provided tables of the NLO calculations of the SDCs.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: Data sharing not applicable to this article as no datasets were generated or analysed during the current study.] Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .