arxiv: v2 [hep-lat] 24 Nov 2016

Transcription

1 esons in Lattice QCD with N = twisted-ass erions arxiv: v2 [hep-at] 24 Nov 206 Vittorio Lubicz Dip. di Mateatica e Fisica, Università Roa Tre & INFN, Sezione di Roa Tre, Roe, Itay E-ai: ubicz@is.uniroa3.it Aurora Meis Dipartiento di Mateatica e Fisica, Università Roa Tre, Roe, Itay E-ai: e.aurora.6@gai.co Sivano Siua Istituto Nazionae di Fisica Nuceare, Sezione di Roa Tre, Roe, Itay E-ai: siua@roa3.inn.it or the ETM Coaboration We present a attice cacuation o the decay constants and asses o D (s) and B (s) esons using the gauge conigurations produced by the European Twisted Mass Coaboration (ETMC) with N = dynaica quarks and at three vaues o the attice spacing a Pion asses are siuated in the range π MeV, whie the strange and char quark asses are cose to their physica vaues. We coputed the ratios o vector to pseoscaar decay constants or asses or various vaues o the heavy-quark ass h in the range 0.7c phys h 3c phys. In order to reach the physica b-quark ass, we expoited the HQET prediction that, in the static iit o ininite heavy-quark ass, a the considered ratios are equa to one. We obtain: D / D =.078(36), D / D =.0769(79), D s / Ds =.087(20), D s Ds =.075(56), B / B = 8(22), B / B =.0078(5), B s / Bs = 74(0) and B s / Bs =.0083(0). Cobining the with the corresponding experienta asses ro the PDG and the pseoscaar decay constants cacuated by ETMC, we get: D = 223.5(8.4) MeV, D = 203(4) MeV, D s = 268.8(6.6) MeV, D s = 26() MeV, B = 85.9(7.2) MeV, B = (7.6) MeV, B s = 223.(5.4) MeV and B s = 54.36(5.3) MeV. 34th annua Internationa Syposiu on Lattice Fied Theory Juy 206 University o Southapton, UK c Copyright owned by the author(s) under the ters o the Creative Coons Attribution-NonCoercia-NoDerivatives 4.0 Internationa License (CC BY-NC-ND 4.0).

2 . Introduction The decay constants o D (s) and B (s) esons are iportant ingredients in the phenoenoogica description o various processes, ike seieptonic and non-eptonic decays o heavy hadrons. It is we known that in the iit o ininite heavy-quark ass the Heavy Quark Eective Theory (HQET) predicts that the ratios o vector (V) to pseoscaar (PS) decay constants or asses are equa to one, i.e. i h (M H /M H ) = and i h ( H / H ) =. When the heavy quark is either the char or the beauty, the spin-avor syetry is broken and the above ratios deviate ro one because o power corrections in / h. Ti now there are ony ew attice cacuations o the D (s) and B (s) decay constants using gauge conigurations with N = 2 [, 2] and N = 2+(+) [3, 4] dynaica quarks. These resuts exhibit a surprisingy non-negigibe dependence on N. In this contribution we present the resuts obtained or the V to PS ratios o decay constants and asses in the char and beauty sectors using the gauge ensebes generated by the European Twisted Mass Coaboration (ETMC) with N = 2++ dynaica quarks [5, 6, 7]. In the ETMC set-up the guon interactions are described by the Iwasaki action, whie the erions are reguarised in the axiay twisted-ass (Mt) Wison attice oruation. We considered three vaues o the attice spacing, naey a = (36), 0.085(30) and 0.069(8), with the owest siuated pion ass being equa to 20 MeV. The vaence quark asses are chosen to be in the ranges: 3 phys 2 phys, 0.7phys s s.2s phys and 0.7c phys c c phys. To extrapoate up to the b-quark sector we have aso considered higher vaues o the vaence heavyquark ass in the range c phys h 3c phys 0.7 phys b. The attice scae was deterined using the experienta vaue o π + [8], whie the physica up/down, strange, char and botto quark asses were obtained [8, 9] by using the experienta vaues or π, K, D and B, respectivey. In Re. [8] eight branches o the anaysis were adopted. They dier in: i) the continuu extrapoation adopting or the scae paraeter either the Soer paraeter r 0 or the ass o a ictitious PS eson ade up o strange(char)-ike quarks; ii) the chira extrapoation perored with itting unctions chosen to be either a poynoia expansion or a Chira Perturbation Theory (ChPT) ansatz in the ight-quark ass; iii) the choice between the ethods M and M2, which dier by O(a 2 ) eects, used to deterine in the RI -MOM schee the ass renoraization constant (RC) Z = /Z P. In the present anaysis we ade use o the input paraeters corresponding to each o the eight branches o Re. [8]. 2. Extraction o asses and decay constants The decay constants o V and PS esons are deined in ters o the atrix eeents 0 hγ µ H ( p,λ) = H H ε λ µ, (2.) 0 hγ 5 H ( p) = 0 µ ( hγµ γ 5 ) H ( p /( h + ) = H 2 H /( h + ), (2.2) where h and are the heavy- and ight-quark asses with h = {c,b} and = {,s}, and ε λ µ is the vector eson poarization. Ground-state asses and decay constants can be deterined by stying two-point correation unctions at arge tie distances, viz. C V (t) = 3 i, x V i ( x,t)v i (0,0) t t in i 0 V i (0) H (λ) 2 cosh[ H (T /2 t)] e H T /2, (2.3) 3 H

3 C P (t) = x P( x,t)p (0,0) 0 P(0) H 2 cosh[ H (T /2 t)] e H T /2, (2.4) t t in H where t in stands or the iniu tie at which the ground state can be considered we isoated. In Eq. (2.3) V i Z A hγ i is the oca vector current, which in our Mt setup renoraizes utipicativey with the RC Z A, whie in Eq. (2.4) P ( h + )Z P hγ 5 = (µ h + µ )hγ 5, where µ h and µ are bare quark asses, is the pseoscaar interpoating ied, which in our Mt setup is renoraization group invariant and does not require any RC. The eson asses H and H are extracted ro the pateaux o the eective ass at arge t t in, whie the correation unctions (2.3) and (2.4) or t t in contain the required atrix eeents. In Eqs. ( ) we considered oca source and sink operators, but we anayzed aso the whoe set o our correation unctions given by the cobinations o oca interpoating operators with those obtained ro a Gaussian searing procedure in both the sink and the source, naey and CSS P,V, where L and S denote oca and seared operators, respectivey. It is straightorward to check that the required oca atrix eeents in Eqs. (2.-2.2) can be extracted ro the LL correation unctions as we as ro an appropriate cobination o the SL and SS ones, that is CP,V SL (t)/ C SS C LL P,V,CLS P,V,CSL P,V P,V (t). For the reasons expained in the Introduction we have considered the oowing ratios R H = H / H and R H = H / H. (2.5) In Fig. we show an exape o the above quantities by coparing the extraction ro Gaussianseared and/or oca correation unctions. The searing techniques aows pateaux to start at earier tie distances, and the SL correation unctions exhibit the best signa to noise ratio. The vaue o t in corresponds to the saest tie distance where the eective asses obtained ro SL and SS correation unctions intercept each other. R H β =.90, L = 32, aµ = aµ ~ s phys, aµh ~ c phys R H β =.90, L = 32, aµ = aµ ~ s phys, aµh ~ c phys Seared-Seared Loca-Loca Seared-Loca t/a oca seared cobination t/a Figure : The ratios R H and R H (see Eq. (2.5)) using ony oca (red points) or aso Gaussian seared operators (bue points). The vaue o t in is shown as the vertica dot-dashed ine. The grey bands are the ass and the decay constant, obtained as constant its in the pateaux regions. 2. The D (s)-eson asses and decay constants We peror a sooth interpoation o the attice data or the ratios R H and R H to the vaues o the physica strange and char quark asses s phys (MS, 2 GeV) = 99.6(4.3) MeV and 2

4 c phys (MS, 2 GeV) = 76(39) GeV [8]. The dependence o R D (s) and R D (s) on the renoraized up/down quark ass = aµ /(az P ) and the attice spacing a is investigated by peroring a cobined chira and continuu extrapoation, based on a poynoia expansion o the or R it D (s) (,a) = P 0 + P + P 2 a 2 + P P 4 a 4, (2.6) where we have taken into account that or our Mt setup the autoatic O(a)-iproveent ipies that discretization eects invove ony even powers o the attice spacing. The resuts obtained with quadratic 2 and quartic a4 ters have not been inced in the ina average (which thereore corresponds to P 3 = P 4 = 0), but they have been considered to estiate the uncertainty reated to the chira and continuu extrapoation, respectivey. The atter ones are shown in Fig. 2, where the physica point corresponds to phys (MS, 2 GeV) = 3.70(7) MeV [8] R D (, a) β =.90 L = 24 β =.90 L = 32 β =.95 L = 24 β =.95 L = 32 β = 2.0 L = 48 Physica point R D (, a) β =.90 L = 24 β =.90 L = 32 β =.95 L = 24 β =.95 L = 32 β = 2.0 L = 48 Physica point Figure 2: Chira and continuu extrapoations o R D and R D based on the poynoia it (2.6) with P 3 = P 4 = 0. The green points represent the vaues at the physica point phys (MS, 2 GeV) = 3.70(7) MeV [8]. Siiar resuts hod as we in the case o the ratios R D s and R D s. In this way at the physica point we get D / D =.0769(7) stat (30) input (3) tin (8) disc (5) chir [79], (2.7) D s / Ds =.075(49) stat (27) input (8) disc (4) tin (2) chir [56], (2.8) D / D =.078(3) stat (9) chir (8) disc (6) tin (5) input [36], (2.9) D s / Ds =.087(6) stat (7) disc (6) input (6) tin (5) chir [20], (2.0) where tota uncertainty (in the square brakets) is the su in quadrature o the statistica and a the systeatic uncertainties, which have been written in order o reevance or each o the ratios. The various systeatic uncertainties are estiated in the oowing way: i) the uncertainty abeed tin is coputed by repeating the anaysis with a vaue o t in shited by two units and taking the ha dierence with the ina estiates; ii) the chira and discretization uncertainties, abeed respectivey as chir and disc, are obtained by considering either P 3 0 or P 4 0 in Eq. (2.6) and taking again ha o the dierence with the ina resuts; iii) the uncertainty abeed input is given by the spread o the resuts over the eight branches o the input paraeters o Re. [8]. By cobining R D (s) with the experienta vaues o the D (s) -eson asses [0] we obtain D = 203 (4) MeV and D s = 26 () MeV, (2.) 3

5 that copare we with the experienta eson asses exp D 22.(4) MeV [0]. = (5) MeV and exp D = s As or the decay constants, existing attice cacuations or R D (s) have been carried out ony with N = 2 + and N = 2 dynaica quarks. The N = 2 + estiate D s / Ds = 0(2) [3] is in good agreeent with our resut, whie the N = 2 resuts D / D =.208(27) [2] and D s / Ds =.26(3)[] are 0% arger than our predictions. Using the pseoscaar decay constants cacuated by ETMC in Re. [] we get D = (8.7) MeV and D s = (6.5) MeV. (2.2) 2.2 The B (s)-eson asses and decay constants We have coputed the ratios R (k) H and R (k) H or a series o asses { (k) h } c with k =,...,8. The resuts are extrapoated to the chira and continuu iits, as shown in Fig. 3. R H (3) (, a) β =.90 L = 24 β =.90 L = 32 β =.95 L = 24 β =.95 L = 32 β = 2.0 L = 48 Continuu iit R H (3) (, a) Figure 3: Chira and continuu extrapoations o R (k) H β =.90 L = 24 β =.90 L = 32 β =.95 L = 24 β =.95 L = 32 β = 2.0 L = 48 Continuu iit and R (k) H (2.6) with P 3 = P 4 = 0. Siiar resuts hod as we in the case o R (k) H s or k = 3, based on the poynoia it and R (k) H s. The HQET predicts that the ratios R H and R H R H /C W ( h ), where C W ( h ) is the perturbative atching correction between u QCD and HQET (coputed in Re. [2] up to next-to-nexteading order), are equa to one in the static heavy-quark iit, viz. i h R H = and i R h H = i R h H /C W ( h ) =. (2.3) Thus, we peror correated poynoia its in / h iposing the static iit constraint, naey R H it phys = + D 2/ 2 h + D 3 / 3 h + D 4/ 4 h, R H it phys = + D / h + D 2 / 2 h + D 3 / 3 h, (2.4) where we have taken into account that, according to HQET, the inear ter is absent in the case o the ass ratio (i.e., D = 0). In Fig. 4 the interpoations o the various ratios in the inverse heavy-quark ass are shown together with the resuts at the b-quark physica point obtained using the vaue phys b ro Re. [9]. 4

6 R H phys / Cw (h) R Hs phys / Cw (h) R H phys / Cw (h) /h (GeV - ) /h (GeV - ) R Hs phys / Cw (h) /h (GeV - ) /h (GeV - ) Figure 4: The dependence o R H (s) phys and R H (s) phys, extrapoated to the chira and continuu iits, on the inverse heavy-quark ass / h (MS, 2 GeV). The its are based on Eqs. (2.4) and the correation atrix aong the data points is accounted or. The vertica bands correspond to / phys b ro Re. [9]. Our ina resuts or B (s) esons are B / B =.0078(8) stat (8) chir (7) tin (5) disc (2) input [4], (2.5) B s / Bs =.0083(6) stat (7) chir (6) disc (3) tin (2) input [], (2.6) B / B = 8(8) stat (0) disc (6) chir (5) tin (2) input [22], (2.7) B s / Bs = 74(7) stat (6) disc (3) tin (2) input () chir [0], (2.8) where the error bget accounts or the sae sources o uncertainties aready considered or the char sector in Sec. 2.. Mass ratios can be cobined with the experienta vaues o B (s) -eson asses [0] to obtain B = (7.6) MeV and B s = 54.8 (6.2) MeV, (2.9) that copare nicey with the experienta vaues exp B = (32) MeV and exp B = 545.4(.6) s MeV [0]. As or the decay constant ratios, we can copare our resuts with a recent coputation [4] obtained ro N = 2++ siuations (ike the ones considered in this work), B / B = 4(26) and B s / Bs = 3(23), as we as with a recent deterination based on the QCD su rue approach [3], B / B = 44(23) and B s / Bs = 47(30). A these estiates are nicey consistent with our resuts. On the contrary we ind again a 0% dierence with the N = 2 deterination B / B =.05(7) ro Re. [2]. 5

7 Eventuay, cobining our resuts or the ratios with the pseoscaar decay constants cacuated by ETMC in Re. [9] yieds 3. Concusions B = 86.4 (7.) MeV and B s = 223. (5.6) MeV. (2.20) We have coputed the asses and the decay constants o vector heavy-ight esons using ETMC gauge conigurations with N = dynaica quarks. Our resuts reproduce very we the experienta vaues o both D (s) - and B (s)-eson asses. We have ound that D (s) / D(s) > and B (s) / B(s) < with a spin-avor syetry breaking eect o +8% in the char sector and 4% in the beauty sector. Our resuts or the decay constant ratio exhibit a tension with the corresponding attice deterinations obtained by ETMC at N = 2 [, 2], whie they are consistent with the indings o Res. [3, 4] obtained by HPQCD with N = 2 + (+) dynaica quarks. Since our present anaysis oow aost the sae steps o the previous ETMC anayses at N = 2, the observed 0% tension ay be due to a dependence on the nuber o sea quarks, and in particuar to the incusion o the strange quark. The possibiity that the observed dierence can be attributed to a quenching eect o the strange quark is a quite interesting issue, because its size wo be arger than what typicay expected. Further investigations at dierent N vaues are thereore required. Acknowedgents We grateuy acknowedge the CPU tie provided by PRACE under the project PRA067 on the BG/Q syste Juqueen at JSC (Gerany) and by CINECA under the speciic initiative INFN-LQCD23 on the BG/Q syste Feri at CINECA (Itay). Reerences [] D. Becirevic et a., JHEP 202 (202) 042 [arxiv: [hep-at]]. [2] D. Becirevic et a., arxiv: [hep-ph]. [3] G.C. Donad et a. [HPQCD], Phys. Rev. Lett. 2 (204) [arxiv: [hep-at]]. [4] B. Coquhoun et a. [HPQCD], Phys. Rev. D 9 (205), 4509 [arxiv: [hep-at]]. [5] R. Baron et a. [ETMC], JHEP 006 (200) [arxiv: [hep-at]]. [6] R. Baron et a. [ETMC], Coput. Phys. Coun. 82 (20) 299 [arxiv: [hep-at]]. [7] R. Baron et a. [ETMC], PoS LATTICE 200 (200) 23 [arxiv:0.058 [hep-at]]. [8] N. Carrasco et a. [ETMC], Nuc. Phys. B 887 (204) 9 [arxiv: [hep-at]]. [9] A. Bussone et a. [ETMC], Phys. Rev. D 93 (206) no., 4505 [arxiv: [hep-at]]. [0] K. A. Oive et a., Chin. Phys. C 40 (206) no.0, [] N. Carrasco et a. [ETMC], Phys. Rev. D 9 (205) 5, [arxiv: [hep-at]]. [2] D. J. Broadhurst and A. G. Grozin, Phys. Rev. D 52 (995) 4082 [hep-ph/940240]. [3] W. Lucha et a., Phys. Rev. D 9 (205) no., 6009 [arxiv: [hep-ph]]. 6