arxiv: v1 [hep-ph] 15 Mar 2016

Transcription

1 Transition Form Factors of χ b (1P) B c lν in QCD arxiv: v1 [hep-ph] 15 Mar 16 K. Azizi 1, H. Sundu, J. Y. Süngü, N. Yinelek 1 Department of Physics, Dogus University, Acibadem-Kadikoy, 347 Istanbul, Turkey Physics Department, Kocaeli University, 4138 Izmit, Turkey kazizi@dogus.edu.tr hayriye.sundu@kocaeli.edu.tr jyilmazkaya@kocaeli.edu.tr neseyinelek@gmail.com Abstract The form factors of the semileptonic χ b (1P) B c lν decay are calculated using theqcdsumruleapproach. Theresults obtained arethenusedto estimate thedecay widths of this transition in all lepton channels. The orders of decay rates indicate that this transition is accessible at LHC for all lepton channels. PACS number(s):11.55.hx, 14.4.Pq, 13..He

2 1 Introduction Quarkonia are flavorless bound states composed of combinations of quarks and their antiquarks. Since discovery of J/ψ meson in 1974, many new quarkonia states have been detected. The quarkonia systems consist of charmonium and bottomonium. Due to the large mass there are no toponium bound states and no light quark-antiquark states because of the mixture of the light quarks in experiments. The large mass difference between charm and bottom quarks prevents them from mixing. The heavy quark bound states may provide key tools for understanding the interactions between quarks, new hadronic production mechanisms and transitions, the magnitude of the CKM matrix elements and also analyzing the results of heavy-ion experiments. The J/Ψ suppression in ultra relativistic heavy-ion collisions was first suggested as a signal of the formation of a quark-gluon plasma (QGP). However, most recently, attentions have shifted to the bottomonium states due to that they are more massive than charmonium states. The bottom quarks and anti-quarks are relatively rare within the plasma, so the probability for reproduction of the bottomonium states through recombination is much smaller than for charm quarks. Consequently, it is expected that the bottomonium system to be a cleaner probe of the QGP than the charmonium system. Hence investigations on the properties of bottomonium systems can help us get useful information not only about the nature of the b b systems, but also on the existence of QGP. Since the discovery of quarkonium states, QCD sum rule technique as one of the most powerful nonperturbative tools to hadron physics [1, ] has played an important role in understanding the quarkonia spectrum. In order to find missing states one should know their physical properties to develop a successful search strategy. Clearly significant progress in understanding of quarkonium production cannot be reached without detailed measurements of the cross sections and fractions of the quarkonia. Thereby, form factors and decay widths of quarkonia become significant for completing quarkonia spectrum. Completing the bottomonium spectrum is a crucial validation of theoretical calculations and a test of our understanding of bottomonium states in the context of the quark model. Bottomonium states are considered as great laboratories to search for the properties od QCD at low energies. In this connection, we investigate the decay properties of tensor χ b (1P) meson as one of the important members of the bottomonia to the heavy B c meson in the present work. The B c meson with J P = is the only meson consists of two heavy quarks with different flavors. Yet, other possible B c states (the scalar, vector, axial-vector and tensor) have not been observed, however many new B c species are expected to be produced at the Large Hadron Collider (LHC) in the near future. Taking into account the two-gluon condensate corrections, the transition form factors of the semileptonic χ b (1P) B c lν decay channel are calculated within the three-point QCD sum rule. We use the values of transition form factors to estimate the decay rate of the transition under consideration at all lepton channels. The interpolating current of χ b (1P) with quantum numbers I G (J PC ) = + ( ++ ) contains derivatives with respect to space-time. So, we start our calculations in the coordinate space then we apply the Fourier transformation to go to the momentum space. To suppress the contributions of the higher states and continuum, we apply a double Borel transformation. 1

3 The article is organized as follows. We derive the QCD sum rules for the transition form factors in Section II. Last section is devoted to the numerical analysis of the obtained sum rules, estimation of the decay rates at all lepton channels and concluding remarks. QCD Sum Rules for χ b B c lν Transition Form Factors The semileptonic χ b B c lν decay is based on b c lν transition at quark level whose effective Hamiltonian can be written as H eff (b c lν l ) = G F V cb cγ µ (1 γ 5 )blγ µ (1 γ 5 )ν, (1) whereg F isthefermiweakcoupling constant andv cb iselement ofthecabibbo-kobayashi- Maskawa (CKM) mixing matrix. After sandwiching the effective Hamiltonian between the initial and final states, the amplitude of this transition is obtained in terms of transition matrix elements. These matrix elements will be parameterized in terms of transition form factors later. In order to start our calculations, we consider the three-point correlation function Π µαβ (p,p,q) = i d 4 xe ip x d 4 ye ip y T j Bc (y)jµ tr,v A ()j χ b(x), () where T is the time-ordering operator and j tr,v A µ () = c()γ µ (1 γ 5 )b() is the transition current. To proceed we also need the interpolating currents of the initial and final mesons in terms of the quark fields, which are given as j χ b αβ (x) = i [ b(x)γα Dβ (x)b(x)+ b(x)γ β Dα (x)b(x)] αβ (3) and j Bc (y) = b(y)γ 5 c(y), (4) where the covariant derivative Dβ (x) denotes the four-derivative with respect to x acting on two sides simultaneousely and it is defined as with and Dβ (x) = 1 [ Dβ (x) Dβ (x)], (5) Dβ (x) = β (x) i g λa A a β (x) (6) Dβ (x) = β (x)+i g λa A a β(x). (7) Here, λ a are the Gell-Mann matrices and A a β (x) denote the external gluon fields.

4 According to the method used, the correlation function in Eq.() is calculated in two different ways. In physical or phenomenological side we obtain it in terms of hadronic parameters such as masses and decay constants. In QCD or theoretical side we evaluate it in terms of QCD degrees of freedom like quark masses as well as quark and gluon condensates via operator product expansion (OPE). The QCD sum rules for form factors are obtained by equating the above representations to each other. After applying a double Borel transformation, the contributions of the higher states and continuum are suppressed. The hadronic side of the correlation function is obtained by inserting complete sets of intermediate states into Eq.(). After performing the four-integrals over x and y, we get µαβ (p,p,q) = j B c () B c (p ) B c (p ) jµ tr,v A χ b (p,ε) χ b (p,ε) j χ b αβ () (p m B c )(p m χ b ) Π PHYS +, where denotes the contribution of the higher states and continuum. To go further, we need to know the following matrix elements: and B c (p ) jµ tr,v χ b (p,ε) = h(q )ǫ µνθη ǫ νλ P λ P θ q η, (9) { } B c (p ) jµ tr,a χ b (p,ε) = i K(q )ǫ µν P ν +ǫ θη P θ P η [P µ b + (q )+q µ b (q )], (1) χ b (p,ε) j χ b αβ = f χb m 3 χ b ε αβ (11) j Bc B c (p ) = i f B c m B c m c +m b, (1) where h(q ), K(q ), b + (q ) and b (q ) are transition form factors, ǫ αβ is the polarization tensor associated with the χ b tensor meson; f χb and f Bc are leptonic decay constants of χ b and B c mesons, respectively, P µ = (p+p ) µ and q µ = (p p ) µ. Putting all matrix elements given in Eqs. (9), (1), (11) and (1) into Eq.(8), the final representation of the correlation function on physical side is obtained as where and Π PHYS µαβ (p,p,q) = { K(q )q α g βµ f χb f Bc m χb m B c 8(m b +m c )(p m χ b )(p m B c ) [ ] b (q )+ K(q ) q µ g αβ 3 [ ] b + (q )+K(q )( +4m χ 3 b ) P µ g αβ } i( 4m χ b )h(q )ε ληβµ P λ P η q α +other structures +..., (13) (8) = m B c +3m χ b q (14) = m 4 B c m B c (m χ b +q )+(m χ b q ). (15) 3

5 Note that we represented only the structures which we will use to find the corresponding form factors. Meanwhile, we used the following summation over the polarization tensors to obtain Eq. (13) ε λ µν ε λ αβ = 1 η µαη νβ + 1 η µβη να 1 3 η µνη αβ, (16) where λ η µν = g µν + p µp ν m χ b. (17) The next step is to calculate the QCD side of the correlation function in deep Euclidean region, where p and p via OPE. Placing the explicit expressions of the interpolating currents into the correlation function and contracting out all quark pairs via Wick s theorem, we obtain Π QCD µαβ (p,p,q) = i3 { Tr d 4 x d 4 ye ip x e ip y [ Sb ia (x y)γ 5Sc aj (y)γ µ(1 γ 5 ) ] Dβ (x)s ji b ( x)γ α +[β α] where S is the heavy quark propagator and it is given by [3] { SQ ai (x) = i d 4 ke ik (x) δ ai g sg ψϕ ai σ ψϕ ( k +m Q )+( k +m Q )σ ψϕ (π) 4 k m Q 4 (k m Q } ) + π 3 α sgg π δ k +m Q k aim Q +, (19) (k m Q )4 where Q = b or c quark. Replacing the explicit expression of the heavy quark propagators in Eq. (18) and applying integrals over x and y, we find the QCD side as ( ) ( ) Π QCD µαβ (p,p,q) = Π pert 1 (q )+Π nonpert 1 (q ) q α g βµ + Π pert (q )+Π nonpert (q ) q µ g βα ( ) ( ) + Π pert 3 (q )+Π nonpert 3 (q ) P µ g βα + Π pert 4 (q )+Π nonpert 4 (q ) ε λνβµ P λ P α q ν }, (18) + other structures. () Here Π pert i (q ) with i = 1,,3,4 are the perturbative parts which are expressed in terms of double dispersion integrals as Π pert i (q ) = ds ds ρ i (s,s,q ) +subtracted terms, (1) (s p )(s p ) where the spectral densities are defined as ρ i (s,s,q ) = 1 π Im[Πpert i ]. The spectral densities corresponding to four different structures shown in Eq. () are obtained as ρ 1 (s,s,q ) = 1 1 x dx dy [ ] 3 m c ( 3+4x+y)+m b ( 5+8x+4y), 16π 4

6 and ρ (s,s,q ) = 1 ρ 3 (s,s,q ) = 1 1 x dx dy dx 1 x dy [ ] 3 m c (3 4x y)+m b ( 1+4x+y), 8π [ ] 3 m c (1 y)+m b (1+y) () 8π ρ 4 (s,s,q ) =. (3) The function Π nonpert i (q ) in nonperturbative parts are calculated from a similiar manner and by considering the two-gluon condensates contributes. Having calculated both the physical and OPE sides of the correlation function, now, we match them to find QCD sum rules for form factors. In the Borel scheme we get K(q ) = 8(m b +m c ) χb f χb f Bc m B c m χb em /M e m Bc /M { s s } ds ds e s M e s M ρ 1 (s,s,q )Θ[L(s,s,q )]+ B M BM Π nonpert 1 (q ), (m b +m c) 4m b b (q ) = 1(m b +m c ) f χb f Bc m B c m χb em χb /M e m Bc /M { s s ds ds e s M e s M ρ (s,s,q )Θ[L(s,s,q )]+ B M BM Π nonpert (q ) 4m b (m b +m c) } ), K(q b + (q ) = 1(m b +m c ) f χb f Bc m B c m χb em χb /M e m Bc /M { s s ds ds e s M e s M ρ 3 (s,s,q )Θ[L(s,s,q )]+ B M BM Π nonpert 3 (q ) 4m b (m b +m c) } +4m χ b K(q ) and h(q 8(m b +m c ) ) = i χb f χb f Bc m B c m χb (4m χ b ) em /M { s 4m b ds s e m Bc /M (m b +m c) ds e s M e s M ρ 4 (s,s,q )Θ[L(s,s,q )]+ B M BM Π nonpert 4 (q ) where M and M are Borel mass parameters; and s and s are continuum thresholds in the initial and final channels. Here, Θ is the step function and L(s,s,q ) is given by L(s,s,q ) = s y(1 x y) sxy +m c (x+y 1) m b (x+y)+q x(1 x y). (6) 5 } (5) (4)

7 The functions B M BM Π nonpert i (q ) are written as B M BM Π nonpert i (q ) = α sgg π 1 dx e m b (1+ M x M )+M x M (m c q x) M x(1+( 1+ M M )x) f i (q ), (7) where f i (q ) are very lengthy functions and we do not present their explicit expressions here. 3 Numerical results To numerically analyze the sum rules obtained for the form factors, we use the meson masses from PDG [4]. Considering the fact that the results of sum rules considerably depend on the quark masses, decay constants and gluon condensate, we use the values of these parameters from different sources. For the quark masses we take into account all the pole values and those obtained at MS scheme from the Refs. [5 15]. For f Bc, we consider all the values predicted using different methods in Refs. [16 19]. In the case of f χb, we use the only value exists in the litreature, i.e. f χb = (.1±.7) []. For the gluon condensate, we also use its value from different sources [1, 1 8] calculated via different approaches. From the sum rules for the form factors it is also clear that they contain extra four auxiliary parameters, namely the Borel parameters M and M as well as continuum thresholds s and s. The general criteria is that the physical quantities like form factors should be independent of these parameters. Therefore, we need to determine their working regions such that the form factors weakly depend on these parameters. To find the Borel windows, we require that the higher states and continuum contributions are sufficiently suppressed and the perturbative parts exceed the nonperturbative contributions and the series of OPE converge. As a result we get the windows: 14 GeV M GeV and 8 GeV M 1GeV. The continuum threshold s ands arenotcompletely arbitrary but theyarerelatedtotheenergy ofthefirst excited stateswiththesamequantumnumbers as the interpolating currents of the initial and final channels. In this work the continuum thresholds are chosen in the intervals 14 GeV s 18 GeV and 43 GeV s 45 GeV. The dependence of form factors K and b +, as examples, on Borel parameters M and M at q = are plotted in figures 1 and. From these figures we see that form factors show overall weak dependence on the Borel mass parameters. The behaviors of the form factors K, b +, b and h in terms of q are shown in figures (3-6). In these figures, the red triangles show the QCD sum rule predictions, yellow-solid line denotes the prediction of fit function obtained using the central values of the input parameters and the green band shows the uncertainty due to errors of input parameters. Note that to obtain the central values, we consider the average values of input parameters discussed above, however, to calculate the uncertainties we consider all errors of these parameters from different sources previously quoted. As it is seen from these figures the sum rules results are truncated at some points. Hence, to enlarge the region to whole physical region we need to find some fit functions such that their results coincide well with the QCD sum rules predictions at reliable regions. For this reason we show the q dependence of form factors including both 6

8 K(q =) - - K(q =) s =14 GeV, s ' =43 GeV, M' =1 GeV -3 s =16 GeV, s ' =44 GeV, M' =1 GeV s =18 GeV, s ' =45 GeV, M' =1 GeV M (GeV ) -3 s =14 GeV, s ' =43 GeV, M =17 GeV s =16 GeV, s ' =44 GeV, M =17 GeV -3 s =18 GeV, s ' =45 GeV, M =17 GeV M' (GeV ) Figure 1: Left: K(q = ) as a function of the Borel mass parameter M at fixed values of s, s and M. Right: K(q = ) as a function of the Borel mass parameter M at fixed values of s, s and M. the sum rules and fit results in figures (3-6). Our numerical calculations reveal that the following fit function well defines the form factors under consideration: [ ( q f(q ) ( q ) ] ) = f exp a +b, (8) m χ b m χ b where the values of the parameters, f, a and b obtained at M = 17 GeV and M = 1 GeV for χ b B c lν transition are presented in table 1. Ourfinalpurposeinthissectionistoobtainthedecaywidthoftheχ b B c lν transition at all lepton channels. The differential decay width for this transition is obtained as f a b K(q ).871± ± ±.48 b (q ).134±.43 GeV 8.973± ±18.68 b + (q ).34±.97 GeV 1.54± ±17.55 h(q ) (.594±.83) 1 4 GeV 5.4± ±1.45 Table 1: Parameters appearing in the fit function of the form factors. dγ G F = V { cb dq 1 3 m 7 χ b π 3 q 6(m l q ) 3/ b (q ) m l q4 [ ] + b + (q ) (m B c m χ b ) m l +(m B c m χ b ) q (m B c +m χ b )q 4 +q 6 7

9 b + (q =) GeV , s =14 GeV, s ' =43 GeV, M' =1 GeV 1,,8 s =16 GeV, s ' =44 GeV, M' =1 GeV s =18 GeV, s ' =45 GeV, M' =1 GeV,6,4,8,6,4 b + (q =) GeV , 1, s =14 GeV, s ' =43 GeV, M =17 GeV,8 s =16 GeV, s ' =44 GeV, M =17 GeV s =18 GeV, s ' =45 GeV, M =17 GeV,6,4,8,6,4,,,,, , M (GeV ), , 1 M' (GeV ) Figure : Left: b + (q = ) as a function of the Borel mass parameter M at fixed values of s, s and M. Right: b + (q = ) as a function of the Borel mass parameter M at fixed values of s, s and M. ] + Re[K(q )b +(q )] [ q 4 +m B c (m l +q ) m χ b (m l +q ) Re[b (q )b +(q )] m lq (m B c m χ b ) + K(q ) [m 4 B c (m l +q )+m 4 χ b (m l +q )+q 4 (m l +q ) m B c ] (m χ b +q )(m l +q )+m χ b q (m l +5q ) + 3 h(q ) m χ b q (m l +q ) Re[K(q )b (q )] m l }. q (9) After performing integration over q in Eq.(9) in the interval m l q (m χb m Bc ), we obtain the total decay widths as presented in table for different leptons. The errors belong Γ(GeV) χ b B c ēν e χ b B c µν µ χ b B c τν τ (1.54±.56) 1 13 (1.41±.5) 1 13 (.398±1.175) 1 14 Table : Numerical results of decay widths at different lepton channels. to the uncertainties coming from the determination of the working regions for auxiliary parameters as well as those of the other input parameters. The orders of decay rates at all lepton channels show that these transitions are accessible at LHC in near future. 8

10 -,5 -, K(q ) -1,5-1, -,5, q (GeV ) Figure 3: K(q ) as a function of q at s = 16 GeV, s = 44 GeV, M = 17 GeV and M = 1 GeV. The red triangles show the QCD sum rule predictions, yellow-solid line denotes the prediction of fit function obtained using the central values of the input parameters and the green band shows the uncertainty due to errors of input parameters. In summary, we have calculated the transition form factors for the semileptonic χ b B c lν transition using QCD sum rule technique. We took into account the two-gluon condensate contributions as nonperturbative effects. We used these form factors to estimate the order of decay widths at all lepton channels. The order of decay width reveal that these transitions can be seen at LHC in near future. Any comparison of the experimental results with our predictions can provide us with essential knowledge on the nature of the tensor χ b (1P) state. 4 Conflict of Interests The authors declare that there is no conflict of interest regarding the publication of this paper. References [1] M. A. Shifman, A. I. Vainshtein, V. I. Zakharov, QCD and resonance physics theoretical foundations, Nucl. Phys. B 147, 385 (1979). [] M. A. Shifman, A. I. Vainshtein, V. I. Zakharov, QCD and resonance physics applications, Nucl. Phys. B 147, 448 (1979). 9

11 -,8 h(q ) GeV - -,6 -,4 -, q (GeV ) Figure 4: The same as figure 3 but for h(q ). 5 4 b + (q ) GeV q (GeV ) Figure 5: The same as figure 3 but for b + (q ). 1

12 -, -1,5 b - (q ) GeV - -1, -,5, q (GeV ) Figure 6: The same as figure 3 but for b (q ). [3] L. J. Reinders, H. Rubinstein, S. Yazaki, Hadron properties from QCD sum rules, Phys. Rept. 17, 1 (1985). [4] K. A. Olive et al. (Particle Data Group), Review of particle physics, Chin. Phys. C 38, 91 (14). [5] H. Abramowicz et al.[h1 and ZEUS Collaborations], Combination and QCD Analysis of Charm Production Cross Section Measurements in Deep-Inelastic ep Scattering at HERA, Eur.Phys.J. C 73, 311 (13), arxiv: [hep-ex]. [6] H. Abramowicz et al. [ZEUS Collaboration], Measurement of beauty and charm production in deep inelastic scattering at HERA and measurement of the beauty-quark mass, JHEP 149, 17 (14), arxiv: [hep-ex]. [7] S. Narison, Summary on m c,b (m c,b ) and precise f D(s),B(s) from heavy-light QCD spectral sum rules, Nucl. Phys. B - Proceedings Supplements 34, 187 (13), 16th Inter. Conf. in QCD, arxiv:19.95 [hep-ph]. [8] K. G. Chetyrkin et al., Charm and Bottom Quark Masses: an Update, Phys. Rev. D 8, 741 (9), arxiv:97.11 [hep-ph]. [9] A. J. Lee et al., Mass of the b-quark from lattice NRQCD and lattice perturbation theory, Phys. Rev. D 87, 7418 (13), arxiv: [hep-lat]. [1] N. Carrasco et al. [European Twisted Mass Collaboration], Up, down, strange and charm quark masses with N f = +1+1 twisted mass lattice QCD, Nucl. Phys. B 887, 19 (14), arxiv: [hep-lat]. 11

13 [11] Y. Kiyo, G. Mishima and Y. Sumino, Determination of m c and m b from quarkonium 1S energy levels in perturbative QCD, Phys. Lett. B 75, 1 (16), arxiv: [hep-ph]. [1] B. Dehnadi, A.H. Hoang and V. Mateu, Bottom and Charm Mass Determinations with a Convergence Test, JHEP 158, 155 (15), arxiv: [hep-ph]. [13] A. Ali, A.Y. Parkhomenko, A.V. Rusov, Precise Calculation of the Dilepton Invariant- Mass Spectrum and the Decay Rate in B ± π ± µ + µ in the SM, Phys. Rev. D 89, 941 (14), arxiv: [hep-ph]. [14] K. G. Chetyrkin, M. Steinhauser, Short-Distance Mass of a Heavy Quark at Order α 3 s, Phys. Rev. Lett. 83, 41 (1999), [15] P. Colangelo, A. Khodjamirian, QCD sum rules, a modern perspective, At the Frontier of Particle Physics/Handbook of QCD, edited by M. Shifman (World Scientific, Singapore), 3, 1495, (1), arxiv:hep-ph/1175. [16] M. J. Baker, J. Bordes, C. A. Dominguez, J. Penarrocha, K. Schilcher, B meson decay constants f Bc, f Bs and f B from QCD sum rules, JHEP 147, 3 (14) arxiv: [hep-ph]. [17] B. Colquhoun et al., B-meson decay constants: a more complete picture from full lattice QCD, Phys. Rev. D 91, (15), arxiv: [hep-lat]. [18] S. Narison, Improved f D (s),f B (s) and f Bc from QCD Laplace sum rules, Int. J. Mod. Phys. A 3,, (15), arxiv: [hep-ph]. [19] E. Bagan, H. G. Dosch, P. Gosdzinsky, S. Narison and J. M. Richard, Hadrons with Charm and Beauty, Z. Phys. C 64, 57 (1994), arxiv:hep-ph/9438. [] T. M. Aliev, K. Azizi, M. Savci, Heavy χ Q tensor mesons in QCD, Phys. Lett. B 69, 164 (1) arxiv:1.767 [hep-ph]. [1] B. L. Ioffe, QCD at low energies, Prog. Part. Nucl. Phys. 56, 3 (6), arxiv:hepph/ [] C. A. Dominguez, L. A. Hernandez and K. Schilcher, Determination of the gluon condensate from data in the charm-quark region, JHEP 157, 11 (15), arxiv: [hep-ph]. [3] R. Horsley et al., Wilson loops to th order numerical stochastic perturbation theory, Phys. Rev. D 86, 545 (1), arxiv: [hep-lat]. [4] B. Chakraborty et al., High-precision quark masses and QCD coupling from n f = 4 lattice QCD, Phys. Rev. D 91, 5458 (15), arxiv: [hep-lat]. [5] C. A. Dominguez, L. A. Hernandez, K. Schilcher and H. Spiesberger, Chiral sum rules and vacuum condensates from tau-lepton decay data, JHEP 153, 53 (15), arxiv: [hep-ph]. 1

14 [6] B. V. Geshkenbein, Calculation of gluon and four-quark condensates from the operator expansion, Phys. Rev. D 7, 747 (4), arxiv:hep-ph/391. [7] H. G. Dosch, M. Jamin, S. Narison, Baryon masses and flavour symmetry breaking of chiral condensates, Phys. Lett. B, 51 (1989). [8] V. M. Belyaev, B. L. Ioffe, Determination of baryon and baryonic masses from QCD sum sules. Strange baryons, Sov. Phys. JETP 57, 716 (1983). 13