On x-independence of RQ = FLQ/F2Q ratio at low x

We give predictions for the ratio RQ(x,Q2)=FLQ(x,Q2)/F2Q(x,Q2) at small values of Bjorken variable x in the first three orders of perturbation theory. We show an approximate x-independence of RQ(x,Q2) at low x and non-large Q2 values (Q2<8--10 m2Q), irrespectively on the gluon density in a proton used in the calculations. This observation could be useful in subsequent phenomenological studies of the heavy flavor production at future lepton-hadron and hadron-hadron colliders.


Introduction
An important new data on the cross sections for the open charm and beauty production in deep inelastic electron-proton scattering (DIS) have appeared [1] by combining the results of research from the H1 and ZEUS collaborations at HERA.Measurements have shown that the production of heavy flavor Q in DIS occurs predominantly due to the photon-gluon fusion process γg → QQ and therefore depend strongly on gluon density in the proton and mass m Q of produced heavy quark 1 .Theoretical studies usually serve to confirm that available HERA data can be described by perturbative heavy flavor generation in QCD (see, for example, review [2] and references therein).Further investigations are planned at future lepton-hadron and hadron-hadron colliders [3], such as LHeC, FCCeh and FCC-hh (see also [4] for review), where the measurements can be performed with much increased precision and extended to much smaller x and high Q2 values.
In our previous consideration [5,6] we have analyzed latest experimental data [1] taken by the H1 and ZEUS Collaborations at HERA.In particular, we have studied the heavy quark contributions to the proton structure function (SF) F 2 (x, Q 2 ) and reduced charm and beauty cross sections σ cc red (x, Q 2 ) and σ bb red (x, Q 2 ) measured mostly at small values of the Bjorken variable x in a wide region of Q 2 [7].An additional important result of the evaluations [6] is that the compact analytical expressions for these DIS coefficient functions have been presented up to next-to-next-to-leading order (NNLO) accuracy.These expressions were used very recently to investigate the top quark production in the FCC-he kinematical regime [8].Here we continue our study [5,6].Using the previously derived expressions for DIS coefficient functions, we study the ratio in the first three orders of perturbation theory We demonstrate an approximate x-independence of this ratio for non-large Q 2 values, namely, Q .Moreover, we show a very slow ratio's dependence on the choice of used gluon density.

Basic formulas
It is well known that the differential cross section σ Q Q of heavy flavor production in DIS can be presented in the simple form: where x and y are the usual Bjorken scaling variables.Suggesting that the SFs ) are driven at small x by gluons (see [2]), it is easily to obtain the following results (see, for example, [9,10]): where the symbol ⊗ marks the Mellin convolution: So, the ratio R Q (x, Q2 ) can be presented as where f g (x, Q 2 ) is the gluon density in a proton 2 .In fact, the ratio R Q (x, Q 2 ) depends slowly on non-perturbative input f g (x, Q 2 ), which contribute to the both numerator and denominator of the ratio As we already noted in Introduction, we study the ratio R Q (x, Q 2 ) at small values of the Bjorken variable x.Using the results [13][14][15], below we give formulas for the high energy asymptotics of collinear coefficient functions of the heavy quark production process in the first three orders of the perturbation theory.

Wilson coefficients
Taking the Wilson coefficient in the m order of perturbation theory, we have ) where e Q is the charge of heavy quark Q and a s (µ 2 ) = α s (µ 2 )/(4π) is a strong coupling constant.Moreover, µ is the renormalization scale, which is usully taken in the following form [2,4] 1.The LO Wilson coefficients have the following form: where So, for the ratio R Q LO (x, Q 2 ) we have: 2. The NLO coefficient functions B k,g (x, a, µ) of photon-gluon fusion subprocess are rather lengthy and not published in print; they are only available as computer codes [16].Following [9,10], it is sufficient to work in the high energy regime, defined by x 1, where they assume the compact form3 [13][14][15]: with and where being the dilogarithmic function.We would like to note that B So, at NLO we have: 3. In the high energy regime the coefficient B k,g (x, a, µ) has the compact form: with where (14), respectively, and where t is given in (14) and are the trilogarithmic function Li 3 (x) and Nilsen Polylogarithm S 1,2 (x) (see [17]).The results for K(a) in the form of harmonic Polylogarithms [18] can be found in [15].So, at NNLO we have:

Gluon density in a proton
In the previous paper [6], the gluon density f g (x, Q 2 ) coming from the generalized double asymptotic scaling (gDAS) approach [19,20] has been applied.The usage of the gluon density leads to a good description of experimental data for the SFs F 2 (x, Q 2 ) and F c 2 (x, Q 2 ) (see, for example, [20,21] and [9,10], respectively).Moreover, the parton densities of the generalized DAS approach have been used successfully [5,6,22] to generate the TMD parton densities (see the recent paper [23]), applied in the k T -factorization approach [11,12], which is very useful in partially to study F Q 2 (x, Q 2 ) and the corresponding (reduced) charm and beauty cross sections σ cc red and σ bb red .Here, with a purpose to implement more information from experimental data, additionally we will use the gluon density, corresponding to H1 parametrization [24] (see also [25,26]) for the SF F 2 (x, Q 2 ): with the following parameters [24] C H = 0.18, a H = 0.048, Λ H = 292 MeV 2 .( Assuming the basic contribution of f g (x, Q 2 ) to the SF F 2 (x, Q 2 ) at small x values, we can suggest it in the following form with some unknown C g H , which value is not so important for the present investigations.

Resume
The predictions are shown in Figs. 1 -3.Our calculations show a flat (x-independent) behavior of R Q (x, Q 2 ) for small values of x (x ≤ 10 −2 ) and at not very large For larger values of Q 2 , some dependence on x appears, especially in the NNLO case, which has strongly singular coefficient functions ∼ 1/(n − 1) 2 and leads to the strengthening ∼ ln x in the NNLO contributions at low x (see Appendix B in Ref. [6]).In the case of c and b quarks, such flat behavior is in agreement with one observed earlier [6,27] within the framework of the k T -factorization approach [11,12].In the case of t quark, this flat behavior was recently discussed in the second paper of Ref. [8].Very large NNLO contribution to R t (x, Q 2 ) for Q 2 4m 2 t indicates a strong need for resummation for small x (see, for example, [13,14,28,29] and discussions therein), but that is beyond the scope of this short paper.
As one can see in Figs. 1 -3, the obtained results for R Q (x, Q 2 ) are very similar for both gluon densities even though they have very different low x asymptotics.So, this shows a very slow dependence of the ratio R Q (x, Q 2 ) on used gluon density.It is interesting that such flat behavior observed for R Q (x, Q 2 ) at the first three orders of perturbation theory strongly contradicts the behavior of the ratio R = F L /(F 2 − F L ), which becomes strongly x-dependent in higher orders of perturbation theory (see e.g., [2]) An attempt to resum the large terms in the perturbation theory coefficients leads [30] to large values of the argument of the strong coupling constant.
We hope that the observed approximate x-independence of R Q (x, Q 2 ) at low x will be useful in further phenomenological studies of the (reduced) cross sections σ QQ red and SFs F Q 2 (x, Q 2 ) at future lepton-hadron and hadron-hadron colliders, such as LHeC, FCC-eh and FCC-hh.

Figure 1 :
Figure 1: The ratio R c (x, Q 2 ) as a function of x at different Q 2 values.The purple (1), azure (3) and yellow (5) curves correspond to the results obtained with AKL gluon density function[23] at LO, NLO and NNLO, respectively.The green (2), orange (4) and blue (6) curves correspond to the results obtained with H1 parametrization from Eq. (25) at LO, NLO and NNLO, respectively.

Figure 2 :
Figure 2: The ratio R b (x, Q 2 ) as a function of x at different Q 2 values.Notation of curves is the same as in Fig. 1.

Figure 3 :
Figure 3: The ratio R t (x, Q 2 ) as a function of x at different Q 2 values.Notation of curves is the same as in Fig. 1.