Isovector axial vector and pseudoscalar transition form factors of in QCD A. Kucukarslan, U. Ozdem, and A. Ozpineci Physics Department, Canakkale Onsekiz Mart University, 7 Canakkale, Turkey Physics Department, Middle East Technical University, 65 Ankara, Turkey (Dated: December, 4 We investigate the isovector axial vector and pseudoscalar form factors of baryon by employing light-cone QCD sum rules. Numerical calculations show that the form factors can be well fitted by the exponential form. We make a comparison with the predictions of lattice QCD, chiral perturbation theory and quark model. Keywords: Baryon axial and pseudoscalar form factors,, light-cone QCD sum rules arxiv:47.75v [hep-ph Dec 4 I. INTRODUCTION Form factors are important properties of hadrons. They describe how hadrons interact with each other and give information about the internal structure of the hadrons such as the measure of the charge and current distributions. The ( resonance is the lightest excitation of the nucleon. Solving the mysterious structure of the resonance has a significant effect to nuclear phenomenology. The is a rather broad resonance close to the πn threshold. Studying the structure of the resonance theoretically of using lattice QCD s an important input to phenomenology that cannot be directly extracted from experiments. This is because the decays strongly with a lifetime of seconds [, and resists experimental probes. Measurements of the magnetic moment exist despite a large experimental uncertainty. Moreover, the was studied within the framework of the lattice [. The decomposition of the ( matrix elements into the appropriate Lorentz invariant form factors was carried out and techniques to calculate form factors were developed and tested using quenched configurations. Besides, the axial charge g A is an important parameter for low-energy effective theories. g A can also be viewed as an indicator of the phenomenon of spontaneous breaking of chiral symmetry of nonperturbative QCD [4. Form factors are nonperturbative objects. Therefore, in order to study form factors, one needs to use a nonperturbative method. QCD sum rules is one of the nonperturbative methods. It is an analytical method which is directly based on the QCD Lagrangian. It is a powerful tool to extract hadron properties [5 8. Using this method static (e.g. baryon masses and dynamic (e.g. magnetic moments and coupling constants hadronic parameters can be determined. An alternative to the traditional QCD sum rules is the light cone QCD sum rules (LCSR. LCSR is especially suitable to study the interactions of hadrons at large momentum transfer [9. In LCSR, the hadronic properties are expressed in terms of the properties of the vacuum and the light cone distribution amplitudes of the hadrons in the process. Since the form factors are expressed in terms of the properties of the QCD vacuum and the distribution amplitudes, any uncertainty in these parameters reflects in the uncertainty of the predictions of the form factors. This method has been rather successful in determining hadron form factors at high Q (see e.g. [ 8. Using the LCSR, the isovector axial vector form factors of octet baryons have been calculated [7. For baryon isovector axial vector and pseudoscalar form factors have been studied using lattice QCD [, chiral perturbation theory [9 and quark models [4,. Our aim in this work is to study the axial vector and pseudoscalar form factors in the framework of light-cone QCD sum rules for the baryon. We give the formulation of the baryon form factors on the light cone and derive our sum rules, then we present our numerical results. In the last section, we conclude our work with a discussion on our results. II. THE AXIAL AND PSEUDOSCALAR FORM FACTORS In the LCSR approach, the procedure begins with the following two-point correlation function: Π µν (p,q = i d 4 xe iqx T[ηµ(A ν(x (p,s, ( akucukarslan@comu.edu.tr ozpineci@metu.edu.tr
where η (x is an interpolating current for the. Here we choose the interpolating current as follows [, η µ ( = ǫ abc [(u at (Cγ µ d b (u c (+(u at (Cγ µ u b (d c ( ( where a, b, c are the color indices and C denotes charge conjugation. The matrix elements of the axial-vector transition are parametrized in terms of four invariant form factors as follows [; (p,s A ν (x (p,s = i ( υα (p,s [g αβ g A (q γ ν γ 5 +g A (q q νγ 5 M + qα q β ( 4M h A (q γ ν γ 5 +h A (q q νγ 5 υ β (p,s ( M wherea ν (x = [(ū(xγ νγ 5 u(x d(xγ ν γ 5 d(x isthe isovector-axialvectorcurrent,q = p p, M isdeltamass, υ α is a Rarita-Schwinger spinor describing spin / fermions. The axial charge of the baryon is defined as g A = g A ( where the factor is due to the spin-/ nature of the baryon. Inthederivationofthe sumrules, summationsoverthespinsofthe baryonwillbeperformedusingthesummation formula for the Rarita-Schwinger spinors: υ µ (p,sυ ν (p,s = ( p +m {g µν γ µγ ν p µp ν m + p µγ ν p νγ µ m } (4 Note that, isoscalar axial vector current differs from the isovector axial vector current by the relative sign of the two terms. This current has an anomaly which would contribute to the form factors. For this reason, in this work we will not study the isoscalar current. The correlationfunction given in Eq. ( can be calculatedin terms ofhadronicpropertiesif p > and (p+q >, and also in terms of QCD parameters and several distribution amplitudes (DAs of the baryon in the deep Euclidean region p and (p+q. In order to construct the sum rules for the form factors, we calculate the correlation function in terms of hadron and quark-gluon parameters. Inserting the interpolating currents for the into the correlation function in Eq. (, we determine Π µν = i 6 d 4 xe iqx { (Cγ µ αβ (γ ν γ 5 ρσ 4ǫ abc a b c q σ (q θ (xq φ ( (p,s [δ ηαδ θσδ φβ S( x λρ +δ ηλ δθσδ φβ S( x αρ +δ ηαδ σθ δ φλ S( x βρ +δ ηβ δθσδ λφs( x αρ 4ǫ abc q a σ (q b θ (q c φ (x (p,s [ } δα η δθ λ δφ σ S( x βρ +δα η δθ β δφ σ S( x λρ (5 where q i (i =,, denote the quark fields. For the case of the baryon, q = q = u and q = d. Setting m u = m d =, the light quark propagator S(x in an external gluon field can be written as: S q (x = ix/ π x 4 q q (+ m x 6 ig s [ dυ x/ 6π x 4G µνσ µν υx µ G µν γ ν i 4π x. The first term of this expression describes the hard-quark propagator. The second term gives the contributions from the nonperturbative structure of the QCD vacuum, that is, the quark and quark-gluon condensates. After applying Borel transformation these terms do not contribute to the sum rules. The same holds for any other contribution that contains higher order condensates which are multiplied by positive powers of x. The last term is due to the correction in the background gluon field and gives rise to four-particle and five-particle baryon distribution amplitudes, which are not yet known. These distribution amplitudes contain contributions from operators of higher conformal spin than the ones used in the estimation of the distribution amplitudes [4. Hence, their inclusion would not be consistent with the use of the distribution amplitudes calculated in [5. Furthermore, in [6, it has been argued that for the nucleon, the contribution of the higher Fock states that will contribute to the four parton DAs is small compared to the valence quark DAs. It would be a reasonable approximation to ignore such higher Fock state contribution also for the baryon. In this work these contributions shall not be taken into account, which leaves us with only the first term in propagator to consider.
The pseudoscalar current is defined as P(x = ( ū(xγ 5 u(x d(xγ 5 d(x (6 The matrix elements of the pseudoscalar current is decomposed in terms of two form factors as follows [, 7; (p,s P(x (p,s = υ σ(p,s [g στ ( g P (q γ 5 + qσ q τ ( h P (q γ 5 υ τ (p,s (7 Then, inserting the interpolating current into the correlation function we determine; Π µν = i 6 d 4 xe iqx { (Cγ µ αβ (γ 5 ρσ 4ǫ abc a b c q σ (q θ (xq φ ( (p,s [δ ηαδ θσδ φβ S( x λρ +δ ηλ δθσδ φβ S( x αρ +δ ηαδ σθ δ φλ S( x βρ +δ ηβ δθσδ λφs( x αρ 4M 4ǫ abc q a σ (q b θ (q c φ (x (p,s [ } δα η δθ λ δφ σ S( x βρ +δα η δθ β δφ σ S( x λρ (8 The matrix elements of the local three-quark operator 4ǫ abc q a α (a xq b β (a xq c γ (a x (p,s (a,, are real numbers denoting the coordinates of the valence quarks 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. We refer the reader to Refs. [5 for a detailed analysis on DAs of, which we employ in our work to extract the axial-vector and pseudoscalar form factors. The long-distance side of the correlation function is obtained using the analyticity of the correlation function, then we can write the correlation function in terms of a dispersion relation of the form Π µ (p,q = π ImΠ µ(s (s p ds (9 The ground state hadron contribution is singled out by utilizing the zero-width approximation ImΠ µ = πδ(s m η (p (p (A ν (x,p(x (p +πρ h (s ( The matrix element of the interpolating current between the vacuum and baryon state is defined as η µ (p,s = λ υ µ (p,s where λ is the baryon overlap amplitude and υ µ (p,s is the baryon spinor. Physical part of the correlation function can be obtained as follows for axial vector current; [{ λ Π µν = i (M p +g A (q (M +4q 4M 5 h A (q ga (q M { 4 q µγ 5 ν +p ν γ 5 µ +γ ν q/γ 5 µ { (M + +q } 6M 4 h A (q + ga (q M q µ γ ν q/γ 5 q β β + ha (q q µ p ν q/γ 5 q β β ga (q { q ν q/γ 5 µ + M { + M ( M +5q M 4 for pseudoscalar current } q µ q ν q/γ 5 q β β } { g A (q +g A (q } q ν γ 5 µ q µ q/γ 5 ν M { (M 4q h A (q + ga (q 6M 4 (M +q M h A (q h A (q (5M +q 6M 4 h A (q 4gA (q M g (q M M } q µ γ ν γ 5 q β β }q µ q ν γ 5 q β β } q µ p ν γ 5 q β β (
4 λ [ g P (q M Π µν = (M γ 5 µ + h P (q q µ γ 5 q τ τ p 6 4M ( In these expressions, only the contribution of the spin-/ baryon is shown. In principle, the correlation function can also include the contribution from spin-/ particles. The overlap of the spin-/ particles with the η µ current can be written as /(p η µ = ( Ap µ +Bγ µ u(p ( where u(p is the spinor describing the spin-/ particle. Hence, when the gamma matrices are put into the order γ µ γ ν q p in the related correlation function, spin-/ states contribute only to the structures which have γ µ at the beginning or which are proportional to p µ. Then the contributions of the spin-/ states in the correlation function are eliminated by ignoring the structures proportional to p µ and the structures that contain a γ µ at the beginning. In this way, only the contributions from spin-/ states are kept [8, 9. The QCD sum rules are obtained by matching the short-distance calculation of the correlation function with the longdistance calculation. As can be seen from Eq. (, the coefficients of the structures p µ γ 5 ν and (q ν γ 5 µ p µ γ 5 ν are proportional to the form factors g A and g A respectively. Similarly, coefficient of the structure γ 5 µ in Eq. ( is proportional to the form factor g P. Hence, in this work we will use the above mentioned structures to determine the g A, g A and g P form factors. The following expressions for axial vector form factors are obtained: x [ g A (q λ M p = f M dx (q px d x 4V(x,x, x x x dx (q px d x [ T +A V(x, x x,x x [ g A λ (q M p = f M dx (q px d x [T +4A+8V(x,x, x x x + dx (q px d x [ T +A+V(x, x x,x, for pseudoscalar form factors (4 x g P (q λ M p = 6f [ M x dx (q px d x T(x,x, x x x x + dx (q px d x T(x, x x,x In order to eliminate the subtraction terms in the spectral representation of the correlation function, the Borel transformation is performed. After the transformation, contributions from excited and continuum states are also exponentially suppressed. Then, using the quark-hadron duality and subtracted, the contributions of the higher states and the continuum can be modeled. Both of the Borel transformation and the subtraction of higher states are achieved by using following substitution rules (see e.g. [5: (5 dx ρ(x (q xp x dx x ρ(xe s(x/m, dx ρ(x (q xp 4 M x dx x ρ(xe s(x/m + ρ(x Q +x m Be s/m, dt ρ(t (q tp dx 6 M 4 x x ρ(xe s(x/m + x Q +x m B [ M d ρ(x dx x (Q +x m B ρ(x x (Q +x m B
5 where m B = m, s(x = ( xm B + x x Q, M is the Borel mass and x is the solution of the quadratic equation for s = : x = [ (Q + m B +4m B Q (Q + m B /(m B, where is the continuum threshold. III. NUMERICAL ANALYSIS In this section, we present our numerical predictions of the axial vector and pseudoscalar form factors of. To obtain the numerical results, we need the expressions of the baryon DAs which are studied in [5. In this work, we obtain our results using the explicit expressions of DAs in Ref.[5. In order to calculate the form factors, we need to determine the value of the residue of the baryon, λ. We choose the value as λ =.8 GeV from analysis of the mass sum rules [,,,, and use the following numerical values of the parameters; M =. GeV [, f =.±. GeV [5. In the traditional analysis of sum rules, the spectral density of the higher states and the continuum are parametrized using quark hadron duality. In this approach, the spectral density corresponding to the contributions of the higher states and continuum is parametrized as ρ h (s = ρ QCD (sθ(s. The predictions for the form factors depend on two auxiliary parameters: the Borel mass M, and the continuum threshold. The continuum threshold signalsthe scale at which, the excited states and continuum start to contribute to the correlation function. Hence it is expected that (m +. GeV =.4 GeV. One approachto determine the continuumthresholdandthe workingregionoftheborelparameterm istoplot thedependence ofthe predictions on M for a range of values of the continuum threshold and determine the values of for which there is a stable region with respect to variations of the Borel parameter M. For this reason, in Figs., we plot the dependence of the form factors g A(Q, g A(Q and g P (Q on M for two fixed values of Q and for various values of in the range GeV < < 4 GeV. As can be seen from these figures, for =.5±.5 GeV, the predictions are practically independent of the value of M for the shown range. The uncertainty due to variations of in this range is much larger than the uncertainty due to variations with respect to M. Note that the determined range of is in the range that one would expect from the physical interpretation of. In Figs., we plot the form factors as a function of Q. The qualitative behavior of the form factors agree with our expectations. The form factors g A and gp are the dominant axial-vector, pseudoscalar form factor and the only one that can be extracted directly from the matrix element at Q =, determining the axial charge and pseudoscalar form factors of the. g A(P, is generally assumed to have a dipole form; g A(P (q = g A(P ((+Q /M A(P (6 Using this fit function we obtain M A =.85 GeV and M P =.74 GeV. Note that, in the vector meson dominance (VMD model, the pole of the form factors is given by the mass of the meson that couples to the current. The lightest axial vector meson has a mass of m A =. GeV [; therefore our results also are smaller from the predictions of the VMD model. We have also tried exponential form g A(P (q = g A(P (exp[ Q /M A(P (7 Our results are shown in Table I. The fit from. GeV Q GeV region produces the empirical value of g A quite well for the baryon. Besides we observe that the axial and pseudoscalar masses are very close to each other for all regions. In Table II we present the different numerical results predicted from other theoretical models. As seen from table, our result is slightly larger than the result obtained from Lattice QCD [, approximately two times smaller compared to the predictions of ChPT [9 and quark models [.
6 Fit Region (GeV g A ( M A (GeV [..6.5 [.5.88.4 [..7. Fit Region (GeV g P( M P (GeV [..7. [.5.5.6 [... TABLE I: The values of exponential fit parameters, g A(P and M A(P for axial and pseudoscalar form factors. The results include the fits from three region. [ [9 [ [ [ This Work g A.9±. 4.5 4.47 4.48 4..7±.6 TABLE II: Different results from theoretical models which are Lattice QCD [, ChPT [9, quark models [ and also our model. IV. CONCLUSION In this study, we have calculated the isovector axial vector and pseudoscalar form factors for baryon within the light cone QCD sum rules method. We extract axial g A, g A and pseudoscalar form factors g P. The form factors (h A,hA and h P cannot obtain from our approach since the necessary DAs have not been known yet. Finally, we have found that our results for the g A and g P form factors compare well with latticeqcd, chiral perturbation theory and quark models results in the high-q region and with those from other theoretical approaches when extrapolated to low-q regions via an exponential fit. The result obtained from the axial charge exponential fit different from other theoretical approaches. ACKNOWLEDGMENTS This work has been supported by The Scientific and Technological Research Council of Turkey (TUBITAK under Project No. T45. The work of A.K. and A. O. is also partially supported by the European Union (HadronPhysics project Study of strongly interacting matter. [ G. Lopez Castro and A. Mariano. Phys.Lett., B 57:9(. [ M. Kotulla, J. Ahrens, J.R.M. Annand, R. Beck, G. Caselotti, et al. Phys.Rev.Lett., 89:7,. [ C. Alexandrou, E.B. Gregory, T. Korzec, G. Koutsou, J.W. Negele, T. Sato, and A. Tsapalis. Phys.Rev., D87:45,. [4 Ki-Seok Choi, W. Plessas, and R.F. Wagenbrunn. Phys.Rev., D8:47,. [5 M. A. Shifman, A.I. Vainshtein, and Valentin I. Zakharov. Nucl.Phys., B47:85 447, 979. [6 M. A. Shifman, A.I. Vainshtein, and Valentin I. Zakharov. Nucl.Phys., B47:448 58, 979. [7 L.J. Reinders, H. Rubinstein, and S. Yazaki. Phys.Rept., 7:, 985. [8 B.L. Ioffe and Andrei V. Smilga. Nucl.Phys., B:9, 984. [9 V. M. Braun and I.E. Filyanov. Z. Phys., C44:57, 989. [ I.I. Balitsky, V. M. Braun, and A.V. Kolesnichenko. Nucl. Phys., B:59 55, 989. [ V.L. Chernyak and I.R. Zhitnitsky. Nucl. Phys., B45:7 7, 99. [ T.M. Aliev and A. Ozpineci. Nucl. Phys., B7:9, 6. [ T.M. Aliev and M. Savci. Phys.Lett., B656:56 66, 7. [4 Zhi-Gang Wang, Shao-Long Wan, and Wei-Min Yang. Eur.Phys.J., C47:75 84, 6. [5 V.M. Braun, A. Lenz, and M. Wittmann. Phys. Rev., D7:949, 6. [6 G. Erkol and A. Ozpineci. Phys.Lett., B74:55 558,.
[7 G. Erkol and A. Ozpineci. Phys.Rev., D8:4,. [8 T.M. Aliev, K. Azizi, and M. Savci. Phys.Rev., D84:765,. [9 Fu-Jiun Jiang and Brian C. Tiburzi. Phys.Rev., D78:754, 8. [ L. Theussl, R.F. Wagenbrunn, B. Desplanques, and Willibald Plessas. Eur.Phys.J., A:9,. [ L. Ya. Glozman, Willibald Plessas, K. Varga, and R.F. Wagenbrunn. Phys.Rev., D58:94, 998. [ K. Glantschnig, R. Kainhofer, Willibald Plessas, B. Sengl, and R.F. Wagenbrunn. Eur.Phys.J., A:57 55, 5. [ B.L. Ioffe. Nucl.Phys., B88:7 4, 98. [4 V. Braun, R.J. Fries, N. Mahnke, and E. Stein. Nucl. Phys., B589:8 49,. [5 Carl E. Carlson and J.L. Poor. Phys.Rev., D 8:758, 988. [6 M. Diehl, T. Feldmann, R. Jakob, and P. Kroll. Eur. Phys. J., C8:49 44, 999. [7 C. Alexandrou, E. B. Gregory, T. Korzec, G. Koutsou, J. W. Negele, T. Sato, and A. Tsapalis. Phys.Rev.Lett., 7:46,. [8 V.M. Belyaev. arxiv:hep-ph/957, 99. [9 V.M. Belyaev and B.L. Ioffe. Sov.Phys.JETP, 57:76 7, 98. [ Frank X. Lee. Phys. Rev., C57: 8, 998. [ W.Y.P. Hwang and Kwei-Chou Yang. Phys.Rev., D49:46 465, 994. [ J. Beringer et al. Phys.Rev., D86:,. 7
8 g (Q = GeV.8.7.6.5.4. =. GeV =. GeV =.5 GeV = 4. GeV g (Q = 4 GeV.4.5..5..5. =. GeV =. GeV =.5 GeV = 4. GeV..5. 4 6 8 M (GeV (a 4 6 8 M (GeV (b g (Q = GeV g ~ (Q = GeV.4..8.6.4..8.6.4....9.8.7.6.5.4.. =. GeV =. GeV =.5 GeV = 4. GeV 4 6 8 M (GeV =. GeV =. GeV =.5 GeV = 4. GeV (c 4 6 8 M (GeV (e g (Q = 4 GeV g ~ (Q = 4 GeV.7.6.54.45.6.7.8.9.5.4... =. GeV =. GeV =.5 GeV = 4. GeV 4 6 8 M (GeV =. GeV =. GeV =.5 GeV = 4. GeV (d 4 6 8 M (GeV (f FIG. : The dependence of the form factors; on the Borel parameter squared MB for the values of the continuum threshold = GeV, = GeV, =.5 GeV and = 4 GeV and Q = and 4 GeV, (a and (b for g A axial form factor, (c and (d for ga axial form factor, (e and (f for gp (Q pseudoscalar form factor.
9 g (Q.8.6.4 =. GeV, M =. GeV, M =. GeV =. GeV, M =. Ge =. GeV, M = 6. GeV, M = 6. GeV =. GeV, M = 6. GeV. 4 6 8 Q (GeV (a g (Q.5.5 =. GeV, M =. GeV, M =. GeV =. GeV, M =. GeV =. GeV, M = 6. GeV, M = 6. GeV =. GeV, M = 6. GeV.5.5 4 6 8 Q (GeV (b g ~ (Q.5 =. GeV, M =. GeV, M =. GeV =. GeV, M =. GeV =. GeV, M = 6. GeV, M = 6. GeV =. GeV, M = 6. GeV 4 6 8 Q (GeV (c FIG. : The dependence of the form factors for the values of the continuum threshold = GeV, = GeV, =.5 GeV and M = and 6 GeV (a for g A axial form factors, (b for ga axial form factor, (c for gp (Q pseudoscalar form factor.