CGC/saturation approach for soft interactions at high energy: long range rapidity correlations

In this paper we continue our program to construct a model for high energy soft interactions that is based on the CGC/saturation approach. The main result of this paper is that we have discovered a mechanism that leads to large long range rapidity correlations and results in large values of the correlation function Ry1,y2≥1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R\left( y_1,y_2\right) \ge 1$$\end{document}, which is independent of y1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$y_1$$\end{document} and y2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ y_2$$\end{document}. Such a behavior of the correlation function provides strong support for the idea that at high energies the system of partons that is produced is not only dense but also has strong attractive forces acting between the partons.

In this paper we continue our effort [16][17][18] to comprehend such interactions, by constructing a model that incorporates the advantages of two theoretical approaches to high energy QCD.
On the other hand, there exists a different approach to high energy QCD: the BFKL Pomeron [55][56][57][58] and its interactions [19,, which is suitable to describe diffractive physics. The BFKL Pomeron calculus turns out to be close to the old Reggeon theory [80], so, for calculating the inclusive characteristics of multiparticle production, we can apply the Mueller diagram technique [81]. The relation between these two approaches has not yet been established, but they are equivalent [74][75][76][77] for the rapidities (ln (s/s 0 )), such that where BFKL denotes the intercept of the BFKL Pomeron.
As we have discussed [16], the parameters of our model are such that for Y ≤ 36, we can trust our approach, based on the BFKL Pomeron calculus. This paper is the next step in our program to build a model for high energy soft scattering, based on an analytical calculation, without using a Monte Carlo simulation. We discuss the correlation function: where Y denotes the total rapidity (Y = ln (s/s 0 ) and s = W 2 , W is the energy in c.m.f.) and y 1 and y 2 are the rapidities of the produced hadrons. σ N SD = σ tot − σ el − σ sd − σ dd where σ tot (σ el , σ sd , σ dd ) are the total, elastic, single, and double diffraction cross sections. The paper is organized as follows. In the following section we discuss the main features of our approach, concentrating on the description of diffractive processes. In Sect. 3, we derive the main formulas for the correlation functions in our approach, while in Sect. 4 we compare our predictions with the available experimental data.

Our model: generalities, elastic amplitude, and inclusive production
In this section we briefly review our model, which successfully describes diffractive [16,17] and inclusive cross sections [18]. The main ingredient of our model is the BFKL Pomeron Green's function, which we determined using the CGC/saturation approach [16,78]. We determined this function from the solution of the non-linear Balitsky-Kovchegov equation [25,26,29], using the MPSI approximation [79,[82][83][84][85] to sum enhanced diagrams shown in Fig. 1a. It has the following form: In these formulas we take a = 0.65; this value was chosen so as to obtain the analytical form for the solution of the BK equation. The parameters λ and φ 0 can be estimated in the leading order of QCD, but due to large next-to-leading order corrections we treat them as parameters of the fit. m is a nonperturbative parameter, which characterizes the large impact parameter behavior of the saturation momentum, as well as the typical sizes of the dipoles that take part in the interactions. The value of m = 5.25 GeV in our model justifies our main assumption that BFKL Pomeron calculus based on a perturbative QCD approach is able to describe soft physics, since m μ soft , where μ soft denotes the natural scale for soft processes (μ soft ∼ QCD and/or the pion mass).
Unfortunately, since the confinement problem is far from being solved, we assume a phenomenological approach for the structure of the colliding hadron. We use a two channel model, which also allows us to calculate the diffractive production in the region of small masses. In this model, we replace the rich structure of the diffractively produced states by a single state with the wave function ψ D , à la Good and Walker [86]. The observed physical hadronic and diffractive states are written in the form where α 2 + β 2 = 1. (2. 3) The functions ψ 1 and ψ 2 form a complete set of orthogonal functions {ψ i } which diagonalize the interaction matrix T The unitarity constraints take the form where G in i,k denotes the contribution of all non-diffractive inelastic processes, i.e. it is the summed probability for these final states to be produced in the scattering of a state i off a state k. In Eq. (2.5) √ s = W denotes the energy of the colliding hadrons, and b the impact parameter. A simple solution to Eq. (2.5) at high energies, has the eikonal form with an arbitrary opacity ik , where the real part of the amplitude is much smaller than the imaginary part. We have , is the probability that the initial projectiles (i, k) reach the final state interaction unchanged, regardless of the initial state rescatterings.
Note that there is no factor 1/2; its absence stems from our definition of the dressed Pomeron. Table 1 Fitted parameters of the model. The values are taken from Ref. [17] Model  In the eikonal approximation we We propose a more general approach, which takes into account new small parameters, which come from the fit to the experimental data (see Table 1 and Fig. 1): The second equation in Eq. (2.9) leads to the fact that b in Eq. (2.8) is much smaller than b and b ; therefore, Eq. (2.8) can be re-written in a simpler form: (2.10) Selecting the diagrams using the first equation in Eq. (2.9) indicates that the main contribution stems from the net diagrams shown in Fig. 1b. The sum of these diagrams [17] leads to the following expression for i,k (s, b): where T (s, b) is given by Eq. In all previous formulas, the value of the triple BFKL Pomeron vertex is known: G 3P = 1.29 GeV −1 .
To simplify the further discussion, we introduce the notation with a = 0.65. It should be stressed that Eq. (2.13) with this value of a is an analytical approximation to the numerical solution for the BK equation [16,17,78]. Therefore, a is not a fitting parameter, as it has been determined from CGC/saturation effective theory. We have G i P (Y, b) = g i (b)G dressed (T ). We recall that the BK equation sums the 'fan' diagrams shown in Fig. 2.
For the elastic amplitude we have (2.14) To determine the correlation function (given in Eq. (1.2)), we need to know the single inclusive cross sections. We have discussed these cross sections in Ref. [18]; for the sake of completeness we give the formula that describes the Mueller diagram of Fig. 3.
We have where Y denotes the total rapidity of the colliding particles, and y is the rapidity of produced hadron. I n (i) is given by a PP is a fitted parameter that was determined in Ref. [18] (see Table 1).

Correlations between two parton showers
The Mueller diagram for the correlations between two parton showers is shown in Fig. 4. Examining this diagram, we see that the contribution to the double inclusive cross section differs from the product of two single inclusive cross sections. There are two reasons for this; the first is that in the expression for the double inclusive cross section, we integrate the product of the single inclusive cross sections, over b, at fixed b. The second is that the summation over i and k for the product of single inclusive cross sections is for fixed i and k. Introducing the following new function enables us to write the analytical expression: Using Eq. (3.1) we can write the double inclusive cross section in the form d 2 σ 2 parton showers dy 1 dy 2 2) with Eq. (2.15) squared, we note the different powers of α and β, which reflect the different summation over i and k, as well as different integration over b.

Correlations in one parton shower: semi-enhanced diagrams
The main theoretical assumption that we make in calculating the correlation in a one parton shower, is that the Mueller diagram technique [81], and the AGK cutting rules [87] are valid. We should note, however, that even if the Mueller diagrams provide the correct description of inclusive processes in QCD, the AGK cutting rules are not valid for calculations of the correlations in QCD [88][89][90]. Nevertheless, we believe that we can neglect the AGK cutting rules violating contributions since, first, they do not lead to long range rapidity correlations, which are the main subject of our concern, and second, as we will show below, the correlations in one parton shower turn out to be negligibly small. It is instructive to write the expression for the first Mueller diagram in the following form (see Fig. 6a): The expression for the first Mueller diagram for two parton showers correlation (see Fig. 4a) has the form Therefore, we could expect that the semi-enhanced diagrams give a larger contribution to the double inclusive cross section than the production from two parton showers. However, Fig. 5 shows that both the value and the increase turn out to be small in the kinematic region accessible to experiment. Even at ultra high energies, shown in Fig. 5b, On the other hand, the contribution of Eq. (3.5) is small and is proportional to G 3P /g i 1. Bearing in mind that G 3P = 1.29 in our approach, one can see that the maximum value for max (H) ≈ 0.012 and we have the values of Therefore, we expect that the contribution of the correlations due semi-enhanced diagrams is negligibly small.
The general expression for the double inclusive cross section (see Fig. 6b) can be written using two new functions S i (y 1 , y 2 ) and S k (y 1 , y 2 ) defined as It takes the form The first Mueller diagram for the correlations from the enhanced diagram is shown in Fig. 7a, and it has the following form: (3.10) this can be re-written as d 2 σ one parton shower enh dy 1 dy 2 An example of typical diagrams is shown in Fig. 7. The formula, summing all diagrams shown in Fig. 7b where G dressed (Y, b) is determined by Eqs. (2.1) and (2.2). The contributions of enhanced diagrams are proportional to the square of the ratios given by Eq. (3.6) and, therefore, they are negligibly small.

Correlation function
The correlation function R (y 1 , y 2 ) is defined as

Kinematic corrections
In all our previous equations we assumed that Y = ln W 2 /W 2 0 with W 0 = 1 GeV. This assumption appears natural for the elastic amplitude and the cross section of the single inclusive production, but it should be re-examined for the correlation function. For this observable, the definition of Y has to be modified to account for the fact that the energy of the parton shower is not equal to W = √ s (s = W 2 ) but smaller than or equal tos =W 2 = x 1 x 2 s = x 1 x 2 W 2 (see Fig. 8, where we show the diffractive cut of the Mueller diagram of Fig. 4a). The simplest way to find x 1 and x 2 is to assume that both p 2 1 = p 2 2 = −Q 2 μ 2 soft , where μ soft is the scale of the soft interactions, μ soft ∼ QC D . In our approach the scale of the hardness for the BFKL Pomeron is  The red line indicates the moment C 2 = |n 2 | / |n 2 | = 2 for the window in rapidity −0.5 < y < +0.5, measured by the CMS experiment [13] Q = m μ soft . Bearing this in mind, the energy variable x 1 (x 2 ) for gluon-hadron scattering is equal to 0 = (x 1 P 1 + p 1 ) 2 = −Q 2 + x 1 2 p 1 · P 1 ; For the two channel model, it is the mass of the diffractive state. We can use the quark structure function to estimate the typical value of x 1 = x 2 , as is suggested in Ref. [91,92]. Using the structure functions at Q 2 ≈ 25 GeV 2 , one finds that |x 1 | ≈ 0.3 ÷ 0.5. In Fig. 9 the values of R (Y, 0, 0) are plotted for Y = ln x 2 1 W 2 /W 2 0 as a function of x 1 .

Correlation in one parton shower: emission from one BFKL Pomeron
In addition to the sources of correlation that have been discussed above, we need to take into account the correlation between two gluons emitted from one BFKL Pomeron (see Fig. 10). At large y 12 = |y 1 − y 2 | the diagram of Fig. 10 induces long range correlations in rapidity, however, at small y 12 this emission is suppressed, and we do not expect a large contribution from this source. The contribution of this diagram can be written in the form For small values of the argument i (y) . (3.18), has no dependence on y 12 , leading to long range rapidity correlations. However, it turns out that the exact computation, leads to very small values of C BFKL (y 1 , y 2 ): approximately 0.2 ÷ 0.4 % of the contributions from the sources discussed above.

Short range rapidity correlation
Besides long range rapidity correlations, the emission from one BFKL Pomeron, as well as the hadronization in one gluon jet, can lead to short range correlations in rapidity. Unfortunately, at present, this contribution cannot be treated on pure theoretical grounds, as it involves confinement effects. To estimate this contribution, we introduce the Mueller diagram shown in Fig. 11, where we describe this correlation by the phenomenological constant a PR , and introduce the correlation length ≈ 2. In the diagram of Fig. 11 for the zigzag line we have a 2 PR exp − y 12 (y 12 = |y 1 = y 2 |). Our estimate for stems from Reggeon phenomenology, in which the zigzag line describes the contribution of the secondary Reggeon, with a propagator exp (−(1 − α R (0))y 12 ) and α R (0) ≈ 0.5. At first sight such a parametrization of the short range correlations neglects the contributions of resonances, which are an essential part of these correlations. Due to the duality between the Reggeon exchange and the resonance contributions, we believe that the Reggeon description (3.20) Figure 9, shows that the correlation function increases with energy and becomes rather large (of the order of 1) at W = 7 TeV. This qualitative feature is in agreement with the experimental data from the LHC. The first set of data is the multiplicity distribution measured by the CMS collaboration [13]. In particular, the value of C 2 = |n 2 | / |n 2 | turns out to be very close to 2, for the window in rapidities −0.5 < η < 0.5. Since C 2 = R (0, 0) + 1 + 1/dN /dη| η=0 where dN /dη η=0 denotes the multiplicity at η = 0, and at W = 7 TeV, it is equal to 5.8, while R (0, 0) = 0.82. 1 The second set of the data is the measurement of the double parton interaction (DPI) [93][94][95][96][97][98]. In the LHC experiments, the double inclusive cross sections of two pairs of back-toback jets with momenta p T,1 and p T,2 were measured with rapidities of two pairs (y 1 and y 2 ) close to each other (y 1 ≈ y 2 ). These pairs can only be produced from two different parton showers. The data were parameterized in the form dσ dy 1 d 2 p T,1 dy 2 d 2 p T,2 = m 2σ eff dσ dy 1 d 2 p T,1 dσ dy 2 d 2 p T,2 (4.1) where m = 2 for pairs of different jets, and m = 1 for identical pairs. One can calculate the rapidity correlation function using Eq. (4.1),

Predictions and comparison with the experiment
For the above the estimates we use σ eff = 12 − 15 mb (see Refs. [93][94][95][96][97][98]) and σ in = σ tot − σ el − σ sd − σ dd ≈ 50 mb for Hence, at fixed n/n, the data at W = 13 TeV with higher multiplicity is above the data at W = 7 TeV the energy W = 7 TeV (see Ref. [17] and references therein). These data confirm that at high energies we are dealing with a system of partons that have a large mutual attraction. The fact that we predict a smaller correlation than we have in this experiment does not discourage us, since the correlation function in Eq. (4.2) differs from the one that we calculate (see Eq. (1.2)).
The forward-backward correlation has been measured by the ATLAS collaboration in Ref. [14]. The observable that was used in Ref. [14] differs from the correlation function R (η 1 , η 2 ), and it can be re-written as The value of ρ n f b ∼ 0.666 [14] indicates large correlations, but it is difficult to compare ρ n f b with our estimates, since ATLAS introduced a specific selection: the p T of all produced particles should be larger than 100 MeV, while R (η 1 , η 2 ) is defined as integrated over all momenta.
Using R (0, 0), we estimate the values for C 2 for W = 13 TeV, and using the formula for the negative binomial distribution, we obtain the multiplicity distribution in the rapidity window −0.5 < η < 0.5 shown in Fig. 12. In Eq. (4.4) n = dN /dη| η=0 , which was calculated in Ref. [18]. From Fig. 12 we expect a violation of the KNO scaling behavior [99]. Accordingly to KNO scaling σ n /σ = F (n/n) with n = 0.5 −0.5 dη dN /dη = dN /dη| η=0 . It turns out that R (y 1 , y 2 ) at fixed energy depends neither on y 1 nor on y 2 , giving a perfect example of long range correlations in rapidity. To understand why we have these features, it is instructive to start from Fig. 4a at small values of Y . In this kinematic region we can replaceG dressed → G bare =T (Y ) and N BK (Y, b) → g i (b)G bare . After simple algebra, the correlation function is equal to (4.5) In Eq. (4.5) we use the fact that in our model g 1 g 2 . Equation (4.5) leads to a correlation function that does not depend on y 1 and y 2 .
On the other hand, at very large is a step function. Plugging in this simple expression, we obtain R (Y ) in Eq. (4.6) denotes a typical impact parameter at large Y , which is proportional to Y . 2 Recall that at high energies, all components of the wave functions in the two channel model, give the same contribution. This is the reason why we do not have an extra factor, depending on α and β. Equation (4.6) shows the logarithmic dependence on y 1 . Using Eq. (4.6) we can estimate the y 1 (y 2 ) dependence of R (y 1 , y 2 )), calculating .  Fig. 13, in which we plot the results of our calculation, one can see that only at large Y > 11 does |b 2 (Y ) | start showing visible Y dependence. Two vertical dotted lines mark the widow in rapidity, which is essential in the calculation of the correlation function at W = 7 TeV for −2 < y 1 , y 2 < 2. We can expect a change of R by 2 %. The actual calculation gives even less.
Using this relation, we estimate R (y 1 , y 2 ) as At W = 7 TeV from Eq. (4.8) we find that R (y 1 , y 2 ) = 1.647, while the exact calculation give 1.437 (see Fig. 9). At W = 13 TeV this simple formula leads to R (y 1 , y 2 ) = 1.72, versus 1.64 from the exact calculations (see Fig. 9). The correlations, measured by the ATLAS collaboration [14], at first sight contradict both our estimate and the CMS data, regarding the multiplicity distribution. We first check Eq. (4.3). The measured observable has the form [14] The numerator of Eq. (4.9) can be written as R (y 1 , y 2 ) y 1 y 2 , where y i is the interval of rapidities where the hadrons are measured. However, at the same value of rapidity R (y 1 , y 1 ) ( y 1 ) 2 corresponds to n(n − 1) − n 2 . Therefore, we have the expression which leads to the following formula for ρ n f b : (4.10) Taking y 1 = 0.5, we see that the first element of the Table  2 is equal to 0.7, which is in good agreement with the experimental value 0.666 ± 0.011. To describe our results given in Table 2, we also need to take short range rapidity correlations into account. In this table, in parentheses we have our estimates, which we obtain describing the correlation function in the form R (y 1 , y 2 ) = R long range (y 1 , y 2 ) with a = 0.7 and = 2. In Eq. (4.11) y 12 = |y 1 − y 2 | and we restrict ourselves to the contribution of the state "1" in Eq. (3.20), since g 1 g 2 . We assumed that C short range (0, 0) = 0 since at y 12 = 0 Eq. (3.20) leads to the long range correlations which we have calculated in Sect. 3.4. One can see that the agreement is not perfect, but it demonstrates that the ATLAS data can be reproduced, by including the short range correlations.

Conclusions
The main result of this paper, is that in our model, which is based on the CGC/saturation approach, we have discovered a mechanism that produces large, long range rapidity correlations at high energies. The large values of the correlation function R (y 1 , y 2 ) ≥ 1 at high energies, lends strong support to the idea that at high energies the system of partons that is produced, is not only dense, but also has strong attractive forces acting between the partons.
Another interesting result is that this mechanism, leads to long range rapidity correlations that are independent of y 1,2 , and stems from the correlation between two parton showers, i.e., we found that in our approach the contribution from enhanced and semi-enhanced diagrams that describe the interaction inside a single parton shower, turns out to be negligibly small. It should be stressed that the contribution of the same diagrams in the process of diffraction dissociation in the region of large mass is essential, as they are responsible for the experimental measured diffraction in the region Table 2 Multiplicity correlations for events at √ s = 7 TeV for events with a minimum of two charged particles in the kinematic interval p T > 100 MeV and |η| < 2.5 for different combinations of forward and backward pseudorapidity interval. The data is take from Ref. [14]. The numbers in parentheses are our estimates of large mass [17]. This feature, is in direct contradiction to the soft Pomeron-based models for the soft scattering amplitude, suggest should cancel (see Ref. [92]). In these models, the correlations from one parton shower are larger than from two parton showers, and they led to the y 1,2 dependence. In other words, in the soft Pomeron based model, the contribution of semi-enhanced and enhanced Pomeron diagrams that describe large mass diffraction, is sufficiently large to induce considerable (actually dominant) contributions to the two particle correlation. Scrutinizing our formulas, we found that the main reason for the small size of the correlation in one parton showers that we observed in our approach, stems from the most theoretically reliable part of our model: from the expression for the dressed BFKL Pomeron Green function.
Another result, our model predicts is the probability of two parton showers production without introducing new phenomenological parameters, 3 which is able to describe the long range rapidity correlations.
We demonstrated that our model is able to describe the LHC data, eminating from the CMS and ATLAS collaborations. These data are certainly insufficient for a thorough analysis of the details of our approach, but they confirm that the long range rapidity correlations are large at high energies. Our prediction for W = 13 TeV is shown in Fig. 9. The correlations should increase with the energy, and the measurements at W = 13 TeV should clarify the situation.
It worthwhile mentioning that the main contributions to the correlation function in terms of Mueller diagrams, has been known from 1970s, as well as the fact that non-enhanced diagrams (production of two and more parton showers) led to long range rapidity correlations, and were needed to describe the correlation at ISR and TEVATRON energies. However, the problem of the contributions of semi-enhanced and enhanced diagrams in the framework of the soft Pomeron calculus, has not been solved theoretically. This, as well as the estimates of the value of the contribution of these diagrams from the experimental data, are very dependent on the model of the Pomeron interactions. From the first paper on CGC/saturation approach [19], it became clear that hard processes could give a considerable contribution to the correlation function. However, it was not clear whether such correlations will be enough to reproduce the experimental data.
It was already appreciated in 1987 that taking into account the perturbative QCD sub-collisions in the Monte Carlo models is not sufficient to reproduces the correlation data from the UA5 experiment [100], and that multiple sub-collisions should be added. In the language of Ref. [100], our main result can be formulated as the statement that only hard and soft processes from two sub-collisions contribute to the correlation function.
In closing, we suggest that the investigation of short range rapidity correlations in MC models, which do not depend on Pomeron interactions and on the contributions of the enhanced and semi-enhanced diagrams, can provide more insight, as such estimates will be complementary to our Reggeon motivated description of short range rapidity correlations.