Bottomonia production and polarization in the NRQCD with $k_T$-factorization. I: $\Upsilon(3S)$ and $\chi_b(3P)$ mesons

The $\Upsilon(3S)$ production and polarization at high energies is studied in the framework of $k_T$-factorization approach. Our consideration is based on the non-relativistic QCD formalism for bound states formation and off-shell production amplitudes for hard partonic subprocesses. The transverse momentum dependent (TMD, or unintegrated) gluon densities in a proton were derived from the Ciafaloni-Catani-Fiorani-Marchesini (CCFM) evolution equation as well as from the Kimber-Martin-Ryskin (KMR) prescription. Treating the non-perturbative color octet transitions in terms of the multipole radiation theory and taking into account feed-down contributions from radiative $\chi_b(3P)$ decays, we extract the corresponding non-perturbative matrix elements for $\Upsilon(3S)$ and $\chi_b(3P)$ mesons from a combined fit to $\Upsilon(3S)$ transverse momenta distributions measured by the CMS and ATLAS Collaborations at the~LHC energies $\sqrt s = 7$ and $13$~TeV and central rapidities. Then we apply the extracted values to describe the CDF and LHCb data on $\Upsilon(3S)$ production and to investigate the polarization parameters $\lambda_\theta$, $\lambda_\phi$ and $\lambda_{\theta\phi}$, %and $\tilde\lambda$, which determine the $\Upsilon(3S)$ spin density matrix. Our predictions have a good agreement with the currently available data within the theoretical and experimental uncertainties.

χ b (3P ) and χ b (2P ) decays and contributions from Υ(3S) decays are taken into account. The transverse momentum dependent gluon densities in a proton were derived from the Ciafaloni-Catani-Fiorani-Marchesini evolution equation as well as from the Kimber-Martin-Ryskin prescription. Treating the non-perturbative color octet transitions in terms of the mulitpole radiation theory, we extract the corresponding non-perturbative matrix elements for Υ(2S) and χ b (2P ) mesons from a combined fit to Υ(2S) transverse momenta distributions measured by the CMS and ATLAS Collaborations at the LHC energies √ s = 7 and 13 TeV and from the relative production rate R χ b (2P ) Υ(2S) measured by the LHCb Collaboration at √ s = 7 and 8 TeV. Then we apply the extracted values to investigate the polarization parameters λ θ , λ φ and λ θφ , which determine the Υ(2S) spin density matrix. Our predictions have a good agreement with the currently available data within the theoretical and experimental uncertainties.

Introduction
Since it was first observed, the production of charmonia and bottomonia in hadronic collisions remains a subject of considerable theoretical and experimental studies [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16]. The theoretical framework for the description of heavy quarkonia production and decays provided by the non-relativistic QCD (NRQCD) factorization [17,18]. This formalism implies a separation of perturbatively calculated short-distance cross-sections for the production of QQ pair in an intermediate Fock state 2S+1 L (a) J with spin S, orbital angular momentum L, total angular momentum J and color representation a from long-distance non-perturbative matrix elements (NMEs), which describe the transition of that intermediate QQ state into a physical quarkonium via soft gluon radiation.
However, NRQCD meets some difficulties in simultaneous description of the charmonia and bottomonia production cross section and polarization data, as we have explained in our previous paper [19]. So, for example, having the NMEs fixed from fitting the charmonia transverse momentum distributions, one disagrees with the polarization observables: if the dominant contribution comes from the gluon fragmentation into an octet QQ pair, the outgoing meson must have strong transverse polarization. The latter disagrees with the latest data [20][21][22][23][24], which show the unpolarized or even longitudinally polarized particles (so called "the polarization puzzle"). Moreover, the NMEs, obtained from the collider data, dramatically depend on the minimal transverse momentum used in the fits [25] and are incompatible with each other when obtained from fitting the different data sets.
A potential solution to this problem was proposed [26] in the framework of a model that interprets the soft final state gluon radiation as a series of color-electric dipole transitions. In this way the NMEs are represented in an explicit form inspired by the classical multipole radiation theory, that leads to unpolarized or only weakly polarized mesons either because of the cancellation between the 3 P contributions or as a result of two successive color-electric E1 dipole transitions in the chain 3 S 1 . This scenario was already successfully applied to describe the recent data on charmonia production and polarization [27,28]. Of course, it is important to investigate the bottomonia production within the same framework.
The data on Υ(nS) and χ b (mP ) mesons have been reported recently by the CMS [29,30], ATLAS [31] and LHCb [32,33] Collaborations at √ s = 7, 8 and 13 TeV. As it was shown [34][35][36][37][38], these data can be explained within the NRQCD, both in polarization and yield. Our present study continues the line started in the previous paper [19] and here we consider the production and polarization of Υ(2S) mesons. To preserve the consistency with our studies [27,28], we apply the k T -factorization QCD approach [39,40] to describe the perturbative production of the bb pair in the hard scattering subprocess. This approach is based on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) [41] or Ciafaloni-Catani-Fiorani-Marchesini (CCFM) [42] evolution equations, which resum large logarithmic terms proportional to ln s ∼ ln 1/x, important at high energies (or, equivalently, at low longitudinal momentum fraction x of proton carried by gluon). Resummation of the terms α n s ln n 1/x, α n s ln n µ 2 /Λ 2 QCD and α n s ln n 1/x ln n µ 2 /Λ 2 QCD up to all orders in the perturbative expansion results in Transverse Momentum Dependent (TMD) gluon distributions, that generalize the factorization of hadronic amplitudes beyond the conventional (collinear) DGLAP-based approximation. For the different aspects of the k T -factorization approach the reader may consult the reviews [43]. We determine the NMEs for Υ(2S) and χ b (2P ) mesons from the Υ(2S) transverse momentum distributions measured by the CMS [29,30] and ATLAS [31] Collaborations at √ s = 7 and 13 TeV and from the relative production ratio R Υ(2S) measured by the LHCb Collaboration at √ s = 7 and 8 TeV [44]. In the calculations we take into account the feed-down contributions from χ b (3P ), χ b (2P ) and Υ(3S) decays. Then, we make predictions for polarization parameters λ θ , λ φ , λ θφ (and frame-independent parameterλ), which determine the Υ(2S) spin density matrix and compare them to the currently available data [20,24].
The outline of our paper is the following. In Section 2 we briefly recall the basic steps of our calculations. In Section 3 we perform a numerical fit and extract the NMEs from the LHC data. Then we test the compatibility of the extracted NMEs with the available LHCb data on Υ(2S) transverse momentum distributions and Tevatron data on the Υ(2S) transverse momentum distributions and polarization. Our conclusions are collected in Section 4.

Theoretical framework
In the present paper we follow mostly the same steps as in our previous paper [19]. Our consideration is based on the off-shell gluon-gluon fusion subprocesses that represents the true leading order (LO) in QCD: where we listed all intermediate color states, J = 0, 1 or 2 and the four-momenta of all particles are indicated in the parentheses. The respective cross sections for 2 → 2 and 2 → 1 subprocesses are: where Q is the Υ and/or χ b meson, φ 1 and φ 2 are the azimuthal angles of initial off-shell gluons having the longitudinal momentum fractions x 1 and x 2 , p T and y are the transverse momentum and rapidity of produced mesons, y g is the rapidity of outgoing gluon and F is the off-shell flux factor [45]. The initial off-shell gluons have non-zero transverse momenta k 2 1T = 0, k 2 2T = 0 and an admixture of longitudinal component in the polarization vectors. So, the gluon spin density matrix is taken in the form where k T is the component of the gluon momentum perpendicular to the beam axis [39,40]. In the collinear limit k T → 0 this expression converges to the ordinary one In all other respects, we follow the standard QCD Feynman rules. As usual, the hard production amplitudes contain spin and color projection operators [46] that guarantee the proper quantum numbers of the state under consideration (see, for example, [19,27,28] for more details). The formation of bb bound states need additional explanation. We employ the mechanism 1 proposed in [26] and used previously [19,27,28]. A soft gluon with a small energy E ∼ Λ QCD is emitted after the hard interaction is over, bringing away the unwanted color and changing other quantum numbers of the produced CO system. In our calculations such soft gluon emission is described by a classical multipole expansion, in which the electric dipole (E1) transition dominates [47]. Only a single E1 transition is needed to transform a P -wave state into an S-wave state and the structure of the respective 3 P 1 + g amplitudes is given by [47]: where p The mechanism [26] is not connected to the choice of factorization scheme (k T or collinear), but represents a completely independent issue. corresponding particles and e µναβ is the fully antisymmetric Levi-Civita tensor. The transformation of color-octet S-wave state into the color-singlet S-wave state is treated as two successive E1 transitions 3 S states with J = 0, 1, 2. For each of these transitions we apply the same expressions (6) (8). The amplitudes (6) (8) lead to the fact that the final state bottomonia come unpolarized [26], either because of the cancellation between the 3 P contributions or as a result of two successive E1 transitions. This property remains true irrespectively of the numerical values of NMEs and only follows from the spin algebra. The expressions (6) (8) can be applied for both gluons and photons (up to an overall color factor) and can be used to calculate the polarization variables in radiative decays in feed-down processes.
As we did in our previous paper [19], we have tested several sets of TMD gluon densities in a proton. Two of them (A0 [48] and JH'2013 set 1 [49]) were obtained from CCFM equation where all input parameters were fitted to the proton structure function F 2 (x, Q 2 ).
Besides that, we have tested a parametrization obtained within the Kimber-Martin-Ryskin (KMR) prescription [50], which provides a method to construct the TMD quark and gluon densities from conventional (collinear) distributions. For the input, we have used recent LO NNPDF3.1 set [51].
The parton level calculations were performed using the Monte-Carlo event generator pegasus [52].

Numerical results
In the present paper we set the masses m Υ ( (1) 0 ] = 2.61 GeV 5 as obtained from the potential model calculations [54]. All the NMEs for Υ(3S) and χ b (3P ) mesons were derived in [19].

Fit of color octet NMEs
We performed a global fit to the Υ(2S) production data at the LHC and determined the corresponding NMEs for both Υ(2S) and χ b (2P ) mesons. We have included in the fitting procedure the Υ(2S) transverse momentum distributions measured by the CMS [29,30] and ATLAS [31] Collaborations at √ s = 7 and 13 TeV. To determine NMEs for χ b (2P ) mesons, we also included into the fit the recent LHCb data [44] on the radiative χ b (2P ) → Υ(2S) + γ decays collected at √ s = 7 and 8 TeV. We have excluded from our fit low p T region and consider only the data at p T > p cut T = 10 GeV, where the NRQCD formalism is believed to be most reliable.
We would like to mention here a few important points. First of all, we found that the 1 ] contributions to Υ(2S) production is almost the same in all kinematical regions probed at the LHC. Thus, the ratio 1 ]/dp T (9) can be well approximated by a constant for a wide Υ(2S) transverse momentum p T and rapidity y ranges at different energies, as seen in Fig. 1. We estimate the mean-square average r = 0.98 ± 0.005, which is practically independent on the TMD gluon density in a proton. So that, we construct the linear combination which can be only extracted from the measured Υ(2S) transverse momentum distributions. Then we use recent LHCb data [44] on the ratio of Υ(2S) mesons originating from the χ b (2P ) radiative decays measured at √ s = 7 and 8 TeV: From the known M r and R 1 ] and, therefore, reconstruct full map of color octet NMEs for both Υ(2S) and χ b (2P ) mesons.
The fitting procedure was separately done in each of the rapidity subdivisions (using the fitting algorithm as implemented in the commonly used gnuplot package [55]) under the requirement that all the NMEs be strictly positive. Then, the mean-square average of the fitted values was taken. The corresponding uncertainties are estimated in the conventional way using Student's t-distribution at the confidence level P = 80%. The results of our fits are collected in Table 1. For comparison, we also presented there the NMEs obtained in the conventional NLO NRQCD by other authors [37]. The corresponding χ 2 /d.o.f. are listed in Table 2, where we additionally show their dependence on the minimal Υ(2S) transverse momenta involved into the fit p cut T . As one can see, the χ 2 /d.o.f. tends to decrease when p cut T grows up and best fit of the LHC data is achieved with A0 gluon, although other gluon densities also return reliable χ 2 /d.o.f. values. We note that including into the fit the latest CMS data [30] taken at √ s = 13 TeV leads to 2 3 times higher values of χ 2 /d.o.f. We have checked that this is true for both the k T -factorization and collinear approaches 2 and, therefore, it could be a sign of some inconsistency between these CMS data and all other measurements.
All the data used in the fits are compared with our predictions in Figs. 2 4. The shaded areas represent the theoretical uncertainties of our calculations, which include the scale uncertainties, uncertainties coming from the NME fitting procedure and uncertainties connected with the choice of the intermediate color-octet mass, added in quadrature. To estimate the scale uncertainties the standard variations µ R → 2µ R or µ R → µ R /2 were introduced with replacing the A0 and JH'2013 set 1 gluon densities by A0+ and JH'2013 set 1+, or by A0− and JH'2013 set 1− ones. This was done to preserve the intrinsic correspondence between the TMD gluon set and the scale used in the CCFM evolution [48,49]. To estimate the uncertainties connected with the intermediate color-octet mass we have varied amount of energy E emitted in the course of transition of unbound color octet bb pair into the observed bottomonium by a factor of 2 around its default value E = Λ QCD . One can see that we have achieved a reasonably good description of the CMS [29,30] and ATLAS [31] data in a whole p T range within the experimental and theoretical uncertainties for the Υ(2S) transverse momentum distributions. The ratios R Υ(2S) measured by the LHCb Collaboration [44] at √ s = 7 and 8 TeV are also reproduced well, see Fig. 4. Finally, we have checked our results with the data, not included into the fit procedure: namely, rather old CDF data [56] taken at the √ s = 1.8 TeV and LHCb data [32,33] taken in the forward rapidity region 2 < y < 4.5 at √ s = 7, 8 and 13 TeV (see Fig. 5). As one can see, we acceptably describe all the data above. Moreover, we find that the KMR gluon is only one which is able to reproduce the measurements in the low p T region.
The CMS Collaboration has measured all of these polarization parameters for Υ(2S) mesons as functions of their transverse momentum in three complementary frames: the Collins-Soper, helicity and perpendicular helicity ones at √ s = 7 TeV [20]. The frameindependent parameterλ = (λ θ + 3λ φ )/(1 − λ φ ) has been additionally studied. The CDF Collaboration has measured λ θ andλ parameters in the helicity frame at √ s = 1.96 TeV [24]. As it was done previously [19], to estimate λ θ , λ φ , λ θφ andλ we generally follow the experimental procedure. We collect the simulated events in the kinematical region defined by the experimental setup, generate the decay lepton angular distributions according to the production and decay matrix elements and then apply a three-parametric fit based on (12).
Our predictions are shown in Figs. 6 9. The calculations were done using the A0 gluon density which provides a best description of the measured Υ(2S) transverse momenta distributions. As one can see, we find only weak or zero polarization in the all kinematic regions, that agrees with the CMS and CDF measurements. The similar results we have obtained earlier for charmonia (J/ψ, ψ ) and Υ(3S) polarization [19,27,28]. Thus, we conclude that the approach [26], which is a corner stone of our consideration, results in a self-consistent and simultaneous description of the entire charmonia family, Υ(2S) and Υ(3S) production data and therefore can provide an easy and natural solution to the long-standing quarkonia production and polarization puzzle.

Conclusion
We have considered the Υ(2S) production at the Tevatron and LHC in the framework of k T -factorization approach. Our consideration was based on the off-shell production amplitudes for hard partonic subprocesses (including both color-singlet and color-octet contributions), NRQCD formalism for the formation of bound states and TMD gluon densities in a proton. The latter were derived from the CCFM evolution equation and KMR scheme. Treating the nonperturbative color octet transitions in terms of multipole radiation theory and taking into account feed-down contributions from the radiative χ b (3P ) and χ b (2P ) decays and contribution from Υ(3S) decays, we extracted long-distance non-perturbative NRQCD matrix elements for Υ(2S) and χ b (2P ) mesons from a fit to Υ(2S) transverse momentum distributions measured by the CMS and ATLAS Collaborations at √ s = 7 and 13 TeV and from the relative production rates R χ b (2P ) Υ(2S) measured by the LHCb Collaboration at √ s = 7 and 8 TeV. Then we estimated polarization parameters λ θ , λ φ , λ θφ and frameindependent parameterλ which determine the spin density matrix of Υ(2S) mesons. We show that treating the soft gluon emission as a series of explicit color-electric dipole transitions within the NRQCD leads to unpolarized Υ(2S) production at moderate and large transverse momenta, that is in agreement with the LHC data.
1 ] /GeV 3 0.0181 ± 0.0007 0.0117 ± 0.0007 0.0074 ± 0.0004 0.0109 ± 0.0014      Fig. 2. The experimental data are from CMS [29,30]. Figure 6: The polarization parameters λ θ , λ φ , λ θφ andλ of Υ(2S) mesons calculated in the CS frame as function of its transverse momentum at √ s = 7 TeV. The A0 gluon density is used. The blue and red histograms correspond to the predictions obtained at |y| < 0.6 and 0.6 < |y| < 1.2, respectively. The experimental data are from CMS [20]. Figure 7: The polarization parameters λ θ , λ φ , λ θφ andλ of Υ(2S) mesons calculated in the helicity frame as function of its transverse momentum at √ s = 7 TeV. Notation of all histograms is the same as in Fig. 6. The experimental data are from CMS [20]. Figure 8: The polarization parameters λ θ , λ φ , λ θφ andλ of Υ(2S) mesons calculated in the perpendicular helicity frame as function of its transverse momentum at √ s = 7 TeV. Notation of all histograms is the same as in Fig. 6. The experimental data are from CMS [20]. Figure 9: The polarization parameters λ θ andλ of Υ(2S) mesons calculated in the helicity frame as function of its transverse momentum at √ s = 1.96 TeV. Notation of all histograms is the same as in Fig. 6. The experimental data are from CDF [24].