# On multiple gluon exchange in hadroproduction

###### Abstract

We consider a contribution to hadroproduction in which the meson production is mediated by three-gluon partonic state, with two gluons coming from the target and one gluon from the projectile. This mechanism involves double gluon density and thus it enters at a non-leading twist, but it is enhanced at large energies due to large double gluon density at small . The three-gluon contribution to hadroproduction is calculated within perturbative QCD in the factorization framework, and it is found to provide a significant correction to the standard leading twist cross-section at the energies of the Tevatron or the LHC. The results are given as differential -dependent cross-sections for polarization components.

On multiple gluon exchange in hadroproduction

Mariusz Sadzikowski

Institute of Physics, Jagiellonian University, Reymonta 4, 30-059 Kraków, Poland

E-mail: mariusz.sadzikowski@uj.edu.pl

\abstract@cs

## 1 Introduction

Understanding of charmonia and bottomonia production in high energy and collisions has been a challenge for the theory of strong interactions since first Tevatron measurements were published of prompt production [1]. The Tevatron data showed huge excess of the measured cross-section over the simplest LO QCD predictions, reaching almost two orders of magnitude at large transverse momenta [2]. This naïve picture of the process is given by standard, color singlet mechanism (CSM) in the collinear approximation. In this mechanism it is assumed, that all the vector quarkonium quantum numbers and kinematic properties are generated already at the partonic level. The production process of the heavy quark – anti quark pair is then driven by gluons, , and the final state gluon emission of the gluon is necessary to produce the pair with the quantum numbers of the vector meson. Due to presence of large scales — the quark mass and large transverse momentum — the process may be computed within perturbative QCD. In the collinear approximation, the emitted gluon recoil generates the transverse momentum of the meson. As said above, however, this simple description fails badly to describe the data and needed to be modified. In particular, the transverse momentum dependence of the collinear CSM is much too steep.

### 1.1 Color octet model

Since it become clear that the collinear color singlet approximation is insufficient, a few mechanisms have been proposed that are able to explain reasonably well the measured heavy quarkonia production. Currently two most successful approaches to prompt vector meson hadroproduction are based on perturbative QCD. The first one is based on the so-called color octet mechanism (COM), that assumes that the pair may be produced with quantum numbers different from the final state meson, in particular in the color octet representation [3, 4]. Originally the COM has been motivated by the heavy quarkonium effective theory in which higher Fock states of the quarkonium wave function may contain an octet pair with some amplitude of the order of . Probably, it may be also viewed as a fragmentation process of partonic state into the meson, somewhat similar to or quark fragmentation into or mesons. What is crucial for effective predictivity, in the COM one assumes existence of universal transition amplitudes from the pair with given quantum numbers to the mesons, provided matching of quark and meson kinematic variables. The COM has been developed up to the next-to-leading order (NLO) within collinear approximation [5, 6, 7], and it was shown to provide a good global fit of prompt quarkonium production data [6], except of the quarkonium polarization [7]. The COM model does not provide a satisfactory description of prompt polarization as a function of transverse momentum neither at the leading order nor at the NLO. Therefore, however sound and sophisticated, the COM may be still incomplete.

### 1.2 factorisation approach

The other rather successful approach to vector quarkonia hadroproduction attempts to explain the wide transverse momentum distribution by taking into account the non-zero transverse momentum of initial state gluons assuming the -factorization (or high energy factorization) [8, 9, 10, 11, 12]. The most interesting quarkonia production at the Tevatron and the LHC occur in the kinematic domain where the incoming parton (gluon) momentum fractions are small, typically . In this region radiation in the QCD evolution may generate sizable transverse momentum, , of the gluons, resulting then with a modified (broader) distribution of quarkonium transverse momentum. In the -factorization approach the standard, collinear parton densities, e.g. for the gluon, are replaced by unintegrated parton densities, e.g. , with a LL constraint . The -factorization approach was applied to prompt quarkonia production at the LO, assuming both color singlet and color octet scenarios. The emerging picture is quite encouraging, but not fully clear yet. A recent description of the recent LHC data within the -factorization CSM approach was quite successful, including even meson polarization [11], but the older description of the Tevatron data required also the color octet contributions [9]. Therefore, in this moment the -factorizing approach also suffers from some deficiencies in providing the consistent global picture of prompt quarkonia production.

## 2 The triple gluon correction

### 2.1 Motivation

As discussed in the previous section, the existing approaches to vector quarkonia hadroproduction, although reproducing well many features of the data, do not provide complete and fully satisfactory global description. Thus it may be still necessary to include yet other production mechanism. One of the potentially important corrections to the standard CSM and COM contributions may be driven by multiple scattering effects. At the lowest order (assuming matching of quark and meson quantum numbers) the correction is driven by a partonic amplitude of three-gluon fusion into heavy quark-antiquark pair, [13]. Since one needs to drag two glouns out of one hadron, this correction employs the double gluon distribution and enters, therefore, beyond the leading twist. It implies a power suppression of this correction w.r.t. the standard, two-gluon cross-sections. The double gluon distribution, however, at small gluon , provides a large enhancement factor , that may well reach about 20 in the relevant kinematic domain. The first estimates suggested that this multiple scattering effect may be even the dominant contribution of the total vector quarkonia cross-section, leaving only little space for the color octet contributions [13]. Moreover, the multiple scattering effects should be also important in vector meson productions in high energy collisions with nuclei, where the nucleus mass number provides an additional enhancement factor of multiple gluon distributions. Therefore we performed a detailed estimate of this rescattering (or triple gluon) correction for the Tevatron and the LHC energies beyond the collinear, leading logarithmic accuracy of Ref. [13].

### 2.2 The framework

The calculations were performed within the -factorization approach. The relevant diagram for an inclusive process, is shown in Fig. 1a. The total amplitude includes amplitudes of this diagram and the ‘upside-down’ diagram in which the two gluons originate from projectile . When squaring the amplitude, one finds diagrams with topology shown in Fig. 1b, with four and two gluons in the -channel, and (not shown) diagrams with three gluon -channel states originating from both and . The latter ‘interference’ diagrams may be neglected as the three gluon state evolution is known to have no energy enhancement, contrary to the two and four gluon states that drive the contributions shown in Fig. 1b. After performing the phase space integrals over remnants and in this contribution one recovers double and single unintegrated gluon distributions originating from and respectively. The three gluon-fusion amplitude into the meson (the middle part of the diagram shown in Fig. 1a) is described in terms of a known impact factor, [14], with polarization vectors in the helicity frame. The proper normalization of the amplitudes is obtained by matching to collinear cross-sections of single and double parton scattering. We assume factorization of the double gluon density and factorization of the impact parameter dependence. With these assumptions the final formula for the triple gluon contribution to production takes the form:

(2.0) |

where and are meson longitudinal momentum fractions of the projectiles, rapidity , p is the meson transverse momentum, , is the meson polarization, k, l and are the gluons transverse momenta, and we combine numerical constants, color factors etc. into a normalization constant .

Note the emergence of , the effective parameter describing parton transverse density and probability of multiple scattering. It appears as a result of forcing two partons from one hadron to be in the same position in the transverse space. Since, in the double gluon density the values of gluon are different, we used off-diagonal gluon densities and included the Shuvaev factor, [15].

### 2.3 Results

In the numerical evaluations unintegrated gluon densities were used derived from the CT10 collinear gluon density [16] using Kimber-Martin-Ryskin approach [17] with the hard scale given by the transverse meson mass, . The running strong coupling constants of a gluon with virtuality was evaluated at the scale , with . We set the multiple scattering parameter value mb, in accordance with Refs. [19, 20].

In Fig. 2a the results of numerical evaluation of the triple gluon correction to the differential meson production cross section (including the branching ratio of the meson decay to muons) are shown for the Tevatron energy ( GeV). For a reference, we show in Fig. 2b a plot taken from Ref. [6], in which results of CSM and of COM fits are presented at LO and NLO accuracy together with CDF data. From the comparison of Fig. 2a and Fig. 2b it is clear, that the triple gluon correction is similar in the magnitude and shape to the standard CSM contribution at the NLO. In more detail, the triple-gluon contribution is found to be larger than the CSM contribution in the whole range of transverse momentum and it reaches about 20 – 25% of the measured cross-section at low transverse momenta. At larger transverse momenta this contributions the relative importance of this correction diminishes and already at GeV, this contribution is negligible.

In Fig. 2a we also show polarized components of the triple gluon correction in the helicity frame. Clearly, at low transverse momentum the longitudinal and total transverse cross-sections are similar, and with increasing the longitudinal component becomes dominant, saturating the total cross-section.

The results for the LHC exhibit very similar pattern, as the results for the Tevatron, so we do not show them here, leaving the broader and more detailed presentation to the forthcoming paper [21].

## 3 Conclusions

We estimated corrections to prompt production at the Tevatron and the LHC coming from a triple gluon incoming partonic state. They are found to be sizable at small transverse momenta, up to 20 – 25 % of the measured cross-section and negligible at large transverse momenta. The polarization from this mechanism varies with , from negligible polarization at small momenta to longitudinal polarization at large . The correction is large enough to affect the fits to the existing, precise data on hadroproduction. The results may be also relevant for production in proton–nucleus collisions, where the multiple gluon exchanges are enhanced by the nucleus mass number.

## References

- [1] J. P. Lansberg, Int. J. Mod. Phys. A 21 (2006) 3857 [hep-ph/0602091].
- [2] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 69 (1992) 3704; Phys. Rev. Lett. 79 (1997) 578.
- [3] G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 51 (1995) 1125 [Erratum-ibid. D 55 (1997) 5853].
- [4] P. L. Cho and A. K. Leibovich, Phys. Rev. D 53 (1996) 150; P. L. Cho and A. K. Leibovich, Phys. Rev. D 53 (1996) 6203.
- [5] M. Butenschoen and B. A. Kniehl, Phys. Rev. Lett. 104 (2010) 072001; Phys. Rev. Lett. 106 (2011) 022003; Phys. Rev. Lett. 107 (2011) 232001.
- [6] M. Butenschoen and B. A. Kniehl, Phys. Rev. D 84 (2011) 051501.
- [7] M. Butenschoen and B. A. Kniehl, Phys. Rev. Lett. 108 (2012) 172002.
- [8] P. Hagler, R. Kirschner, A. Schafer, L. Szymanowski and O. V. Teryaev, Phys. Rev. Lett. 86 (2001) 1446; Phys. Rev. D 63 (2001) 077501.
- [9] S. P. Baranov, Phys. Rev. D 66 (2002) 114003.
- [10] S. P. Baranov and N. P. Zotov, JETP Lett. 88 (2008) 711.
- [11] S. P. Baranov, A. V. Lipatov and N. P. Zotov, Phys. Rev. D 85 (2012) 014034.
- [12] B. A. Kniehl, D. V. Vasin and V. A. Saleev, Phys. Rev. D 73 (2006) 074022.
- [13] V. A. Khoze, A. D. Martin, M. G. Ryskin and W. J. Stirling, Eur. Phys. J. C 39 (2005) 163.
- [14] A. Bzdak, L. Motyka, L. Szymanowski and J. -R. Cudell, Phys. Rev. D 75 (2007) 094023.
- [15] A. G. Shuvaev, K. J. Golec-Biernat, A. D. Martin and M. G. Ryskin, Phys. Rev. D 60 (1999) 014015.
- [16] H. L. Lai, M. Guzzi, J. Huston, Z. Li, P. M. Nadolsky, J. Pumplin and C.-P. Yuan, Phys. Rev. D 82 (2010) 074024.
- [17] M. A. Kimber, A. D. Martin and M. G. Ryskin, Phys. Rev. D 63 (2001) 114027.
- [18] D. Acosta et al. [CDF Collaboration], Phys. Rev. D 71 (2005) 032001.
- [19] F. Abe et al. [CDF Collaboration], Phys. Rev. D 56 (1997) 3811.
- [20] V. M. Abazov et al. [D0 Collaboration], Phys. Rev. D 81 (2010) 052012.
- [21] L. Motyka and M. Sadzikowski, in preparation.