Advances in High Energy Physics

Volume 2013 (2013), Article ID 829803, 10 pages

http://dx.doi.org/10.1155/2013/829803

## Comparison of Analytical Solution of DGLAP Equations for at Small by Two Methods

^{1}Department of Physics, Gauhati University, Guwahati 781014, India^{2}Department of Physics, Physics Academy of North East, Gauhati University, Guwahati 781014, India^{3}Department of Physics, Cotton College, Guwahati 781001, India

Received 26 December 2012; Revised 20 March 2013; Accepted 24 April 2013

Academic Editor: Luca Stanco

Copyright © 2013 Neelakshi N. K. Borah 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.

#### Abstract

The DGLAP equation for the nonsinglet structure function at LO is solved analytically at low by converting it into a partial differential equation in two variables: Bjorken and and then solved by two methods: Lagrange’s auxiliary method and the method of characteristics. The two solutions are then compared with the available data on the structure function. The relative merits of the two solutions are discussed calculating the chi-square with the used data set.

#### 1. Introduction

The structure function of the nucleon has played a pivotal role in our understanding of the internal composition of the proton and the neutron. Experimentally, structure functions are extracted from the cross-section measurement in deep inelastic scattering experiments, and theoretically in QCD, they have a simple interpretation in terms of the quark and gluon momentum distribution, a set of universal functions known as parton distribution functions (PDFs). In DIS kinematics, the PDFs are expressed in terms of two DIS variables: the Bjorken and , the four momentum transferred squared. These functions are not calculable in QCD but their evolution in is predicted by a set of integrodifferential equations known as DGLAP equations [1–4]. The standard and most widely used practice of studying these structure functions is through the numerical solution of these equations [5, 6]. The dependence of the parton structure function at some initial scale , which is small but not so small to break down the perturbative picture, is assumed on various physical grounds, and the parameters are then obtained from a global fit of DIS data. It has been shown by numerous analyses that these equations are in good agreement with the DIS data over a wide kinematic region in and . They can explain data down to very small value of traditionally explained by soft processes and also for very small value of , where these equations are not expected to work. These give credibility to the DGLAP approach as a proper one to study the parton distribution functions. However, apart from the numerical solution, there is the alternative approach of studying analytically these equations at small and there are many analytical solutions available in the literature [7–10], and the present authors have also pursued such an approach with reasonable phenomenological success [11–15]. The analytical approach, though not possible to carry out to higher order in () space due to the complex nature of the splitting functions involved, is intuitive one in the sense that the solutions obtained allow us to visualize their dependence on the variables. In this paper, we study some analytical solutions of the nonsinglet structure functions that is, the flavour dependent contributions to the structure functions considering the corresponding DGLAP evolution equations. As is well known, the nonsinglet structure functions, in DIS plays an important role for precise description of the quark densities; it is comparatively easy because it is not coupled to the singlet and the gluon and can be regarded as a starting ground for the analysis of the other structure functions. We convert the LO DGLAP equation which is an integrodifferential equation into a partial differential equation in the two variables () by a Taylor series expansion valid to be at low . The resulting equation is then solved analytically by two different methods: Lagrange’s auxiliary method [16] and the method of characteristics [17, 18].

In earlier work, the DGLAP equation was solved approximately by using either Lagrange’s method [11–13] or by the method of characteristics [14, 15]. Besides that, the levels of approximation were also different. The aim of the paper is to make a detailed comparison of the predictions of the two methods with two different levels of approximations. We have highlighted above how the present work differs from the earlier works.

The paper is organized as follows: in Section 2 we give the formalism, Section 3 is devoted to discussion of the solutions by comparison with the available data, and in Section 4 we give our conclusion.

#### 2. Formalism

The nonsinglet flavour dependent contribution is defined as where is the density of quark of th flavour and the sum run over all quark flavours. The corresponding nonsinglet structure function is the charged weighed sum of these densities multiplied by . Since the gluon and the singlet quark are flavour independent, the nonsinglet contribution evolves independently in the DGLAP approach and is given by [21] where and are the splitting functions which give the probability of radiating a parton with momentum fraction from a parton with higher momentum fraction . The symbol stands for the convolution integral with respect to the first variable defined as Using the explicit form of the splitting function from [21] the evolution equation for the nonsinglet structure function can be written as Here and are QCD beta function at LO.

To simplify and reduce the integrodifferential equation to a partial differential equation, we introduce a variable and expand the argument as Using (5) we expand in Taylor series as The series (6) is convergent [11–15], and hence at small , we can approximate these by Now putting (7) in (4) we get Carrying out the integration in , we can write (8) as Equation (9) is a partial differential equation for the nonsinglet structure function with respect to the variables and . Beyond its traditional use in evolution, it gives also evolution at small . There are various methods for solving the partial differential equation in two variables. We here adopt the two different methods as mentioned in Section 1.

While performing the integration in and neglecting terms and higher, we can also express (8) as This is what we get by considering during integration.

We solve both the PDEs (9) and (10) with the two formalisms described here, the Lagrange’s method and method of characteristics. Though both these PDEs are obtained from the same equation (4), the levels of approximations are different for both (9) and (10). We continue our discussion for comparing the formalisms with the solutions derived by solving (9) and (10).

##### 2.1. Solution by the Method of Characteristics

Being a differential equation in two variables, it requires two boundary conditions for its solution. But usually we have only one given boundary condition which is the nonperturbative distribution of the function at some initial scale . So, the solutions obtained are not unique, but they give only a range of solutions. This limitation can be avoided by adopting the method of characteristics [17, 18]. From the theory of differential equation, we know that most of the important properties of the solution of (9) depend on the principal part of the equation, that is, the left hand side in (9). But this part is actually a total derivative along the solution of the characteristic equation which gives the characteristic curve of (9). That is, along the characteristic curve, which is a solution of (12), the partial differential equation becomes an ordinary differential equation, and then we can solve it with only one boundary condition.

On using (12) in (9), the left hand side becomes an ordinary derivative with respect to , and the equation becomes an ordinary differential equation: where If the characteristic curve passes through a point , that is, in the space, we get the solution of (12) to be If the characteristic curve cuts the initial curve at a point , then (15) gives so that, Since is any point on the space, therefore dropping the bars over and , the equation of the characteristic is Integrating (13) over from to along the characteristic curve (18), we get the solution for the nonsinglet structure function as Equation (19) is the analytical solution for the nonsinglet structure function of (9) within the present formalism. We can also express it as Using the same formalism for the PDE (10), we get the characteristic equation as Substituting (21) in (10), we get an ordinary differential equation with respect to , where Solving (22) along the characteristic curve, in a similar way we get a different form of solution for the nonsinglet structure function as which can be expressed also as where Equations (19) and (24) are the two analytical solutions of (9) and (10), which are obtained from the same evolution equation (4) with different level of approximations.

##### 2.2. Solution by the Lagrange’s Auxiliary Method

To solve (9) by the Lagrange’s auxiliary method [16], we write the equation in the form where with The general solution of (27) is obtained by solving the following auxiliary system of ordinary differential equations If and are the two independent solutions of (30), then in general, the solution of (27) is where is an arbitrary function of and .

Solving the auxiliary system, we get The functions and are defined as The explicit analytical form of in the leading approximation is [22, 23] In this approach we try to find a specific solution that satisfies some physical conditions on the structure function. Such a solution can be extracted from the combination of and linear in , where and are two quantities to be determined from the boundary conditions on . The two physically plausible boundary conditions are The first one corresponds to a nonperturbative input at some low momentum scale, and the second one corresponds to large behaviour of any structure function at any momentum transfer consistent with constituent quark counting rules. These two boundary conditions leads us to The term does not appear in (42) due to (40). With these boundary conditions, the solution of (27) takes the form From (37) we know Hence, we get This gives the evolution for , nonsinglet structure function in LO at small .

In a similar way when we solve (10) by the Lagrange’s method we write that equation in a form where with Solution of this ODE (46) then leads us to a solution for the nonsinglet structure function as given below Here, as defined by (35) using gives But in this case, Defining as Equation (49) becomes, with for . Here measures deviation of (53) from solution given by (45).

We note that the apparent absence of dependence in the solution (53) is due to algebraic cancellation and boundary condition (40).

Equations (45) and (53) are the solutions of (9) and (10). They are obtained from the same evolution equation (4) as noted earlier. In the next section, we consider their phenomenological utility with respect to each other vis-a-vis the available experimental data. Then we perform a chi-square test to test their compatibility with the data.

#### 3. Results and Discussion

We have obtained two sets of analytical solutions of the same DGLAP evolution equation by applying two different methods, having different levels of approximations. The question we have to address is: which solution is valid in which kinematic region? As well as which set of solution is more compatible with data and why? We test the validity of the solutions by comparing them directly with the available data on the nonsinglet structure function. We also compare our solutions with the MSTW 08 numerical solutions to test the analytical methods. Experiments usually publish the data only for the proton and the neutron structure functions and . So the nonsinglet structure function data is to be extracted from these data using the formula . While deriving the data for the nonsinglet, the statistical and systematic errors of the individual data are added in quadrature, so that maximum possible errors come out. For our analysis we use the data from [19, 20, 24]. To evolve our solutions, we use the MRST2004 [25] input and MSTW2008 [26] input for two different representations of our solutions.

For comparison at first we take the solutions (19) and (45). In Figure 1, we plot the solutions given by (19) and (45) as functions of at some representative fixed , where the data in the said references are given. Though we explore the range of and for CERN-WA25 experiment and and for EMC collaboration data, due to very few experimental data points available, we plot here the data within the range only. In the given range, we see that in the small region the solution by the Lagrange’s method (see (45)) gives a better description of the experimental data than the solution by the method of characteristics (see (19)), whereas in the larger region the opposite is true.

In Figure 2, we plot the solutions (i.e., (19) and (45)) as functions of at different fixed and compare them with the same set of data. Here also we notice the same pattern consistent with the observation in the first graph. The solution by the Lagrange’s method (45) is good for low values of , particularly for the data from the CERN-WA25 [19], whereas the second solution by the method of characteristics (19) represents correctly the data which fall comparatively on the high side, more particularly the data from EMC [20]. Comparing both the graphs we see that this observation, that is, validity of the solution (45) in the low region and that of (19) in the region , is true for all values of explored here.

In Figure 3, we plot the solutions given by (24) and (53) as functions of at some fixed values with the same set of data. In the given range of data we see that both our analytical solutions are no more compatible with the experimental data beyond for the entire range explored here. Similarly in Figure 4, the solutions (24) and (53) are plotted as function of against some fixed values. From Figure 4 also we see that our analytical solutions (24) and (53) do not follow the trend of data beyond for the entire range explored.

In Figures 5 and 6, we plot our set of solutions (19), (45), (24), and (53), respectively, with the BCDMS data where we explore a relatively high range, , and the corresponding range is also high . Here we use the recent MSTW08 input to evolve our solutions and we compare our solutions with the numerical solution by MSTW08 [27]. We compare our solutions as a function of for different values. From the figure we observe that though our both sets of solutions follow the general trend of data and agree with the numerical solutions towards low range, as we approach the high range our solutions overshoot both data and the numerical solutions.

From the above observations we can conclude that for both the set of analytical solutions (19), (45), (24), and (53) obtained from the two PDEs (9) and (10), towards low they converge to the same limit, that is, they predict the same behaviour at low value of . In case of the solutions obtained by Lagrange’s method given by (45) and (53), while (45) shows logarithmic growth with the increasing values, the other solution given by (53) remains almost constant for increasing for fixed values. The solutions by method of characteristics (19) and (24), also show very slow growth with increasing for fixed values.

We note that (10) is less accurate than (9), since the infinite sum in (5) is approximated by only two terms as shown in (11) to derive it and so are the solutions (24) and (53). It is also to be noted that (9) and (10) can be considered equivalent to the accuracy for . Similarly the solutions (19) and (45) are also equivalent to the accuracy for . But for the solutions (24) and (53), due to the neglect of few terms in (10) they do not require to be so. The solutions (19) and (24) appear to be different. The reason of difference is due to the fact that the solution (24) was obtained after neglecting a few terms in (10). Similar is the case for (45) and (53).

For a quantitative estimate of the goodness of fit between the solutions of the two analytical methods with the experimental data, we also do an testing using the formula , where the theoretical and experimental values are for the same th data point with estimated uncertainty . Usually if d.o.f is not much larger than one, the theoretical calculations are considered as being statistically consistent with data. In Table 1. we show that the d.o.f values are consistent with the above discussion for the first set of solutions given by (19) and (45). The numbers of degrees of freedom for the experiments CERN WA25 and EMC are 20 and 32, respectively.

We also compare (24) and (53), the two solutions of (10), obtained by method of characteristics and Lagrange’s method, respectively. d.o.f. values for the two solutions obtained by solving the PDE (10) are shown in Table 2. The analysis also supports our above discussion.

Let us discuss the compatibility or otherwise of the present work with those of López and Ynduráin [28] and Martin [29]. For small , Lopez-Yndurain and Martin-like analysis leads to the following behaviour of nonsinglet structure function: where and is anomalous dimension given as This is to be compared with (19), (24), (45), and (53) above. Using the standard MRST [25] PDF, where . Thus except for the factor , which originates from the input [25], the present work differs from that of López and Ynduráin [28] and Martin [29]. In a sense it is close to the work of Vovk et al. [30], where the effective is dependent contrary to the expectations of López, Ynduráin and Martin.

#### 4. Conclusion

The Taylor approximated DGLAP equation for the nonsinglet structure function, which turns out to be a partial differential equation in two variables, is solved analytically by two different methods: the Lagrange’s auxiliary method and the method of characteristics. However, further approximations in the partial differential equation give us a rather unsatisfactory result. The solutions are valid vis-a-vis the data and numerical solution in two different region of and : the Lagrange’s method solution for the lower and range and the method of characteristics solution for the low and comparatively high range. Considering the solutions together, they are valid in a wide range of as discussed. This demonstrated that two powerful methods of solving differential equations can be applied in the DGLAP framework to obtain analytical solutions. Results of these methods to the polarized structure function have been reported elsewhere [31].

#### References

- V. N. Gribov and L. N. Lipatov, “Deep inelastic ep scattering in perturbation theory,”
*Soviet Journal of Nuclear Physics*, vol. 15, p. 438, 1972. View at Google Scholar - L. N. Lipatov, “The parton model and perturbation theory,”
*Soviet Journal of Nuclear Physics*, vol. 20, p. 94, 1975. 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,”*Journal of Experimental and Theoretical Physics*, vol. 46, pp. 641–653, 1977. View at Google Scholar - 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 - S. Forte, L. Garrido, J. I. Latorre, and A. Piccione, “Neural network parametrization of deep-inelastic structure functions,”
*Journal of High Energy Physics*, vol. 2002, article 062, 2002. View at Publisher · View at Google Scholar - L. del Debbio, S. Forte, J. I. Latorre, A. Piccione, and J. Rojo, “Unbiased determination of the proton structure function ${F}_{2}^{p}$ with faithful uncertainty estimation,”
*Journal of High Energy Physics*, vol. 2005, article 080, 2005. View at Google Scholar - R. D. Ball and S. Forte, “Double asymptotic scaling at HERA,”
*Physics Letters B*, vol. 335, no. 1, pp. 77–86, 1994. View at Publisher · View at Google Scholar - R. D. Ball and S. Forte, “A direct test of perturbative QCD at small
*x*,”*Physics Letters B*, vol. 336, no. 1, pp. 77–79, 1994. View at Google Scholar - A. V. Kotikov and G. Parente, “Small
*x*behavior of parton distributions with soft initial conditions,”*Nuclear Physics B*, vol. 549, no. 1-2, pp. 242–262, 1999. View at Publisher · View at Google Scholar - L. Mankiewick, A. Saalfeld, and T. Weigl, “On the analytical approximation to the GLAP evolution at small
*x*and moderate ${Q}^{2}$,”*Physics Letters B*, vol. 393, no. 1-2, pp. 175–180, 1997. View at Google Scholar - D. K. Choudhury and P. K. Dhar, “An improved analysis of deep inelastic neutrino data with the DGLAP equations at low and high
*x*,”*Indian Journal of Physics*, vol. 81, no. 2, pp. 259–274, 2007. View at Google Scholar · View at Scopus - D. K. Choudhury and P. K. Dhar, “
*x*-distribution of non-singlet structure function at low and high*x*,”*Indian Journal of Physics*, vol. 83, no. 12, pp. 1699–1709, 2009. View at Google Scholar · View at Scopus - J. K. Sarma, D. K. Choudhury, and G. K. Medhi, “x-distribution of deuteron structure function at low-x,”
*Physics Letters B*, vol. 403, no. 1-2, pp. 139–144, 1997. View at Publisher · View at Google Scholar - D. K. Choudhury and P. K. Sahariah, “A solution of the DGLAP equation for gluon at low
*x*,”*Pramana—Journal of Physics*, vol. 58, no. 4, pp. 599–610, 2002. View at Google Scholar · View at Scopus - D. K. Choudhury and P. K. Sahariah, “The next-to-leading order (NLO) gluon distribution from DGLAP equation and the logarithmic derivatives of the proton structure function ${F}_{2}(x,{Q}^{2})$ at low
*x*,”*Pramana—Journal of Physics*, vol. 65, no. 2, pp. 193–213, 2005. View at Publisher · View at Google Scholar - I. N. Sneddon,
*Elements of Partial Differential Equations*, McGraw-Hill, New York, NY, USA, 1957. View at Zentralblatt MATH · View at MathSciNet - W. E. Williams,
*Partial Differential Equations*, Clarendon Press, Oxford, UK, 1980. - E. C. Zachmanoglou and D. W. Thoe,
*Introduction to Partial Differential Equations with Applications*, The Williams & Wilkins Co., Baltimore, Md, USA, 1976. View at Zentralblatt MATH · View at MathSciNet - D. Allasia, C. Angelini, A. Baldini et al., “${Q}^{2}$ Dependence of the Proton and Neutron Structure Functions from Neutrino and Antineutrino Scattering in Deuterium,”
*Zeitschrift für Physik C*, vol. 28, pp. 321–333, 1985. View at Google Scholar - J. J. Aubert, G. Bassompierre, K. H. Becks et al., “Measurements of the nucleon structure functions ${F}_{2}^{N}$ in deep inelastic muon scattering from deuterium and comparison with those from hydrogen and iron,”
*Nuclear Physics B*, vol. 293, pp. 740–786, 1987. View at Google Scholar - L. F. Abbot, W. B. Wood, and R. M. Barnett, “Quantum-chromodynamic analysis of eN deep-inelastic scattering data,”
*Physical Review D*, vol. 22, no. 3, pp. 582–593, 1980. View at Publisher · View at Google Scholar - A. Deshmukhya and D. K. Choudhury, “Non-singlet structure function at high
*x*and neutrino scattering,”*Indian Journal of Physics A*, vol. 75, no. 2, p. 175, 2001. View at Google Scholar - A. Deshmukhya and D. K. Choudhury, “Asymptotics and corrections to asymptotics at low
*x*,”*Indian Journal of Physics A*, vol. 75, no. 5, p. 573, 2001. View at Google Scholar - A. C. Benvenuti, D. Bollini, G. Bruni et al., “A comparison of the structure functions ${F}_{2}$ of the proton and the neutron from deep inelastic muon scattering at high ${Q}^{2}$,”
*Physics Letters B*, vol. 237, no. 3-4, pp. 599–604, 1990. View at Publisher · View at Google Scholar - A. D. Martin, R. G. Roberts, W. S. Stirling, and R. S. Thorne, “Parton distributions incorporating QED contributions,”
*The European Physical Journal C*, vol. 39, no. 2, pp. 155–161, 2005. View at Publisher · View at Google Scholar - A. D. Martin, W. S. Stirling, R. S. Thorne, and G. Watt, “Parton distributions for the LHC,”
*The European Physical Journal C*, vol. 63, no. 2, pp. 189–285, 2009. View at Publisher · View at Google Scholar - The Durham HEP Data, http://durpdg.dur.ac.uk/hepdata/online/f2/structindex.html.
- C. López and F. J. Ynduráin, “Behaviour of deep inelastic structure functions near physical region endpoints from QCD,”
*Nuclear Physics B*, vol. 171, pp. 231–252, 1980. View at Google Scholar · View at Scopus - F. Martin, “Results on nucleon structure functions in quantum chromodynamics,”
*Physical Review D*, vol. 19, no. 5, pp. 1382–1397, 1979. View at Publisher · View at Google Scholar · View at Scopus - V. I. Vovk, A. V. Kotikov, and S. I. Maksimov, “QCD parametrization for structure functions of deep inelastic scattering,”
*Theoretical and Mathematical Physics*, vol. 84, no. 1, pp. 744–751, 1990. View at Publisher · View at Google Scholar · View at Scopus - N. N. K. Borah, D. K. Choudhury, and P. K. Sahariah, “Non-singlet spin structure function ${g}_{1}^{NS}(x,t)$ in the DGLAP approach,”
*Pramana—Journal of Physics*, vol. 79, no. 4, pp. 833–837, 2012. View at Google Scholar