arxiv: v1 [hep-ph] 13 Jan 2015

Transcription

1 A QCum rules alulation of the J/ψD s strong oupling onstant B. Osório Rodrigues, 1, 2 M. E. Brao, 3 M. Chiapparini, 1 and A. Cerqueira Jr. 3 1 Instituto de Físia, Universidade do Estado do Rio de Janeiro, Rua São Franiso Xavier 524, 255-9, Rio de Janeiro, RJ, Brazil. 2 Centro Brasileiro de Pesquisas Físias - CBPF, arxiv: v1 [hep-ph] 13 Jan 215 Rua Xavier Sigaud 15, , Rio de Janeiro, RJ, Brazil. 3 Fauldade de Tenologia, Universidade do Estado do Rio de Janeiro, Rod. Presidente Dutra Km 298, Pólo Industrial, , Resende, RJ, Brazil. Abstrat In this work, we alulate the form fators and the oupling onstant of the strange-harmed verte J/ψD s in the framework of the QCum rules b studing their three-point orrelation funtions. All the possible off-shell ases are onsidered,, D s and J/ψ, resulting in three different form fators. These form fators are etrapolated to the pole of their respetive off-shell mesons, giving the same oupling onstant for the proess. onstant is g J/ψD s = GeV 1. Our final result for the J/ψD s oupling 1

2 I. INTRODUCTION In the past few deades, the Quantum Chromodnamis Sum Rules (QCDSR) ommunit has been ativel working in the alulation of strong oupling onstants of a great variet of verties. Partiularl, our group has alulated a large number of harmed verties as the D Dπ [1], D Dρ [2], J/ψD D [3] among man others that were revisited in the review of ref. [4]. Some reent works in QCDSR target verties with the heav beaut mesons (B, B, B s, B s,...) as the B sbk [5], B sb s η ( ) [6], B ( ) B Υ [7, 8],B ( ) B J/ψ [8] and/or strangeharmed mesons (, D s, 1, D s,...) as the J/ψ [9, 1], J/ψD sd s [1], J/ψD s1, φd s1 [11], D s η ( ) [6] and man others. Regarding the strange-harmed verties, whih are a subjet of major interest nowadas, due to the detetion of various new eoti states, it is notieable that some of them have masses that lie just above the ( ) ( ) mass threshold. B eoti, we mean that the usual quark-model desription as q q pairs does not hold for these new states. An eample of suh eoti state is the Y (414), first observed b the CDF ollaboration [12], whih has an observed dea into the pair J/ψφ. Some interpretations for its internal struture are the tetraquark state [13], D s moleular state [14], D s D s moleular state [15 17] and hbrid state [15]. Among them, we highlight the moleular D s D s hpothesis, so that the dea Y (414) J/ψφ an be understood with the help of an intermediate dea as Y (414) D s D s J/ψφ. In this dea, the verte J/ψD s is present, thus a more preise knowledge of its oupling onstant and form fators ma help in the understanding of the fundamental onstitution of the Y (414) meson. Also, it ma be useful in the studies of other eoti states as the X(435) [18 21], Y (4274) [22] and the Y (466) [23, 24], for eample. In this work, we investigate the J/ψD s verte within the QCDSR formalism, following the development of our previous work for the J/ψ verte [9]. This development onsists in alulating the form fators from the three point orrelation funtion in two was: one with quarks and gluons degrees of freedom, the OPE side, and the other with hadroni degrees of freedom, the phenomenologial side. After that, a double Borel transform is performed on eah side and using the quark-hadron dualit priniple, we equate both OPE and phenomenologial sides. If a window of stabilit is obtained, the form fator and oupling onstant an be etrated. In the following sets., we present the QCDSR formalism and our results for this alu- 2

3 lation. II. FORMALISM The starting point to perform the sum rule, is the alulation of the three point orrelation funtion of the verte. In this ase, the J/ψD s is a vetor-vetor-pseudosalar (V V P ) verte. Eah meson will be represented b an interpolating urrent, whih ontains the quantum numbers of the meson. We an alulate three orrelation funtions [4]: one with the vetor meson J/ψ off-shell, another with the vetor meson D s off-shell and a third one with the pseudo-salar meson off-shell: Γ (M) (p, p ) = d 4 d 4 e ip e iq Γ (M) (, ) (1) where q = p p is the transferred momentum, M is the off-shell meson (M =, J/ψ, D s) and Γ (M) (, ) is presented in eqs. (2)-(4) below, Γ (J/ψ) µν (, ) = T {j Ds 5 ()j J/ψ µ ()j D s ν ()}, (2) Γ (D s ) µν (, ) = T {j Ds Γ (Ds) 5 ()j D s ν ()j J/ψ µ ()}, (3) µν (, ) = T {j D s ν ()j Ds 5 ()jµ J/ψ ()}, (4) where T is the time ordered produt and jµ J/ψ, j D s ν the mesons J/ψ, D s and respetivel. and j Ds 5 are the interpolating urrents of Aording to the QCDSR, these orrelation funtions an be alulated in two different was: using hadron degrees of freedom, alled the phenomenologial side, or using quark degrees of freedom, alled the OPE side. In the following subsets., we alulate the orrelations funtions for the phenomenologial and OPE sides, and the sum rule. A. The phenomenologial side The effetive Lagrangian, whih represents the proess in the J/ψD s verte, is [25, 26]: ( ) L J/ψD s = g J/ψD s ε αβγδ α ψ β γ D sδ D+ s + γ D + sδ D s, (5) where ε 123 = +1 is the Levi-Civita totall antismmetri tensor. 3

4 From this Lagrangian, we an obtain the verties of the hadroni proess neessar to the alulation of the phenomenologial side of the QCDSR. For the J/ψ, D s and off-shell ases the are, respetivel: D s(p)j/ψ(q) (p ) = ig (J/ψ) J/ψD s (q 2 )ɛ β (q, λ)ɛ δ (p, λ)q α p γ ε αβγδ, (6) J/ψ(p)D s(q) (p ) = ig (D s ) J/ψD s Ds(q2 )ɛ β (p, λ)ɛ δ (q, λ)p α q γ ε αβγδ, (7) J/ψ(p) (q) D s(p ) = ig (Ds) J/ψD s (q 2 )ɛ β (p, λ)ɛ δ(p, λ)p α p γε αβγδ, (8) where g (M) J/ψD s Ds(q2 ) is the form fator of the J/ψD s verte with meson M off-shell (M = J/ψ,, D s). In the alulation of the phenomenologial side, we also make use of the following hadroni matri elements: j P 5 P (q) = P (q) j P 5 = f P m 2 P m q1 + m q2, (9) V (q) j V = f V m V ɛ (q), (1) j V V (q) = f V m V ɛ (q), (11) where P is the pseudo-salar meson ( ), V is a vetor meson (J/ψ or D s), q is the fourmomentum of the respetive meson, m 2 P,V is the squared mass of the meson, f P,V are the dea onstants, m q1 and m q2 are the onstituent quarks of meson P and ɛ is the polarization vetor, with = µ(ν) for J/ψ(D s). Finall, the orrelation funtions (2)-(4) of the phenomenologial side for eah off-shell ase read: Γ phen(j/ψ) µν = Γ phen(d s ) µν = Γ phen(ds) µν = Cg (J/ψ) J/ψD s Ds(q2 )p λ p σ ε µνλσ (p 2 + m 2 D s )(q2 + m 2 J/ψ )(p 2 + m 2 ) + h.r., (12) Cg (D s ) J/ψD s Ds(q2 )p λ p σ ε µνλσ (p 2 + m 2 J/ψ )(q2 + m 2 D s )(p 2 + m 2 ) + h.r., (13) Cg (Ds) J/ψD s Ds(q2 )p λ p σ ε µνλσ (p 2 + m 2 J/ψ )(q2 + m 2 )(p 2 + m 2 D s ) + h.r., (14) where h.r. stands for the ontributions of higher resonanes and ontinuum states of eah meson and C is defined as: C = f f D s f J/ψ m 2 m D s m J/ψ (m + m s ). (15) From eqs. (12)-(14), it is lear that this verte has onl one tensor struture to work within the formalism of the QCDSR. 4

5 B. The OPE side The OPE side is obtained from eqs. (2)-(4), with the interpolating urrents written in terms of the quark fields: jµ J/ψ = γ µ, j D s ν = γ ν s and j D s 5 = i γ 5 s. B onstrution, the OPE side is given b an epansion known as Wilson s Operator Produt Epansion, dominated b the perturbative term and followed b non-perturbative ontributions: OP E(M) Γµν = Γ pert(m) µν + Γ non-pert(m) µν, (16) where Γ pert(m) µν is the perturbative term and Γ non-pert(m) µν are the non-perturbative ontributions to the orrelation funtion. When alulating the form fator, this epansion ehibits a rapid onvergene and an be trunated after a few terms [4]. Considering the similarities between the J/ψD s and both J/ψ and J/ψD D verties, we epet, from the former, a similar behavior regarding the OPE series as seen in the two latter [3, 9]. Therefore, it should be adequate to onsider non-perturbative ontributions up to the mied quark-gluon ondensate: Γ non-pert µν = Γ qq µν + Γ mq qq µν + Γ g2 G 2 µν + Γ qgσgq µν + Γ mq qgσgq µν. (17) Calulating eq. (16) using the non-perturbative ontributions of eq. (17) orresponds to alulating the diagrams of fig. 1. As a onsequene of the use of the double Borel transform, onl the J/ψ off-shell ase has ontributions from all the non-perturbative terms of eq. (17). All the non-perturbative terms, eept the gluon ondensates (fig. 1d-1i), are suppressed for both and D s off-shell ases. This has been taken into aount in fig. 1, wherefore these ontributions were omitted from the beginning. The off-shell ase is also omitted in fig. 1, as it is just an interhange between the mesons D s and of the D s off-shell ase. Using dispersion relations, the perturbative term (fig. 1a) for a given meson M off-shell an be written in the following form: Γ pert(m) µν (p, p ) = 1 4π 2 ρ pert(m) µν (s, u, t) dsdu, (18) (s p 2 )(u p 2 ) where the spetral densit ρ pert(m) µν (s, u, t) an be obtained from the Cutkosk s rules. The quantities s = p 2, u = p 2 and t = q 2 are the Mandelstam variables. In general, eq. (18) is the main ontributing term of the OPE series in a QCDSR alulation. For a V V P verte, the spetral densit an be parametrized as: ρ pert(m) µν (s, u, t) = 3 λ F (M) (s, u, t)p λ p σ ε µνλσ, (19) 5

6 where λ = (u + s t) 2 4us and F (M) is an invariant amplitude. For the J/ψ and off shell, the invariant amplitude an be written as: F (J/ψ) = (m s m )(A + B) m s, (2) F (Ds) = F (D s ) = (m m s )B m, (21) where A = [ k p ] k os θ s p, B = k os θ, s p os θ = 2p k u + ( 1) (ɛ) (m 2 s m 2 ) 2 p k k = k2 ɛm 2 s (1 ɛ)m 2, k = s + ɛ(m2 s m 2 ) 2, s p = s + u t λ 2, p s = 2 s,, and ɛ = (1) for ( ) (J/ψ) off-shell. The first non-perturbative ontribution to the orrelation funtion (16) is the quark ondensate showed in fig. 1b: Γ ss (J/ψ) µν = ss p λ p σ ε µνλσ (p 2 m 2 )(p 2 m 2 ). (22) Figure 1 is the first order mass term for the strange quark ondensate ontribution to eq. (16), numeriall less important and usuall negligible in the three point QCDSR: ) Γ ms ss (J/ψ) µν = m m s ss ( p 2 + p 2 2m 2 2 (p 2 m 2 ) 2 ( ) p 2 m 2 2 p λ p σ ε µνλσ. (23) Full epressions for the gluon ondensates ( g 2 G 2 ) of fig. 1d-1i and mied quark-gluon ondensates ( sgσ Gs ) of fig. 1j-1o for the J/ψ off-shell ase an be found in Appendi A. C. The QCD Sum Rule Two more steps are neessar in order to alulate the QCDSR to obtain the form fators. First, we make the hange of variables p 2 P 2, p 2 P 2 and q 2 Q 2 followed b a double Borel transform [27, 28] to both sides of the QCDSR in eqs. (12)-(14) and (16), whih involves the transformation: P 2 M 2 and P 2 M 2, where M and M are the 6

7 Borel masses. After that, we equate the phenomenologial and OPE sides, invoking the quark-hadron dualit from whih the sum rule is obtained. The seond step is to eliminate the h.r. terms from the phenomenologial side in eqs. (12)- (14). This is ahieved b the introdution of the ontinuum threshold parameters s and u in the OPE side and b taking advantage of the quark-hadron dualit. These parameters fulfill the following relations: 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 inoming and outoming mesons respetivel and m is the mass of the first eited state of these mesons. After performing these steps, we obtain the QCDSR epressions for the form fators: where BB [ g (J/ψ) J/ψDs (Q 2 ) = 3 s u λ 1 4π 2 s inf u inf F (J/ψ) e s M 2 e u M 2 dsdu + BB [Γ non-pert ], (24) C (Q 2 +m 2 D /M 2 J/ψ )e m2 s e m2 Ds /M 2 g (D s ) J/ψD s Ds(Q2 ) = g (Ds) J/ψD s (Q 2 ) = 3 4π 2 s s inf u u inf 1 λ F (D s ) e s M 2 e u M 2 dsdu + BB [ ] Γ g2 G 2, (25) C (Q 2 +m 2 J/ψ /M 2 e m2 Ds /M 2 Ds )e m2 [ ] Γ g2 G 2 3 s u λ 1 4π 2 s inf u inf F (Ds) e s M 2 e u M 2 dsdu + BB, (26) C (Q 2 +m 2 J/ψ /M 2 e m2 D /M 2 Ds )e m2 s [ ] ] stands for the double Borel transform and BB [Γ non-pert ] = BB Γ g2 G 2 the D ( ) s off-shell ases, as alread mentioned. The onstant C is defined in eq. (15). As in previous works, the definition for the oupling onstant g J/ψD s is given b [29]: lim g (M) Q 2 m 2 J/ψDs Ds(Q2 ). (27) M In order to alulate this, it is neessar to etrapolate the results for the form fator to the region of Q 2 <. From eqs. (24)-(26) it is lear that we an obtain the oupling onstant from three different form fators, one for eah meson off-shell. However, the oupling onstant must be the same regardless the form fator used for the etrapolation. This ondition is used to minimize the unertainties in the alulation of the oupling onstant, whih will beome quite lear in the following setion. for 7

8 III. RESULTS AND DISCUSSION Equations (24)-(26) show the three different form fators of the J/ψD s verte. The numerial alulation of these form fators give results that must be fitted to an analtial funtion of Q 2. In order to minimize the unertainties regarding the fitting proedure, it is required that these three form fators lead to the same oupling onstant [29]. This ondition also helps to redue the errors of finding the Borel masses and ontinuum thresholds neessar to the alulation. Regarding the Borel masses, whih satisf the relation M 2 = m2 o m 2 i M 2, the an assume an value within the Borel window of stabilit, in whih the phenomenologial and OPE side should be ompatible. The Borel window an be redued b requiring that the pole ontribution must be bigger than the ontinuum ontribution b at least 5% and that the perturbative term ontributes with more than 5% to the total orrelation funtion. We have worked with the average values of the form fators within the Borel window. Hene, it is not neessar to use an speifi Borel mass value, minimizing the unertainties assoiated to the hoie of this quantit [9]. The proedure is briefl skethed as follows: the average value of the form fator is alulated in the Borel window for eah value of Q 2 used. The standard deviation is then used to automatize the analsis of the stabilit of the form fator regarding the Borel masses and ontinuum threshold parameters. This stabilit riterion then indiates the optimal values of the ontinuum threshold parameters and the Borel windows we have to work with. Thus, not onl is a good stabilit ensured in the Borel window but also in the whole Q 2 interval, without the need to hose a speifi value for the Borel mass. The ontinuum threshold parameters, s and u, are defined as s = (m i + i ) 2 and u = (m o + o ) 2, where the quantities i and o have been determined imposing the most stable Borel window, as mentioned before. In order to inlude the pole and to elude the ontributions from higher resonanes and ontinuum states, the values for J/ψ, D s and Ds annot be far from the eperimental value (when available) of the distane between the pole and the first eited state [3, 31]. Our analsis has found that the best values are Ds =.6 GeV and D s = J/ψ =.5 GeV, whih leads to remarkabl stable Borel windows for the three off-shell ases, as it an be seen in fig. 2 for J/ψ and off-shell. Figures for Ds off-shell are omitted as the are ver similar to the off-shell ase. 8

9 Also in fig. 2, we an hek that the perturbative term is the leading one in the OPE series, followed b the quark ondensate (for the J/ψ off-shell ase) and the gluon ondensates. The summed ontribution of the non-perturbative terms, m s ss, sgσgs and m s sgσgs, is small ( 2%), whih means we ould easil neglet these terms without hanging our results signifiantl. Likewise, the gluon ondensates have a ver small ontribution to the total orrelation funtion in the and D s off-shell ases and the ma also be negleted. Besides the values of the Borel masses and the ontinuum thresholds, we also need to know the values of dea onstants, ondensates and hadrons and quark masses. The hadron masses used are m Ds = GeV, m D s = GeV and m J/ψ = 3.97 GeV [3]. The other parameters are presented in table II. In table I, we present the Q 2 and Borel windows found for eah form fator g (M) J/ψD s (Q 2 ) (where M is the meson off-shell), together with its parametrization and the orresponding oupling onstant g J/ψD s. Its assoiated error σ is alulated using the method that will be eplained later in this work. TABLE I: Parametrization of the form fators and numerial results for the oupling onstant of this work. Quantit J/ψ off-shell D s off-shell off-shell Q 2 (GeV 2 ) [.5, 2.] [1., 3.] [1., 4.] M 2 (GeV 2 ) [7., 7.9] [5.3, 7.3] [4.9, 5.9] g (M) J/ψD s Ds(Q2 ) A B+Q 2 Ae Q2 /B Ae Q2 /B A GeV 2.3 GeV GeV 1 B (GeV 2 ) g (M) J/ψD ± σ s Ds (GeV 1 ) The Borel windows presented in table I (M 2 row) satisf the alread mentioned onditions regarding the pole and ontinuum ontributions, as shown in fig. 3. The form fator obtained for the J/ψ off-shell ase was well adjusted b a monopolar urve, while for the Ds and off-shell ases, the form fators were well adjusted b eponential urves (fig. 4). The oupling onstants of eah off-shell ase are shown in table I. The present different values among them, however, when the unertainties are taken into aount, the are all 9

10 ompatible, as suggested b fig. 4, where the error bars stand for the unertainties (σ) of the oupling onstants also presented in table I. The unertainties were alulated following the same proedures used in our previous works [2, 5, 9]. The are estimated b studing the behavior of the oupling onstants when eah parameter is varied individuall between its upper and lower limits. Eperimental parameters have their own errors (shown in table II), while for the QCDSR parameters, we varied the thresholds s and u in ±.1 GeV and the Q 2 window in ±2%. For the error due to the Borel mass (M 2 ), we used the standard deviation of the form fator within the Borel window. After that, we alulated the average of all these ontributions and their respetive standard deviations to obtain the oupling onstant and its unertaint for a given off-shell ase. In table II, we an see how the variane of eah parameter affets the final value of the oupling onstant. TABLE II: Perentage deviation of the oupling onstants related to eah parameter. Parameter g (J/ψ) J/ψD s (%) g (D s ) J/ψD s (%) g (Ds) J/ψD s (%) f Ds = ± 6.1(MeV)[3] f D s = 31 ± 13(MeV) a f J/ψ = 416 ± 6(MeV)[3] m = (GeV)[3] m s = (MeV)[3] M 2 (GeV 2 ) b s ±.1(GeV), u ±.1(GeV) Q 2 ± 2% (GeV 2 ) b ss = (29 ± 15) 3 (MeV 3 )[32] g 2 G 2 =.88 ±.16(GeV 4 )[33] sgσ Gs = (.8 ±.2) ss (GeV 5 )[34] Fitting parameters (A and B) a Value obtained using f D s /f Ds = 1.17 from [35, 36]. b The intervals for these quantities are those of table I. Aording to table II, the variation of most of the parameters (the dea onstants, ondensates and the Borel mass for eample) has a small impat on the value of the oupling 1

11 onstant. On the other hand, the unertaint in the quark harm mass is the one with the biggest propagation over to the final value of the oupling onstant. The fat that the sensitivit related to the Borel mass is small was alread epeted, as we have ver stable Borel windows (fig. 2). IV. CONCLUSION In this work we have alulated the g J/ψD s oupling onstant b three different QCD sum rules: one with the meson off-shell, another with the D s meson off-shell and a third one with the J/ψ meson off-shell. This proedure allowed us to redue the unertainties related to the method, leading to ompatible oupling onstants, as seen in fig. 4. Taking the mean value between the numbers presented in table I, we obtain the following final result for g J/ψD s : g J/ψD s = GeV 1. (28) This oupling onstant was obtained from sum rules that respet the riteria of the pole dominane over the ontinuum, the perturbative ontribution of the OPE being the dominant one and the form fator stabilit regarding the Borel mass in the whole Borel window, as shown in figs. 2 and 3. Regarding the form fators, we observe some similarities with our previous work for the J/ψ verte [9]. Among suh similarities, there is the form fator given b a monopolar parametrization when the heaviest meson (J/ψ) is off-shell, while an eponential one is the ase when one of the lightest mesons ( or D s) is off-shell. Furthermore, the OPE series presents, in both verties, the same hierarh for the ontribution of eah term, as well as omparable ontributions among terms of the same dimension. We an also ompare our result for the g J/ψD s results of previous QCDSR works, presented in table III. oupling onstant (eq. (28)) with the The omparison with ref. [11] is straightforward, its result differs from ours b 3%, both being ompatible within 2σ. Similar onlusions an be obtained from heav quark effetive theor (HQET) [39], from where the relation g J/ψD s = g J/ψDs /m Ds stands, whih allows us to ompare the value of eq. (28) with the result of our previous work [9]. Using m Ds = GeV, we have g J/ψDs /m Ds = 3.4 GeV 1, a different result ( 3%) 11

12 TABLE III: Values of oupling onstants obtained using different methods. Method and referene QCDSR [11] Coupling onstant g J/ψD s = 3.3 ±.62GeV 1 QCDSR [3] g J/ψD D = 4. ±.6GeV 1 QCDSR [9] g J/ψDs = VMD [37, 38] g J/ψDD = 7.44 from the g J/ψD s in this work, but in agreement with the also sizable differene ( 23%) found using the relation g J/ψD D = g J/ψDD /m D in [3]. A third omparison an be made invoking the SU(3) smmetr [39]: g J/ψD s = g J/ψD D, from whih we an ompare the result of eq. (28) with the g J/ψD D of ref. [3]. From this smmetr, we an see that these two oupling onstants are ompatible with eah other within 1σ, with a differene between them of 7%. This differene is omparable with the one found using the SU(3) relation g J/ψDsDs = g J/ψDD in ref. [9]. Finall, we an ompare our result with the oupling onstant g J/ψDD from the vetor meson dominane (VMD) [37, 38] using the SU(3) and HQET relations, whih leads to g J/ψD s = 3.84 GeV 1. This oupling is also ompatible with the one of eq. (28) within 2σ. Aknowledgements This work has been partiall supported b CNPq and CAPES. [1] F. S. Navarra, M. Nielsen, and M. E. Brao, Phs. Rev. D 65, 3752 (22). [2] B. O. Rodrigues et al., Nul Phs A 852, 127 (211). [3] R. D. Matheus et al., Int J Mod Phs E 14, 555 (25). [4] M. E. Brao et al., Prog Part Nul Phs 67, 119 (212). [5] A. Cerqueira Jr., B. O. Rodrigues, and M. Brao, Nul Phs A 874, 13 (212). [6] E. Yazii et al., Eur Phs J Plus 128, 113 (213). [7] A. Abbasi and V. Bashir, Eur Phs J A 49, 13 (213). [8] Z.-G. Wang, Phs. Rev. D 89, 3417 (214). 12

13 [9] B. O. Rodrigues, M. Brao, and M. Chiapparini, Nul Phs A 929, 143 (214). [1] R. Khosravi and M. Janbazi, Phs. Rev. D 87, 163 (213). [11] R. Khosravi and M. Janbazi, Phs. Rev. D 89, 161 (214). [12] T. Aaltonen et al. (CDF Collaboration), Phs. Rev. Lett. 12, 2422 (29). [13] F. Stanu, J Phs G Nul Part Phs 37, 7517 (21). [14] Qiao, Cong-Feng and Tang, Liang, EPL 17, 311 (214). [15] N. Mahajan, Phs. Lett. B 679, 228 (29). [16] R. M. Albuquerque, M. E. Brao, and M. Nielsen, Phs. Lett. B 678, 186 (29). [17] X. Liu and S.-L. Zhu, Phs. Rev. D 8, 1752 (29). [18] X. Liu, Z.-G. Luo, and Z.-F. Sun, Phs. Rev. Lett. 14, 1221 (21). [19] Z.-J. Zhao and D.-M. Pan, , URL [2] Z. Jian-Rong and H. Ming-Qiu, Commun Theor Phs 54, 175 (21). [21] Y.-L. Ma, Phs. Rev. D 82, 1513 (21). [22] T. Aaltonen et al. (CDF Collaboration), , URL [23] X. L. Wang et al. (Belle Collaboration), Phs. Rev. Lett. 99, 1422 (27). [24] D.-Y. Chen, X. Liu, and T. Matsuki, , URL [25] X. Liu, B. Zhang, and S.-L. Zhu, Phs. Lett. B 645, 185 (27). [26] R. S. Azevedo and M. Nielsen, Phs. Rev. C 69, 3521 (24). [27] A. Khodjamirian, in Continuous Advanes in QCD 22, edited b K. A. Olive (22), hap. 4, pp [28] P. Colangelo and A. Khodjamirian, in At the Frontier of Partile Phsis, edited b M. Shifman (21), hap. 2, pp [29] M. Brao et al., Phs. Lett. B 521, 1 (21). [3] K. Nakamura et al., J Phs G Nul Part Phs 37, 1422 (21). [31] A. M. Badalian and B. L. G. Bakker, Phs. Rev. D 84, 346 (211). [32] C. MNeile et al., Phs. Rev. D 87, 3453 (213). [33] S. Narison, Phs. Lett. B 76, 412 (212). [34] B. Ioffe, Prog Part Nul Phs 56, 232 (26). [35] A. M. Badalian, B. L. G. Bakker, and Y. A. Simonov, Phs. Rev. D 75, 1161 (27). [36] D. Beirevi et al., Phs. Rev. D 6, 7451 (1999). [37] Y. Oh, T. Song, and S. H. Lee, Phs. Rev. C 63, 3491 (21). 13

14 [38] Y. Oh, W. Liu, and C. M. Ko, Phs. Rev. C 75, 6493 (27). [39] X. Liu, X.-Q. Zeng, and X.-Q. Li, Phs. Rev. D 74, 743 (26). Appendi A: Full epressions for g 2 G 2 and qgσgq orrelators - J/ψ off-shell ase. OPE. Equation (A1) is the double Borel transform of the gluon ondensate ontribution to the BB [ Γ g2 G 2 ] = g2 G 2 96π 2 =d,e,f,g,h,i dα 1 e 1 M 2 α 1 M 2 t t M 2 +M 2 + [ ] BB Γ g2 G 2 () = ( α 1 +α M 2 ) 1 M [ 2 m 2 m2 s α2 1 m2 M 2 ] α 1 M 2 1 p λ p σ ε µνλσ =d,e,f,g,h,i F () (α 1 ), (A1) where the inde = d to i orresponds to the diagrams of fig. 1d-1i, with F (d) -F (i) given b eqs. (A2)-(A7) respetivel. F (d) = m (2α 2 1m 2 M 2 3α 1 M 2 + 2α 1 m m s 2α 1 m 2 + 3)/4 (A2) F (e) = m (α 3 1tM 4 2α 2 1M 4 2α 2 1tM 2 + 4α 2 1m m s M 2 + 6α 1 M 2 + α 1 t + 4α 1 m 2 s 4α 1 m m s 4)/(8(α 1 M 2 1)) F (f) = (α 1 m M 2 + m s m )/2 F (g) = 2m s (2α1m 2 m s M 2 3α 1 M 2 + 2α 1 m 2 s 2α 1 m m s + 3)/(α 1 M 2 1) 3 F (h) = (2α 1 m M 2 + 3m s 2m )/(α 1 M 2 1) F (i) = (α 1 m M 2 + m s m )/(α 1 M 2 1) (A3) (A4) (A5) (A6) (A7) The double Borel transforms of the mied quark-gluon ondensate ontributions are given b eqs. (A8) and (A9). B M 2B M 2 [ Γ sgσgs ] = sgσgs 12M 4 M 4 ( 3 [ Q 2 M 2 M 2 + m 2 M 4 + 2m 2 M 2 M 2 + m 2 M 4] M 2 M 4 +3M 4 M 2) p λ p σ ε µνλσ e m 2 M 2 e m2 M 2. (A8) B M 2B M 2 [ Γ m s sgσgs ] = m s sgσgs 24M 6 M 6 ( 2m 3 [ M 6 + 3M 4 M 2 + 3M 2 M 4 + M 6] + m M 2 M 2 [ 12M 4 + 6M 2 M 2 + M 4]) p λ p σ ε µνλσ e m 2 M 2 e m2 M 2. (A9) 14

15 J/ψ(D s) J/ψ J/ψ (s) Ds(J/ψ) s( ) D s ss D s ss (a) (b) () J/ψ(D s) J/ψ(D s) J/ψ(D s) (s) (s) (s) D s(j/ψ) s( ) D s(j/ψ) s( ) D s(j/ψ) s( ) (d) (e) (f) J/ψ(D s) J/ψ(D s) J/ψ(D s) (s) (s) (s) D s(j/ψ) s( ) D s(j/ψ) s( ) D s(j/ψ) s( ) (g) (h) (i) J/ψ J/ψ J/ψ D s ss D s ss D s ss (j) (k) (l) J/ψ J/ψ J/ψ D s ss D s ss D s ss (m) (n) (o) FIG. 1: Contributing OPE diagrams for J/ψ(D s) off-shell. The off-shell ase orresponds to a permutation of the and D s mesons of the D s off-shell ase. 15

16 (Q 2 =1.GeV 2 ) Total Perturbative ss +ms ss g 2 G 2 sgσgs +ms sgσgs (Q 2 =1.GeV 2 ) Total Perturbative g 2 G 2 g (J/ψ) J/ψD s 2 1 g ( ) J/ψD s M 2 (GeV 2 ) M 2 (GeV 2 ) (a) (b) FIG. 2: OPE ontributions for J/ψ off-shell (panel (a)) and off-shell (panel (b)). Perentual ontribution Pole Continuum Perentual ontribution Pole Continuum M 2 (GeV 2 ) M 2 (GeV 2 ) (a) (b) FIG. 3: Pole and ontinuum ontributions for the J/ψ off-shell (panel (a)) and for off-shell (panel (b)), both at Q 2 = 1 GeV 2. 16

17 6 5 J/ψ off-shell D s off-shell off-shell (Q 2 ) g (M) J/ψD s Q 2 (GeV 2 ) FIG. 4: Form fators of the J/ψD s verte. 17