A model for strong interactions at high energy based on the CGC/saturation approach

We present our first attempt to develop a model for soft interactions at high energy, based on the BFKL Pomeron and the CGC/saturation approach. We construct an eikonal-type model, whose opacity is determined by the exchange of the dressed BFKL Pomeron. The Green function of the Pomeron is calculated in the framework of the CGC/saturation approach. Using five parameters we achieve a reasonable description of the experimental data at high energies (W≥0.546\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W\ge 0.546$$\end{document} TeV) with overall χ2/d.o.f.≈2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2/d.o.f. \approx 2$$\end{document}. The model results in different behavior for the single- and double-diffraction cross sections at high energies. The single-diffraction cross section reaches a saturated value (about 10 mb) at high energies, while the double-diffraction cross section continues growing slowly.


Introduction
The strong interaction at high energies is one of the most difficult and unrewarding problems of high energy physics. The reason for this is the embryonic stage of our understanding of non-perturbative QCD. Traditionally, we consider the strong interaction at high energy as a typical example of processes that occur at long distances, where the unknown force confining quarks and gluons plays a crucial role, making all our theoretical efforts to treat these processes fruitless. The description of these processes which we need for practical purposes is the field of high energy phenomenology, based on Pomeron calculus [1][2][3][4]. The LHC data [5][6][7][8][9][10] shows that in many cases models based on this phenomenology failed to agree with the results of the classical set of soft interaction data: total, elastic, and diffraction cross section as well as elastic slope and the inclusive production of the secondary hadrons [11][12][13][14][15][16][17][18]. However, there is a glimpse of hope due to the following two facts: models that fit the LHC data have been proposed [19][20][21]; and after two decades of experience in high energy phenomenology we have learned that the more theoretically based the phenomenological input is, the more appropriate and apprehensible, the resulting description of the data we obtain.
In Ref. [22] we reviewed our model which describes successfully all high energy data, including those at the LHC, Fig. 1 The structure of the parton cascade for a fast dipole and its relation to the Pomeron interaction. Helical lines denote gluons. The wavy lines describe BFKL Pomerons. The blobs stand for triple Pomeron vertices 2 and which incorporates theoretical ingredients from N = 4 SYM [23][24][25][26][27][28][29][30] and from perturbative QCD [31][32][33][34][35][36][37][38][39]. In the present paper we improve this approach, by including more pertinent theoretical input. First, we introduce a more constructive meaning to our old idea [15,40] that there is only one Pomeron that describes both soft and hard interactions. In perturbative QCD the BFKL Pomeron at high energy takes the following form [31,32] denotes the BFKL Pomeron Green function,ᾱ S the QCD coupling, r and R are the sizes of the two interacting dipoles. Y = ln s, where s = W 2 . W denotes the energy of the interaction, and b the impact parameter for the scattering amplitude of the two dipoles.
From Eq. (1.1) it is obvious that the BFKL Pomeron is not a pole in angular momentum but a branch cut, since its Y-dependence has an additional ln s term; it does not reproduce the exponential decrease at large b, which follows from the general properties of analyticity and unitarity [41,42]; as the exchange of the BFKL Pomeron depends on the sizes of the dipoles, consequently, the BFKL Pomeron does not factorize.
The Pomeron that appears in N = 4 SYM [23][24][25] corresponds to the BFKL Pomeron in QCD with the following glossary: Glossary ≡ AdS-CFT correspondence [26]: gives the Pomeron contribution at high energies for ln 2 ww * ≤ 4DY . For a wider kinematic region the contribution of the BFKL Pomeron can be found in Ref. [33].
denotes the QCD-like coupling.
Hereby we generalize our approach, dealing with the BFKL Pomeron instead of the simple Regge pole that was used in our previous model.
The second innovation is related to the Pomeron interaction. The LHC data supports the assumption that dense systems of partons (gluons) are produced in the proton-proton interactions at high energy. Such a system of partons appears naturally in the CGC/saturation approach [47][48][49][50][51][52][53][54][55][56][57][58][59][60][61][62][63][64][65][66], and it provides a successful description of the general properties of the average event at the LHC [67][68][69][70][71], and of the long range rapidity angular correlations [72][73][74]. In this paper, we use the CGC/saturation approach to describe the Pomeron interactions, replacing the Pomeron calculus. This strategy allows us not only to treat the Pomeron interactions, but also to include the saturation phenomenon, which was beyond the scope of our previous model. The parton cascade, which originates from the decay of one gluon to two gluons in QCD, can be described equally well in two ways. The first, using the QCD expression for this decay, describes the change of probability to have n gluons at rapidity Y , due to the decay of one gluon to two. The equation in this approach is a linear functional equation for the generating functional (see Refs. [53,54,[58][59][60][61][62][63][64][65][75][76][77]). The alternate way is to sum the Pomeron fan diagrams (see Fig. 1) in the framework of the BFKL Pomeron calculus [78][79][80][81][82][83].
In this paper we use the solution of the functional equation which was proposed and discussed in Ref. [84]. For completeness of presentation, we repeat the main ideas of the solution, and explain the physical meaning of the phenomenological parameters that we have introduced in our model.
The first simplification arises when we consider the interaction of the parton cascade with a large target (say with a heavy nucleus). In this case the functional equation reduces where Q s denotes the saturation scale [47][48][49][50][51][52]. The non-linear corrections are small, and the solution is the BFKL Pomeron; 2. r 2 Q 2 s (Y, b) ∼ 1 (vicinity of the saturation scale). The scattering amplitude has the following form [85][86][87]: where φ 0 is a constant, and where the critical anomalous dimension γ cr can be found from 1 (deeply inside the saturation domain). The amplitude approaches unity [88]: In spite of our understanding of all qualitative features of the solution, we do not have an equation for an analytical solution [66], which we need to reconstruct the parton cascade. The parton cascade can be described as the amplitude for the production of dipoles of size r i at impact parameters b i . This amplitude can be written as The solution to the non-linear equation is of the following general form: Unfortunately, we cannot find the coefficient C n , for the general non-linear equation. For the case of the simplified BFKL kernel (see Refs. [84,88]) the solution can be found, and we can suggest a simple formula that provides a very accurate solution of Eq. (2.5) (see Ref. [84]): with a = 0.65. This formula allows us to find C n (φ 0 , r ), and to reconstruct the amplitude of Eq. (2.4).

Summing Pomeron loops (MPSI approximation)
It was shown in Ref. [89] that in the BFKL Pomeron calculus for the parton cascade (see Fig. 1), the integration over rapidities of the triple Pomeron vertices, suggests that the value of the typical rapidity is of the order of Y − Y i ∼ 1/ BFKL . Consequently, only large Pomeron loops with rapidity of order Y , contribute at high energies [90][91][92][93]. To sum such loops we use the MPSI approximation developed in Ref. [90][91][92][93]. The essence of this approximation is to use the tchannel unitarity constraint, which is satisfied by the one BFKL Pomeron exchange. Indeed, at any value of Y , the BFKL Pomeron has the following property from t-channel unitarity [47,94] (see Fig. 2b): The MPSI approximation is illustrated in Fig. 2a, where the first non-trivial loop diagram is presented. This approximation enables us to evaluate the Pomeron loops, using the fan diagram structure of the parton cascade. The general MPSI equation for the sum of enhanced Pomeron diagrams, has the form which leads to a new Pomeron Green function (dressed Pomeron): In the last equation we used Eq. (2.6) with Since, for the proton-proton scattering, r = R, z > 0 (see Eq. (2.10)), we are dealing with parton cascades in the saturation domain. We recall that the saturation domain corresponds to z > 0, (r 2 Q 2 s ≥ 1), while z < 0(r 2 Q 2 s ≤ 1) characterizes the perturbative QCD region.
In Ref. [95][96][97] the MPSI approximation, as well as the equivalence of the CGC/saturation approach and the BFKL Pomeron calculus, was proven for a wide range of rapidities: For larger Y the MPSI approximation does not give the exact answer, since we have not introduced the vertex of the four Pomeron interaction, which violates the simple structure of the parton cascade shown in Fig. 1. The errors that stem from neglecting the four Pomeron interaction, have been evaluated in Ref. [22,[95][96][97].

Model: main formulas and parameters
In this section we describe our model. Its main ingredient is the sum of the Pomeron loops, which leads to a new dressed Pomeron Green function.

Dressed Pomeron
The resulting Green function of the Pomeron is given by Eq. (2.9). Using Eq. (2.7) and Eq. (2.9) (see Fig. 3a) we obtain the following expression (see Refs. [22,84] for more details): where (x) is the incomplete Euler gamma function (see 8.35 of Ref. [98]). The function T (Y − Y 0 , r, R, b) can be found from Eq. (2.9) and has the form where we used two inputs: r = R and Q 2 s = 1/(m 2 R 2 ) S (b) exp (λY ). The parameter λ =ᾱ S χ (γ cr ) /(1−γ cr ) in leading order of perturbative QCD. From phenomenology λ turns out to have the value λ = 0.2 ÷ 0.3 [99,101,102]. S (b) is a pure phenomenological profile function, which we choose to be of the form The parameter m represents the inverse size of the dipole m ∼ 1/r = 1/R. Unfortunately, we have no theoretical estimate for this mass. It may be large, reflecting the masses of glueballs and the small size of the typical dipole in a hadron [103]. Note that S (b) has a correct, exponential decrease at large b. This is an advantage of our approach, as it enables us to introduce a non-perturbative scale in a physical motivated way, for the observable that characterizes the principal property of the parton cascade. Therefore, in the framework of our approach we do not face the theoretical problem of large b behavior, which is the main unsolved problem in the CGC/saturation approach [104,105] (Fig. 3).

Interaction of dressed Pomerons
The interaction of a dressed Pomeron with a hadron is a non-perturbative problem, which cannot be solved for the moment. From the microscopic point of view this interaction depends on the size of a typical dipole in a hadron, on the probability of finding such a dipole, and on the interaction coupling. Since this interaction originates at long distances, we cannot calculate it even in the CGC/saturation approach. Introducing two phenomenological constants, g and m 1 , we describe the vertex of the hadron-Pomeron interaction as follows: To account for the interaction of the dressed Pomerons with hadrons, we use the strategy that has been suggested in Ref. [14], and which is based on the fact that we anticipate the value of g in Eq. (3.4) will be large. In this case we can evaluate the scattering amplitude in the following kinematic region of rapidities: The difference with our previous model reviewed in Ref. [22] lies in the value of G 3IP which was a phenomenological parameter, while now we are able to estimate it from the CGC/saturation approach. Finally, the opacity has the form (see Fig. 3b) In Eq. (3.6) we assumed that m m 1 . The factor 1.29 stem from estimates of the triple Pomeron vertex in the CGC/saturation approach.

Elastic amplitude
The elastic amplitude is (3.8)

Single diffraction
The cross section for single diffraction can be written as whereN SD (Y ) = d 2 b N SD Y, b and N SD (Y, b) has been calculated in Ref. [84]. It has the form . This definition of T is valid only in the region where T < 1. A more general formula is given in Ref. [84]. Equation (3.10) which sums the diagrams of Fig. 4a, where the double wavy lines crossed by the dashed one, denote the dressed Pomeron structure in terms of the produced particle. In Reggeon calculus it is referred to as the cut Pomeron. We would like to emphasis that in our approach this contribution is the solution to the equation for singlediffractive production of Ref. [108], which is given in Ref. [84]. The profile function for the single-diffraction production is taken from Eq. (3.25) of Ref. [15].

Double diffraction
The double diffraction cross section has the form can be determined from the simple expression where T = T (Y, b) of Eq. (3.2) with the same comments as for Eq. (3.10). Equation (3.11) sums the diagrams shown in Fig. 4b. 2 The profile S DD (b) is given by (3.13) 2 We thank our referee for pointing out that in our treatment of double diffraction, we have not included the contribution arising from the superposition of two (projectile and target) single diffraction processes, which may provide corrections to our present results. We will include this process in our planned two channel treatment.

Phenomenological parameters
In this section we summarize our phenomenological parameters and provide theoretical estimates for them. Altogether, we have five parameters: g, φ 0 , λ, m and m 1 .
• λ in the CGC/saturation approach, can be calculated in the leading order of perturbative QCD. It characterizes the energy dependence of the saturation scale in proton-proton collisions. Theoretical estimates give λ = α S χ (γ cr ) /(1−γ cr ) ≈ 4.88ᾱ S in leading order of perturbative QCD, where γ cr = 0.37. However, the estimates with a running QCD coupling, as well as CGC/saturation phenomenology, lead to λ = 0.2 ÷ 0.3. • In the vicinity of the saturation line G IP = φ 0 r 2 Q 2 s (1−γ cr ) (see Eq. 2.1). φ 0 denotes the value of the Pomeron Green function on the saturation line (at z = 0, r Q s = 1). The exact value of φ 0 cannot be determined without specifying the Pomeron-hadron interaction in more detail than we have. However, φ 0 ∝ᾱ 2 S so we expect φ 0 to be small.
• m 1 and m are pure phenomenological parameters, in our formulas we assumed that m m 1 . We make this assumption in order to simplify the formula.
• g is a pure phenomenological parameter which we assumed to be larger than G 3IP .

Cross sections and the values of the parameters
We determine the five parameters that our model depends on, by fitting to experimental data for the following set of observables: total, inelastic, and elastic cross sections, for single and double-diffractive production cross sections, and for the slope of the forward elastic differential cross section. We fit to the high energy data with W ≥ 0.546 TeV. We choose the minimal energy W = 0.546 TeV in our fit, as starting from this energy the CGC/saturation approach is able to describe the data on inclusive production in proton-proton collisions (see Ref. [100]). On the other hand the energy W = 0.2 TeV is too low, as at this energy saturation occurs in ion-ion and proton-ion collisions, but not in proton-proton collisions [101]. The quality of the fit can be seen from Fig. 5 and the values of parameters are presented in Table 1.
Our first observation is that the values of all parameters are in agreement with our expectations given in Sect. 3.4. Second, the overall fit has χ 2 /d.o. f. ≈ 2 of which 40% is due to our failure to reproduce the TOTEM value of B el at W = 7 TeV, and therefore, except for this point, we have a reasonably good description of the data. been taken from Ref. [109] for energies less than the LHC energy. At the LHC energy for total and elastic cross section we use TOTEM data [9,10] and for single-and double-diffraction cross sections are taken from Ref. [5]  We discuss some regularities in our fit, which could be useful for further and more profound understanding of the microscopic physics.
• We obtain a good description both of the values and of the energy dependence for total, inelastic, and elastic cross sections in a wide energy range: W = 0.546 ÷ 57 TeV. • At lower energies the values of B el are rather close to the experimental ones, but a glance at Fig. 5d shows that the energy behavior in our model is milder than that of the experimental one. The LHC value of B el is considerably higher than our prediction. A natural con-  jecture would be that this behavior is a direct consequences of long standing and unsolved problem in the CGC/saturation approach: i.e. the large impact parameter behavior of the BFKL Pomeron [104,105], however, we do not think that this is correct. Indeed, the problem noticed in Ref. [104,105] is one of the principal problems of the CGC/saturation approach, and has not yet been solved. Several theoretical models [106,107] show that the correct (exponential decreasing exp(−μb)) large impact parameter behavior of the scattering amplitude, does not influence the BFKL Pomeron, and does not even produce shrinkage of the diffraction peak (i.e. energy behavior of B el ). We can also see indications of the weak influence of the large b tail of the profile function in the data. Indeed, we were able to describe both σ tot and σ el which could be possible, only if one profile function describes the typical b behavior correctly.
Since there is no reason to expect that large b behavior could be crucial in the description of the available data, we assume that the CGC/saturation approach will not be changed by the (incorrect) large b behavior of the BFKL Pomeron, and that all non-perturbative corrections have to be included in Q s [87]: the only dimensional parameter of the theory. It should be stressed that this assumption leads to correct exponential fall off of the amplitude in the framework of our approach. We believe that the weak energy dependence of the elastic slope follows mostly from rather small shadowing corrections (see Fig. 6b), which we hope will be alleviated in our future planned two channel treatment. • Our single-diffractive production cross section has values within the experimental error, but our model predicts a lower cross section at W = 7 TeV, than the one given by ALICE [5]. The model's results are shown in Table 2 and support the idea that the single diffraction is saturated at high energy. 3 In our model the saturation, and oscillatory behavior, of the single-diffraction cross section stems from two sources: the simple one channel parametrization and the robust features of CGC/saturation approach. We plan to develop a more sophisticated (say two channel) model to disentangle these two effects, and we trust that this will also cure the present oscillatory behavior. • The striking feature of our model and, perhaps of the data, is that the double diffractive cross section increases with energy (see Table 2).

Partial amplitudes
We believe that the information contained in the impact parameter dependence of partial amplitude is instructive for understanding the nature of strong interactions at high energy. It is also useful for illustrating the strong and weak features of the model. Our present model, which is a single channel model has only one elastic amplitude, A el (s, b) which is shown in Fig. 6a.
One can consider our proton as a gray disk (A el (s, b = 0) < 1), even at energies as high as W = 57 TeV. This behavior stems mostly from the b-dependence of the Green function of the dressed Pomeron (see Fig. 6b).
The typical b increases with energy. Note that such an increase has been seen in the behavior of B el versus energy (see Fig. 5d). This behavior is due to strong saturation in dipole-dipole scattering, as the slope of the Pomeron trajectory for the BFKL Pomeron, is equal to zero.
In Fig. 7 dσ SD (s, b) /db 2 Eq. (3.9) and dσ DD (s, b) /db 2 (Eq. 3.11) are plotted. Note that the single-and doublediffraction production have quite different distributions in b. dσ DD (s, b) /db 2 , as is expected, has a peripheral form having a minimum at b = 0, the maximum and the width of the b-distribution of dσ DD (s, b) /db 2 grows considerably with energy. On the other hand, the peripheral nature of the single-diffractive production starts to appear only at high energies. The typical distribution has two maxima, at b = 0 and b ≈ 1 f m, both decrease with energy, while the width of the distribution slowly increases with energy. Such unexpected behavior stems mostly from the rather transparent dipole-dipole interaction, which is due to the values of the fitted parameters in our model.
The feedback from the impact parameter behavior of the scattering amplitude is that there is a need to increase the shadowing corrections in the dipole-dipole scattering. From our experience with the description of soft interactions at high energies, we know that one way of achieving this is by using a two channel model.

Comparison with other models on the market
In brief, this model is a one channel eikonal-type model, with a dressed Pomeron whose form has been adapted from the CGC/saturation approach. It differs from our previous model (see review [22]) which is a two-component model having three different partial amplitudes (see Fig. 8a). The striking feature of the two-component model is that two amplitudes become black at low energies. The resulting elastic amplitude in the two-component model is shown in Fig. 8b, and is similar to that of our present model. This is not surprising, as both models provide a good description of the data. However, the fact that in our present approach none of the ingredients, elastic amplitude and Green function of the dressed Pomeron, reach the black limit, looks surprising, especially so, since the recent model proposed by KMR (see Ref. [111]; Fig. 9), supports the fact that two of the partial amplitudes are black.
We wish to emphasize that our two models: the present one and the two-component model, lead to completely different predictions for single diffraction: in the first the cross section is saturated, while in the second it grows with energy.
Over the past few years a number of models have been constructed [18,111,112] based on Reggeon Field Theory  Im Re whose results for energies below that of the LHC are similar i.e. they adhere to the general trend of the experimental data, in that their results for σ tot , σ el , B el , σ SD and σ DD increase with increasing energy. The same applies to the Monte Carlo program MBR [113] (which is an "event generator" based on an enhanced PYTHIA8 simulation) and QGSJET-II [21].
Following the appearance of the preliminary results for single-and double-diffraction cross sections by the TOTEM Collaboration [114,115] and the CMS Collaboration [116] which suggests that the growth of σ SD and σ DD may be leveling off (or even decreasing) at W = 7 TeV. KMR [20] have modified their model by including energy dependent couplings, so as to be in accord with the TOTEM results. We would like to stress that the published results of the ALICE Collaboration [117] have the single-and double-diffractive cross sections still increasing at LHC energies.

Conclusions
In this paper we present a first attempt to develop a consistent approach based on the BFKL Pomeron and the CGC/saturation approach for soft interactions at high energy. We follow our general strategy for constructing models for strong interactions at high energy i.e. to maximize the theoretical ingredients, and to minimize the number of phenomenological parameters.
We construct an eikonal-type model whose opacity is determined by the exchange of the dressed BFKL Pomeron. The Green function of the Pomeron is calculated in the framework of the CGC/saturation approach. Having only five parameters we obtain a reasonable description of the experimental data at high energies (W ≥ 0.546 TeV). One of these five parameters λ, determines the energy dependence of the saturation scale, its value λ = 0.323 is a bit higher than the values that have been found from the description of the DIS and heavy ion scattering data, but it is close to them.
Using the value of λ from the fit we can estimate the value of the intercept of the BFKL Pomeron since λ = 4.88ᾱ S while BFKL = 2.8ᾱ S ≈ 0.2. From Eq. (2.11) we see that we can trust the MPSI approximation for Y ≤ 36, and therefore, the MPSI approximation provides the exact answer for the entire kinematic region of energies quoted.
In our model we find different behavior for the singleand double-diffraction cross sections at high energies. The single diffraction reaches a saturated value (about 10 mb) at high energies, while the double-diffraction cross section grows steadily. The reason for this is the different energy and impact parameter dependences, of the diagrams describing σ SD (Fig. 4a) and σ dd (Fig. 4b).
It turns out that in the model, all ingredients are far from being a black disc, in contradiction to our previous model. This illustrates how important it is to find a theoretical approach for soft processes.
We consider this paper as a first attempt to expand the CGC/saturation approach to describe soft processes at high energy. We plan to include more details of the CGC/saturation theory in our model and, in particular, to account for the running QCD coupling and to develop a two channel model to disentangle the two effects: the simple eikonal approach and the CGC/saturation features, which are included in our model. This paper provides an illustration of how the LHC data has stimulated our thinking.