Beyond DGLAP improved saturation model

We present a modification of the DGLAP improved saturation model with respect to the nonlinear correction (NLC). The GLR-MQ improved saturation model is considered by employing the parametrization of proton structure function due to the Laplace transforms method, which preserves its behavior success in the low and high Q2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q^{2}$$\end{document} regions. We show that the geometric scaling holds for the GLR-MQ improved model in a wide kinematic region rQs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$rQ_{s}$$\end{document}. These results are comparable with other models in a wide kinematic region rQs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$rQ_{s}$$\end{document}. The behavior of the dipole cross sections, with respect to the GLR-MQ improved saturation model, are comparable with the Color Glass Condensate (CGC) model. The model describes the dipole cross sections in the inclusive and diffractive processes. We also compare the nonlinear corrections to the impact-parameter dependent saturation (IP-Sat) model with the impact-parameter dependent color glass condensate (b-CGC) dipole model. Finally, we consider the linear and nonlinear corrections to the IP Non-Sat model. These results provide a benchmark for further investigation of QCD at small x in future experiments such as the Large Hadron Collider and Future Circular Collider projects.


I. Introduction
The color dipole picture (CDP) [1] has been introduced to study a wide variety of small x inclusive and diffractive processes at HERA.The dipole approach, at small values of Bjorken x, gives a clear interpretation of the highenergy interactions.This regime of QCD is characterized by high gluon densities because the proton structure is dominated by dense gluon systems [2][3][4] and predicts that the small x gluons in a hadron wavefunction should form a Color Glass Condensate [5].The gluon saturation effects are observable at very small x values and characterized by a hard saturation momentum Q s (x).The saturation scale is a border between dense and dilute gluonic systems as where xg(x, Q 2 ) is the gluon distribution function and πR 2 is the target area where R is the correlation radius between two interacting gluons.Indeed the parameter R controls the strength of the nonlinearity.The saturation scale rises with decreasing x and at small enough x, Q s ≫Λ QCD where Λ QCD is the QCD cut-off parameter at each heavy quark mass threshold (i.e., Λ n f QCD ).Since nonlinear dynamics are known to become sizable only at small-x, so the nonlinear contribution to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution [6] leads to an equation of the form * Electronic address: boroun@razi.ac.ir where α s ≡α s C A /π and the value of R is order of the proton radius (R ≃ 5 GeV −1 ), if the gluons are distributed through the whole of proton, or much smaller (R ≃ 2 GeV −1 ) if gluons are concentrated in hot spot within the proton.This was a vast subject initiated by Gribov, Levin, Ryskin, Mueller and Qiu (GLR-MQ) [7], as the second nonlinear term in ( 2) is responsible for gluon recombination.This term arises from perturbative QCD diagrams which couple four gluons to two gluons.So that two gluon ladders recombine into a single gluon ladder.It leads to saturation of the gluon density at low Q 2 with decreasing x.A closer examination of the small x scattering is resummation powers of α s ln(1/x) where leads to the k T -factorization form [8]. In the k T -factorization approach the large logarithms ln(1/x) are relevant for the unintegrated gluon density in a nonlinear equation.Solution of this equation develops a saturation scale where tame the gluon density behavior at low values of x and this is an intrinsic characteristic of a dense gluon system [9].
The main goal of this paper is to consider the nonlinear corrections to the DGLAP improved saturation model.In fact, the DGLAP improved saturation model will be modified to the GLR-MQ improved saturation model.This is based on the nonlinear evolution of the gluon density at small values of x.These results will be compared with the nonlinear saturation dynamics which is explicitly incorporated into the CGC model.One of the well-known impact-parameter dependent saturation models is the IP-Sat model [10,11].This is a simple dipole model that incorporates the physics of saturation and all known properties of the gluon saturation.In this case the saturation boundary is approached via the DGLAP evolution, that is, by the eikonalization of the gluon distribution, which effectively represents higher twist contributions.The b-CGC and the IP-Sat models are easily generalized from DIS off protons to DIS off nuclei [12].This paper is organized as follows.In Sec.II, we introduce the color dipole model for calculating the dipole cross sections in the GBW, the DGLAP improved saturation, the b-CGC dipole , the IP-Sat models and also the exclusive diffractive processes.In Sec.III, we present the GLR-MQ improved saturation model to consider the color dipole cross section at low values of x.Then in Sec.IV, we present a detailed numerical analysis and our main results.We summarize our main results in Sec.V.

II. Dipole cross section
Dipole representation provides a convenient description of DIS at small x.There, the scattering between the virtual photon γ * and the proton is seen as the color dipole where the transverse dipole size r and the longitudinal momentum fraction z with respect to the photon momentum are defined.The amplitude for the complete process is simply the production of these subprocess amplitudes, as the DIS cross section is factorized into a lightcone wave function and a dipole cross section.Using the optical theorem, this leads to the following expression for the γ * p cross-sections where subscripts L and T referring to the transverse and longitudinal polarization state of the exchanged boson.
Here Ψ L,T are the appropriate spin averaged light-cone wave functions of the photon and σ dip (x, r) is the dipole cross-section which related to the imaginary part of the (qq)p forward scattering amplitude.The variable z, with 0 ≤ z ≤ 1, characterizes the distribution of the momenta between quark and antiquark.The square of the photon wave function describes the probability for the occurrence of a (qq) fluctuation of transverse size with respect to the photon polarization [1,2].The dipole hadron cross section σ dip contains all information about the target and the strong interaction physics with where b is a particular impact parameter (IP) as and S(b) is the S-matrix element of the elastic scattering.The cross section at a given impact parameter b is proportional to the dipole area, the strong coupling, the number of gluons in the cloud and the shape function by the following form [10] where the hard scale is assumed to have the form and the parameters C and µ 2 0 are obtained from the fit to the DIS data [1].For multi Pomeron exchange, the eikonalised dipole scattering amplitude of Eq.( 6) can be expanded as where dσ dip /d 2 b = 2N (x, r, b) and the n-th term in the expansion corresponds to n-Pomeron exchange [10].Eq.( 6) is known as the Glauber-Mueller dipole cross section [13] and can also be obtained within the McLerran-Venugopalan model [14].The exponential form of the function T (b) is determined from the fit to the data as where the parameter B G was found [10] to be 4.25 GeV −2 .
In the original Golec-Biernat-Wüsthoff (GBW) model [1], the dipole cross section was proposed to have the eikonallike form where Q s (x) plays the role of the saturation momentum, parametrized as Q 2 s (x) = Q 2 0 (x/x 0 ) −λ .Parameters Q 0 and x 0 set dimension and absolute value of the saturation scale and exponent λ governs x behavior of Q 2 s .This model was updated in [2,15] to improve the large Q 2 description of the proton structure function by a modification of the small r behavior of the dipole cross section to include the DGLAP evolved gluon distribution.Since the energy dependence in large Q 2 region is mainly due to the behavior of the dipole cross section at small dipole size r, therefore authors in Refs.[2,15] investigated the DGLAP evolution for small dipoles.Bartels-Golec-Bienat-Kowalski (BGBK) improved the dipole cross section by adding the collinear DGLAP effects.Indeed the BGBK model is the implementation of QCD evolution in the dipole cross section which depends on the gluon distribution.The following modification of the DGLAP improved saturation model [1] proposed for the dipole cross section as Indeed BGBK model is successful in describing dipole cross section at large values of r as the two models (GBW and BGBK) overlap in this region but they differ in the small r region where the running of the gluon distribution starts to play a significant role.Indeed the DGLAP improved model of σ dip significantly improves agreement at large values of Q 2 without affecting the physics of saturation responsible for transition to small Q 2 .As expected, geometrical scaling is true for the DGLAP improved model curve for the scaling variable rQ s ≥1 and for the GBW model curve for the whole region [1].
The saturated version of the dipole model may in principle be derived from the Color Glass Condensate effective theory for QCD according to Eq.( 6) where at small r this expression (i.e., Eq.( 6)) becomes Eq.( 6) is referred to as the IP-Sat model, while Eq. ( 11) is referred to as the IP Non-Sat model.The BGBK and CGC models considered only the dipole cross section integrated over the impact parameter b [16].The BGBK model was modified to include the impact parameter dependence as denoted by the IP-Sat model and the CGC model was also modified to include the impact parameter dependence as denoted by the b-CGC model.The dipole cross section can be calculated in the CGC approach from the relation where σ 0 = 2πR 2 p and where Y = ln(1/x) and k = χ ′′ (γ s )/χ ′ (γ s ) where χ is the LO BFKL [17] characteristic function.The scattering amplitude N (x, r) can vary between zero and one, where N = 1 is the unitarity limit.To introduce the impact parameter dependence into the CGC model, the b-CGC model for the dipole cross section is defined by the following form [16] where the impact parameter dependence of the saturation scale Q s was introduced by where the parameter B CGC , instead of σ 0 in the CGC dipole model, is a free parameter and is determined by other reactions, namely the t distribution of the exclusive diffractive processes at HERA.The parameters were fixed by a combination of theoretical constraints [11] and a fit to DIS data.Another one of the main advantages of dipole models is the description of the diffractive process [2,18].The cross section for the diffractive qq production reads [8] dσ where t = ∆ 2 , and ∆ is the four-momentum transferred into the diffractive system from the proton.In Eq.( 16), the generalised optical theorem is applied in the framework of the dipole picture.At small values of the diffractive mass M 2 ∼ Q 2 the elastic scattering of the qq pair dominates, while at larger values of the mass M 2 ≫Q 2 , the qqg contribution dominates (due to gluon production in the final diffractive state).The treatment of the qqg component goes beyond the saturation model since this is not present in the inclusive analysis [2,18].This component was computed in the two gluon exchange approximation with an additional assumption of strong ordering of transverse momenta of the qq pair and the gluon.
In the transverse coordinate representation, the qqg system is treated as a color octet dipole 88 where the coupling of two t-channel gluons is relative by a weight factor where N C is the number of colors.Thus, the color dipole cross section for exchange of a two gluon system for octet dipole reads [2,18] In the next section, we consider the color dipole cross sections due to the behavior of the linear and nonlinear gluon density and compare with the other models.
The linear gluon densities are obtained with respect to the Laplace transform technique by employing the parametrization of proton structure function, then applied the GLR-MQ evolution equation for the nonlinear gluon densities.Some approximated analytical solutions in the color dipole model have been reported in recent years [19,20] with considerable phenomenological success due to a parametrization of the deep inelastic structure function for electromagnetic scattering with protons.

III. GLR-MQ improved saturation model
We will present an approach to the description of the color dipole cross section at small x, alternative to that based on the DGLAP improved saturation model.From a more theoretical viewpoint it is known that in the low x, low Q 2 region gluon recombination effects are not negligible and reduce the growth of the gluon parton distribution function.The GLR-MQ equation for the gluon density, where an extra non-linear term, quadratic in the gluon density, was added to the linear DGLAP evolution equation by the following form where χ = x x0 and x 0 is the boundary condition that the gluon distribution (i.e., G(x, µ 2 ) = xg(x, µ 2 )) joints smoothly onto the linear region.We note that at x≥x 0 (= 10 −2 ) the non-linear corrections are negligible.The nonlinear shadowing term, ∝ −[g] 2 , arises from perturbative QCD diagrams.In this regime the gluons in the proton form a dense system with mutual interaction and recombination which also leads to the saturation of the total cross section.Other early works on this topic can be found in [21,22].In what follows, the hard scale is assumed to have the form µ 2 = C/r 2 + µ 2 0 as for light quarks the gluon distribution is evaluated at x = x BJ = µ 2 /(µ 2 + W 2 ) and for the charm quark the gluon structure function is evaluated at where m c is the charm quark mass and W refers to the photon-proton center-of-mass energy.The non-linear equation (i.e., Eq.( 18)) shows that the strong rise that is corresponding to the linear QCD evolution equations at small-x and Q 2 can be tamed by screening effects.The first iteration of Eq.18 reads where the nonlinear correction to the gluon distribution function (i.e., G NLC (x, µ 2 ) ) is obtained by the following form Here G(x, µ 2 ) and G(x, µ 2 0 ) are the linear gluon distributions, and obtained from the parametrization F 2 using the Laplace transform techniques [23,24], at µ 2 and µ 2 0 scales respectively.At the initial scale µ 2 0 , the low x behavior of the non-linear gluon distribution is assumed to be [25] Therefore the non-linear correction to the gluon distribution at µ 2 scale for x < x 0 reads The gluon distribution due to the non-linear corrections can be analytically solved at small x with respect to the linear gluon distribution behavior.The linear gluon distributions (i.e., G(x, µ 2 ) and G(x, µ 2 0 )) in Eq.( 23) are defined with respect to the most parametrization suggested in Refs.[23] and [24].The authors in Ref. [23] have an expression for the asymptotic part of F 2 (no-valence) as for x≤0.09.In Ref. [24], the authors obtained two quadratic expressions in ln 2 (1/x) using second-order linear differential equation as well as Laplace transforms for the leading-order (LO) gluon distribution function, respectively.In the first method in Ref. [24], the LO DGLAP equation for the evolution of the proton structure function F 2 (x, Q 2 ) is rearranged into an inhomogeneous second-order differential equation by the following form Eq.( 25), with the new variable υ = ln(1/x) becomes a linear 2 nd order inhomogeneous equation, as and the definition G(υ, Q 2 ) = G(e −υ , Q 2 ).In Ref. [24], the authors have found the parametrization of the gluon distribution G 4 (υ, Q 2 ) which is calculated as a second degree polynomial in υ whose coefficients are quadratic polynomials in ln(Q 2 ) for x≤0.09 as Therefore where the gluon distribution in x-space reads as a simple quadratic polynomial in ln(1/x) with quadratic polynomial coefficients in ln(Q 2 ) by the following form In the second method, the authors [24] have suggested a new parametrization based on Laplace transforms.The DGLAP evolution is written as follows where w = ln(1/z) and The function h(υ) in Eq.( 29) is where P gq is the gluon-quark splitting function.The function F 2 (x, Q 2 ) in Eq.( 30) is sum of the proton structure function F 2 -dependent terms in the DGLAP evolution equation by By making a Laplace transform in υ, we can factor Eq.( 29), since the Laplace transform of a convolution is the product of the Laplace transform of the factors, so that Solving Eq.( 29) for g in s-space, we have Thus, inverting the Laplace transform of the factors, then the gluon distribution is defined by Therefore, the gluon distribution in x-space reads for 0 < x≤0.06.The standard representation for QCD coupling in LO approximation is defined by where β 0 is the one loop correction to the QCD βfunction and t = ln Q 2 Λ 2 , Λ is the QCD cut-off parameter with α s (M 2 z ) = 0.118.The ln 2 (1/x) behavior of the DIS proton structure function (i.e., Eq.( 24)) at small values of x is compatible with saturation of the Froissart bound at each value of Q 2 .The authors, in Ref. [23], have shown that this behavior may be the signal for the saturation or gluon recombination processes at high parton densities.The gluon distribution in Eqs.(26)(27)(28)(29)(30)(31)(32)(33)(34)(35)(36), according to the results in Ref. [24], is determined from the DGLAP evolution equation for the proton structure function.Thus in Eq.( 20), the nonlinear corrections to the gluon behavior at low x and Q 2 values are considered, where it is compatible with ln 2 (1/x) behavior of parton densities at very small x in the QCD evolution framework.Now, we can estimate the non-linear corrections to the gluon distribution (i.e., Eq.( 23)) due to the linear gluon distributions (i.e., Eq.( 28) and (36)) for small x and we will use the non-linear corrections to the dipole cross sections, and in the next section, the accuracy of the results will be discussed in comparison with the CGC model.

IV. Results
The linear and nonlinear methods are presented based on the solutions of the DGLAP and GLR-MQ evolution equations at the leading-order accuracy in perturbative QCD, respectively.The dipole cross-sections (Eqs.6, 10, 11 and 17) require the gluon density G(x, µ 2 ) for all scales µ 2 .These gluon distributions [24] are obtained directly in terms of the parameterization of the structure function F 2 (x, µ 2 ) and its derivative.The resulting linear and nonlinear gluon distribution functions for various dipole sizes for x≤10 −2 are shown in Fig. 1.The dipole size determines the evolution scale µ 2 .In this figure, we plot the r dependence of the nonlinear corrections to the gluon distribution for R = 2 GeV −1 .Nonlinear corrections play an important role on gluon distribution as x and µ 2 decrease.A depletion occurs at x < 10 −2 where these results show that the nonlinear behavior of the gluon distribution function is tamed.This taming behavior of nonlinear gluon distribution function toward low x at low µ 2 values becomes significant when considering the color dipole cross section at the hot spot point (i.e., R = 2 GeV −1 ).We have calculated the linear and nonlinear corrections to the ratio σ dip /σ 0 in a wide range of the dipole size at the LO approximation.Results of  28) and ( 36)) and nonlinear (Eq.( 23)) gluon distribution functions for various dipole sizes.
calculations and comparison with the GBW [1] and CGC [5] models for x = 10 −4 are presented in Fig. 2. The linear corrections to the ratio of color dipole cross sections at LO approximation are comparable with the GBW model at low and high r values.The nonlinear corrections to the ratio of color dipole cross sections are comparable with the GBW model for r 10 −2 and r≥10 0 and also are comparable with the CGC model for 10 −2 r≤10 0 .Indeed the nonlinear corrections tame the behavior of the dipole cross section at r 10 −1 .The effective parameters in the GBW model have been extracted from a fit of the HERA data as, λ = 0.288, x 0 = 3.04×10 −4 , C = 0.38 and µ 2 0 = 1.73 [1].Parameters of the CGC dipole model fixed at the LO BFKL according to the original CGC fit [5] with respect to the values γ s = 0.63, k = 9.9, N 0 = 0.7, λ = 0.177 and x 0 = 2.70×10 −7 [16].The dipole cross sections are evaluated according to the four active flavors, which take into account charm quark mass.The quark mass, in the CGC model, was taken to be 1.4 GeV although in our calculations it is 1.29 GeV [26].
An important property of the saturation formalism is the geometric scaling phenomenon, which means that the scattering amplitude and corresponding cross sections can scale on the dimensionless scale rQ s .A particular interests present the linear and nonlinear ratio σ dip /σ 0 defined by the scaling variable rQ s .In Fig. 2 (right hand), we observe that the nonlinear corrections to the ratio σ dip (rQ s (x))/σ 0 merge into the GBW curve for rQ s 10 −1 .The results of the GLR-MQ improved saturation model due to the parametrization of the proton structure function have become a function of a single variable, rQ s , for almost all values of r at LO approximation.
a qq pair in the color singlet state as the diffractive γ * p→qqp ′ cross section is proportional to σ 2 (x, r) by Eq.( 16).We have calculated the ratio σ 2 dip /σ 2 0 for the diffractive qq production into r and rQ s respectively and compared the ratio with the GBW and CGC models in Fig. 3.The linear corrections to the ratio σ 2 dip /σ 2 0 are comparable with the GBW model in a wide rang of r although the nonlinear corrections are comparable with the CGC model for r < 1 and with the GBW model for r≥1.In Fig. 3 (left hand), the geometrical scaling of the nonlinear corrections to the ratio σ 2 dip /σ 2 0 is visible in comparison with the linear curve.The nonlinear curve merges into one solid line in the right plot where the dipole cross section is plotted as a function of the scaling variable rQ s .This is a reflection of geometric scaling in the nonlinear corrections in comparison with the GBW model for the diffractive qq production.
In addition to the contributions of the qq states, it is
important to include the contributions of the qqg final states of the diffractive processes in the nonlinear corrections to the ratio σ 2 dip /σ 2 0 .In Fig. 4 the ratio of the dipole cross sections are determined by the qqg component and are compared with the GBW and CGC models.The linear and nonlinear ratio of the dipole cross sections, in comparison with the results in Fig. 3, are deviated from the GBW and CGC models respectively.The reason for this deviation is because the qqg component, interacting with the proton with the same dipole cross section as the qq system, goes beyond the saturation model [18].Indeed, in Fig. 4, the linear and nonlinear cross sections are modified due to the weighted factor C A /C F although this component is not present in the inclusive analysis.Now we consider the nonlinear corrections to the the qq differential cross section dσ dip /d 2 b.In Fig. 5, the linear and nonlinear corrections to the impact parameter dependent dipole cross section due to the GLR-MQ equation are considered and compared with the b-CGC model for x = 10 −4 .In this figure (i.e., Fig. 5) the linear corrections to the IP-Sat (i.e., b-Sat) model are comparable with the b-CGC model in a wide range of the impact pa- The linear and nonlinear corrections to the impact parameter dependent dipole cross section versus the impact parameter b (Eq.( 6)) compared with the b-CGC model (Eq.( 14)) for the dipole sizes r = 0.1, 1 and 2 fm at x = 10 Indeed, the GLR-MQ improved saturation model illustrates unitarity with an increase of r as b decreases.At small r, the IP-Sat model (Eq.6) becomes the IP Non-Sat model (Eq.11)where the interaction between the dipole and the hadron is described by the exchange of one gluon.The linear and nonlinear behavior of dσ dip /d 2 b in the IP Non-sat model are considered in Fig. 6.In this model, the behavior of the dσ dip /d 2 b is directly dependent on the gluon distribution function.Saturation effects are not visible in this model as b decreases.However this behavior tamed due to the nonlinear corrections to the gluon density.For small dipole sizes the distributions are almost similar but they differ significantly as r becomes large.A comparison of the resulting dσ dip /d 2 b according to the linear as well as nonlinear behavior of the dipole cross section for b = 0 at x = 10 −4 presented in Fig. 7.The resulting dipole cross-sections in linear and nonlinear corrections are shown in Fig. 7.We observe that, in this figure, the nonlinear corrections suppress the behavior of the large dipoles in the IP Non-Sat model.Indeed, this behavior tamed at large r for b = 0 where with the increase r, µ decreases to the value of µ 0 .We also note that adding nonlinear corrections to the IP-Sat model decreases the dipole cross-section for 0.2 fm < r < 1 fm at b = 0.The linear and nonlinear corrections to the dipole cross section to the IP-Sat model reach the saturation plateau at r > 1 fm.It is interesting to increase the impact parameter from for the IP-Sat (Eq.( 6)) as well as IP Non-Sat (Eq.( 11)) versus r at x = 10 −4 for b = 0 and b = 1 fm, respectively.b = 0 to fm for the nonlinear behavior of the dipole cross sections to the IP Non-Sat in Fig. 8.The proton dipole cross section at different impact parameters with and without nonlinear corrections are shown in Fig. 8 for the IP-Sat as well as IP Non-Sat versus r at x = 10 −4 .The IP-Sat and IP Non-Sat dipole cross sections are very similar in the range 0≤r≤4 fm for b = 1 fm.Consequently, for large impact parameter sizes the distributions are almost similar but they differ significantly as b becomes small due to the nonlinear corrections.Indeed, the nonlinear corrections become stronger at larger impact parameters for the IP Non-Sat model [3].In Fig. 9 we show the IP Non-Sat to IP-sat cross-section ratios as a function of r for x = 10 −4 .We depict the ratio as a function of r for b = 0 and 1fm.Note that the ratio increases much faster as a function of r for b = 0 than for b = 1fm.We further note that at larger r, the ratio remains near almost unity for b = 1fm.The large difference between IP Non-Sat and IP-Sat comes from the decreases in the impact parameter values.

V. Conclusions
In this paper we proposed a modification of the saturation model which takes into account the GLR-MQ evolution of the gluon distribution.We have presented a certain theoretical model at LO approximation to describe the color dipole cross sections based on the Laplace transforms method at small values of x (the Bjorken variable x is fixed to be x = 10 −4 ).We have used a nonlinear correction to the dipole cross sections from a parametrization of the proton structure function with a rescaled variable m c .The nonlinear corrections to the dipole cross sections in the description of inclusive and diffractive DIS at small x, according to the saturation scale and geometric scaling, are consistent with analytical saturation models in a wide range of r and rQ s , respectively.We find that the ratio σ dip /σ 0 due to the DGLAP improved saturation model is consistent with the GBW saturation model, although the nonlinear corrections to this ratio with respect to the GLR-MQ improved saturation model is consistent with the CGC saturation model especially in the range 0.05 < r≤1.In the simplest case of the qq system for the ratio σ 2 dip /σ 2 0 in the diffractive processes, the linear and nonlinear corrections show good agreement with the GBW and CGC models in a wide range of r and rQ s .The linear and nonlinear corrections to the ratio σ 2 dip /σ 2 0 in the diffractive processes due to the component qqg deviates from the GBW and CGC models, because the qqg system goes beyond the saturation models.We developed nonlinear corrections to the impact parameter dependent dipole cross sections, dσ dip /d 2 b.The nonlinear corrections to the IP-Sat model are comparable with the b-CGC model in a wide range of the impact parameter b and the dipole size r.The linear and nonlinear corrections considered in the IP-Sat and IP Non-Sat models for the impact parameters b = 0 and b = 1 fm in the range 0≤r≤4.The behavior of the nonlinear corrections to the IP Non-Sat model tamed in a wide range of r.This behavior for the IP-Sat model shows that the dipole cross section saturated early for b = 0 in comparison with b = 1 fm for r > 1 fm.The nonlinear corrections to the IP-Sat and IP Non-Sat models show that those behaviors are comparable in the range 0≤r≤4 for the impact parameter b = 1 fm.In conclusion, by considering the statistical errors due to the effective parameters, the nonlinear corrections to the dipole cross sections give a reasonable data description in comparison with the other models.Indeed, the GLR-MQ improved saturation model tames the DGLAP improved model behavior when the results compared to models described based on the recombination of gluons at low x.
FIG.5:The linear and nonlinear corrections to the impact parameter dependent dipole cross section versus the impact parameter b (Eq.(6)) compared with the b-CGC model (Eq.(14)) for the dipole sizes r = 0.1, 1 and 2 fm at x = 10 −4 .

− 4 .
FIG.5:The linear and nonlinear corrections to the impact parameter dependent dipole cross section versus the impact parameter b (Eq.(6)) compared with the b-CGC model (Eq.(14)) for the dipole sizes r = 0.1, 1 and 2 fm at x = 10 −4 .