arxiv: v1 [hep-ph] 12 Jul 2016

Transcription

1 Analysis of the strong coupling form factors of Σ b NB and Σ c ND in QCD sum rules uo-liang Yu, Zhi-ang Wang, and Zhen-Yu Li 2 Department of Mathematics and Physics, North China Electric power university, Baoding 73, People s Republic of China 2 School of Physics and Electronic Science, uizhou Normal College, arxiv:68.346v [hep-ph] 2 Jul 26 uiyang 558, People s Republic of China (Dated: August 2, 26 In this article, we study the strong interaction of the vertexes Σ b NB and Σ cnd using the threepoint QCD sum rules under two different dirac structures. Considering the contributions of the vacuum condensates up to dimension 5 in the operation product expansion, the form factors of these vertexes are calculated. Then, we fit the form factors into analytical functions and extrapolate them into time-like regions, which giving the coupling constant. Our analysis indicates that the coupling constant for these two vertexes are Σb NB =.43±.eV and ΣcND = 3.76±.5eV. PACS numbers: 3.25.Ft; 4.4.Lb Introduction By this time, many heavy baryons have been observed by BaBar, Belle and CLEO Collaborations[ 4] such as + 2 antitriplet states(λ + c,ξ + c,ξ c, the 2+ and 3+ 2 sextet states(ωc,σ c,σ c and (Ω c,σ c,σ c [2]. Besides, several S-wave bottom baryon states such as Λ b,σ b,σ b,ξ b and Ω b have also been observed by CDF and LHCb Collaborations[5, 6]. The SELEX collaboration even have reported the observation of the signal for the doubly charmed baryon state Ξ + cc[7, 8]. Since then, people have showed great interest in studying the properties of these heavy baryons which contains at least a heavy quark[9 3]. The charm and bottom baryon states which contain one (two heavy quark(s are particularly interesting for studying dynamics of the light quarks in the presence of the heavy quark(s, and serve as an excellent ground for testing predictions of the quark models and heavy quark symmetry. The properties of the heavy baryons such as mass spectrum, radiative and strong decays have been studied by many researches, which is very important for us to further understand the heavy flavor physics[4 2]. In this regards, the strong coupling constants associating with heavy baryons play an important role in describing the strong interaction among the heavy baryons and other participated hadrons. In addition, the properties of B and D mesons in nuclear medium are closely related with Electronic address: yuguoliang2@63.com Electronic address: zgwang@aliyun.com

2 2 their interactions with the nucleons[22 24], i.e. D +p or n Λ + c,σ + c or Σ c B +p or n Λ b or Σ b From these processes, we can see that it is significant to know the values of the related strong couplingconstants Σb NB and ΣcND which isessentialto determinethe modificationson the masses, decay constants and other parameters of the B and D mesons in nuclear medium. Up to now, only a few works on the strong coupling constants of the heavy baryons with the nucleon and heavy mesons have been reported[9, 2, 25, 26]. On the other hand, QCD sum rules is one of the most powerful non-perturbative methods, which is also independent of model parameters. In recent years, numerous research articles have been reported about the precise determination of the strong form factors and coupling constants via QCDSR[27 4, 44 46]. In this work, we use the QCDSR formalism to obtain the coupling constants of the strong vertexs Σ b NB and Σ c ND, where the contributions of vacuum condensates up to 5 in the OPE are considered. Theoutlineofthispaperisasfollows. InSect.2, westudythestrongvertexsσ b NB andσ c ND using the three-point QCDSR under two different dirac structures p/γ 5 and q/γ 5. Besides of the perturbative contribution, we also consider the contributions of qq, and qq at OPE side. In Sect.3, we present the numerical results and discussions, and Sect.4 is reserved for our conclusions. 2 QCD sum rules for Σ b NB and Σ c ND The three-point correlation functions of these two vertices Σ b NB and Σ c ND can be written as: Π(p,p,q = i 2 d 4 x d 4 ye ip.x e ip.y τ(j N (yj B[D] (J Σb [Σ c](x ( Where τ is the time ordered product and J Σb [Σ c](x,j N (y and J B[D] ( are the interpolating currents of the hadrons Σ b [Σ c ], N and B[D] respectively: J Σb [ c ] (x = ǫ ijk (u it (xcγ µ d j (x γ 5 γ µ b[c] k (x (2 J N (y = ǫ ijk (u it (ycγ µ u j (y γ 5 γ µ d k (y (3 J B[D] ( = u(γ 5 b[c]( (4 where C is the charge conjugation operator, and i, j and k are color indices. According to the QCD sum rules, the three-point correlation function can be calculated in two different ways. In the first way, the calculation is carried out in hadron degrees of freedom, called the phenomenological side. Secondly, it is called OPE side which is calculated in quark degrees of freedom. Then, invoking the quark-hadron duality, we equate the phenomenological and OPE sides from which the QCD sum rules for the strong coupling form factors is attained.

3 3 2. The phenomenological side We insert a complete set of intermediate hadronic states with the same quantum numbers as the operators J Σb [Σ c](x,j N (y and J B[D] ( into the correlation function Eq( to obtain the phenomenological representations. After isolating the ground-state contributions, the correlation function is written as: Π HAD (p,p,q = J N N(p J B[D] B[D](q Σ b [Σ c ](p J Σb [Σ c] (p 2 m 2 Σ b [Σ c] (p 2 m 2 N (q2 m 2 B[D] N(p B[D](q Σ b [Σ c ](p + (5 Where h.r. stands for the contributions of higher resonances and continuum states. And the matrix elements appearing in the above equation can be parameterized as the following formulas: J N N(p = λ N u N (p,s (6 J B[D] B[D](q = i m2 B[D] f B[D] m u +m b[c] (7 Σ b [Σ c ](p J Σb [Σ c] = λ Σb [Σ c]u Σb [Σ c](p,s (8 N(p B[D](q Σ b [Σ c ](p = Σb NB[Σ cnd]u N (p,s iγ 5 u Σb [Σ c](p,s (9 Where λ N and λ Σb [Σ b ] are residues of N and Σ b [Σ b ] baryons, f B[D] is the leptonic decay constant of B[D] meson and Σb NB[Σ cnd] is the strong coupling form factor of the vertices Σ b NB and Σ c ND. Considering these parameters, Eq.(5 can be written as: Π HAD (p,p,q = i 2m2 B[D] f B[D] λ N λ Σb [Σ c]g Σb NB[Σ cnd] m b[c] +m u (p 2 m 2 Σ b [Σ c] (p 2 m 2 N (q2 m 2 B[D] } (m N m Σb [Σ c] m 2 Σb [Σc] γ 5 +(m Σb [Σ c] m N p/γ 5 +q/p/γ 5 m Σb [Σ c]q/γ The OPE side + ( Now, we briefly outline the operator product expansion(ope for the three-point correlation Eq.(. Firstly, we contract the quark fields in the correlation with Wich s theorem. Π(p,p,q OPE = i 2 d 4 x d 4 ye ip.x e ip.y ǫ abc ǫ ijk γ 5 γ ν S cj d (y xγ µcsu bit (y xcγ νsu ah (yγ 5Sb[c] hk ( xγ µγ 5 γ 5 γ ν S cj d (y xγ µcs ait u (y xcγ νs bh u (yγ 5S hk b[c] ( xγ µγ 5 } Secondly, we replace the heavy and light quark propagators with the following full propagators[45 47], Sb[c] mn (x = i (2π 4 d 4 ke ik.x δ mn g s αβ mn σ αβ (k/+m b[c] +(k/+m b[c] σ αβ k/ m b[c] 4 (k 2 m 2 Q 2 + π2 αs k δ 2 +m b[c] k/ mn m b[c] + } 3 π (k 2 m 2 (2 b[c] 4 (

4 4 Su[d] mn (x = i x/ 2π 2 x 4δ mn m u[d] 4π 2 x 2δ mn qq 2 ig sλ ij A A [ ] θη 32π 2 x 2 x/σ θη +σ θη x/ + ( i m u[d] 4 ( x/ x2 92 m2 qq i m u[d] 6 x/ (3 Where m,n are the color indices, and qq is the uu and dd in Eq(3. After these above substitutions in Eq(, we carry out Fourier transformation in D=4 dimensions using the following formulas: [(y x 2 ] n = d D t (2π De it.(y x i( n+ 2 D 2n π D/2Γ(D/2 n ( D/2 n Γ(n t 2 (4 d D [y 2 ] n = t.y i( n+ 2 D 2n π D/2Γ(D/2 n ( D/2 n (2π De it Γ(n t 2 (5 Before the preformation of four-x and four-y integrals, the replacements x µ i p µ and y µ i p µ are carried out. After these processes, the integrals turn into Dirac delta functions which are used to simplifythefour-integralsoverk andt. ThefollowingstepistoperformtheFeynmanparametrization, after which the following function is used to carried out the remaining four-integral over t. d D t (2π D [t M 2 ] α = i( α Γ(α D/2 (4π D/2 Γ(α (M 2 α D/2 Where M 2 = m 2 b[c] x+p2 x(x+y +p 2 y(x+y q 2 xy. After further simplification, the three-point correlation in OPE side show the following Dirac structures: Π OPE (p,p,q = Π (q 2 γ 5 +Π 2 (q 2 p/γ 5 +Π 3 (q 2 q/p/γ 5 +Π 4 (q 2 q/γ 5 (6 Where eachπ i denotes contributionscoming fromperturbativeandnonperturbativeparts. In general, we expect that we can choose either dirac structure Π i (with i =,2,3,4 of the correlations Π(p,p,q to study the hadronic coupling constants. In our calculations, we observe that the structure p/γ 5 and q/γ 5 are the pertinent dirac structures. After taking its imaginary parts ofπ i, we get the spectral densities ρ i (s,s,q 2 ofthe corresponding Dirac structure. Using dispersion relations, each Π i can be written as: Π OPE i (Q 2 = ds ds ρpert i (s,s,q 2 +ρ non pert i (s,s,q 2 (s p 2 (s p 2 Where s = p 2, s = p 2 and Q 2 = q 2. As examples, we give the perturbative and nonperturbative (7

5 5 parts of the spectral densities for the two Dirac structures p/γ 5 and q/γ 5 p/γ 5 (s,s,q 2 [ ] = dx dy 32π 4 2m b[c] (x+y+m d (3x+3y m u (x+y (x+y [ ] [ x(m 2 b[c] +Q2 y+sx(x+y +s y(x+y +s(x+y m b[c] x(2x+2y ] ( +3m d (x (x+y m u (x 2 +x(y 4 2y +3 +s (x+y [m b[c] x(2y ( ] +2(y y +y 3m d (x+y m u (x+y 2 +9m b[c] m d m u x+9m b[c] m d m u y ρ pert x 6m b[c] m d m u 3m b[c] m 2 u x 3m b[c]m 2 u y +2m b[c]q 2 x 2 y m b[c] Q 2 x 2 +2m b[c] Q 2 xy 2 m b[c] Q 2 xy 6m d m 2 u x 6m dm 2 u y +6m dm 2 u +3m dq 2 x 2 y +3m d Q 2 xy 2 3m d Q 2 xy } 3m d Q 2 y 2 +3m d Q 2 y m u Q 2 x 2 y m u Q 2 xy 2 +2m u Q 2 xy +2m u Q 2 y 2 3m u Q 2 y Θ[H 2 (s,s,q 2 ] x ρ pert q/γ 5 (s,s,q 2 [ ( = dx dy 32π 4 (x+y 2 (x+y s m b[c] xy(2x+2y +3m d (x 2 (y +x(y 2 3y + (y y m u (x 2 (y 3+x(y 2 7y+3+(3 2yy ( +s y m b[c] (x(2y +2(y y+3m d (y (x+y m u (x(y 3+y 2 5y +3 +9m b[c] m d m u y 6m b[c] m d m u 3m b[c] m 2 uy +2m b[c] Q 2 xy 2 m b[c] Q 2 xy 6m d m 2 uy +3m d Q 2 xy 2 ] 3m d Q 2 xy 3m d Q 2 y 2 +3m d Q 2 y m u Q 2 xy 2 +3m u Q 2 xy +2m u Q 2 y 2 3m u Q 2 y [ ] m b[c] (x(6y +6(y y+3m d (3y 2(x+y +m u ( 3x(y 2 3y 2 +y 6 [ ]} x(m 2 b[c] +Q2 y+sx(x+y +s y(x+y Θ[H 2 (s,s,q 2 ] x ρ non pert p/γ 5 (s,s,q 2 = dx dy 3( uu dd ( 2π 2 3x+3y Θ[H 2 (s,s,q 2 ] uu 48π 2 (m 2 b[c] +Q2 2 6m 3 b[c] m d 2m 3 b[c] m u 6m 2 b[c] m dm u +6m 2 b[c] m2 u +s [ 2m 2 b[c] m b[c]m u +2Q 2] 2m 2 b[c] Q2 +6m b[c] m d Q 2 +s [m b[c] (m u 2m b[c] 2Q 2] m b[c] m u Q 2 3m d m u Q 2 +3m 2 u Q2 2Q 4} Θ[H [s,s,q 2 ]] x 2 (3x+3y 2m b[c] dx dy α s π x3 6π 2 δ[h 2 [s,s,q 2 ]] (x+y qqg + 8 4π 2 9(m 2 b[c] +Q2 4 8m 6 b[c] 8m5 b[c] m d 36m 5 b[c] m u 8m 4 b[c] m dm u +8m 4 b[c] m2 u [ +39m 4 b[c] Q2 36m 3 b[c] m dq 2 32m 3 b[c] m uq 2 +3m 2 b[c] Q4 3s m 3 b[c] (3m b[c] 2m u +4m 2 b[c] Q2 +Q 4] +3s [ 3m 4 b[c] 2m3 b[c] m u +4m 2 b[c] Q2 +Q 4] 8m b[c] m d Q 4 2m b[c] m u Q 4 +9Q 6} Θ[H [s,s,q 2 ]]

6 x ρ non pert q/γ 5 (s,s,q 2 3 } = dx dy 6 2π 2 (48y 32 dd +(32 48y uu Θ[H 2 (s,s,q 2 ] uu ] 48π 2 (m 2 b[c] +Q2 2 3m b[c] [2m 2 b[c] (m d 2m u +m b[c] m u (m u 5m d +2m d m 2 u +Q 2[ 2m 2 b[c] 6m b[c]m d +m b[c] m u +2m d m u ]+s [m b[c] (2m b[c] m u +2Q 2] 2Q 4} x Θ[H [s,s,q 2 2 (8xy x+8y 2 8y+2m b[c] ]]+ dx dy α s π x3 96π 2 (x+y 2 δ[h 2 [s,s,q 2 ]] qqg 8 4π 2 9(m 2 b[c] +Q2 4 27m 6 b[c] 54m5 b[c] m d 54m 5 b[c] m u 72m 4 b[c] m dm u +8m 4 b[c] m2 u [ 6m 4 b[c] Q2 +36m 3 b[c] m dm 2 u +s m 3 b[c] ( 5m b[c] +36m d +m u 4m b[c] Q 2 (6m b[c] 9m d m u 9Q 4] +2s [ m 3 b[c] (3m b[c] 72m d 2m u 2m b[c] Q 2 ( 3m b[c] +9m d +m u +3Q 4] 72m 3 b[c] m dq 2 64m 3 b[c] m uq 2 72m 2 b[c] m dm u Q 2 +45m 2 b[c] Q4 8m b[c] m d Q 4 6m b[c] m u Q 4 8m d m u Q 4 +2Q 6} Θ[H [s,s,q 2 ]] Where Θ denotes the unit-step function, and H [s,s,q 2 ], H 2 [s,s,q 2 ] are defined as: H [s,s,q 2 ] = s (8 H 2 [s,s,q 2 ] = x(m 2 b[c] +Q2 y+sx(x+y +s y(x+y (9 2.3 The strong coupling constant We perform a double Borel transformation[48] to the physiological as well as the OPE sides. Then, we equate these two sides, invoking the quark-hadron duality from which the sum rule is obtained. As an example, the form factors for the structure p/γ 5 is: m 2 p/γ5 Σ b [Σc] Σ b NB[Σ cnd] (Q2 = e M 2 e m2 N (m b[c] +m u (Q 2 +m 2 B[D] M2 2 m 2 B[D] f B[D]λ Σb [Σ c]λ N (m Σb [Σ c] m N s ds (m b[c] +m u+m d 2 u (2m u+m d 2 ds e s M 2 e s M 2 [ρ pert p/γ 5 (s,s,q 2 +ρ non pert p/γ 5 (s,s,q 2 Where M and M 2 arethe Borelparameters, s and u are twocontinuum threshold parameterswhich areintroduced to eliminate the h.r. terms. These parametersfulfill the followingrelations:m 2 i <s <m 2 i and m 2 o <u <m 2 o, where m i and m o are the masses of the incoming and out-coming hadrons respectively and m is the mass of the first excited state of these hadrons. 3 The results and descussions Present section is devoted to the numerical analysis of the sum rules for the coupling constants. The decay constants parameters used in this work are taken as f B =(248±23 exp ±25 Vub MeV[49], f D =(25.8 ± 8.5 ± 2.5MeV[5], λ N =(. ±.5eV 6 [5], λ Σb =(.62 ±.8eV 3 [52] and λ Σc =(.45 ±.5eV 3 [52]. We take the masses of the hadronic from reference[53], ]} 6 (2

7 7 where m B =( ±.7MeV, m D =(864.84±.7MeV, m N =( ±.2MeV, m Σb =(58.3±.9MeV,m Σc =(2452.9±.4MeVand of quarkm b =(4.8±.3eV,m c =(.275±.25eV, m d =( MeV, m u =( MeV. The vacuum condensates are taken to be the.5 standard values uu = dd = (.8 ±. (.24 ±.ev 3 [54], sg s σs = m 2 ss [54], m 2 = (.8 ±.ev 2, gs 2 = (.22 ±.4eV 4 [55]. From Eq.(2, we also know that the value of the form factor Σb NB[Σ cnd] is the function of the input parameters, including the Borel parameters M 2 and M2, 2 continuum threshold s and u, the momentum Q 2. The workingregionsfor the M 2 and M2 2 aredetermined by requiringnot only that the contributions of the higher states and continuum be effectively suppressed, but also that the contributions of the higher-dimensional operators are small. In other words, we should find a good plateau which will ensure OPE convergence and the stability of our results[48]. The plateau is often called Borel window. Considering these factors, the Borel windows are chosen as 7(3eV 2 M 2 < 4(7eV 2 and 3(2eV 2 M2 2 < 7(6eV 2 for the strong vertex Σ b NB(Σ c ND(see Figures-4. From these figures, we can see that the values are rather stable with variations of the Borel parameters, it is reliable to extract the form factors. In addition, the continuum parameters, s = (m i + i 2 and u = (m o + o 2 are employed to include the pole and to suppress the h.r. contributions. The values for i and o can not be far from the experimental value of the distance between the pole and the first excited state[48]. In general, these two continuum thresholds s and u are determined by the relations s (m i +.5eV 2 and u (m o +.5eV 2. According to these considerations, we take s = 37.4(7.6eV 2 and u =.99(.99eV 2 for the strong vertex Σ b NB(Σ c ND..2.2 structure pγ 5 M2=5eV 2 structure pγ 5 M=eV 2 structure qγ 5 M2=5eV 2 structure qγ 5 M=eV 2 Σb NB (ev (Q 2 =ev Σb NB (ev (Q 2 =ev M 2 (ev M2 2 (ev 2 FI. : Σb NB as a function of M 2 at average values of the continuum thresholds. FI. 2: Σb NB as a function of M 2 2 at average values of the continuum thresholds. However, in order to obtain the coupling constants, it is necessary to extrapolate these results into physical regions(q 2 <, which is realized by fit the form factors into suitable analytical functions. It is indicated that we should get the same values for the coupling constants for the different dirac structure p/γ 5 or q/γ 5 when we take Q 2 = m 2 B[D]. This above procedure can help us minimizing the

8 Σc ND (ev (Q 2 =ev structure pγ 5 M2=4eV 2 structure qγ 5 M2=4eV 2 Σc ND (ev (Q 2 =ev structure pγ 5 M=5eV 2 structure qγ 5 M=5eV M 2 (ev M2 2 (ev 2 FI. 3: ΣcND as a function of M 2 at average values of the continuum thresholds. FI. 4: ΣcND as a function of M 2 2 at average values of the continuum thresholds..5 QCD sum rules for structure pγ QCD sum rules for structure qγ 5 Fit function for structure pγ 5 6 Fit function for structure qγ 5 Σb NB (ev QCD sum rules for structure pγ 5 QCD sum rules for structure qγ 5 Fit function for structure qγ 5 Fit function for structure pγ 5 Σc ND (ev Q 2 (ev Q 2 (ev 2 FI. 5: Σb NB as a function of Q 2 at average values of the continuum thresholds and Borel mass parameters. FI. 6: ΣcND as afunction of Q 2 at average values of the continuum thresholds and Borel mass parameters. uncertainties in the calculation of the coupling constant, which will be quite clear in the following section. From our analysis, we observe that the dependence of the form factors on Q 2 can be well described by the following fit function(see Figures5-6 : Σb NB[Σ cnd](q 2 2 = C exp Q C 2 +C 3 exp Q2 C 4 (2 Where the values of C,C 2,C 3 and C 4 for different dirac structures are presented in Table for these two strong vertexes Σ b NB and Σ c ND. The fit function is used to determine the value of the strong coupling constant at Q 2 = m 2 B[D] for different structures, and the results are also presented in Table. The errors existing in these results arise from the uncertainties of the input parameters together with the uncertainties coming from the determination of the working regions of the auxiliary parameters. It is indicated from Figure 5, Figure 6 and Table that different dirac structure can give compatible

9 9 TABLE I: Parameters appearing in the fit function of the coupling form factor for Σ b NB and Σ cnd. Structure C (ev C 2(eV 2 C 3(eV C 4 (ev 2 Σ b NB p/γ 5.4± ±.2.9± ±4..55±. q/γ 5.78±. 4.99±.2.74± ±4..3±. Σ cnd p/γ 5 4.2±.2 2.5± ±.2 q/γ ± ± ± ±.4 results for each strong coupling constants when we take Q 2 = m 2 B[D] in the fit function(eq.2. For example, the results of the coupling constant for vertex Σ b NB are.55 and.3 for p/γ 5 and q/γ 5 structure respectively. Thus, we can take the average of the coupling constants from two different diracstructure for eachvertex, which are.43±.and 3.76±.5for Σ b NB and Σ c ND respectively. 4 Conclusion In this article, wehavecalculatedthe formfactorsofthe vertexesσ b NB andσ c ND in the space-like regions by three-point sum rules. Then we fit the form factors into analytical functions, extrapolated them into the time-likeregions,and obtainedthe strongcouplingconstants Σb NB and ΣcND. These calculated results can be used to analyze the related experimental results at LHC as well as the heavy ion collision experiments like PANDA at FAIR. Acknowledgment This work has been supported by the Fundamental Research Funds for the Central Universities, rant Number 26M S33. [] B. Aubert et al, Phys. Rev. Lett. 97, 232 (26 [2] K. Nakamura et al, J. Phys. 37,752(2 [3] T. Lesiak, hep-ex/6242. [4] J. L. Rosner, J. Phys. 34, S27 (27. [5] M. Paulini, arxiv: [6] E. Klempt and J. M. Richard, Rev. Mod. Phys. 82, 95 (2 [7] M. Mattson et al, Phys. Rev. Lett. 89, 2(22 [8] A.Ocherashvili et al, Phys. Lett. B 628, 8 (25 [9] A.Faessler, Th. utsche, M.A.Ivanov, J.. Korner, V.E. Lyubovitskij, D.Nicmorus, andk.pumsa-ard, Phys. Rev. D 73, 943 (26 [] B. Patel, A. K. Rai, and P. C. Vinodkumar, J. Phys. 35, 65 (28; J. Phys. Conf. Ser., 22 (28 [] Xiang Liu, H.X.Chen, Y.R.Liu et al. Phys. Rev. D 77, 43(28

10 [2] Chun Mu, Xiao-Wang, Xiao-Lin Chen etal. Chin. Phys. C38, 3 (24 [3] J.Z. Zhang, M.Q.Huang,Phys.Rev.D78, 945(28; Phys. Lett. B 674,28(29; Chin. Phys. C 33, 385(29 [4] M.Karliner, H.J.Lipkin,Phys. Lett. B 66,539(28 [5] M.Karliner, J.L.Rosner,Phys. Rev. D 9, 947(24 [6] F.S.Navarra, M.Nielsen, K.Tsushima,Phys. Lett. B 66, 335(25 [7] A.Khodjamirian, Th.Mannel, N.Offen et al.phys. Rev. D83,943 (2 [8] T.M.Aliev, K.Azizi, M.Savci,J. Phys. 4, 653(23; 4, 653(24; J. Phys.: Conf. Ser. 556, 26 (24 [9] K.Azizi, H.Sundu, S. Sahin, PHYS. REV. D9,4(24 [2] K.Azizi, Y.Sarac, H.Sundu, Nucl. Phys. A 943, 59 (25 [2] Z..Wang,Phys. Rev. C85,4524 (22 ;Eur.Phys.J.C7,86(2;Eur.Phys.J.C73,2533(23 [22] Z..Wang, T.Huang,Phys. Rev. C84,482 (2 [23] A. Kumar, Adv. High Energy Phys. 24, (24 [24] A. Hayashigaki, Phys. Lett. B 487, 96 (2 [25] F. S. Navarra, M. Nielsen, Phys. Lett. B 443, 285 (998 [26] A. Khodjamirian, Ch. Klein, and Th. Mannel, et al. J. High Energy Phys. 9,6 (2 [27] M.E. Bracco, M. Chiapparini, F.S. Navarra, M. Nielsen, Phys. Lett. B 659, 559 (28 [28] M.E. Bracco, M. Nielsen, Phys. Rev. D 82, 342 (2 [29] Z.. Wang, Phys. Rev. D 89, 347 (24 [3] A. Khodjamirian, Th Mannel, N. Offen, Y.M. Wang, Phys. Rev. D 83, 943 (2 [3] A. Khodjamirian, Ch. Klein, Th. Mannel, Y.M. Wang, arxiv:8.297 [hep-ph] [32] T.M. Aliev, M. Savci, arxiv: [hep-ph] [33] T.M. Aliev, M. Savci, arxiv: [hep-ph] [34] T. Doi, Y. Kondo, M. Oka, Phys. Rep. 398, 253 (24 [35] R. Altmeyer, M. oeckeler, R. Horsley et al., Nucl. Phys. Proc. Suppl. 34, 373 (994 [36] Z.. Wang, Phys. Rev. D 74, 47 (26 [37] A. Cerqueira Jr, B.O. Rodrigues, M.E. Bracco, Nucl. Phys. A 874, 3 (22 [38] B.O. Rodriguesa, M.E. Braccob, M. Chiapparinia, Nucl. Phys. A 929, 43 (24 [39] E. Yazici et al., Eur. Phys. J. Plus. 28(, 3 (23 [4] R. Khosravi, M. Janbazi, Phys. Rev. D 87, 63 (23 [4] R. Khosravi, M. Janbazi, Phys. Rev. D 89, 6 (24 [42] L.J. Reinders, H. Rubinstein, S. Yazaki, Phys. Rep. 27, (985 [43] P. Pascual, R. Tarrach, Lect. Notes Phys. 94, (984 [44] Z.. Wang, Z.Y. Di, Eur. Phys. J. A 5, 43 (24 [45] L.J. Reinders, H. Rubinstein, S. Yazaki, Phys. Rep. 27, (985 [46] P. Pascual, R. Tarrach, Lect. Notes Phys. 94, (984 [47] Z.. Wang, Z.Y. Di, Eur. Phys. J. A 5, 43 (24 [48] P.Colangelo, A. Khodjamirian, At the frontier of particle physics. in Handbook of QCD, vol. 3. (World Scientific, Singapore, 2, p arxiv:hep-ph/75

11 [49] A. Khodjamirian, B and D Meson Decay Constant in QCD, in Proceeding of 3rd Belle Analysis School, 22, 2 (KEK, Tsukuba, Japan, 2 [5] B. I. Eisenstein et al. (CLEO Collab., Phys. Rev. D78, 523 (28 [5] K. Azizi, N. Er, Eur. Phys. J. C 74, 294 (24 [52] K. Azizi, M. Bayar, A. Ozpineci, Phys. Rev. D 79, 562 (29 [53] K. A. Olive et al. (Particle Data roup, Chin. Phys. C 38, 9 (24 [54] B. L. Ioffe, Prog. Part. Nucl. Phys. 56, 232 (26 [55] V. M. Belyaev, B. L. Ioffe, Sov. Phys. JETP 57, 76 (983; Phys. Lett. B 287, 76 (992