arxiv:hep-ph/ v2 13 Mar 2001

Transcription

1 TTP00-13 NIKHEF Some higher moments o deep inelastic structure unctions at next-to-next-to leading order o perturbative QCD arxiv:hep-ph/ v2 13 Mar 2001 A. Rétey a and J.A.M. Vermaseren b a Institut ür Theoretische Teilchenphysik, Universität Karlsruhe, D Karlsruhe, Germany b NIKHEF Theory Group Kruislaan 409, 1098 SJ Amsterdam, The Netherlands Abstract We present the analytic next-to-next-to-leading QCD calculation o some higher moments o deep inelastic structure unctions in the leading twist approximation. We give results or the moments N=1,3,5,7,9,11,13 o the structure unction F 3. Similarly we present the moments N=10,12 or the lavour singlet and N=12,14 or the non-singlet structure unctions F 2 and F L. We have calculated both the three-loop anomalous dimensions o the corresponding operators and the three-loop coeicient unctions o the moments o these structure unctions. 1

2 1 Introduction The determination o the next-to-next-to-leading NNL) order QCD approximation or the structure unctions o deep inelastic scattering has become important or the understanding o perturbative QCD and necessary or an accurate comparison o perturbative QCD with the increasing precision o experiments. Such calculations however are rather complicated and hence a complete NNL result does not exist as o yet. The one-loop anomalous dimensions were calculated in 1, 2. In 3 see also the reerences therein) the complete one-loop coeicient unctions were obtained. Anomalous dimensions at 2-loop order were obtained in 4, 5, 6, 7, 8, 9, 10, 11, 12. The 2-loop coeicient unctions were calculated in 13, 14, 15, 16, 17, 18, 19, 20, 21. Analytical results o the 3-loop anomalous dimensions and coeicient unctions o the moments o F 2 and F L are only known or the moments N=2,4,6,8 in both the singlet and non-singlet case and additionally or N =10 in the non-singlet case rom 22 and 23. In addition the Gross- Llevellyn Smith sum rule, which corresponds to the irst moment o F νp+ νp 3 has been calculated at this order 24. For a complete reconstruction o the x-dependence o the structure unctions via an inverse Mellin-transormation one would need the moments or all N that is either all even or all odd integer values). Additionally one needs them both or F 2 and F 3 in order to untangle the various quark and gluon contributions. The determination o the NNL approximation or generic N is work in progress 25, but probably will not be inished in the near uture. The available moments o F 2 have been used by a number o authors to make a reconstruction o the complete structure unctions at NNL by a variety o means 23, 26, 27, 28. Additionally they can be used to obtain a better value o α S 29 It should be clear that it is important to have as large a number o moments as possible. First, these results can be immediately used to increase the precision o phenomenological investigations o deep inelastic scattering 30. Second, the moments will be a very important check or the new methods and programs needed or the determination o the 3-loop results or arbitrary N. Unortunately it is not very easy to increase the number o moments, because each new moment requires roughly ive times the computer resources that its predecessor needs. With the advent o better computers this means that by now it has been possible to obtain two more moments or the singlet case and one additional moment or the non singlet F 2 case. This should allow or instance a somewhat better determination o α S. More important however is the determination o the irst seven odd moments o F νp+ νp 3 to three loops. To this end we used the the same programs as in 23, state o the art computers and a new version o the symbolic manipulation program FORM 31 that supports now 64-bit architectures and to some extend parallel computers see also 32). We could push the limit in these calculations to include two new moments N =10,12) in the calculation o the lavour singlet structure unctions F L and F 2, and two new moments N = 12,14) or the lavour nonsinglet structure unctions F L and F 2. Additionally we have computed the moments N=3,5,7,9,11,13 o the structure unction F 3. We do not expect more moments to become available beore the complete results or all N will be presented. 2 The ormalism This calculation ollows the one presented in 23 see also 21) in every detail, so we only will give a very short review on the methods used. We need to calculate the hadronic part o the amplitude or unpolarized deep inelastic scattering 2

3 which is given by the hadronic tensor W µν x, Q 2 ) = 1 d 4 ze iq z p, nucl J µ z)j ν 0) nucl, p 4π = g µν q µ q ν ) 1 2x F Lx, Q 2 ) 4x + g 2 µν p µ p ν q 2 p µq ν + p ν q µ ) 2x ) 1 q 2 2x F 2x, Q 2 ) +iǫ µνρσ p ρ q σ p q F 3x, Q 2 ) 1) where the J µ are either electromagnetic or weak hadronic currents and x = Q 2 /2 p q) is the Bjorken scaling variable with 0 < x 1. Q 2 = q 2 is the transered momentum and p, nucl is the nucleon state with momentum p. In these equations spin averaging is assumed. The longitudinal structure unction F L is related to the structure unction F 1 by F L = F 2 2xF 1. For electronnucleon scattering J µ is the electromagnetic quark current and F 3 vanishes. For neutrino-nucleon scattering J µ is an electroweak quark current which has an axial vector contribution and F 3, which describes parity violating eects that arise rom vector and axial-vector intererence, will not vanish. Using the dispersion relation technique one can relate the hadronic tensor to the ollowing 4-point Green unctions: W µν p, q) = 1 2π ImT µν p, q), T µν p, q) = i d 4 ze iq z p, nucl T J µ z)j ν 0) nucl, p. Applying a ormal operator product expansion in terms o local operators to the time-ordered product o two quark currents leads to: i d 4 ze iq z T J ν1 z)j ν2 0) = ) N 1 Q 2 g ν1ν 2 q ) ν 1 q ν2 q 2 q µ1 q µ2 C jl,n Q2 µ 2, a s) N,j g ν1µ 1 g µ2ν 2 q 2 ) g ν1µ 1 q ν2 q µ2 g ν2µ 2 q ν1 q µ1 + g ν1ν 2 q µ1 q µ2 C j 2,N Q2 µ 2, a s) +i ǫ ν1ν 2µ 1ν 3 g ν3ν4 q ν4 q µ2 C j 3,N Q2 µ 2, a s) q µ3...q µn O j,{µ1,...µn } 0) + higher twists, j = α, ψ, G 2) Here we have introduced the notation a s = α s /4π) = g 2 /4π) 2 and everything is assumed to be renormalized. The sum over N runs over the standard set o the spin-n twist-2 irreducible lavour non-singlet quark operators and the singlet quark and gluon operators: O α,{µ1,...,µn } = ψλ α γ {µ1 D µ2 D µn } ψ, α = 1, 2,...,n 2 1) O Ψ,{µ1,...,µN } = ψγ {µ1 D µ2 D µn } ψ O G,{µ1,...,µN } = G {µµ1 D µ2 D µn 1 G µnµ} Application o this OPE to Eq. 1) leads to an expansion or the unphysical values x. From the proper analytical continuation to the physical region 0 < x 1 one inds or the moments o the structure unctions F 2, F L and F 3 : M k,n 2 = M 3,N 1 = dxx N 2 F k x, Q 2 ) = dxx N 1 F 3 x, Q 2 ) = i=α i=α,ψ,g C i k,n Q2 µ 2, a s)a i nucl,n, k = 2, L 3) C i,n Q2 µ 2, a s)a i nucl,n 4) 3

4 with the spin averaged matrix elements p, nucl O j,{µ1 µn } p, nucl = p {µ1 p µn } A j nucl,n p2 µ 2 ) 5) In the derivation o 3) one needs the symmetry properties o T µν under x x. This is why one can only ind either even or odd moments rom these equations, dependent on the process under consideration. For F 3 we will only consider the lavor non-singlet contributions, due to the properties o the operators O ψ and O G under charge conjugation there should not be a singlet contribution see e.g. 33). The scale-dependence o the coeicient unctions is then covered by the renormalization group equations: k=ψ,g µ 2 µ 2 + βa s) { µ 2 µ 2 + βa s) a s γn ns a s } δ jk γ jk N Q Ci,N ns 2 Q 2 C k i,n ) µ 2, a s = 0, i = 2, 3, L 6) ) µ 2, a s = 0, i = 2, L; j = ψ, G 7) The non-singlet coeicient unctions and anomalous dimensions don t depend on the index α, and we have adopted the conventional collective denotation ns or them. 3 The even moments o F 2 and F L Equation 2) is a relation between operators and does not depend on the hadronic states to which the OPE is applied. Using the method o projectors 34 one can ind both the coeicient unctions and the anomalous dimensions or the even moments o F 2 and F L as deined in Eq. 6) rom the ollowing 4-point Green unctions: Tµν qγqγ = i d 4 ze iq z p, quark T J µ z)j ν 0) quark, p 8) Tµν gγgγ = i d 4 ze iq z p, gluon T J µ z)j ν 0) gluon, p 9) Applying to Eqs. 8) the projectors q {µ1 q µn } N P N N! p µ1 p µn and, to project out the dierent Lorentz projections these projectors are valid in D = 4 2ǫ and or the leading twist approximation): p=0 10) q2 P L = p q) 2 pµ p ν 3 2ǫ q 2 P 2 = 2 2ǫ p q) 2 pµ p ν + 1 ) 2 2ǫ gµν, as well as the corresponding lavour projections results in the equations: T qγqγ,s k,n Q2 µ 2, a s, ǫ) = C ψk,n Q2 µ 2, a s, ǫ)z ψψ N a s, 1 ǫ ) + CG k,n Q2 µ 2, a s, ǫ)z Gψ N a s, 1 ) ǫ ) T gγgγ,s k,n Q2 µ 2, a s, ǫ) = C ψk,n Q2 µ 2, a s, ǫ)z ψg N a s, 1 ǫ ) + CG k,n Q2 µ 2, a s, ǫ)zn GG a s, 1 ) ǫ ) T qγqγ,ns k,n Q2 µ 2, a s, ǫ) = Ck,N ns Q2 µ 2, a s, ǫ)zn ns a s, 1 ǫ )Ans,tree quark,n ǫ) A ψ,tree quark,n ǫ) A G,tree gluon,n ǫ) k = 2, L 11) 4

5 From these equations the coeicient unctions and rom the Z ij N the anomalous dimensions can be calculated in the usual way. It should be mentioned that in the Eqs. 11) on the let hand side ater applying the projectors 10) we are let with only diagrams o the massless propagator type, a problem solved at 3-loop order long ago 35 and implemented in an eicient way in the FORM package MINCER 36. On the right hand side only the tree level diagrams contributing to the Matrix elements survive. It turns out that much computing time can be saved when calculating additionally Green unctions with external ghosts to get rid o the unphysical polarization states o the external gluons instead o using the very complicated projection onto physical states. Also, rom Eqs. 11) one can determine the ZN GG and ZGψ N only to order α2 s. To obtain the α3 s-contributions one can calculate additionally Greens unctions with external scalar ields φ that couple to gluons only at tree level. Altogether, to obtain the coeicient unctions and anomalous dimensions or the even moments with N=10,12 o F 2 and F L the ollowing diagrams had to be calculated q=quark, g=gluon, γ=photon, h=ghost, φ=scalar ield): tree 1-loop 2-loops 3-loops Lorentz projections qγqγ qφqφ gγgγ hγhγ gφgφ hφhφ Total The q {µ1 q µn } in Eq. 10) are the harmonic i.e. symmetrical and traceless) part o the tensor q µ1 q µn. The number o terms in these harmonic tensors explodes as N increases and this is the real limitation in these calculations considering the computing time as well as disk-space usage. In spite o a very eicient implementation o these tensors see 37) or N = 12, singlet and N = 14, nonsinglet, individual diagrams had a disk space usage up to and over 100 GB. Altogether the calculation o all the above diagrams or N = 10, 12 took approximately 5 weeks on a Compaq Server with 8 Alpha processors running at 700 MHz, 4 GB o RAM and GB o disk-space. The N = 14 nonsinglet calculation took comparable resources. 4 The odd moments o F 3 The coeicient unctions and anomalous dimensions o the odd moments o the structure unction F 3 can be obtained in the same way as the non-singlet part o F 2 and F L but now considering the time-ordered product o one vector current V µ and one axial vector current A ν. The axial current introduces the appearance o a γ 5 and some care has to be taken to treat it correctly within the ramework o dimensional regularization. We adopt the deinition used in 24: γ µ γ 5 = i 1 6 ǫ µνρσγ µ γ ν γ ρ γ σ Projecting out the lavour non-singlet part and the corresponding Lorentz structure with: 1 P 3 = i p αq β 1 2ǫ)2 2ǫ) ǫµναβ p q one inds products o metric tensors which have to be considered as D-dimensional objects. Since this deinition o γ 5 in D dimensions violates the axial Ward identity one needs to renormalize A µ with a renormalization constant Z A and additionally apply a inite renormalization with Z 5, both o these constants are given to 3-loop order in 24. Combining all this inally leads to Z A a s, 1 ǫ )Z 5a s, ǫ)t ns 3,N Q2 µ 2, a s, ǫ) = C 3,N Q2 µ 2, a s, ǫ)z 1 N a s, 1 ǫ )Ans N,tree ǫ) 5

6 Due to the γ 5 insertion at one o the vertices, some o the symmetries that were used to minimize the number o diagrams could not be applied in this case and to determine the T3,N ns 1076 = ) diagrams had to be evaluated, which took about 6 weeks or the moments N =1,3,5,7,9,11,13. 5 Results Using the strategies sketched in the previous section we ind the ollowing results or the coeicient unctions and anomalous dimensions. Again ollowing Re. 23, we present the combined singlet and non-singlet results or F 2 and F L in terms o lavour actors which are deined in the ollowing table or n number o lavours: l 2 l 11 l 02 l g 2 l g 11 non-singlet 1 3 n n =1 e singlet 1 1 =1 e ) 2 n n n =1 e =1 e ) 2 n n n =1 e2 The numerical values o the anomalous dimensions or the spin-even operators contributing to F 2 and F L that now are known to NNL order are: γ ψψ 2 = a s + a 2 s n l 02 n ) s n n 2 +l n n 2 )) γ ψψ 4 = a s + a 2 s n l 02 n ) s n n 2 +l n n 2 )) γ ψψ 6 = a s + a 2 s n l 02 n ) s n n 2 +l n n 2 )) γ ψψ 8 = a s + a 2 s n l 02 n ) s n n 2 +l n n 2 )) γ ψψ 10 = a s + a 2 s n l 02 n ) s n n 2 +l n n 2 )) γ ψψ γ ψψ,ns 12 = a s + a 2 s n l 02 n ) s n n 2 +l n n 2 )) 14 = a s + a 2 s n ) s n n 2 ) 6

7 γ ψg 2 = a s n a 2 s n s n n 2 ) γ ψg 4 = a s n a 2 s n s n n 2 ) γ ψg 6 = a s n a 2 s n s n n 2 ) γ ψg 8 = a s n a 2 s n s n n 2 ) γ ψg 10 = a sn a 2 s n s n n 2 ) γ ψg 12 = a sn a 2 s n s n n 2 ) γ Gψ 2 = a s + a 2 s n ) s n n 2 ) γ Gψ 4 = a s + a 2 s n ) s n n 2 ) γ Gψ 6 = a s + a 2 s n ) s n n 2 ) γ Gψ 8 = a s + a 2 s n ) s n n 2 ) γ Gψ 10 = a s + a 2 s n ) s n n 2 ) γ Gψ 12 = a s + a 2 s n ) s n n 2 ) γ GG 2 = a s n a 2 s n s n n 2 ) γ GG 4 = a 2 s n ) + a s n ) s n n 2 ) γ GG 6 = a 2 s n ) + a s n ) s n n 2 ) γ GG 8 = a 2 s n ) + a s n ) s n n 2 ) γ GG 10 = a 2 s n ) + a s n ) s n n 2 ) γ GG 12 = a 2 s n ) + a s n ) s n n 2 ) The corresponding coeicient unctions read: 7

8 C ψ 2,2 = a s + a 2 s n l 02 n ) s n l 11 n n 2 +l n n 2 )) C ψ 2,4 = a s + a 2 s n l 02 n ) s n l 11 n n 2 +l n n 2 )) C ψ 2,6 = a s + a 2 s n l 02 n ) s n l 11 n n 2 +l n n 2 )) C ψ 2,8 = a s + a 2 s n l 02 n ) s n l 11 n n 2 +l n n 2 )) C ψ 2,10 = a s + a 2 s n l 02 n ) s n l 11 n n 2 +l n n 2 )) C ψ 2,12 = a s + a 2 s n l 02 n ) s n l 11 n n 2 +l n n 2 )) C ψ,ns 2,14 = a s + a 2 s n ) s n l 11 n n 2 ) C2,2 G = 0.5 a s n a 2 s n s n n l g 11 ) n2 C2,4 G = a s n a 2 s n s n n l g 11 ) n2 C2,6 G = a s n a 2 s n s n n lg 11 ) n2 C2,8 G = a sn a 2 s n s n n lg 11 ) n2 C2,10 G = a sn a 2 s n s n n l g 11 ) n2 C2,12 G = a s n a 2 s n s n n l g 11 ) n2 8

9 C ψ L,2 = a s + a 2 s n l 02 n ) s n l 11 n n 2 +l n n 2 )) C ψ L,4 = a s + a 2 s n l 02 n ) s n l 11 n n 2 +l n n 2 )) C ψ L,6 = a s + a 2 s n l 02 n ) s n l 11 n n 2 +l n n 2 )) C ψ L,8 = a s + a 2 s n l 02 n ) s n l 11 n n 2 +l n n 2 )) C ψ L,10 = a s + a 2 s n l 02 n ) s n l 11 n n 2 +l n n 2 )) C ψ L,12 = a s + a 2 s n l 02 n ) s n l 11 n n 2 +l n n 2 )) C ψ,ns L,14 = a s + a 2 s n ) s n l 11 n n 2 ) CL,2 G = a sn a 2 s n s n n l g 11 ) n2 CL,4 G = a s n a 2 s n s n n l g 11 ) n2 CL,6 G = a s n a 2 s n s n n lg 11 ) n2 CL,8 G = a s n a 2 s n s n n lg 11 ) n2 CL,10 G = a sn a 2 s n s n n lg 11 ) n2 CL,12 G = a sn a 2 s n s n n l g 11 ) n2 The numerical values o the anomalous dimensions or the odd moments o F 3 read: 9

10 γ1 ns = 0 γ3 ns = a s + a 2 s n ) s n n 2 ) γ ns 5 = a s + a 2 s n ) s n n 2 ) γ ns 7 = a s + a 2 s n ) s n n 2 ) γ ns 9 = a s + a 2 s n ) s n n 2 ) γ ns 11 = a s + a 2 s n ) s n n 2 ) γ ns 13 = a s + a 2 s n ) s n n 2 ) and the coeicient unctions are: C ns 3,1 = 1 4 a s + a 2 s n ) s n n 2 ) C ns 3,3 = a s + a 2 s n ) s n n 2 ) C ns 3,5 = a s + a 2 s n ) s n n 2 ) C ns 3,7 = a s + a 2 s n ) s n n 2 ) C ns C ns 3,9 = a s + a 2 s n ) s n n 2 ) 3,11 = a s + a 2 s n ) s n n 2 ) C ns 3,13 = a s + a 2 s n ) s n n 2 ) 6 Acknowledgements A.R. would like to thank Sven Moch, Timo van Ritbergen and Thomas Gehrmann or many useul discussions and Denny Fliegner or technical support. J.V. would like to thank the University o Karlsruhe and the DFG or repeated hospitality. This work was supported by the DFG under Contract Ku 502/8-1 DFG-Forschergruppe Quanteneldtheorie, Computeralgebra und Monte-Carlo-Simulationen ) and the Graduiertenkolleg Elementarteilchenphysik an Beschleunigern at the Univerität Karlsruhe. 10

11 Reerences 1 D. J. Gross and F. Wilczek, Phys. Rev. D8, ). 2 D. J. Gross and F. Wilczek, Phys. Rev. D9, ). 3 W. A. Bardeen, A. J. Buras, D. W. Duke, and T. Muta, Phys. Rev. D18, ). 4 E. G. Floratos, D. A. Ross, and C. T. Sachrajda, Nucl. Phys. B129, ). 5 Nucl. Phys. B139, ). 6 E. G. Floratos, D. A. Ross, and C. T. Sachrajda, Nucl. Phys. B152, ). 7 A. Gonzalez-Arroyo, C. Lopez, and F. J. Yndurain, Nucl. Phys. B153, ). 8 A. Gonzalez-Arroyo and C. Lopez, Nucl. Phys. B166, ). 9 C. Lopez and F. J. Yndurain, Nucl. Phys. B183, ). 10 W. Furmanski and R. Petronzio, Phys. Lett. B97, ). 11 G. Curci, W. Furmanski, and R. Petronzio, Nucl. Phys. B175, ). 12 R. Hamberg and W. L. van Neerven, Nucl. Phys. B379, ). 13 D. W. Duke, J. D. Kimel, and G. A. Sowell, Phys. Rev. D25, ). 14 A. Devoto, D. W. Duke, J. D. Kimel, and G. A. Sowell, Phys. Rev. D30, ). 15 D. I. Kazakov and A. V. Kotikov, Nucl. Phys. B307, ). 16 D. I. Kazakov and A. V. Kotikov, Phys. Lett. B291, ). 17 D. I. Kazakov, A. V. Kotikov, G. Parente, O. A. Sampayo, and J. S. Guillen, Santiago de Compostela Univ. - US-FT ,rec.Jul.) 6 p. 18 W. L. van Neerven and E. B. Zijlstra, Phys. Lett. B272, ). 19 E. B. Zijlstra and W. L. van Neerven, Phys. Lett. B273, ). 20 E. B. Zijlstra and W. L. van Neerven, Nucl. Phys. B383, ). 21 S. Moch and J. A. M. Vermaseren, Nucl. Phys. B573, ), hep-ph/ S. A. Larin, T. van Ritbergen, and J. A. M. Vermaseren, Nucl. Phys. B427, ). 23 S. A. Larin, P. Nogueira, T. van Ritbergen, and J. A. M. Vermaseren, Nucl. Phys. B492, ), hep-ph/ S. A. Larin and J. A. M. Vermaseren, Phys. Lett. B259, ). 25 J. A. M. Vermaseren and S. Moch, 2000), hep-ph/ G. Parente, A. V. Kotikov, and V. G. Krivokhizhin, Phys. Lett. B333, ), hepph/ A. Kataev, A. V.Kotikov, G. Parente and A. V. Sidorov, Phys. Lett. B388, ), hep-ph/ W. L. van Neerven and A. Vogt, Nucl. Phys. B568, ), hep-ph/ J. Santiago and F. J. Yndurain, Nucl. Phys. B563, ). 11

12 30 W. L. van Neerven and A. Vogt 2000), hep-ph/ J. A. M. Vermaseren, Symbolic Manipulation with Form version 2, Tutorial and Reerence Manual Computer Algebra Nederland, Amsterdam, 1991). 32 D. Fliegner, A. Retey, and J. A. M. Vermaseren, 1999), hep-ph/ A. J. Buras, Rev. Mod. Phys. 52, ). 34 S. G. Gorishnii and S. A. Larin, Nucl. Phys. B283, ). 35 K. G. Chetyrkin and F. V. Tkachev, Nucl. Phys. B192, ). 36 S. A. Larin, F. V. Tkachev, and J. A. M. Vermaseren, NIKHEF-H S. A. Larin, T. van Ritbergen, and J. A. M. Vermaseren, Prepared or 3rd International Workshop on Sotware Engineering, Artiicial Intelligence and Expert systems or High-energy and Nuclear Physics, Oberammergau, Germany, 4-8 Oct J. A. Gracey, Phys. Lett. B322, ), hep-ph/ A Conventions Here we give the complete expressions or the newly computed moments and coeicient unctions. The notation o the color actors is as usual: The Casimir operators o the undamental and adjoint representation are denoted by C F and C A and their values or the color group SU3) are 4 3 and 3, respectively. Additionally we are using the symmetric structure constants o SU3) or which d abc d abc = For the trace normalization o the undamental representation we have inserted T F = 1 2. For a generic SUn) group the number o generators equals N A = n 2 1. The values o the Riemann ζ unction are written as ζ n = ζn). The only check or the newly computed 3-loop results is rom a large n -expansion o 38 where the n 2 terms are calculated. These terms are in agreement with ours. To test the new parallel version o FORM we have also recalculated the lower moments or F 2 and F L and ound complete agreement with 23 except or one missing term. In C G 2,2 there should be the extra term a3 s n C 2 F ζ 5). The numerical values given in this reerence are correct. It should be noted that the terms in d abc d abc enter or the irst time at the three loop level and help in the determination o P S qq P S qq. B Results or the moments o F 2 and F L γ10 GG = a s C A a s n a 2 sca a 2sC F n a 3 s C3 A a 3 sc F C A n a 3 s C2 F n ζ ζ 3 + a 3 sc 2 An a 3 s C Fn ζ a 3 s C An a 2sC A n

13 γ12 GG = a s C A a s n a 2 s C An a 3 s C3 A a 2 3 sca 2 + a 2s C F n a 3 sc F C A n ζ ζ 3 + a 3 s C2 F n + a 3 scan ζ 3 + a 3 s C F n a 3 sc A n γ Gψ 10 = a s C F a 2 sc F C A a 3 s C FC 2 A + a 2 sc 2 F ζ 3 + a 3 sc 2 FC A a 3 s C3 F ζ ζ 3 + a 3 sc F C A n a 3 s C2 F n ζ ζ 3 + a 3 s C F n 2 + a 2sC F n

14 γ Gψ 12 = a s C F 79 + a 2sC F C A a 2 s C Fn a 3 s C FC 2 A a 2 sc 2 F ζ 3 + a 3 sc 2 FC A a 3 s C3 F ζ ζ 3 + a 3 sc F C A n ζ 3 + a 3 s C2 F n ζ 3 + a 3 sc F n γ ψg 10 = a s n 28 + a 2s 165 C F n a 2s C An + a 3 sc F C A n ζ 3 + a 3 s C2 F n a 3 sc 2 An a 3 s C F n ζ ζ a 3 s C An γ ψg 12 = a s n a 2 s C Fn a 2s C An a 3 sc F C A n ζ ζ 3 + a 3 s C2 F n + a 3 scan ζ 3 + a 3 s C Fn a 3 sc A n

15 γ ψψ 10 = a s C F a 2 s C FC A a 2 s C2 F a 2 s C Fn a 2s l 02C F n a 3 sc F CA ζ 3 + a 3 s C2 F C A ζ 3 + a 3 scf ζ 3 + a 3 s C FC A n ζ 3 + a 3 scfn ζ 3 + a 3 s C Fn a 3 sl 02 C F C A n ζ 3 + a 3 s l 02CF 2 n ζ 3 + a 3 sl 02 C F n

16 γ ψψ 12 = a s C F a 2 s C F C A a 2 s C2 F a 2 s C F n a 2s l 02C F n a 3 sc F C 2 A + a 3 s C2 F C A + a 3 sc 3 F ζ ζ ζ 3 + a 3 s C F C A n ζ ζ a 3 sc 2 F n a 3 s C F n 2 + a 3 sl 02 C F C A n ζ 3 + a 3 s l 02CF 2 n ζ 3 + a 3 sl 02 C F n γ14 ns = a s C F a 2 sc F C A a 2 s C F n a 2 scf a 3 s C F C A n ζ 3 + a 3 sc F C 2 A + a 3 s C F n ζ 3 + a 3 scf 2 C A ζ 3 + a 3 scf 2 n ζ 3 + a 3 s C3 F ζ 3 16

17 C ψ 2,10 = 1 + a s C F a 2sC F C A a 2 s C F n a 2 sl 02 C F n a 3 s C F C A n a 3 sc F C 2 A ζ ζ 5 + a 3 sc F n ζ 3 + a 3 s C2 F C A ζ 3 + a 2 s C2 F ζ ζ ζ ζ ζ ζ ζ 5 + a 3 s C2 F n ζ ζ 4 + a 3 scf ζ ζ ζ 5 + a 3 s l 02C F C A n ζ ζ 4 + a 3 sl 02 C F n ζ 3 + a 3 s l 02CF 2 n ζ a 3 s l 11n d abc d abc n ζ ζ ζ 4 17

18 C ψ 2,12 = 1 + a s C F a 2 s C F C A ζ 3 + a 2 s C2 F ζ 3 + a 2 sc F n a 2sl 02 C F n a 3 s C F C A n ζ 4 + a 3 sc F C 2 A ζ ζ ζ ζ 5 + a 3 sc F n ζ 3 + a 3 s C2 F C A a 3 s C2 F n ζ ζ ζ ζ 4 + a 3 sc 3 F ζ ζ ζ ζ 5 + a 3 sl 02 C F C A n ζ ζ 4 + a 3 s l 02C F n ζ 3 + a 3 sl 02 CF 2 n ζ ζ 4 + a 3 d abc d abc sl 11 n n ζ ζ 5 18

19 C2,14 ns = 1 + a s C F a 2 s C FC A ζ 3 + a 2 s C Fn a 2 scf ζ 3 + a 3 s C FC A n ζ a 3 s C FC 2 A + a 3 s C Fn ζ ζ ζ ζ 3 + a 3 scfc 2 A ζ ζ 5 + a 3 scfn a 3 s C3 F ζ ζ 5 + a 3 s l 11n d abc d abc n ζ ζ ζ ζ ζ ζ ζ 5 19

20 C ψ L,10 = a s C F a 2 s C F C A ζ 3 + a 2 s C F n a 2 scf a 2 s l 02C F n ζ 3 + a 3 sc F C A n ζ 3 + a 3 s C F CA ζ ζ 5 + a 3 sc F n a 3 s C2 F C A ζ ζ 5 + a 3 scf 2 n ζ 3 + a 3 s C3 F ζ ζ 5 + a 3 sl 02 C F C A n ζ 3 + a 3 s l 02C F n a 3 sl 02 CFn ζ 3 + a 3 s l d abc d abc 11n n ζ ζ 5 20

21 C ψ L,12 = a s C F a 2sC F C A 13 + a 2 s C Fn ζ a 2 s C2 F + a 2 s l 02C F n a 3 s C Fn ζ 3 + a 3 sc F C A n ζ 3 + a 3 s C FC 2 A ζ ζ 5 + a 3 sc 2 FC A a 3 s C2 F n ζ ζ ζ a 3 scf ζ ζ 5 + a 3 s l 02C F C A n ζ 3 + a 3 sl 02 C F n a 3 s l 02CF 2 n ζ 3 + a 3 d abc d abc sl 11 n n ζ ζ 5 CL,14 ns = a s C F a 2 sc F C A ζ 3 + a 2 s C Fn a 2 scf ζ 3 + a 3 s C FC A n ζ 3 + a 3 s C FCA ζ ζ 5 + a 3 sc F n a 3 s C2 F C A ζ ζ 5 + a 3 scfn ζ 3 + a 3 s C3 F ζ ζ 5 + a 3 s l d abc d abc 11n n ζ ζ 5 21

22 C2,10 G = a s n a 2 s C An ζ 3 + a 2 s C F n ζ 3 + a 3 sc A n ζ 3 + a 3 s C2 A n ζ ζ ζ 5 + a 3 s C F C A n ζ ζ ζ 5 + a 3 s C F n ζ 3 + a 3 scf 2 n ζ ζ a 3 sl g d abc d abc 11 n2 N A 33 ζ ζ ζ 5 22

23 C2,12 G = a s n a 2 s C An a 2 s C F n a 3 sc A n 2 + a 3 s C2 A n 13 ζ ζ ζ ζ ζ ζ 3 + a 3 s C F C A n ζ ζ ζ 5 + a 3 s C F n ζ 3 + a 3 scf 2 n ζ ζ a 3 sl g d abc d abc 11 n2 N A CL,10 G = a s n a 2 s C An a 2 s C Fn a 3 sc A n ζ 3 + a 3 s C2 A n 91 ζ ζ ζ ζ ζ ζ ζ 5 + a 3 sc F C A n ζ 3 + a 3 s C Fn 2 + a 3 scfn ζ ζ 5 + a 3 sl g d abc d abc 11 n N A ζ ζ 5 23

24 CL,12 G = a s n a 2sC A n a 2sC F n a 3 s C An ζ 3 + a 3 s C2 A n ζ ζ 5 + a 3 sc F C A n ζ ζ 5 + a 3 s C F n ζ 3 + a 3 scf 2 n ζ ζ 5 + a 3 sl g d abc d abc 11 n N A ζ ζ 5 C Results or the moments o F 3 γ3 ns = a s C F a 2sC F C A a 2sC F n a scf a 3 s C FC A n ζ 3 + a 3 s C F CA ζ 3 + a 3 sc F n a 3sC 2FC A ζ 3 + a 3 s C2 F n ζ 3 + a 3 s C3 F ζ 3 + a 3 d abc d abc sn n 288 γ5 ns = a s C F a 2s 15 C F C A a 3 sc F C A n ζ 3 + a 3 s C F n a 3 sc 2 F n a 3 sn d abc d abc n ζ 3 + a 2s C F n a 2 s C2 F + a 3 sc F CA ζ 3 + a 3 s C2 F C A ζ 3 + a 3 scf ζ

25 γ7 ns = a s C F a 2 s C2 F + a 3 s C FC 2 A + a 2sC F C A a 3 sc 2 FC A a 3 sc 3 F γ9 ns = a s C F a 2sC F n a 3s C FC A n ζ ζ 3 + a 3 s C F n ζ 3 + a 3sC 2F n ζ ζ 3 + a 3sn d abc d abc n a 2s C F C A a 2s C F n a 3sC F C A n a 2 scf ζ 3 + a 3 s C F CA ζ 3 + a 3 s C F n a 3 scf 2 C A ζ 3 + a 3sC 2F n ζ 3 + a 3 s C3 F ζ 3 + a 3 s n d abc d abc n γ11 ns = a s C F a 2 s C2 F + a 3 sc F C 2 A + a 3 s C2 F C A + a 2s C F C A a 3s C FC A n + a 2s C Fn ζ ζ 3 + a 3 sc F n ζ ζ 3 + a 3 sc 2 Fn a 3 s C3 F + a 3 s n d abc d abc n ζ

26 γ13 ns = a s C F a 2sC F C A a 2 s C2 F + a 3 s C FC A n ζ 3 + a 3 sc F CA ζ 3 + a 3 s C Fn a 3 scfc 2 A ζ 3 + a 3 s C2 F n ζ 3 + a 3 scf ζ 3 + a 3 d abc d abc sn n a 2sC F n C ns 1 = 1 + a s C F 3 + a 2 s C FC A 23 + a 2 s C Fn +4 + a 2 s C2 F a 3 sc F C A n ζ ζ 5 + a 3 sc F CA 2 + a 3 s C Fn 2 + a 3 sc 3 F a 3 s C2 F C A + a 3sn d abc d abc n ζ ζ 3 C3 ns = 1 + a s C F a 2 sc F C A ζ 3 + a 2sC F n a 3 s C F C A n ζ ζ 4 + a 3 sc F C 2 A + a 3 s C F n 2 + a 3 sc 2 F n a 3 s C3 F a 3 s C2 F n + a 2 sc 2 F ζ ζ ζ ζ 3 + a 3 s C2 F C A ζ ζ ζ ζ 4 40ζ ζ ζ ζ ζ ζ 4 260ζ 5 + a 3 s n d abc d abc n ζ 3 26

27 C5 ns = 1 + a s C F a 2 s C FC A ζ 3 + a 2s C F n a 2 s C2 F ζ 3 + a 3 sc F C A n ζ ζ 4 + a 3 s C FCA ζ ζ ζ 5 + a 3 sc F n ζ 3 + a 3 s C2 F C A ζ ζ 4 180ζ 5 + a 3 scfn ζ ζ 4 + a 3 s C3 F ζ ζ 4 376ζ 5 + a 3 s n d abc d abc n ζ 3 C7 ns = 1 + a s C F a 2 s C FC A a 2 sc 2 F ζ 3 + a 2s C Fn ζ 3 + a 3 s C FC A n a 3 s C FCA a 3 sc F n ζ 3 + a 3 s C2 F C A ζ ζ a 3 scfn a 3 s C3 F + a 3 s n d abc d abc n ζ ζ ζ ζ ζ ζ ζ ζ ζ ζ ζ ζ 5 27

28 C9 ns = 1 + a s C F a 2 s C F C A a 2 s C2 F ζ ζ 3 + a 3 sc F C A n a 2s C F n ζ ζ 4 + a 3 s C F CA ζ ζ ζ 5 + a 3 sc F n ζ 3 + a 3 s C2 F C A ζ ζ ζ 5 + a 3 scf 2 n ζ ζ 4 + a 3 s C3 F ζ ζ ζ 5 + a 3 s n d abc d abc n ζ 3 C11 ns = 1 + a s C F a 2 sc F C A a 2 s C2 F ζ 3 + a 2sC F n ζ 3 + a 3 sc F C A n a 3 s C FC 2 A ζ 5 + a 3 s C Fn ζ 3 + a 3 s C2 F C A a 3 s C2 F n + a 3 sc 3 F ζ ζ ζ ζ ζ ζ ζ ζ ζ ζ 5 + a 3 s n d abc d abc n ζ ζ ζ 4

29 C13 ns = 1 + a s C F a 2 s C F C A ζ 3 + a 2 s C F n a 2 sc 2 F ζ 3 + a 3 s C F C A n a 3 sc F C 2 A ζ ζ 5 + a 3 sc F n ζ 3 + a 3 s C2 F C A ζ ζ ζ ζ ζ ζ 5 + a 3 s C2 F n ζ ζ 4 + a 3 scf ζ ζ ζ 5 + a 3 d abc d abc sn n ζ 3 29