arxiv: v1 [hep-ph] 8 May 2015

Transcription

1 , Gravitational form factors and transverse spin sum rule in a light front quark-diquark model in AdS/QCD Dipankar Chakraarti, 1 Chandan Mondal, 1 and Asmita Mukherjee arxiv:155.13v1 hep-ph 8 May 15 1 Department of Physics, Indian Institute of Technology Kanpur, Kanpur 816, India Department of Physics, Indian Institute of TechnologyBomay, Powai, Mumai 476, India (Dated: July, 18) Astract The gravitational form factors are related to the matrix elements of the energy-momentum tensor T µν. Using the light front wave functions of the scalar quark-diquark model for nucleon predicted y the soft-wall AdS/QCD, we calculate the flavor dependent A(Q ), B(Q ) and C(Q ) form factors. We also present all the matrix element of the energy-momentum tensor in a transversely polarized state. Further, we evaluate the matrix element of Pauli-Luanski operator in this model and show that the intrinsic spin sum rule involves the form factor C. The longitudinal momentum densities in the transverse impact parameter space are also discussed for oth unpolarized and transversely polarized nucleons. PACS numers: t, 14..Dh, x 1

2 I. INTRODUCTION Understanding the spin structure of the proton, which means how the proton spin (1/) is distriuted among its constituent quarks and gluons is one of the fundamental prolems in hadron physics. Most of the studies, oth theoretical and experimental, mainly aim at the longitudinal spin or helicity. Understanding the transverse spin and transverse angular momentum of the proton is a much more involved prolem. The complications associated with the transverse angular momentum are est understood in light-front framework, in which one gets an intuitive picture of deep inelastic electron-proton scattering processes. The longitudinal angular momentum operator is kinematical on the light-front, whereas the transverse angular momentum and rotation operators are dynamical. This implies that the partonic structure of the transverse spin is different from that of the longitudinal spin structure of the proton. Recently, several sum rules have een proposed in the literature aout transverse spin. In 1, a sum rule was derived in terms of the intrinsic transverse spin operator on the lightfront. Unlike the transverse rotation operator or Pauli-Luanski operator, matrix elements of the intrinsic spin operator are frame independent. This sum rule was explicitly verified in perturation theory for a dressed quark at one loop. Another sum rule was proposed in 3 involving the transversity distriution. A new transverse polarization sum rule was proposed in 4 6 which was interpreted at the partonic level. This was partially motivated y 7 where a relation etween the expectation values of equal time transverse rotation operators and the gravitational form factors is derived using delocalized states in the rest frame of the nucleon. Authors of 4, 5 analyzed the matrix elements of the transverse component of the Pauli- Luanski operator for a transversely polarized state and related it to the gravitational form factors A() and B(). In 8 1 it was pointed out that the aove result is frame dependent. In fact the only frame independent result is otained in terms of the intrinsic spin operators on the light-front; the corresponding relation not only involves A() and B() ut also the higher twist term C(), and the contriution from C() is not suppressed. In this work, we verify the statements made in 1 in a model calculation. Here, we evaluate the GFFs for a transversely polarized proton from the energy momentum tensor and verify the sum rule for the transverse spin in a light front quark-diquark model.

3 For this work, we take a phenomenological light front quark-diquark model recently proposed y Gutsche et al. 11 with the corrected parameters given in Ref.1. In this model, the diquark is considered to e scalar(i.e. scalar diquark model) and the light front wave functions for the proton are constructed from the wave functions otained in light front AdS/QCD correspondence13. The parameters in this model are fixed y fitting to the electromagnetic form factors of the nucleons. Using the overlap formalism of light front wave functions, we calculatethegffsfromtheenergymomentumtensor (T µν )foratransversely polarizedproton. The intrinsic spin operators which can e derived from the transverse components of the Pauli- Luansky operator are shown to satisfy the sum rule consistent with 1. The Fourier transform of the gravitational form factor in the impact parameter space has interesting interpretations14, 15. The Fourier transform of the form factor A(Q ) gives the longitudinal momentum density (p + density) in the transverse impact parameter space. We have evaluated this momentum density in our model. For unpolarized nucleon the momentum density is axially symmetric whereas for a transversely polarized nucleon, the deviation from the axially symmetric distriution is found to e dipolar in nature. In Sec.II, we descrie the model very riefly efore providing the results for the GFFs in Sec.III. In Sec.IV, we derive the matrix elements of the energy momentum tensor and the Pauli-Luanski operator for a transversely polarized proton state with transverse momentum P =. The sum rule for the intrinsic transverse spin otained from the Pauli-Luanski operator involves the GFFs A(), B() and C(). In Sec.VI, we summarize our main results. In Sec.V, the longitudinal momentum densities in the transverse impact parameter space for oth unpolarized and transversely polarized nucleons are discussed. The detail expressions of the matrix elements of T µν are provided in the appendix. II. LIGHT-FRONT QUARK-DIQUARK MODEL FOR THE NUCLEON In quark-scalar diquark model, the nucleon with three valence quarks is considered as an effectively composite system of a fermion and a neutral scalar ound state of diquark ased on one loop quantum fluctuations. The generic ansatz for the massless LFWFs as proposed in 11 3

4 is ψ + +q (x,k ) = ϕ (1) q (x,k ), ψ + q (x,k ) = k1 +ik xm n ϕ () q (x,k ), ψ+q (x,k ) = k1 ik ϕ () q xm (x,k ), n ψ q (x,k ) = ϕ (1) q (x,k ), (1) where ψ λ N λqq (x,k ) are the LFWFs with specific nucleon helicities λ N = ± and the struck quark q has a spin λ q = ±, where plus and minus correspond to + 1 and 1 respectively. For the nucleons, q can e either up (u) or down(d) quark. The functions ϕ (1) q (x,k ) and ϕ () q (x,k ) are the wave functions predicted y soft-wall AdS/QCD13 ϕ q (i) (x,k ) = N q (i) 4π log(1/x) (i) κ 1 x xa q (1 x) (i) q exp The normalizations of the Dirac and Pauli form factors are fixed as Parameters u d a (1).35. (1).8 1. a () () N (1) N () k log(1/x). () κ (1 x) TABLE I: List of the parameters used in the LF quark-diquark model for κ = 46.6 MeV F q 1(Q ) = n q I q 1 (Q ) I q 1(), Fq (Q ) = κ q I q (Q ) I q (), (3) so that F q 1() = n q and F q () = κ q where n u =, n d = 1 and the anomalous magnetic moments for the u and d quarks are κ u = and κ d =.33. The structure integrals, I q i (Q ) otained from the LFWFs have the form as I q 1 (Q ) = 1 dxx a(1) q (1 x) 1+(1) q R q (x,q )exp Q 4κ log(1/x), (4) 4

5 I q (Q ) = 1 dxx a(1) q 1 (1 x) +(1) q σ q (x)exp Q 4κ log(1/x), (5) with R q (x,q ) = 1+σq κ (x)(1 x) x Mn log(1/x) 1 Q 4κ log(1/x), (6) σ q (x) = N() q N (1) q x a() q a (1) q (1 x) () q (1) q. (7) In this model, the value of the AdS/QCD parameter κ is taken to e 46.6 MeV and the other parameters are fixed y fitting to the electromagnetic properties for the proton and neutron such as form factors, magnetic moments and charge radii1. For completeness, the parameters are listed in the Tale I. Here, we should mention that the value of the parameter κ depends on exact AdS/QCD model, here we use the value of κ = 46.6 MeV as determined y fitting the nucleon form factors with experimental data in Ref.16. III. GRAVITATIONAL FORM FACTORS The gravitational form factors(gffs) which are related to the matrix elements of the the stress tensor (T µν )play animportant roleinhadronicphysics. Foraspin1/composite system, the matrix elements of T µν involve four gravitational FFs 4, 1 P,S T µν i () P,S = Ū(P,S ) B i (q ) P µ Pν M n + (A i (q )+B i (q )) 1 (γµ Pν +γ ν Pµ ) + C i (q ) qµ q ν q g µν M n + C i (q )M n g µν U(P,S), (8) where P = (P +P )/ and q = P P and A(q ), B(q ), C(q ) and C(q ) are the GFFs. The spin-nonflip form factor A is analog of the Dirac form factor F 1. A(q ) allows us to measure the momentum fractions carried y each constituent of a hadron. According to Ji s sum rule17, J q = A q ()+B q (). Thus, one has to measure the spin-flip form factor B to find the quark contriutions to the nucleon spin. B(q ) is analogous to the Pauli form factor F for the vector current. In the light-front representation, one can easily compute spin-nonflip and spin-flip 5

6 GFFs y calculating the ++ component of the matrix elements of the the stress tensor 18 as Here we consider the Yukawa Lagrangian P +q, T++ i () (P + ) P, = A i(q ), (9) P +q, T++ i () (P + ) P, = (q1 iq ) B i(q ) M. (1) L = i ψγ µ ( µ ψ) ( µ ψ)γ µ ψ m ψψ + 1 ( µ φ)( µ φ) 1 λ φφ+gφ ψψ, (11) which leads to the corresponding energy momentum tensor as T µν = i ψγ µ ( ν ψ) ψγ µ ν ψ+( µ φ)( ν φ) g µν L. (1) Using the two particle Fock states for J z = + 1 and Jz = 1 and the light-front wave functions given in Eq.(1), we evaluate the GFFs A(q ) and B(q ) depending on different flavors (struck quark) as A(q ) = A q (q )+A (q ) = I 1q(q )+I 1 (q ) I q 1 (), (13) B(q ) = B q (q )+B (q ) = M n I q (q ) I (q ) I q (), (14) where A q/ (q ) and I q/ i (q ) are the GFFs and structure integrals corresponding to quark/scalar diquark. The explicit expressions of the structure integrals are listed in the Appendix. The integrals, I q i (Q ) in the denominators in the right hand side of Eq.(13) and Eq.(14) have the form as I q 1 (Q ) = I q (Q ) = 1 1 dxx a(1) q (1 x) 1+(1) q R q (x,q )exp Q 4κ log(1/x), (15) dxx a(1) q 1 (1 x) +(1) q σ q (x)exp Q 4κ log(1/x), (16) with R q (x,q ) and σ q (x) as defined in Eqs.(6,7). In Fig. 1(a) and (), we show total A(q ) and B(q ) depending on different struck quarks. The contriutions of quark and the diquark to the total spin-nonflip and spin-flip GFFs are shown in Fig. 1(c)-(f). One notices that at zero momentum transfer, A() = A q ()+A () = 1 and B() = B q ()+B () = as expected. 6

7 IV. MATRIX ELEMENTS OF THE ENERGY-MOMENTUM TENSOR Here we consider the transversely polarized state to calculate the matrix elements of T µν. The transversely polarized state (polarized along +ve x direction) is given y P,S (1) = 1 ) ( Ψ p(p +,P ) + Ψ p(p +,P ), (17) where Ψ p (Ψ p ) represents the two-particle Fock state corresponding to Jz = + 1 (J z = 1 ). For the transversely polarized state we calculate the matrix elements of T µν for Ψ p going to Ψ p and Ψ p going to Ψ p.to evaluate the right hand side of Eq.(8), we use the matrix elements of the different γ matrices listed in the appendix of the Ref. 1. Here we list only the final expressions of all the matrix elements for zero skewness and the detailed calculations are given in the appendix. Ψ p(p ) T ++ Ψ p(p) + Ψ p(p ) T ++ Ψ p(p) = (P + ) (I q I ) I q () (iq ) = B(Q ) (P+ ) (iq ), (18) M n Ψ p(p ) T +1 Ψ p(p) + Ψ p(p ) T +1 Ψ p(p) = (I 4q I 4 ) I q () P + (iq 1 q ) = B(Q ) P+ (iq 1 M q ). (19) n Ψ p(p ) T + Ψ p(p) + Ψ p(p ) T + Ψ p(p) = (I 6q I 6 ) I q () P + i(q ) = B(Q ) P+ M n i(q ). () Ψ p(p ) T + Ψ p(p) + Ψ p(p ) T + Ψ p(p) = (I 8q I 8 ) I q () (iq ) = A(Q )(M n ) B(Q ) (q ) +C(Q ) 4(q ) + M n M C(Q )(4M n )(iq ). (1) n We keep only terms linear in q, which are relevant for evaluation of the matrix elements of transverse spin. The matrix elements of T µν up to O(q) are: P,S (1) T ++ P,S (1) = (P + ) (I q I ) I q () (iq ) = B(Q ) (P+ ) (iq ), () M n P,S (1) T +1 P,S (1) =, (3) P,S (1) T + P,S (1) =, (4) P,S (1) T + P,S (1) = (I 8q I 8 ) I q () (iq ) = M n A(Q )+ C(Q )(iq ). (5) 7

8 1 struck quark: u struck quark: d.8 struck quark: u struck quark: d.8.6 A (Q ).6 B (Q ).4.4. (a) Q GeV () Q GeV.35.3 struck quark: u struck quark: d.3.5 struck quark: u struck quark: d.5. A q (Q )..15 B q (Q ) (c) Q GeV (d) Q GeV.8.7 struck quark: u struck quark: d A (Q ).5 B (Q ) struck quark: u struck quark: d (e) Q GeV (f) Q GeV FIG. 1: (Color online) Plots of gravitational form factors A(Q ) and B(Q ) for oth struck quarks u and d. The contriutions of quark and diquark to the nucleon are shown in (c)-(f). We see from the Eq.() that the matrix element of T ++ in a transversely polarized state doesnot depend onthe formfactora(q ). It depends onlyontheformfactor B(Q ). Whereas, Eq.(5)impliesthatthematrixelement oft + dependsonoth C(Q )anda(q ). Fornonzero skewness, the matrix element of T + has a term proportional to q + as shown in 1. But the 8

9 .1. C (Q ) struck quark: u struck quark: d (a) Q GeV Cq (Q )..3.4 C (Q ) struck quark: u.5 struck quark: u struck quark: d struck quark: d () Q GeV (c) Q GeV FIG. : (Color online) Plots of gravitational form factor C(Q ) for oth struck quarks u and d. The contriutions of quark and diquark to the nucleon are shown in () and (c). matrix element of T + in this quark-diquark model is zero when we consider only the term linear in q. The main reason is that the LFWFs are independent of quark mass in this model. For zero skewness, the results of this quark-diquark model are consistent with Ref. 1. Using Eq.(13) and (5), we evaluate the C(Q ) form factor. In Fig. we show the form factor C(Q ) for different struck quark. The quark and the diquark contriutions are shown in Fig. () and (c) respectively. A. Matrix element of the Pauli-Luanski operator The Pauli-Luanski operator is defined as 1 W 1 i = 1 F i P+ + K 3 P 1 Ẽ i P, (6) 9

10 where F i and E i are the light front transverse rotation and transverse oost operators. K 3 is the longitudinal oost operator. The matrix elements of the operators F i, K3 and Ẽ i in a transversely polarized state are given y PS (1) Fi PS(1) = i(π) 3 δ 3 () P S (1) T + i () PS (1) and P S (1) T + i () PS (1) q=, (7) PS (1) K i PS 3 (1) = i (π)3 δ 3 () P S (1) T ++ i () PS (1), (8) q= PS (1) Ẽ i PS(1) = i(π) 3 δ 3 () P S (1) T ++ i () PS (1), (9) q= where = P P. Using the results of the individual matrix element in Eqs. (7), (8) and (9), the matrix element of the total Pauli-Luanski operator W 1 can e written in this quark-diquark model as PS (1) W 1 PS (1) PS (1) PS (1) = PS(1) W 1 PS (1) (π) 3 P + δ 3 () = 1 P + { (I9q I +I9 I ) } P + ζ I() q q= { + P+ I 8q() I 8 () { }+ P I q I() q (P+ ) () I() q I () } I() q = 1 P + P + M n{a()+b()+ C()}+ P (P+ ) B() M n = M n A()+B()+ C(). (3) One can notice that as B() = (Fig. 1-), only the matrix element of T + make contriution to the matrix element of total W 1 operator in a transversely polarized state. For different struck quarks, the matrix element of total W 1 operator will give different values due as C() is different for u and d quarks. The non-vanishing contriution of C to the matrix element of Pauli-Luanski operator W 1 has een reported previously in 1, 19,. For a massive particle like nucleon, the intrinsic spin operators can e related with the Pauli-Luanski operators through the following relations 1,, 1, M n J i = W i P i J 3 = ǫ ij ( 1 Fj P + +K 3 P j 1 Ej P ) P i J 3, 1

11 J 3 = W+ P + = J3 + 1 P +(E1 P E P 1 ), (31) where J 3 is the helicity operator. The matrix element of the intrinsic spin operator in a transversely polarized dressed quark state has een explicitly demonstrated in. The matrix element of the intrinsic spin operator in a transversely polarized state in the quark-diquark model is given y PS (1) J 1 PS (1) PS (1) PS (1) = 1 M n PS (1) W 1 PS (1) PS (1) PS (1) = 1 A()+B()+ C(). (3) In 1, the authors have descried in detail whether the matrix elements of the intrinsic spin operator are frame independent or not whereas the authors in 4, 5 have claimed that the results are frame independent though they have calculated only in a frame where P = (a) ρ u () fm fm fm () hm fm (c) ρ T u () fm fm fm (d) hm fm FIG. 3: (Color online) The longitudinal momentum densities for the active u quark in the transverse plane, upper panel for unpolarized nucleon, lower panel for nucleon polarized along x-direction. () and (d) are the top view of (a) and (c) respectively. 11

12 (a) ρ d () fm fm fm () hm fm (c) ρ T d () fm fm fm (d) hm fm FIG. 4: (Color online) The longitudinal momentum densities for the active d quark in the transverse plane, upper panel for unpolarized nucleon, lower panel for nucleon polarized along x-direction. () and (d) are the top view of (a) and (c) respectively. V. LONGITUDINAL MOMENTUM DENSITY IN TRANSVERSE PLANE According to the standard interpretation 14, 1 4 the charge and anomalous magnetization densities in the transverse plane can e identified with the two-dimensional Fourier transform(ft) of the electromagnetic form factors in the light-cone frame with q + = q + q 3 =. Similar to the electromagnetic densities, one can evaluate the gravitomagnetic density in transverse plane y taking the FT of the gravitational form factor 14, 15. Since the longitudinal momentum is given y the ++ component of the energy momentum tensor P + = dx d x T ++, (33) 1

13 (a) ρ u N () fm fm fm () hm fm (c) N ρ () fm Tu fm fm (d) hm fm 6 4 FIG. 5: (Color online) The longitudinal momentum densities of nucleon in the transverse plane for the active quark u, upper panel for unpolarized nucleon, lower panel for nucleon polarized along x-direction. () and (d) are the top view of (a) and (c) respectively. it is possile to interpret the Fourier transform of the gravitational form factor A(Q ) as the longitudinal momentum density in the transverse plane 15. For a unpolarized nucleon the momentum density can e defined as q ρ() = )e iq. = (π) A(Q dq π QJ (Q)A(Q ), (34) where = represents the impact parameter and J is the cylindrical Bessel function of order zero and Q = q. Under isospin symmetry, the momentum density is same for oth proton and neutron. Due to polarization, the density gets modified y a term which involves the spin flip form factor B(Q ). For transversely polarized nucleon, the momentum density is 13

14 (a) ρ d N () fm fm fm () hm fm (c) N ρ () fm Td fm fm (d) hm fm FIG. 6: (Color online) The longitudinal momentum densities of nucleon in the transverse plane for the active quark d, upper panel for unpolarized nucleon, lower panel for nucleon polarized along x-direction. () and (d) are the top view of (a) and (c) respectively. given y15 ρ T () = ρ()+sin(φ φ s ) dq π Q M n J 1 (Q)B(Q ), (35) where M n is the mass of nucleon. The transverse polarization of the nucleon is given y S = (cosφ sˆx+sinφ s ŷ) and the transverse impact parameter is denoted y = (cosφ ˆx+sinφ ŷ). Without loss of generality, the polarization of the nucleon is chosen along x-axis ie., φ s =. The second term in Eq.(35), provides the deviation from circular symmetry of the unpolarized density. Results for the momentum density ρ() for the active quark u for oth unpolarized and the transversely polarized nucleon are shown in Fig.3. Similar plots for the active quark d are shown in Fig.4. The plots show that the unpolarized densities are axially symmetric and have 14

15 hm hm (a) fm () fm hm. hm (c) fm 1 (d) fm 1 FIG. 7: (Color online) The momentum density asymmetry (ρ T () ρ()) in the transverse plane for a nucleon polarized in x-direction, (a) for u quark () for d quark (c) for nucleon when the active quark is u and (d) for nucleon when the active quark is d. the peak at the center of the nucleon( = ). For the nucleon polarized along x-direction, the peak of densities gets shifted towards positive y-direction and the densities no longer have the symmetry. It can also e noticed that the width of the density for d quark is larger ut the height of the peak is sufficiently small compare to u quark. In Fig.5 and Fig.6, we show the totalnucleonp + densities (fermionicplusosonic) foroththeunpolarizedandthetransversely polarized nucleon for different active quark u and d respectively. The total angular-dependent part of the densities are small. So the shifting of the densities in ρ T () are also very small and it is practically invisile in Fig.5(d) and Fig.6(d). However, removing the axially symmetric part of the density from ρ T () i.e, if we look at (ρ T () ρ()), one can find that the total angular-dependent part of the density(i.e. distortion from the symmetry) displays a dipole 15

16 pattern (Fig.7(c) and (d)). The angular-dependent part of the densities for active quarks u and d are shown in Fig.7(a) and () respectively. Both plots show the dipole pattern ut it is roader for d quark than u quark. Though the dipolar distortion for individual quark is quite large when oth quark and diquark contriutions are added together for the nucleon, the distortions ecome small irrespective of struck quark flavor(fig.7(c) and (d)). VI. SUMMARY The main result of this work is to show that the sum rules for the intrinsic spin for a transversely polarized proton involve the form factors A q, B q and C q in agreement with the claim in Ref.1. We demonstrated this in a recently proposed light front quark-diquark model where the LFWFs are modeled from the wave functions otained from light front AdS/QCD. We have also shown explicit Q ehavior of the gravitational form factors in this model. We have evaluated the longitudinal momentum density (p + density) in the transverse plane for oth unpolarized and polarized nucleon. For transversely polarized nucleon, the asymmetries in the distriutions for an individual quark is quite large, ut when the contriutions from the quark and the osonic diquark are considered, the overall asymmetries in the nucleon ecome small ut are shown to e dipolar in nature. Appendix A: Matrix elements of T µν 1. T ++ : up going to up matrix element T ++ = i ψγ + ( + ψ) ψγ + + ψ+( + φ)( + φ) = iψ +( + ψ + ) ψ + + ψ+ +( + φ)( + φ). (A1) Ψ p(p +,P = q) T ++ q Ψ p(p +,P = ) = Fq 1 (Q ) I q 1(), Ψ p(p +,P = q) T ++ Ψ p(p +,P = ) = F 1 (Q ) I q 1(), (A) 16

17 where F1(Q q ) = (P + ) k dx 16π x 3 ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 ) ( = (P + ) dxx x a 1 (1 x) 1+1 N ) x a + (1 x) +3 1 N 1 Mn ( κ ) log(1/x) Q exp log(1/x) Q 4 κ 4 = (P + ) I 1q, (A3) F1 (Q ) = (P + ) k dx (1 x) ψ k 16π3 + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 ) = (P + ) dx(1 x) x a 1 (1 x) ( κ log(1/x) x Q ) 4(1 x) = (P + ) I 1. exp ( N N 1 ) x a (1 x) +3 1 M n log(1/x) κ x Q 4(1 x) (A4) Ψ p(p +,P = q) T ++ Ψ p(p +,P = ) = (P + ) I1q (Q ) I q 1 () + I 1(Q ) I q 1 (). (A5) Using the matrix elements Eq.(8), Ψ p(p +,P = q) T ++ Ψ p(p +,P = ) = (P + ) A(Q ). (A6) So, A(Q ) = I1q (Q ) I q 1() + I 1(Q ) I q. (A7) 1(). T ++ : up going to down plus down going to up matrix elements Ψ p(p ) T ++ q Ψ p(p) + Ψ p(p ) T ++ q Ψ p(p) = Fq (Q ) I q (), Ψ p(p ) T ++ Ψ p(p) + Ψ p(p ) T ++ Ψ p(p) = F (Q ) I q (), (A8) 17

18 where F q (Q ) = (P + ) k dx 16π x {ψ k 3 + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} +{ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} = (P + ) (iq ) dx N 1 x a 1+a (1 x) 1+ + exp log(1/x) Q N 1 M n κ 4 = (P + ) (iq )I q, (A9) F (Q ) = (P + ) k dx 16π (1 x) {ψ k 3 + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} +{ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} = (P + ) (iq ) dx N 1 x a 1+a (1 x) 1+ + exp log(1/x) x Q N 1 M n κ (1 x) 4 = (P + ) (iq )I. (A1) Ψ p (P ) T ++ Ψ p (P) + Ψ p (P ) T ++ Ψ p (P) Using the matrix elements from Eq.(8), So, = (P + ) Iq (Q ) I q () I (Q ) I q () Ψ p(p ) T ++ Ψ p(p) + Ψ p(p ) T ++ Ψ p(p) = B(Q ) (P+ ) (iq ). M n (iq ). (A11) (A1) B(Q Iq (Q ) ) = M n I q () I (Q ) I q (). (A13) 3. T +1 : up going to up matrix element T +1 = i ψγ + ( 1 ψ) ψγ + 1 ψ+( + φ)( 1 φ) = iψ +( 1 ψ + ) ψ + 1 ψ+ +( + φ)( 1 φ). (A14) Ψ p(p +,P = q) T +1 q Ψ p(p +,P = ) = Fq 3 (Q ) I q 1(), Ψ p(p +,P = q) T +1 Ψ p(p +,P = ) = F 3 (Q ) I q 1 (), (A15) 18

19 where F3(Q q ) = P + k dx 16π 3 ( k 1 )ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 ) { = P + dx x a 1 (1 x) 1+ + ( κ log(1/x) Q 4 exp log(1/x) Q κ 4 )} q 1 + ( N ( N N 1 ) x a (1 x) +4 1 N 1 ) x a (1 x) +4 1 M n M n κ log(1/x) (iq ) = P + (q 1 I 3q +iq I 3q ), (A16) F3(Q ) = P + k dx 16π 3 k 1 ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 ) { = P + dx x a1+1 (1 x) ( κ log(1/x) x Q )} ( q 1 4(1 x) + N exp log(1/x) x Q κ (1 x) 4 ( N ) x a 1 (1 x) +3 1 N 1 Mn ) x a 1 (1 x) +3 1 N 1 M n κ log(1/x) (iq ) = P + (q 1 I 3 +iq I 3 ), (A17) Ψ p (P ) T +1 Ψ p (P) = (I 3q +I 3 ) I q 1() Using the matrix elements from Eq.(8), P + q 1 + (I 3q +I 3 ) I q P + (iq 1() ). (A18) Ψ p(p ) T +1 Ψ p(p) = A(Q )P + (q 1 )+ 1 (A(Q )+B(Q ))P + ( iq ). (A19) Comparing Eq. A18 and Eq. A19, A(Q ) = (I 3q +I 3 ) I q, 1() A(Q )+B(Q ) = (I 3q +I 3 ) I q 1(). (A) 4. T +1 : up going to down plus down going to up matrix elements Ψ p(p ) T +1 f Ψ p(p) + Ψ p(p ) T +1 f Ψ p(p) = Fq 4 (Q ) I q (), 19

20 Ψ p(p ) T +1 Ψ p(p) + Ψ p(p ) T +1 Ψ p(p) = F 4 (Q ) I q (), (A1) where F4(Q q ) = P + k dx 16π 3 ( k 1 ) {ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} +{ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} = P + (iq 1 q ) dx N 1 x a 1+a 1 (1 x) exp log(1/x) Q N 1 M n κ 4 = P + (iq 1 q )I 4q, (A) F4 (Q ) = P + k dx 16π 3 k 1 {ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} +{ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} = P + (iq 1 q ) dx N N 1 1 M n x a 1+a +1 (1 x) exp log(1/x) x Q κ (1 x) 4 = P + (iq 1 q )I 4. (A3) Ψ p (P ) T +1 Ψ p (P) + Ψ p (P ) T +1 Ψ p (P) = (I 4q I 4 ) I q P + (iq 1 () q ). (A4) Using the matrix elements from Eq.(8), Ψ p (P ) T +1 Ψ p (P) + Ψ p (P ) T +1 Ψ p (P) = B(Q ) P+ M n (iq 1 q ). (A5) We keep only terms which are linear in q, then ignoring the term q 1 q, Ψ p (P ) T +1 Ψ p (P) + Ψ p (P ) T +1 Ψ p (P) =. (A6) 5. T + : up going to up matrix element T + = i ψγ + ( ψ) ψγ + ψ+( + φ)( φ) = iψ +( ψ + ) ψ + ψ+ +( + φ)( φ). (A7)

21 Ψ p (P+,P = q) T + q Ψ p (P+,P = ) = Fq 5(Q ) I q 1(), Ψ p (P+,P = q) T + Ψ p (P+,P = ) = F 5 (Q ) I q 1(), (A8) where F5(Q q ) = P + k dx 16π 3 ( k )ψ k + (x, )ψ { = P + dx x a 1 (1 x) 1+ + ( κ log(1/x) Q 4 exp log(1/x) Q κ 4 )} q + ( N (x, k )+ψ k (x, )ψ (x, k 1 1 ) ( N N 1 ) x a (1 x) +4 1 N 1 ) x a (1 x) +4 1 M n M n κ log(1/x) ( iq1 ) = P + (q I 5q iq 1 I 5q ), (A9) F5 (Q ) = P + k dx 16π 3 k ψ k + (x, )ψ { = P + dx x a1+1 (1 x) ( κ log(1/x) x Q )} q 4(1 x) + exp log(1/x) x Q κ (1 x) 4 ( N (x, k )+ψ k (x, )ψ (x, k 1 1 ) N 1 ) x a 1 (1 x) +3 1 M n ( N N 1 ) x a 1 (1 x) +3 1 M n κ log(1/x) ( iq1 ) = P + (q I 5 iq 1 I 5 ), (A3) Ψ p (P ) T + Ψ p (P) = (I 5q +I 5 ) I q 1 () P + q (I 5q +I 5 ) I q 1 () P + (iq 1 ). (A31) Using the matrix elements from Eq.(8), Ψ p(p ) T + Ψ p(p) = A(Q )P + q + 1 (A(Q )+B(Q ))P + (iq 1 ). (A3) Comparing Eq. A31 and Eq. A3, A(Q ) = (I 5q +I 5 ) I q 1 (), A(Q )+B(Q ) = (I 5q +I 5 ) I q 1() 1. (A33)

22 6. T + : up going to down plus down going to up matrix elements where Ψ p(p ) T + f Ψ p(p) + Ψ p(p ) T + f Ψ p(p) = Fq 6 (Q ) I q (), Ψ p(p ) T + Ψ p(p) + Ψ p(p ) T + Ψ p(p) = F 6 (Q ) I q (), F6(Q q ) = P + k dx 16π 3 ( k ) {ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} (A34) +{ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} = P + i(q ) dx N 1 x a 1+a 1 (1 x) exp log(1/x) Q N 1 M n κ 4 = P + i(q ) I 6q, (A35) F6(Q ) = P + k dx 16π 3 k {ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} +{ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ 1 = P + i(q ) (x, k )} dx N 1 x a 1+a +1 (1 x) exp N 1 M n log(1/x) x Q κ (1 x) 4 = P + i(q ) I 6. (A36) Ψ p(p ) T + Ψ p(p) + Ψ p(p ) T + Ψ p(p) = (I 6q I 6 ) I q () P + i(q ). (A37) Using the matrix elements from Eq.(8), Ψ p (P ) T + Ψ p (P) + Ψ p (P ) T + Ψ p (P) = B(Q ) P+ M n i(q ). Ignoring the term (q ), (A38) Ψ p (P ) T + Ψ p (P) + Ψ p (P ) T + Ψ p (P) =. (A39) 7. T + : up going to up matrix element T + = ψ +m + ψ i ( φ) + 1 λ φ +interaction terms. (A4)

23 where F q 7 (Q ) = = Ψ p (P+,P = q) T + q Ψ p (P+,P = ) = Fq 7(Q ) I q 1(), Ψ p (P+,P = q) T + Ψ p (P+,P = ) = F 7 (Q ) I q 1(), k dx(k +m ) ψ k 16π 3 x + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 ) x x a 1 (1 x) 1+1 {m + κ (1 x) x log(1/x) + (1 x) Q } ( N ) x a + 4 N 1 (1 x) +3 1 {( κ )m log(1/x) Q Q (1 x) ( κ ) 4 4 log(1/x) + Q 4 + κ log(1/x) = I 7q, M n ( κ (1 x) log(1/x) + Q 4 )} exp log(1/x) Q κ 4 (A41) (A4) F7 k dx 1 ( (Q ) = k 16π 3 q ) Q 1 x 4 +λ ψ k + (x, )ψ (x, k ) +ψ k (x, )ψ (x, k 1 1 ) = dx x a 1 (1 x) 1 {(λ Q 4 )+ κ (1 x) log(1/x) + (1 x) Q } ( N ) x a + 4 N 1 (1 x) 1 { κ 4 (1 x) 4 Mn (log(1/x)) + κ (1 x) log(1/x) (λ xq ) x Q ( Q )} 4 4 (x x)+λ exp log(1/x) x Q κ (1 x) 4 = I 7, (A43) Ψ p(p +,P = q) T + Ψ p(p +,P = ) = I 7q +I 7 I q 1(). (A44) The diquark mass λ appears only in C(Q ) and is taken to e.6gev. Using the matrix elements from Eq.(8), Ψ p(p +,P = q) T + Ψ p(p +,P = ) = A(Q )(M n + (q ) If we ignore the (q ) dependent term, ) B(Q ) (q ) +C(Q )4(q ) + C(Q )(4M n). (A45) Ψ p(p +,P = q) T + Ψ p(p +,P = ) = A(Q )(M n)+ C(Q )(4M n). (A46) 3

24 So A(Q )+ C(Q ) = 1 I 7q +I 7 M n I q. (A47) 1() 8. T + : up going to down plus down going to up matrix elements where F q 8 (Q ) = F 8(Q ) = Ψ p(p ) T + q Ψ p(p) + Ψ p(p ) T + q Ψ p(p) = Fq 8 (Q ) I q (), Ψ p(p ) T + Ψ p(p) + Ψ p(p ) T + Ψ p(p) = F 8 (Q ) I q (), k dx(k +m ) {ψ k 16π 3 x + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} +{ψ k + (x, )ψ dx N 1 N 1 = iq exp log(1/x) = iq I 8q, (x, k )+ψ 1 (x, k )ψ 1 (x, k )} x a 1+a (1 x) 1+ + κ (1 x) M n Q κ 4 log(1/x) + Q (1 x) 4 k dx 1 ( k 16π 3 q ) Q 1 x 4 +λ {ψ k + (x, )ψ (x, k ) +ψ = iq k (x, )ψ 1 1 +m (A48) (A49) (x, k )}+{ψ k + (x, )ψ (x, k )+ψ k 1 (x, )ψ (x, k 1 )} dx N 1 x a 1+a (1 x) 1+ κ (1 x) N 1 M n log(1/x) + Q (1 x) +(λ Q 4 4 ) exp log(1/x) x Q κ (1 x) 4 = P + i(q )I 8. (A5) Ψ p(p ) T + Ψ p(p) + Ψ p(p ) T + Ψ p(p) = (I 8q I 8 ) I q () (iq ). (A51) Using the matrix elements from Eq.(8), Ψ p(p ) T + Ψ p(p) + Ψ p(p ) T + Ψ p(p) = A(Q )(M n ) B(Q ) (q ) M n +C(Q ) 4(q ) M n 4 + C(Q )(4M n )( iq ). (A5)

25 Ignoring the (q ) dependent term, So Ψ p (P ) T + Ψ p (P) + Ψ p (P ) T + Ψ p (P) = A(Q )(M n )+ C(Q )(4M n )( iq ). (A53) A(Q )+ C(Q ) = 1 M n I 8q I 8 I q (). (A54) The interaction terms will not contriute in the process. It will contriute only when we consider the higher order corrections. Appendix B: Matrix elements of T + for nonzero skewness To calculate the matrix element of Pauli-Luanski operator, first we need to evaluate the matrix elements of T + for nonzero skewness. For this section we use the following frame as ) P = (P +,P,P ) = (P +,, M, P + P = (P +,P,P ) = q = P P = where t = ζ M +q 1 ζ and q = Q. ( (1 ζ)p +, q, q ) +M, (1 ζ)p + ( ) ζp +,q, t+q, (B1) ζp + 1. T + : up going to down plus down going to up matrix elements Ψ p(p ) T + f Ψ p(p) + Ψ p(p ) T + f Ψ p(p) = Fq 9(ζ,Q ) I q (), Ψ p (P ) T + Ψ p (P) + Ψ p (P ) T + Ψ p (P) = F 9(ζ,Q ) I q (), (B) where F q 9(ζ,Q ) 5

26 = P + k dx ( 1 x ) 1/ 1 ζ( k 16π 3 1 x ) {ψ + (x, k )ψ (x, k ) (x, k )} + ψ (x, k )ψ (x, k 1 1 )}+{ψ + 1 (x, k )ψ (x, k + 1 )+ψ 1 (x, k )ψ 1 = 4iP +N 1 ( 1 x ) 1/ logxlogx 1/ { dx 1 ζ x a 1 (1 x) 1 N 1 κ M n 1 x (1 x)(1 x ) x a 1 (1 x ) x a 1 (1 x ) 1 x a 1 (1 x) 1 A +(q ) (logx ) 4κ 4 (1 x ) A 3 q + x a 1 (1 x) 1 x a 1 (1 x ) +1 (q logx } Q logx ( logx ) ) exp +1 κ (1 x )A q κ κ (1 x ) A q = ip + I9q(ζ,Q I )+ip + (q ) I9q(ζ,Q II ), (B3) F9 (ζ,q ) = P + k dx ( 1 x ) 1/ 1 ζ k 16π 3 {ψ 1 x + (x, k )ψ (x, k ) + x a 1 (1 x) 1 x a (1 x ) (q ) x logx κ (1 x ) A = ip + I9(ζ,Q I )+ip + (q ) I9(ζ,Q II ), (x, k )} + ψ (x, k )ψ (x, k 1 1 )}+{ψ + 1 (x, k )ψ (x, k + 1 )+ψ 1 (x, k )ψ 1 = 4iP +N 1 ( 1 x ) 1/ logxlogx 1/ { dx 1 ζ x a 1 (1 x) 1 N 1 κ M n 1 x (1 x)(1 x ) x a 1 (1 x ) x a 1 (1 x ) 1 x a 1 (1 x) 1 +(q x (logx ) A ) 4κ 4 (1 x ) 4 A 3 } x logx exp with Q logx κ ( ) +1 κ (1 x ) 4 A (B4) A q (x) = log(1/x ) κ (1 x ) + log(1/x) κ (1 x), A (x) = log(1/x ) κ (1 x ) + log(1/x) κ (1 x), (B5) where x = x ζ 1 ζ and k = k 1 x 1 ζ q for the struck quark and x = x 1 ζ and k = k + x 1 ζ q for the scalar diquark. So, Ψ p(p ) T + Ψ p(p) + Ψ p(p ) T + Ψ p(p) = (II 9q +II 9 ) I q () (ip + )+ (III 9q +III I q () 9 ) (ip + (q ) ). (B6) 6

27 Using the matrix elements from Eq.(8): Ψ p(p ) T + Ψ p(p) + Ψ p(p ) T + Ψ p(p) = 1 ζ( ζ) A(Q )+B(Q ) M n (ip + ) 1 ζ + B(Q ζ ) +C(Q ζ ) ip + (q 4M n 1 ζ M n 1 ζ ). Ignoring the term (q ) and comparing Eqs.(B6) and (B7), (I I 9q +I I 9 ) I q () = 1 ζ( ζ) A(Q )+B(Q ) M n. 1 ζ (B7) (B8) So ( Ψ p q (P ) T + Ψ p (P) + Ψ p (P ) T + Ψ p (P) ). q= (I9q I +I9 I = i ) ζ I q q= () = im n A()+B(). (B9) 1 A. Harindranath, A. Mukherjee and R. Rataole, Phys. Lett. B 476, 471 () A. Harindranath, A. Mukherjee and R. Rataole, Phys. Rev. D 63, 456 (1). 3 B. Bakker, E. Leader and T. L. Trueman, Phys. Rev. D 7, 1141 (4). 4 X. Ji, X. Xiong and F. Yuan, Phys. Lett. B 717, 14 (1). 5 X. Ji, X. Xiong and F. Yuan, Phys. Rev. Lett. 19, 155 (1). 6 X. Ji, X. Xiong, F. Yuan, Phys. Rev. Lett 111, 3913 (13). 7 M. Burkardt, Phys. Rev. D 7, 94 (5). 8 E. Leader, C. Lorce, Phys. Rev. Lett. 111, 3911 (13). 9 A. Harindranath, R. Kundu, A. Mukherjee, R. Rataole, Phys. Rev. Lett. 111, 391 (13). 1 A. Harindranath, R. Kundu and A. Mukherjee, Phys. Lett. B 78, 63 (14). 11 T. Gutsche, V. E. Lyuovitskij, I. Schmidt and A. Vega, Phys. Rev. D 89, 5433 (14). 1 C. Mondal and D. Chakraarti, arxiv: hep-ph. 13 S. J. Brodsky and G. F. de Teramond, Phys. Rev. D. 77, 567 (8). 14 O.V. Selyugin and O.V. Teryaev, Phys. Rev. D79, 333 (9). 7

28 15 Z. Aidin and C. E. Carlson, Phys. Rev. D 78, 715 (8). 16 D. Chakraarti, C. Mondal, Phys. Rev. D 88, 736 (13); D. Chakraarti, C. Mondal, Eur. Phys. J. C 73, 671 (13). 17 X. D. Ji, Phys. Rev. Lett. 78, 61 (1997); Phys. Rev. D 58, 563 (1998). 18 S. J. Brodsky, D. S. Hwang, Nucl. Phys. B543, 39 (1999); S. J. Brodsky, D. S. Hwang, B. -Q. Ma and I. Schmidt, Nucl. Phys. B593, 311 (1); S. J. Brodsky, M. Diehl and D. S. Hwang, Nucl. Phys. B596, 99 (1). 19 Y. Hatta, K. Tanaka, S. Yoshida, JHEP 13, 3 (13). E. Leader, Phys. Lett. B 7, 1 (13). 1 G. A. Miller, Phys. Rev. Lett. 99, 111 (7). C. E. Carlson and M. Vanderhaeghen, Phys. Rev. Lett. 1, 34 (8). 3 C. Granados, C. Weiss, JHEP 141, 9 (14). 4 D. Chakraarti and C. Mondal, Eur. Phys. J. C 74, 96 (14). 8