Triple-parton scatterings in proton–nucleus collisions at high energies

A generic expression to compute triple-parton scattering (TPS) cross sections in high-energy proton–nucleus (pA) collisions is derived as a function of the corresponding single-parton cross sections and an effective parameter encoding the transverse parton profile of the proton. The TPS cross sections are enhanced by a factor of about 9A≃2000\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$9\,A\simeq 2000$$\end{document} in pPb as compared to those in proton–nucleon collisions at the same center-of-mass energy. Estimates for triple charm (cc¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{c}\overline{\mathrm{c}}$$\end{document}) and bottom (bb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{b}\overline{\mathrm{b}}$$\end{document}) production in pPb collisions at LHC and FCC energies are presented based on next-to-next-to-leading-order calculations for cc¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{c}\overline{\mathrm{c}}$$\end{document} and bb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{b}\overline{\mathrm{b}}$$\end{document} single-parton cross sections. At sNN=8.8\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{_{\textsc {nn}}}}= 8.8$$\end{document} TeV, about 10% of the pPb events have three cc¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{c}\overline{\mathrm{c}}$$\end{document} pairs produced in separate partonic interactions. At sNN=63\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{_{\textsc {nn}}}}= 63$$\end{document} TeV, the pPb cross sections for triple-J/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi $$\end{document} and triple-bb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{b}\overline{\mathrm{b}}$$\end{document} are O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathscr {O}}}$$\end{document}(1–10 mb). In the most energetic collisions of cosmic rays in the upper atmosphere, equivalent to sNN≈\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{_{\textsc {nn}}}}\approx $$\end{document} 400 TeV, the TPS cc¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{c}\overline{\mathrm{c}}$$\end{document} cross section equals the total p–Air inelastic cross section.


Introduction
The extended nature of hadronic systems and their growing parton density when probed at increasingly higher collision energies, makes it possible to produce multiple particles with large transverse momentum and/or mass ( p 2 T + m 2 3 GeV) in independent multiparton interactions (MPIs) in high-energy proton-(anti)proton (pp, pp) collisions [1][2][3][4][5]. Many experimental final states -involving the concurrent production of heavy-quarks, quarkonia, jets, and gauge bosons -have been found to be consistent with double parton scatterings (DPS) processes at Tevatron (see e.g. [6]) and the LHC (see e.g. [7][8][9] for a selection of the latest a e-mail: david.d'enterria@cern.ch b e-mail: snigirev@lav01.sinp.msu.ru results). Multiple hard parton interaction rates depend chiefly on the transverse overlap of the matter densities of the colliding hadrons, and provide valuable information on (1) the badly known 3D parton profile of the proton, (2) the unknown energy evolution of the parton density as a function of impact parameter (b), and (3) the role of many-parton correlations in the hadronic wave functions [10]. In our previous work [11][12][13][14], we highlighted the importance of studying DPS also in proton-nucleus (pA) and nucleus-nucleus (AA) collisions, as a complementary means to improve our understanding of hard MPIs in pp collisions. The larger transverse parton density in a nucleus (with A nucleons) compared to that of a proton, results in enhanced DPS cross sections coming from interactions where the two partons of the nucleus belong to the same or to two different nucleons, providing thereby useful information on the underlying multiparton dynamics [11][12][13][14][15][16][17][18][19][20][21].
The possibility of triple-parton scatterings (TPS) in hadronic collisions has also been considered in the literature [22][23][24][25][26], and estimates of their expected cross sections have been recently provided for pp collisions [27]. In this paper, we extend our latest work and derive for the first time quantitative estimates of the cross sections for observing three separate hard interactions in a pA collision through a factorized formula which depends on the underlying single-parton scattering (SPS) cross sections normalized by the square of an effective cross section σ eff,tps , characterizing the transverse area of triple-partonic interactions, that is closely related to the DPS-equivalent σ eff,dps parameter [27]. The paper is organized as follows. In Sect. 2, we review the theoretical expression for TPS cross sections in generic hadron-hadron collisions, first expressed as a convolution of SPS cross sections and generalized parton densities dependent on parton fractional momentum x, virtuality Q 2 , and impact parameter b, and then in their factorized form as a function of σ eff,tps . In Sect. 3, a generic expression for TPS cross sections in pA collisions is presented based on realistic parametrizations of the nuclear transverse profile. As concrete numerical examples, Sect. 4 provides estimates for triple charm (cc) and bottom (bb) cross sections from independent parton scatterings in proton-lead (pPb) collisions at the LHC and future circular collider (FCC) [28] energies, as well as in protonair collisions at the highest energies observed so far, based on next-to-next-to-leading-order (NNLO) calculations of the corresponding SPS cross sections. The main conclusions are summarized in Sect. 5.

Triple-parton-scattering cross sections in hadron-hadron collisions
In a generic hadronic collision, the inclusive TPS cross section from three independent hard parton scatterings (hh → abc) can be written as a convolution of generalized parton distribution functions (PDF) and elementary cross sections summed over all involved partons [22][23][24][25][26][27]29] σ tps Here, Γ i jk 3 ) encode all the parton structure information of relevance for TPS.
Without any loss of generality, any TPS cross section can always be expressed in a very economical form as a triple product of single inclusive cross sections, where σ eff,tps encodes all the unknowns related to the generalized PDFs and their correlations in space, color, flavor, p T ,... The value of σ eff,tps can be estimated theoretically making a few common approximations. First, the triple-PDF are commonly assumed to be factorizable in terms of longitudinal and transverse components, i.e.
describes the transverse parton density of the hadron, often considered a universal function for all types of partons, from which the corresponding hadron-hadron overlap function is derived: Making the further assumption that the longitudinal components reduce to the product of independent single PDF, D i jk , the effective TPS cross section has a very simple geometric interpretation as the inverse of the cube of the integral of the hadron-hadron overlap function: which is closely related to the similar quantity, determined in DPS measurements. In the proton-proton case, making use of Eqs. (5), (6) and (7), for a wide range of proton transverse parton profiles f (b), we found a simple relationship between the effective DPS and TPS cross sections: which, for the typical σ eff,dps = 15 ± 5 mb values extracted from a broad range of DPS measurements at Tevatron [6] and LHC [4,[6][7][8][9], translates into This data-driven numerical value allows for the computation of any TPS cross section in pp collisions via Eq. (3), once the corresponding SPS cross section, Eq. (2), is known.
In the next section, we extend and exploit these results for the pA case. We note that, since we start off with a numerical value for σ eff,dps that is directly obtained from data, our approach effectively takes into account effects which go beyond pure geometrical quantities computed in terms of transverse overlap of parton densities. Indeed, the fact that the experimentally extracted σ eff,dps = 15 ± 5 mb values are about a factor of 2 smaller (i.e., that the DPS cross sections are about twice larger) than expected from Eq. (7) for a "naive" proton profile, is indicative of the presence of "beyond mean-field" effects, including e.g. parton correlations in momentum, flavor, color, x,... Perturbative partonic correlations due to 1-to-2 parton splittings have been discussed for the DPS case in [30,31] in pp collisions and in [21] for pA collisions. The first attempt to estimate those for TPS in pp collisions was done in [26]. As we will see later, however, any such potential effects are of less relevance in pA than in pp collisions, since the majority of the TPS yields for the former arise from interactions among different "target" nucleons for which any potential partonic correlations are of smaller magnitude than in the pp case.

Triple-parton-scattering cross sections in proton-nucleus collisions
The first theoretical analysis of double-and triple-parton scattering cross sections in pA collisions was done in [15,16], where large enhancements of the DPS and TPS yields were anticipated. Here we focus on the derivation of a compact "pocket formula" for TPS scattering in pA collisions, particularly useful for phenomenological and numerical applications. The starting point of our approach is to remember that the parton flux in pA compared to pp is enhanced by the nucleon number A and, modulo shadowing effects in the nuclear PDF [32], the single-parton cross section for any hard process is that of proton-nucleon (pN) collisions (with N = p,n including their appropriate relative fraction in the nucleus) scaled by the factor A [33], Here dz is the nuclear thickness function given by the integral of the nuclear parton density function (commonly parametrized in terms of a "Woods-Saxon" Fermi-Dirac distribution [34]) over the longitudinal direction with respect to the impact parameter b between the colliding proton and nucleus, normalized to d 2 b T A (b) = A. In order to obtain a TPS "pocket formula" of the form of Eq. (3) for pA collisions, we follow the approach developed in our previous work for the DPS case [11][12][13][14]. The TPS pA cross section is thus obtained from the sum of three contributions: -A "pure TPS" cross section, given by Eq. (3) for pN collisions scaled by A, namely: -A second contribution, involving interactions of partons from two different nucleons in the nucleus, depending on the square of T A , where the factor (A−1)/A is introduced to account for the difference between the number of nucleon pairs and the number of different nucleon pairs. -A third term, involving interactions among partons from three different nucleons, depending on the cube of T A , σ tps,3 pA→abc = σ tps pN→abc · σ 2 eff,tps · C pA , with The factor (A−1)(A−2)/A 2 is introduced to take into account the difference between the total number of nucleon TPS and that of different nucleon TPS.
The inclusive TPS cross section for three hard parton subprocesses a, b, and c in pA collisions is thus obtained from the sum of the three terms (11), (12), and (14): which is enhanced by the factor in parentheses compared to the corresponding TPS cross section in pN collisions scaled by A. The analytical (unintegrated) expression for TPS cross sections in proton-nucleus collisions was first derived in [15], but our compact expression (16) provides a more useful formula to easily derive numerical TPS estimates for any pA system for phenomenological purposes. The value of this enhancement factor, as well as the relative role of each one of the three TPS components, can be obtained for pPb evaluating the integrals (13) and (15) using the standard Fermi-Dirac spatial density for the lead nucleus (A = 208, radius R A = 6.36 fm, and surface thickness a = 0.54 fm) [34]. The first integral is identical to the overlap function at zero impact parameter for the corresponding AA collision, F pA = (A − 1)/A T AA (0) = 30.25 mb −1 [11][12][13][14]. The second one can be obtained by means of a Glauber Monte Carlo (MC) [33] and amounts to C pA = 4.75 mb −2 . From the relationship (8) between the effective DPS and TPS cross sections, and the experimental σ eff,dps = 15 ± 5 mb value [4,[6][7][8][9], we can finally determine the relative importance for pPb of the three TPS terms of Eq. (16): σ tps,1 pA→abc : σ tps,2 pA→abc : σ tps,3 pA→abc = 1 : 4.54 : 3.56. Namely, in pPb collisions, 10% of the TPS yields come from partonic interactions within just one nucleon of the lead nucleus, 50% involve scatterings within two nucleons, and 40% come from partonic interac-tions in three different Pb nucleons. The fact that, for heavy nuclei, the second and third terms of Eq. (16) are much larger than the first one was first pointed out in Ref. [15], and we provide here an exact numerical evaluation for the pPb case. The sum of the three contributions in Eq. (16) amounts to 9.1, namely the TPS cross sections in pPb are about nine times larger than the naive expectation based on A-scaling of the corresponding pN TPS cross sections, Eq. (11). We note that in the DPS case the equivalent pA enhancement factor was [1 + σ eff,dps F pA /A] 3 [11][12][13][14]. The final formula for TPS in the proton-nucleus case reads where the effective TPS pA cross section in the denominator depends on the effective pp one and on pure geometric quantities directly derivable from the well-known nuclear transverse profile: σ 2 eff,tps,pA = A/σ 2 eff,tps + 2.46 F pA /σ eff,tps + C pA where the latter equality is obtained using Eqs. (8)- (9). The effective TPS cross section in the pPb case amounts thereby to σ eff,tps,pA = 0.29 ± 0.05 mb. This value is very robust with respect to the parametrization of the underlying proton and nucleus transverse profiles. Indeed, by using simplified Gaussian proton and nucleus transverse densities, all relevant factors in Eq. (16) can be analytically calculated, and the effective TPS pA cross section can be simply written as a function of the proton and nucleus radii: σ 2 eff,tps,pA = 3/4 σ 2 eff,dps /{A[1+9/2A (r p /R A ) 2 +4A 2 (r p /R A ) 4 ]}, which amounts to σ eff,tps,pA 0.28 mb (fixing r p so as to σ eff,dps = 15 mb), in perfect agreement with our more accurate estimate above.

Triple cc and bb production cross sections in pA collisions
As concrete numerical examples of our calculations, following our previous similar pp study [27], we compute the charm (pPb → cc + X) and bottom (pPb → bb + X) TPS cross sections first at the LHC and FCC center-of-mass (c.m.) energies, and then also in proton-air collisions of relevance for ultra-high-energy cosmic-ray collisions in the atmosphere. These processes are dominated by gluon-gluon scattering gg → qq at low parton fractional momentum x, for which the DPS and TPS mechanisms have a growing contribution to the total inclusive production at increasing c.m. energy. This expectation has been discussed for the DPS case in [35], and we extend those studies to the TPS case here. The TPS heavyquark cross sections are computed via Eq. (17) for m = 1, i.e. σ tps pPb→cc,bb = (σ sps pN→cc,bb ) 3 /(6 σ 2 eff,tps,pA ) with σ eff,tps,pA given by (18), and σ sps pN→cc,bb is calculated via Eq. (2) at NNLO accuracy using a modified version [36] of the Top++ (v2.0) code [37]. Top++ is run with N f = 3, 4 light flavors, charm and bottom pole masses set to m c,b = 1.67, 4.66 GeV [38], default renormalization and factorization scales set to μ R = μ F = 2 m c,b , and using the NNLO ABMP6 PDF of the proton [39] and the nuclear PDF modification factors of the Pb nucleus given by EPS09-NLO [32]. The PDF uncertainties include those from the proton and nucleus, as obtained from the corresponding 28 (30) eigenvalues of the ABMP16 (EPS09) sets, combined in quadrature. The dominant uncertainty is that linked to the theoretical scale choice, which is estimated by modifying μ R and μ F within a factor of 2. In the pp case, such a theoretical NNLO setup yields SPS heavyquark cross sections which are larger by up to 20% compared to the NLO [40,41] predictions at the LHC, reaching a better agreement with the experimental data, and showing a much reduced scale uncertainty (±50%, ±15% for cc, bb) [36]. In the pPb case, the inclusion of EPS09 nuclear shadowing reduces moderately the total charm and bottom cross sections in pN compared to pp collisions, by about 10% (13%) and 5% (10%) at the LHC (FCC). Since the TPS pPb cross section go as the cube of σ sps pN→qq , the impact of shadowing is amplified and leads to 15-35% reductions with respect to the result obtained if one used the pp (instead of the pN) SPS cross section in Eq. (17). At √ s = 5.02 TeV, our theoretical SPS prediction (σ sps,nnlo pPb→cc = 650 ± 290 sc ± 60 pdf mb) agrees well with the ALICE total D-meson measurement [42] extrapolated using [40] to a total charm cross section (σ alice pPb→cc = 640 ± 60 stat +60 −110 syst mb, Fig. 1 left). Table 1 collects the heavy-quark cross sections and associated uncertainties predicted in pPb collisions at the nominal LHC and FCC c.m. energies. The large SPS cc cross section at the LHC (∼1 b) results in triple-cc cross sections from independent parton scatterings amounting to about 20% of the inclusive charm yields. Since the total inelastic pPb cross sections are σ pPb = 2.2, 2.4 b at √ s nn = 8.8 and 63 TeV [33], charm TPS takes place in about 10% of the pPb events at 8.8 TeV. At the FCC, the theoretical TPS charm cross section even overcomes the inclusive charm one. Such a seemingly "unphysical" result indicates that quadruple, quintuple,... parton-parton scatterings are expected to produce extra cc pairs with non-negligible probability in pPb at √ s nn = 63 TeV. We recall that inclusive cross sections can be related to the factorial moments of the multiplicity distribution and, thus, SPS, DPS, TPS... cross sections are not bounded by the inelastic cross section [43]. The huge TPS cc cross sections at the FCC will make triple-J/ψ production     observable. Indeed, the SPS J/ψ cross section corresponds to a few percent of the cc one [11][12][13][14], which translates 1 into σ (3 × J/ψ + X) ≈ 1 mb. Triple-bb cross sections remain comparatively small, in the 0.1 mb range, at the LHC but 1 We note that the large production of double quarkonia seen in the data seem to indicate values of σ eff,dps smaller than the (15 ± 5 mb) world-average used here [44]. Therefore, one would also expect that the corresponding TPS yields for J/ψ could well be above our 1 mb order-of-magnitude estimate.
reach ∼10 mb (i.e. 3% of the total inclusive bottom cross section) at the FCC. Figure 1 shows pPb cross sections over √ s nn ≈ 40 GeV-100 TeV for SPS (solid bands), TPS (dashed bands) for charm (left) and bottom (right) production, and total inelastic (dotted curve, in both plots). The TPS cross sections are small at low energies but rise fast with √ s, as the cube of the SPS cross section evolution. Whenever the theoretical central value of the TPS cross section overcomes the inclusive charm cross section, indicative of multiple (beyond three) cc-pair production, we equalize it to the latter. Above √ s nn ≈ 25 TeV, the total charm and inelastic pPb cross sections are equal implying that the average number of charm pairs produced in pPb collisions is larger than one. In the bb case, such a situation only occurs at much higher c.m. energies, above 500 TeV.
The most energetic hadronic collisions observed in nature occur in collisions of O(10 20 eV) cosmic rays, at the socalled "GZK cutoff" [45], with N and O nuclei at rest in the upper atmosphere. To study the amount of triple heavyquark production produced in such collisions at equivalent c.m. energies of √ s nn ≈ 430 TeV, we show in Fig. 2 similar curves as those in Fig. 1 [33]. Around the GZK cutoff, the cross sections for inclusive as well as TPS charm production equal the total inelastic proton-air cross section, σ pAir ≈ 0.61 b, indicating that the average number of cc-pairs produced in p-Air collisions is larger than one. In the bb case, about 20% of the p-Air collisions produce bottom hadrons but only about 4% of them have TPS production. These results are clearly of relevance for the hadronic models commonly used for the simulation of the interaction of ultrarelativistic cosmic rays with the atmosphere [46] which, so far, do not include any heavy-quark production. Given that the measurements [47] observe unexplained excesses in the number of muons compared to the model predictions, and that charmed and bottom mesons feed less the non-muonic component of the air-shower, it is worth to explore the impact of properly including all such multiple heavy-quark production in the MC generators used in high-energy cosmic-ray physics.

Summary
We have derived for the first time estimates of the cross sections for triple-parton scattering (TPS) cross sections in proton-nucleus collisions as a function of the corresponding single-parton cross sections and an effective σ eff,tps,pA parameter characterizing the transverse densities of partons in the proton and nucleus. Using NNLO predictions for single heavy-quark production, we have shown that three cc-pairs are produced from separate parton interactions in ∼10% of the pPb events at the LHC. At FCC energies, more rare processes such as triple-J/ψ and triple-bb production have cross sections reaching the 1-10 mb range. At even higher energies, of a few hundred TeV reachable in the highest-energy collisions of cosmic rays with the nuclei in the atmosphere, the average number of cc-pairs produced in p-Air collisions is larger than one. The quantitative results presented here are of relevance for a proper description and understanding of final states with multiple hard particles in heavy-ion collider and cosmic-ray physics at very high energies.