Abstract
The observation in small size collision systems, pp and pA, of strong correlations with long range in rapidity and a characteristic structure in azimuth, the ridge phenomenon, is one of the most interesting results obtained at the large hadron collider. Earlier observations of these correlations in heavy ion collisions at the relativistic heavy ion collider are standardly attributed to collective flow due to strong final state interactions, described in the framework of viscous relativistic hydrodynamics. Even though data for small size systems are well described in this framework, the applicability of hydrodynamics is less well grounded and initial state based mechanisms have been suggested to explain the ridge. In this review, we discuss particle correlations from the initial state point of view, with focus on the most recent theoretical developments.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
While the focus of the physics programme at the large hadron collider (LHC) is the discovery and understanding of the properties of the previously missing piece in the Standard Model – the Higgs boson – and the search for its eventual failure, it has also shown very surprising and unexpected aspects of quantum chromodynamics (QCD), particularly in small collisions systems, pp and pA. One of the most exciting observations made in high multiplicity pp collisions by the CMS collaboration during the first LHC run is the discovery of the correlations between produced particles over large intervals of rapidity, peaking at zero relative azimuthal angle [1]. This phenomenon was dubbed ridge due its shape in the azimuthal anglerapidity plot, and constitute one of the key findings at the LHC (see Fig. 1).
Later on, this structure was found by other collaborations and for smaller multiplicities [2,3,4,5] and in association with Z boson production [6]. A similar ridge structure was also observed in pPb collisions at the LHC by the four large collaborations [7,8,9,10]. A maximum in the correlations also appears at azimuthal angle \(\pi \), called the away side ridge in contrast to the near side ridge peaked at zero azimuthal angle. They have also been observed in PbPb collisions, see e.g. [11,12,13] for PbPb results and a comparison with those in pPb. Similar correlations were observed in AuAu, dAu and \(^3\)HeAu collisions at the relativistic heavy ion collider (RHIC) [14,15,16,17,18,19]. They have also been observed in photoproduction on Pb in ultraperipheral collisions (UPCs) at the LHC [20]. Their existence in smaller systems like \(e^+e^\) collisions [21] at the large electronpositron collider and deep inelastic scattering (DIS) events in ep at the HadronElektronRinganlage [22] has been scrutinised, but the results are not conclusive. The ridge is the most striking feature in the long list of similarities between small and large collision systems in the observed results for many observables [23,24,25,26,27].
The standard explanation for such azimuthal asymmetries in heavy ion collisions (HICs) is the existence of strong final state interactions that lead to a situation where viscous relativistic hydrodynamics can be applied, see the reviews [28, 29]. The dynamics leading to such a situation, called hydrodynamisation, is unclear [30] and both strong and weak coupling explanations have been proposed, see e.g. [31] and refs. therein. Furthermore, the hydrodynamic description seems to hold for large anisotropies, i.e. rather far from local equilibrium.
The hydrodynamic description of the azimuthal asymmetries in pp and pPb collisions at the LHC is successful [29, 32]. There is a hot ongoing discussion on the explanation of such success and of the fact that hydrodynamics seems to be the effective description for long wavelength modes in any field theory, see e.g. [33] and refs. therein. But it demands a very careful choice of initial conditions, specifically that the proton is modelled as a collection of constituent quarks or hot spots. The description seems to be pushed to the limit of small collision areas and low particle densities where nonhydrodynamic modes play a very important role, as seen in both hydrodynamic studies [29] and in those that consider a weak coupling quasiparticle picture in transport frameworks [34].
Therefore, the hydrodynamic explanation for the azimuthal asymmetries in small systems looks tenuous. Besides, causality arguments show that long range correlations in rapidity must come from the very early stages of the collision [35]. And hydrodynamic calculations demand initial conditions that contains long range rapidity correlations, initial energy and particle density and flow profiles and, unless they are assumed to be completely washed out by final state interactions, correlations. So, beyond addressing the obvious fundamental question: is the strong interaction dynamics capable to lead to collectivity through final state interactions even in small systems, or is the origin of the ridge correlations different in pp and pPb collisions than in HICs?, searching for correlations coming from the initial state is linked to the understanding of the dynamics prior to the use of hydrodynamics and to the provision of well grounded initial conditions for hydrodynamic calculations.
This contribution is devoted to the description and discussion of those frameworks which lead to correlations among partons in the initial stage that, if not washed out by strong final state interactions and hadronisation – that we will assume in the following, may lead to azimuthal asymmetries as observed in data. We start by those based on the weak coupling but nonperturbative description of dense partonic systems offered by the the Color Glass Condensate (CGC) effective theory, see the reviews [36,37,38] and the book [39]. This will be the subject of Sects. 2 and 3. We will then review other explanations inspired in QCD in Sect. 4, to end with a summary and discussions in Sect. 5.
Our focus will be on recent formal developments and we will base the presentation in our own works and formalism, trying to make connection with the other formalisms which differ in notation. We will comment briefly on the status of the comparison to experimental data in the summary.
2 Two particle correlations from the CGC
The observation of the ridge correlations in small size systems has triggered a lot of efforts to understand whether the structure of the initial state itself can lead, in pp and pA collisions, to such correlations without resourcing to final state interactions. Over the last decade, several mechanisms have been suggested to explain the ridge correlations in the CGC framework. The two most successful ones are the “domain structure of the target” developed in [40,41,42] and the “glasma graph approach” introduced in [35, 43, 44].^{Footnote 1}
The underlying mechanism for the domain structure of the target can be summarised as follows: the hadronic target is assumed to contain domains of oriented chromoelectric fields in the transverse plane. When two partons (normally assumed to be gluons when the scattering takes place at high energies and the probed values of momentum fraction of the partons, x, is small) from the projectile are close enough to scatter on the same domain, they get a common final momentum that reflects the correlated structure of the fields in the target. As gluons belong to the adjoint, thus real, representation of the \(SU(N_c)\) group, the correlation holds for both parallel and antiparallel momenta, thus justifying the near and away side structures. The size of the domain in the target is assumed to be of order \(1/Q_s\), with \(Q_s\) being the saturation momentum which is the characteristic transverse momentum for the partons in the saturated target wave function described by the CGC [36,37,38]. Projectile partons lying closer than \(1/Q_s\) contribute mainly to particle production in the region of transverse momentum \(p_T\gtrsim Q_s\). Therefore, this mechanism should mainly be applicable in that transverse momentum region.
Note that this model implies a nontrivial target structure that goes beyond the usual isotropic averages employed in CGC calculations, see below. While still lacking justification from first principles (although indeed CGC numerical calculations indicate that field correlations in the hadron wave functions are characterised by length scales \(\sim 1/Q_s\) [49,50,51]), this explanation is often used for qualitative discussion and understanding of numerical results, and may have further implications on e.g. spin or Transverse Momentum Distributions (TMDs) physics. Numerical studies based on models containing this domain structure have been performed in [52,53,54]. They show correlations that go beyond leading number of colours, see the discussion below, and lead to odd harmonics, see Sect. 3.
On the other hand, the glasma graph approach to two particle correlations is very successful to describe many features of the data as shown in [55,56,57,58,59,60], but the physics behind this approach was not clear. This issue has been studied in [61] and it has been shown that a genuine quantum effect, Bose enhancement of the gluons in the projectile wave function, leads to final state correlations in the glasma graph approach.^{Footnote 2}
The concept of Bose enhancement for a generic quantum system can be understood by considering a state with fixed occupation number, \(\{n_i(p)\}\), of N species of bosons at different momenta which, up to some normalisation factor, can be written as
with \(a^\dagger _i(p)\) the creation operator of the boson and \(i=1,2,\ldots , N\). The mean particle density \({\tilde{n}}\) is defined as the expectation value of the number operator in this state:
The two particle correlator in momentum space C(p, k) is defined in a similar way and can be calculated in a trivial manner:
The first term on the right hand side of Eq. (3) is the square of the mean particle density and the second term is the Bose enhancement term. It vanishes when the momenta of the two bosons are different and gives an enhancement when the momenta of two bosons coincide which is \({\mathcal O}(1/N)\), due to the fact that it contains a single sum over the species index. The physics behind this is the fact that only bosons of the same species are correlated with each other.
Let us now describe how Bose enhancement arises in the CGC and leads to final state correlations by considering the double inclusive gluon production within the glasma graph approach. In this approach each gluon is assumed to come from a different colour charge density in the projectile wave function that is rapidity invariant.^{Footnote 3} For our purposes, these colour charge densities can be conveniently represented in terms of gluon creation and annihilation operators in the incoming projectile wave function. After averaging over the target fields the glasma graphs can be written as sum of three types of diagrams (see Fig. 2).
Type A diagrams describe the case when two gluons with transverse momenta \(\mathbf{k}_1\) and \(\mathbf{k}_2\) scatter independently on the target, acquiring transfer of transverse momentum \(\mathbf{pk}_1\) and \(\mathbf{qk}_2\) so that the outgoing gluons have transverse momenta \(\mathbf{p}\) and \(\mathbf{q}\). Type B and Type C diagrams include interference contributions which are also interesting to study but, for now, let us focus on the Type A contribution and show how the Bose enhancement effect can be observed by studying these diagrams alone. The Type A contribution to the double inclusive gluon production can be written as
where \(\mathrm{in}\rangle \) is the wave function of the incoming projectile and \(N(\mathbf{pk})\) is the dipole scattering amplitude – the scattering amplitude for a two gluon system to scatter on the target. Moreover, the rapidity dependence of the gluon creation and annihilation operators is integrated over. The explicit dependence on rapidity becomes important only when the rapidity difference between the observed particles is parametrically large, \(\varDelta \eta \gtrsim 1/\alpha _s\).
The evaluation of the expectation value of any operator in the incoming projectile state requires a two averaging procedure in the CGC. In [61], averaging over the valence colour charge density is performed first which leads to the density matrix operator \(\hat{\rho }\) on the soft gluon Hilbert space. Then, the second averaging over the soft gluons is performed using this density matrix operator. The two particle correlator that appears in the Type A contribution calculated with this procedure leads to the following result:
where \(S_\perp \) is the transverse area of the projectile. The first term on the right hand side of Eq. (5) is the classical term which corresponds to the square of the number of gluons, while the second term is the typical Bose enhancement term, relatively suppressed by the number of states in the adjoint colour representation.
If we consider a situation where the incoming projectile has intrinsic saturation momentum \(Q_s\) and the momenta of the produced gluons are also \(\sim Q_s\), i.e. \(\mathbf{p}\sim \mathbf{q}\sim Q_s\), then the production amplitude is dominated by the contributions \(\mathbf{k}_1\sim \mathbf{k}_2\sim Q_s\). The initial state correlations are encoded in the Bose enhancement terms in Eq. (5), which are delta functions. The interaction with the target is obtained by convoluting the two particle correlator with the dipole amplitudes \(N(\mathbf{pk_1})N(\mathbf{qk_2})\). Since, in this kinematics, the momentum transfers from the target ( \(\mathbf{pk_1}\sim \mathbf{qk_2}\ll Q_s\) ) are small and the Bose enhancement terms involve delta functions, these initial state correlations naturally transform into angular correlations between the produced gluons in the final state. In a more general case, the delta functions, which are an artefact of considering a translationally invariant projectile, are smeared when convoluted with the dipole scattering amplitudes but this should not completely destroy the final state angular correlations.
The immediate question that arises after the study of gluons is whether quarks are subject to correlations in the CGC. This question has been posed in [64] where the correlations between the produced quarks were studied. The results in [61] show that the origin of the correlations between the produced gluons is the Bose enhancement of the projectile gluons. Due to their fermionic nature, one expects quarks to experience Pauli blocking which effectively amounts to a suppression of the probability of finding two quarks with the same quantum numbers in the CGC state. Therefore, one should expect a negative correlation between the final state quarks that originate from the initial state ones. On the other hand, the correlation between the gluons is found to be long range in rapidity since the CGC wave function is dominated by the rapidity integrated soft gluon field. Thus, another important question to answer is: are the (anti)correlations between the final state quarks long or short range in rapidity? The answer to this question is not obvious a priori. In the projectile wave function, quarks are produced via splitting of the rapidity invariant gluons into quarkantiquark pairs. However, the splitting amplitude itself depends on the rapidity of the quark and antiquark. Moreover, due to this splitting in the projectile wave function the expression for the production cross section of quarks is much more complicated compared to the one for gluons. These questions are answered in [64] where it was shown that the initial state correlations between the quarks in the projectile wave function are not distorted by the small momentum transfer from the target in specific kinematics. In these kinematics, the rapidity difference between the produced quarks is relatively large, i.e. \(\eta _1\eta _2\gg 1\). Moreover, a large contribution comes from the situation where the transverse momenta of the produced quarks \(\mathbf{p}\) and \(\mathbf{q}\) are of the same order and much larger than the saturation scale of the projectile \(Q_s\), and the saturation scale of the projectile is much larger than saturation scale of the target \(Q_T\), i.e. \(\mathbf{p}\sim \mathbf{q}\gg Q_s\gg Q_T\). Then the contribution to the production cross section that is sensitive to correlations has the following behaviour (see Eqs. (3.19) and (3.20) in [64] for the full expressions):
The negative sign of this contribution shows that it suppresses the classical term as opposed to the gluon case. This is the result of the Pauli blocking effect in quarkquark production. Moreover, this effect decays exponentially with the rapidity difference between the two produced quarks, which shows that it is short range in rapidity. However, this exponential decrease is tempered by two powers of the rapidity difference.
Besides, there is another physical effect present in the glasma graph approach which is referred to as the Hanbury–Brown–Twiss (HBT) correlations between the produced gluons^{Footnote 4} [67]. The diagrams in the glasma graph approach that lead to HBT correlations are those in Fig. 3 after performing pair wise contraction of the colour charges in the projectile wave function. Assuming a translationally invariant projectile wave function, the contribution from Type B and Type C diagrams to the production cross section is
If the translational invariance condition is relaxed, then the delta functions are smeared over a scale of the inverse size \(R^{1}\) of the projectile: \(\mathbf{p\pm q}\sim R^{1}\). This size R represents the radius of the gluon cloud inside the proton and its inverse is smaller than the saturation scale, \(R^{1}< Q_s\). Moreover, it is also shown in [67] that the HBT correlations are long range in rapidity just as the Bose enhancement effect. Thus, the strength of the HBT correlations is equal when the rapidities of the two produced gluons are similar (\(\eta _1\simeq \eta _2\)) or when the difference between them is large (\(\eta _1\eta _2\gg 1\)).
To sum up, the correlation function \(C(\mathbf{p,q})\), formally defined as the ratio of double inclusive gluon production cross section to the square of the single one, in the glasma graph approach contains two physical effects and can be written as follows:
The first term on the right hand side of Eq. (8) is the classical contribution which originates from the square of single inclusive production. \(C(\mathbf{p,q})\big _{\mathrm{BE}}\) represents the effect of Bose enhancement of the gluons in the projectile wave function. As described above, this effect leads to a correlation of the final state gluons. On the other hand, \(C(\mathbf{p,q})\big _\mathrm{HBT}\) represents the HBT correlations in the glasma graph approach which directly introduces correlations between the final state gluons. Both \(C(\mathbf{p,q})\big _{\mathrm{BE}}\) and \(C(\mathbf{p,q})\big _{\mathrm{HBT}}\) are rapidity independent, therefore long range in rapidity. The Bose enhancement contribution is suppressed by the transverse area of the projectile with respect to the HBT contribution (actually by the number of “sources” \(Q_s^2 S_\perp \)). However, it leads to correlations whose width in momentum space is determined by the saturation momentum \(Q_s\). On the other hand, the HBT contribution is not suppressed but it gives a narrow peak in momentum space with width \(R^{1}\). This comparison is shown in Fig. 4.
In the explicit calculations in the glasma graph approach to double inclusive particle production,^{Footnote 5} the averaging over the target configurations that leads to the dipole scattering amplitude is performed expanding this amplitude to the lowest non trivial order in the target field strength, corresponding to two gluon exchange between the gluons in the projectile wave function and the target. The dipole scattering amplitudes \(N(\mathbf{pk_1})\) and \(N(\mathbf{qk_2})\) introduced in Eq. (4) are assumed to originate from single pairs of target fields and therefore this approach does not take into account the effects of multiple scatterings in a dense target. Therefore, this approach is only valid for pp collisions^{Footnote 6}. In [70], the inclusive production of two and three gluons is computed beyond the glasma graph approach by including the multiple scattering effects, which extends the validity of the glasma graph approach from pp to pA collisions.^{Footnote 7}
Apart from taking into account the multiple scattering effects in [70], a systematic way to identify each term in the double inclusive gluon production cross section and to characterise whether it is a Bose enhancement or HBT contribution is introduced. This identification is performed by adopting the following strategy. When calculating the double inclusive gluon production, one has to average over four colour charges (two in the amplitude and two in the complex conjugate amplitude) in the projectile wave function: \(\langle \rho ^{a_1}(\mathbf{x_1})\rho ^{a_2}(\mathbf{x_2})\rho ^{b_1}(\mathbf{y_1})\rho ^{b_2}(\mathbf{y_2})\rangle _P\). Here, \((\mathbf{x_i},a_i)\) and \((\mathbf{y_i},b_i)\) stand for the transverse position and colour indices of the colour charge densities in the amplitude and in the complex conjugate amplitude, respectively. The averaging over the colour charge distributions in the projectile is commonly performed by using a generalised McLerran–Venugopalan (MV) model [73, 74] where the weight functional is Gaussian. Then, the average of any product of colour charge densities factorises into a product of all possible pair, Wicklike contractions. The correlator of two colour charge densities in momentum space can be defined as
The function \(\mu ^2(\mathbf{k, p})\) characterises the structure of the projectile. It can be written as
where \(F\big [ (\mathbf{k+p})R\big ]\) is a soft form factor with maximal value F(0), and R is the radius of the projectile. Function T defines the transverse momentum dependent distribution of the valence charges. The soft form factor identifies whether a term is a contribution to the Bose enhancement of the projectile gluons or a contribution to the HBT correlations of the produced gluons. For example, in our set up the produced gluons have momenta \(\mathbf{p}\) and \(\mathbf{q}\), while the projectile gluons carry transverse momenta \(\mathbf{k_1}\) and \(\mathbf{k_2}\). In this case, \(\mu ^2(\mathbf{p,q})\) gives a maximal contribution when \(\mathbf{p+q}=0\) which clearly can be identified as the HBT correlations of the produced gluons. \(\mu ^2(\mathbf{k_1,k_2})\) is peaked when \(\mathbf{k_1+k_2}=0\) which is a contribution to the Bose enhancement of the gluons in the projectile wave function.
On the other hand, multiple scattering effects on the dense target are taken into account by introducing the standard Wilson lines in the CGC framework. In this framework, the interaction between the projectile and the target is assumed to be eikonal which amounts to the situation where each parton produced by the projectile colour charge scatters on the target by picking up a colour rotation described by a Wilson line which is defined as an exponential of the target field ordered in the \(x^+\) coordinate:
at the amplitude level. Here, \(T^a_{{\mathcal {R}}}\) is the \(SU(N_c)\) generator in the representation \({{\mathcal {R}}}\) which can be the fundamental one for a quark and the adjoint one for a gluon. This leads to the appearance in the cross section of double dipole and quadrupole amplitudes (in the adjoint representation) of the type
which have to be averaged over the target field distributions. The dipole and the quadrupole operators are defined as
The cross section has to be integrated over four transverse coordinates. In principle, the maximal contribution should come from the area in coordinate space, i.e. when all the four coordinates are far away from each other. However, all four points cannot be far away from each other since the target field ensemble has to be colour neutral, and colour neutralisation in the CGC happens on scales of order \(1/Q_s\). Therefore, the maximal contribution to the integral must come from the configurations where the four points are combined into pairs, such that each pair is a singlet and the distance between the pairs is large. This is the leading contribution in \(1/(Q_s^2 R^2)\) to the integral on transverse coordinates^{Footnote 8} – not, by any means, a good representation of the target averages of ensembles of Wilson lines by themselves. Taking into account only such configurations is equivalent to calculating the target averages of products of any number of Wilson lines by factorising them into averages of pairs with basic Wick contraction. In this case, target averaging of the double dipole and the quadrupole amplitudes can be written as
where we have defined \(d(\mathbf{x,y})\equiv \langle s(\mathbf{x,y})\rangle _T\). Then, by using the function \(\mu ^2(\mathbf{k,p})\) given in Eq. (10) for the projectile colour charge density correlators and using the factorisation ansatz described above for the double dipole and quadrupole amplitudes, the double inclusive gluon production cross section is computed and the nature of all the terms is identified. Moreover, it is also shown that the contributions to final state correlations comes from the quadrupole terms in two gluon production.^{Footnote 9}
3 Odd azimuthal harmonics from the CGC
To describe one key existing problem in usual CGC calculations, namely the absence of odd harmonics, let us briefly discuss double inclusive particle production in more depth. Within the approximations described in Sect. 2, double inclusive gluon production is computed in pA collisions in [70]. The production cross section of two gluons with rapidities \(\eta _1\) and \(\eta _2\), and transverse momenta \(\mathbf{k}_1\) and \(\mathbf{k}_2\) (see Fig. 5) reads
where
Here, function \(\mu ^2\) defines the structure of the projectile, Eq. (10). On the other hand, function \(L^i(\mathbf{k},\mathbf{q})\) is the eikonal Lipatov vertex which is defined as
Note that, as discussed in Sect. 2, all correlations are subleading in \(1/N_c\) which is due to the use of Gaussian averages – the MV model. Moreover, the double gluon inclusive production cross section is written in terms of the dipole averages assuming translational invariance of the target (a standard approximation which is reasonable in pA for large nuclei):
A convenient way to study two particle correlations is through a Fourier decomposition into harmonics defined in for the azimuthal angle \(\varDelta \phi \) between the produced gluons with transverse momenta \(\mathbf{k}_1\) and \(\mathbf{k}_2\). When Fourier expanded, the double inclusive gluon spectrum Eq. (17) can be written as
where
One way of defining the \(p_T\) dependence of the Fourier coefficients is by fixing one of the momenta (\(k_1=p^{ref}_T\)), and treating the other one as a free variable (\(k_2=p_T\)). With this choice, the azimuthal harmonics are defined as
The key theoretical problem for the description of the two particle correlations within the CGC is the absence of the odd harmonics. This problem is analyzed in [41], and it is shown to be due to the symmetry \((\mathbf{k}_2\rightarrow \mathbf{k}_2)\) of the double inclusive gluon production cross section Eq. (17) which is also referred to as the “accidental symmetry of the CGC”.
Three ways^{Footnote 10} have been proposed to solve this problem.^{Footnote 11} On the one hand, the projectile and target can be characterised by a more involved structure than that considered in the usual MV averages [40,41,42, 52,53,54].
On the other hand and as discussed previously, within the CGC framework each produced gluon originates from a separate colour charge density in the projectile wave function as shown in Fig. 5. The contributions to the projectile wave function that emerge from merging of the gluons before the interaction with the target or splitting of a gluon into two gluons emitted from the same colour charge density, are not taken into account in the standard CGC calculations. Recently, in [80, 81] it is shown that the accidental symmetry of the CGC can be broken by including such corrections to the projectile wave function. The even and odd parts of the double inclusive gluon production cross section under the accidental symmetry are computed separately, and finally, the azimuthal harmonics are calculated. The corresponding numerical studies and a comparison with data are performed in [82, 83].
Finally, inclusive gluon production is usually studied within the eikonal approximation in the CGC framework. In recent studies [84, 85], it is shown that the accidental symmetry can be broken by going beyond this eikonal approximation. In the next subsection, we introduce a systematic way to include subeikonal corrections in CGC calculations and show how these corrections give rise to nonvanishing odd harmonics.
3.1 Subeikonal corrections in the CGC
In inclusive gluon production at central rapidity in pA collisions, both the projectile and the target are highly energetic since they are boosted from their initial rapidity to the central rapidity where the collision occurs. Therefore, in this case both colliding objects can be treated in the CGC framework. This corresponds to defining the projectile by the colour charge \(J^\mu _a(x)\),
and the target by the colour field \(A^\mu _a(x)\) that is given as
Let us recall that these expressions of the colour charge of the projectile and the colour field of the target are defined within the eikonal approximation, which is justified by the large energy of both colliding objects. If for the dilute projectile the eikonal approximation can be trusted at a given energy, the same approximation for a large target can be true only for larger energies. The eikonal approximation for the target amounts to the following three conditions:

1.
\(A^{\mu }_a(x)\simeq \delta ^{\mu }A^_a(x)\): Neglecting the (+) and transverse components of the colour field of the target.

2.
\(A^{\mu }_a(x)\simeq A^{\mu }_a(x^+,\mathbf{x})\): Neglecting the \(x^\) dependence in the colour field of the target.

3.
\(A^\mu (x)\propto \delta (x^+)\): Assuming that the target field is peaked around \(x^+=0\) due to Lorentz contraction, which is also known as the shockwave approximation.
In realistic kinematical conditions under which the experiments are performed, the energies are not asymptotic and the eikonal approximation is not always justified. While for a dilute projectile it is usually valid even for high energy collisions, this is not necessarily true for a large nucleus. Relaxing any of the above approximations accounts for corrections to the eikonal limit. In [86, 87], a systematic method to compute the corrections to the eikonal limit by relaxing the third approximation is developed. This corresponds to treating the colour field of the target with a finite longitudinal support \(L^+\) along the \(x^+\) direction, thus replacing Eq. (27) by
Such subeikonal corrections are thus subleading with respect to the infinite Lorentz contraction of the target.
Before discussing the results, let us give a brief sketch of the method employed to derive the noneikonal corrections. Let us consider the production of a single gluon with transverse momenta \(\mathbf{k}\) and longitudinal momenta \(k^+\) in pA collisions at central rapidity. The dilute projectile is still treated in the eikonal approximation and defined with the charge density \(J^\mu _a(x)\) given in Eq. (26). On the other hand, the eikonal approximation is relaxed for the dense target that is defined by the colour field \(A^\mu _a(x)\) given in Eq. (28) with a finite support from 0 to \(L^+\) in the longitudinal direction. In this case, the production cross section can be written as the square of the gluon production amplitude averaged over the projectile and target distributions and integrated over impact parameter \(\mathbf{B}\):
Here, \(\lambda \), a and \({\underline{k}}=(k^+, \mathbf{k})\) are the polarization, colour and momentum of the produced gluon^{Footnote 12}. For a target with finite longitudinal width, the gluon production amplitude \(\mathcal{M}^a_\lambda (\underline{k}, \mathbf{B})\) is composed of three different contributions: gluon production before, while and after the projectile propagates through the target. At leading order, it is possible to relate the total gluon production amplitude and the background retarded gluon propagator by using the LSZ reduction formula and the perturbative expansion of the colour field of the target [88]. In the light cone gauge \(A^+=0\), the total gluon production amplitude can be written in terms of the \((i)\) component of the background retarded gluon propagator \(G_R^{\mu \nu }(x,y)\) as
Since the colour field of the target is independent of \(x^\), one can introduce the onedimensional Fourier transform of the background retarded gluon propagator and write it in terms of of the background scalar propagator \({{\mathcal {G}}}_{k^+}^{\mu \nu }({\underline{x}}, {\underline{y}})\). Then, the \((i)\) component of the retarded background gluon propagator reads
The background scalar propagator \({{\mathcal {G}}}_{k^+}^{ab}(\underline{x}, \underline{y})\) satisfies the scalar Green’s equation whose solution formally can be written as a path integral
with the Wilson line
following the Brownian trajectory \(\mathbf{z}(z^+)\). In the limit of vanishing longitudinal width, \(x^+y^+\rightarrow 0\), the background scalar propagator \({{\mathcal {G}}}_{k^+}^{ab}(\underline{x}, \underline{y})\) reduces to the standard Wilson line introduced in Eq. (11) and one recovers the eikonal limit. Therefore, it can be safely concluded that all the noneikonal effects that are due to the finite longitudinal width of the target are encoded in the background scalar propagator. This also means that an expansion of \({{\mathcal {G}}}_{k^+}^{ab}(\underline{x}, \underline{y})\) can be performed in terms of an eikonal parameter, with the first term in this expansion corresponding to the eikonal limit and higher order terms to the corrections to this limit.
In order to perform an eikonal expansion of the background scalar propagator \({{\mathcal {G}}}_{k^+}^{ab}(\underline{x}, \underline{y})\), one should first discretise the scalar background propagator. In the eikonal limit, \(k^+/(x^+y^+)\) is much larger than any transverse scale in the problem. In the large \(k^+\) limit, it is natural to consider a generic path as a perturbation around the classical free path,
where the transverse positions at step n are on the straight line
between the initial and final points, and the perturbation \(\mathbf{u}_n\) satisfies the boundary conditions \(\mathbf{u}_0=\mathbf{u}_N=0\) with N being the number of discretised steps. Once the expansion around the free classical path is performed for fixed initial and final positions, one should perform another expansion for small \(\mathbf{x}\mathbf{y}\), since \(\mathbf{x}\mathbf{y}\) is parametrically small in the large \(k^+\) limit. After performing these two expansions up to second order in \((x^+y^+)\) – the finite longitudinal width of the target, the scalar background propagator \(\mathcal{G}_{k^+}^{ab}(\underline{x}, \underline{y})\) reads
The first term on the right hand side of Eq. (36) is the standard Wilson line defined in Eq. (11). \({{\mathcal {O}}}\big [(x^+y^+)/k^+\big ]\) terms are the first order corrections to the strict eikonal limit which we refer to as nexttoeikonal (NEik) corrections. Similarly, \({{\mathcal {O}}}\big [ (x^+y^+)^2/(k^+)^2\big ]\) terms are the second order corrections and they are referred to as nexttonexttoeikonal (NNEik) corrections. The terms that are denoted as \(U_{[\alpha ,\beta ]}(x^+,y^+; \mathbf{y})\) are the decorated Wilson lines which only appear beyond strict eikonal order. The first subscript \(\alpha \) in the decorated Wilson lines stands for the order of expansion around the classical path while the second subscript \(\beta \) stands for the order of the expansion around the initial transverse position \(\mathbf{y}\). The reason why these objects are referred to as decorated Wilson lines is related with their structure. These objects involve a background field insertion into the standard Wilson lines along the \(+\)direction in a given \(+\)coordinate. For example, the first decorated Wilson line is defined as
The other decorated Wilson lines have similar structure with one or more background field insertions. We do not present the structure of all the decorated Wilson lines due their complexity and lengthy expressions (see [87]). One can easily get the expression for the gluon production amplitude at NNEik accuracy given in Eq. (30) by using the expression of the retarded background gluon propagator Eq. (31) and the expression derived for the background scalar propagator Eq. (36).
As discussed above, the retarded background gluon propagator \(G_R^{\mu \nu }(x,y)_{ab}\) and, therefore, the scalar background propagator \({{\mathcal {G}}}_{k^+}^{ab}(\underline{x}, \underline{y})\) are the main building blocks for computing cross sections in high energy pA collisions. In [86, 87], these propagators were used to calculate the single inclusive gluon production cross section in pA collisions at NNEik accuracy. The same formalism can be adopted to compute double inclusive gluon production and hence the azimuthal harmonics in pA collisions [84, 85].
In [89], the results of the single inclusive gluon production cross section at NNEik accuracy in pA collisions are used to study the weak field limit (i.e. glasma graph approximation) of this result which corresponds to single inclusive production in pp collisions. In this limit, the decorated Wilson lines are expanded to first order in the background field of the target \(A_a^(z^+,\mathbf{y})\). For example, the first decorated Wilson line given in Eq. (37) reduces to
This simplification allows us to calculate the Lipatov vertex at NNEik accuracy. After expanding the eikonal and noneikonal terms to first order in powers of the background field, which corresponds to the glasma graph approach in usual CGC calculations, the Lipatov vertex at NNEik accuracy can be written as
The first term on the right hand side of Eq. (39) corresponds to the strict eikonal limit, thus it gives the eikonal Lipatov vertex defined in Eq. (21). The second and the third terms are the NEik and NNEik corrections respectively. The structure of the vertex suggests that the corrections to the amplitude due to finite width of the target may exponentiate.
This observation is further studied recently in [84] where it was shown that indeed the noneikonal corrections due to finite longitudinal width of the target in the weak field limit exponentiate and can be written as modified Lipatov vertex. By computing the corresponding three diagrams (see Fig. 6) at the amplitude level and keeping the phase \(e^{ik^x^+}\) which is taken to be one in the eikonal limit, the noneikonal Lipatov vertex that accounts for all order corrections to the eikonal limit due to finite longitudinal width of the target in the weak field limit reads
where \(k^\equiv \mathbf{k}^2/(2k^+)\). This structure was observed in the context of jet quenching in [90,91,92] previously, however the identification of this building block for its use to include noneikonal corrections in CGC calculations is done in [84] for the first time, further illustrating the close relation between CGC and jet quenching calculations [88].
It is now straightforward to compute the noneikonal single inclusive gluon production in pp (i.e. dilutedilute) collisions which formally reads
An additional modification that is needed to account for the finite longitudinal width of the target is to adopt a modified expression for the correlator of two target fields. Motivated by the non zero longitudinal extent of the target, the fields can be located at different positions that are separated by the colour correlation length in the target \(\lambda ^+\). In this case, the two target field correlator reads
where function \(n(x^+)\) defines the onedimensional target density that we take as constant, \(n_0=n(x^+)\) for \(0\le x^+\le L^+\), and 0 elsewhere. Moreover, function \(a(\mathbf{q})\) is the potential in momentum space which can be taken to be of Yukawa type, i.e. \(a(\mathbf{q})^2=m^2/(\mathbf{q}^2+m^2)^2\) with m being the Debye screening mass or inverse colour correlation length. In the eikonal limit, when \(\lambda ^+\rightarrow 0\) for a constant potential and constant onedimensional target density, one recovers the standard MV expression for the two target field correlator. Using Eq. (42) one can integrate over the longitudinal coordinates that appear in the phases in noneikonal Lipatov vertices. The final result of the noneikonal single inclusive gluon production cross section in pp collisions then reads
where we assume that the longitudinal width of the target is much larger than the colour correlation length, \(\lambda ^+\ll L^+\). In the cross section, Eq. (43), all the noneikonal effects are encoded in the function \({{\mathcal {G}}}_1^\mathrm{NE}(k^;\lambda ^+)\) which is defined as
with \(k^\equiv \mathbf{k}^2/2k^+\). In the limit of vanishing \((k^\lambda ^+)\) we have
and the well known eikonal limit for the single inclusive gluon production in the dilute target limit is recovered. Therefore, function \({{\mathcal {G}}}_1^\mathrm{NE}(k^;\lambda ^+)\) can be interpreted as the function that accounts for the relative importance of the noneikonal effects with respect to the eikonal limit of the single inclusive gluon production in the dilute target limit.
In Fig. 7, the ratio of the noneikonal to eikonal single inclusive gluon production cross sections, i.e. function \({{\mathcal {G}}}_1^\mathrm{NE}(k^;\lambda ^+)\), is plotted as a function of rapidity for different values of the transverse momenta of the produced gluon at correlation length \(\lambda ^+=0.5\) fm. The results show that with increasing rapidity of the produced gluon, the effects of the noneikonal corrections vanish as expected from analytical predictions. Up to rapidity \(\eta =2.5\), the relative importance of the corrections varies between 15 and 2% depending on the value of the transverse momenta.
Noneikonal double inclusive gluon production cross section in pp scattering can be computed in a similar manner. The main difference between the single and double inclusive production is that one needs to compute the target average of the four field correlator. This can be accomplished by factorising the average of the colour fields of the target into all possible Wick contractions which can be written
where each two target field correlator is defined in Eq. (42). The integrals over the longitudinal coordinates can be performed using this definition and the final result for the noneikonal double inclusive gluon production cross section in pp scattering can be organised in the following way:
where the subscripts denote the single trace terms (\(I^{(i)}_\mathrm{1tr}\)) or the double trace term (\(I^{(i)}_{\mathrm{2tr}}\)) which are analogue to double dipole and quadrupole operators in pA scattering discussed previously. The explicit expressions for these contributions read
and, finally,
In this setup (see Fig. 5), \(\mathbf{k}_1\mathbf{q}_1\) and \(\mathbf{k}_2\mathbf{q}_2\) are the transverse momenta of the two gluons in the projectile, \(\mathbf{k}_1\) and \(\mathbf{k}_2\) are the transverse momenta of the produced gluons in the final state, and \(\mathbf{q}_1\) and \(\mathbf{q}_2\) are the transverse momenta transferred from the target to the projectile during the interaction. By using the definition of function \(\mu ^2(\mathbf{p},\mathbf{k})\) given in Eq. (10) and the behaviour of the soft form factor, one can easily identify each term in the noneikonal double inclusive gluon production cross section given in Eqs. (48), (49) and (50). Clearly, the term in \(I^{(0)}_{\mathrm{2 tr}}\) corresponds to the square of the single inclusive production and does not give any contribution to the correlations. The terms in \(I^{(1)}_{\mathrm{2 tr}}\) corresponds to the Bose enhancement of the target gluons since the soft form factor is peaked around \(\mathbf{q}_1=\mathbf{q}_2\). Finally, the terms in \(I^{(1)}_{\mathrm{1 tr}}\) contribute to the HBT correlations of the produced gluons and to Bose enhancement of the projectile gluons.
In the noneikonal double inclusive gluon production cross section, two functions appear that account for the noneikonal effects: \({{\mathcal {G}}}_1^\mathrm{NE}(k_1^;\lambda ^+)\) presented in Eq. (44) and a new function \({{\mathcal {G}}}_2^\mathrm{NE}(k_1^,k_2^;L^+)\). This new function is defined as
which in the eikonal limit, i.e. \(L^+\rightarrow 0\), goes to unity,
Different from the eikonal double inclusive gluon production cross section, in the noneikonal expression the mirror images are given by \({\underline{k}}_2\rightarrow {\underline{k}}_2\) where \({\underline{k}}_2\equiv (k_2^+,\mathbf{k}_2)\). The mirror images of the terms that are accompanied by the function \({{\mathcal {G}}}_2^\mathrm{NE}(k_1^,k_2^;L^+)\) are now accompanied by \({{\mathcal {G}}}_2^\mathrm{NE}(k_1^,k_2^;L^+)\). However, as obvious from the definition given in Eq. (51), this function is not symmetric under this transformation. Moreover, in certain kinematic regimes the behaviour of \({{\mathcal {G}}}_2^\mathrm{NE}(k_1^,k_2^;L^+)\) differs completely from \({{\mathcal {G}}}_2^\mathrm{NE}(k_1^,k_2^;L^+)\). Particularly, in a kinematic region where \(k_1^\sim k_2^\), one gets
This creates an asymmetry between the terms (\({\underline{k}}_1,{\underline{k}}_2\)) and their partners (\({\underline{k}}_2\rightarrow {\underline{k}}_2\)). This asymmetry which comes from the noneikonal effects reminds the asymmetry between the forward and backward peaks of the ridge structure observed in two particle production. Therefore, noneikonal corrections break the accidental symmetry present in usual CGC calculations and can give rise to odd harmonics. In the remaining of the section we briefly examine the numerical relevance of this effect.
A detailed numerical analysis of the azimuthal structures in two particle correlations based on the noneikonal double inclusive gluon production cross section given in Eq. (47) with Eqs. (48), (49) and (50), is performed in [85] where it is assumed that:

1.
the colour sources inside the projectile have a Gaussian distribution such that \(\mu ^2(\mathbf{k}, \mathbf{q})=\mu ^2(2\pi )^2\delta ^{(2)}(\mathbf{k}+\mathbf{q})\), with \(\mu \) being the width of the Gaussian;

2.
the Yukawa type potential that defines the target field correlators is given by \(a(\mathbf{q})^2=\mu _T^2/(\mathbf{q}^2+\mu ^2_T)^2\), with \(\mu _T\) being an infrared regulator analogous to a Debye mass;

3.
the transverse area of the projectile \(S_\perp \) is defined through \((2\pi )^2\delta ^{(2)}(\mathbf{q}\mathbf{q})\rightarrow S_{\perp }\).
In this analysis, function \({{\mathcal {G}}}_2^\mathrm{NE}(k_1^,k^_2;L^+)\) that encodes the noneikonal effects defined in Eq. (51) is rewritten as
using \(k^=k^2/2k^+\), \(k^+=ke^\eta /{\sqrt{2}}\), \(k=\mathbf{k}\).^{Footnote 13}
In Fig. 8 the azimuthal harmonics up to \(v_5\) are computed by using the definition given in Eq. (25). \(p_T^{ref}\) is taken to be \(1\, \mathrm{GeV}\) for different values of \(\sqrt{s_{NN}}\) and \(\eta _1=\eta _2=\eta \). The plot shows that the value of the odd harmonics decreases with increasing centreofmass energy at fixed rapidity \(\eta \). This behaviour is the natural outcome of the fact that noneikonal corrections become smaller with the increasing Lorentz gamma factor. Therefore, one can conclude that the noneikonal corrections can be negligible for collisions at high centreofmass energy such as the ones at the LHC but they can be important for collisions at RHIC with \(\sqrt{s_{NN}}\le 200\, \mathrm{GeV}\). On the other hand, at any fixed energy the value of odd azimuthal harmonics decreases with increasing rapidity \(\eta \). This behaviour is also expected, since the value of the odd harmonics is directly related to the noneikonal corrections. The eikonal expansion parameter can be written as \(p_TL^+e^{\eta }\) in terms of rapidity and, therefore, noneikonal corrections (and thus the value of odd harmonics) decrease with increasing rapidity and vanish completely in the strict eikonal limit.^{Footnote 14}
In Fig. 9, the ratio \(v_n(L^+)/v_n(1.5\, \mathrm{fm})\) is plotted as a function \(L^+\) which reveals a very interesting feature of the effects of noneikonal corrections on azimuthal harmonics: odd harmonics depend strongly on the size of the target while even ones are almost independent of it. Even though the explicit relation between multiplicity and \(L^+\) requires a more dedicated study, the scaling of the odd harmonics with \(L^+\) in Fig. 9 qualitatively resembles the results of the analysis performed in [82] where it is shown that the value of \(v_3\) increases with the increasing multiplicity.
4 NonCGC explanations
Besides explanations to the ridge phenomenon based on the CGC, there are others that address its origin in the initial state of the collision or, at least, do not demand hydrodynamics or transport at work. It must be noted that the existence of long range rapidity correlations was discussed long ago as a consequence of multiple scattering, see [93, 94].
This approach was pushed forward in string models for multiparticle production, see e.g. [95]. Later on, several models that consider string interactions were argued to lead to azimuthal asymmetries: string percolation [96,97,98] with the creation of azimuthally anisotropic strong chromoelectric fields, colour reconnection [99] that is able to produce some of the QGPlike features observed in pp, and string repulsion [100, 101]. It is not clear whether the dynamics contained in these approaches can be considered as pure initial state but they offer a mechanism to produce the ridge in collisions between small systems that does not require any explicit final state rescattering, see [102] for a model that explicitly requires parton and/or hadron rescattering to build azimuthal asymmetries even with just two strings. In all these approaches, particle production from a single string is still isotropic and the anisotropy is built after string breaking.
A stringbased model is also proposed in [103, 104]. There, valence diquarkquark flux tubes or strings in the incoming protons overlap and produce more particles in the transverse than in the longitudinal direction of the flux tube. Such anisotropic particle production leads to azimuthal asymmetries and the prediction has been made that it should also be visible in photoproduction, with particle production becoming maximal in the plane of the deflected electron (in ep collisions) or proton (in UPCs).
It should be noted that all these approaches are inspired in the string behaviour of the QCD interaction in the nonperturbative domain, in contrast to the CGC that relies on perturbation theory for a small coupling constant. Indeed already in the framework of Reggeon field theory, some ideas have been pushed [105] on the spatial variation of the transverse density in the hadron that resemble those in the CGC. Or CGC arguments have been extended to the soft physics domain and applied to describe azimuthal correlations, see [106, 107] and subsequent papers of this group.^{Footnote 15}
Finally, let us indicate that azimuthal asymmetries arise in several processes when the nucleon is studied and characterised beyond collinear parton densities, as in the framework of Wigner distributions and TMD parton densities, see [110, 111]. Azimuthal asymmetries then arise in final observables like dijet production in DIS [112,113,114]. But, although these calculations are often performed in a framework close to that of the CGC which is related with the TMD framework at small x [115], it goes beyond the standard CGC context to link with other physics like spin.
5 Summary and discussions
In this manuscript we have discussed the explanations that are currently proposed to describe particle correlations, the ridge, observed in experimental data in small collision systems, pp and pA, from the initial state point of view. Our main focus has been those weak coupling explanations based on the CGC. We have assumed that correlations among partons in the initial stage leave an imprint on those among particles in the final state, i.e. they are not washed out by final strong final state interactions and hadronisation.
First, we have focused on the standard eikonal treatment within the glasma graph approximation which is valid for collisions between dilute objects – pp. We have reviewed the studies which have shown that this approximation encodes two different type of contributions, namely the Bose enhancement of both projectile and target gluons and also HBT correlations of the produced final gluons.
We have summarised the procedure that should be adopted to extend the validity of glasma graph approximation from dilutedilute to dilutedense (i.e. from pp to pA) collisions by taking into account the multiple scattering effects in the dense target. We have shown that the structure of the double inclusive gluon production cross section is symmetric under (\(\mathbf{k}_2\rightarrow \mathbf{k}_2\)), which is known as the accidental symmetry of the CGC. Since this symmetry is the reason for vanishing odd harmonics in the CGC framework, we have discussed the suggested mechanisms to break this accidental symmetry.
In particular, we have focused on a specific mechanism to break this symmetry which is based on going beyond the standard eikonal approximation and including the subeikonal corrections that are due to the finite longitudinal width of the target. We have argued that such noneikonal corrections, when included in the glasma graph approach to two particle correlations, successfully generate nonzero odd harmonics in specific kinematics. We would like to emphasise here that we make no attempt to compare the results with experimental data but only address the existence and size of the noneikonal effects on the azimuthal structure. As expected from noneikonal corrections, their value and thus that of the odd harmonics decrease rapidly with increasing centreofmass energy. This decrease is strong since the analysis is performed for a dilute target – a slower decrease of the size of the odd harmonics with increasing energy could be expected in a dilutedense collision. Besides, the treatment of such nonleading eikonal corrections shows explicitly the link of the formalisms used in CGC and jet quenching calculations.
At this point, we should comment briefly on the comparison with experimental data. The main characteristics concerning azimuthal asymmetries in small systems observed in experiment [1,2,3,4,5, 7,8,9,10, 13, 16,17,19] are:

The even and odd harmonics extracted using correlations between two and more particles, are of similar size to those found in larger systems.

They show the same dependence on the mass of the measured hadron as found in larger systems (see [116] for an approach in the CGC).

Even harmonics show a much weaker dependence on the multiplicity in the event than odd harmonics.

\(v_2\) and \(v_3\) found in pAu, dAu and \(^3\)HeAu collisions at RHIC show, for central (headon) collisions, the ordering \(v_2^{p\mathrm {Au}}<v_2^{d\mathrm {Au}}\approx v_2^{^3\mathrm {HeAu}}\), \(v_3^{p\mathrm {Au}}\approx v_3^{d\mathrm {Au}}< v_3^{^3\mathrm {HeAu}}\).

Measurements of many particle cumulants show evidence of collectivity. For example, four particle cumulants \(c_2\{4\}=\langle e^{i2(\phi _1+\phi _2\phi _3\phi _4)}\rangle 2\langle e^{i2(\phi _1\phi _2)} \rangle ^2\) (\(v_2\{4\}=[c_2\{4\}]^{1/4}\)) change sign from positive to negative with increasing associated multiplicity, with a smooth behaviour from small to large systems and from smaller to larger energies.
In the glasma graph approximation, several studies were done that describe pp data in a reasonable manner [55,56,57,58]. Later on, these studies were extended to pA with diverse degree of modelling [59, 60]. Then, odd azimuthal harmonics were introduced following [80, 81], with the corresponding numerical studies and a comparison with data performed in [82, 83]. In these latter studies a successful comparison with RHIC and LHC data was initially claimed, which was later corrected after the criticism in [117]. Nevertheless, it must be stated that none of the numerical calculations can be considered as a full implementation of the theoretical framework and that some results are still to be clarified from an analytical point of view, e.g. those in [59, 60] about the second Fourier coefficient defined through four particle correlations (\(v_2\{4\}\)) where the mentioned change of sign in \(c_2\{4\}\) is attributed to multiple scattering beyond glasma graphs.
Finally, we have also shortly commented on approaches that are not based on, or go beyond, CGC ideas, to study the two particle correlations from the initial state.
To conclude, let us indicate that future experimental programmes and facilities [27, 118,119,120,121] will address the physics of small systems and the transition from small to large, particularly the onset and understanding of collectivity which is a central question in QCD at high energies.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: This review focuses on the most recent theoretical developments in the study of particle correlations from the initial state point of view. Therefore, there is no experimental data used in this review.]
Notes
Apart from these two approaches, there are also other CGCbased mechanisms to describe the two particle correlations from the initial state. In [45, 46], it is argued that long range rapidity correlations can be explained by the spatial variation of the partonic density in the target. On the other hand, in [47], the origin of two particle correlations is explained through the rapidity evolution of dipole operators by breaking the mean field approximation. Correlations in the hadron wave functions as described in the CGC have been recently considered in [48] but with the aim of providing initial conditions for hydrodynamic evolution beyond simple energy, flow and particle density transverse profiles.
In [62, 63] a collective behaviour and azimuthal asymmetries are obtained for small systems in a perturbative framework as a consequence of quantum interference and colour correlations. Spatial anisotropies result into momentum anisotropies via multipole radiation patterns. This approach, albeit formulated in a rather different language, shows similarities with the glasma graph approach.
These two assumptions are justified at high enough energy or at small x where the colour charge density is high so gluons can be treated semiclassically, and they have evolved radiation tails that populate phase space uniformly in \(\ln 1/x\).
The previously discussed conditions in the projectile of large colour charge and rapidity independence of the gluon distribution are assumed to hold in pp.
This extension is studied numerically in [71], and also analytically in [72]. The main difference between the latter and [70] is the computation framework. Two gluon correlations in [72] are computed within the \(k_T\)factorised approach which is difficult to generalise to three or more particles. Moreover, the results of [72] are valid only in the large \(N_c\) limit as opposed to the results in [70] which are valid for finite \(N_c\). In this sense, [72] can be considered as the first attempt to generalise the glasma graph approach to two gluon production from pp to pA.
Besides, the role of the centrality or multiplicity event selection for the breaking of the accidental symmetry and the appearance of odd azimuthal harmonics has been analysed in [79].
Hereafter, we use the underline notation to indicate that for coordinates \(\underline{x}=(x^+,\mathbf{x})\) and for momentum \(\underline{k}=(k^+,\mathbf{k})\).
Assuming that L is the size of the target in its rest frame, in the centreofmass frame \(L^+\) is taken as
where A is the mass number of the nucleus and \(\gamma \simeq \sqrt{s_{NN}}/(2m_N)\) accounts for the Lorentz contraction in the centreofmass frame. Moreover, for the numerics the gluonic size of the projectile is taken to be \(B_p=4\, \mathrm{GeV^{2}}\), the transverse size of the projectile is assumed to be \(S_{\perp }=2\pi B_p\approx 9.8\, \mathrm{mb}\) and the size of the target in its rest frame for a Pb nucleus is taken to be \(L=12\, \mathrm{fm}\). Finally, the number of colours is taken to be \(N_c=3\) and the colour correlation is set to be \(\lambda ^+=0\).
In Fig. 8, the unrealistic peaks that account for the HBT contributions are due to the use of \(\mu ^2(\mathbf{k},\mathbf{q})\propto \delta ^{(2)}(\mathbf{k}+\mathbf{q})\). In a more realistic treatment, \(\mu ^2\) can be chosen as Gaussian which would peak around \(\mathbf{k}+\mathbf{q}=0\) and show a bell shape behaviour.
References
V. Khachatryan et al. (CMS), JHEP 09, 091 (2010)
M. Aaboud et al. (ATLAS), Phys. Rev. C 96, 024908 (2017)
V. Khachatryan et al. (CMS), Phys. Rev. Lett. 116, 172302 (2016)
G. Aad et al. (ATLAS), Phys. Rev. Lett. 116, 172301 (2016)
V. Khachatryan et al. (CMS), Phys. Lett. B 765, 193 (2017)
M. Aaboud et al. (ATLAS), Eur. Phys. J. C 80, 64 (2020)
S. Chatrchyan et al. (CMS), Phys. Lett. B 718, 795 (2013)
B. Abelev et al. (ALICE), Phys. Lett. B 719, 29 (2013)
G. Aad et al. (ATLAS), Phys. Rev. Lett. 110, 182302 (2013)
R. Aaij et al. (LHCb), Phys. Lett. B 762, 473 (2016)
S. Chatrchyan et al. (CMS), Phys. Lett. B 724, 213 (2013)
B.B. Abelev et al. (ALICE), Phys. Rev. C 90, 054901 (2014)
M. Aaboud et al. (ATLAS), Eur. Phys. J. C 77, 428 (2017)
B. Alver et al. (PHOBOS), Phys. Rev. Lett. 104, 062301 (2010)
B.I. Abelev et al. (STAR), Phys. Rev. C 80, 064912 (2009)
A. Adare et al. (PHENIX), Phys. Rev. Lett. 114, 192301 (2015)
L. Adamczyk et al. (STAR), Phys. Lett. B 747, 265 (2015)
A. Adare et al. (PHENIX), Phys. Rev. Lett. 115, 142301 (2015)
C. Aidala et al. (PHENIX), Nat. Phys. 15, 214 (2019)
The ATLAS collaboration (ATLAS), Twoparticle azimuthal correlations in photonuclear ultraperipheral Pb+Pb collisions at 5.02\(\sim \)TeV with ATLAS (2019)
A. Badea, A. Baty, P. Chang, G.M. Innocenti, M. Maggi, C. Mcginn, M. Peters, T.A. Sheng, J. Thaler, Y.J. Lee, Phys. Rev. Lett. 123, 212002 (2019)
I. Abt et al. (ZEUS), (2019)
S. Schlichting, P. Tribedy, Adv. High Energy Phys. 2016, 8460349 (2016)
C. Loizides, Nucl. Phys. A 956, 200 (2016)
B. Schenke, Nucl. Phys. A 967, 105 (2017)
J.L. Nagle, W.A. Zajc, Ann. Rev. Nucl. Part. Sci. 68, 211 (2018)
Z. Citron et al., Report from Working Group 5: Future Physics Opportunities for HighDensity QCD at the LHC with Heavyion and Proton Beams, vol. 7, pp. 1159–1410 (2019)
S. Jeon, U. Heinz, Introduction to Hydrodynamics, pp. 131–187 (2016)
P. Romatschke, U. Romatschke, Relativistic Fluid Dynamics In and Out of Equilibrium, Cambridge Monographs on Mathematical Physics (Cambridge University Press, Cambridge, 2019) (ISBN 9781108483681, 9781108750028)
P. Romatschke, Eur. Phys. J. C 77, 21 (2017)
L. Keegan, A. Kurkela, P. Romatschke, W. van der Schee, Y. Zhu, JHEP 04, 031 (2016)
B. Schenke, P. Tribedy, R. Venugopalan, Phys. Rev. Lett. 108, 252301 (2012)
A. Kurkela, W. van der Schee, U.A. Wiedemann, B. Wu, Phys. Rev. Lett. 124, 102301 (2020)
A. Kurkela, U.A. Wiedemann, B. Wu, Eur. Phys. J. C 79, 965 (2019)
A. Dumitru, F. Gelis, L. McLerran, R. Venugopalan, Nucl. Phys. A 810, 91 (2008)
E. Iancu, A. Leonidov, L. McLerran, The color glass condensate: an introduction, in QCD Perspectives on Hot and Dense Matter. Proceedings, NATO Advanced Study Institute, Summer School, Cargese, August 6–18, 2001, pp. 73–145 (2002)
L. McLerran, The color glass condensate and glasma (2008). arXiv:0804.1736 [hepph]
F. Gelis, E. Iancu, J. JalilianMarian, R. Venugopalan, Ann. Rev. Nucl. Part. Sci. 60, 463 (2010)
Y.V. Kovchegov, E. Levin, Quantum Chromodynamics at High Energy, vol. 33 (Cambridge University Press, Cambridge, 2012). http://www.cambridge.org/de/knowledge/isbn/item6803159(ISBN 9780521112574, 9780521112574, 9781139557689)
A. Kovner, M. Lublinsky, Int. J. Mod. Phys. E 22, 1330001 (2013)
A. Kovner, M. Lublinsky, Phys. Rev. D 83, 034017 (2011)
A. Kovner, M. Lublinsky, Phys. Rev. D 84, 094011 (2011)
N. Armesto, L. McLerran, C. Pajares, Nucl. Phys. A 781, 201 (2007)
A. Dumitru, K. Dusling, F. Gelis, J. JalilianMarian, T. Lappi, R. Venugopalan, Phys. Lett. B 697, 21 (2011)
E. Levin, A.H. Rezaeian, Phys. Rev. D 84, 034031 (2011)
E. Iancu, A.H. Rezaeian, Phys. Rev. D 95, 094003 (2017)
V.A. Abramovsky, N.V. Prikhod’ko, Dipole amplitude correlation in saturation model beyond mean field approximation (2005). arXiv:hepph/0512343
J.L. Albacete, P. GuerreroRodríguez, C. Marquet, JHEP 01, 073 (2019)
A. Krasnitz, R. Venugopalan, Phys. Rev. Lett. 84, 4309 (2000)
A. Krasnitz, R. Venugopalan, Phys. Rev. Lett. 86, 1717 (2001)
T. Lappi, Phys. Rev. C 67, 054903 (2003)
A. Dumitru, V. Skokov, Phys. Rev. D 91, 074006 (2015)
A. Dumitru, L. McLerran, V. Skokov, Phys. Lett. B 743, 134 (2015)
A. Dumitru, A.V. Giannini, V. Skokov, Anisotropic particle production and azimuthal correlations in highenergy pA collisions (2015). arXiv:1503.03897 [hepph]
K. Dusling, R. Venugopalan, Phys. Rev. Lett. 108, 262001 (2012)
K. Dusling, R. Venugopalan, Phys. Rev. D 87, 051502 (2013)
K. Dusling, R. Venugopalan, Phys. Rev. D 87, 054014 (2013)
K. Dusling, R. Venugopalan, Phys. Rev. D 87, 094034 (2013)
K. Dusling, M. Mace, R. Venugopalan, Phys. Rev. Lett. 120, 042002 (2018)
K. Dusling, M. Mace, R. Venugopalan, Phys. Rev. D 97, 016014 (2018)
T. Altinoluk, N. Armesto, G. Beuf, A. Kovner, M. Lublinsky, Phys. Lett. B 751, 448 (2015)
B. Blok, C.D. Jakel, M. Strikman, U.A. Wiedemann, JHEP 12, 074 (2017)
B. Blok, U.A. Wiedemann, Phys. Lett. B 795, 259 (2019)
T. Altinoluk, N. Armesto, G. Beuf, A. Kovner, M. Lublinsky, Phys. Rev. D 95, 034025 (2017)
Y.V. Kovchegov, D.E. Wertepny, Nucl. Phys. A 906, 50 (2013)
Y.V. Kovchegov, D.E. Wertepny, Nucl. Phys. A 925, 254 (2014)
T. Altinoluk, N. Armesto, G. Beuf, A. Kovner, M. Lublinsky, Phys. Lett. B 752, 113 (2016)
S. Ozonder, Phys. Rev. D 91, 034005 (2015)
S. Ozonder, Turk. J. Phys. 42, 78 (2018)
T. Altinoluk, N. Armesto, A. Kovner, M. Lublinsky, Eur. Phys. J. C 78, 702 (2018)
T. Lappi, B. Schenke, S. Schlichting, R. Venugopalan, JHEP 01, 061 (2016)
T. Altinoluk, N. Armesto, D.E. Wertepny, JHEP 05, 207 (2018)
L.D. McLerran, R. Venugopalan, Phys. Rev. D 49, 2233 (1994)
L.D. McLerran, R. Venugopalan, Phys. Rev. D 50, 2225 (1994)
M. Martinez, M.D. Sievert, D.E. Wertepny, JHEP 07, 003 (2018)
M. Martinez, M.D. Sievert, D.E. Wertepny, JHEP 02, 024 (2019)
M.K. Davy, C. Marquet, Y. Shi, B.W. Xiao, C. Zhang, Nucl. Phys. A 983, 293 (2019)
C. Zhang, C. Marquet, G.Y. Qin, S.Y. Wei, B.W. Xiao, Phys. Rev. Lett. 122, 172302 (2019)
E. Gotsman, E. Levin, Eur. Phys. J. C 78, 792 (2018)
A. Kovner, M. Lublinsky, V. Skokov, Phys. Rev. D 96, 016010 (2017)
Y.V. Kovchegov, V.V. Skokov, Phys. Rev. D 97, 094021 (2018)
M. Mace, V.V. Skokov, P. Tribedy, R. Venugopalan, Phys. Lett. B 788, 161 (2019). [Erratum: Phys. Lett. B 799, 135006 (2019)]
M. Mace, V.V. Skokov, P. Tribedy, R. Venugopalan, Phys. Rev. Lett. 121, 052301 (2018). [Erratum: Phys. Rev. Lett. 123, 039901 (2019)]
P. Agostini, T. Altinoluk, N. Armesto, Eur. Phys. J. C 79, 600 (2019)
P. Agostini, T. Altinoluk, N. Armesto, Eur. Phys. J. C 79, 790 (2019)
T. Altinoluk, N. Armesto, G. Beuf, M. Martinez, C.A. Salgado, JHEP 07, 068 (2014)
T. Altinoluk, N. Armesto, G. Beuf, A. Moscoso, JHEP 01, 114 (2016)
Y. MehtarTani, Phys. Rev. C 75, 034908 (2007)
T. Altinoluk, A. Dumitru, Phys. Rev. D 94, 074032 (2016)
U.A. Wiedemann, Nucl. Phys. B 588, 303 (2000)
M. Gyulassy, P. Levai, I. Vitev, Nucl. Phys. B 594, 371 (2001)
Y. MehtarTani, C.A. Salgado, K. Tywoniuk, JHEP 04, 064 (2012)
A. Capella, A. Krzywicki, Phys. Rev. D 18, 4120 (1978)
A. Capella, A. Krzywicki, E.M. Levin, Phys. Rev. D 44, 704 (1991)
N.S. Amelin, N. Armesto, M.A. Braun, E.G. Ferreiro, C. Pajares, Phys. Rev. Lett. 73, 2813 (1994)
M.A. Braun, C. Pajares, Eur. Phys. J. C 71, 1558 (2011)
C. Andres, A. Moscoso, C. Pajares, Phys. Rev. C 90, 054902 (2014)
I. Bautista, L. Cunqueiro, J.D. de Deus, C. Pajares, J. Phys. G 37, 015103 (2010)
A. Ortiz Velasquez, P. Christiansen, E. Cuautle Flores, I. Maldonado Cervantes, G. Paic, Phys. Rev. Lett. 111, 042001 (2013)
C. Bierlich, G. Gustafson, L. Lonnblad, Phys. Lett. B 779, 58 (2018)
C.B. Duncan, P. Skands, Fragmentation of two repelling lund strings (2019). arXiv:1912.09639 [hepph]
J.L. Nagle, R. Belmont, K. Hill, J. Orjuela Koop, D.V. Perepelitsa, P. Yin, Z.W. Lin, D. McGlinchey, Phys. Rev. C 97, 024909 (2018)
J.D. Bjorken, S.J. Brodsky, A. Scharff Goldhaber, Phys. Lett. B 726, 344 (2013)
S.D. Glazek, S.J. Brodsky, A.S. Goldhaber, R.W. Brown, Phys. Rev. D 97, 114021 (2018)
K. Boreskov, A. Kaidalov, O. Kancheli, Eur. Phys. J. C 58, 445 (2008)
E. Gotsman, E. Levin, U. Maor, S. Tapia, Phys. Rev. D 93, 074029 (2016)
E. Gotsman, E. Levin, U. Maor, Eur. Phys. J. C 76, 607 (2016)
C.Y. Wong, Phys. Rev. C 84, 024901 (2011)
T.A. Trainor, J. Phys. G 40, 055104 (2013)
A. Dumitru, T. Lappi, V. Skokov, Phys. Rev. Lett. 115, 252301 (2015)
Y. Hatta, B.W. Xiao, F. Yuan, Phys. Rev. Lett. 116, 202301 (2016)
T. Altinoluk, N. Armesto, G. Beuf, A.H. Rezaeian, Phys. Lett. B 758, 373 (2016)
A. Dumitru, V. Skokov, Phys. Rev. D 94, 014030 (2016)
H. Mantysaari, N. Mueller, B. Schenke, PoS DIS2019, 060 (2019)
T. Altinoluk, R. Boussarie, JHEP 10, 208 (2019)
B. Schenke, S. Schlichting, P. Tribedy, R. Venugopalan, Phys. Rev. Lett. 117, 162301 (2016)
J.L. Nagle, W.A. Zajc, Phys. Rev. C 99, 054908 (2019)
S. Vigdor, The Case for Continuing RHIC Operations (2012). https://www.bnl.gov/npp/docs/The%20Case%20for%20Continuing%20RHIC%20Operations%20_draft%205.pdf
A. Accardi et al., Eur. Phys. J. A 52, 268 (2016)
J. Abelleira Fernandez et al. (LHeC Study Group), J. Phys. G 39, 075001 (2012)
A. Abada et al. (FCC), Eur. Phys. J. ST 228, 755 (2019)
Acknowledgements
We would like to Ian Balitsky, Boris Blok, Giovannni Chirilli, Adrian Dumitru, Cyrille Marquet, Mauricio Martinez, Larry McLerran, Alfred Mueller, Jamie Nagle, Carlos Pajares, Matt Siervert, Vladimir Skokov, Raju Venugopalan, Douglas Wertepny, Urs Wiedemann and William Zajc, and the participants of the 2018 Benasque workshop on Collectivity and correlations in highenergy hadron and nuclear collisions and the COST workshop on collectivity in small systems, for collaborations and many fruitful discussions. We would like to specially thank Pedro Agostini, Guillaume Beuf, Alex Kovner and Michael Lublinsky for a longstanding collaboration. NA is supported by Ministerio de Ciencia e Innovación of Spain under projects FPA201783814P and Unidad de Excelencia María de Maetzu under project MDM20160692, by the European research Council under project ERC2018ADG835105 YoctoLHC, by Xunta de Galicia under project ED431C 2017/07, by Consellería de Educación, Universidade e Formación Profesional as Centro de Investigación do Sistema universitario de Galicia (ED431G 2019/05), and by FEDER. TA is supported by Grant No. 2018/31/D/ST2/00666 (SONATA 14  National Science Centre, Poland). This work has been performed in the framework of COST Action CA 15213 “Theory of hot matter and relativistic heavyion collisions” (THOR), MSCA RISE 823947 “Heavy ion collisions: collectivity and precision in saturation physics” (HIEIC) and has received funding from the European Union’s Horizon 2020 research and innovation programme under Grant Agreement No. 824093.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Laura Tolos
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Altinoluk, T., Armesto, N. Particle correlations from the initial state. Eur. Phys. J. A 56, 215 (2020). https://doi.org/10.1140/epja/s10050020002256
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epja/s10050020002256