Investigation of proton structure function at HERA in light of an analytical solution to the Balitsky-Kovchegov equation

Ranjan Saikia ,Pragyan Phukan and Jayanta Kumar Sarma

1 HEP laboratory,Department of Physics,Tezpur University,Tezpur,Assam-784028,India

2 Department of Physics,Moran College,Moranhat,Assam-785670,India

Abstract In this paper,the proton structure functionat small-x is investigated using an analytical solution to the Balitsky–Kovchegov(BK)equation.In the context of the color dipole description of deep inelastic scattering (DIS),the structure functionis computed by applying the analytical expression for the scattering amplitude N(k,Y) derived from the BK solution.At transverse momentum k and total rapidity Y,the scattering amplitude N(k,Y)represents the propagation of the quark-antiquark dipole in the color dipole description of DIS.Using the BK solution we extracted the integrated gluon density xg(x,Q2) and then compared our theoretical estimation with the LHAPDF global data fits,NNPDF3.1sx and CT18.Finally,we have investigated the behavior ofin the kinematic region of 10-5 ≤x ≤10-2 and 2.5 GeV2 ≤Q2 ≤60 GeV2.Our predicted results forwithin the specified kinematic region are in good agreement with the recent high-precision data forfrom HERA(H1 Collaboration) and the LHAPDF global parametrization group NNPDF3.1sx.

Keywords: quantum chromodynamics (QCD),perturbative QCD (pQCD),proton structure function,small-x,Balitsky–Kovchegov equation

1.Introduction

Understanding the substructure of the nucleon within matter is one of the fundamental research topics in high-energy particle physics.It is crucial to understand the structure of the nucleon in order to comprehend the fundamental structure of matter.With the introduction of high-energy accelerator facilities,we have been able to understand the substructure of the nucleon within the framework of quantum chromodynamics(QCD).Understanding the substructure of the nucleon has relied heavily on the structure functions of the nucleon.Deep inelastic scattering (DIS) experiments on hadrons by leptons have yielded significant data on the distribution of partons within hadrons in terms of the quark and gluon distributions.The DIS cross section is associated with the structure functions of the nucleon in relation to parton distributions.The measurements of the proton’s structure functionsat HERA have begun a new era of parton density measurement within the nucleon [1–4].These structure functions can be correlated with the momentum distributions of partons within the nucleon,and thus the parton distribution functions (PDFs).At high energies,or equivalently at small-x(Bjorkenx) values,the gluon density dominates among partons,and hence the dominant contribution inobservations comes exclusively from gluons.Therefore,measurements of these structure functions at small-xare essential for calculating gluon distribution functions and visualising the overall hadronic wave function in high-density QCD.

The DGLAP (DokshitzerGribovLipatovAltarelliParisi)evolution equation [5–9],the famous and well-established evolution equation,acts as a basic tool for the theoretical investigation of DIS structure functions.This equation has been successfully applied to address available HERA data in the moderate kinematical region of Bjorken’sx(x≥0.01).At small-x,HERA data shows a steep rise in the behaviour of gluons in the region.This steep behaviour of gluons at small-xis well described by the famous BFKL (Balitsky–Fadin–Kuraev–Lipatov) evolution equation [10,11].However,the rapid growth of gluons at small-xcannot continue indefinitely;otherwise,the physical cross-sections will violate unitarity and the Froissart–Martin bound [12].As a result,both the DGLAP and BFKL evolution equations fail to account for implicit physics in high-density QCD.In order to maintain the unitarity of the theory,the infinite growth of gluons must be slowed down by certain processes.The phenomena of gluon recombination and saturation have provided the solution to the problems faced at small-xby the linear evolution equations [13–17].The nonlinear phenomena of gluon recombination and saturation lead to nonlinear terms in the DGLAP and BFKL equations.The Balitsky–Kovchegov(BK) equation [18–21],nonlinearization of the BFKL equation,and mean-field approximation to the Jalilian–Marian–Iancu–McLerran–Weigert–Leonidov–Kovner(JIMWLK)equation [22–25] describe the behaviour of gluon density in the small-xregion.The BK equation is an integro-differential equation in coordinate space that can be transformed to momentum space to yield a partial differential equation useful for phenomenological studies at various high-energy accelerator facilities.This equation has been studied and solved numerically;however,due to its complexity,we have not seen an exact analytical solution to this equation to date.We recently proposed an approximate analytical solution to the BK equation that can be useful for phenomenology in the context of gluon saturation and hence small-xphysics [26].

The collaborations at HERA contributed high-precision data for measuring the proton structure functionat different kinematical regions ofxandQ2[27,28].Some recent works on the measurement of proton structure functions using various evolution equations and approaches are described in[29–33],and the results correlate well with the experimental data.In [34],for the first time,the reproduction of measurements of the DIS proton structure function at high energy from the color dipole description in momentum space was investigated.Using the knowledge of asymptotic solutions of the BK equation,the authors of [34] measured the charm structure functionas a function ofxfor various values ofQ2and compared their results favourably with those of the HERA experiment.At small-x,it has been shown at HERA that the contributions in PDFs come exclusively from gluons.Thus,the determination of gluon density is essential in order to understand the overall hadronic wave function at small-x.We do not see any theoretical measurement of gluon density at small-xusing BK evolution theory,which motivates the current work.We previously presented an analytical expression for the scattering amplitudeN(k,Y),which describes the propagation of the quark-antiquark dipole through the target color field in the color dipole model at various rapidities,Y[26].In this work,we want to extract the gluon densityxg(x,Q2) and analyze the proton structure functionat small-xby combining experimental data with QCD evolution theory using the analytical solution of the BK equation in momentum space.We compute the gluon densityxg(x,Q2)and the proton structure functionat various kinematical regions of HERA and analyze their behavior at small-x.

The plan of the paper is as follows:in section 2,we relate the dipole-proton cross-section with forward scattering amplitude within the color dipole description of DIS in QCD.Then,using the analytical solution of the BK equation,we discuss how the structure functionof the proton can be obtained from the color dipole description.In section 3,we present thexevolution of the integrated gluon densityxg(x,Q2) for twoQ2values,viz.35 GeV2and 100 GeV2using the BK evolution theory.Our results are compared with the LHAPDF global data fits NNPDF3.1sx [35] and CT18[36].Both have LHC and HERA data included in their analyses.Finally,the obtained numerical results on the proton structure functionusing the analytical solution of the BK equation and comparison with the data from the H1 Collaboration and LHAPDF global parameterization group NNPDF3.1sx [35] are shown.Section 4 follows the discussion of the summary and the conclusion of our work.

2.Proton structure functions in the color dipole description

2.1.Color dipole description

In QCD,the color dipole description of DIS[37–39]has been a useful tool for various perturbative QCD (pQCD) calculations at small-xof high density QCD.It is especially convenient to work with the pQCD dipole description of DIS when it comes to the energy dependence of the scattering amplitude of any DIS event.The scattering process in the color dipole description of DIS can be factored into several steps where an incoming high-speed virtual photon (γ*) after fluctuation from QCD vacuum changes to a quark-antiquarkdipole.The formed dipole then scatters off the target proton,eventually forming the final-state particles.In the color dipole description,following [40],the total γ*pcross section for the γ*pscattering can be factored as

ΨT,Lrepresents the wave function for the virtual photon (γ*)to fluctuate into the quark-antiquarkdipole with a transverse size ofr.T,Lrepresent the transverse and longitudinal polarisation states of the virtual photon,respectively,andzrepresents the longitudinal momentum fraction of the virtual photon carried by the quark (or antiquark).Q2represents the virtuality of the photon,whileYrepresents the total rapidity.The dipole-proton cross section(r,Y),which is the imaginary component of the forward scattering amplitude of the dipole on the target proton,contains all the information on the hadronic interactions.

The dipole-proton cross section,(r,Y),can be obtained from the elastic dipole-proton scattering amplitude A (r,Y)using the optical theorem as [41,42]

N(r,Y,b) is the value for the imaginary part of the forward elastic dipole-proton scattering amplitude.Now,if one considers the target proton as a homogeneous disk of radiusRp,then the dipole-proton cross section can be related to the forward scattering amplitudeN(r,Y)by the following relation[34] :

The scattering amplitudeN(r,Y) will come from the solution of the BK equation.

After obtaining the γ*pcross-section,one can obtain directly the proton structure functionfrom the γ*pcrosssection through the relation

Now,we see how the structure functionFp2 can be obtained using the scattering amplitudeN(r,Y) obtained from the analytical solution of the BK equation.The BK equation describes the high-energy evolution of the dipole-target proton scattering amplitudeN(r,Y) in the dipole description,which in turn describes the propagation of the quark-antiquark dipole through the target color field.We express the γ*pcross-section inN(k,Y)and hence the structure function of the proton in momentum space,as the asymptotic behaviour of the solution of the BK equation can be expressed naturally in momentum space.For that,let us transformN(r,Y)toN(k,Y)by the following simple Fourier transform

The proton structure functionin momentum space related toN(k,Y) can be expressed as follows using the discussion from above and some algebraic calculations [34]:

where the photon wave function ˜Ψ is now expressed in the momentum space and can be found in [34] andNcare the number of colors.The scattering amplitudeN(k,Y) comes from the solution of the BK equation.The scattering amplitudeN(k,Y) at total rapidityYand transverse momentumkobeys the BK equation in the momentum space as [20]

Both the equations belong to the same universality class,and their solutions are of the same nature,i.e.travelling wave nature.We suggested an approximate analytical solution of the BK equation(7)in connection with the FKPP equation(8)using the homotopy perturbation method (HPM) in [26].Following [26],the solution of the BK equation (7) in connection with the FKPP equation (8) is given by

Nis the initial condition atY=0 i.e.N(k,0)=N.This solution of the BK equation gives the scattering amplitudeN(k,Y) at any given rapidityY＞0,once the initial conditionNis known to us.For the initial conditionN,we use the GBW(Golec–Biernat and Wusthoff) initial condition given by [48]

Qis the squared value of the initial saturation momentum of gluons that can be fitted from the existing HERA data,and its value is 0.24 GeV2[49].As we are dealing with the the BK equation in momentum space,the use of the GBW initial condition would be helpful as this can be transformed into momentum space simply,resulting in

Therefore,we can write the GBW initial condition (10) in momentum space as

Now,we use the above equation in equation (9) and replace the initial conditionNwith the GBW initial condition in momentum space,which gives

The scattering amplitude expression given above is an approximate analytical solution to the BK equation (7).

Figure 1.The probability distribution of the virtual photon to emit a quark-antiquark pair as a function of k with mq →0 and z for various Q2.

In light of the above discussions,we will now investigate the proton structure functionFp2 by combining experimental data from the HERA experiment with QCD evolution theory.To do this,we will use the analytical solution of the BK equation,which is expressed above.

3.Results and discussions

To investigate the proton structure function,we use the expression equation(6)with the analytical solution of the BK equation (14).We handle the expression equation (6)numerically for various inputs ofxandQ2.In equation (6),)z;Q2represents the probability of a virtual photon emitting a quark-antiquark pair with the momentum fractionz(quark) and (1-z) (antiquark) of the virtual photon in the momentum space.Its expression is given by [34]

The solution (14) gives the propagation of thedipole in the color dipole description of QCD,which in turn gives the gluon content of the proton,i.e.the unintegrated gluon density.To calculate conventional gluon densityxg(x,Q2),we integrate the BK solution(14)over transverse momentum with the relation

The extracted results of the gluon density from the BK solution are compared with the global data fits NNPDF3.1sx and CT18.The results are shown in figure 2.

For the proton structure functionwe plug the expression given in equation (15) (withmq→0 andz=1/2)together with the BK solution given in equation (14) into equation(6).We setNc=3 andRp≈0.831fm≈4.22 GeV-1in the expression equation (6) and solve the expression numerically for variousQ2.The value of theRpis taken from the recent study on the proton radius [50].We analyse the expression equation (6) numerically and compare our results with the HERA measurements of the proton structure function from the H1 Collaboration [28] with constraints to the kinematic region: 10-5≤x≤10-2and 2.5 ≤Q2≤60 GeV2.We also compare our results with the LHAPDF global parameterization group NNPDF3.1sx [35].The reason for choosing the particular kinematic regionx≤10-2is to describe the small-xbehaviour of high-energy amplitude as the BK equation is only applicable at small-x.For theQ2range which is too high,we would need corrections from the DGLAP equation,which we cannot skip atQ2which is too high.The results are shown in figure 3.

4.Summary and conclusion

In this work,we have investigated the proton structure functionin light of an analytical solution of the BK equation.Proton structure functions have been investigated by different collaborations at experimental facilities such as HERA and the LHC.On the phenomenological side,the proton structure functions have been investigated by studying different QCD evolution equations and well tested at existing experimental facilities.To investigate the proton structure functions at small-x,the BK equation is most suitable for testing the experimental data on observables.In this work,we have stated our discussion on how we can obtain the proton structure functionfrom the color dipole description of the DIS in QCD.In the color dipole approach,the proton structure functionis obtained directly from the virtual photon-proton cross-section,which can be expressed in the dipole-proton scattering amplitudeN(k,Y).The dipole-proton scattering amplitudeN(k,Y) is obtained from the solution of the BK equation.Using the BK solution,we extracted thexevolution of the integrated gluon densityxg(x,Q2) and compared our results with those of LHAPDF global data fits NNPDF3.1sx and CT18.Both NNPDF3.1sx and CT18 have HERA and LHC data included in their analyses.Finally,we computed the proton structure functionand compared it with the measurement ofat HERA from the H1 Collaboration and the LHAPDF global parameterization group NNPDF3.1sx.Our predicted results are in good agreement with the experimental results within the kinematic region we have constrained.Beyond that,we have to consider corrections to the DGLAP equation.We have also shown in this work how the virtual photon wave function Ψ would behave against transverse momentumkwithmq→0,z=1/2 at variousQ2.We found that the probability of the virtual photon emitting a quark-antiquark pair increases askincreases to a certain peak value,after which it starts to fall askincreases.We have seen that the maximum probability is the same for differentQ2values,irrespective of the value ofkat the maximum peak for differentQ2values.

Figure 2.x evolution of gluon density xg(x,Q2)extracted from the BK solution compared with the global data fits NNPDF3.1sx and CT18.

Figure 3.The results of the proton structure function as a function of x at various Q2 obtained in this work are compared with the data from the H1 Collaboration [28] and global data fit NNPDF3.1sx [35].

In this work,we have seen the ability of our BK solution to describe physics at small-x.We have successfully applied our BK solution to calculate the gluon densityxg(x,Q2) and the proton’s structure functionat small-xwithin the kinematic region we have constrained.Also,our results are testable at future experimental facilities such as the LHeC(Large Hadron electron Collider) [51,52],EIC (Electron Ion Collider) [53],and the FCC-eh (Future Circular Collider electron-hadron) [54].In these future experimental facilities,the measurement of the proton structure function will be performed at much lower values ofxwith increased precision.Nevertheless,we could investigate the proton structure functionusing an analytical solution of the BK equation within the constrained region.We conclude that the analytical solution of the BK equation can serve as a convenient tool for further studies at small-xand high-density QCD.We hope that the BK equation,together with future experimental facilities,will help us understand and explore phenomena inside hadrons at small-xin the near future.

ORCID iDs