arxiv: v2 [hep-ph] 7 Oct 2011

Transcription

1 Tensor form factors of nucleon in QCD arxiv:7.4584v2 hep-ph 7 Oc Güray Erkol and Altug Ozpineci 2 Laboratory for Fundamental Research, Ozyegin University, Kusbakisi Caddesi No:2 Altunizade, Uskudar Istanbul 34662, Turkey 2 Physics Department, iddle East Technical University, 653 Ankara, Turkey (Dated: June 5, 28) We extract the isovector tensor nucleon form factors, which play an important role in understanding the transverse spin structure of the nucleon when related to the quark helicity-flip generalized parton distributions via their first moments. We employ the light-cone QCD sum rules to leading order in QCD and include distribution amplitudes up to twist 6 in order to calculate the three tensor form factors H T, E T and H T. Our results agree well with those from other approaches in the low and high momentum-transfer regions. PACS numbers: 3.75.Gx, 4.2.Dh, e, 2.38.Gc Atthetwist-twolevelthequarkstructureofthenucleoncanbedescribedbythreequarkdistributionfunctionsf q (x), g q (x) and hq (x), where q is the corresponding quark. The unpolarized distribution, fq (x), and the spin-dependent helicity distribution, g q (x), give a measure of quark contribution to longitudinal momentum and the net helicity of a nucleon. These two chiral-even distribution functions are well-known quantities, which can be extracted from inclusive deep-inelastic scattering data. The transversity distribution, h q (x), measures the net density of quarks having parallel polarization to that of the transversely polarized nucleon. h q (x), being a chiral-odd quantity, is difficult to measure experimentally because it requires a chiral-odd probe to couple in the cross section. The transversely polarized Drell-Yan process, 2 and electro- and photo-production of mesons off the polarized nucleons 3 are suggested as appropriate ways to measure h q (x). Recently, hq (x) has been extracted for the first time by an analysis of the experimental data from Belle 4, HERES 5 and COPASS 6 collaborations on azimuthal single spin asymmetry in semi-inclusive deep-inelastic scattering 7. At the leading twist, a more general description of the nucleon structure is given by eight generalized parton distributions (GPD): Two chiral-even spin-independent distributions, H and E, two chiral-even spin-dependent distributions, H q and Ẽq and four chiral-odd spin-dependent distributions, H q T, Eq T, Hq T and Ẽq T 8. Important information about the quark-gluon structure of the nucleon is encoded in these quantities. They describe, e.g., how partons are distributed in the transverse plane with respect to motion of the nucleon or the contribution of quark orbital angular momentum to total spin of the nucleon (see, e.g., the reviews in Refs. 9). Takingthe moments ofthe GPDs one canobtain the nucleonform factorsofthe localcurrents. The first momentsof the helicity-conserving GPDs, H q, E q and H q, Ẽq, are constrained by the nucleon form factors of the electromagnetic and axial currents, respectively. Similarly, the helicity-flip GPDs, H q T, Eq T, Hq T and Ẽq T, can be related to the nucleon form factors of the tensor current by taking the first moment. ore explicitly, we have dx X(x,ξ,q 2 ) = X(q 2 ), for X = H q T, Eq T and H q T. Note that the first moment of Ẽq T (x,ξ,q2 ) vanishes therefore we have Ẽq T (q2 ) =. Here q = p p is the momentum transfer and ξ = n q/2 denotes the longitudinal momentum transfer, where n is a light-like vector. The nucleon matrix elements of the tensor isovector current are parameterized in terms of above form factors as follows: N(p ) T µν N(p) = ū(p ) iσ µν HT I= (q 2 )+ γ µq ν γ ν q µ ET I= (q 2 )+ P µq ν P ν q µ H I= 2m N 2m 2 T (q 2 ) u(p), () N where T µν = ūiσ µν u diσ µν d is the isovectortensor current, P = p +p, m N is the nucleon mass, and X I= X u X d for any of the form factors, X = H T, E T or H T (hereafter, the I = superscript will be dropped since we will be dealing with the I = case exclusively). The tensor form factors at Q 2 = q 2 = give what are called the tensor charges. Previously, they have been investigated in the framework of QCD sum rules, 2, relativistic quark models 3, axial-vector dominance model 4 and chiral-quark soliton model(χqs) 5. The dependence of the guray.erkol@ozyegin.edu.tr ozpineci@metu.edu.tr

2 tensor form factors on Q 2 have been extracted up to GeV 2 from χqs 6 and up to a few GeV 2 from lattice QCD 7, 8. Our aim in this work is to study the nucleon tensor form factors in the framework of light-cone QCD sum rules method (LCSR) 9. LCSR have proved to be rather successful in extracting the values of the hadron form factors at high-momentum transfers, e.g., the electromagnetic and the axial form factors of the nucleon have been calculated to leading order and with higher-twist corrections in Refs. 2, 2. To our knowledge, while various estimations of H T (Q 2 ) at low- and high-momentum transfers exist, with this study it is the first time in the literature E T (Q 2 ) and H T (Q 2 ) are extracted. In this method one starts with the following two-point correlation function: Π N µν(p,q) = i d 4 xe iqx Tη N ()T µν (x) N(p), (2) whereη N (x) is the nucleoninterpolatingfield. Thereareseverallocaloperatorswith the quantum numbersofspin-/2 baryons one can choose from. Here we work with the general form of the nucleon interpolating field parameterized as follows: 2 η N = 2ǫ abc (u at (x)cj l db (x))j2 l uc (x), (3) l= with J = I, J 2 = J 2 = γ 5 and J 2 2 = β, which is an arbitrary parameter that fixes the mixing of two local operators. We would like to note that when the choice β = is made the interpolating field above gives what is known as Ioffe current for nucleon. Here u(x) and d(x) denote the u- and d-quark fields, respectively, a, b, c are the color indices and C denotes charge conjugation. The short-distance physics corresponding to high momenta p 2 and q 2 is calculated in terms of quark and gluon degrees of freedom. Inserting the interpolating field in Eq. (3) into the correlation function in Eq. (2), we obtain 2 2 { Π N µν = d 4 xe iqx (CJ l 2 ) αγ J l 2 S u ( x)iσ µν ρβ 4ǫabc u a α ()ub β (x)dc γ () N l= +(J2 l ) ρα (CJ l ) T S u ( x)iσ µν γβ 4ǫabc u a α (x)ub β ()dc γ () N +(J2 l ) } ρβ CJ l S d ( x)iσ µν αγ 4ǫabc u a α ()ub β ()dc γ (x) N, (4) where u and d denote the quark fields. S u,d represent the quark propagator S q (x) = ix/ 2π 2 x 4 q q (+ m2 x 2 ) ig s dυ 2 6 x/ 6π 2 x 4G µνσ µν υx µ G µν γ ν i 4π 2 x 2. (5) Here the first term gives the hard-quark propagator. The second term represents the contributions from the nonperturbative structure of the QCD vacuum, namely, the quark and quark-gluon condensates. These contributions are removed by Borel transformations. The last term is due to the correction in the background gluon field and gives rise to four-particle (N = 4 with an additional gluon) and five-particle (N = 5 with an additional quark anti-quark pair) baryon distribution amplitudes. In Ref. 22, higher Fock states (N > 3) of the nucleon wavefunction have been considered and it has been shown by an application to electromagnetic form factors of the nucleon that the valence Fock states dominate while N = 4, 5 Fock states give only small contributions. Following the common practice, in this work we shall not take into account such contributions, which leaves us with the first term in Eq. (5) to consider. The matrix elements of the local three-quark operator 4ǫ abc u a α (a x)u b β (a 2x)d c γ (a 3x) N (6) (a,2,3 are real numbers in the interval, denoting the position of the quarks along the ray pointing from the origin to the space-time point x) can be expanded in terms of DAs using the Lorentz covariance, the spin and the parity of the baryon. Based on a conformal expansion using the approximate conformal invariance of the QCD Lagrangian up to -loop order, the DAs are then decomposed into local nonperturbative parameters, which can be estimated using QCD sum rules or fitted so as to reproduce experimental data. The long-distance side of the correlation function is obtained using the analyticity of the correlation function, which allows us to write the correlation function in terms of a dispersion relation of the form Π N µ (p,q) = π ImΠ N µ (s) (s p 2 ) ds

3 3 The ground-state hadron contribution is singled out by utilizing the zero-width approximation Im Π N µ = πδ(s m2 N ) ηn N(p ) N(p ) T µν N(p) expressingthe correlationfunction as a sharp resonanceplus continuum with the continuum threshold, s. The matrix element of the interpolating current between the vacuum and baryon state is defined as η N N(p,s) = λ N υ(p,s) where λ N is the nucleon overlap amplitude and υ(p,s) is the nucleon spinor. The QCD sum rules are obtained by matching the short-distance calculation of the correlation function with the long-distance calculation. Using the most general decomposition of the matrix element in Eq. (6) we obtain the following sum rules for the tensor form factors: H T = e m 2 N {m /2 N Y(,)+ 4λ N E T = e m 2 N {m /2 N Y(,)+ 2λ N H T = e m 2 N {m /2 N Y(,)+ 2λ N where we define Y(k,n) = } 2 m 3 2 N m3 2Y(2,4i+)+ N 2Y(,3)+ m 3 NZ(4i+)+m 3 Nx Z(3), i= i= } 2 m 3 2 N m3 2Y(2,4i+)+ N 2Y(,3)+ m 3 NZ(4i+)+m 3 Nx Z(3), i= i= } 2 m 3 2 N m3 2Y(2,4i+)+ N 2Y(,3)+ m 3 NZ(4i+)+m 3 Nx Z(3), i= d x ( ) ke s(t2)/ 2 {( β)f n +(+β)f n+ +( β)f n+2 +(+β)f n+3 }, Z(n) = e s/2 Q 2 +x 2 ( β)(f n +F n+2 )+(+β)(f n+ +F n+3 ), m2 N for H T form factor at the structure iq/σ µν, for E T form factor at the structure γ µ q ν γ ν q µ and for H T form factor at the structure P µ q ν P ν q µ. The explicit form of the functions that appear in the above sum rules are given in terms of DAs as follows: for H T : F = 2 F 2 = F 3 = 2 F 4 = 2 F 5 =, F 6 =, F 7 = F 8 = t2 t2 t2 t2 t2 t2 F 9 = 4 F = 4 F = dt P +T T 2 +T 7 (t,, t ), dt A +A 2 +2A 3 +V V 2 +2V 3 (t,, t ), dt A +A 3 V V 3 (t, t, ), dt P +S 2T T 3 T 7 (t, t, ), dt A +V (t, t, ), dt 3T (t, t, ), ρ ρ ρ i= dt T T 2 T 5 +T 6 2T 7 2T 8 (t,ρ, t ρ), dt A A 2 +A 3 +A 4 A 5 +A 6 V +V 2 +V 3 +V 4 +V 5 V 6 (t,ρ, t ρ), dt A A 2 +A 3 +A 4 A 5 +A 6 +V V 2 V 3 V 4 V 5 +V 6 (t, t ρ,ρ), (7) (8)

4 4 F 2 = F 3 =, F 4 =, F 5 =, F 6 =, ρ dt T +3T 2 2T 3 2T 4 +3T 5 T 6 +4T 7 +4T 8 (t, t ρ,ρ), for E T : t2 F = 2 F 2 = 2 F 3 = 2 F 4 = 4 F 5 = 2 F 6 = 2 F 7 =, F 8 = 4 F 9 = 2 F =, F = 2 t2 t2 t2 t2 t2 t2 F 2 = 4 F 3 = 2 F 4 = 2 F 5 = 2 F 6 = 2 dt dt P +S +2T T 3 T 7 (t,, t ), dt A +V (t,, t ), dt A +A 2 V +V 2 (t, t, ), dt T (t, t, ), dt T (t,, t ), dt A +V (t,, t ), dt T (t, t, ), dt t ρ ρ t ρ ρ ρ dt T +T 2 +T 5 T 6 +2T 7 +2T 8 (t,ρ, t ρ), dt A A 2 +A 3 +A 4 A 5 +A 6 +V V 2 V 3 V 4 V 5 +V 6 (t, t ρ,ρ), dt T 2 T 3 T 4 +T 5 +T 7 +T 8 (t, t ρ,ρ), P P 2 +S S 2 2T 2T 2 +T 3 +T 4 +2T 5 +3T 7 T 8 (t, t ρ,ρ), A +A 2 +2A 3 A 4 2A 5 +V +V 2 2V 3 +V 4 2V 5 (t, t ρ,ρ), dt A A 2 +A 3 V +V 2 +V 3 (t, t ρ,ρ), dt P P 2 S +S 2 +2T T 3 T 4 T 7 T 8 (t, t ρ,ρ), and for H T : F =, F 2 = t2 F 3 =, dt V A (t,, t ),

5 5 t2 F 4 = 2 dt T (t, t, ), F 5 =, F 6 = F 7 =, F 8 = 2 F 9 =, F =, F =, t2 t2 F 2 = 2 F 3 = F 4 = F 5 = 2 dt F 6 = dt dt V A (t,, t ), dt T (t, t, ), t t ρ ρ ρ dt T 2 T 3 T 4 +T 5 +T 7 +T 8 (t, t ρ,ρ), P P 2 +S S 2 4T 2T 2 +3T 3 +T 4 +2T 5 +5T 7 T 8 (t, t ρ,ρ), A +A 2 +2A 3 A 4 2A 5 +V +V 2 2V 3 +V 4 2V 5 (t,ρ, t ρ), dt A A 2 +A 3 V +V 2 +V 3 (t, t ρ,ρ), dt P P 2 S +S 2 +2T T 3 T 4 T 7 T 8 (t, t ρ,ρ). TheDA sp i,a i, S i andt i appearinthemostgeneraldecompositionofthenucleonmatrixelement23andthetwist-5 DA s A, V and T havebeen calculated in Ref. 2. In obtaining the sum rules, Borel transformationis applied to eliminate the subtraction terms in the spectral representation of the correlation function and to exponentially suppress the contributions from excited and continuum states. The contributions of the higher states and the continuum are modeled using the quark-hadron duality and subtracted. Both of the Borel transformation and the subtraction of the higher states are carried out using the following substitution rules: ρ(x) dx (q xp) 2 dx, x x ρ(x)e s(x)/2 ρ(x) dx (q xp) 4 (9) dx ρ(x) 2 + x x 2ρ(x)e s(x)/2 Q 2 +x 2 e s/2, m2 N where s(x) = ( x)m 2 N + x x Q2. is the Borel mass and x is the solution of the quadratic equation for s = s : x = (Q 2 +s m 2 N )2 +4m 2 N (Q2 ) (Q 2 +s m 2 N) /(2m 2 N), where s is the continuum threshold. To obtain a numerical prediction for the form factors, the residue, λ N, is also required. The residue can be obtained from the mass sum rule, which is given by: λ 2 Ne m2 N /2 = 6 24π 2(5+2β +5β2 )E 2 (x) m2 96 2( +β)2 qq 2 m2 8 2( +β2 ) qq ( +β2 ) qq 2, where x = s / 2, and E n (x) = e x n i= x i i!.

6 H T Exp. Fit. QCDSR CQS Lattice E T Exp. Fit. QCDSR ~ H T Exp. Fit QCDSR Q 2 (GeV 2 ) Q 2 (GeV 2 ) Q 2 (GeV 2 ) FIG.. Tensor form factors H T(Q 2 ), E T(Q 2 ) and H T(Q 2 ) as a function of Q 2 above Q 2 = GeV 2. We also show our fitted form factors (see the text for details) and give the results for H T as obtained from χqs 6 and from lattice QCD 7 We use the following parameter values: qq = (.243) 3 GeV 3, m 2 =.8 GeV2 and m N =.94 GeV. The DAs of the nucleon are given in Ref. 2 and they are expressed in terms of some nonperturbative parameters which are calculated using QCD sum rules or phenomenological models. In this work, we use Chernyak-Zhitnitsky-like model of the DAs (see Ref. 2 for details) to calculate the tensor form factors. Note that unlike other form factors such as axial ones, the tensor form factors are renormalization-scale, µ, dependent. Our results correspond to µ 2 = GeV 2 as we are using the numerical values of the DAs at this scale 24. There are several parameters in the sum rules that need to be determined. We use the continuum threshold value s 2.25 GeV 2, which is pretty much fixed in the literature from nucleon spectrum analysis. We also give our numerical results at 2 =2 GeV 2. The sum rules and hadron properties are expected to be independent of the mixing parameter β. In order to see if we can achieve such an independence, we make a reparameterization using β = tanθ and observe that a stability region with respect to a change in the mixing parameter can be found around cosθ, which we use to produce our results. We plot the tensor form factors H T (Q 2 ), E T (Q 2 ) and H T (Q 2 ) as a function of Q 2 in Fig.. We show the form factors above Q 2 = GeV 2. We also give the results for H T as obtained from χqs 6 and from lattice QCD 7, which have been calculated at a renormalization scale of µ 2 =.36 GeV 2 and µ 2 = 4 GeV 2, respectively. We use the following relation 25 to relate these results of the form factors to those at µ 2 = GeV 2 : X(µ 2 ) = ( αs (µ 2 ) 4/27 ) α s (µ 2 i ) π α s(µ 2 i ) α s(µ 2 ) X(µ 2 i ), () where µ i is the initial renormalization scale and α s (µ 2 ) = 4π 9ln(µ 2 /Λ 2 QCD ) 64 8 lnln(µ 2 /Λ 2 QCD ) ln(µ 2 /Λ 2 QCD ) () with Λ QCD =.25 GeV and N c = N f = 3. In order to extrapolate the form factors to low-momentum region, we have fitted QCDSR results to various forms including monopole and dipole, which unfortunately fail to give a reasonable description of data with a two-parameter fit. On the other hand, an exponential form, viz., F T (Q 2 ) = F T ()exp Q 2 /m 2 T, (2) allows a plausible description of data with a two-parameter fit. Note that the LCSR are reliable for Q 2 larger than a few GeV 2. We have, however, tried three fit regions, namely, Q 2 > GeV 2, Q 2 >.5 GeV 2 and Q 2 > 2 GeV 2. Our results are given in Table I. H T () values are to be compared with those from experiment as H T () = , from χqs as H T () =.23 6, from external-field QCD sum rules method as H T () =.29±.5 2, and from lattice QCD as H T () =.5±.5 7 at µ 2 = GeV 2. We find that our results for H T agree well with those from other theoretical approaches, when we make a fit from a more reliable high-q 2 region, namely from Q 2 > 2 GeV 2 (all theoretical methods overestimate H T when compared to the phenomenological result in 7). In Fig. we also show our fitted form factor as well as those from other theoretical approaches. In the high-q 2 region we observe that our H T results compare well with lattice-qcd results. With the finding that an exponential fit above 2 GeV 2 gives more

7 7 TABLE I. The values of exponential fit parameters, namely F T() and m T, of tensor form factors. We give the results of the fits from three regions. Form Factor Fit Region (GeV 2 ) F T() m T (GeV) H T E T H T accurate results for H T, we quote our final results for the other two form factors at Q 2 = as extrapolated from this region as E T () =.96 and H T () =.43. In conclusion, we have extracted the isovector tensor form factors of nucleon, which play an important role in understanding the transverse spin structure of the nucleon when related to the quark helicity-flip generalized parton distributions via their first moments. To this end, we have employed LCSR method to leading order in QCD and included distribution amplitudes up to twist 6. We have found that our results for the H T form factor compare well with lattice-qcd results in the high-q 2 region and with those from other theoretical approaches when extrapolated to low-q 2 regions via an exponential fit. We have also calculated the tensor form factors E T and H T for the first time in the literature. We will extend our work with a detailed analysis of tensor form factors with respect to each quark contribution so as to extract isoscalar, octet and singlet form factors as well as those of strange baryons, in a future publication. ACKNOWLEDGENTS This work has been supported by The Scientific and Technological Research Council of Turkey (TÜBİTAK) under project number T245. The work of A. O. is also partially supported by the European Union (HadronPhysics2 project Study of strongly interacting matter). R. Jaffe and X.-D. Ji, Phys.Rev.Lett. 67, 552 (99); Nucl.Phys. B375, 527 (992). 2 J. Cortes, B. Pire, and J. Ralston, Z.Phys. C55, 49 (992). 3 D. Ivanov, B. Pire, L. Szymanowski, and O. Teryaev, Phys.Lett. B55, 65(22); R. Enberg, B. Pire, and L. Szymanowski, Eur.Phys.J. C47, 87 (26);. El Beiyad, B. Pire,. Segond, L. Szymanowski, and S. Wallon, Phys.Lett. B688, 54 (2). 4 K. Abe et al., Phys.Rev.Lett. 96, 2322 (26). 5 A. Airapetian et al., Phys.Rev.Lett. 94, 22 (25). 6 E. Ageev et al., Nucl.Phys. B765, 3 (27). 7. Anselmino,. Boglione, U. D Alesio, A. Kotzinian, F. urgia, et al., Phys.Rev. D75, 5432 (27); Nucl.Phys.Proc.Suppl. 9, 98 (29); A. Bacchetta, A. Courtoy, and. Radici, Phys.Rev.Lett. 7, 2 (2). 8 X.-D. Ji, Phys.Rev.Lett. 78, 6 (997); A. Radyushkin, Phys.Rev. D56, 5524 (997); P. Hoodbhoy and X.-D. Ji, D58, 546 (998);. Diehl, Eur.Phys.J. C9, 485 (2). 9. Diehl, Phys.Rept. 388, 4 (23); X. Ji, Ann.Rev.Nucl.Part.Sci. 54, 43 (24). P. Hagler, Phys.Lett. B594, 64 (24). H.-x. He and X.-D. Ji, Phys.Rev. D52, 296 (995). 2 H.-x. He and X.-D. Ji, Phys.Rev. D54, 6897 (996). 3 I. Schmidt and J. Soffer, Phys.Lett. B47, 33 (997); B. Pasquini,. Pincetti, and S. Boffi, Phys.Rev. D72, 9429 (25). 4 L. P. Gamberg and G. R. Goldstein, Phys.Rev.Lett. 87, 242 (2). 5 H.-C. Kim,. V. Polyakov, and K. Goeke, Phys.Rev. D53, 475 (996); Phys.Lett. B387, 577 (996); C. Lorce, Phys.Rev. D79, 7427 (29). 6 T. Ledwig, A. Silva, and H.-C. Kim, Phys. Rev. D82, 3422 (2). 7. Gockeler et al., Phys.Lett. B627, 3 (25).

8 8 P. Hagler et al., Phys. Rev. D77, 9452 (28). 9. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl. Phys. B47, 385 (979); B47, 448 (979); L. J. Reinders, H. Rubinstein, and S. Yazaki, Phys. Rept. 27, (985); B. L. Ioffe and A. V. Smilga, Nucl. Phys. B232, 9 (984); V.. Braun and I. E. Filyanov, Z. Phys. C44, 57 (989); I. I. Balitsky, V.. Braun, and A. V. Kolesnichenko, Nucl. Phys. B32, 59 (989); V. L. Chernyak and I. R. Zhitnitsky, B345, 37 (99). 2 V.. Braun, A. Lenz, and. Wittmann, Phys. Rev. D73, 949 (26). 2 G. Erkol and A. Ozpineci, Phys.Rev. D83, 422 (2); T.. Aliev, K. Azizi, and A. Ozpineci, Nucl. Phys. A799, 5 (28). 22. Diehl, T. Feldmann, R. Jakob, and P. Kroll, Eur.Phys.J. C8, 49 (999). 23 V. Braun, R. J. Fries, N. ahnke, and E. Stein, Nucl. Phys. B589, 38 (2). 24 V. Chernyak and I. Zhitnitsky, Nucl.Phys. B246, 52 (984). 25. Gluck, E. Reya, and A. Vogt, Z.Phys. C67, 433 (995); V. Barone, A. Drago, and P. G. Ratcliffe, Phys.Rept. 359, (22). 8