arxiv: v2 [hep-ph] 25 Sep 2012

Transcription

1 Order-v 4 Relativistic Corrections to Υ Inclusive Decay into a Charm pair Wen-Long Sang, Hai-Ting Chen, 2 and Yu-Qi Chen 2 arxiv: v2 [hep-ph] 25 Sep 202 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 00049, People s Republic of China 2 Key Laboratory of Frontiers in Theoretical Physics, The Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 0090, People s Republic of China Abstract In this work, we determine the short-distance coefficients for Υ inclusive decay into a charm pair through relative order v 4 within the framework of NRQCD factorization formula. The shortdistancecoefficient oftheorder-v 4 color-singlet NRQCDmatrixelement isobtainedthroughmatching the decay rate of b b( 3 S [] c cgg in full QCD to that in NRQCD. The double and single IR divergences appearing in the decay rate are exactly canceled through the next-to-next-to-leadingorder renormalization of the operator O( 3 S [8] and the next-to-leading-order renormalization of the operators O( 3 P [8] ). To investigate the convergence of the relativistic expansion arising from the color-singlet contributions, we study the ratios of the order-v 2 and -v 4 color-singlet short-distance coefficients to the leading-order one. Our results indicate that though the order-v 4 color-singlet short-distance coefficient is quite large, the relativistic expansion for the color-singlet contributions in the process Υ c c+x is well convergent due to a small value of v. In addition, we extrapolate the value of the mass ratio of the charm quark to the bottom quark, and find the relativistic corrections rise quickly with increase of the mass ratio. PACS numbers: t, 2.38.Bx, 2.39.St, 3.20.Gd

2 I. INTRODUCTION NRQCD factorization formula[] provides a systematical approach to express the quarkonium inclusive decay rate (cross section) as the sum of product of short-distance coefficients and the NRQCD matrix elements. The short-distance coefficients can be expanded as perturbation series in coupling constant α s at the scale of the heavy quark mass m. The long-distance matrix elements can be expressed in a definite way with the typical relative velocity v of the heavy quark in the quarkonium. The relativistic corrections to the quarkonium decay and production have been widely studied. In some processes, the next-to-leading-order (NLO) relativistic corrections are sizable, some of which even surpass the leading-order (LO) contributions. For the quarkonium production, typical examples are the double charmonia /ψ+η c production at the B factories [2, 3] and the double η c production through η b decay [4]. For the quarkonium annihilate decay, it happens in the /ψ inclusive decay and /ψ(υ) inclusive decay into a lepton pair or charm pair [5 7]. It may arouse worries about the relativistic expansion in the NRQCD approach. In Refs. [8, 9], the authors resumed a collect of relativistic corrections to the process e + e /ψ+η c and found the relativistic expansion is of well convergence. Similarly, the authors of Ref. [4] studied the process η b η c (ms)η c (ns) and found the relativistic expansion is convergent very well, though the order-v 2 corrections are both negative and large. In Ref. [0], the authors considered the order-v 4 corrections to /ψ inclusive decay and found the contribution from the color-singlet matrix element in this order is not as large as the LO and NLO contributions. The order-v 6 corrections to this process were even calculated in Ref. [], where the convergence is furthermore confirmed. It has been shown that the order-v 2 corrections to the process Υ c c+x are extremely huge [6]. The ratio of the short-distance coefficient of the order-v 2 matrix element to that of the LO one approaches to 2. It seriously spoils the relativistic expansion. So it urges us to calculate the order-v 4 corrections and investigate the convergence of the relativistic expansion for this process. Moreover, since the momenta of some gluons can be simultaneously soft, there exist complicated IR divergences in the calculations. It is technically challenging to cancel the IR divergences and obtain the IR-independent short-distance co- 2

3 efficients through the color-octet mechanism, and it also provides an example to examine the NRQCD factorization formula. The remainder of this paper is organized as follows. In Sec. II, we describe the NRQCD factorization formula for the Υ inclusive decay into a charm pair. In Sec. III, we list our definitions and present the techniques used to compute the decay rates. We elaborate the calculations on determining the short-distance coefficients at relative order v 4 in Sec. IV. The QCD corrections to the color-octet operators are also calculated in this section. Sec. V is devoted to discussions and summaries. II. NRQCD FACTORIZATION FORMULA FOR Υ c c+x According to the NRQCD factorizationformula, throughrelative order v 4, thedifferential decay rate for Υ inclusive decay into a charm pair can be expressed as []: dγ[υ c c+x] = df ( 3 S [] O( 3 S [] m 2 Υ + df 2( 3 S [] P( 3 S [] m 4 Υ + df( S [8] 0 ) O( S [8] m 2 + df 3( 3 S [] Q m 6 ( 3 S [] Υ + df 4( 3 S [] Q m 6 2 ( 3 S [] + df(3 S [8] m 2 O( 3 S [8] Υ + =0,,2 Υ 0 ) Υ df( 3 P [8] ) m 4 O( 3 P [8] ) Υ, () where m indicates the mass of the bottom quark b, O( 2S+ L [c] ) indicates the NRQCD operator with the spectroscopic state 2S+ L [c] with the spin S, orbital angular momentum L, total angular momentum, and color quantum number c, and O( 2S+ L [c] ) H H O( 2S+ L [c] ) H indicates the NRQCD matrix element averaged over the spin states of H. Thecolorquantum number c = and8inthenrqcdoperatorsdenotethecolor-singlet and the color-octet respectively. The NRQCD operators in the factorization formula () are When we were doing this work, a work about order-v 4 corrections to gluon fragmentation appeared at arxiv recently [2]. In that work, the authors applied the similar techniques to cancel IR divergences. 3

4 defined by 2 O( 3 S [] = ψ σχ χ σψ, (2a) P( 3 S [] = [ ψ σχ χ ( ) i D σψ +ψ ( ] ) 2 i 2 D σχ χ σψ, (2b) O( 3 S [8] = ψ σt a χ χ σt a ψ, (2c) O( S [8] 0 ) = ψ T a χ χ T a ψ, (2d) O( 3 P [8] 0 ) = d ψ ( ) 2 i D σ T a χ χ ( i ) 2 D σ T a ψ, (2e) O( 3 P [8] = ψ ( i 2 D [i, σ j]) T a χ χ ( i 2 D [i, σ j]) T a ψ, (2f) O( 3 P [8] 2 ) = ψ ( i 2 D (i, σ j)) T a χ χ ( i 2 D (i, σ j)) T a ψ, (2g) Q ( 3 S [] = ψ ( ) 2 i 2 D σχ χ ( ) 2 i 2 D σψ, (2h) Q 2 ( 3 S [] = [ ] ψ σχ χ ( ) i D σψ +ψ ( ) 2 i 4 D σχ χ σψ, (2i) where ψ and χ are Pauli spinor fields that annihilates the bottom quark and creates the bottomantiquark, respectively, σ i denotesthepaulimatrix, D indicatesthegauge-covariant derivative, and d = 4 2ǫ represents the space-time dimensions. In (2), A [i, B j] and A (i, B j) indicate the antisymmetric tensor and the symmetric traceless tensor respectively, which are defined via [2] A [i, B j] 2 (Ai B j A j B i ), A (i, B j) 2 (Ai B j +A j B i ) d δij A k B k. (3a) (3b) In (), we omit the term associated with the matrix element of the operator [ ψ σχ χ ( ) D ge +ge ( ] D σψ ψ ) D ge +ge D σχ χ σψ, (4) 2 which is shown to be dependent and be a linear combination of the matrix elements of the operators (2h) and (2i) accurate up to relative order v 4 [0]. In order to get the decay rate, both the NRQCD matrix elements and the short-distance coefficients in () should be determined. The NRQCD matrix elements have been studied through many nonperturbative approaches, such as lattice QCD [5], nonrelativistic quark model [6], and fitting the experimental data [7, 8]. 2 The dimensional regularization in quarkonium calculations including definitions of operators was first given in Refs.[3, 4] 4

5 On the other hand, based on the factorization, the short-distance coefficients can be perturbatively determined through matching the decay rates of the relevant processes at parton level in full QCD to these in NRQCD. At the leading order in α s, the coefficients df( S [8] 0 ), df( 3 S [8], and df( 3 P [8] ) can bedetermined through the processes b b( S [8] 0 ) c cg, b b( 3 S [8] c c, and b b( 3 P [8] ) c cg respectively. The color-singlet short-distance coefficients can be determined through the process b b( 3 S [] c cgg. Among these short-distance coefficients, F( 3 S [8] has been calculated up to the next-to-leading order in α s in Ref. [9]. F ( 3 S [] and F 2( 3 S [] are obtained in Refs. [20 22] and Refs. [6, 7], respectively. As mentioned in Ref. [0], since there exists the relation Q 2 ( 3 S [] Υ = Q ( 3 S [] Υ (+ O(v 2 )), to order v 4, we can only determine the sum of the short-distance coefficients: 3 df( 3 S [] df 3( 3 S [] +df 4( 3 S []. (5) As we will see, in order to get the IR-independent coefficient df( 3 S [], we are required to calculate df( 3 S [8] and df(3 P [8] ), and make use of the color-octet mechanism to cancel the IR divergences appearing in the decay rate of the process b b( 3 S [] c cgg. III. KINEMATIC DEFINITIONS AND PHASE-SPACE DECOMPOSITIONS A. Kinematic definitions We assign m c to the mass of charm quark, and assign p and p 2 to the momenta of the incoming bottom quark b and antiquark b. The momenta satisfy the on-shell relations: p 2 = p2 2 = m2. p and p 2 can be expressed as linear combinations of their total momentum P and half their relative momentum p: p = P/2+p, p 2 = P/2 p. (6) We assign l,l 2 to the momenta of the final charm pair. Therefore the momentum of the virtual gluon yields to Q = l +l 2. In addition, we take k to the momentum of the gluon in the process b b c cg, and take k,k 2 to the momenta of the two gluons in the process b b c cgg. 3 In Refs. [24 26], the authors may provide a potential approach to distinguish df 3 ( 3 S [] and df 4( 3 S [] in (5), nevertheless, it is enough for us to determine df( 3 S [] in this work. 5

6 To facilitate the evaluation, we introduce a set of dimensionless variables x = 2k P, x P 2 2 = 2k 2 P, z = Q2 P 2 P 2, r = 4m2 c M, (7) 2 where M denotes the mass of the bottomonium. At the leading order in v, there is M = 2m. All the Lorentz invariant kinematic quantities can be expressed in terms of these new variables. B. Phase-space decompositions In this subsection, we present the techniques for phase-space calculations of the relevant processes. We decompose each phase-space integral into two parts, which is proved to significantly simplify the calculations in the following section.. b b( 3 S [] c cgg The process b b( 3 S [] c cgg involves four-body phase-space integral, which can be expressed as dφ 4 = d d k d d k 2 d d l d d l 2 (2π) d δ d (P k (2π) d 2k 0 (2π) d 2k2 0 (2π) d 2l 0 (2π) d 2l2 0 k 2 l l 2 ).(8) As the treatment in Ref. [6], we decompose the space-space integral into dφ 4 = dφ 3 dφ 2, (9) 2π where dφ 3 and dφ 2 are defined via dφ 3 = dφ 2 =P 2 d d Q d d k (2π) d 2Q 0 (2π) d 2k 0 d d l d d l 2 d d k 2 (2π) d δ d (P k (2π) d 2k2 0 k 2 Q), (2π) d δ d (Q l (2π) d 2l 0 (2π) d 2l2 0 l 2 ). (0a) (0b) On the other hand, we can also separate the squared amplitude of the process into two parts: the charm part and the bottom part, i.e., pol,col M( 3 S [] 2 = L (ab)µν H (ab) µν (3 S [], () 6

7 where polarizations and colors of the initial and final states are summed, a,b denote the color indices, and the charm part L (ab)µν is given by L (ab)µν = δab 2 g 2 s Q 4Tr[( l +m c )γ µ ( l 2 m c )γ ν ], (2) and the bottom part H µν (ab) ( 3 S [] ) accounts for the remainder. According to the current conservation, we have dφ 2 L (ab)µν = δab 2 ) ( g µν + Qµ Q ν L, (3) where the Lorentz invariance L is explicitly calculated L = α s 3z (2+ r z ) r z. (4) As we will see, since L is a common factor in the whole calculations, we get (4) in 4- dimensions. As a consequence, the decay rate yields Γ( 3 S [] = 2 2π L δ ab dφ 3 2 where a symmetry factor 2 Q 2 ( (ab) H µν (3 S [] g µν + Qµ Q ν Q 2 ), (5) is included to account for the indistinguishability of the two gluons. The second term in the parenthesis of (5) does not contribute, due to the current conservation. It is accustomed to reduce dφ 3 into the integration over two dimensionless variables: { } dφ 3 = P2 fǫe 2 2ǫγE ǫ dx 28π 3 dx 2 (x +x 2 +z)[( x ( x 2 ) z]. (6) Γ(2 2ǫ) where f ǫ ( 4πµ2 M 2 e γ E ) ǫ. Here µ in f ǫ represents the dimensional-regularization scale, and γ E denotes the Euler constant. The boundaries of z,x, and x 2 are readily inferred x z x 2 x z x, 0 x z, r z. (7) It is convenient to make a further change of variables: x = ( z)x, x 2 = ( z)( x)[ ( z)xy]. (8) ( z)x This change of variables is particularly useful in extracting the relativistic corrections. Through this transformation, the boundaries of the new variables are simplified to 0 x, 0 y. (9) 7

8 Finally, the decay rate reduces to Γ( 3 S [] = f2 ǫe 2ǫγ [ E dxdy Γ(2 2ǫ) where ( z)x y( y)x 2 ( x) 2 H( 3 S [] P 2 ( z) 3 x( x) δ ab 2 (4π) 4 ( z)x 2 ( gµν ) H (ab) ] ǫ ( z) 4ǫ L H( 3 S [], (20) µν (3 S [] ). (2) 2. b b( 3 P [8] ) c cg Analogously, we separate the phase-space integral of this process into dφ 3 2 = dφ 2 dφ 2 2, (22) 2π where dφ 2 is given in (0b) and dφ 2 2 indicates d d Q d d k dφ 2 2 = (2π) d δ d (P k (2π) d 2Q 0 (2π) d 2k 0 Q). (23) dφ2 2 can be straightforward integrated out as dφ 2 2 = ( z) 2ǫ 8π Γ( ǫ) Γ(2 2ǫ)Γ(+ǫ) f ǫ, (24) Similarly, we separate the squared amplitude into the charm part and the bottom part Hµν( ab 3 P [8] where ), and obtain the decay rate Γ = Γ( ǫ)f ǫe ǫγ E Γ(2 2ǫ) ( z) 2ǫ L H( 3 P [8] ), (25) H( 3 P [8] ) δab ( z)( g µν )Hµν ab(3 P [8] ). (26) 2 6π 2 3. b b( 3 S [8] c c The phase-space integral of this process is quite simple. Therefore we directly present the decay rate where H( 3 S [8] is defined via Γ( 3 S [8] = δ( z)h( 3 S [8] L, (27) H( 3 S [8] δab ( g µν )Hµν ab(3 S [8]. (28) 2 P 2 8

9 IV. DETERMINING THE SHORT-DISTANCE COEFFICIENTS In this section, we determine the differential short-distance coefficients of Υ include decay into a charm pair at relative order v 4. df( 3 S [8], df( 3 P [8] ), and df(3 S [] can be determined through calculating the decay rates of the perturbative processes b b( 3 S [8] c c, b b( 3 P [8] ) c cg, and b b( 3 S [] c cgg, respectively. A. S-wave color-octet We first determine the differential short-distance coefficient df( 3 S [8] through matching the decay rate of b b( 3 S [8] c c in full QCD to that in NRQCD. We also carry out the computations for the renormalization of the operator O( 3 S [8], which will produce mixing with the operator O( 3 P [8] ) at the next-to-leading order in α s and with the operator O( 3 S [] at the next-to-next-to-leading order in α s. Moreover we consider the renormalization of the operators O( 3 P [8] ), which will induce mixing with the operator O(3 S []. The renormalized operators will be utilized to cancel the IR divergences through the color-octet mechanism when we determine the differential short-distance coefficients df( 3 P [8] ) and df(3 S [].. df( 3 S [8] The corresponding factorization formula for b b( 3 S [8] c c is expressed as dγ( 3 S [8] = df(3 S [8] O( 3 S [8] m 2 b b( 3 S [8] ). (29) Through calculations, we obtain the expression of H( 3 S [8] defined in (28) as H( 3 S [8] = πα s 2m2. (30) Inserting (30) into (27), we immediately deduce the differential decay rate in full QCD. By making use of O( 3 S [8] =, where the quark pair state is non-relativistically b b( 3 S [8] normalized, we readily get df( 3 S [8] = δ( z) πα s 2 L. (3) 9

10 (a) (b) (c) (d) FIG. : The Feynman graphs for the NLO QCD corrections to the operator O( 3 S [8]. We suppress the graphs which do not give rise to operator mixing S [8] 3 P [8] In this subsection, we consider the NLO QCD corrections to the operator O( 3 S [8]. There are four diagrams illustrated in Fig. 4. The vertex in the middle signifies the operator O( 3 S [8]. Since the diagrams involve UV divergences, the operator needs to be renormalized. In this work, we uniformly use MS scheme to carry out renormalization. We express the renormalized operator as O( 3 S [8] MS = O(3 S [8] +δ O+δ 2 O +O(αs 3 ), (32) where we truncate the expansion up to order α 2 s, which is enough to current work. δ O and δ 2 O arethecorrespondingnloandthenext-to-next-to-leading-order(nnlo)counterterms respectively. For convenience, we can expand an operator O in terms of α s : O = O (0) +O () +O (2) +. (33) where the superscript (n) represents order-α n s contribution. In the following calculations, we will also utilize the color decompositions [27] T a T b T a T b = N2 c T a T a, 4Nc 2 N c T a T b T b T a = N2 c + N2 c 2 T a T a. 4Nc 2 2N c (34a) 4 Here weconsideronly the diagramswhich arerelevantto our currentwork. Otherdiagramswill takeeffect when one considers the NLO QCD corrections to the short-distance coefficients. The similar calculations can also be found in Refs. [27 30] 0

11 With the NRQCD Feynman rules [3], the color-octet contribution from Fig. (a) reads I a = (ig s) 2 T a T a 4m 2 N c d d k (2π) d p 0 k 0 ( p k) 2 +iǫ 2m 4i(p p p kp k/k 2 ) k 2 +iǫ = g2 s T a T a p p d 2 d d k 2m 2 N c d (2π) d k 3 = α ( s T a T a fǫ f ǫ 3πm 2 N c ǫ UV ǫ IR i p 0 k 0 ( p k) 2 2m i +iǫ ) p p, (35) where f ǫ (4πe γ E ) ǫ. Since the coefficient is proportional to ǫ UV ǫ IR, the factor d coming from the loop integral in (35) is replaced with 4 in the MS scheme. Moreover, we add the factor f ǫ which is always associated with the MS scheme. It is not hard to find that Fig. (b-d) give the same contributions as Fig. (a) up to a color factor. Summing all the contributions and employing (32), we get O( 3 S [8] MS = O( 3 S [8] (0) 5α s 9πm 2 fǫ δ O = 5α s 9πm 2 fǫ ǫ UV ǫ IR O( 3 P [8] ). O( 3 P [8] )+O(α2 s), (36a) (36b) 3. 3 P [8] 3 S [] and In the similar way, we get O( 3 P [8] )() = 8α ( s p 4 N fǫ f ) ǫ O( 3 S [] 27πm 2 (d ) 3 +O(α ǫ UV ǫ s), 2 (37) IR O( 3 P [8] ) MS = O(3 P [8] )(0) 8α s 27πm 2 p 4 N (d ) 3 fǫ ǫ IR O( 3 S [] +O(α 2 s), (38) where N =, (d )(d 2) 2, and (d 2)(d+) 2 for = 0,,2. In (37) and (38), we explicitly keep the factor d stemming from the S-wave extraction, i.e., p i p j p2 d δij. The same treatment will be implemented in the following section 5. 5 We thank G. T. Bodwin for communications on this point.

12 (a) (b) (c) (d) (e) (f) (g) (h) (i) (j) (k) (l) (m) (n) (o) (p) (q) (r) FIG.2: TheFeynmangraphsfortheNNLOQCDcorrections totheoperatoro( 3 S [8]. Wesuppress the graphs which do not give rise to operator mixing with O( 3 S [] S [8] 3 S [] We proceed to deal with the NNLO QCD corrections to the operator O( 3 S [8]. The relevant Feynman graphs are illustrated in Fig. 2. There are totally eighteen Feynman graphs which do contribution. We are able to project out the color-singlet contribution by employing the color decompositions T a T b T c T a T b T c = others, (39a) T a T b T c T b T a T c = others, (39b) where we merely keep the color-singlet part. Through simple analysis, we find the color factors are 7 27 Fig. 2(b-d,h,j,l,n,o,r). 2 for the diagrams Fig. 2(a,e-g,i,k,m,p,q) and for the diagrams 27 2

13 With the Feynman rules, the amplitude of Fig. 2(a) reads I 2a = (ig s) 4 d d k d d l i i 6m 4 (2π) d (2π) d p 0 k 0 ( p k) 2 +iǫ p 2m 0 k 0 ( p k) 2 +iǫ 2m i i p 0 k 0 l 0 ( p k l) 2 +iǫp 2m 0 k 0 l 0 ( p k l) 2 +iǫ 2m 4i(p p p kp k/k 2 ) 4i[(p k) (p k) (p k) l(p k) l/l 2 ] k 2 +iǫ l 2 +iǫ = g4 s (d 2) 2 d d k d d l (p p ) 2 4m 4 (d ) 2 (2π) d (2π) d k 3 l ( k + l ) 2, (40) where the color factor is suppressed. Analogous to (40), we are able to get all others I 2b = g4 s (d 2) 2 d d k d d l (p p ) 2 4m 4 (d ) 2 (2π) d (2π) d k 2 l 2 ( k + l ) 2, I 2c = g4 s (d 2) 2 d d k d d l (p p ) 2 4m 4 (d ) 2 (2π) d (2π) d k 3 l 2 ( k + l ), I f = g4 s (d 2) 2 d d k d d l (p p ) 2 4m 4 (d ) 2 (2π) d (2π) d k 3 l, 3 I h = I j = I m = I a, I i = I k = I n = I b, I d = I e = I g = I o = I p = I q = I r = I c, I l = I f. (4) Including the color factors, and summing over all these contributions, we get I = 5 2gs 4 (d 2) 2 d d 54 m 4 (d ) k d d l (p p ) 2 2 (2π) d (2π) d k 3 l 3 ( 20αs 2 p 4 fǫ = f ) 2 ǫ, (42) 243π 2 (d ) m 4 ǫ UV ǫ IR where we use p 2 = p 2. In addition, we also need to calculate the NLO QCD corrections to the counterterm δ O. Employing (37), we readily obtain δ O () = 5α s fǫ 9πm 2 ǫ UV O( 3 P [8] 40α )() s 2 p 4 fǫ = 243π 2 (d ) m 4 Combining (42) and (43), we get ( O( 3 S [8] (2) = 20α2 s p 4 243π 2 (d )m 4 3 ǫ UV ( fǫ f ) ǫ O( 3 S [] ǫ UV ǫ ).(43) IR f ǫ 2 + f ) ǫ 2 O( 3 S [] ǫ 2 UV ǫ 2 +δ 2 O, (44) IR

14 where the UV divergence can be canceled through the counterterm. Finally, we present the renormalized operator and the corresponding counterterm O( 3 S [8] MS = O(3 S [8] (0) 5α s 9πm 2 fǫ ǫ IR O( 3 P [8] )+ 20α2 s p 4 f2 ǫ O( 3 S [] 243π 2 (d )m 4 ǫ 2,(45a) IR δ 2 O = 20α2 s p 4 f2 ǫ O( 3 S [] 243π 2 (d )m 4 ǫ 2. (45b) UV B. P-wave color-octet In this subsection, we determine the short-distance coefficients df( 3 P [8] ) through calculating the decay rates of the perturbative processes b b( 3 P [8] ) c cg. The factorization formulas for these processes are expressed as dγ( 3 P [8] ) = df(3 P [8] ) O( 3 P [8] m 4 ) + df(3 S [8] O( 3 S [8] b b( 3 P [8] ) m 2 b b( 3 P [8] ). (46) We can utilize (25) to calculate the decay rates in full QCD. First, we need to obtain H( 3 P [8] ) defined in (26). To understand the IR structure and show the IR cancelation, here we separate H( 3 P [8] ) into two parts: H(3 P [8] ) = H d( 3 P [8] )+H s( 3 P [8] ). H d( 3 P [8] ) include the terms proportional to /( z), which contribute the whole IR divergences to the decay rates from the region with the real gluon being soft. H s ( 3 P [8] ) take in charge of the remainder, which is absent of any singularity. Correspondingly, we use the subscripts d, s in both dγ( 3 P [8] ) and df(3 P [8] ) to denote the contributions from the two parts. The advantage of this classification will be recognized when we determine the short-distance coefficient df( 3 S [].. H d ( 3 P [8] ) We directly present the expressions of H d ( 3 P [8] ) which are the same for = 0,,2. H d ( 3 P [8] ) = It is useful to make the following expansion = ( r) 2ǫ δ( z)+ ( z) +2ǫ 2ǫ 5α 2 s ( ǫ)p2 3m 4 ( z)(3 2ǫ), (47) ( z 4 ) + ( ) ln( z) 2ǫ z + +O(ǫ 2 ), (48)

15 where the plus functions are defined via r (a) + f(z) r [ ] a f(z) f(). (49) Inserting (47) into (25) and employing (48), we are able to get the differential decay rates dγ d ( 3 P [8] ), which embrace IR divergences. Incorporating the factorization formulas (46) and the expressions (3) and (36a), we find the IR divergences appearing on the left hand side of (46) are exactly canceled by the renormalized S-wave color-octet matrix element on the right hand side. It renders the short-distance coefficients df( 3 P [8] ) free of any singularity: df d ( 3 P [8] ) = 5α2 s Lf {[ ǫ 8 ǫ IR ] π ) ǫ ( f ǫ )+ 53 ( f 2ln( r)+ ǫ δ( z) ( ) [ 2+ 0ǫ ] +4ǫ z + 3 2ln 2 ( r) 0 3 ln( r) ( ) } ln( z). (50) z In (50), we keep the terms linearly dependent on ǫ, which will induce finite contributions to the short-distance coefficient df( 3 S [] H s ( 3 P [8] ) We present the expressions of H s ( 3 P [8] ) as H s ( 3 P [8] 0 ) = 5α2 s (5 z)( ǫ)p2, 2m 4 (3 2ǫ) (5a) H s ( 3 P [8] = s[2+z 2(+z)ǫ]p 2 5α2, 6m 4 (3 2ǫ) (5b) H s ( 3 P [8] 2 ) = 5α2 s[4+z (+z)ǫ]p 2 6m 4 (5 6ǫ+4ǫ 2 ). (5c) The corresponding differential decay rates can be achieved by multiplying a constant factor: df s ( 3 P [8] 0 ) df s ( 3 P [8] df s ( 3 P [8] 2 ) = 5(5 z)α2 slf ǫ 36 = 5α2 slf ǫ 8 = α2 slf ǫ 8 [ 2ǫln( z)+ 5ǫ ], (52a) 3 [ z +2 2ǫ(z +2)ln( z)+ ] ǫ, (52b) 2(z +5) 3 [z +4 2ǫ(z +4)ln( z)+ 5 ] (3z +69)ǫ. (52c) 5

16 C. S-wave color-singlet In this section, we determine the short-distance coefficient df( 3 S [] through calculating the decay rate of the process b b( 3 S [] c cgg. The corresponding factorization formula is expressed as dγ( 3 S [] = df ( 3 S [] O( 3 S [] m 2 H + df 2( 3 S [] P( 3 S [] m 4 H + df 3( 3 S [] Q m 6 ( 3 S [] H + df 4( 3 S [] Q m 6 2 ( 3 S [] + df(3 S [8] m 2 O( 3 S [8] H + H df( 3 P [8] ) m 4 O( 3 P [8] ) H, (53) where the subscript H in the matrix elements represents b b( 3 S []. As mentioned in Sec. II, the two matrix elements in the second line of (53) are equal at relative order v 4. Therefore, we will determine the combined short-distance coefficient df( 3 S [] defined in (5) in this subsection. There are six diagrams for this process. The formula for the decay rate is given in (20). Firstly, we need to calculate H( 3 S [], which is defined in (2). To separate the relativistic corrections 6, we expand H( 3 S [] in powers of p2 : H( 3 S [] = H (0) ( 3 S [] +H (2) ( 3 S [] p2 m 2 +H(4) ( 3 S [] p4 m 4 +O(p6 ), (54) where the first two orders have been considered in Ref. [6]. Our remainder task is to calculate H (4) ( 3 S [] and the corresponding decay rate. At order v4, the decay rate involves IR divergences. Analogous to the treatment in the previous subsection, we separate H (4) ( 3 S [] into three parts: H (4) ( 3 S [] = H (4) d (3 S [] +H (4) s ( 3 S [] +H (4) r ( 3 S []. H (4) d (3 S [] proportional to /( z) contributes the whole IR divergences to the decay rate in the region where the two real gluons are simultaneous soft. In the following, we will demonstrate that the IR divergences can be thoroughly canceled by the renormalized S-wave color-octet matrix element (45a) with the short-distance coefficient (3), together with the renormalized P-wave matrix elements (38) with the short-distance coefficients df d ( 3 P [8] ) given in (50). H(4) s ( 3 S [] contributes the whole IR divergence to the decay rate in the region where only one of the 6 Since in the region z r, the decay rate will develop a logarithmic dependence on r, i.e., lnr, which is sensitive to the value of r, we will not expand the r appearing in L. For further explanations, we refer the reader to Ref. [6]. 6

17 real gluons is soft, as a result, it should be proportional to either /x or /( x). We will show that the IR divergence can be thoroughly canceled by the renormalized P-wave matrix elements (38) with short-distance coefficients df s ( 3 P [8] ) given in (52). H(4) r ( 3 S [] takes in charge of the remainder and therefore corresponds to a finite contribution to the decay rate. Similarly, we use the subscripts d, s, r in both dγ( 3 S [] and df(3 S [] to denote the contributions from the three parts.. H (4) d (3 S [] In this subsection, we calculate the decay rate and short-distance coefficient related to H (4) d (3 S []. First, we present the expression H (4) d (3 S [] = 40α 3 s 243m 2 πx( x)( z) [3(2y 2 2y +)+(8y 2 8y +)ǫ 8y( y)ǫ 2 ].(55) This term is proportional to /( z). In addition, we notice it also contains the factor /x( x). Therefore, we expect it will contribute a double IR pole to the decay rate at the endpoints of z and x. Employing the expansion = ( r) 4ǫ δ( z)+ ( z) +4ǫ 4ǫ ( ) 4ǫ z + ( ) ln( z) z + +O(ǫ 2 ), (56) and integrating out the variables x,y, we are able to obtain the differential decay rate dγ d ( 3 S [] {[ = 20α3 s Lf2 ǫ p4 243πm 6 ǫ + 4 4ln( r) +8ln 2 ( r) 6ln( r) 7π2 2 ǫ ] 3 ( ) [ ] ( ) } 4 ln( z) δ( z) +2lnz (57) z ǫ z + Inserting (57) into the factorization formula (53), we see the IR divergences on the left hand side are exactly canceled by the S-wave color-octet contribution, together with the renormalized P-wave color-octet matrix elements with the short-distance coefficients df d ( 3 P [8] ) on the right hand side. Straightforwardly, we get {[ df d ( 3 S [] = 0α3 s L 729π 2π2 3 ln 2µ2 ] δ( z) where µ indicates the factorization scale. M + 0 2ln( r) ln µ2 2 3 M 2 +4ln2 ( r) ln( r)+ 3 9 ( ) [4ln µ2 z M +2lnz + 20 ] ( ) } ln( z) +8, (58) 2 3 z

18 2. H (4) s ( 3 S [] The expression of H s (4) ( 3 S [] reads H s (4) ( 3 S [] = 20α3 s 243πm 2 x { 8y 2 +8y z 3 + 2(z 5)(y2 y)+z ǫ x ( y +yz) 4 x [ 2( z) 3 y 6 +2(z 5)(z ) 2 y 5 6(z 3 z 2 +3z 2)y 4 24(z 2 5z +6)y 3 +(z 4 +z 3 +9z 2 27z +96)y 2 +2(z 3 2z 2 3z 8)y ( +z 2 +3z +6+ǫ 6( z) 3 y 6 +6(z 5)(z ) 2 y 5 2(z ) (z 2 34z +8)y 4 8(z 2 4z +2)y 3 (z 4 z 3 9z 2 3z 92)y 2 )]} 2(3z 2 +7z +2)y z 2 +z +2. (59) This term is proportional to either /x or /( x), therefore, it contributes single poles to the decay rate at the endpoints of x, i.e., x 0 and x. After integrating out the variables x,y, we have dγ s ( 3 S [] { [ = 20α3 s Lf2 ǫ p4 3z z 5 43z 4 0z 3 +29z 2 248z πm 6 ǫ 3( z) 4 ] } +3(3z 5 4z 4 +4z 3 +2z 2 9z +6)lnz 4(3z +5)ln( z). (60) ItisnothardtofindthattheIRdivergence in(60)isexactly canceledbytherenormalized P-wave color-octet matrix elements with the short-distance coefficients df s ( 3 P [8] ), as we expected. We present the final result df s ( 3 S [] [ = 0α3 sl {(3z +5)ln µ2 287π M + 6z 5 6z 4 59z 3 +53z 2 3z ( z) 4 ] } +3(3z 5 4z 4 +4z 3 +2z 2 9z +6)lnz 2(3z +5)ln( z). (6) 3. H (4) r ( 3 S [] Finally, we deal with H r (4) ( 3 S [] ), which will produce a finite contribution to the decay rate. Since the expression of H (4) r ( 3 S [] is both tedious and cumbersome, here we merely 8

19 present the differential short-distance coefficient as df r ( 3 S [] { ( )[ α 3 = sl z 5tan 6(86z 5 279z z z π( z) 4 z ( ) z 5694z +609)tan + z( z)( z z 2 z ] 7296z 3 476z 4 ) 6(299z z z z 2 205z +700) } lnz ( z)(6896z z z z ). (62) D. Summarizing the differential short-distance coefficients We summarize the differential short-distance coefficients calculated in the above: df( 3 S [8] df( 3 P [8] ) df( 3 S [] = πα sl 2 = 5α2 s L 8 = df d( 3 S [] δ( z), (63a) {[ ln µ2 M + 5 ] ( ) 2 3 2ln( r) δ( z) 2 +A }, (63b) z + df s( 3 S [] where A are 5 z 4+z, 2+z, and for = 0,,2 respectively dfr( 3S [], (63c) + V. DISCUSSIONS AND SUMMARIES A. Discussions Applying the formulas (63) for the differential short-distance coefficients obtained in the last section, we now make some discussions.. df( 3 S [] / in the limit of z 0 We first discuss the short-distance coefficient df( 3 S [] / in the limit of z 0. It is not hard to derive L df( 3 S [] = 35α3 s z 0 287π lnz 70α3 s 287π ln µ2 M πα3 s α3 s 4374π. (64) In (64), we notice that the limitation bears the logarithmic divergence lnz. It is no surprise, owing to the fact we actually do not regularize the singularity when z approaches to 0. 9

20 0 t z t 2 z 0 z t i z z FIG. 3: Ratios of the differential short-distance coefficients. The blue dashed line represents the distribution t (z), while the red dot-dashed line represents the distribution t 2 (z). We specify the factorization scale µ = M in t 2 (z). Moreover, we find that the coefficient of ln z equals exactly to the corresponding coefficient of IR divergence in the decay rate of the process b b( 3 S [] 3g up to a factor 2π 2, as our expectation (The constant factor originates in from (5).). On the other hand, when either of the two real gluons becomes soft, there exists IR divergence which is regularized in dimensional regularization and canceled by the renormalized P-wave color-octet matrix elements. The term proportional to ln µ2 M 2 is related to the divergence. We are delight with that the coefficient of ln µ2 M 2 is exactly double that of lnz. 2. Color-singlet differential short-distance coefficients In (63), the two types of plus functions ( z ) + and ( ln( z) z ) + diverge as z. Since these singularities arise when the momenta of the real gluons in the final states go to 0, the distributions are actually unreliable in this region. Nevertheless, the singularities in the distributions are smeared when one integrates out z and so the integrated short-distance coefficients are well behaved. In order to investigate the dependence of the relativistic corrections on the virtuality of the intermediate gluon, it is intriguing and enlightening to study two ratios: t (z) df 2 ( 3 S [] /df ( 3 S [] and t 2(z) df( 3 S [] /df ( 3 S [], where df ( 3 S [] and df 2 ( 3 S [] are defined in () and have been obtained in (53) of Ref. [6]. To see it clearly, we plot Fig. 3 to illustrate the two ratios. In the figure, we observe that t (z) which reflects the NLO relativistic corrections is negative and its magnitude rises quickly 20

21 TABLE I: Ratios of the short-distance coefficients for Υ c cgg. The charm mass is selected to that of the D meson [6, 22]. The masses of the D meson and bottomonia are taken from Ref. [32]. r F 2 ( 3 S [] /F ( 3 S [] F(3 S [] /F ( 3 S [] Υ(S) ln 2 ( µ2 M 2 )+0.9ln( µ2 M 2 ) Υ(2S) ln 2 ( µ2 M 2 )+0.6ln( µ2 M 2 ) Υ(3S) ln 2 ( µ2 M 2 )+0.4ln( µ2 M 2 ) with increase of the virtuality of the immediate gluon. However t 2 (z) which reflects the order-v 4 relativistic corrections is positive in small values of z and turns to negative in large values. The magnitude of t 2 (z) is sizable in most values of z. 3. Integrated color-singlet short-distance coefficients Finally, we integrate out the variable z and investigate the integrated short-distance coefficients. In Table. I, we list the ratios of the order-v 2 and the order-v 4 color-singlet short-distance coefficients to the LO one for the processes of Υ(nS) inclusive decay into a charm pair. We learn from the table that the ratio F( 3 S [] /F ( 3 S [] is both positive and sizable, nevertheless, the relativistic expansion is of well convergence due to a small value of v (e.g., v 2 0..). The situation is quite similar to the case for the process /ψ ggg [0]. To further study the relation between the relativistic corrections and r, we extrapolate the value of r and investigate the ratios: g (r) F 2 ( 3 S [] (r)/f ( 3 S [] (r) and g 2 (r) F( 3 S [] (r)/f ( 3 S [] (r). We illustrate the two functions in Fig. 4. From the figure, we find the NNLO relativistic corrections become more and more important with r increase. B. Summaries In this work, we determine the short-distance coefficients within the framework of NRQCD factorization formula for Υ inclusive decay into a charm pair through relative order v 4. The order-v 4 color-singlet differential short-distance coefficient df( 3 S [] is obtained through matching the decay rate of b b( 3 S [] c cgg in full QCD to that in NRQCD. The double and single IR divergences appearing in the decay rate are exactly canceled through the NNLO renormalization of the operator O( 3 S [8] and the NLO renormalization of the 2

22 80 60 g r g 2 r g i r r r FIG. 4: Ratios of the short-distance coefficients as functions of r. The blue dashed line represents the distribution g (r), while the red dot-dashed line represents the distribution g 2 (r). We specify the factorization scale µ = M in g 2 (r). operators O( 3 P [8] ). To investigate the magnitude of the relativistic corrections and the convergence of the relativistic expansion, we show both the ratios of the differential shortdistance coefficients t i (z) and the ratios of the short-distance coefficients g i (r). Our results indicate that though F( 3 S [] is quite large, the relativistic expansion from the color-singlet contributions in the process Υ c c+x are well convergent due to a small value of v. In addition, we extrapolate g i (r) to a large range of r, and find the relativistic corrections rise quickly with increase of r. Acknowledgments We thank G. T. Bodwin for helpful discussions. The research of W. S. was supported by China Postdoctoral Science Foundation under Grant No. 202M The research of H. C. and Y. C. was supported by the NSFC with Contract No [] G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 5, 25 (995) [Erratum-ibid. D 55, 5853 (997)] [arxiv:hep-ph/ ]. [2] E. Braaten and. Lee, Phys. Rev. D 67, (2003) [Erratum-ibid. D 72, (2005)] [hep-ph/02085]. 22

23 [3] Z. -G. He, Y. Fan and K. -T. Chao, Phys. Rev. D 75, 0740 (2007) [hep-ph/ [HEP- PH]]. [4] W. -L. Sang, R. Rashidin, U-R. Kim and. Lee, Phys. Rev. D 84, (20) [arxiv: [hep-ph]]. [5] W. -Y. Keung and I.. Muzinich, Phys. Rev. D 27, 58 (983). [6] H. -T. Chen, Y. -Q. Chen and W. -L. Sang, Phys. Rev. D 85, (202) [arxiv: [hep-ph]]. [7] H. -T. Chen, W. -L. Sang and P. Wu, Commun. Theor. Phys. 57, 665 (202). [8] G. T. Bodwin, D. Kang, T. Kim,. Lee and C. Yu, AIP Conf. Proc. 892, 35 (2007) [hep-ph/06002]. [9] G. T. Bodwin,. Lee and C. Yu, Phys. Rev. D 77, (2008) [arxiv: [hep-ph]]. [0] G. T. Bodwin and A. Petrelli, Phys. Rev. D 66, 0940 (2002) [hep-ph/020520]. [] Hee Sok Chung, une-haak Ee, ungil Lee and Wen-Long Sang, work in progress. [2] G. T. Bodwin, U-R. Kim and. Lee, arxiv: [hep-ph]. [3] E. Braaten and Y. -Q. Chen, Phys. Rev. D 55, 2693 (997) [hep-ph/96040]. [4] E. Braaten and Y. -Q. Chen, Phys. Rev. D 55, 752 (997) [hep-ph/970242]. [5] G. T. Bodwin, S. Kim and D. K. Sinclair, Nucl. Phys. Proc. Suppl. 34, 434 (994). G. T. Bodwin, D. K. Sinclair and S. Kim, Phys. Rev. Lett. 77, 2376 (996) [arxiv:hep-lat/ ]. [6] E.. Eichten and C. Quigg, Phys. Rev. D 52, 726 (995) [hep-ph/ ]. G. T. Bodwin, D. Kang, T. Kim,. Lee and C. Yu, in Quark Confinement and the Hadron Spectrum VII: 7th Conference on Quark Confinement and the Hadron Spectrum-QCHS7, edited by. Emilio, F. T. Ribeiro, N. Brambilla, A. Vairo, K. Maung, and G. M. Prosperi, AIP Conf. Proc. 892 (AIP, New York, 2007), p. 35 [arxiv:hep-ph/06002]. Y. Q. Ma, K. Wang and K. T. Chao, Phys. Rev. Lett. 06, (20) [arxiv: [hep-ph]]. [7] H. -K. Guo, Y. -Q. Ma and K. -T. Chao, Phys. Rev. D 83, 4038 (20) [arxiv: [hep-ph]]. [8] G. T. Bodwin, H. S. Chung, D. Kang,. Lee and C. Yu, Phys. Rev. D 77, (2008) [arxiv: [hep-ph]]. [9] Y. -. Zhang and K. -T. Chao, Phys. Rev. D 78, (2008) [arxiv: [hep-ph]]. [20] H. Fritzsch and K. H. Streng, Phys. Lett. B 77 (978) 299. [2] A. Y. Parkhomenko and A. D. Smirnov, Mod. Phys. Lett. A 9, 5 (994) 23

24 [arxiv:hep-ph/ ]. [22] D. Kang, T. Kim,. Lee and C. Yu, Phys. Rev. D 76, 408 (2007) [arxiv: [hepph]]. [23] H. S. Chung, T. Kim and. Lee, Phys. Rev. D 78, 4027 (2008) [arxiv: [hep-ph]]. [24] N. Brambilla, E. Mereghetti and A. Vairo, Phys. Rev. D 79, (2009) [Erratum-ibid. D 83, (20)] [arxiv: [hep-ph]]. [25] N. Brambilla, E. Mereghetti and A. Vairo, HEP 0608, 039 (2006) [Erratum-ibid. 04, 058 (20)] [hep-ph/060490]. [26]. P. Ma and Q. Wang, Phys. Lett. B 537, 233 (2002) [hep-ph/ ]. [27] A. Petrelli, M. Cacciari, M. Greco, F. Maltoni and M. L. Mangano, Nucl. Phys. B 54, 245 (998) [hep-ph/ ]. [28] M. Beneke, F. Maltoni and I. Z. Rothstein, Phys. Rev. D 59, (999) [hep-ph/ ]. [29] G. T. Bodwin, E. Braaten, D. Kang and. Lee, Phys. Rev. D 76, (2007) [arxiv: [hep-ph]]. [30] Z.-G. He, Y. FanandK.-T.Chao, Phys.Rev.D8, (200) [arxiv: [hep-ph]]. [3] G. T. Bodwin and Y. -Q. Chen, Phys. Rev. D 60, (999) [hep-ph/ ]. [32] K. Nakamura et al. [Particle Data Group],. Phys. G 37, (200). 24