Rajkhowa, R

Sarma, J K

IACS Institutional Repository

(2004)/<%U PParticular solutions of GLDAP evolution equations in leading order and structurefunctions at low-jrR Rajkhowa and | K Sarma*D e p a r tm e n t o f  P h y sic s , T e/pU r U n iv ersity . N apaam .T c z p u r-7 8 4  02H. A ssam . IndiaE -m aii . jk .s@ teE u .cm ct.inReceived 26 June 200J, accru'd 3 November 2003\l)s trac t W e p re se n t p a r t ic u la r  s o lu tio n s  o f  s in g le t and  n o n -s in g ie t G h b o v -L ip a lo v -D o k sh ilz c i-A lta rc lli  P a iis i (G L D A P ) ev o lu tio n  equation.s m leading order (LO ) a t low-x. W e o b ta in  r-evo lu lions o f  p ro ton  and  neutron  structure functions and r-evo lu lions o f  deu tcron  structure functions at low- \ from GLDAP evo lu tion  equations. T h e  re su lts  o f  /-evo lu tions are com pared  w ith H ERA  Io w -a and \ow Q^  data and those o f  v-cvolutions arc com pared  with NM(' low-A and low-Q^ data.Kc\w(ird.s P a rtic u la r  so lu tio n , c o m p le te  so lu tio n , A lte re lli-P a r is i eq u a tio n , s tru c tu re  func tion  PALS Nos. . 1 2 .3 8 .B x , 12 3 9 .-x , 1 3 .6 0 .H b1. IntroductionI he Gnbov-Lipalov-Dokshitzer-Altarelli-Parisi (GLDAP) evolution equations! 1-4j are fundamental tools to study thef ln ( ( ) ' /  )) and x- evolutions of structure functions, where\ and are Bjorken scaling and four momenta transfer in a deep in e la s tic  scattering (DIS) Process [51 respectively and A is the QCD cut-off parameter. On the other hand, the study of structure functions at 1o w -a  has become topical in view [6 ] of high energy  collider and supercollider experiments (7]. Solutions ot GLDAP evolution equations give quark and gluon structure lu n c iio n s  which produce ultimately proton, neutron and deutcron structure functions. Though numerical solutions are iivailablc m the literature [8], the explorations of the possibility obtaining analytical solutions of GLDAP evolution equations are a lw a y s  interesting. In this connection, some particular solutions computed from general solutions of GLDAP evolution equations at low-jc in leading order have already been obtained applying Taylor expansion method [9] and f-evolutions [ 10] and .r- evolutions [11] o f structure functions have been presented.T'he present paper reports particular solutions of GLDAP ‘^ '^ olution equations computed from complete solutions in leading order at low-jc and calculation of t and x-evolutions for(-orresponding A u th o rsinglet and non-singlet structure functions and hence /-evolutions of proton and neutron structure functions and A-evolulions of deutcron structure functions. In some instance, we can deal with particular solutions more conveniently than with the general solutions [ 12). In calculating structure functions, input data points have been taken from the experimental data directly, unlike the usual practice of using an input distribution function introduced arbitrarily. Results of proton and neutron structure functions are compared with the HERA Iow-a low-Q^ data and those of deutcron structure functions are compared with the NMC 1ow-a Iow-Q^ data. Comparisons are also made with the results of earlier solutions flO, 11, 131 of GLDAP evolution equations. In Section 2, necessary theory has been discus.sed. Section 3 gives results and discussion.2. TheoryThough the basic theory has been discussed elsewhere [ 10, 11,13], the essential steps of the theory have been presented here for clarity. The GLDAP evolution equations for singlet and non- singlet structure functions in the standard forms are f l4 ] :d F ^ jx j )dt- ^ [ { 3  + 4 ln (l-  jc)} Fi{x,t)+ l f ix , t )  + ll(x,t)] = 0 ( 1)©20041ACS368anddFj ' \ x , t) Af  dtwhere.R Rajkhowa and J K SarmaUsing eqs. (8)-(10) in eqs. (3)-{5) and performmo, integrations, we get ^f  (,r,r) = 2 j - ^ {  (l + )F2'^(x / w, i ) - 2 R ^ jc ,r ) |. (3)X ^/ |U / )  = ^  /V^  J{w- + { l - w f ] G ( x / w , l ) d w ,  (4)and=  2 / ^  {(I +  w^) F , ^ \ x / w. 0 - 2 /s '^ ^ U , / ) } . (5)l~ wLet us introduce the variable 1-w and note that [15]X  Xw^  V1— • (6)The series (6) is convergent for |w| < 1. Since jc o v  < I , so 0 < w < l~jc and hence the convergence criterion is satisfied. Now, using Taylor expansion method [9], wc can rewrite G(.vAv, t) as/G(x/w,  r) = G X  +  X  ^ 1 4 ^  J\  A-1 y.  C( .^ 0  + *S »* ^i  . ' (  V  „• V  •S-C(Ar. <)d x ,  +- -  (7)which covers the whole range of m, 0 < « < 1-a-. Since jc is small in our region of discussion, the terms containing x^ and higher powers of jr can be neglected as the first approximation as discussed in our earlier works [11,12,14] and G(x/w, t) can be approximated for small-x asG[xt  w, f )BG(x,r)  + x * dG(x,i)t=i(8)Similarly, F ^ ix /w ,  t) and F^^{x / w, t) can be approximated forsroall-xasF / ( x  / w ,  t )  =  f / ( x , /) +  x ^ M *  ,t=i ^  ^andF^^{x I w , t )B  F^^ (x, t) + x ^  ud  X(9)(10)dFi(x , t )d xn  = Ny ~  (1 ~ jc)(2 ~ + 2 JC" )G( JC,/)+ I - —. (^1 ~ jc) ( 5~4^  + 2jc^ ) + — ln(l / jc ) |  — - (i^ ,1 2  2 J rA "andHere, / = In and Aj = 4/(33-2A^^), Nj being thenumber of flavours and A is the QCD cut off parameter.i“j2jc ln (l/jc)+<?JCNow using eqs. (11) and (12) in eq. (1), we havedFj^x j)dt tA{x)Fiix, t)  + B{x)-^^-^ ' ' -  dx+C(x) G (x ,/)+£X x)dG(x. t)dx= 0 .Let us assume fur simplicityG(x,t) = K{x) F2 (x , t ) ,  where K(x) is a function of x. Now, eq. (14) givesd F h x j )  Afd tL(x) Fi ix ,T)+  M(x)-*^\^'*’'^ - 0 .  lit" a^whereA{x) = 3 + 4 ln ( l '" jc ) - ( l  “'A*)(3 + a), 5 (a) = A(1 -  a'" )  + 2a ln(l/A ),C (x )  =  l / 2 / V ^ ( l - x ) ( 2 - x  +  2x ^ ) ,- i ( l - x ) ( 5 - 4 x  + 2x^) + [ ^ |] ln ( j^ ^ ( x )D{x) = NjX L(x) = /\(x )+  /f(x)C (x)+  D(x)d xand Af (x) = B{x) + K{x) D (x ).Secondly, using eq. (13) in eq. (2), we havedF^HxJ)dFj'HxJ)d t  tP{x)F^Hx, t )  + Q(x)'-d x■0, (17where Differentiating eq. (20) with respect to a , we getParticular solutions ofGLDAP evolution equations etc 369P ( x )  = 3 + 4 In(l -  X) -  (1 -  jc) (jf + 3).anda  = ~ —t expM{x)- d xQ(x) = x { l - x ^ ) ~ 2 x  \nx .The general solutions of eqs. (16) is [9, 12] F (U, V) = 0, where F is an arbitrary function andPutting the value of a  in eq. (20). we obtain the envelopef’>V/ s 1h  (w. r) = —j t~ exp 1 Ux)AjM(x)  Mix)(L (21)antiv[x, t ,F^) = C^ torm a solution of equationdx dt dF^ix, t)AfM(x) - t  -A j  L(x )F2 (x,t) Solving eq.( 18) we obtain.y ( , . , .K ’ ) = ,cx pa n dV ( .v . r . / s ^ > f = ; ; 'U - ,O e x p r j^ a r(18)which is merely a particular solution of the general solution. Unlike the case of ordinary differential equations, the envelope is not a new locus. It is to be noted that when /3 is an arbitrary function of a , then the elimination of a  in eq. (20) is not possible. Thus, the general solution can not be obtained from the complete «olution |9 |, Actually, the general solution of a linear partial differential equation of order one is the totality of envelopes of all one-parameter families (21) obtained from a complete solution.Now, definingr^ S/ \ A 2h  = expUx)A f M { x )  M i x )  ^dxat r = where /q ~ (G j ttt any lower value Q = we get from eq. (21)/s (A ,r)^ f '2 ( A ' o , r ) [ ] , (22)which gives the t-evolution of singlet structure function F/(jr, t) - Proceeding exactly in the same way, and defining= exp II f  U and V are two independent solutions of eq. ( 18) and if a  a nd  13 arc arbitrary constants, then V =  ail  +  P  may be taken as a complete solution of eq. (18). We take this form as this is the simplest form of a complete solution which contains both the arbitrary constants a  and )3. Earlier [10,11,13], we considered a solution A (/ -f = 0, where A and B are arbitrary constants.But that is not a complete solution having both the arbitrary vve get for non-singlet structure function constants as this equation can be transformed to the form V =C(y, where C = -A/B, i e. the equation contains only one arbitrary constant.P(x)AjQ(x) Q{x)jdx(23)Now, the complete solution [9,12]f-'ziwt) exp Ux)Mix)dx = a re x p  . d x  +p  (19) ^ Mix)which gives the /-evolution of non-singlet structure functionF.!^hx,t).Again defining•s a two-parameter family of surfaces, which does not have an envelope, since the arbitrary constants enter linearly [91. i^ifferentiatingeq. (19) with respect to )8, w eget0  = 1 which is iibsurd. Hence, there is no singular solution. The one-parameler l^ amily determined by taking has equationUx)A. Mix) Mix)dxl^iix. /)exp Ux)L' M(x)dx = a t  exp J L f_ L _ d U : + a^ (2 0 )A.  j M(>'ix)we obtain from cq. (21)F2*(x.r)= F/(xo.r)exp| JX,\Ux)AfMix)  Mix)dx (24)372 R Rajkhowa and J K Sarmarespectively, with the HERA low-;c, low-Q^ data [16]. Here, proton structure functions , z) m easured in the range2 < 0 ^  < 50  GeV- , 0.73 < ::<  0.88 and neutron structure functions measured in the range 2<Q^ <50GeV^ , 0.3 < 2 < 0.9 have been used. Moreover, here Pj < 200MeV, where P^ is the transverse momentum of the final state baryon and2= w h e r e ^  are the four momentaof the incident proton and the exchanged vector boson coupling to the positron and p ' is the four-momentum of the final statebaryon. Though we compare our results with y = 2 in ^  relation with data, our results with y maximum, which are equivalent to our earlier results are equally valid. For /-evolutions of deuteron, proton and neutron structure functions, the results will be the range-bounded by our new and old results. But for cvolutions of deuteron structure function, new and old results have not any significance difference.In Figure 1, we present our results of /-evolutions of proton structure functions (solid lines) for the representative values of X given in the figure. Data points at lowest-0“ values in the figure are taken as input to test the evolution equation (33). Agreement is found to be excellent. In the same figure, we alsoFigure 1. /-e v o lu tio n s  o f  p ro to n  s tru c tu re  fu n c tio n s  F {  ( so lid  lin es) fo r th e  re p re se n ta tiv e  v a lu e s  o f  x. D a ta  p o in ts  a t low cst-Q ^  v a lu es  a re  tak en  a s  in p u t to  tes t the  e v o lu tio n  e q u a tio n  (3 3 ). W e a lso  p lo t the  re su lts  o f  /- e v o lu tio n s  o f  p ro to n  s tru c tu re  fu n c tio n s  (d ash ed  lines) fo r o u r  ea r lie r  so lu tio n s  fro m  eq . (3 5 )  o f  G L D A P  ev o lu tio n  eq u a tio n s . F o r co n v en ien ce , v a lu e  o f  e ach  d a ta  p o in t is in c re a se d  by  ad d in g  0 .2 /, w h ere  / = 0 , 1 ,2 ,3 ,... a rc  the  n u m b erin g s  o f  c u rv e s  co u n tin g  fro m  th e  b o tto m  o f  th e  lo w erm o st c u rv e  as th e  0 -lh  o rd e r.plot the results of /-evolutions of proton structure functmn^  ^Ff  (dashed lines) for our earlier solutions from eq. (35)ofGU),\p evolution equations. We observe that our new results are better agreement with data than the old ones.In Figure 2, we present our results of/-evolutions of ncutrun structure functions F? (solid lines) for the representative value* of X given in the figure. Data points at lowest- Q- values m the figure are taken as input to test the evolution eq. (34). Agrecnicn' is found to be excellent. In the same figure, we also phu results of /-evolutions of neutron structure functions (dasfied lines) for our earlier solutions from eq. (36) of GLDAP evolution equations. We observe that in this case also, our new resuKs are in better agreement with data than the old ones.1 8 1 7 1.6 1 5 1 4 1 3 1 2 1 1 1.^109LL0 8  0 7  0 6  0.5 0 4  0 3  0.2 0 1 0ExcO 00033. z^ O 9x*0 00033, z=0 7t i - r i x=.0 00033.x=0 00033, ZsO 0310 15 20Q^ (GeV^ )25 30F ig u r e  2 . /-evo lu tion .s o f  p ro to n  s tru c tu re  fu n c tio n s  F^ (so lid  line".) lor th e  re p re se n ta tiv e  v a lu e s  o f  x. D a ta  p o in ts  a t l o w e s t - v a l u c . s  arc taken as in p u t to  te s t th e  ev o lu tio n  eq u a tio n  (34). W e a lso  p lo t the  re.sults ol / ev o lu tio n s o f  n eu tro n  stru c tu re  fu n c tio n s F2 (d ash ed  lines) for out cailiti so lu tio n s fro m  cq. (3 6 ) o f  G L D A P  e v o lu tio n  eq u a tio n s . F o r convenience value o f  each  d a ta  po in t is in c reased  by  ad d in g  0 .2 /, w here  i = 0* 1 . a rc  the  n u m b erin g s  o f  c u rv e s  co u n tin g  fro m  the bo tto m  o f  the lovvc rrnosi cu rv e  as th e  0 -th  o rder.For a quantitative analysis of ;c-distribution.s of structure functions, we calculate the integrals that occurred in eq. (30) loi N j. -  4. In Figure 3, we present our results of jr-distribution ot deuteron structure functions F2 for K(x) = constant (solid lines). K{x) = (dashed lines) and for AT(.x) = (dotted lines), where a, b, c and d  are constants and for representative values of given in each figure, and compare them with NMC deuteron low-x low- Q- dala[ 17], In each, the data point for x- value ju.st below 0.1 has been taken as input F2 (x^, /).Particular solutions o f GLDAP evolution equations etc 373\{ we take K{x) = 4.5 in cq. (30), then agreement of the result vith experimental data is found to be excellent. On the other h'lnd if vve take K(x) = then agreement of the results with exfierimental data is found to be good at a = 4.5, h = 0.01. Again ,1 we take K(x) = then agreement of the results with experimental data is found to be good at c = 5 .6  = 1.In Figure 4, we present the sensitivity of our results for dillerent constant values of K(x), we observe that at Kix) = 4.5, agreement of the results with experimental data is found to be excellent. If value of Kt\) is increased, the curve goes upward direction and if value of K(x) is decreased, the curve goes downward direction. But the nature of the curve is similar.In Figure 5, we present the sensitivity of our results for different values of a at fixed value of h. Here, we lake h = 0.01. We observe that at a = 4.5, agreement of the results with experimental data is found to be excellent. If value of a is increased, the curve mewes upward and if value of r/ is decreased, curve g(Ks downward. But the nature of the curve is similar.Figure 3. \ d is t i ib u lio n s  o f  d c u le ro n  s t ru c tu re  fu n c tio n s  / ' f  fo r K(x) = CDihtaiit (solid lines), K{.\) = ax^ (d a sh ed  lines) an d  to r  L(x) ~ re  (do tted  liiiLs;, wheic a, b, r  a n d  d  a rc  c o n s ta n ts , a n d  c o m p a re  th e m  w ith  N M C  dcuicion low-x loW ‘ d a ta  In e a c h , th e  d a ta  p o in t  fo r  x- v a lu e  ju s t  klow  0 I has been  taken  as in p u t t) F o r c o n v en ien ce , value o ft.‘di diia point is increa.sed by ad d in g  0 .2 /, w h ere  / =  0 , 1 ,2 ,  3. ... a re  the mimhcnngs o f  cu rv e s  c o u n tin g  fro m  th e  b o tto m  o f  th e  lo w e rm o st cu rv e  as lht‘ O-th Older0 06 0.08ABgure 4. S ensitiv ity  o f  o u r  r e su lts  fo r  d if fe re n t c o n s ta n t v a lu es  o f  K{x).F ig u re  5. S ensitiv ity  o f  o u r resu lts for d iffe ren t values o f  a at fixed  value o f  h = 0 01In Figure 6, we present the sensitivity of our results for different values ofh  at fixed value of a. Here, we take a = 4.5, we observe that at h = 0.01, agreement of the results with experimental data is excellent. If value of b is increased, then the curve gt)es downward and if value of h is decreased, the curve goes upward. But we observe that difference of the curves for h = 0.01, 0.001, 0.0001 is very small and all these curves arc overlapped. Here also the nature of the curves is similar.In Figure 7, we present the sensitivity of our results for different values of V at fixed value of'cf. Here, we lake r/ = 1. We observe that at c = 5, agreement of the results with experimental data is excellent. If value of V is increased, the curve goes upward and if value of V is decreased, the curve goes downward direction. But the nature of the curves is similar.In Figure 8, we present sensitivity of our results for different values of d  at fixed value of c. Here, we take r  = 5, we observe that at = 1, agreement of the results with experimental data is excellent. If value o f 'd' is increased, then the curve goes

Particular Solutions of GLDAP Evolution Equations in Leading Order and Structure Functions at Low-χ

http://arxiv.iacs.res.in:8080/jspui/bitstream/10821/7205/1/Particular%20Solutions%20of%20GLDAP%20Evolution_By%20R%20Rajkhowa.pdf

Particular Solutions of GLDAP Evolution Equations in Leading Order and Structure Functions at Low-χ

Abstract

Similar works

Full text

Available Versions

IACS Institutional Repository