# Quark jets scattering from a gluon field: from saturation to high

###### Abstract

We continue our studies of possible generalization of the Color Glass Condensate (CGC) effective theory of high energy QCD to include the high (or equivalently large ) QCD dynamics as proposed in [1]. Here we consider scattering of a quark from both the small and large gluon degrees of freedom in a proton or nucleus target and derive the full scattering amplitude by including the interactions between the small and large gluons of the target. We thus generalize the standard eikonal approximation for parton scattering which can now be deflected by a large angle (and therefore have large ) and also lose a significant fraction of its longitudinal momentum (unlike the eikonal approximation). The corresponding production cross section can thus serve as the starting point toward derivation of a general evolution equation that would contain DGLAP evolution equation at large and the JIMWLK evolution equation at small . This amplitude can also be used to construct the quark Feynman propagator which is the first ingredient needed to generalize the Color Glass Condensate (CGC) effective theory of high energy QCD to include the high dynamics. We outline how it can be used to compute observables in the large (high ) kinematic region where the standard Color Glass Condensate formalism breaks down.

## 1 Introduction

Twist expansion and collinear factorization approach [2] to particle production in QCD is a powerful and extremely useful formalism for particle production in high energy hadronic/nuclear collisions at high . However it is not expected to be valid at high energy and/or for large nuclei where twist expansion breaks down due to high gluon density effects (gluon saturation). The Color Glass Condensate (CGC) formalism (see [3] for reviews) is an effective field theory approach to QCD at high energies which relies on the fact that at high energy (or at small ) a hadron or nucleus wave function contains many gluons, referred to as gluon saturation [4, 5], and hence is a dense many-body system which is most efficiently described via semi-classical methods [6]. The most significant aspect of CGC is perhaps the emergence of a dynamical scale, called the saturation scale , which grows with energy (or ) and hence can be semi-hard. The CGC formalism thus can be used to compute quantities such as gluon multiplicities,… which are not amenable to the standard perturbative methods. Even though applications of the CGC formalism to hadronic/nuclear processes in a limited range of kinematics at RHIC and the LHC have been quite successful [7], the CGC formalism has its shortcomings, namely it is not valid when one probes large modes of the target proton or nucleus. In case of particle production in hadronic collisions this happens when high particles are produced since and are kinematically related, . This is specially important for particle production in mid rapidity as well as the Large Hadron Collider (LHC) where due to the large center of mass energy of the collision a large range of transverse momentum becomes accessible. Furthermore, the large region will be the dominant part of kinematics covered in the proposed Electron Ion Collider (EIC), at least in the earliest stages [8]. Therefore it is desirable to devise a formalism that not only incorporates the physics of saturation but also has the correct high and/or large physics encoded.

Generalizing the CGC formalism to include high physics would have significant ramifications not only for saturation physics and the work to determine its domain of applicability, it would also enable one to describe a wide range of phenomena using the same formalism. For example, saturation physics is commonly employed to provide the initial conditions for the hydrodynamic evolution of the medium, the Quark Gluon Plasma, created in high energy heavy ion collisions but is not applicable to jet (radiative) energy loss and the interactions between high partons and the produced medium. Elastic energy loss and loss (shift) of rapidity are also not present in the current formulation of saturation physics but must be included in a more general description which could be specially significant for cold matter energy loss and broadening.

Toward this goal we proposed a more general approach in [1]. We considered a high energy quark scattering from a target proton or nucleus whereas the quark not only scatters from the small gluons of the target, represented by a soft color field, but also from the large gluons in the target. We resummed the multiple scatterings of the quark from the soft classical fields to all orders in the number of soft scatterings but kept only the first scattering from the large modes. However we did not consider the large modes themselves scattering from the small gluons. Here we continue this approach and proceed to include interactions between the large gluons and the soft background field representing small gluons of the target.

We start by a brief overview of the approximations used in high energy (eikonal) scattering and how it is used in saturation physics when viewed in the target rest frame. We then give a brief summary of the approximations and the methods used in our approach in [1]. We again consider scattering of a quark from a proton or nucleus target including both small and large gluon modes. We then proceed to calculate and resum multiple interactions between the small color fields and the large gluons in the target. We then briefly outline how the calculated scattering amplitude may be used to extract the quark propagator in this more general setting and how the quark propagator may be used to compute physical observables in full range of (and/or ).

## 2 Eikonal approximation, multiple scattering at small and beyond

Here we remind the reader of the approximations involved in high energy (eikonal) scattering. As this is standard and already covered in detail in [1] we will be brief here. We define the light cone coordinates as

(1) |

and similarly for momenta and fields. The small gluons of the target are modeled as a soft color (background) field . One can either work in the frame where both the projectile and the target are moving fast, or in the frame where the projectile quark is fast and the target is at rest. In either frame the projectile quark moving to the right (along the direction to be specific) will have a large component of momentum and will couple to the conjugate component of the target color field . Furthermore, the target color field is independent of so that . We also define a light-like vector

(2) |

with which can be used to extract the Lorentz index of the soft color field and express it as so that which will help keep the expressions compact.

The -th order of the scattering of a quark, with momentum , from the color field of the target is depicted in Fig. 1 (target is shown as an ellipse) where label the coordinate positions of the field in the target (one should think of this as the projectile quark multiply scattering while going right through the target so that there are no propagators between the quark line and the target field). Diagrams of this type resum into a path-ordered infinite Wilson line provided one neglects the transverse momenta of the intermediate quark lines and the phases one picks up after integrating over the of each intermediate quark propagator . The integration over the minus component of the intermediate propagators forces a path ordering such that the scattering is sequential along the longitudinal direction, i.e., and so on (see [1, 9, 10] for details).

The scattering amplitude can then be written as

(3) |

where the infinite Wilson line is defined as

(4) |

and depicted in Fig. 2.

This infinite Wilson line resums multiple scatterings of a high energy quark (moving along the positive axis) on a soft background color field to all orders in the soft field . Due to the eikonal approximation (see also [11] where the first energy suppressed terms are investigated) the transverse position of the quark does not change during the scattering, i.e., the projectile quark does not get a significant deflection (small angle scattering). Color matrices are in the fundamental representation and soft color field represents small gluon modes of the target.

In [1] we went beyond eikonal approximation by including scattering from a large gluon field denoted which, unlike the soft field , carries large longitudinal momentum and can therefore cause a large deflection of the projectile quark. Due to the possibility of this large angle deflection (so that the final state quark has a large transverse momentum) it was necessary to introduce a rotated frame, denoted bar-ed frame, where the scattered quark is moving along a new longitudinal direction . The bar-ed coordinates are related to the original coordinates (projectile quark is moving along the direction) via the rotation matrix in dimensions.

(5) |

The elements of this -d rotation matrix are expressed in terms of the -momentum of the scattered quark. We also defined a new light cone vector which projects out the plus component in the new frame, so that . We then managed to resum all the multiple scatterings of the projectile quark from the soft field and one scattering from the large (sometimes referred to as the hard field where hard refers to a large longitudinal momentum). This is shown in Fig. 3 (target is not explicitly drawn),

Diagrams of this type resum into [1]

(6) | |||||

with ,
and the semi-infinite, anti path-ordered Wilson lines in the fundamental representation are now defined
as ^{1}^{1}1See [1] regarding the use of path vs. anti path-ordered
Wilson lines in propagators vs. amplitudes.

(7) |

and

(8) |

where anti path-ordering (AP) in the amplitude means fields with the largest argument appear to the left.

### 2.1 Multiple scatterings of the large gluon

We now proceed to consider interactions of the hard (large ) gluon with the soft background field. First let us consider the case when only the hard gluon interacts with the soft field but not the initial or final state quark, as shown in Fig. 4

The amplitude can be written as

where the free gluon propagator is

(9) |

and the triple gluon vertex is

(10) |

The large gluon field is denoted while the soft field is . It is straightforward to simplify the Lorentz structure of the amplitude by repeated use of the gauge condition as well as the null condition and the fact that the soft field has only a component and therefore is proportional to . We get

(11) | |||||

As before the soft field is independent of the coordinate which allows us to do the integration over coordinate which gives which in turn is used to do the integration setting upon integrating . We also note that the only dependence left is in the phase factors since kills all the other possible dependent terms. Then one can carry out the and integrations which lead to delta functions which are then used to perform the remaining integrations over setting and . In other words the soft field and the large (hard) field are at the same space coordinate except that the soft field does not depend on , unlike the hard field . After performing the above integrations we are left with

(12) | |||||

with the most important point being the soft and hard field are at the same point. We now go ahead and consider one more soft scattering of the hard gluon as shown in Fig. 5,

We can now repeat the same steps as before, using the gauge choice and the null vector condition to simplify the Lorentz structure, performing the integration over puts the fields and at the same point () as before, however there is now an integration over . Let us now consider the integration,

(13) |

This integration can be done using the standard contour integration techniques realizing that the pole is always below the real axis since . This integral is then

(14) |

The most essential point here is the reappearance of the theta function which forces a path ordering of the soft scatterings. This is not surprising since the soft multiple scatterings are eikonal. The rest of the analysis goes through as before and we get

(15) | |||||

Including more soft scattering on the hard gluon line is straightforward and proceeds as usual, with each extra soft scattering path ordered. This allows one to resum all the soft scatterings of the hard gluon and write it as

(16) | |||||

where the derivative acts on the coordinate of the adjoint Wilson line (anti path-ordered in the amplitude)
and arises from the fact that
one can write the soft field at the last scattering point as a derivative on the Wilson line.
This amplitude is symbolically shown in
Fig. 6 where the thick solid line attached to the hard gluon line depicts a
semi-infinite (and anti path-ordered) Wilson line in the adjoint representation
(analog of eq. (8)) ^{2}^{2}2Recall that Adjoint representation is real so
that ..

Finally we note that this amplitude vanishes in the soft (eikonal) limit, i.e. when and .

### 2.2 Multiple scatterings of the large gluon and the final state quark

We now consider the case when both the large gluon and the final state quark multiply scatter from the soft background field. The first diagram not included so far is when both the large gluon and the final state quark scatter once as shown in Fig. 7.

This scattering amplitude is given by

Many of the steps involved in simplifying this expression are identical to the previous ones, for example we use the fact that the soft field is independent of the coordinate to perform the integration over the minus components of their coordinates leading to delta functions relating the components of the momenta. The Lorentz structure can also be simplifies as before by repeated use of the gauge condition as well as the fact that and that the soft fields (in their respective frames) have only minus components which allows extraction of their Lorentz index by the use of the light-like vectors and so that and . As in the case of only the hard gluon scattering we just considered one gets and . The most important part is the integration over momentum of the intermediate quark line which can be written as

(17) |

keeping in mind that both and that we see that the integral above has two poles which are on the opposite side of the real axis. This is completely different from the eikonal scattering where the intermediate quark propagators have poles which are all on the same side of the real axis which leads to path ordering along the direction. This integral can be evaluated using the standard contour integration techniques and gives

(18) |

To proceed further and to stay consistent with the approximations made for strict eikonal scattering where one neglects terms of the order we will ignore the phase factors above. We then see that the two different path orderings corresponding to the two theta functions add to unity and path ordering disappears. This can be understood as the following, integration over any of the poles forces the other propagator to go off-shell and to become space-like in which case there is no absolute ordering between the interaction vertices at and . However, if we consider further soft scatterings of the hard gluon and the final state quark they will be path ordered with respect to and respectively. This is straightforward but long and we will just quote the final result as the calculations proceeds as earlier and there are no further subtle points. Resuming all the soft scatterings of the hard gluon and the final state quark then gives

(19) | |||||

and is depicted in Fig. 8 below where the thick solid lines denote semi-infinite (and anti path-ordered) Wilson lines in fundamental (attached to the final state quark line) and adjoint (attached to the hard gluon line) representations. Also, we have . We again note that this amplitude also vanishes in the soft limit.

### 2.3 Multiple scatterings of the initial state quark and the large gluon

We now consider the next class of diagrams in which both the large gluon and the initial state quark scatter from the soft background field. The lowest order diagram, not included so far, is shown in Fig. 9 below

This amplitude can be written as

Again most of the steps are identical to before; one integrates over component of the soft fields coordinates, eventually setting some components of momenta equal to each other. Lorentz structure is simplified using the gauge condition and the null vector as well as extracting the Lorentz index of the soft field . Integration over again sets . The crucial step now is the integration over which we focus on,

(20) |

The two poles are now both below the real axis so that the integration gives a non-zero value only if , we get

(21) |

now, unlike in the case of scattering from the large gluon and the final state quark, the relative sign between the two phase factors is negative (recall that the poles were on the opposite side in that case). Hence ignoring the phases again consistent with strict eikonal approximation, one gets a cancellation between the two terms so that , therefore this amplitude identically vanishes! It is straightforward to include any number of further soft scatterings from the initial or final state quark or the large gluon. However, it can be shown these further scatterings do not affect this null result. Therefore we conclude that one can not have simultaneous soft scatterings from the initial state quark and the large gluon. There are no other diagrams to consider, therefore this completes our derivation of the amplitude for scattering of a quark from the small and large gluon fields of the target.

The result for the full scattering amplitude at all (or any ) can thus be written as

(22) |

where , , and are given by eqs. (3,6,16,19) respectively. This is our main result. We further recall [1] that one needs to use the covariant coupling in the large terms so that the above amplitude smoothly reduces to the eikonal amplitude in the soft scattering limit, i.e. when . To extract the quark propagator one defines the effective vertex as

(23) |

in terms of which the propagator can be written as

(24) |

A final remark is in order here, we have treated the target as consisting of gluons only and have totally ignored the contribution of quarks at large . We expect sea quarks to appear when one performs a one loop correction to our result and can therefore be included in principle. On the other hand inclusion of valence quarks at large is an open problem at this point and will require a detailed investigation which is beyond the scope of this work.

## 3 Discussion and summary

We have derived the amplitude for scattering of a high energy quark on the gluon field of a proton or nucleus target including both small and large gluon modes of the target. This generalizes the standard expressions for eikonal scattering and is, to the best of our knowledge, the first gluon saturation-based calculation which includes large gluons (in the target proton or nucleus) as dynamical degrees of freedom. As such it allows one to investigate many important high and/or large phenomena which are not accessible to the standard gluon saturation formalism. This is specially essential for a proper and quantitative understanding of the outcome of the experiments at the proposed Electron Ion Collider and the Large Hadron Collider.

The derived scattering amplitude is a ”tree-level” expression which can already be used to investigate several phenomena, for example broadening and elastic energy loss as well as the nuclear modification factor [12]. Due to the presence of large gluon fields in the target the scattered quark can undergo an arbitrarily large deflection and pick up large transverse momenta. Furthermore, it can lose longitudinal momentum and undergo a potentially large rapidity loss which is not contained in saturation formalism. One can also extract the quark propagator from this scattering amplitude and use it to calculate the tree-level production cross sections for particle production in high energy collisions for any transverse momentum (at small or large ). One would also need to calculate the gluon propagator [13] in this formalism which is a straightforward extension of the present work. This would allow one to investigate cold matter radiative energy loss including both the fully coherent, present in saturation formalism, as well as the partially coherent/incoherent energy loss which is not present in the saturation formalism. Examples of where our results in the present form can be used are single inclusive or di-jet production in DIS as well as in high energy proton-proton and proton-nucleus collisions. The present work generalizes the saturation formalism and enlarges the transverse momentum range (recall and are kinematically related) where gluon saturation-based models are applied and improve their quantitative accuracy [14]. In addition one would also be able to investigate forward-backward (in rapidity) correlations in our framework.

Naturally one expects that our tree-level expression will be renormalized when one considers radiative (one-loop) corrections. In analogy with renormalization of product of Wilson lines [14, 15, 16] in small QCD which leads to the JIMWLK/BK evolution equation [17, 18], we expect the renormalization of the scattering cross section here to lead to a more general evolution equation which incorporates both the DGLAP [19] and JIMWLK evolution equations; due to the presence of the large gluon field (not present in saturation formalism) one expects the one loop corrections to result in the DGLAP evolution equation in the large limit. On the other hand and due to presence of the eikonal term one would expect to recover the JIMWLK evolution equation in the small limit. Therefore it will be enormously beneficial to calculate the one loop corrections to our result. It may also be possible to reformulate this as an effective action approach, analogous to the McLerran-Venugopalan model [6]. If so this would make it possible to treat both the early stages in the formation of a Quark Gluon Plasma and the high jet energy loss in high energy heavy ion collisions using the same formalism, at least in the earliest times after the collision [20]. In summary, the present work takes the first step toward deriving a formalism that generalizes the Color Glass Condensate framework by including the physics of high and large .

## Acknowledgments

We acknowledge support by the DOE Office of Nuclear Physics through Grant No. DE-FG02-09ER41620 and by the Idex Paris-Saclay though a Jean d’Alembert grant. We would like to thank the staff of Pedro Pascual Science Center in Benasque, Spain for their kind hospitality during the completion of this work. We also thank T. Altinoluk, N. Armesto, F. Gelis, E. Iancu, Yu. Kovchegov, A. Kovner, C. Lorcé, C. Marquet, A.H. Mueller, S. Munier, B. Pire, C. Salgado, G. Soyez, R. Venugopalan, D. Wertepny and B. Xiao for critical questions, illuminating discussions and helpful suggestions.

## References

- [1] J. Jalilian-Marian, Phys. Rev. D 96, no. 7, 074020 (2017) doi:10.1103/PhysRevD.96.074020 [arXiv:1708.07533 [hep-ph]].
- [2] R. Brock et al. [CTEQ Collaboration], Rev. Mod. Phys. 67, 157 (1995). doi:10.1103/RevModPhys.67.157
- [3] E. Iancu and R. Venugopalan, In *Hwa, R.C. (ed.) et al.: Quark gluon plasma* 249-3363 doi:10.1142/97898127955330005 [hep-ph/0303204]. F. Gelis, E. Iancu, J. Jalilian-Marian and R. Venugopalan, Ann. Rev. Nucl. Part. Sci. 60, 463 (2010) doi:10.1146/annurev.nucl.010909.083629 [arXiv:1002.0333 [hep-ph]]. J. Jalilian-Marian and Y. V. Kovchegov, Prog. Part. Nucl. Phys. 56, 104 (2006) doi:10.1016/j.ppnp.2005.07.002 [hep-ph/0505052]. H. Weigert, Prog. Part. Nucl. Phys. 55, 461 (2005) doi:10.1016/j.ppnp.2005.01.029 [hep-ph/0501087].
- [4] L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rept. 100, 1 (1983). doi:10.1016/0370-1573(83)90022-4
- [5] A. H. Mueller and J. w. Qiu, Nucl. Phys. B 268, 427 (1986). doi:10.1016/0550-3213(86)90164-1
- [6] L.D. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994); ibid. 3352 (1994); ibid. 50, 2225 (1994); S. Jeon and R. Venugopalan, Phys. Rev. D 71, 125003 (2005) doi:10.1103/PhysRevD.71.125003 [hep-ph/0503219]; A. Dumitru, J. Jalilian-Marian and E. Petreska, Phys. Rev. D 84, 014018 (2011) doi:10.1103/PhysRevD.84.014018 [arXiv:1105.4155 [hep-ph]].
- [7] T. Lappi, EPJ Web Conf. 137, 07013 (2017) doi:10.1051/epjconf/201713707013 [arXiv:1611.07668 [hep-ph]].
- [8] E. C. Aschenauer et al., arXiv:1602.03922 [nucl-ex].
- [9] A. J. Baltz, F. Gelis, L. D. McLerran and A. Peshier, Nucl. Phys. A 695, 395 (2001) doi:10.1016/S0375-9474(01)01109-5 [nucl-th/0101024]; F. Gelis and A. Peshier, Nucl. Phys. A 697, 879 (2002) doi:10.1016/S0375-9474(01)01264-7 [hep-ph/0107142].
- [10] A. Kovner and U. A. Wiedemann, In Hwa, R.C. (ed.) et al.: Quark gluon plasma 192-248 doi:10.1142/97898127955330004 [hep-ph/0304151]; J. Casalderrey-Solana and C. A. Salgado, Acta Phys. Polon. B 38, 3731 (2007) [arXiv:0712.3443 [hep-ph]].
- [11] T. Altinoluk and A. Dumitru, Phys. Rev. D 94, no. 7, 074032 (2016) doi:10.1103/PhysRevD.94.074032 [arXiv:1512.00279 [hep-ph]]; T. Altinoluk, N. Armesto, G. Beuf and A. Moscoso, JHEP 1601, 114 (2016) doi:10.1007/JHEP01(2016)114 [arXiv:1505.01400 [hep-ph]]; T. Altinoluk, N. Armesto, G. Beuf, M. Martínez and C. A. Salgado, JHEP 1407, 068 (2014) doi:10.1007/JHEP07(2014)068 [arXiv:1404.2219 [hep-ph]]; Y. V. Kovchegov and M. D. Sievert, arXiv:1808.09010 [hep-ph].
- [12] J. Jalilian-Marian, Nucl. Phys. A 748, 664 (2005) doi:10.1016/j.nuclphysa.2004.12.001 [nucl-th/0402080]. J. Jalilian-Marian, Y. Nara and R. Venugopalan, Phys. Lett. B 577, 54 (2003) doi:10.1016/j.physletb.2003.09.097 [nucl-th/0307022].
- [13] A. Ayala, J. Jalilian-Marian, L. D. McLerran and R. Venugopalan, Phys. Rev. D 52, 2935 (1995) doi:10.1103/PhysRevD.52.2935 [hep-ph/9501324].
- [14] G. A. Chirilli, B. W. Xiao and F. Yuan, Phys. Rev. D 86, 054005 (2012) doi:10.1103/PhysRevD.86.054005 [arXiv:1203.6139 [hep-ph]], Phys. Rev. Lett. 108, 122301 (2012) doi:10.1103/PhysRevLett.108.122301 [arXiv:1112.1061 [hep-ph]]; B. W. Xiao and F. Yuan, arXiv:1407.6314 [hep-ph]; Z. B. Kang, I. Vitev and H. Xing, Phys. Rev. Lett. 113, 062002 (2014) doi:10.1103/PhysRevLett.113.062002 [arXiv:1403.5221 [hep-ph]]; E. Iancu, A. H. Mueller and D. N. Triantafyllopoulos, JHEP 1612, 041 (2016) doi:10.1007/JHEP12(2016)041 [arXiv:1608.05293 [hep-ph]]; A. M. Stasto and D. Zaslavsky, Int. J. Mod. Phys. A 31, no. 24, 1630039 (2016) doi:10.1142/S0217751X16300398 [arXiv:1608.02285 [hep-ph]].
- [15] A. Dumitru and J. Jalilian-Marian, Phys. Rev. Lett. 89, 022301 (2002) doi:10.1103/PhysRevLett.89.022301 [hep-ph/0204028], Phys. Lett. B 547, 15 (2002) doi:10.1016/S0370-2693(02)02709-0 [hep-ph/0111357]; A. Dumitru, A. Hayashigaki and J. Jalilian-Marian, Nucl. Phys. A 765, 464 (2006) doi:10.1016/j.nuclphysa.2005.11.014 [hep-ph/0506308], Nucl. Phys. A 770, 57 (2006) doi:10.1016/j.nuclphysa.2006.02.009 [hep-ph/0512129].
- [16] A. Dumitru, J. Jalilian-Marian, T. Lappi, B. Schenke and R. Venugopalan, Phys. Lett. B 706, 219 (2011) doi:10.1016/j.physletb.2011.11.002 [arXiv:1108.4764 [hep-ph]]. A. Dumitru and J. Jalilian-Marian, Phys. Rev. D 82, 074023 (2010) doi:10.1103/PhysRevD.82.074023 [arXiv:1008.0480 [hep-ph]].
- [17] J. Jalilian-Marian, A. Kovner, L. D. McLerran and H. Weigert, Phys. Rev. D 55, 5414 (1997); J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert, Nucl. Phys. B 504, 415 (1997), Phys. Rev. D 59, 014014 (1999), Phys. Rev. D 59, 014015 (1999), Phys. Rev. D 59, 034007 (1999), A. Kovner, J. G. Milhano and H. Weigert, Phys. Rev. D 62, 114005 (2000); A. Kovner and J. G. Milhano, Phys. Rev. D 61, 014012 (2000); E. Iancu, A. Leonidov and L. D. McLerran, Nucl. Phys. A 692, 583 (2001), Phys. Lett. B 510, 133 (2001); E. Ferreiro, E. Iancu, A. Leonidov and L. McLerran, Nucl. Phys. A 703, 489 (2002).
- [18] I. Balitsky, Nucl. Phys. B 463, 99 (1996); Yu.V. Kovchegov, Phys. Rev. D 61, 074018 (2000).
- [19] G. Altarelli and G. Parisi, Nucl. Phys. B 126, 298 (1977); V.N. Gribov, L.N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972); ibid. 675 (1972); Yu. Dokshitzer, Sov. Phys. JETP 46, 641 (1977).
- [20] A. Kovner, L. D. McLerran and H. Weigert, Phys. Rev. D 52, 6231 (1995) doi:10.1103/PhysRevD.52.6231 [hep-ph/9502289]; A. Kovner, L. D. McLerran and H. Weigert, Phys. Rev. D 52, 3809 (1995) doi:10.1103/PhysRevD.52.3809 [hep-ph/9505320].