# Nucleon 3D structure from double parton scattering: a Light-Front quark model analysis

###### Abstract

Double parton distribution functions (dPDFs) represent a tool to explore the 3D partonic structure of the proton. They can be measured in high energy proton-proton and proton nucleus collisions and encode information on how partons inside a proton are correlated among each other. dPFDs are studied here in the valence region by means of a constituent quark model scenario within the relativistic Light Front approach, where two particle correlations are present without any additional prescription. Furthermore, a study of the QCD evolution at high energy scale, of the model results, has been completed in order to compare our predictions with future data analyses. In closing, results on the evaluation of the so called , crucial ingredient for the description of double parton scattering, where dPDFs can be accessed, are presented and discussed.

## 1 Introduction

The observation of multiple parton interactions (MPI), occurring in high energy hadron-hadron collisions, will be an important goal of the experimental and theoretical analyses of the processes studied, e.g., at the LHC. In particular, in these kind of events, more than one parton of a hadron interact with partons of the other colliding hadron. Naturally, MPI contribution to the total cross section is suppressed with respect to the single parton interaction. However, the measurement and estimate of MPI cross sections is an important challenge being MPI events a background for the search of new Physics. Furthermore, from a hadronic point of view, we are interested on MPI due to the possibility of accessing new fundamental information on the partonic proton structure. We focus our studies on the the double parton scattering (DPS), the most simple case of MPI, which can be observed, in principle, in several processes, e.g., with dilepton productions and double Drell-Yan processes (see, Refs. [1, 2, 3, 4, 5] for recent reviews). At the LHC, DPS, has been observed some years ago [6] and represents a background for the Higgs production in several channels. From a theoretical point of view, the DPS cross section is written in terms of a new quantity, the so called double parton distribution function (dPDF), , which describe the joint probability of finding two partons of flavors with longitudinal momentum fractions and separation in the transverse plane inside the hadron, see Ref. [7]. Here is the renormalization scale. However, to date, being dPDFs poorly known, in order to qualitatively estimate the magnitudo of the DPS cross section, the following approximation is usually adopted:

(1) |

i.e., dPDFs are evaluated in a fully factorized ansatz in terms of standard one-body parton distribution functions (PDF), , and , the function encoding all parton correlations in the transverse plane. The spirit of the latter assumption is to neglect all possible double parton correlations (DPCs) between the two interacting partons, being the latter almost unknown, see e.g. Ref. [8] for updates. Moreover dPDFs are non perturbative quantities in QCD so that they cannot be easily evaluated from the theory. Nevertheless, as already deeply discussed in Refs. [9, 10, 11, 12], quark model calculations of dPDFs could help to grasp the basic feature of such quantities. This strategy has been largely used in the past to study unknown distributions. However, since in this scenario dPDFs are calculated at the low hadronic scale of the model, , in order to compare the obtained outcomes with future data taken at high energy scales, , it is then necessary to perform the perturbative QCD (pQCD) evolution of the model calculations, using the dPDF evolution equations, see Refs. [13, 14]. The idea supporting our analysis is that, thanks to this procedure, future data analyses of the DPS processes could be guided, in principle, by model calculations. In the first part of the present paper, we focus our attention on the study of the role of DPCs in dPDFs in order to verify, as a first step, the validity of the approximation Eq. (1), often used for data analyses. In particular, it is here recall that in all model calculations of dPDFs, see Refs. [9, 10, 11], the assumption Eq. (1) is violated. See also Ref. [15] for details on the violation of the factorized ansatz due to model independent relativistic effects in the calculated dPDFs. Moreover, as already pointed out, it is fundamental to realize to which extend DPCs survive at very high energy scales, where, due to the large population of partons, the role of DPCs could be less important than at the low scale of the used model. To this aim in Refs. [11, 16], DPCs in dPDFs have been studied at the energy scale of the experiments. In particular in Ref. [16] an extension of the approach used in Refs. [11] has been provided to include sea quarks and gluon degrees of freedom. Here, the most important results of the these analyses will be summarized. In the last section of the present paper, an investigation of the so called will be discussed. Usually, indeed, DPS cross section, in processes with final state , is written through the following ratio (see e.g. Ref. [17]):

(2) |

where is combinatorial factor depending on the final states and ( for or for ) and is the single parton scattering cross section with final state . Expressing the cross section in Eq. (2) in terms of product of , one assumes that, as a first approximation, double parton distributions can be written as in Eq. (1). The present knowledge on DPS cross sections has been condensed in the experimental and model dependent extraction of [17, 18, 19, 20, 21, 22, 23, 24]. To date, 15 mb, compatible, within errors, with a constant, irrespective of centre-of-mass energy of the hadronic collisions, the final state and . Let us stress again that the latter result is obtained considering the fully uncorrelated scenario shown in Eq. (1) where all DPCs are neglected. Also in this case non perturbative methods have been used to study , through the calculation of dPDFs, in order to characterize “signals” of DPCs in . On top of that, in Ref. [25], the possible dependence of the latter quantity on the longitudinal momentum fraction carried by the interacting partons has been addressed. Recent results on this topic will be here discussed, including also new intriguing outcomes addressed in Refs. [15, 26]. In Ref. [15] the contribution of relativistic effects to have been estimated by means of the Light-Front approach (LF), already used for the evaluation of dPDFs in Ref. [11]. Furthermore, in Ref. [26], has been calculated in a AdS/QCD framework, showing, also here, a strong dependence of in the valence region.

## 2 Calculation of dPDFs within LF CQMs

In this section, details of the calculations of dPDFs, within the LF approach (see Refs. [27, 28] for fundamental details), will be discussed verifying if the factorized ansatz Eq. (1) would work in this relativistic scenario. In this framework, among all positive features of the LF approach, let us remark that one can obtain a fully Poincaré covariant description of relativistic strongly interacting systems with a fixed number of on-shell constituents, LF boosts and the “plus” components of momenta () are kinematical operators. Furthermore, being the LF hypersurface, the one where the initial conditions of the system are fixed, tangent to the light-cone, the kinematics of DIS processes is naturally obtained. For these reasons, this strategy has been adopted in the past for the evaluation unknown distributions, see e.g. Refs. [29, 30, 31]. All the details on the present calculations are deeply discussed and reported in Ref. [11] and will not be repeated here. The final expression of the dPDF in momentum space for two unpolarized quarks of flavor , reads:

being the relative transverse momentum of one of the parton in the amplitude and its complex conjugate and the momentum of the quark.

The canonical proton wave function is embedded in the function here above, which can be written as follows:

(3) |

where and are the canonical parton helicity and the isospin, respectively. Here the short notation is introduced. In Eq. (3), the Melosh operators, which allow to rotate the LF spin into the canonical one are introduced:

(4) |

being the total free energy mass of the system and and generic canonical spins. In this scenario, being , appearing in delta functions in Eq. (2), a kinematical quantity, dPDFs are well defined in the region, at variance with what happens in the instant form calculations of PDFs and dPDFs (see, e.g., Ref. [10]). Let us remark that the main ingredient in Eq. (2) is the canonical proton wave function which has been calculated by means of a relativistic CQM, the so called hyper-central CQM described in Ref. [32]. The choice of this model is motivated by its simplicity and capability to basically reproduce the spectrum of light-hadrons. Since actually for dPDFs there are not yet available data, model calculations of these quantities could provide essential information on their relevant features. Before discussing in the next section the results of the calculations of dPDFs, one should notice that in principle the role of the model independent Melosh rotations can be important. Indeed, from Eq. (4), one can realize that such operator con introduce, in the calculation of dPDFs, non trivial correlations between and . As deeply discussed in Ref. [15], in order to taste possible effects of these objects, it is instructive to study the following quantity:

(5) |

being the intrinsic transverse component of the quark momentum. In the latter quantity use has been made of the commonly adopted SU(6) symmetry in order to evaluate the spin part of the proton wave function. In this scenario, the quantity (5) is rather model independent. Here we consider two fast partons (FF) with , one slow and and one fast parton (SF) with and two slow partons (SS) with . The calculation of Eq. (5) in these kinematic regions is presented in Fig. (1) where one may identify three distinct regions as a function of . For the Melosh’s in all kinematic configurations reduce to unity. In an intermediate region of the curves show a dip whose depth depends on the chosen kinematic configuration and at larger the curves flattens out with different asymptotics. This complicated pattern, generated by Melosh’s rotations, affects the calculation of dPDFs, which, in general, are distributions evaluated also at .

## 3 Results of the calculations of dPDFs at the hadronic scale

In this section the main results of the calculations of dPDFs at the hadronic scale, GeV, will be presented. In particular, as already mentioned, here the emphasis of the analysis is focused on testing the validity of the approximation Eq. (1). As already discussed in the previous section, in general, Melosh rotations introduce correlations between and so that, as also found in Ref. [11], at the hadronic scale such factorization is violated in all model calculations of dPDFs within the LF approach, see Ref. [15]. Furthermore, in order to study the validity of the approximations Eq. (1), in the and dependence of dPDFs, the following ratio have been evaluated:

(6) |

As one can see on the left panel of Fig. (2), being the ratio , Eq. (6), different from the unity in all the kinematical range, also the factorized form of dPDFs in terms of the product of single PDFs is not supported by the present approach. In closing, at the hadronic scale, in all model calculations of dPDFs (see also Refs. [9, 10]), the assumption Eq. (1) is violated due to DPCs.

## 4 pQCD evolution of the calculated dPDFs

A fundamental point, discussed e.g. in Ref. [11], is the analysis of the pQCD evolution of model calculations of dPDFs. This procedure is essential to relate CQM predictions with present and future experimental analyses. For the moment being, the pQCD evolution of dPDFs is known only for the longitudinal momentum dependence, which means , and using the same energy scale for both the acting partons. In these case, the evolution equations are obtained as a proper generalization of the usual Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) ones (see Refs. [13, 14] for details), defined for the evolution of PDFs. In the first part of this section, results will be shown for the non-singlet sector. In particular, the ratio , Eq. (6), has been calculated using the dPDFs evaluated at a generic high energy scale, e.g., GeV, using pQCD evolution. As one can see on the right panel of in Fig. 2, for small values of , e.g, close to the LHC kinematics, . This means that dynamical correlations, for valence quarks, are suppressed after the evolution. Nevertheless, in order to complete this analysis, in Ref. [16], new ratios sensitive to longitudinal correlations have been defined and calculated at high energy scale including perturbative and non perturbative degrees of freedom. To this aim let us generalize the expression of for different partonic species:

(7) |

where here , and are the single PDFs for two given partons of flavor and and GeV. In such ratio, dPDFs in the numerator and PDFs in the denominator evolve by means of different evolution equations. Due to this feature, is sensitive to non perturbative correlations, encoded in the proton wave function used to calculate dPDFs and PDFs and to the perturbative ones due to the difference in the evolution equations of dPDFs and PDFs. In order to disentangle these two different contributions to understand which among theme could affect the dPDF evaluations, the following ratios have been defined and calculated within the LF CQM approach:

(8) |

where here:

(9) |

The latter quantity is calculated by evolving the product of PDFs by means of dPDF evolution. Due to this feature, the numerator and denominator of evolve differently in pQCD while the input at the initial scale is the same, so that is sensitive to perturbative correlations. On the contrary, since in ratio the numerator and denominator evolve within the same scheme, but using different input at the initial scale, it is sensitive to non perturbative correlations. Due to these features, the ratio , shown in Fig. 3 (left panel) for two , is particularly emblematic. The full ratio , Eq. (7), (dashed line), clearly influenced by both perturbative (dot-dashed line) and non-perturbative (continuous line) effects, is compared with those where perturbative and non-perturbative correlations are disentangled, contributing to the behavior of dPDFs at low values of and . As one can see in the left panel of Fig. 3, for the gluon-gluon distribution such components coherently interfere. Here we report results only for gluons, being the highest partonic component at LHC kinematics. Furthermore, in order to show how much these conclusions do not depend on the used model, a semi factorized ansatz has been adopted in order to include non perturbative see quarks and gluons at a given initial scale. All details of the procedure are discussed in Ref. [16]. We assume that at a given initial scale one has:

(10) | ||||

where in the above expression, and are the dPDFs and PDFs respectively calculated by means of the LF CQM approach obtained evolving the same quantities from to and is the PDF evaluated through to the MSTW2008 parametrization [33], see Ref. [16] for further details. In particular, within this choice, GeV. As one can see in the right panel of Fig. 3, where has been calculated starting the evolution from and , conclusions arisen from the precedent analysis, see left panel of Fig. 3, are confirmed. For gluon-gluon the correlations induced at low- still contain a specific sign of the correlations introduced in the valence sector and this is due to the presence of the valence component in the quark-singlet sector in the evolution procedure. The strength of the correlations seems to become smaller but they are still sizable and should not be neglected.

## 5 Calculation of

As previously mentioned, an important quantity, relevant for the experimental analyses of DPS is the effective cross section, whose expression in terms of PDFs and dPDFs, has been presented in Ref. [25]. In order to discuss the main features of , we restrict the analysis to the zero rapidity region and show results, within the LF approach, in the left panel of Fig. 4. Here the latter quantity has been calculated at GeV for a sea and a valence quarks. It is worth to notice that the three old experimental extractions of from data, Refs. [34, 35, 36], lie in the obtained range of values of the calculated . It is fundamental to remark that, at the hadronic scale, as demonstrated in Ref. [15], if Melosh rotations would neglected, the average value of would change by a factor 2, making the impact of relativistic effects in very strong. In the same figure, as one can see, DPCs generate a strong dependence of in the valence region. Such feature, as already discussed, is related to the combined effect of pQCD evolution correlations and non perturbative ones. In order to provide a fully model independent analysis of the dependence of , in Ref. [26], such quantity has been calculated through an evaluation of dPDFs based on the AdS/QCD correspondence, see Refs. [37, 38, 39] for fundamental details. In particular the authors found that, at the hadronic scale, the dependence of is comparable with the one found within the LF approach. Moreover, in order to emphasize such feature, in the right panel of Fig. 4, the ratio has been plotted. As one can see, such ratio strongly depends on in the valence region, at the variance of the case in which longitudinal correlations were neglected.

## 6 Conclusions

In this work, dPDFs have been calculated within a fully Poincaré covariant constituent quark model approach, the Light-Front one, fulfilling the essential properties of dPDFs, see Ref. [11]. The main goal of this first analysis was to establish the role of DPCs in dPDFs. In particular, the factorization of dPDFs, at the hadronic scale of the model, in the and dependences is violated, as already found in previous analyses c.f. Refs. [9, 10]. In order to compare our results with actual and future experimental studies, the pQCD evolution of the calculated dPDFs has been deeply investigated as discussed in Refs. [11, 16] where, at very high energy scales, like the experimental ones, it was found that longitudinal correlations may survive also at small values of . Furthermore, as also explained in Ref. [16], where also non perturbative degrees of freedom have been taken into account in the analysis, there are different source of longitudinal correlations, the perturbative one, induced by the pQCD evolution scheme and the non perturbative ones due to the model used to describe the proton wave function, used to calculate the dPDFs. For example, in the gluon gluon case, such correlations coherently interfere in the small region at very high energy scales making theme sizable. However, since dPDFs can not be directly measured in DPS, within this model calculation approach, the so called , a quantity experimentally studied to estimate the DPS cross section, has been calculated. The latter one could be very important from the theoretical and experimental point of view to obtain information on partonic structure of the proton being sensitive to DPCs. In particular, signals of DPCs in , have been found in its dependence in the valence region also at high energy scales, for details see Ref. [25]. Furthermore has been also calculated within the AdS/QCD correspondence and also in this case, an important dependence has been found, in particular in the valence region. Since, the role of correlations is crucial to grasp new information on the partonic structure of the proton, relativistic effects in the calculation of dPDFs have been studied in order to identify model independent correlations in the dependence of dPDFs. In particular, as deeply discussed in Ref. [15], the role of Melosh rotations, essential feature of the LF approach, introduce non negligible correlations in such dependence. Furthermore new studies for the extraction of the proton dPDFs from collisions are under investigations, see Ref. [40]. We conclude that the investigation of dPDFs, also trough the analysis of , could unveil new and interesting information on the 3D partonic structure of the proton. I thank Sergio Scopetta, Marco Traini, Vicente Vento and Federico Alberto Ceccopieri for a nice and fruitful collaboration on this subject. This work was supported in part by the Mineco under contract FPA2013-47443-C2-1-P and SEV-2014-0398.

## References

## References

- [1] Gaunt J R and Stirling W J 2010 JHEP 03 005
- [2] Diehl M, Ostermeier D and Schafer A 2012 JHEP 1203 089
- [3] Manohar A V and Waalewijn W J 2012 Phys. Rev. D 85 114009
- [4] Bartalini P and Fanò arXiv:1003.4220 [hep-ex]
- [5] Bansal S et al., arXiv:1410.6664 [hep-ph]
- [6] Aad G et al. [ATLAS Collaboration] 2013 New J. Phys. 15 033038
- [7] Paver N and Treleani D 1982 Nuovo Cim. A 70 215
- [8] Diehl M, Kasemets T and Keane S 2014 JHEP 05 118
- [9] Chang H M, Manohar A V and Waalewijn W J 2013 Phys. Rev. D 87 034009
- [10] Rinaldi M, Scopetta S and Vento V 2013 Phys. Rev. D 87 11
- [11] Rinaldi M, Scopetta S, Traini M and Vento V 2014 JHEP 12 028
- [12] Broniowski W and Arriola E R 2014 Few Body Syst. 55 381
- [13] Kirschner R 1979 Phys. Lett. B 84 266
- [14] Shelest V P, Snigirev A M and Zinovev G M 1982 Phys. Lett. B 113 325
- [15] Rinaldi M and Ceccopieri F A 2017 Phys. Rev. D 95 no. 3 034040
- [16] Rinaldi M, Scopetta S, Traini M C and Vento V 2016 JHEP 10 063
- [17] Jung H et al. Proceedings, 7th International Workshop on Multiple Partonic Interactions at the LHC (MPI@LHC 2015) https://bib-pubdb1.desy.de/record/297386, DESY-PROC-2016-01
- [18] Akesson T et al. [Axial Field Spectrometer Collaboration] 1987 Z. Phys. C 34 163
- [19] Alitti J et al. [UA2 Collaboration] 1991 Phys. Lett. B268 145
- [20] Abe F et al. [CDF Collaboration] 1997 Phys. Rev. D 56 3811
- [21] Abazov V M et al. [D0 Collaboration] 2010 Phys. Rev. D 81 052012
- [22] Aad G et al. [ATLAS Collaboration] 2013 New J. Phys. 15 033038
- [23] Chatrchyan S et al. [CMS Collaboration] 2014 JHEP 03 032
- [24] Aaij R et al. [LHCb Collaboration] 2014 JHEP 04 091
- [25] Rinaldi M, Scopetta S, Traini M and Vento V 2016 Phys. Lett. B 752 40
- [26] Traini M, Scopetta S, Rinaldi M and Vento V 2017 Phys. Lett. B 768 270
- [27] Keister B D and Polyzou W N 1991 Adv. Nucl. Phys. 20 225
- [28] Brodsky S J, Pauli H C and Pinsky S S 1998 Phys. Rept. 301 299
- [29] Boffi S, Pasquini B and Traini M 2003 Nucl. Phys. B 649 243
- [30] Pasquini B, Traini M and Boffi S 2005 Phys. Rev. D 71 034022
- [31] Traini M 2014 Phys. Rev. D 89 no. 3 034021
- [32] Faccioli P, Traini M and Vento V 1999 Nucl. Phys. A 656 400
- [33] Martin A D, Stirling W J, Thorne R S and Watt G 2009 Eur. Phys. J. C 63 189
- [34] Akesson T et al. [Axial Field Spectrometer Collaboration] 1987 Z. Phys. C 34 163
- [35] Alitti J et al. [UA2 Collaboration] 1991 Phys. Lett. B 268 145
- [36] Abe F et al. [CDF Collaboration] 1997 Phys. Rev. D 56 3811
- [37] Maldacena J M 1998 Adv. Theor. Math. Phys. 2 231 Gubser S S, Klebanov I R and Polyakov A M 1998 Phys. Lett. B 428 105 Witten E 1998 Adv. Theor. Math. Phys. 2 253
- [38] Gursoy U and Kiritsis E 2008 JHEP 0802 032 Gursoy U, Kiritsis E and Nitti F 2008 JHEP 0802 019
- [39] Brodsky S J, de Teramond G F, Dosch H G and Erlich J 2015 Phys. Rept. 584 1
- [40] Salvini S, Treleani D and Calucci G 2014 Phys. Rev. D 89 016020