W otrzymanych równaniach ewolucji każdy moment n sprzęga się jedynie z wyższymi momentami (n+p, p > 0), które szybko maleją do 0 wraz z rosnącym p. Stąd, z (5.7) można pozostawić zamknięty układ M + 1 równań dla M + 1 kolejnych obciętych momentów, począwszy od najniższego Nq:
dfNo(xmin,Q2) _ as(Q2) r^(M),
dln Q2 |
2tt |
+c[UNl(xmin)fN°+1(xr, |
™,q2) + . |
df No+1 (xmin, Q2) |
MQ2), |
dlnQ2 |
27r 1 |
+c{MN^\(xmin)fN°+2(xm, |
in ,Q2) + - |
dfN°+M(xmin,Q2) |
<*s(Q2)r |
d\n Q2 |
2tt C |
Rozwiązaniem powyższych równań jest w | |
f'(Xmin,Q2) = (/'(Z. |
j nin,Q20)~ |
xexp \2tt |
1)
fc=i+l
AT0+M
N°+M(Xmin,Q2'.
[«;+!(^n)/yV°+1(^min,Q2)
f Nq+M t,
(5.9)
(5.10)
/ _ «0+M _ \
fi^mimQ2)= /‘(a:rom,Qo)- Aik(Xmin)fk(Xmin,QŹ)]
\ k=i+1 /
<eXP 111 MOIM5)) + ? Aik(Xmin)fk(Xmin,k,Q2), (5.11)
w przypadku biegnącej o,,. Elementy macierzowe D^\xmin) oraz Aij(xmir) podano w dodatku B pracy [H2].
Celem pracy [H2] było otrzymanie równań (5.5)-(5.11), a następnie porównanie rozwiązań (5.10)-(5.11) z wynikami uzyskanymi metodą rozwinięcia na wielomiany Cze-byszewa w przestrzeni zmiennej x. Podejście wielomianowe jest jedną z technik rozwiązywania równań DGLAP, sprowadzając pierwotne różniczkowo - całkowe równanie ewolucji do układu liniowych równań różniczkowych. Wielomiany Czebyszewa [24] były z powodzeniem stosowane przez J. Kwiecińskiego celem dyskretyzacji wielu zagadnień QCD (zob. np. [25]). Szczegółowy opis metody można znaleźć w dodatku A [H2]. Rozwiązując równanie ewolucji z użyciem wielomianów Czebyszewa, otrzymaliśmy w [H2] wyniki dla funkcji rozkładu partonów, które następnie posłużyły do wyliczeń obciętych momentów (3.21). W granicy xmin —» 0, rozwiązania pokrywały się z analitycznymi przewidywaniami (3.18), potwierdzając dobrą dokładność metody, opartej o wielomiany Czebyszewa. Porównanie wyznaczonych tym sposobem obciętych momentów, pozwoliło ocenić dokładność półanalitycznych rozwiązań (5.10)-(5.11) niediagonalnego równania ewolucji. Analizę przeprowadziliśmy w przybliżeniu LO dla szerokiego zakresu wartości punktu obcięcia: 10-5 < Xmin < 0.9 oraz wirtualności Q2: 1 < Q2 < 100 GeV2. Przedmiotem badań były funkcje niesingletowe o ogólnej postaci parametryzacji wejściowej
f(x,Ql) =Afxai(l-x)a2, (5.12)
14