the scattering experiments on deep inelastic electron - proton ( dis ) at hera have provided measurements of the inclusive structure function @xmath0 as well as of the @xmath1 and @xmath2 in very small bjorken variable @xmath3 values ( @xmath4 ) . in these collisions
the proton target is analysed through a hard probe with virtuality @xmath5 .
the momentum fraction is @xmath6 , where @xmath7 and @xmath8 are the four - momenta of the incoming proton and of the virtual photon probe . in the kinematical region of small @xmath3
, the gluon is the leading parton driving the behavior of the deep inelastic observables .
the standard qcd evolution @xcite gives a powerlike growth for the gluon distribution and related quantities , and this feature leads , in principle , to the unitarity violation at asymptotic energies , requiring a control of the gluon distribution in high energies . in the partonic language , in the infinite momentum frame , the small momentum fraction region corresponds to the high parton density domain , which is connected with the black disk limit of the proton target and with the parton recombination phenomenon .
these issues can be addressed through a non - linear dynamics beyond the usual dglap formalism ( for a review , see @xcite ) .
the complete knowledge about the non - linear dynamical regime plays an important role in the theoretical description of the reactions in the forthcoming experiments rhic and lhc .    in the breit frame
, the qcd factorization theorem allows to calculate the scattering processes through the convolution of the partonic subprocess with the parton distribution functions ( pdf s ) .
there , the degrees of freedom of the theory are the quasi - free partons ( quarks and gluons ) . on the other hand ,
more recently the theoretical description of the small @xmath3 physics has been widely analyzed in the target rest frame , which is a powerful tool concerning an unified picture for both the inclusive and the diffractive scatterings , including vector meson production @xcite .
now , the degrees of freedom are the color dipoles , which are the most simple configurations considering the virtual photon fock states expansion .
its main appeal is a quite simplified picture for the different mentioned processes , based on general properties of quantum mechanics .
the description of dis in the color dipole picture is quite intuitive , allowing a simple representation instead of the involved one from the breit ( infinite momentum ) frame , and such framework was proposed by gribov many years ago @xcite .
considering small values of the bjorken variable @xmath3 , the virtual photon fluctuates into a @xmath9 pair ( dipole ) with fixed transverse separation @xmath10 at large distances upstream of the target and interacts in a short time with the proton .
more complicated configurations should be considered for larger transverse size systems , for instance the photon fock state @xmath11 .
an immediate consequence of the lifetime of the pair ( @xmath12 ) to be bigger than the interaction s one is the factorization between the photon wavefunction and the cross section dipole - proton in the @xmath13 total cross section .
the wavefunctions are perturbatively calculable , namely through qed for the @xmath9 configuration @xcite and from qcd for @xmath14 @xcite .
the effective dipole cross section should be modeled and encloses both perturbative and non - perturbative content .
however , since the interaction strength relies only on the configuration of the interacting system the dipole cross section turns out to be universal and may be employed in a wide variety of small @xmath3 processes @xcite .
currently , there are several models for the dipole cross section based either in pure phenomenological parametrizations or in a more theoretical ground ( for a review , see ref .
the main feature in those models is the description of the energy dependence of the interaction , namely taking into account the interplay between hard and soft domains .
the dipole cross section should be consistent with the sharp growth on energy at small transverse separation @xmath10 ( large gluon density ) and a softer behavior for larger @xmath10 ( regge - like phenomenology ) .
an additional ingredient is the expectation for saturation effects in high energies as a consequence of the unitarity requirements .
indeed , the growth of the gluon distribution should be tamed at very small @xmath3 and it has been found that the corrections are important already in the present hera kinematics @xcite . moreover , such effects are associated with higher twist contributions concerning the standard linear evolution equation , i.e. the dglap formalism @xcite .    here , we take into account a formalism providing the unitarity corrections to the deep inelastic scattering at small @xmath3 , namely the glauber - mueller approach in qcd .
it was introduced by a. mueller @xcite , who developed the glauber formalism to study saturation effects in the quark and gluon distributions in the nucleus considering the heavy onium scattering .
afterwards , the authors of ref .
@xcite extended that approach that has as a limit the glr results @xcite .
it is obtained an evolution equation taking into account the unitarity corrections ( perturbative shadowing ) , generating a non - linear dynamics which is related with higher twist contributions .
its main characteristic is to provide a theoretical framework for the saturation effects , relying on the multiscattering of the pqcd pomeron .
the latter is represented through the usual gluonic ladder in the double logarithmic approximation .    summarizing the glauber - mueller approach , the gluonic content of the nucleus or nucleons
is obtained in the following way : in the rest frame , a virtual probe ( gluon ) decays into a gluon pair interacting with the nucleon inside the nucleus .
the multiple scatterings of the pair give rise to the unitarization of the corresponding cross section .
the calculations are performed in the double logarithmic approximation ( dla ) , corresponding to the condition @xmath15 , where @xmath16 and @xmath17 are the qcd coupling constant and the gluon anomalous dimension , respectively . in this approximation ,
the transverse separation of the pair remains fixed allowing an eikonal description for the interaction gluon pair - nucleon through incoherent multiscatterings @xcite .
the cross section for the interaction can be expressed in terms of the nucleon gluon distribution @xmath18 and of the transverse separation @xmath10 of the gluon pair .
the procedure for an initial state quark - antiquark pair is similar , up to proper color coefficients .
such a formulation has produced comprehensive phenomenological applications : the inclusive structure function @xmath19 @xcite , the longitudinal @xmath1 and charmed @xmath2 ones @xcite , the logarithmic slope @xmath20 @xcite and other related quantities have been calculated in the breit system .
the respective nuclear case has been investigated in refs . @xcite and
the equivalence with other high density qcd approaches has been reported in refs .
moreover , the asymptotic limit of the inclusive structure function and its logarithmic slope are estimated in @xcite . regarding the rest frame ,
glauber - mueller has been also used to estimate the saturation effects for dis and diffractive dissociation @xcite as well as it is considered as initial condition for a high energy evolution equation ( for a review , see @xcite ) .
therefore , the glauber - mueller approach gives a good framework for the unitarity effects ( saturation ) in the nucleon and nuclear sectors , providing the dynamics of the observables in a quantitative level .    in this work
we make use of the parton saturation formalism to study the description of the observables driven by the gluonic content of the proton in the color dipole picture .
the inclusive structure function @xmath19 is calculated properly , disregarding approximations commonly considered in previous calculations @xcite .
the structure functions @xmath1 and @xmath2 are presented for the first time using the glauber - mueller approach and the rest frame in comparison with the experimental data .
the saturation effects are included in the effective dipole cross section corresponding to the small size dipole contributions .
the large dipole sizes are taken into account through an anstz for the non - perturbative region . here , we choose to freeze the gluon distribution at large distances @xcite .
the dipole framework provides a clear identification in the transverse distance @xmath10 range where the perturbative and non - perturbative sectors ( soft domain ) contribute to the calculated quantities .
large dipoles correspond to soft domain and the small ones are connected with the hard sector .
as we will see , the photon wavefunctions play the role of a weight function selecting small dipole sizes , with a non - negligible contribution coming from the large dipole sizes .
however , when considering the production of heavy quark pairs ( charm and bottom ) , their masses largely diminish the soft contribution .
these issues are addressed throughout this paper .
this work is organized as follows . in the next section we shortly review the deep inelastic scattering in the rest frame , introducing the main formulae for the further analysis . in sec .
( 3 ) one addresses the effective dipole cross section considering saturation effects ( unitarity corrections ) encoded in the glauber - mueller approach , pointing out its main properties and discussing the large transverse separation contribution .
( 4 ) is devoted to calculate the theoretical estimates for the hera observables , focusing on those ones dominated by the gluon content . in the last section
we draw our conclusions and comments .
in the reference frame where the target ( proton ) has infinite momentum the usual description of the dynamics evolution is the following : the emitted partons from the target remain quasi - free for enough time in such a way that the virtual probe ( photon ) detects them as real particles ( asymptotic freedom ) . on the other hand , in the small @xmath3 region
it is suitable to describe the evolution in a system where the target is at rest : in this situation the evolution is related to the partonic fluctuations of the probe and their interactions with the nucleon .
the rest frame physical picture is advantageous since the lifetimes of the photon fluctuation and of the interaction process are well defined @xcite .
the simplest case is the quark - antiquark state ( color dipole ) , which is the leading configuration for small transverse size systems .
its life time can be estimated by the uncertainty principle through the energy fluctuation associated with the emerging pair . the well known coherence length is expressed as @xmath21 , where @xmath3 is the bjorken variable and @xmath22 the proton mass .
for instance , in deep inelastic at hera kinematics reaching at @xmath23 , the coherence length is about @xmath24 fm , which is a distance larger than the radius of any atomic nuclei .
an immediate consequence of such a picture is the factorization between virtual photon wavefunctions and the interaction cross section on the corresponding amplitude in the impact parameter space representation .
a striking consequence of the formulation above is that the photoabsortion cross section can be derived from the expectation value of the interaction cross section for the multiparticle fock states of the virtual photon weighted by the light - cone wave functions of these states @xcite .
the scattering matrix has an exact diagonalization in the @xmath25 fock states representation , where the partial - wave amplitudes subjected to the @xmath26-channel unitarization are identified properly .
then , the photoabsortion cross section can be cast into the quantum mechanical factorized form , @xmath27    the formulation above is valid even beyond perturbation theory , since it is determined from the space - time structure of the process .
the @xmath28 are the photon wavefunctions ( for transverse , @xmath29 , and longitudinal , @xmath30 , polarizations ) describing the pair configuration , where @xmath31 and @xmath32 are the fraction of the photon s light - cone momentum carried by the quark and antiquark of the pair , respectively .
the transverse separation of the pair is @xmath33 .
the precise normalizations of the wavefunctions can be determined through the fock expansion @xmath34 , where @xmath35 is the @xmath36-wavefunction renormalization constant , @xmath37 denotes the bare state and @xmath38 is the coefficient determining the probability of the @xmath9 pair fluctuation in the photon @xcite . considering completely normalized states , then @xmath39 and the normalization for transverse photons
is obtained from @xmath40 .
the remaining normalizations , i.e. , longitudinal component and cross section , are consequently fixed @xcite .
the explicit expressions are well known , @xmath41\:\varepsilon^2\,k_1 ^ 2\,(\varepsilon\ , r )   \,+\ , m_q^2\ : k_0 ^ 2\,(\varepsilon\ , r)\,\big\ } \label{eq9}\end{aligned}\ ] ] and @xmath42    clarifying the notation , the auxiliary variable @xmath43 , with @xmath44 the light quark mass , and @xmath45 and @xmath46 are the mc donald functions of rank zero and one , respectively .    the quantity @xmath47 is interpreted as the cross section of the scattering of the effective dipole with fixed transverse separation @xmath48 @xcite .
it is directly dependent on the unintegrated gluon distribution , @xmath49    the unintegrated gluon distribution , @xmath50 , vanishes at the gluon transverse momentum @xmath51 ( and similarly for the factor @xmath52 ) due to gauge invariance , and therefore the dipole cross section has to be infrared finite .
the most important feature of the dipole cross section is its universal character , namely it depends only on the transverse separation @xmath48 of the color dipole .
the dependence on the external probe , i.e. the photon virtuality , is included in the wavefunctions .
the main technical difficulty in eq . ( [ eq11 ] ) is to model the unintegrated gluon distribution function in a suitable way , mainly in the small transverse momentum ( @xmath53 ) region ( infrared sector ) . in particular , to obtain these distributions one should solve numerically an evolution equation , which makes the procedure cumbersome for practical use ( instead , for a prompt parametrization see @xcite ) . in general
, this is avoided introducing an anstz for the effective dipole cross section and analyzing the process in the impact parameter space .
the main feature of the current models in the literature is to interpolate the physical regions of small transverse separations ( qcd - parton improved model picture ) and the large ones ( regge - soft picture ) .
below , we quote two of them , which have connections with the work performed here .
the phenomenological saturation model of golec - biernat and wsthoff ( hereafter gbw ) gives a good description of dis data @xcite .
the corresponding cross section interpolates between color transparency , i.e. @xmath54 at small @xmath10 , and constant cross section @xmath55 at large @xmath10 ( confinement ) .
such a procedure ensures @xmath56-saturation , while parton saturation at low @xmath3 is obtained with an eikonal - inspired shape for the dipole cross section , @xmath57 where the parameters @xmath58 mb , @xmath59 and @xmath60 are fitted to the hera dis inclusive data with   @xmath61 , whereas @xmath62 gev@xmath63 sets the dimension .
the @xmath3-dependent saturation radius @xmath64 scales the pair separation @xmath10 in the cross section and is associated with the mean separation between partons in the nucleon .
this approach was used to describe diffractive dissociation in a parameter - free way , considering also the required @xmath14 configuration . some criticism to the gbw model , mainly concerning a better knowledge of the parton ( gluon ) distribution at low @xmath10 and its saturation at small @xmath3
are postponed to the next section .
the most of them can also to be found in refs .
@xcite .
our work in this paper is closer to the mcdermott et al .
one @xcite ( hereafter mcdfgs ) , where perturbative qcd relates the dipole cross section to the leading logarithmic gluon distribution function in the proton , at lla(@xmath56 ) accuracy , @xmath65 where the identification @xmath66 is allowed at leading - log level .
a study of the dipole cross section is carried for all transverse separations @xmath10 in @xcite .
the small @xmath10 region is described by the usual gluon pdfs , evoluting with the scale @xmath67 , while for large separations ( @xmath68 fm ) the dipole cross section is driven by the pion - proton contribution with the typical soft energy behavior from the hadronic sector . moreover
, the taming of the parton density is implemented by hand starting at the named critical transverse separation @xmath69 , stated when the dipole cross section reaches one half of its maximum value labeled by the pion - proton cross section . in connection with the present work ,
the glauber - mueller approach gives the eq .
( [ eq14 ] ) at the born level . instead of using an ad hoc control of the gluon distribution , the glauber - mueller provides corrections required by unitarity in an eikonal expansion . for the large @xmath10 region ,
we choose to follow a similar procedure from the gbw model , namely saturating the dipole cross section ( @xmath10-independent constant value ) .
having defined the notation and reviewed the main properties settled by the rest frame representation of deep inelastic process , in the next section we address the unitarity corrections formalism contained in the glauber - mueller approach , that we will consider in the calculation of @xmath19 , @xmath1 and @xmath2 .
the glauber formalism concerns mainly interactions with a target nucleus , allowing to calculate the amount of unitarity corrections to the nuclear cross section .
however this approach can be extended to take into account the evolution of the partonic densities ( saturation ) through the multiple scatterings .
below , we review shortly the main properties of the glauber formalism in qcd , either in the nuclear case as its application to the nucleon case .
we indicate the original papers @xcite for a complete presentation . it should be stressed that saturation effects mean unitarity corrections to the observables .
indeed , the asymptotic calculations have produced an unified @xmath70 pattern for the cross section and gluon function instead of a effectively saturated one @xcite .
since the small @xmath3 limit is driven by the gluonic interactions , we consider a virtual probe @xmath71 with invariant mass @xmath56 which decays into a gluon pair @xmath72 having a transverse separation @xmath48 and transverse momentum @xmath73 .
then , the pair interacts with the target through a gluonic ladder at fixed transverse separation . following the discussion from the previous section , the photoabsortion cross section for the probe particle in terms of @xmath3 and virtuality @xmath56 in the nuclear case is @xmath74 where the variables have the same identification as in the @xmath9 pair discussion .
the quantity @xmath75 is the light - cone wavefunction for the gluon pair .
the glauber s multiscattering theory employs the phase shift method to describe processes at high energy for an incident particle undergoing into successive scatterings .
when @xmath3 is small , the coherence length @xmath76 is bigger than the mean radius @xmath77 for all nuclei and the interaction of the initial parton stands through the entire nuclear path , providing coherent interactions with all target partons along the distance @xmath76 .
these scatterings are coherent in an interaction of a hadronic fluctuaction , i.e. quark or gluon pairs , with the nucleus .
there are interference effects among them , generating a reduction in the nuclear cross section , @xmath78 . otherwise , when @xmath79 , full incoherent scatterings occur leading to the expectation that the nuclear cross section equals to @xmath80 .
the well known glauber formula for the total cross section of a hadronic state with the nucleus is @xmath81 where @xmath82 is a profile function containing the dependence on the impact parameter @xmath83 , which is the conjugate variable to the momentum transfer @xmath84 .
it is related to the nucleon distribution inside the nucleus and encodes the information about the angular distribution of the scattering .
we discuss its particular shape later on .
( [ glauber ] ) is quite general and allows to describe the hadronic fluctuations of the photon or any virtual probe , as saw in the last section .
considering the scattering amplitude dependent on the usual mandelstan variables @xmath26 and @xmath84 , now written in the impact parameter representation @xmath83 , @xmath85 the corresponding total and elastic cross sections ( from optical theorem ) are rewritten in the impact parameter representation ( @xmath86 ) as @xmath87 a very important property when treating the scattering in the impact parameter space is the simple definition for the unitarity constraint @xcite . at fixed @xmath83
, the constraint can be expressed in the following way , @xmath88 with @xmath89 denoting the sum of contributions from all the inelastic channels .
the constraint above has a simple solution .
if the real part of the scattering amplitude vanishes at the high energy limit , corresponding to small @xmath3 values , the solution is @xmath90\,,\\ \sigma_{tot } & = & 2\,\int d^2{\bf b}\ , \left[\ , 1-\rm{e}^{-\frac{1}{2}\,\omega\,(s,\,{\bf b})}\ , \ , \right]\ , , \label{sigeik } \end{aligned}\ ] ] where the opacity @xmath91 is an arbitrary real function and
it should be determined by a detailed model for the interaction .
the opacity function has a simple physical interpretation , namely @xmath92 corresponds to the probability that no inelastic scatterings with the target occur . to provide the connection with the glauber formalism
, the opacity function can be written in the factorized form @xmath93 , considering @xmath94 normalized as @xmath95 ( for a detailed discussion , see i.e. @xcite ) .    from eqs .
( [ sigeik ] ) and ( [ glauber2 ] ) , we identify the opacity @xmath96 .
moreover , it has been found that the same formalism for multiple scatterings can be applied to the nucleon case . to proceed ,
we should determine the @xmath72 cross section .
the gluon pair cross section is equivalent to the quark pair one , up to a color factor ( @xmath97 ) .
the ( @xmath9 pair ) dipole - proton cross section is well known @xcite , calculated starting from eq .
( [ eq11 ] ) , and in double logarithmic approximation ( dla ) has the following form @xmath98 with the @xmath10-dependent scale @xmath99 .
considering eq .
( [ dipdla ] ) one can connect directly the dipole picture with the usual parton distributions ( gluon ) , since they are solutions of the dglap equations . in our case
, we follow the calculations in ref .
@xcite and consider the effective scale @xmath100 .
such a value differs from @xcite , where it is @xmath101 , which is obtained by an averaging procedure on the transverse size integral of @xmath1 .
however , in further studies in vector mesons it was found that @xmath102 ranges from 4 - 15 , and @xmath19 and @xmath1 are not sensitive to those variations .
thus , these values are consistent in leading logarithmic @xmath56 approximation .    from the above expression
, one obtains a dipole cross section satisfying the unitarity constraint and a framework to study the unitarity effects ( saturation ) in the gluon dglap distribution function .
hence , hereafter we use the glauber - mueller dipole cross section given by @xmath103    here , some comments are in order .
the glauber - mueller approach is valid in the small @xmath3 region , and the gluon emission is described in the double logarithmic approximation of perturbative qcd .
the interaction of the quark pair with the nucleon ( proton ) occurs through ladder diagrams exchange , satisfying the dglap equations in the dla limit .
the high energy limit allows to treat successive scatterings as independent collisions , meaning the process described by the eikonal picture of a relativistic particle crossing the nucleus .
moreover , as a consequence of no correlation among nucleons inside the nucleus ( in the nuclear case ) , there are no correlations among partons from different partonic cascades , stressing that only the fastest parton interacts with the target . the corrections coming from the slowest partons in the cascade ( emitted by the pair ) lead to the agl non - linear evolution equation @xcite , and they have been considered recently to describe diffractive dis in ref .
@xcite . regarding criticisms to the glauber - mueller approach
, we indicate the recent paper @xcite for a complete discussion about the eikonal - like models , concerning their advantages and limitations as well as pointing out the improvements to be taken into account to introduce the proper corrections .
now , we proceed to calculate numerical estimates of the dipole cross section using the glauber - mueller approach through eqs .
( [ dipdla]-[dipolgm ] ) .
then we are calculating saturation effects in the color dipole picture .
firstly , we need to discuss the profile function @xmath104 .
this function contains information about the angular distribution in the scattering , namely the @xmath84 dependence ( quark pair - ladder and proton - ladder couplings ) .
both of them can be approximated by an exponential parametrization , leading to a simple gaussian shape in the impact parameter space , @xmath105 , where @xmath106 is the atomic number and @xmath77 is the target radius .
we will keep this notation although we are only concerned with the nucleon case .
the @xmath107 value should be determined from data , ranging between @xmath108 gev@xmath109 for the proton case ( see discussions in @xcite ) . for nuclear reactions , a more suitable shape for the profile
should be taken into account , since the gaussian approximation is no longer appropriate to describe the nuclear profile for large @xmath106 . here
, we have used the value ( @xmath110 gev@xmath109 ) obtained from a good description of both inclusive structure function and its derivative @xcite .
such a value corresponds to significative unitarity corrections to the standard dglap input even in the current hera kinematics .
now , we discuss in a detailed way the main characteristics emerging from the dipole cross section eq .
( [ dipolgm ] ) . in order to do this , in fig .
( [ dipsigmas ] ) one shows the glauber - mueller dipole cross section as a function of dipole transverse size @xmath10 at fixed momentum fraction @xmath3 . for a better illustration on the expected partonic saturation effects
, we run the bjorken @xmath3 down to a quite small value @xmath111 ( thera region ) .
hereafter , we are using the grv gluon distribution at leading order @xcite in the input eq .
( [ dipdla ] ) , whose choice we justify below .
the solid lines correspond to the dipole cross section calculation , eq .
( [ dipolgm ] ) , whereas the dashed lines are the gbw model @xcite presented for comparison .
the general shape in terms of the dipole size comes from the balancing between the color transparency @xmath112 behavior and the gluon distribution .    here
some comments about the large transverse separation are in order : although perturbative qcd provides reliable results at small distances ( small dipole sizes ) , the nonperturbative sector is not still completely understood .
the usual pdf s are evoluted from a perturbative initial scale @xmath113 gev@xmath63 , and one has little information about the behavior at @xmath114 . in general
one makes use of regge phenomenology to estimate those contributions ( see , for instance @xcite ) .
thus , extrapolating to lower virtuality regions ( large dipole sizes ) one needs a parametrization regarding the nonperturbative sector .
this is the main justification of the use the grv94 parametrization @xcite in our calculations . bearing in mind that @xmath115 , its evolution initial scale is @xmath116 gev@xmath63 allowing to scan dipole sizes up to @xmath117 gev@xmath118 (= 0.62 fm ) .
for the most recent parametrizations , where @xmath119 gev@xmath63 ( @xmath120 fm ) the amount of nonperturbative contribution in the calculations should increase .
an additional advantage is that grv94 does not include non - linear effects to the dglap evolution since the parametrization was obtained from rather large @xmath3 values .
this feature ensures that the parametrization does not include unitarity corrections ( perturbative shadowing effects ) in the initial scale .
now , we should introduce an ansatz for the large transverse separation region .
a more phenomenological way is to match the pqcd dipole cross section with the typical hadronic one @xmath121 at @xmath122 as performed in @xcite .
however , due to the significant growth of the pqcd dipole cross section at high energies , we choose an alternative ansatz : the gluon distribution is frozen at scale @xmath122 , namely @xmath123 .
then , for the large distance contribution @xmath124 the gluon distribution reads as @xmath125 leading to the correct behavior @xmath126 as @xmath127 . in a more sophisticated case
, one can substitute the freezed scale @xmath128 by the saturation scale @xmath129 to take into account a realistic value of the gluon anomalous dimension in all kinematic region ( see correlated issues in @xcite ) .
recently , the phenomenological model of ref .
@xcite has produced a good description of hera data in both inclusive and diffractive processes .
it is constructed interplaying the color transparency behavior @xmath130 at small dipole sizes and a flat ( saturated ) behavior at large dipole sizes @xmath131 ( confinement ) .
the expression has the eikonal - like form , @xmath132 \,,\end{aligned}\ ] ] where @xmath133 gev@xmath63 and the three fitted parameters are @xmath134 mb , @xmath135 and @xmath60 and the notation for the saturation radius @xmath136 .
the gbw total cross section lies below the typical hadronic cross section : for instance in the pion - proton case , convoluting the pion wavefunction squared with the gbw dipole cross section we would have a constant cross section at high energies @xmath137 . despite describing data in good agreement
, gbw has some details that deserve some discussions : the approach does not present a dynamical hypothesis for the saturation phenomena and does not match dglap evolution . in gbw
, saturation is characterized by the @xmath3-dependent saturation radius @xmath138 instead of the scale coming from glauber - mueller , @xmath139 , which can be easily extended for the nuclear case @xcite .    in fig .
( [ dipsigmas ] ) one shows the glauber - mueller dipole cross section and the gbw model @xcite .
we choose to compare them due to the fact that gbw is actually a particular case of the glauber - mueller approach , considering an oversimplified profile function .
we should show below that this fact allows to construct an extended saturation model with dglap evolution @xcite ( bgk model ) . in the lower plots , where @xmath140 and @xmath141 , the gm cross section underestimates the gbw one .
however , as @xmath3 decreases the gluon distribution in the proton rises producing a bigger dipole cross section .
this feature is clear in the upper plots , for smaller @xmath142 and @xmath111 , where gm overestimates gbw by a significant factor mostly at intermediate @xmath10 values .
finally , we show the connection between the glauber - mueller approach with the saturation model with dglap evolution @xcite , observing that it is timely since the latter is quite efficient in describing data and its qualitative success corroborates quantitative qcd studies .
considering the particular case of central collisions , namely scattering at impact parameter @xmath143 ( @xmath144 ) , the glauber - mueller approach produces , @xmath145    the integration over @xmath146 can be promptly carryied out , and introducing the notation for the proper normalization for the dipole cross section , @xmath147 the eq .
( [ dipolgmb0 ] ) recovers the simple expression for the saturation model dglap evoluted @xcite , @xmath148\ , .
\end{aligned}\ ] ]    for a phenomenological analysis , the parameter @xmath55 and the scale @xmath149 are determined from data in @xcite . in our case
, it assumes the well defined value @xmath150 mb , using @xmath151 gev@xmath109 . for a larger radius , for instance @xmath152 gev@xmath109 ,
one obtains a value @xmath153 mb , closer to the gbw one . here , the virtuality scale is @xmath100 ) , whereas bgk choose the parametric form @xmath154 and a two - parameter initial condition for the gluon distribution function .
this section is devoted to the study and estimate of the gluon driven observables measured at hera kinematical domain in the rest frame .
the first one is the inclusive structure function @xmath156 , the main quantity testing the small @xmath3 physics .
the unitarity corrections are well established for this observable considering glauber - mueller approach @xcite as well as its high energy asymptotics @xcite , namely the black disk limit .
we review these issues considering the dipole picture ( rest frame ) , using a more complete analysis similar to @xcite , but mostly , discuss in detail the role played by the nonperturbative physics needed to describe the structure function , and where in the transverse separation @xmath10 range it starts to be important .
the longitudinal structure function @xmath157 is also addressed , verifying the frame invariance in comparison with previous laboratory frame calculations @xcite .
the longitudinal wavefunction strongly suppress large @xmath10 contributions , thus selecting smaller nonperturbative contribution in comparison with the @xmath19 case .
moreover , @xmath1 is one of the main observables scanning possible higher twist corrections in the standard operator product expansion ( ope ) @xcite .
therefore , a reasonable description of this quantity suggests that the glauber - mueller formalism ( or similar eikonal - like approaches ) take into account the most important contributions to the complete higher - twist corrections at current kinematical regimes .
the structure function @xmath158 gives the charm quark content on the proton and is directly driven by the gluon distribution . therefore it is a powerful observable to scan saturation effects in the small @xmath3 region .
however the current experimental status requires more dedicated measurements and a better statistics .
we verify a consistent description in the rest frame corroborating the similar analysis in the dipole models @xcite and in those ones considering unitarity corrections in the laboratory frame @xcite .
now we perform estimates for the inclusive structure function in the rest frame considering the glauber - mueller dipole cross section . the expression for @xmath19 , with the explicit integration limits on photon momentum fraction @xmath31 and transverse separation @xmath10 is , @xmath159    the notation has been introduced in the previous sections . in the fig .
( [ figf2 ] ) one shows the estimates for the structure function for representative virtualities @xmath56 from the latest h1 collaboration measurements @xcite .
the longitudinal and transverse contributions are shown separately , the longitudinal one being subdominant as is well known .
an effective light quark mass ( @xmath160 quarks ) was taken , with the value @xmath161 mev , and the target radius is considered @xmath110 gev@xmath109 , in agreement with ref .
it should be stressed that this value leads to larger saturation corrections rather than using radius ranging over @xmath162 gev@xmath109 .
the soft contribution comes from the freezing of the gluon distribution at large transverse separation as discussed at the previous section .
the gluon distribution considered is grv94 at leading order @xcite , @xmath163 , whose choice has been justified in the previous section .    from the plots we verify a good agreement in the normalization ,
however the slope seems quite steep .
this fact is due to the modeling for the soft contribution and it suggests that a more suitable nonperturbative input should be taken . indeed , in ref .
@xcite such a question is addressed , claiming that the correct input is the pion - proton cross section parametrized through the donnachie - landshoff pomeron .
it is found that the large transverse separations give a larger contribution at low @xmath56 , whereas it vanishes concerning higher virtualities .    to clarify this issue
, we can calculate the integrand @xmath164 used in the @xmath10-integration for the cross section @xmath165 , that should have significant nonperturbative content .
we plot in the fig .
( [ intpartial ] ) this quantity at low virtuality @xmath166 gev@xmath63 and at higher one @xmath167 gev@xmath63 for momentum fraction ranging on @xmath168 , verifying that the main contribution comes from an asymmetric peak at @xmath169 fm for @xmath166 gev@xmath63 , while it is shifted to @xmath170 at @xmath167 gev@xmath63 . in our calculation
, the perturbative contribution holds up to @xmath171 fm , therefore the region @xmath172 gives a non marginal contribution to the cross section at low virtualities .
indeed , one has found that it reaches about 10 % at @xmath166 gev@xmath63 and that when the virtuality increases the contribution gradually vanishes .
in fact , using the most recent pdf s this situation is more critical since @xmath122 is smaller ( @xmath173 gev@xmath63 ) .
this suggests that the photon piece @xmath174 multiplying the dipole cross section enhances the @xmath10-integrand to smaller @xmath10 at high @xmath56 , corroborating a similar conclusion already found in ref .
@xcite .    to clarify the role played by the soft nonperturbative contribution to the inclusive structure function and to verify the frame invariance of the approach , in the fig .
( [ f2comp ] ) we plot separately the perturbative contribution and parametrize the soft contribution introducing the nonperturbative structure function @xmath175 @xcite , which is added to the perturbative one .
the soft piece normalization is @xmath176 .
such procedure is done in order to compare explicitly the results found in ref .
accordingly , we have used just shadowing corrections for the quark sector , taking into account only the transverse photon wavefunction and zero quark mass .
the integration on the transverse separation is taken for @xmath177 , with @xmath116 gev@xmath63 for leading order grv94 gluon distribution .
this leads to a residual contribution absorbed in the the soft piece coming from transverse separations @xmath178 .
we considered the target radius being @xmath151 gev@xmath63 ( supported by @xcite ) , which produces a correction more important than the value @xmath179 gev@xmath63 .
it is again verified that the soft contribution is important at small virtualities and decreasing as it gets larger . in the plots ,
the dot - dashed lines represent only the perturbative calculations using the particular approximations indicated above , and the solid lines represent the results when we add the soft term .
concluding , we have a theoretical estimate , i.e. no fitting procedure , of the inclusive structure function @xmath0 through the glauber - mueller approach for the dipole cross section , detecting a non negligible importance of a suitable input for the large dipole size region .      as we saw in the sec .
( [ sec2 ] ) , the inclusive structure function can be expressed in terms of the cross sections @xmath180 and @xmath181 for the absortion of virtual photons transversally and longitudinally polarized , @xmath182 . at small @xmath3 ,
the longitudinal structure function is @xmath183 from qed , the longitudinal photons have zero helicity ( @xmath184 ) and therefore they have a virtual character . in the naive parton model ,
the helicity conservation for the electromagnetic vertex implies to the callan - gross relation @xmath185 and consequently a vanishing value for the longitudinal structure function @xmath186 , considering the scattering photon - quarks ( spin 1/2 ) . from qcd theory ,
this quantity has a non - zero value due to the gluon radiation , as is encoded in the altarelli - martinelli expression ( see discussion in @xcite ) @xmath187\ , , \label{martinelli}\end{aligned}\ ] ] where @xmath188 is the inelasticity variable .
therefore , the structure function @xmath1 is an auxiliary observable to detect saturation effects ( unitarity corrections ) in the gluon distribution .    experimentally , the determination of the @xmath1 is quite limited , providing few data points .
most recently , the h1 collaboration has determined the longitudinal structure function through the reduced double differential cross section @xcite , @xmath189 where @xmath190 . for large inelasticity
, the reduced cross section becomes @xmath191 and the contribution of @xmath1 is enhanced for large @xmath192 . the longitudinal structure function should be obtained only in the region of large inelasticity , covered in a large range at hera . in ref .
@xcite , two methods were used to perform the extraction : ( i ) for large @xmath193 gev@xmath63 , @xmath1 is obtained through the extrapolation method , using a nlo dglap qcd fit ( in the restrict kinematic range @xmath194 and @xmath195 gev@xmath63 ) to extrapolate @xmath19 into the high @xmath192 region .
( ii ) at low @xmath196 gev@xmath63 , the behavior of @xmath19 as a function of @xmath197 is obtained using the derivative method , based on the cross section derivative @xmath198 .
the data points obtained are consistent with the previous measurements , however they are more precise and lying into a broader kinematical range .
therefore , in the following we use only the new data points to analyse .    in fig .
( [ flrframe ] ) we present the estimates for the @xmath1 structure function , in representative virtualities as a function of @xmath3 . for the calculations ,
it was considered light quarks ( @xmath199 ) with effective mass @xmath200 mev and the target radius @xmath110 gev@xmath109 .
the large @xmath10 region is considered by the freezening of the gluon distribution at that region .
our expression for the observable is then , @xmath201    the behavior is in agreement with the experimental result , either in shape as in normalization .
a better description can be obtained by fine tunning the target size or the considered gluon distribution function , however it should be stressed that the present prediction is parameter - free and determined using the dipole picture taking into account unitarity ( saturation ) effects in the effective dipole cross section .
we verify that the rest frame calculation , taking into account the dipole degrees of freedom and unitarity effects produces similar conclusions to those ones using the breit system .
for instance , in a previous work @xcite , the unitarity corrections to the longitudinal structure function were estimated in the laboratory frame considering the eq .
( [ martinelli ] ) , with unitarized expressions for @xmath19 and @xmath202 , obtaining that the expected corrections reach to @xmath203 % as @xmath204 , namely on the kinematical sector of the upcoming thera project .
the higher twist corrections to the longitudinal structure function have been pointed out .
for instance , bartels et al .
@xcite have calculated numerically the twist - four correction obtaining that they are large for @xmath205 and @xmath1 , although with opposite signs .
this fact leads to remaining small effects to the inclusive structure function by almost complete cancellation between those contributions .
the higher twist content is analyzed considering the model @xcite as initial condition .
concerning @xmath1 , it was found that the twist - four correction is large and has negative signal , concluding that a leading twist analysis of @xmath1 is unreliable for high @xmath56 and not too small @xmath3 .
the results are in agreement with the simple parametrization for higher twist ( ht ) studied by the mrst group in ref .
@xcite , where @xmath206 .
the second term would parametrize the higher twist content . in our case ,
the unitarity corrections provide an important amount of higher twist content , namely it takes into account some of the several graphs determining the twist expansion ( for recent discussions in these issues , see @xcite ) .      in perturbative qcd ,
the heavy quark production in electron - proton interaction occurs basicaly through photon - gluon fusion , in which the emitted photon interacts with a gluon from the proton generating a quark - antiquark pair .
therefore , the heavy quark production allows to determine the gluon distribution and the amount of unitarity ( saturation ) effects for the observable .
in particular , charmed mesons have been measured at deep - inelastic at hera and the corresponding structure function @xmath158 is defined from the differential cross section for the @xmath207 pair production , @xmath208\,f_2^{c\bar{c}}(x , q^2)\,,\end{aligned}\ ] ] with @xmath192 being the inelasticity variable . in the laboratory frame , the dominant mechanism is the boson - gluon fusion @xmath209 .
hence , the charm structure function is directly driven by the gluon distribution and provides constraints for the gluonic function . in leading order ( lo ) , it is written as @xcite , @xmath210 where @xmath211 .
the mass factorization scale lies in the range @xmath212 .
such a scale introduces an uncertainty of about 10 % , and an additional source of uncertainty is the charm mass , in general ranging on @xmath213 gev .
the standard qcd coefficient function is labeled by @xmath214 .
experimentally , the latest measurements of the charm structure function are obtained by measuring mesons @xmath215 production @xcite . from the theoretical input ,
it was used nlo coefficient functions , considering charm mass @xmath216 gev and factorization - normalization scale @xmath217 .
the function @xmath158 shows an increase with decreasing @xmath3 at constant values of @xmath56 , whereas the rise becomes more intense at higher virtualities .
the data are consistent with the nlo dglap calculations . concerning the ratio @xmath218 , the charm contribution to @xmath19 grows steeply as @xmath3 diminishes .
it contributes less than 10% at low @xmath56 and reaches to about 30 % for @xmath219 gev@xmath63 @xcite .
once more the color dipole picture will provide a quite simple description for the charm structure function in a factorized way .
now , the glauber - mueller dipole cross section is weighted by the photon wavefunction constituted by a @xmath207 pair with mass @xmath220 .
our expression for the charmed contribution in deep inelastic is thus written as @xmath221 where @xmath222 is the probability to find in the photon the @xmath207 color dipole with the charmed quark carrying fraction @xmath31 of the photon s light - cone momentum with @xmath223 polarizations .
for the correspondent wavefunctions , the quark mass in eqs .
( [ eq9],[eq10 ] ) should be substituted by the charm quark mass @xmath220 . here
, we should take care of the connection between the regge parameter @xmath224 and the bjorken variable @xmath225 . for calculations with the light
quarks these variables are equivalent , however for heavier quarks the correct relation is @xcite , @xmath226    in fig .
( [ f2cfig ] ) we show the estimates for the charm structure function as a function of @xmath225 [ eq . ( [ xbj ] ) ] for representative virtualities . in our calculations , it was used the charm mass @xmath227 gev , the target size @xmath110 gev@xmath109 and the frozen gluon distibution at large @xmath10 . we have verified a small soft contribution as in the @xmath1 case , decreasing as the virtuality rises .
there is a slight sensitivity to the value of the charm mass , increasing the overall normalization as @xmath220 diminishes .
such a feature suggests that the charm mass is a hard scale suppressing the non - perturbative contribution to the corresponding cross section , which is in agreement with the recent bfkl color dipole calculations of nikolaev - zoller @xcite and donnachie - dosch @xcite .    regarding the breit system description , in ref .
@xcite it was found strong corrections to the charm structure function , which are larger than those of the @xmath19 ones making use of expression ( [ f2ccms ] ) .
considering the ratio @xmath228 , the corrections predicted by the glauber - mueller approach would reach to 62 % at values of @xmath229 ( thera region ) .
then , an important result is a large deviation of the standard dglap expectations at small @xmath3 for the ratio @xmath218 due to the saturation phenomena ( unitarization ) . with our calculation one verifies that it is obtained a good description of data in both reference systems , suggesting a consistent estimation of the unitarity effects for this quantity .
we study the dipole picture for the description of deep inelastic scattering , focusing on observables driven directly by the gluon distribution . starting from the dipole cross section provided by the glauber - mueller approach in qcd
, we perform estimates for the inclusive structure function @xmath19 , the longitudinal function @xmath1 and the charm structure function on the proton @xmath2 .
for each of the observables discussed , we obtain theoretical estimates , in the rest frame , without further fitting procedure , in good agreement with the updated data from hera .
the resulting calculations corroborate a quite consistent picture for the unitarity corrections from the glauber - mueller approach in both breit and rest reference systems . in the laboratory frame
the unitarity effects are connected with the gluon distribution function , whereas in the color dipole framework the basic block is the dipole cross section which is corrected considering saturation effects .
the small transverse separation @xmath10 region is dominated by the leading log dglap formalism , with the additional ingredient of unitarization phenomenon as the momentum fraction acquires quite small values .
such corrections are associated with the taming of the gluon distribution in the very small @xmath3 region , in general named saturation regime .
however , it should be stressed that the glauber - mueller approach and similar eikonal - like models provide a logarithmic @xmath70 asymptotic behavior for the inclusive structure function and gluon distribution , instead of a constant value for asymptotic energies .
the large transverse separation is described by non - perturbative aspects of qcd . since
this domain is not well determined at the moment , some modelling of the soft region is needed . in this work
we choose the anstz in which the gluon distribution is frozen for virtualities above a cut radius @xmath230 , which corresponds to the region @xmath231 . a convenient choice for the gluon pdf in order to cover the widest possible kinematical window
diminishes the uncertainty coming from the soft sector .
the most appropriated input is the grv94 parametrization , where @xmath232 fm is found , whereas it can take values @xmath233
fm for the more recent pdf s . throughout the paper we used the target size @xmath151 gev@xmath109 , which corresponds to strong unitarity corrections .
when performing a comparison with the phenomenological model gbw , we have found that the glauber - mueller approach underestimate the dipole cross section from gbw at not small @xmath234 .
instead , for very low @xmath235 the glauber - mueller overestimates gbw due to the strong increasing of the gluon function in this region . concerning the saturation model with dglap evolution ( bgk )
, it is a particular case of the gm approach when considering central scattering @xmath143 . despite
that the bgk model matches dglap evolution , glauber - mueller describes more properly the realistic impact parameter dependence of the process .
moreover , in gm the extension to the nuclear case is built in .    when considering the structure function @xmath19 , we have found that it is dominated by small transverse distances contributions .
however , a non - negligible content from the soft sector is present .
moreover , the photon wavefunctions enhances the dipole cross section into smaller dipole sizes , since the weight function selects smaller @xmath10 as the virtuality @xmath56 diminishes .
our estimates here are parameter - free , however a fine tunning of the parameters can improve the data description .
furthermore , we notice that in calculations from @xcite , only the aligned jet dipole configuration @xmath31 , @xmath236 ( and only transverse contribution ) is considered , whereas we take into account all configurations , including the symmetric ones .
thus , all dipole sizes , even those from the non - perturbative region are included in our results .
concerning @xmath1 , the estimates are consistent with the previous calculations in the breit system and are in good agreement with data .
a remarkable feature is that the glauber - mueller approach in the color dipole framework gives important higher twist contributions to the leading twist calculation in a simple way .
as is well known , @xmath1 is the main quantity to study the expected higher twist effects in low virtualities .
the function @xmath2 is directly dependent of the gluon distribution and important unitarity corrections had been found when considering the breit frame . here , we verify consistent results in the rest frame in comparison with the previous ones in the fast proton system .
we verified that the charm mass suppressed soft contributions in comparison with the @xmath19 case , and the results present a slight dependence with the specific value of @xmath220 .    in conclusion ,
the glauber - mueller approach provides a well stablished formalism to take into account the unitarity effects .
it allows to estimate the higher twist contributions to relavant observables in a simplified way , while matching dglap evolution equation at born level and including the impact parameter dependence properly .
mvtm acknowledges martin mcdermott ( liverpool university - uk ) , igor ivanov ( ikp - forschungszentrum juelich , germany ) and victor gonalves ( ifm - ufpel , brazil ) for useful enlightenements .
this work was partially financed by cnpq , brazil .