Spin-Dependent Coupled Altarelli-Parisi Equations in the NLO and the Method of Successive Approximations
The coupled Altarelli-Parisi (AP) equations for polarized singlet quark distribution and polarized gluon distribution, when considered in the small limit of the next to leading order (NLO) splitting functions, reduce to a system of two first order linear nonhomogeneous integrodifferential equations. We have applied the method of successive approximations to obtain the solutions of these equations. We have applied the same method to obtain the approximate analytic expressions for spin-dependent quark distribution functions with individual flavour and polarized structure functions for nucleon.
The study of the evolution of the quark and gluon contributions at small (Bjorken variable) towards the spin of the proton through Altarelli-Parisi equations (AP) [1–3] is an important area of DIS. There are no data below and as a result polarized gluon distribution is basically unconstrained at small . There are theoretical arguments that polarized gluon distribution and the unpolarized gluon distribution are connected through the relation at small , but they cannot be verified due to lack of data. Precise measurement of the polarized structure function and its logarithmic scale dependence can determine at small and thus it can reduce the extrapolation uncertainties of in the integral entering in the proton spin sum rule. EIC [4, 5] will allow for a determination of down to a very small value of and it will eventually give the gluon contribution to the spin of the proton over all to about 10 percent accuracy. The same set of measurements will also provide a significantly better determination of the total quark contribution .
The Jacobi polynomial method [5–9] is one of the important methods for obtaining the solutions of spin-dependent Altarelli-Parisi equations. The main advantage of this method is that it allows us to factorize the and dependence of the structure function in a manner that allows an efficient parameterization and evolution of the structure function. The method of successive approximations is used to solve the integral equation . In such a method one begins with a crude approximation to a solution by using an initial condition and improves it step by step by applying a repeatable operation (Picard's method of successive approximations). In this work, we have applied this method for obtaining the solutions of the spin-dependent integrodifferential coupled Altarelli-Parisi equations at small in the NLO and we begin the process by using the boundary condition that the parton distribution vanishes at [11–13]. We have shown that the application of this method in solving these equations results in solutions which appear as summation of series, each term of which is the product of and -dependent functions.
We have structured our work as follows: in Section 2, we have given a method of solution of a system of two first order linear homogeneous differential equations with variable coefficients under certain conditions. In Section 3, we have shown that in the small limit of the splitting functions and under some reasonable approximation of the coupling constant, the AP equations for polarized parton distributions become two first order simultaneous linear nonhomogeneous integrodifferential equations. By using the method described in Section 2 and the method of variation , we have shown that the solutions can be improved through successive approximations. The same procedure is applied to obtain the approximate analytical expressions for polarized quark distributions with individual flavour and using them we have obtained the expressions for polarized structure functions for proton and as well as for neutron . We have compared our solutions with some numerically obtained solutions.
2. Method of Solving a System of FLHE
A system of two first order linear homogenous differential equations (FLHE) with variable coefficients can be given as where , , , and are known coefficients, and are unknown functions to be determined, and is the independent variable.
Equations (2.1) and (2.2) are analytically solvable if the coefficients and , , are constants, that is, independent of . But, as noted in , there is no general method of solving such a system of equations when the coefficients are not constant.
From (2.12) we have where .
Equation (2.13) implies that we can write where is a function of to be determined.
3. Spin-Dependent AP Equations and Polarized Structure Functions in NLO
3.1. Altarelli-Parisi Equations
, the running coupling constant of QCD in NLO, is defined as where and is the number of active flavours.
We define To proceed further and to apply our formalism, we, as in , use the assumption where is a numerical parameter.
3.2. Solutions of AP Equations for and in NLO
We first solve (3.1) and (3.2) for obtaining the approximate analytic expressions for and in NLO. Using the assumption (3.10) and the small splitting functions in NLO, these equations can be written as , , , are some known constants and , , and , are known functions of .
As described in , we obtain the first approximate solutions of (3.11) by replacing and under the integrals appearing in the right-hand side of these equations by their boundary values at [11–13]: With these substitutions (3.11) become Equations (3.13) are two first order simultaneous linear homogeneous differential equations with variable coefficients. We solve these equations by the method described in Section 2 and find the solutions as where and and are constants of integration. From now on we shall represent by and by .
Now applying the input distributions at , we can find out the constants of integration and . With these the solutions after first approximation become where and subscript of and in (3.17) refers to the first approximate solutions.
The solutions of the homogeneous parts of (3.19), that is, the solutions of the first order linear coupled homogeneous equation (3.13) can be obtained by the method described earlier and the solutions are given as (3.14). Now, to obtain the solutions of the nonhomogeneous coupled equation (3.19) we apply the method of variation . Thus the solutions of (3.19) can be given as , , , and their tilde counterparts are known functions of . Equations (3.20) and (3.21) are the second iterative solutions of (3.11) and in comparison to the first iterative solutions (3.17), they are closer to the numerical results as seen from Figures 1 and 2.
We again substitute and from (3.20), and (3.21) respectively, for the and appearing in the integrals in the right-hand side of (3.11) and the resulting equations can be written as , , and and their tilde versions as appearing in (3.22) are known functions of .
Proceeding in a similar manner we can obtain the solutions of (3.22) as The expressions for and , and their tilde versions are calculable functions of . Equations (3.23) are the third iterative solutions of (3.11). Proceeding in this way we can obtain the solutions after successive approximations, which can be written as where , , , and are calculable functions of . Equations (3.24) and (3.25) are our main results.
3.3. Polarized Structure Functions ,
Now to obtain the first approximate solution of (3.28) we replace under the integrals appearing in the right-hand side of (3.28) by its boundary value at [11–13]: With this substitution, (3.28) becomes The solution of (3.30) can be given as where, is the dependent constant of integration, .
Now applying the input distribution at , that is, the solution (3.31) can be written as where the subscript refers to the first approximate solution. Equation (3.33) is the first approximate solution of (3.28).
Again substituting from (3.36) for appearing under the integrals in the right-hand side of (3.28), we have Proceeding in a similar manner we can obtain the solution of (3.37) as where Equation (3.38) is the solution of (3.28) after third approximation. Similarly, the fourth iterative solution will be where
Proceeding in this way, we have the solution of (3.28) after approximation as where Now using the expression for (3.24), the analytical expression for individual quark distributions after approximation can be given as We now use (3.44) and (3.25) (for ) to obtain the analytical expressions for and in NLO as where subscript indicates approximations and and are called Wilson coefficients  given in the small limit as Equations (3.44) and (3.45) are our main results.
3.4. Results and Discussion
Among the several analyses that have included all or most of the present world data on polarized structure functions [18–23] we have used here LSS'05 NLO (MS) input distributions (set-1) at . We have taken and of (3.10) to be .
Initially, as described in Section 4, we have worked out up to third approximation and obtained the solutions after third approximation (3.23) of the approximate Altarelli-Parisi equations at small region (3.11). However, the -dependent parts of the solutions (3.23), namely, , , , and are two-fold integrations of certain hypergeometric and logarithmic functions (the th approximate solutions involve fold such integrations).
To obtain the analytical forms of the fourth approximate solutions, we parametrize the results of the third iteration in the range (the range of where LSS'05 numerical results for parton distributions are available) by the following effective functional forms: where Using (3.47) we get the following approximate analytic forms after fourth approximation:
The expressions for the -dependent functions appearing in (3.49) are calculable.
We observe that with more and more iterations our analytic solutions come closer to the LSS'05 numerical results.
There may be two sources from where some errors have crept in.(1)As the parametrizations were done only in the range (as in LSS'05), it may not be adequate in calculating the integrated quantities like and which involve integrations of our solutions (3.49) in the range . (2)We obtained the solutions of AP equations with small approximation. The disagreement of the integrated quantities with the experimental values, as observed in Table 1, perhaps indicates that incorporation of high effect is important.
In the case of polarized quark distributions with individual flavour also, we have initially worked up to third iterative solutions (3.38). The fourth iterative solution (3.40) contains terms like which can be obtained by evaluating the three integrals as given by (3.41). As such integrals involve several hypergeometric functions along with logarithmic functions, to proceed for approximate fourth iterative analytical solution, we simplify it by performing the parametrizations in the range (the range of where LSS'05 numerical results for parton distributions are available) to get the following effective expressions for
We have obtained the fourth iterative solutions by using these effective expressions for , where as shown by (3.51).We now compare our work for polarized flavour specific quark distributions with the LSS'05 NLO (MS), set-1 numerical results at (Figures 3, 4, and 5). We have observed that increasing iterations bring our solutions for individual quark densities closer to the numerical results at small . In the high range our results however deviate from the numerical results probably due to the use of small splitting functions. Another observation is that while the numerical (LSS'05) result gives negative , our solution for becomes positive beyond .
We now use our formalism in the next to leading order (NLO) to calculate the structure functions and as given by (3.45) and compare them with the LSS'05 numerical results (Figures 6 and 7). It is observed that our approximate analytical results are compatible with that obtained numerically (LSS'05).
QCD analysis of the quark and gluon contributions towards the spin of the nucleon in the small region is very important for a clear understanding of the spin structure of the nucleon and this is mainly done through the Altarelli-Parisi [1–3] evolution equations. In this work we have given a formalism based on the method of successive approximations, for obtaining analytical solutions in the next to leading order (NLO), valid in small region.
In Section 2 we have given a method for solving a system of two first order linear homogeneous differential equations with variable coefficients.
In Section 3 we have obtained approximate analytical solutions of Altarelli-Parisi equations for the polarized singlet quark density and polarized gluon density in the small limit at NLO by using a method described in Section 2 along with the method of iteration. It is observed that, with increasing number of iterations, the solutions approach the numerical results, specifically in the small region. We have also given the analytical expressions for the individual polarized quark densities and using them we have obtained the expressions for the polarized structure functions and . Our results are found to be compatible with those obtained numerically (LSS'05). It is possible that such agreement will improve if the assumption (3.10) can be removed and it will be our attempt in the future communication.
Consider where , and .
G. Altarelli and G. Parisi, “Asymptotic freedom in parton language,” Nuclear Physics B, vol. 126, no. 2, pp. 298–318, 1977.View at: Google Scholar
Yu. 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, p. 641, 1977.View at: Google Scholar
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.View at: Google Scholar
EIC Collaboration, September 2010, http://web.mit.edu/eicc/Documentation.html.
E. Leader, A. V. Sidorov, and D. B. Stamenov, “Polarized parton densities in the nucleon,” Physical Review D, vol. 58, no. 11, Article ID 114028, 1998.View at: Google Scholar
A. V. Sidorov, E. Leader, and D. B. Stamenov, “NLO QCD analysis of polarized deep inelastic scattering,” International Journal of Modern Physics A, vol. 13, pp. 5573–5592, 1998.View at: Google Scholar
V. G. Krivokhizhin, “QCD analysis of structure functions in terms of Jacobi polynomials,” Zeitschrift für Physik C, vol. 36, no. 1, pp. 51–59, 1987.View at: Google Scholar
S. T. Monfared, A. N. Khorramian, F. Arbabifar, and S. A. Tehrani, “The spin dependent parton distribution functions and their moments,” Acta Physica Polonica B, vol. 41, no. 12, pp. 2921–2927, 2010.View at: Google Scholar
G. F. Simmons, Differential Equations, Tata McGraw-Hill, New Delhi, India, 1979.
R. G. Roberts, The Structure of Proton, Cambridge University Press, Cambridge, UK, 1990.
F. J. Yndurain, The Theory of Quark and Gluon Interactions, Springer, Berlin, Germany, 1992.
G. Ridolfi, R. D. Ball, and S. Forte, “Scale dependence and small-χ behaviour of polarized parton distributions,” Nuclear Physics B, vol. 444, no. 1-2, pp. 287–309, 1995.View at: Google Scholar
B. Lampe and E. Reya, “Spin physics and polarized structure functions,” MPI-PhT/98-23 DO-TH 98/02, 1998.View at: Google Scholar
A. Deshamukhya and D. K. Choudhury, “Non-singlet structure function beyond leading order,” in Proceedings of the 2nd Regional Conference of Physics Academy of North East, 2000.View at: Google Scholar
M. Hirai, S. Kumano, and N. Saito, “Determination of polarized parton distribution functions and their uncertainties,” Physical Review D, vol. 69, no. 5, Article ID 054021, 10 pages, 2004.View at: Google Scholar
J. Blumlein and H. Bottcher, “QCD analysis of polarized deep-inelastic scattering data and parton distributions,” Nuclear Physics B, vol. 636, no. 1-2, pp. 225–263, 2002.View at: Google Scholar
A. V. Sidorov, E. Leader, and D. B. Stamenov, “Role of positivity constraints in determining polarized parton densities,” Journal of High Energy Physics, vol. 6, 2005.View at: Google Scholar
E. Leader, A. V. Sidorov, and D. B. Stamenov, “Impact of CLAS and COMPASS data on polarized parton densities and higher twist,” Physical Review D, vol. 75, no. 7, Article ID 074027, 10 pages, 2007.View at: Google Scholar
V. Y. Alexakhin, Yu. Alexandrov, G. D. Alexeev et al., “Spin asymmetry and the spin-dependent structure function of the deuteron at low values of x and ,” Physics Letters B, vol. 647, no. 5-6, pp. 330–340, 2007.View at: Google Scholar
A. Airapetian and The HERMES Collaboration, “Single-spin asymmetries in semi-inclusive deep-inelastic scattering on a transversely polarized hydrogen target,” Physical Review Letters, vol. 94, no. 1, Article ID 012002, 6 pages, 2005.View at: Google Scholar
A. Airapetian and The HERMES Collaboration, “Precise determination of the spin structure function g1 of the proton, deuteron, and neutron,” Physical Review D, vol. 75, no. 1, Article ID 012007, 48 pages, 2007.View at: Google Scholar
A. Airapetian and The HERMES Collaboration, “Quark helicity distributions in the nucleon for up, down, and strange quarks from semi-inclusive deep-inelastic scattering,” Physical Review D, vol. 71, no. 1, Article ID 012003, 36 pages, 2005.View at: Google Scholar
M. Alekseev, V. Yu. Alexakhin, Yu. Alexandrov et al., “The polarised valence quark distribution from semi-inclusive DIS,” Physics Letters B, vol. 660, no. 5-6, pp. 458–465, 2008.View at: Google Scholar