The ρ-meson longitudinal leading-twist distribution amplitude Hai-Bing Fu, Xing-Gang Wu, Hua-Yong Han, Yang Ma, and Huan-Yu Bi Department of Physics, Chongqing University, Chongqing 433, P.R. China arxiv:49.353v [hep-ph] 5 Sep 4 In the present paper, we suggest a convenient model for the vector ρ-meson longitudinal leadingtwist distribution amplitude φ ;ρ, whose distribution is controlled by a single parameter B ;ρ. By choosing proper chiral current in the correlator, we obtain new light-cone sum rules (LCSR) for the B ρ TFFs A, A and V, in which the δ -order φ ;ρ provides dominant contributions. Then we make a detailed discussion on the φ ;ρ properties via those B ρ TFFs. A proper choice of can make all the TFFs agree with the lattice QCD predictions. A prediction of V ub has also been presented by using the extrapolated TFFs, which indicates that a larger leads to a larger V ub. To compare with the BABAR data on V ub, the longitudinal leading-twist DA φ ;ρ prefers a doubly-humped behavior. PACS numbers: 3.5.Hw,.55.Hx,.38.Aw, 4.4.Be I. INTRODUCTION The meson light-cone distribution amplitude (LCDA) relates to the matrix elements of the nonlocal light-ray operators sandwiched between the hadronic state and the vacuum. It provides underlying links between the hadronic phenomena at the small and large distances. The leading-twist (twist-) LCDA arouses people s great interests e to its role in hard exclusive processes [], i.e. it always provides dominant contributions to those processes. The investigation on the LCDAs has been the subject of numerous studies, especially for the simpler LCDAs of the pseudoscalar mesons pion, kaon and etc. As for the vector ρ meson, its LCDAs are much more involved e to its complex structures: there are chiraleven or chiral-odd LCDAs for the ρ meson, the ρ meson has longitudinal ( ) and transverse ( ) polarization states. Several approaches have been adopted to study the ρ-meson LCDA properties, such as the QCD sum rules [ 5], the light-cone quark model [6], the AdS/QCD model [7] and etc. The properties of the ρ-meson LCDA are also helpful for us to understand other vector mesons LCDAs. That is, by taking the SU f (3)-breaking effect into consideration, one can further study the K meson LCDAs [6, 8, 9]. twist δ δ δ δ 3 φ ;ρ φ ;ρ / / 3 / φ 3;ρ,ψ 3;ρ,Φ 3;ρ, Φ 3;ρ φ 3;ρ,ψ 3;ρ,Φ 3;ρ / 4 / / φ 4;ρ,ψ 4;ρ,Ψ 4;ρ, Ψ 4;ρ φ 4;ρ,ψ 4;ρ TABLE I. The ρ-meson LCDAs up to twist-4 and δ 3 -order. To understand the complex structures of the ρ-meson LCDAs, as suggested in Refs.[, ], it is convenient to arrange them by a parameter δ, i.e. δ m ρ /m b 6%. Email: wuxg@cqu.e.cn A classification of the ρ-meson twist-, twist-3 and twist- 4 LCDAs up to δ 3 -order have been collected in Ref. [], which are rearranged in Table I. Up to twist-4, there are fifteen ρ-meson LCDAs. All those ρ-meson LCDAs, especially the higher-twist DAs, are far from affirmation, thus, it is helpful to find a reliable way to study the properties of a specific ρ-meson LCDA. In the paper, we shall concentrate our attention on the ρ-meson longitudinal leading-twist distribution amplitude φ ;ρ, which is at the δ -order. Forthepurpose, weadopttheb ρtransition form factors (TFFs) to study the φ ;ρ properties. The B ρ TFFs are key components for the semileptonic decay B ρlν. When the leptons are massless, only three TFFs A, A and V provide non-zero contributions to the B ρlν decay [3]. In the following, we shall adopt the QCD light-cone sum rules(lcsr)[4 6] to deal with those non-zero B ρ TFFs. In comparison to the convention SVZ sum rules [7], the LCSR is based on the operator proct expansion near the light cone, and all its non-perturbative dynamics are parameterized into LCDAs of increasing twists. One advantage of LCSR lies in that it allows to incorporate information about high-energy asymptotics of correlation functions in QCD, which is accumulated in the LCDAs. Furthermore, as a tricky point of the LCSR approach, by proper choosing the correlation function (correlator), especially by choosing the chiral correlator, one can highlight the wanted LCDAs contributions while highly suppress the unwanted LCDAs contributions to the LCSR [8 ]. As shall be shown later, by taking a proper chiral correlator, we can highlight the contributions of φ ;ρ to the B ρ TFFs. The obtained LCSRs for those TFFs shall show strong dependence on φ ;ρ. Thus, via the comparison with the data or predictions from other approaches, it shall provides us good opportunities to determine the behavior of φ ;ρ. The remaining parts of the paper are organized as follows. In Sec. II, we present the formulas for the B ρ TFFs under the LCSR approach. In Sec. III,
we present our numerical results and discussions on the B ρ TFFs, the B ρ semi-leptonic decay width and the Cabibbo-Kobayashi-Maskawa(CKM) matrix element V ub. Sec. IV is reserved for a summary. II. CALCULATION TECHNOLOGY To derive the LCSRs for the B ρ TFFs, we suggest to start from the following chiral correlator Π µ (p,q) = } i d 4 xe iq x ρ(p,λ) T { q (x)γ µ ( γ 5 )b(x),j B (), () where the current j B (x) = i b(x)( γ 5 )q (x). This is different from our previous choice of j B (x) = i b(x)( + γ 5 )q (x)[],inwhichtheδ -orderφ ;ρ providedominant contribution. The advantage of our present new choice lies in that one can highlight the contributions from the chiral-even ρ-meson LCDAs such as φ ;ρ, φ 3;ρ, ψ 3;ρ, Φ 3;ρ, Φ 3;ρ, φ 4;ρ and ψ 4;ρ ; while the chiral-odd LCDAs φ ;ρ, φ 3;ρ, ψ 3;ρ, Φ 3;ρ, φ 4;ρ, ψ 4;ρ, Ψ 4;ρ, Ψ 4;ρ provide zero contributions. We also find that within the remaining chiraleven LCDAs, only the LCDAs φ ;ρ, φ 3;ρ and ψ 3;ρ provide dominant contributions to the LCSR, while the remaining LCDAs as φ 4;ρ, ψ 4;ρ, Φ 3;ρ and Φ 3;ρ provide negligible contribution. Moreover, the twist-3 DAs ψ3;ρ and φ 3;ρ (andthelatersimplifieda ρ(u) = u dv[φ ;ρ (v) φ 3;ρ(v)]) can be related to φ ;ρ under the Wandzura-Wilczek approximation [3, ] φ 3;ρ(x,µ) = [ Φ x ρ (x,µ)+φ x ρ (x,µ) ] /, () ψ3;ρ x (x,µ) = [ xφ ρ (x,µ)+xφ x ρ (x,µ) ], (3) A ρ x (x,µ) = [ xφ ρ (x,µ)+xφ x ρ (x,µ) ] /, (4) where Φ x ρ (x,µ) = x v and Φρ x (x,µ) = x dvφ ;ρ (v,µ) v. Thus,theleading-twistφ ;ρ providesdominant contribution to the LCSRs either directly or indirectly. Then, as required, those TFFs can indeed provide us a useful platform for testing the properties of φ ;ρ. dvφ ;ρ (v,µ) The correlator () is an analytic function of q defined at both the space-like and the time-like q -regions. The correlator can be treated by inserting a completed set of intermediate hadronic states in physical region. It can also be treated in the framework of the operator proct expansion in deep Euclidean region, then all its nonperturbative dynamics are parameterized into LCDAs. Those two results can be related by the dispersion relation, and the final LCSR can be achieved by applying the Borel transformation. Following the standard LCSR proceres, we obtain the LCSRs for the B ρ TFFs f B A (q )e m B /M = m b m ρf ρ { [ ] m B (m Θ(c(u,s B +m ρ ) u e s(u)/m ))φ 3;ρ(u) m ρ um Θ(c(u,s ))C ρ(u) m ρ Dα dve } s(x)/m X M Θ(c(X,s ))[Φ 3;ρ (α)+ Φ 3;ρ (α)], (5) f B A (q )e m B /M = m b m ρ(m B +m ρ )f ρ { [ m B u e s(u)/m um Θ(c(u,s ))A ρ (u)+ m ρ um 4 Θ(c(u,s ))C ρ (u) ] } + m b m ρ 4u 4 M 6 Θ(c(u,s ))B ρ(u) f B V(q )e m B /M = m b m ρ(m B +m ρ )f ρ m B +m ρ Dα dve s(x)/m X 3 M 4Θ(c(X,s ))[Φ 3;ρ ( α)+ Φ 3;ρ ( α)],(6) e s(u)/m u M Θ(c(u,s ))ψ 3;ρ(u), (7) where s( ) = [m b (q m ρ)]/ with = ( standsforuorx),andx = a +va 3. f ρ representstheρmeson decay constant. c(u,s ) = us m b +ūq uūm ρ. Θ(c(,s )) is the usual step function: when c(,s ) <, it is zero; otherwise, it is. Θ(c(u,s )) and Θ(c(u,s )) are defined via the integration u M e s(u)/m Θ(c(u,s ))f(u) = u 3 M 4e s(u)/m Θ(c(u,s ))f(u) = u u u f(u)+δ(c(u,s )), (8) M e s(u)/m u 3 M 4e s(u)/m f(u)+ (c(u,s )), (9)
3 where δ(c(u,s )) = e s/m f(u ) C and [ f(u (c(u,s )) = e s/m ) u M u d C C ( ) ] f(u), uc u=u C = m b + u m ρ q and u [,] is the solution of c(u,s ) =. The simplified DAs B ρ(u) = u dvφ 4;ρ (v) and C ρ(u) = u dv v dw [ψ 4;ρ (w)+φ ;ρ (w) φ 3;ρ (w)]. III. NUMERICAL RESULTS A. Input parameters A ;ρ (GeV ) b ;ρ (GeV ) a (GeV) a (.GeV). 8.96.643.8.5. 7.88.68.8.67. 6.4.64 +.6 +. +. 4.5.57 +.4 +.8 +..9.537 +.58 +.7 QCD evolution []. Table II shows a. Another way of constructing the LCWF has been suggested under the AdS/QCD-DA theory [4], ( ψ ;ρ (x,ζ) = N κ exp κ ζ ) ( ) exp m f π κ, which leads to the AdS/QCD-DA model [5] φ ;ρ (x,µ ) = 3m f πf ρ () dζµ J (µ ζ) ψ ;ρ (x,ζ), (3) where κ = m ρ /, m f =.4GeV [6 8]. N can be fixed by the normalization condition, e.g. we have N µ=gev =. and N µ=.gev =.9. TABLE II. The parameters for the ρ-meson leading-twist LCDA ψ ;ρ, and its corresponding second Gegenbauer moment a at GeV and. GeV, respectively..5.5 φ ;ρ (x, µ ) ψ3;ρ (x, µ ) The leading-twist LCDA φ ;ρ can be derived from its light-cone wavefunction (LCWF), i.e. φ ;ρ (x,µ ) = 3 dk f ρ 6π 3ψ ;ρ (x,k ), () k µ where f ρ = f ρ/ 5. One way of constructing the LCWF model has been suggested in by Wu and Huang []. Following the idea, the WH-DA model for φ ;ρ is φ ;ρ (x,µ ) = A ;ρ 3mq 8π 3/ f [+B ρb ;ρ C3/ (ξ)] ;ρ Erf b ;ρ µ +m q Erf ( b ;ρ ) m q,() where µ is the factorization scale, the error function Erf(x) = x e t dt/ π and the constitute light-quark mass m q 3 MeV. Apart from the normalization condition, we take the average of the transverse momentum.37 GeV [3] as another constraint. The DA k / ;ρ parameters for [.,.] are presented in Table II, in which its second Gegenbauer moment a under two typical scales, i.e. µ =.GeV and.gev, are also presented. The value of a at.gev is obtained by.5....4.6.8 A ρ(x, µ )...4.6.8.5.5.5 x WH-DA for [.,.] BB-DA AdS/QCD-DA..4.6.8 WH-DA for =. φ 3;ρ (x, µ ) WH-DA for =...4.6.8 FIG.. A comparison of φ ;ρ (x,µ), φ 3;ρ(x,µ ), ψ 3;ρ(x,µ ) and A ρ(x,µ ) under various models. The shaded band are for the WH-DA with [.,.]. The solid and dashed lines are for BB-DA and AdS/QCD-DA, respectively. In Fig., we present a comparison of various DA models, i.e. the WH-DA with [.,.], the AdS/QCD-DA and the BB-DA. Here the BB-DA stands for the usual Gegenbauer DA expansion with its second Gegenbauer moment a (GeV) =.5(7) as suggested by Ball and Braun [9]. We also present the twist-3 DAs in Fig., where the WH-model is derived from the relations and the AdS/QCD and BB ones are from Ref.[7] and Ref.[9], respectively. Fig. shows that by varying [.,.], the WH-DA shall vary from the single-peaked behavior to
4 the doubly-humped behavior. It is noted that by setting =.4, the WH-DA behaves close to the AdS/QCD-DA; and by setting =., the WH-DA behaves close to the BB-DA. Thus the WH-DA provides a convenient form to mimic the behaviors of various DAs suggested in the literature. Since the LCDA φ ;ρ dominantly determines the B ρ TFFs, inversely, a definite TFFs shall be helpful for fixing the properties of φ ;ρ. In the following, we shall take the WH-DA to study the B ρ semi-leptonic decays. To do the numerical calculation, we take f ρ =.6±.3GeV [9], f B =.6±.9GeV [], and the b-quark pole mass m b = 4.8 ±.5GeV. The ρ and B-meson masses are taken as m ρ =.775GeV and m B = 5.79GeV [3]. As a cross check of the LCSRs (5,6,7), if taking the second Gegenbauer moment a of φ ;ρ as that of Refs.[, 3], we obtain consistent TFFs with those of Refs.[, 3]... +. s A 3.5(5) 3.(5) 3.5(5) M A.4(3).7(3) 3.(3) s A 3.5(5) 3.(5) 3.5(5) M A 3.(3).8(3).6(3) s V 3.5(5) 3.(5) 3.5(5) M V 3.(3) 3.3(3) 3.5(3) TABLE III. The continuum threshold s and the Borel parameter M under various DA models. The central values are for m b = 4.8 GeV and f B =.6 GeV. To set the parameters, such as the Borel window and the continuum threshold s, for the B ρ TFFs LCSR, we adopt the criteria: I) We require the continuum contribution to be less than 3% of the total LCSR. II) We require all the high-twist DAs contributions to be less than 5% of the total LCSR. III) The derivatives of Eqs.(5,6,7) with respect to ( /M ) give three LCSRs for the B-meson mass m B. And we require their predicted B-meson mass to be fulfilled in comparing with the experiment one, e.g. m th B mexp B /mexp B less than.%. We present the values of s and M for the B ρ TFFs in Table III. B. B ρ TFFs and the B ρ semi-leptonic decay By using the chirallcsrs (5,6,7), we present the B ρ TFFs under the WH-DA in Table IV, in which typical values for =.,.,. are adopted. The scale µ It is noted that Refs.[, 3] adopt the usual correlator to derive the LCSR, in which more LCDAs that our present LCSR have to be taken into consideration. A () A () V() =..9+.49.33.3+.33.4.53+.3. =..+.35.4.4+.4.9.58+.7. = +..+.3.7.+.54.37.64+.4.8 TABLE IV. The B ρ TFFs A, A and V at the large recoil region q =. The errors are squared averages of all the mentioned source of errors. is set tobe the typical momentum transferofthe process, µ (m B m b )/.GeV. Table IV shows that the TFFs decreases with the increment of or a. To show how the TFFs A, A and V are affected by the LCDAs, we present the TFFs at the maximum recoil region q = in Table V. Up to twist-4 and δ 3 -order accuracy, we show that only the δ -order LCDAs provide non-zero contributions. And among those δ -order LCDAs, A ρ provides the dominant contribution to A, and ψ 3;ρ provides the solitary contribution to V. By using the the 3-particle DAs Φ 3;ρ and Φ 3;ρ suggested in Refs.[, 3], we show they provide less than.% contributionstoa anda, whichexplainwhythoselcdas have been neglected in many references. As a further step, we can apply the B ρ TFFs to study the properties of the semileptonic decay B ρlν. The LCSRs is applicable in the region q 4GeV, while the physical allowable region for the TFFs is (m B m ρ ) GeV. Thus, certain extrapolation is needed. We adopt the formula F i (q ) = F i ()/[ a i q /m B +b i(q /m B ) ] to do the extrapolation, in which F,,3 stand for A, A and V, respectively. The parameters a i and b i are determined by requiring the quality <, which is defined as [3] t Fi (t) Fi fit (t) = t F, (4) i(t) where t [, 7,,,4]GeV. The fitted parameters for the extrapolation are put in Table VI. We put the extrapolated B ρ TFFs A (q ), A (q ) and V(q ) in Fig., where the lattice QCD predictions [33, 34] are presented for a comparison. It shows that by varying [.,.], the LCSRs predictions agree with the lattice QCD predictions and it shall increasewith the increment of. One can get astrong constraint on the ρ-meson DAs when more precise lattice QCD results are given. The differential decay width for the semileptonic decay B ρlν can be found in Ref.[3]. One can cut off the uncertainty from V ub by calculating the differential decay width / V ub dγ/dq, which is presented in Fig. 3. Moreover, the total decay width (Γ) can be separated as Γ + Γ, where Γ (Γ ) stands for the decay width of the ρ-meson longitudinal (transverse) components. We present the total decay width Γ and the ratio Γ /Γ in Table VII, where the errors in Table VII are squared av-
5 =. B ;ρ =. B ;ρ =. A () A () V() A () A () V() A () A () V() ψ3;ρ / /.53 / /.58 / /.64 φ 3;ρ.3 / /. / /. / / A ρ /. / /.5 / /. / C ρ.8.9 /.8. /.7. / Φ 3;ρ, Φ 3;ρ..7 /..8 /.. / Total.9.3.53..4.58...64 TABLE V. The B ρ TFFs A, A and V at the large recoil region q =, in which the WH-DA with =.,. and., has been adopted. F i a i b i A.36.7.7 =. A.86.956.6 V.97.97.7 A.949.6.3 =. A.684.748.7 V.998.969. A.859..8 =. A.499.55.4 V.899.843. TABLE VI. The fitted parameters a i and b i for the B ρ TFFs. is a measure of the quality of extrapolation. =. B ;ρ =. B ;ρ =. B ;ρ =. Γ/ V ub 4.69 +.34 5.57. +6.4 4.4.99 +4.9 3.39.5 +4.5 3. Γ /Γ.643 +.49.375.68+.348.363.597+.3.36.595+.37.37 TABLE VII. Total decay width Γ/ V ub and the ratio Γ /Γ under the WH-DA model. erage of the mentioned source of errors. C. The CKM matrix element V ub We take two types of semi-leptonic decays as a try to determine theckm matrixelement V ub. Oneisthe B - type via the process B ρ l + ν l with branching ratio and lifetime B(B ρ l + ν l ) = (.34±.8) 4 and τ(b ) =.59 ±.7ps [3]. Another is the B + -type via the process B + ρ l + ν l with branching ratio and lifetime B(B + ρ l + ν l ) = (.7 ±.3) 4 and τ(b + ) =.64 ±.8ps [3]. Our predicted weighted averages of V ub are given in Table VIII, the BABAR predictions are also included as a comparison. The errors come from the choices of the b-quark mass, the Borel window and the threshold parameter s, the experimental uncertainties are from the measured lifetimes V ub =..5±.38 =..76±.36 = +..9±.35 = +. 3.±.35 LCSR [].75 ±.4 BABAR [35] ISGW [37].83±.4 LCSR [].85 ±.4 BABAR [36] ISGW [37].9±.4 TABLE VIII. The weighted average of V ub in unit 3 from both the B + -type and B -type and for the WH-DA with =.,.,. and., respectively. The estimations of the BABAR collaboration [35, 36] are also presented as a comparison. and decay ratios. Table VIII show that V ub increases with the increment of. We observe that by choosing [.,.], its central value is consistent with the BABAR prediction[35, 36]. The BABAR predictions are based on the LCSR [] or ISGW [37] predictions on B ρ TFFs, respectively. IV. SUMMARY We have constructed a convenient model () for the longitudinal leading-twist LCDA φ ;ρ, in which a single parameter a controls its longitudinal behavior. By varying [.,.], the φ ;ρ shall evolve from a single-peaked behavior to a doubly-humped behavior. We have discussed its properties via the B ρ TFFs by using the LCSR approach. Those chiral LCSRs for the TFFs shall be dominated (over 99%) by φ ;ρ either directly or indirectly, then those TFFs do provide us a platform to test the properties of the leading-twist φ ;ρ. After extrapolating the TFFs to their allowable region, we compare them to those of lattice QCD predictions. Those TFFs become smaller with the incre-
6 A(q ) A(q ) V (q ).9.8.7.6.5.4.3..8.6.4..8.6.4. 4.5 4 3.5 3.5.5.5 WH-DA (B ;ρ =.) WH-DA (B ;ρ =.) WH-DA (B ;ρ =.) Lattice (β = 6.) [34] Lattice (β = 6.) [34] Lattice [33] 4 6 8 4 6 8 WH-DA (B ;ρ =.) WH-DA (B ;ρ =.) WH-DA (B ;ρ =.) Lattice (β = 6.) [34] Lattice (β = 6.) [34] Lattice [33] q (GeV ) 4 6 8 4 6 8 WH-DA (B ;ρ =.) WH-DA (B ;ρ =.) WH-DA (B ;ρ =.) Lattice (β = 6.) [34] Lattice (β = 6.) [34] Lattice [33] q (GeV ) 4 6 8 4 6 8 q (GeV ) FIG.. The extrapolated A (q ), A (q ) and V(q ) underthe WH-DA model, where the lattice QCD estimations [33, 34] are included as a comparison. ment of, and either a too smaller or a too larger is not allowed by the lattice QCD estimations. If we have a precise lattice QCD estimation, we can get a more strong constraint on the ρ-meson DA behavior. In comparison to the BABAR prediction on the differential decay width / V ub dγ/dq [35] and also the TFFs from Lattice QCD, we find that smaller B < ;ρ. (or a ;ρ (GeV) <.8) is not suitable. As a further / Vub dγ/dq ( GeV ).8.6.4..8.6.4. LQCD[33] LQCD(β = 6.)[34] LQCD(β = 6.)[34] WH-DA(B ;ρ =.) WH-DA(B ;ρ =.) WH-DA(B ;ρ =.) 4 6 8 4 6 8 q (GeV ) FIG. 3. Differential decay width / V ub dγ/dq for the WH-DA model. The lattice QCD estimations [33, 34] are included as a comparison. comparison with the BABAR prediction on the V ub, we observe that by choosing [.,.], its central values show a better agreement with the BABAR prediction [35, 36]. Thus, we can predict that φ ;ρ prefers a doubly-humped behavior other than the sing-peaked behavior of the transverse leading-twist LCDA φ ;ρ as suggested by Ref.[]. As a final remark, the present approach can also be applied to study other vector meson LCDAs. For example, by taking the SU f (3)-breaking effect into consideration, one can further study the K meson LCDAs. A detailed discussion on the K meson LCDAs under the same approach of the present paper is in preparation. Acknowledgments: This work was supported in part by Natural Science Foundation of China under Grant No.755 and No.758, and by the Fundamental Research Funds for the Central Universities under Grant No.CQDXWL--Z. [] G.P. Lepage and S.J. Brodsky, Phys. Rev. D, 57 (98); S.J. Brodsky and G.P. Lepage, Adv. Ser. Direct. High Energy Phys. 5, 93 (989). [] V.L. Chernyak and A.R. Zhitnitsky, Phys. Rept., 73 (984). [3] A.R. Zhitnitsky, I.R. Zhitnitsky and V.L. Chernyak, Sov. J. Nucl. Phys. 38, 775 (983). [4] P. Ball and V.M. Braun, Phys. Rev. D 54, 8 (996). [5] A.V. Pimikov, S.V. Mikhailov and N.G. Stefanis, Few Body Syst. 55, 4 (4). [6] H.M. Choi, C.R. Ji, Phys. Rev. D 75, 349 (7). [7] M. Ahmady and R. Sandapen, Phys. Rev. D 87, 543
7 (3). [8] P. Ball andm. Boglione, Phys. Rev.D68, 946 (3). [9] P. Ball and R. Zwicky, Phys. Lett. B 633, 89 (6). [] P. Ball and R. Zwicky, Phys. Rev. D 7, 49 (5). [] P. Ball, V.M. Braun, Y. Koike and K. Tanaka, Nucl. Phys. B 59, 33 (998); P. Ball and V.M. Braun, Nucl. Phys. B 543, (999). [] H.B.Fu, X.G.Wu, H.Y.HanandY.Ma, arxiv:46.389 (4). [3] P. Ball and V.M. Braun, Phys. Rev. D 55, 556 (997). [4] L.I. Balitsky, V.M. Braun and A.V. Kolesnichenko, Nucl. Phys. B 3, 59 (989); V.M. Braun and I.E. Filyanov, Z. Phys. C 44, 57 (989). [5] V.L. Chernyak and A.R. Zhitnitsky, Nucl. Phys. B 345, 37 (99). [6] P. Ball, V.M. Braun and H.G. Dosch, Phys. Rev. D 44, 3567 (99) [7] M.A. Shifman, A.I. Vainshtein, and V.I. Zakharov, Nucl. Phys. B47 385 (979); Nucl. Phys. B47, 448 (979). [8] T. Huang and Z.H. Li, Phys. Rev. D 57, 993 (998); T. Huang, Z.H. Li and X.Y. Wu, Phys. Rev. D 63, 94 (). [9] Z.G. Wang, M.Z. Zhou and T. Huang, Phys. Rev. D 67, 946 (3). [] X.G. Wu, T. Huang and Z.Y. Fang, Phys. Rev. D 77, 74 (8); X.G. Wu, Eur. Phys. J. C 57, 665 (8); X.G. Wu and T. Huang, Phys. Rev. D 79, 343 (9). [] S. Wandzura and F. Wilczek, Phys. Lett. B 7, 95 (977). [] X.G. Wu and T. Huang, Phys. Rev. D 8, 344 (). [3] X.H. Guo and T. Huang, Phys. Rev. D 43, 93 (99). [4] S.J. Brodsky and G.F. de Teramond, Phys. Rev. D 77, 567 (8); arxiv:8.54 [5] J.R. Forshaw and R. Sandapen, Phys. Rev. Lett. 9, 86 (). [6] G. Soyez, Phys. Lett. B 655, 3 (7). [7] J.R. Forshaw, R. Sandapen, and G. Shaw, JHEP 6, 5 (6). [8] J.R. Forshaw and G. Shaw, JHEP 4, 5 (4). [9] P. Ball, V.M. Braun, A. Lenz, JHEP 78, 9 (7). [3] P. Ball and V.M. Braun, Phys. Rev.D58, 946 (998). [3] J. Beringer et al., Particle Data Group, Phys. Rev. D 86, (). [3] P. Ball and R. Zwicky, Phys. Rev. D 7, 45 (5); Phys. Lett. B 65, 5 (5). [33] J.M. Flynn et.al. (UKQCD Coll.), Nucl. Phys. B 46, 37 (996). [34] K.C. Bowler, J.F. Gill, C.M. Maynard, J.M. Flynn (UKQCD coll.), JHEP 45, 35 (4). [35] P. del Amo Sanchez, et.al (BABAR Collaboration) Phys. Rev. D 83, 37 (). [36] B. Aubert, et.al (BABAR Collaboration) Phys. Rev. D 7, 5(5). [37] D. Scora and N. Isgur, Phys. Rev. D 5, 783 (995).