arxiv: v2 [hep-ph] 16 Sep 2011

Transcription

1 Strong coupling constants of heavy spin 3/ aryons with light pseudoscalar mesons arxiv:.54v [hep-ph] Sep T. M. Aliev, K. Azizi, M. Savcı Physics Department, Middle East Technical University, 53 Ankara, Turkey Department of Physics, Doğuş University, Acıadem-Kadıköy, 347 Istanul, Turkey Astract The strong coupling constants among memers of the heavy spin 3/ aryons containing single heavy quark with light pseudoscalar mesons are calculated in the framework of the light cone CD sum rules. Using symmetry arguments, some structure independent relations among different correlation functions are otained. It is shown that all possile transitions can e descried in terms of one universal invariant function whose explicit expression is Lorenz structure dependent. PACS numer(s):.55.hx, 3.75.Gx, 3.75.Jz taliev@metu.edu.tr permanent address:institute of Physics,Baku,Azeraijan kazizi@dogus.edu.tr savci@metu.edu.tr

2 Introduction The heavy aryons have een at the focus of much attention oth theoretically and experimentally during the last decade. The heavy quarks inside these aryons provide windows helping us see somewhat further under the skin of the nonperturative CD as compared to thelight aryons. Sinceapartofpolarizationofheavyquarkistransferredtothearyon, so investigation of polarization effects of these aryons can give information aout the heavy quark spin. Moreover, these aryons provide a possiility to study the predictions of heavy quark effective theory (HET). Besides the spectroscopy of heavy aryons, which have een discussed widely in the literature, their electromagnetic, weak and strong decays are very promising tools to get knowledge on their internal structure. In this decade, essential experimental results have een otained in the spectroscopy of heavy aryons. The + [ ] antitriplet states Λ + c, Ξ + c, Ξ c [ Λ + c (593), Ξ + c (79), Ξ c(79)] as well as the + [ 3 + ] sextet states Ω c,σ c,ξ c [Ω c,σ c,ξ c ] have een oserved []. Among the S wave ottom aryons, only the Λ, Σ, Σ, Ξ and Ω have een discovered. It is expected that the LHC will open new horizons in the discovery of the excited ottom aryon sates [] and provide possiility to study the electromagnetic properties of heavy aryons as well as their weak and strong transitions. The experimental progress in this area stimulates intensive theoretical studies (for a review see for instance [3, 4] and references therein). Theoretical calculations of parameters characterizing the decay of the heavy aryons will help us etter understand the experimental results. The strong coupling constants are the main ingredients for strong decays of heavy aryons. These couplings occur in a low energy scale far from the asymptotic freedom region, where the strong coupling constant etween quarks and gluons is large and perturation theory is invalid. Therefore, to calculate such coupling constants a nonperturative approach is needed. One of the most reliale and attractive nonperturative methods is CD sum rules [5]. This method is ased on CD Lagrangian and does not contain any model dependent parameter. The light cone CD sum rules (LCSR) method [] is an extended version of the traditional CD sum rules in which the operator product expansion (OPE) is carried over twists rather than the dimensions of the operators as in the case of traditional sum rules. In the present work, we calculate the strong coupling constants among sextet of the heavy spin 3/ aryons containing single heavy quark with the light pseudoscalar mesons in the framework of the LCSR. Using symmetry arguments, we show that all allowed strong transitions among memers of these aryons in the presence of light pseudoscalar mesons can e expressed in terms of only one universal invariant function. Note that the coupling constants of heavy spin / aryons with pseudoscalar and vector mesons have een recently calculated in [7, 8]. The heavy spin 3/ heavy spin / aryon pseudoscalar meson and heavy spin 3/ heavy spin / aryon vector meson coupling constants have also een calculated in the same framework in [9,]. It should e mentioned here that the couplings of the heavy aryons with mesons is first calculated in [] within the framework of HET. Rest of the article is organized as follows. In section, we derive some structure independent relations among the corresponding correlation functions, and demonstrate how the considered coupling constants can e calculated in terms of only one universal function. In this section, we also derive the LCSR for the heavy spin 3/ aryon light pseudoscalar

3 meson coupling constants using the distriution amplitudes (DA s) of the pseudoscalar mesons. Section 3 is devoted to the numerical analysis of the related coupling constants and discussion. Light cone CD sum rules for the coupling constants of pseudoscalar mesons with heavy spin 3/ aryons In this section, the LCSR for the coupling constants among heavy spin 3/ aryons and light pseudoscalar mesons are derived. We start our discussion y considering the following correlation function: Π µν = i d 4 xe ipx P(q) T {η µ (x) η ν ()}, () where P(q) is the pseudoscalar meson with momentum q, η µ is the interpolating current for the heavy spin 3/ aryons and T denotes the time ordering operator. The correlation function in Eq. () is calculated in two different ways: in terms of hadronic parameters called the physical or phenomenological representation, in terms of CD degrees of freedom y the help of OPE called theoretical or CD representation. Matching then these two representations of the same correlation function, we otain the CD sum rules for strong coupling constants. To suppress contriutions of the higher states and continuum, we apply the Borel transformation with respect to the momentum squared of the initial and final states to oth sides of the sum rules. We start our calculations y considering the physical side. Inserting complete sets of hadrons with the same quantum numers as the interpolating currents and isolating the ground states, we otain Π µν = η µ() B (p) B (p)p(q) B (p+q) B (p+q) η ν () (p m )[(p+q) m ] +, () where B (p+q) and B (p) are the initial and final spin 3/ states, and m and m are their masses, respectively. The dots in Eq. () represent contriutions of the higher states and continuum. It follows from Eq. () that, in order to calculate the phenomenological part of the correlation function, the following matrix elements are needed: η µ () B (p) = λ B u µ (p), B (p+q) η ν () = λ B ū ν (p+q), B (p)p(q) B (p+q) = g B B Pū α (p)γ 5 u α (p+q), (3) where λ B and λ B are the residues of the initial and final spin 3/ heavy aryons, g B B P is the strong coupling constant of pseudoscalar mesons with heavy spin 3/ aryons and

4 u µ is the the Rarita-Schwinger spinor. Using the aove matrix elements and performing summation over the spins of the Rarita Schwinger spinors defined as ( u µ (p,s)ū ν (p,s) = (/p+m) g µν + 3 γ µγ ν + p µp ν 3m p ) µγ ν p ν γ µ, (4) 3m s in principle, one can find the final expression of the correlation function in phenomenological side. However, the following two principal prolems are unavoidale: ) all Lorentz structures are not independent; ) not only the heavy spin 3/, ut also the heavy spin / states contriute to the physical side, i.e., the matrix element of the current η µ, sandwiched etween the vacuum and the heavy spin / states, is nonzero and determined in the following way: η µ B(p,s = /) = A(4p µ mγ µ )u(p,s = /), (5) wheretheconditionγ µ η µ = isimposed. Toremovethecontriutionoftheunwantedheavy spin / aryons and otain only independent structures, we order the Dirac matrices in a specific way and eliminate the ones that receive contriutions from spin / states. Here, we choose the γ µ /p/qγ ν γ 5 ordering of the Dirac matrices and otain the final expression Π µν = λ B λ B g ( B B P (p m )[(p+q) m )] g µν /p/qγ 5 +other structures with γ µ at the eginning and ) γ ν γ 5 at the end, or terms that are proportional to (p+q) ν or p µ, () and to calculate the strong coupling constant g B B P, we choose the structure g µν /p/qγ 5, which is free of the unwanted heavy spin / states. Now, we proceed to calculate the correlation function from CD side. For this aim, we need to know the explicit expression for the interpolating current of the heavy spin 3/ aryons. In constructing the interpolating current for these aryons, we use the fact that the interpolating current for this case should e symmetric with respect to the light quarks. Using this condition, the interpolating current for the heavy aryons with J = 3/ containing single heavy quark is written as } η µ = Aǫ {(q ac Cγ a µ q) c +(qcγ a µ )q c +( a Cγ µ q)q c, (7) where A is the normalization factor, and a, and c are the color indices. In Tale, we present the values of A and light quark content of heavy spin-3/ aryons. Before otaining the explicit form of the correlation functions on the CD side, we first try to otain relations among the correlation functions of the different transitions using some symmetry arguments. Next, we show that all possile transitions can e descried in terms of only one universal invariant function. We follow the approach given in [7 ], where the coupling constants of heavy spin / aryons with light mesons as well as spin 3/ aryons spin / aryons light mesons vertices are calculated (see also [ ] for the couplings among light aryons and light mesons). Here, we should stress that the relations which are presented elow are independent of the choice of Lorentz structures. We start our 3

5 Σ +(++) (c) Σ (+) (c) Σ () (c) Ξ (+) (c) Ξ () (c) Ω () (c) q u u d u d s q u d d s s s A /3 /3 /3 /3 /3 /3 Tale : The light quark content q and q for the heavy aryons with spin 3/ discussion y considering the Σ Σ π transition. The invariant function descriing this strong transition can e written in the following general form: Π Σ Σ π = g πūuπ (u,d,)+g π dd Π (u,d,)+g π Π (u,d,), (8) where g πūu, g π dd and g π show coupling of the π meson to the ūu, dd and states, respectively. The interpolating current of π meson is written as J π = u,d g π qq qγ 5 q, (9) where g = g πūu π dd =, g π =. The invariant functions Π, Π and Π descrie the radiation of π meson from u, d and quarks of heavy Σ aryon, respectively, and they can formally e defined as: Π (u,d,) = ūu Σ Σ, Π (u,d,) = dd Σ Σ, Π (u,d,) = Σ Σ. () Fromtheinterpolating current of Σ aryon, we see that it is symmetric under the exchange u d, so Π (u,d,) = Π (d,u,) and we immediately otain Π Σ Σ π = ] [Π (u,d,) Π (d,u,), () where under the SU () f limit, Π Σ Σ π =. Now, we proceed to otain the invariant function responsile for other transitions containing the π meson. The invariant function for Σ + Σ + π transition can e otained from the Σ Σ π case y making the replacement d u, and using the fact that η Σ µ = η Σ + µ, from which we get, 4Π (u,u,) = ūu Σ + Σ +. () The factor 4 on the left hand side appears due to the fact that each Σ + contains two u quark, hence there are 4 possile ways for radiating π from the u quark. From Eq. (), we get Π Σ + Σ + π = Π (u,u,). (3) 4

6 Similar arguments lead to the following result for the Σ Σ π transition: Π Σ Σ π = Π (d,d,). (4) Consider the strong transition Ξ () Ξ () π. The invariant function for this decay can e otained from the Σ Σ π case using the fact that η Ξ µ = η Σ µ (d s) and η Ξ µ = η Σ µ (u s). As a result, we get Π Ξ Ξ π = Π (u,s,), Π Ξ Ξ π = Π (d,s,). (5) Now, we proceed to find relations among the invariant functions involving charged π ± mesons. We start y considering the matrix element dd Σ Σ, where d quarks from the Σ and Σ form the final dd state, and u and quarks are the spectators. The matrix element ūd Σ + Σ explains the case where d quark from Σ and u quark from Σ + form the ūd state and the remaining u and are eing again the spectators. From this oservations, one expects that these matrix elements e proportional to each other. Our calculations support this expectation. Hence, Π Σ Σ + π = ūd Σ + Σ = dd Σ Σ Making the exchange u d in Eq. (), we otain Π Σ Σ π + = du Σ Similarly, one can easily show that = Π (d,u,). () ūu Σ = Σ Σ = Π (u,d,). (7) Π Σ + Σ π+ = Π (d,u,), Π Σ Σ π = Π (u,d,), Π Ξ Ξ π + = Π (d,s,), Π Ξ Ξ π = Π (u,s,). (8) Calculation of the coupling constants of the memers of heavy spin 3/ aryons to other pseudoscalar mesons can e done in a similar way as we did for the π meson. Here, we shall say that in our calculations, we ignore the mixing etween η and η mesons and only consider η 8 instead of physical η meson. The interpolating current for η 8 meson has the following form J η8 = [ūγ 5 u+ dγ 5 d sγ 5 s], (9) we see that g η8 ūu = g η8 dd =, and g η8 ss =. () 5

7 For instance, consider the Σ Σ η 8 transition. Following the same lines of calculations as in the π meson case, we immediately otain, Π Σ Σ η 8 = [Π (u,d,)+π (d,u,)]. () The invariant function responsile for the Ξ Ξ η 8 transition can e written as: Π Ξ Ξ η 8 = g η8 ūuπ (u,s,)+g η8 ssπ (u,s,)+g η8 Π (u,s,) = [Π (u,s,) Π (u,s,)] = [Π (u,s,) Π (s,u,)]. () For the remaining transitions containing the η 8 meson we otain Π Σ + Σ + η 8 = Π (u,u,), Π Σ Σ η 8 = Π (d,d,), Π Ξ Ξ η 8 = [Π (d,s,) Π (s,d,)], Π Ω Ω η 8 = 4 Π (s,s,). (3) We also find the following relations for transitions involving K mesons: Π Ξ Σ K = Π Σ Ξ K = Π (d,u,), Π Ξ Σ K = Π Σ Ξ K + = Π (u,d,), Π Ξ Σ + K = Π Σ + Ξ K+ = Π (u,u,), Π Ω Ξ K = Π Ξ Ω K + = Π Ω Ξ K = Π Ξ Ω K = Π (s,s,), Π Ξ Σ K = Π Σ Ξ K = Π (d,d,), (4) The expressions for the charmed aryons can easily e otained y making the replacement c and adding to charge of each aryon a positive unit charge. So far we have otained all possile strong transitions among the heavy spin 3/ aryons withpseudoscalarmesonsthataredescriedintermsofonlyoneinvariantfunctionπ. This functioncanecalculatedindeepeuclideanregion,where p + and (p+q) + using the OPE in terms of DA s of the pseudoscalar mesons as well as light and heavy quark propagators. Inotainingtheexpression ofπ incdside, thenonlocalmatrixelements of types P(q) q(x)γq() and P(q) q(x)g µν q() appear, where Γ is any aritrary Dirac matrix. Up to twist 4 accuracy, these matrix elements are determined in terms of the DA s of the pseudoscalar mesons. These matrix elements as well as the explicit expressions of DA s are given in [7 9].

8 In calculation of the invariant function Π, we also need to know the expressions of the light and heavy quark propagators. The light quark propagator, in the presence of an external gluon field, is calculated in []: S q (x) = i/x π x m q 4 4π x qq [ /x ig s du i m q 3π G µνσ µν ( i m ) ( q 4 /x x 9 m qq i m ) q /x π x G µν(ux)σ µν i G 4π x uxµ µν (ux)γ ν ( ( )] x Λ ln )+γ E, (5) 4 where γ E.577 is the Euler constant, and Λ is the scale parameter. In further numerical calculations, we choose it as Λ = (.5 ) GeV (see [,]). The heavy quark propagator in an external gluon field is given as: S (x) = S free (x) ig s d 4 k (π) 4e ikx [ du /k +m (m k ) Gµν (ux)σ µν + ] u m kx µg µν γ ν,() where S free (x) is the free heavy quark operator in coordinate space and it is given y S free (x) = m K (m x ) 4π x i m /x 4π x K (m x ), (7) where K and K are the modified Bessel function of the second kind. Using the explicit expressions of the heavy and light quark propagators and definitions of DA s for the pseudoscalar mesons, we calculate the correlation function from the CD side. Equating the coefficients of the structure g µν /p/qγ 5 from othsides of the correlation function and applying the Borel transformation with respect to the variales p and (p+q) in order to suppress the contriutions of the higher states and continuum, we get the following sum rules for the strong coupling constants of the pseudoscalar mesons with heavy spin 3/ aryons: g B B P = e λ B λ B m M + m M Π, (8) where M and M are the Borel mass parameters correspond to the initial and final heavy aryons, respectively. The contriutions of the higher states and continuum are otained y invoking the duality condition, which means that aove the thresholds s and s the doule spectral density ρ h (s,s ) coincides with the spectral density derived fromcd side of the correlation function. The procedure for otaining doule spectral density from CD side and sutraction of higher states and continuum contriutions are explained in detail in [3], which we use in the present work. The masses of the initial and final aryons are equal to each other, so we take M = M = M and the residues λ B and λ B are calculated in [4]. The explicit expression for Π function is quite lengthy and we do not present its explicit form here. 7

9 3 Numerical results This section is devoted to the numerical analysis of the sum rules for strong coupling constants of pseudoscalar mesons with heavy spin 3/ aryons. The main input parameters of these LCSRs are DA s of the pseudoscalar mesons which are given in [7 9]. Some other input parameters entering to the sum rules are qq( GeV) = ( MeV) 3 and µ π (GeV) = f π m π /(m u +m d ) = (.43±.4)GeV [5], ss =.8 ūu, α π sg = (.±.4) GeV 4, m = (.8±.) GeV [], f π =.3 GeV, f K =. GeV and f η8 =.3 GeV [7] g Ξ Ξ π 3.. s = 45 GeV s = 44 GeV s = 43 GeV s = 4 GeV s = 4 GeV s = 4 GeV M (GeV ) Figure : The dependence of the strong coupling constant for the Ξ Ξ π transition at several fixed values of s. The sum rules for the strong coupling constants include also two auxiliary parameters: Borel mass parameter M and continuum threshold s. These are not physical quantities, hence the result for coupling constants should e independent of them. Therefore, we should look for working regions of these parameters, where coupling constants remain approximately unchanged. The upper limit of M is otained requiring that the contriution of the higher states and continuum is small and constitutes only few percent of the total dispersion integral. The lower ound of M is otained demanding that the series of the light cone expansion with increasing twist should e convergent. These conditions lead to the working region 5 GeV M 3 GeV for the ottom heavy spin 3/ aryons and 4 GeV M GeV for the charmed cases. In these intervals, the twist 4 contriutions does not exceed (4 )% of the total result. Our analysis also shows that contriution 8

10 Bottom Baryons Charmed Baryons g Ξ Ξ π 4±9 (4.±.5) g Ξ + c Ξ + c π ± (5.3±.) g Σ Σ π + 75±5 (8.±.7) g Σ + c Σ c π+ 4± (.±.7) g Ξ Σ + K 5±4 (7.8±.8) g Ξ + c Σ ++ c K 37±8 (.±.) g Ω Ξ K 73± (8.±.8) g Ω c Ξ + c K 37±9 (9.8±.8) g Σ + Σ + η 8 45± (5.8±.) g Σ ++ c Σ ++ c η 8 5± (7.8±.) g Ξ Ξ η 8 9±4 (.7±.3) g Ξ + c Ξ + c η 8.5±.8 (3.±.4) g Ω Ω η 8 95± (.5±.8) g Ω c Ω c η 8 5± (3.±.7) Tale : The values of the strong coupling constants for the transitions among the heavy spin 3/ aryons with pseudoscalar mesons. of the higher states and continuum is less than 5%. The continuum threshold s is not totally aritrary ut it is correlated to the energy of the first excited state with the same quantum numers as the interpolating current. Our calculations show that in the interval (m B +.4) GeV s (m B +.8) GeV, the strong coupling constants weakly depend on this parameter. As an example, let us consider the Ξ Ξ π transition. The dependence of the strong coupling constant for the Ξ Ξ π transition on M at different fixed values of the s is depicted in Fig. (). From this figure, we see that the strong coupling constant for Ξ Ξ π shows a good staility in the working region of M. This figure also depicts that the result of strong coupling constant has weak dependency on the continuum threshold in its working region. From this figure, we deduce g Ξ Ξ π = 4±7. From the same manner, we analyze all considered strong vertices and otain the numerical values as presented in Tale (). Note that, in this Tale, we show only those couplings which could not e otained y the SU() symmetry rotations. The errors in the values of the coupling constants presented in the Tale include uncertainties coming from the variations of the s and M as well as those coming from the other input parameters. Here, we should stress that in the Tale (), we only present the modules of the strong coupling constants, since the sum rules approach can not predict the sign of the residues of the heavy aryons. Our numerical calculations show that the HET are violated approximately 5% (%) for the coupling constants of the heavy aryons containing (c) quark. Finally, we also check the SU(3) f symmetry violating effects and see that they change the results maximally aout 8%. At the end of this section, it should e mentioned that the predictions of the sum rules on light aryon-meson couplings strongly depend on the choice of the structure (for more detail see [7]). In connection with this point, here immediately arises the question whether or not a similar situation occurs for the case of the heavy aryon-light meson couplings. In order to answer this question, we also analyze the coupling constants predicted y the g µν γ 5 structure, which are presented in Tale () (see the values in the rackets). From these results, it follows that the values of the strong coupling constants of the heavy hadrons containing (c) quark with light pseudoscalar mesons decrease y a factor of aout 8(4) 9

11 times for the g µν γ 5 structure. However, sticking on the same criteria as is used in [7], we find that the g µν /p/qγ 5 is a more pertinent Dirac structure. In conclusion, the strong coupling constants of light pseudoscalar mesons with heavy spin 3/ aryons have een studied within LCSR. Using symmetry arguments, the Lorenz structure independent relations among different correlation functions have een otained. It has een shown that all possile transitions can e written in terms of one universal invariant function. Furthermore, it has een oserved that the values of the coupling constants are strongly structure dependent similar to the case of light aryon-meson couplings. The numerical values of those strong coupling constants which could not e otained via the SU() symmetry rotations have een also presented.

12 Appendix : In this appendix, we present some details of our calculations, i.e. how we perform the Fourier and Borel transformations as well as continuum sutraction. For this aim, let us consider the following generic term: T = d 4 x e ipx du e iuqx f(u) K ν(m x ) ( x ) n, (.) where K ν is the modified Bessel function of order ν appearing in the propagator of heavy quark. Using the integral representation of the modified Bessel function K ν (m x ) = Γ(ν +/)ν πm ν ( x dt cos(m t) ) ν (t x ) ν+/, (.) we get T = d 4 x du e ipx Γ(ν +/)ν f(u) πm ν dt cos(m t) ( x ) n ν (t x ) ν+/, (.3) where P = p + uq. For further calculations we go to the Euclidean space. Using the identity we have T = i ν πm ν Γ( n ν ) duf(u) Z = n Γ(n) dt e im t dα α n e αz, (.4) dy y n ν dv v ν e vt d 4 xe i P x y x v x, where means vectors in Euclidean space and we will consider only the real part of the complex exponential function, e im t. After performing Gaussian integral over x, we otain T = i ν π πm ν Γ( n ν ) duf(u) dt e im t dy y n ν The next step is to perform the integration over t. As a result we otain T = iν π m ν Γ(n ν ) Let us define new variales, duf(u) dy y n ν λ = v+y, dv v ν e m P 4(y+v) dv v ν e vt e (y +v). (.) 4v e P 4(y+v) (y +v). (.7) τ = y v +y, (.8) (.5)

13 hence T = iν π m ν Γ(n ν ) duf(u) dλ dτ λ n+ν 3 τ n ν ( τ) ν e m 4λ( τ) e P 4λ (.9) Performing Doule Borel transformation with respect to the p and ( p+ p) using the we get B(M )B(M )T = i ν π m ν Γ(n ν ) δ( M B(M )e αp = δ(/m α), (.) u 4λ )δ( M duf(u) Performing integrals over u and λ, we otain B(M)B(M )T = iν 4 π m ν Γ(n ν) dλ dτ λ n+ν 3 τ n ν ( τ) ν e m 4λ( τ) e u(u ) q 4λ u ). (.) 4λ ( M dτf(u ) 4 )n+ν τ n ν ( τ) ν e m q M ( τ) e M +M, (.) where, u = M M +M and M = M M. Replacing τ = x, we will have M +M B(M )B(M )T = iν+ 4 π m ν Γ(n ν ) ( )n+ν M dxf(u ) x n ν ( x ) ν e m 4 Finally, after changing the variale η = and using q = m x P, we get where B(M )B(M )T = iν+ 4 π m ν Γ(n ν) f(u ) ( ) Ψ α,β, m = M Γ(α) ( M 4 q M ( x ) e M +M. (.3) )n+ν ( ) m P e M +M Ψ α,β, m, (.4) M dηe ηm M η β α (η ) α, (.5) with α = n ν and β = ν. Now, let us discuss how contriution of the continuum and higher states are sutracted. For this aim we consider a generic term of the form ( A = (M ) n f(u )Ψ α,β, m M ). (.) We should find the spectral density corresponding to this term (see also [8]). The first step is to expand f(u ) as f(u ) = Σa k u k. (.7)

14 As a result we get ( ) M A = M n ( M Σa M k +M M +M ) k Γ(α) dηe η m M η β α (η ) α. (.8) Introducing new variales, σ = M and σ =, we have M σ k A = Σa k (σ +σ ) n+k Γ(α) = Σa k σ k Γ(n+k)Γ(α) = Σa k σ k Γ(n+k)Γ(α) = Σa k ( ) k Γ(n+k)Γ(α) dηe ηm (σ +σ ) η β α (η ) α dηe ηm (σ +σ ) η β α (η ) α dξe ξ(σ +σ ) ξ n+k dηη β α (η ) α dξξ n+k e (ξ+ηm )(σ +σ ) dηη β α (η ) α dξξ (( n+k d dξ )k e (ξ+ηm )σ ) e (ξ+ηm )σ. Applying doule Borel transformation with respect to σ s and σ s, we otain the spectral density ( ) k ρ(s,s ) = Σa k Γ(n+k)Γ(α) dηη β α (η ) α (.9) dξξ n+k ( ( d dξ )k δ(s (ξ +ηm )) ) δ(s (ξ +ηm )). (.) Performing integration over ξ, finally we otain the following expression for the doule spectral density: ( ) k ( ρ(s,s ) = Σa k dηη β α (η ) α (s ηm Γ(n+k)Γ(α) ) n+k ( d ) ) k δ(s s ) ds θ(s ηm ), (.) or ( ) k ρ(s,s ) = Σa k Γ(n+k)Γ(α) s /m ( dηη β α (η ) α (s ηm ) n+k ( d ) ) k δ(s s ). ds Using this spectral density, the continuum sutracted correlation function in the Borel scheme corresponding to the considered term can e written as: Π su = s m s ds ds ρ(s,s )e s /M e s /M. (.3) m Defining new variales, s = sv and s = s( v), we get s Π su = ds dv ρ(s,s )(4s)e sv/m e s( v)/m. (.4) m 3 (.)

15 Using the expression for the spectral density, one can get Π su ( ) k = Σa k Γ(n+k)Γ(α) sv/m Integrating over v, finally we otain s m ds dv ( ( d ) k s k dv )k δ(v /) dη η β α (η ) α (sv ηm ) n+k e sv/m e s( v)/m. (.5) s Π su ( ) k ( ) k = Σa k ds Γ(n+k)Γ(α) m k s k [ ( d sv/m ] dv )k dη η β α (η ) α (sv ηm )n+k e sv/m e s( v)/m. v=/ (.) References [] K. Nakamura et al., J. Phys. G 37, 75 (). [] G. Kane, (ed.), A. Pierce, (ed.), Perspectives on LHC physics, (Michigan U.) pp. Hackensack, USA: World Scientific (8) 337 p. [3] N. Roerts, M. Pervin, Int. J. Mod. Phys. A 3, 87 (8). [4] N. Nielsen, F. S. Navarra and S. H. Lee, arxiv: hep ph/ [5] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 47, 385 (979). [] I. I. Balitsky, V. M. Braun, and A. V. Kolesnichenko, Nucl. Phys. B 3, 59 (989) [7] T. M. Aliev, K. Azizi, M. Savci, Phys. Lett. B 9 (). [8] T. M. Aliev, K. Azizi, M. Savci, Nucl. Phys. A 85 () 4. [9] T. M. Aliev, K. Azizi, M. Savci, Eur. Phys. J. C 7, 75 (). [] T. M. Aliev, K. Azizi, M. Savci, V. S. Zamiralov, Phys. Rev. D 83, 97 (). [] S. L. Zhu and Y. B. Dai, Phys. Lett. B 49, 7 (998); P. Z. Huang, H. X. Chen, S. L. Zhu, Phys. Rev. D 8, 947 (9). [] T. M. Aliev, A. Özpineci, S. B. Yakovlev, V. Zamiralov, Phys. Rev. D 74, (). [3] T. M. Aliev, A. Özpineci, M. Savcı and V. Zamiralov, Phys. Rev. D 8, (9). [4] T. M. Aliev, K. Azizi, A. Özpineci and M. Savcı, Phys. Rev. D 8, 93 (9). 4

16 [5] T. M. Aliev, A. Özpineci, M. Savcı and V. Zamiralov, Phys. Rev. D 8, 54 (). [] T. M. Aliev, K. Azizi, M. Savcı, Nucl. Phys. A 847 (). [7] P. Ball, JHEP, (999). [8] P. Ball, V. M. Braun, and A. Lenz, JHEP 5, 4 (). [9] P. Ball, R. Zwisky, Phys. Rev. D 7, 45 (5). [] I. I. Balitsky and V. M. Braun, Nucl. Phys. B 3, 54 (989). [] K. G. Chetyrkin, A. Khodjamirian, and A. A. Pivovarov, Phys. Lett. B 5, 5 (8). [] I. I. Balitsky, V. M. Braun, and A. V. Kolesnichenko, Nucl. Phys. B 3, 59 (989). [3] V. M. Belyaev, V. M. Braun, A. Khodjamirian and R. Rükl, Phys. Rev. D 5, 77 (995). [4] T. M. Aliev, K. Azizi, A. Özpineci, Phys. Rev. D 79, 55 (9). [5] A. Khodjamirian, Th. Mannel, N. Offen and Y. M. Wang, Phys. Rev. D 83, 943 (). [] V. M. Belyaev and B. L. Ioffe, Sov. Phys. JETP, 57, 7 (98); Phys. Lett. B 87 (99) 7; Phys. Lett. B 3 (993) 35. [7] T. Doi, Y. Kondo and M. Oka, Phys. Rept. 398, 53 (4). [8] V. A. Beilin, A. V. Radyushkin, Nucl. Phys. B (985). 5