RETRACTED
This article has been retracted as it was found to contain a substantial amount of material from following published articles: (1) “Evolution equations for truncated moments of the parton distributions,” Phys. Lett. B644 (2007) 284-287; [hep-ph 0610282], (2) “Truncated Mellin Moments: Useful Relations and Implications for the Spin Structure Function g₂,” Acta Phys. Polon. B42 (2011) 1231-1246; [arXiv: 1106.3753].

Advances in High Energy Physics
Volume 2014 (2014), Article ID 304369, 7 pages
http://dx.doi.org/10.1155/2014/304369

Research Article

Extracting the QCD Cutoff Parameter Using the Bernstein Polynomials and the Truncated Moments

A. Mirjalili,¹ M. M. Yazdanpanah,² and Z. Moradi¹

¹Physics Department, Yazd University, P.O. Box 89195-741, Yazd, Iran
²Faculty of Physics, Shahid Bahonar University of Kerman, Kerman, Iran

Received 31 August 2014; Revised 10 October 2014; Accepted 13 October 2014; Published 11 November 2014

Academic Editor: Elias C. Vagenas

Copyright © 2014 A. Mirjalili et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP³.

Abstract.

Since there are not experimental data over the whole range of -Bjorken variable, that is, , we are inevitable in practice to do the integration for Mellin moments over the available range of experimental data. Among the methods of analysing DIS data, there are the methods based on application of Mellin moments. We use the truncated Mellin moments rather than the usual moments to analyse the EMC collaboration data for muon-nucleon and WA25 data for neutrino-deuterium DIS scattering. How to connect the truncated Mellin moments to usual ones is discussed. Following that we combine the truncated Mellin moments with the Bernstein polynomials. As a result, Bernstein averages which are related to different orders of the truncated Mellin moment are obtained. These averaged quantities can be considered as the constructed experimental data. By accessing the sufficient experimental data we can do the fitting more precisely. We do the fitting at leading order and next-to-leading order approximations to extract the QCD cutoff parameter. The results are in good agreement with what is being expected.

1. Introduction

Deep inelastic scattering (DIS) provides us with a facility to improve our knowledge for subatomic structure and check the precision of quantum chromodynamic (QCD) theory at the parton level. To do the phenomenological task, it is inevitable for one to use the experimental data, resulting from DIS for parton densities and the nucleon structure functions. To make the calculations more reliable, one is forced to apply the limits for the interval of -Bjorken variable which are imposed by experimental data. This means that we are not able to use the usual and current definition of Mellin moment for parton densities or nucleon structure functions to do the required calculations. This situation automatically leads us to define and use the truncated moments [1–5]. One can find in [1] the relation which exists between the usual moments and the truncated moments that we also use in this paper.

On the other hand, one can employ the Bernstein polynomials to restrict the computations to the regions where there are available experimental data. In [6] one can find that the Bernstein averages of the nucleon structure function can be written as a linear combination of the moments. The truncated moments depend on a Bjorken variable. This variable is in fact the commencing point of an experimental data interval. This point is related to the lower limit of integration to obtain the truncated moments. These moments depend also on the used energy scale. If we go one step further and do another Mellin transformation we can get rid of any dependence on the Bjorken variable. This yields to us a new moment which is called through this paper “moment of moment.” As it was pointed out this quantity can be related to usual moments [1]. So it is possible to use this quantity to construct its Bernstein averages. By changing the order of “moment of moment,” it is possible to get more and more Bernstein averages which provide us with more experimental data dependent input to do a precise fitting.

The organization of this paper is as follwos. In Section 2 we give a review to do the evolution of Mellin moments in the leading order (LO) and next-to-leading order (NLO) approximations. Section 3 is devoted to how the evolution of truncated moments can be done. We introduce the Bernstein polynomials in Section 4. We find the relation between moment of truncated moment (moment of moment) and the related Bernstein averages in Section 5. We use the relation which was obtained in previous section and do a fitting in Section 6 to extract the QCD cutoff parameter at the LO and NLO approximations. Finally we give our conclusion in Section 7.

2. An Overview of Parton Densities Evolution in Mellin Moment Space

The evolution of parton densities with respect to transferred momentum is given by [7–13]

In this equation the symbol is representing a convolution integral so as

In (1), is representing the parton densities and . The is QCD cutoff parameter which is determined by fitting the experimental data. This parameter can appear, for instance, in the renormalized running coupling constant at the NLO approximation which is given by

Here and are the two universal coefficients of QCD- function. In (3) we mean . In (1), is the splitting function which has the following expansion:

where and are representing the splitting coefficients at the LO and NLO approximations, respectively.

Now let us have a brief review on how we can solve the integrodifferential of DGLAP equation in Mellin moment space. The employed method to get the solution will also be used when we intend to get the solution for truncated Mellin moments.

In this regard we should use the Mellin transformation:

By doing the Mellin transformation on (1), we will arrive at

where is the anomalous dimension which is being obtained, doing the Mellin transformation on the splitting function. It has the following expansion:

In order to get the solution of evolution equation in Mellin moment space we need as well to use the QCD -function which has the following form:

and we then get

Considering (3) and (7) at the LO approximation, we can write

Substituting (9), (10), and (11) in (6), we will obtain

which is known as the solution of evolution equation at the LO approximation in the Mellin moment space.

To extend the calculation to the NLO approximation we should do the following. Instead of (9) and (10) we will have

Substituting (14) into (6) will take us to the solution of Mellin moment of the parton densities at the NLO approximation:

Equation (15) can be represented as below, using the exponential and Taylor expansion:

By accessing the parton densities at initial energy scale, , we can obtain its evolved Mellin moments with respect to different values of . We plot in Figure 1 the first and second moment of parton densities versus , using (15). We use the valence parton density in [14] at as the input parton density.

Figure 1: First and second moment of valence parton density versus values at the LO and NLO approximations.

In order to obtain the parton densities in the -Bjorken space, one can use the inverse Mellein transformation as follows:

Different methods exist to get the parton densities in space. One of them can be found in [15]. There are different sets of parton densities in -space. Throughout the whole work we will use the ones of [14].

3. Evolution of the Truncated Mellin Moments

What we did in previous section is based on integration of parton densities over the whole range of variable, starting from . But it does not practically happen, since in the region where the partons do not have any contribution in the carried momentum by nucleon. Furthermore in this region the invariant energy of hadronic system in DIS is raising out to infinity which is physically forbidden. So it is inevitable for one to consider the lower limit of the related integration at the initial value which is different from zero. We should also note that there are not any experimental data for the mentioned region which is again another evidence not to take the lower limit of integration to equal zero. In this case the usual Mellin moment appears in another format which is called truncated Mellin moment [16]:

If we do an integration like the one in (18) on the DGLAP evolution equation, (1), then we obtain

where is representing the related splitting function of the truncated moments and is given by

In order to obtain the solution of (19) we need to get the relation which exists between the usual moments and the truncated moments. Since there is similarity between the evolution equation of usual moments and the truncated moments, that is, (1) and (19), we can use the same technique which is used for (1) to solve the evolution equation for the truncated moments. For this purpose we do Mellin transformation on (19) and call the result as “moment of moment (MM)” which is given by

By doing the Mellin transformation on (19) we get the following equation:

where

Substituting (20) into (23) we achieve a quantity which is analogous to anomalous dimension but it is in fact the Mellin moment of the truncated splitting function and is given by

The “MM” quantity, given by (21), does not have by itself any usage. In order to employ it in practice, it should be written in terms of the usual Mellin moments which appears as

Now we can solve evolution equation (19) as we did for usual evolution equation, (1), just by considering this point that during the calculations we should use the “moment of moment” instead of usual moment.

Therefore following the same strategy which leads us to the solution of evolution equation, (1), we can obtain the following results as the solution of the evolution equation for the truncated moments, (19), at the LO and NLO approximations, respectively:

4. Bernstein Polynomials

The Bernstein polynomials are used as a weight function to get the average of a function, over the interval [17, 18] in which we can write

where is the concerned averaged quantity. The Bernstein polynomials have the following representation:

The required quantities which are sited in the concerned interval are given by

To ensure that the integral in (27) is equivalent to this integral over we should divide this equation to a normalization factor which is given by .

If we do the Taylor expansion for the Bernstein polynomials we then obtain

or equivalently

Substituting (31) into (27), the averaged quantity, , can be represented by

In terms of the Mellim moments, (32) can be represented by

5. Truncated Moments and the Brenstein Plolynomails

The truncated moments of parton densities are functions of the carried momentum fraction and the scale energy which are represented by and (or ), respectively. Therefore we can write

or equivalently (see (21))

where, in (34), is the truncated moment which is defined by (18) and in (35) is its related moment which we called before MM (“momet of moment”). We mean from the order of “MM,” the number in (35). According to (35) the average of the truncated moment is related to a sum over its “MM.” The advantage of this equation is that by increasing the order of “MM,” that is, , we can access more averaged quantities of the related truncated moment. This enable us to construct more data dependent input which will improve the precision of the fitting.

6. Extracting the QCD Cutoff Parameter

Here we are intending to show how we can use (35) to extract the QCD cut-off parameter, . First of all we need to determine the digit amounts of and for the averaged quantity . For this purpose, considering the EMC collaboration data for muon-nucleon DIS scattering [19] and WA25 data for neutrino-deuterium DIS [20] at , we specify the minimum and maximum values of the variable. Then using (29) we determine the amount of and which covers the interval. We can also determine the and by paying attention to the behaviour of Bernstein polynomials with respect to which leads to similar result. As can be seen from Figure 2 the peaks of the plots are moving by changing the and values. We choose the in which its peak is placed in the interval range of variable where the population of the experimental data is large. By specifying the indices of the Bernstein polynomials we can calculate the in (35) as the averaged quantity.

Figure 2: The Bernstein polynomials with respect to for two different and values.

By determining and we should first obtain in (35), based on the two theoretical and experimental methods. To follow the theoretical method we go back to Section 3 which gives us the evolution of truncated Mellin moments. To construct the truncated moment, according to (18) we need the parton density function. We get this function from [14] and take the -valence parton density at as our initial input. We then evolve the truncated moment to energy scale using (26) which gives us the evolved truncated moments at the LO and NLO approximations, respectively. Therefor considering the determined and values at a specified energy scale , using (35), the as a function of the order of truncated moment, , can be obtained. We should note that the mathematical property of Bernstein polynomials compels having always [21].

Now we illustrate how to obtain the truncated Mellin moment, , experimentally. First we fit a phenomenological function to the available experimental data for the -valence quark distribution which are here at . As a result the numerical values of , , and parameters are determined. Now the “MM,” that is, , can be obtained experimentally, using (25). Following that, using (35), the averaged quantity is also determined in terms of the linear combination of “MM” at a specified order while the numerical values of and are determined at given value of .

We can now go to a further step and determine the parameter. Back to (35) we see that we have the from two stands. On the one hand, we achieve this function theoretically where all quantities are specified except the parameter which sits on the definition of renormalized coupling constant through “MM” function in (35) (see (26) and (27)). On the other hand, we also access the experimental values of at specified which is here 15  while the ad are determined. By changing the order of “MM,” that is, , we can construct experimental values for this quantity at different order. Fitting these numerical values to the analytical function of while there is unknown parameter will yield to us the numerical value for this parameter.

To do more precise fitting we need to increase the produced experimental data. This can be done by changing the numerical values for and . For this purpose we use (27) to get the averaged amount of -valence quark distribution at . This can be done by substituting in (27) with the related phenomenological quark distribution as a function of -Bjorken variable and . Now we can change the amount of and indices. If the numerical amount of the Bernstein average for the valence density, , is near to its experimental value in the concerned interval of -variable then the chosen values for the and are acceptable. We can extend this procedure to find more acceptable values for the and . For each value of the and we can repeat the procedure described in the two previous paragraphs to find the produced experimental data, , at the specified order of . It is obvious that by changing the order of “MM” we can construct and produce more experimental data for the at the desired and values.

Therefore we have some experimental data which contain different values of and as well as different numerical values for the order of the “MM. ” By doing a global fit we vary in our fitting the amount of , , and . Following that we are able to extract the numerical values for the at the LO and NLO approximations, respectively. What we got here are and , where the four quark active flavours, , are used. Let us call these numerical values as the “central values of the extracted results for ” since they are obtained using the constructed Bernstein averages rather than directly the experimental data. The result is corresponding to what is expected for the QCD cutoff parameter [22, 23]. We should notice that the results of [23] have been obtained at five quark active flavours, , by taking the Z-boson mass. But what is reported there is comparable with what we got for the at the NLO approximation. We should also note that the NLO value of extracted in [23] can be compared with our results for the central values of the extracted only after transforming the result of [23] from to , using the matching condition for quark threshold on the running of the QCD coupling [24, 25]. Meantime the NNLO result of the analysis in [23] may contain additional theoretical uncertainties. For more discussion of this subject see [26]. We can also see that the numerical values for at the NLO approximation are greater than the LO one which is an acceptable behaviour in all fitting procedure. We depict in Figure 3 the plot of at for different values of and as a function of the MM order, . As can be seen, the agrement with the produced experimental data is acceptable and the NLO result is better than the LO one.

Figure 3: The averaged quantity of the truncated Mellin moment as a function of the truncated moment order, .

7. Conclusion

Since there is not experimental data for the whole range of the -Bjorken variable, one should restrict the limit of integration for the Mellin moment transformation. This will lead us to define the truncated Mellin moments. Based on the DGLAP equation, the evolution of (usual) Mellin moments can be done. Here in this work we were looking to obtain the evolution equations for the truncated moments. The truncated Mellin moments are related to usual moments, using a new quantity which was called moment of moment (MM). In fact this quantity is the Mellin transformation of the truncated moment. “MM” behaves like the usual moments by this difference that contains two indices, , to indicate the order of the concerned moment. So it can be evolved as usual moment. However in practice we only considered one of these indices, , as the order of the moment.

On a further step we obtained the average of truncated moment using the Bernstain polynomials. In this case the average of truncated moment can be written in terms of linear combination of the “MMs”. This made a possibility to do a fitting and obtain the . For this purpose we used the experimental data for the -valence of parton density at  . We got first the phenomenological parameterizations for this data. Then we obtained its truncated moment and in continuation we would also obtain the moment of this truncated moment which was called “moment of moment.” This means that we were able to produce the experimental data for the average of the truncated moment by changing its order.

Two free indices which exist in definition of the Bernstein polynomials can be determined so as to give us the averaged amount for the valence quark densities near the available experimental data over the concerned interval of the -Bjorken variable. In other words, by suitably choosing , we manage to adjust the region where the average is peaked to that in which we have experimental data.

To do the fitting we made a global fit, using the produced experimental data which was related to the average of truncated moments and extracted the cutoff parameter of QCD at the LO and NLO approximations which corresponded to what one was expecting.

The technique which we used in this paper to do the fitting can be extended to the polarized nucleon structure function. The advantage of this method will appear more in this case since there are few experimental data for the polarized nucleon structure function. But our method makes this possibility to construct more experimental data dependent input by changing the order of truncated moment. We hope to report about this issue in our further research task.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

S. Forte and L. Magnea, “Truncated moments of parton distributions,” Physics Letters B, vol. 448, no. 3-4, pp. 295–302, 1999.

S. Forte, L. Magnea, A. Piccione, and G. Ridolfi, “Evolution of truncated moments of singlet parton distributions,” Nuclear Physics B, vol. 594, no. 1-2, pp. 46–70, 2001.

D. Kotlorz and A. Kotlorz, “Evolution equations for truncated moments of the parton distributions,” Physics Letters B, vol. 644, no. 4, pp. 284–287, 2007.

D. Kotlorz and A. Kotlorz, “Evolution equations of the truncated moments of the parton densities. A possible application,” Acta Physica Polonica B, vol. 40, no. 6, pp. 1661–1671, 2009.

D. Kotlorz and S. V. Mikhailov, “Cut moments and a generalization of DGLAP equations,” Journal of High Energy Physics, vol. 2014, no. 6, article 065, 2014.

C. J. Maxwell and A. Mirjalili, “Direct extraction of QCD from moments of structure functions in neutrino-nucleon scattering, using the CORGI approach,” Nuclear Physics B, vol. 645, pp. 298–308, 2002.

Y. L. Dokshitzer, “Calculation of the structure functions for deep inelastic scattering and e+ e- annihilation by perturbation theory in quantum chromodynamics,” Soviet Physics-JETP, vol. 46, pp. 641–653, 1977.

Y. L. Dokshitzer, “Perturbational calculation of the deep inelastic and e+ e- annihilation structure functions in QCD,” Zhurnal Éksperimental'noi i Teoreticheskoi Fiziki, vol. 73, p. 1216, 1977.

V. N. Gribov and L. N. Lipatov, “e+ e− pair annihilation and deep inelastic e p scattering in perturbation theory,” Soviet Journal of Nuclear Physics, vol. 15, pp. 675–684, 1972.

V. N. Gribov and L. N. Lipatov, “Annigilyatsiya e+ e− par i gluboko neuprugoe ep-rasseyanie v teorii vozmushchenii,” Yadernaya Fizika, vol. 15, no. 6, pp. 1218–1237, 1972.

V. N. Gribov and L. N. Lipatov, “Deep inelastic e p scattering in perturbation theory,” Soviet Journal of Nuclear Physics, vol. 15, p. 438, 1972.

V. N. Gribov and L. N. Lipatov, “Gluboko neuprugoe ep-rasseyanie v teorii vozmushchenii,” Yadernaya Fizika, vol. 15, no. 4, pp. 781–807, 1972.

G. Altarelli and G. Parisi, “Asymptotic freedom in parton language,” Nuclear Physics B, vol. 126, no. 2, pp. 298–318, 1977.

M. Glück, E. Reya, and A. Vogt, “Dynamical parton distributions of the proton and small-x physics,” Zeitschrift für Physik C, vol. 67, no. 3, pp. 433–447, 1995.

A. Mirjalili and K. Keshavarzian, “The NLO QCD calculation of sea quark distributions in the CORGI approach, based on the constituent quark model,” International Journal of Modern Physics A, vol. 22, no. 24, pp. 4519–4535, 2007.

A. Cafarella, C. Coriano, and M. Guzzi, “NNLO logarithmic expansions and exact solutions of the DGLAP equations from x-space: new algorithms for precision studies at the LHC,” Nuclear Physics B, vol. 748, pp. 253–308, 2006.

F. J. Yndurain, “Reconstruction of the deep inelastic structure functions from their moments,” Physics Letters B, vol. 74, pp. 68–72, 1978.

J. Santiago and F. J. Ynduráin, “Improved calculation of in electroproduction and in neutrino scattering to NNLO and determination of ,” Nuclear Physics B, vol. 611, no. 1–3, pp. 447–466, 2001.

J. J. Aubert, G. Bassompierre, K. H. Becks, et al., “Measurements of the nucleon structure functions in deep inelastic muon scattering from deuterium and comparison with those from hydrogen and iron,” Nuclear Physics B, vol. 293, pp. 740–786, 1987.

D. Allasia, C. Angelini, A. Baldini, et al., “Measurement of the neutron and proton structure functions from neutrino and antineutrino scattering in deuterium,” Physics Letters B, vol. 135, no. 1–3, pp. 231–236, 1984.

A. N. Khorramian, A. Mirjalili, and S. A. Tehrani, “Next-to-leading order approximation of polarized valon and parton distributions,” Journal of High Energy Physics, vol. 2004, no. 10, article 062, 2004.

A. L. Kataev, A. V. Kotikov, G. Parente, and A. V. Sidorov, “Next-to-next-to-leading order QCD analysis of the revised CCFR data for xF ₃ structure function and the higher twist contributions,” Physics Letters B, vol. 417, no. 3-4, pp. 374–384, 1998.

S.-Q. Wang, X.-G. Wu, and S. Brodsky, “Reanalysis of the higher order perturbative QCD corrections to hadronic Z decays using the principle of maximum conformality,” Physical Review D, vol. 90, Article ID 037503, 2014.

W. Bernreuther and W. Wetzel, “Decoupling of heavy quarks in the minimal subtraction scheme,” Nuclear Physics B, vol. 197, no. 2, pp. 228–236, 1982.

W. Bernreuther and W. Wetzel, “Erratum to “decoupling of heavy quarks in the minimal subtraction scheme”, Nuclear Physics B, vol. 197, pp. 228–236, 1982,” Nuclear Physics B, vol. 513, p. 758, 1998.

A. L. Kataev and S. V. Mikhailov, “Generalization of BLM within the -expansion and the principle of maximal conformality,” In press, http://arxiv.org/abs/1408.0122