arxiv: v1 [hep-lat] 13 Nov 2013

Transcription

1 Peudocalar decay contant / f π, f D and f D with N f = ETMC configuration arxiv: v1 [hep-lat] 13 Nov 2013 P. Dimopoulo (a,b), R. Frezzotti (b,c), P. Lami (d,e), V. Lubicz (d,e), E. Picca (d,e), (d,e), G.C. Roi (b,c), F. Sanfilippo ( f ), S. Simula (e), C. Tarantino (d,e) (a) Centro Fermi - Mueo Storico della Fiica e Centro Studi e Ricerche E. Fermi, Rome, Italy (b) Dipartimento di Fiica, Univerità di Roma Tor Vergata, Rome, Italy. dimopoulo@roma2.infn.it, frezzotti@roma2.infn.it, roig@roma2.infn.it (c) INFN, Sezione di Tor Vergata", Rome, Italy (d) Dipartimento di Matematica e Fiica, Univerità Roma Tre, Rome, Italy. lamipaolo@gmail.com, lubicz@fi.uniroma3.it, e.picca88@gmail.coorenzo.riggio@gmail.com, tarantino@fi.uniroma3.it (e) INFN, Sezione di Roma Tre, Rome, Italy. imula@roma3.infn.it ( f ) Laboratoire de Phyique Théorique (Bât. 210), Univerité Pari Sud, F Oray-Cedex, France. fr.anfilippo@gmail.com For the ETM Collaboration We preent a lattice QCD calculation of the peudocalar decay contant, f D and f D performed by the European Twited Ma Collaboration with N f = dynamical fermion. We imulated at three different value of the lattice pacing, the mallet being approximately 0.06 f m, and with pion mae a mall a 210MeV. Our main reult are: +/ f π + = 1.183(17), + = 154.4(2.1)MeV, f D = 242.1(8.3)MeV, f D = 201.9(8.0)MeV, f D / f D = 1.199(25) and ( f D / f D )/( / f π ) = 1.005(15). 31t International Sympoium on Lattice Field Theory LATTICE 2013 July 29 - Augut 3, 2013 Mainz, Germany Speaker. c Copyright owned by the author() under the term of the Creative Common Attribution-NonCommercial-ShareAlike Licence.

2 1. Introduction and imulation detail An accurate determination of the peudocalar (PS) meon decay contant i a crucial ingredient for the determination of the CKM matrix element, and in turn for teting the Standard Model (SM) and earching new phyic (NP). In thi contribution we preent a lattice QCD calculation of the, f D and f D decay contant uing the enemble of gauge configuration produced by the European Twited Ma (ETM) Collaboration with four flavor of dynamical quark (N f = 2+1+1), which include in the ea, beide two light ma degenerate quark, alo the trange and the charm quark. The imulation were carried out at three different value of the invere bare lattice coupling β that allow for a controlled extrapolation to the continuuimit, and at different lattice volume. For each enemble we ued a ubet of well-eparated trajectorie to avoid autocorrelation. We imulated the pure gauge Iwaaki action [1] for gluon, and the Wilon Twited Ma Action [2] for ea quark, which at maximal twit allow for an automatic O(a)-improvement [3]. To avoid mixing in the trange and charm ector we adopted a non-unitary etup in which valence quark are imulated for each flavor uing the Oterwalder-Seiler action [4]. In order to minimize dicretization effect in the determination of the peudocalar meon mae, the value of the Wilon parameter r are alway choen o that the two valence quark in a meon have oppoite value of r. At each lattice pacing different value of light ea quark mae have been conidered. The light valence and ea quark mae are alway taken to be equal. In the light, trange and charm ector the quark mae were imulated in the range 0.1 m phy 1.2 m phy and 0.7 m phy c µ c 2.0 mc phy, repectively. µ l 0.5 m phy, 0.7 m phy µ We tudied the dependence of the PS meon decay contant on the light quark ma fitting imultaneouly the data at different lattice pacing and volume. The lattice pacing are: a = (36), (30), (18) fm [5], and the lattice volume goe from 2 to 3 fm. The pion mae, extrapolated to the continuum and infinite volume limit, range from 210 to 450MeV. Within our analye we ued: i) the reult for r 0 /a (r 0 i the Sommer parameter [6]) obtained from either a linear or a quadratic extrapolation to the chiral limit of the value computed in [7, 8], and ii) the quark ma renormalization contant Z m = 1/Z P computed in [9] in the RI-MOM cheme uing two different method, labelled a M1 and M2, which differ by O(a 2 ) effect. 2. Calculation of In order to compute the kaon decay contant, we firt performed a mall interpolation of our lattice data to the phyical value of the trange quark ma m determined in [5]. Then we analyzed the dependence of the decay contant a a function of the (renormalized) light-quark ma and the lattice pacing and extrapolated it to the phyical point. The SU(2) ChPT prediction at the next-to-leading order (NLO) for read a follow: ( = P ) 4 ξ l logξ l + P 2 ξ l + P 3 a 2 K FSE f, (2.1) where ξ l = 2B 0 /16π 2 f 2 0 with B 0 and f 0 being the low-energy contant (LEC) entering the LO chiral Lagrangian. The term proportional to a 2 in Eq. (2.1) account for leading dicretization 2

3 effect. The factor K FSE f repreent the correction for the finite ize effect (FSE) in the kaon decay contant, computed following the work of Colangelo, Dürr and Haefeli (CDH) [10]. r Continuum Limit at m ud Continuum Limit at m u /M Continuum Limit at m ud Continuum Limit at m u 0.22 Figure 1: Chiral and continuum extrapolation of r 0 (left) and /M (right) baed on the NLO ChPT fit given in Eq. (2.1). Lattice data have been corrected for FSE uing the CDH approach [10]. The chiral and continuum extrapolation of i hown in Fig. 1 in unit of either r 0 or the ma M of a fictitiou PS meon made of two trange-like valence quark with ma fixed at r 0 m = The impact of dicretization effect uing r 0 a the caling variable i at the level of 3%. In order to keep the extrapolation to the continuuimit well under control we repeated the analyi adopting a different choice for the caling variable, namely intead of r 0 we introduced the ma M, which i affected by cutoff effect imilar to the one of a K-meon without however any ignificant dependence on the light quark ma.. Thu, we performed the chiral and continuum extrapolation alo for the ratio /M. The comparion between the analye performed in unit of r 0 and M clearly how that, when M i choen a the caling variable, the dicretization effect on change from 3% down to 1.5%. For the chiral extrapolation we adopted both the NLO ChPT prediction (2.1) and a polynomial formula in ξ l. Notice in Fig. 1 that after taking the continuuimit the kaon decay contant ha been extrapolated at two different value of the light quark ma. The firt one, labelled by a magenta diamond, i the reult obtained for in the iopin ymmetric limit and correpond to the extrapolation at the average up/down quark ma m ud = 3.70(17)MeV [5]. The econd one, repreented by a violet diamond, i the value extrapolated to the up quark ma m u = 2.47(11)MeV, obtained adopting the ma ratio m u /m d = 0.5 from [11], and it correpond to the quantity + corrected for leading trong iopin breaking effect. In fact, it can be hown that the firt order correction due to the QCD iopin breaking effect depend only on the valence quark. The ea quark effect enter proportionally to the quare of the up/down ma difference, (m d m u ) 2, an effect which i well below the preent preciion. Thu, in order to correct for leading iopin breaking effect it i ufficient to extrapolate to the up quark ma. Electromagnetic iopin breaking correction are a much more challenging iue and have not been conidered in the preent tudy. The variou ource of ytematic uncertaintie are etimated a follow. The difference of the reult obtained uing either r 0 or M a the caling variable i taken a a ytematic uncertainty on aociated with dicretization effect. A ytematic uncertainty related to the chiral extrapol- 3

4 ation i obtained by comparing the reult of the NLO ChPT fit with thoe of the polynomial fit. A for FSE we compared the reult obtained by applying the CDH correction with the one obtained without correcting for FSE. The two method M1 and M2, ued in [9] to calculate the renormalization contant Z P in the RI-MOM cheme, allow u to etimate a ytematic uncertainty due to the ma renormalization contant. Finally, the error on our determination of the trange quark ma repreent another ource of uncertainty, and it ha been included in the tat+fit error, which alo combine the tatitical uncertainty and the error aociated with the fitting procedure (dominated by the continuuimit and by the eparation between the phyical point and the lightet imulated pion ma). Combining the variou ource of uncertainty we find our final reult for +: + = 154.4(1.8) tat+ f it (0.5) Chiral (1.0) Dic. (0.4) FSE (0.2) ZP MeV = 154.4(2.1)MeV. (2.2) Dividing the reult (2.2) by the experimental value of the pion decay contant, which ha been ued a input to et the lattice cale [5], we obtain for the ratio +/ f π + the value +/ f π + = 1.183(14) tat+ f it (4) Chiral (8) Dic. (4) FSE (1) ZP (2) fπ + = 1.183(17), (2.3) which can be compared with the FLAG average +/ f π + = 1.205(18) at N f = 2, +/ f π + = 1.192(5) at N f = and +/ f π + = 1.195(5) at N f = [12]. In the iopin ymmetric limit we get for the value = 155.6(1.6) tat+ f it (0.5) Chiral (1.1) Dic. (0.4) FSE (0.2) ZP MeV = 155.6(2.1)MeV, (2.4) which can be compared with the FLAG average = 158.1(2.5)MeV at N f = 2 and = 156.3(0.8)MeV at N f = [12], and for the ratio / f π the reult / f π = 1.193(13) tat+ f it (4) Chiral (8) Dic. (4) FSE (1) ZP (2) fπ = 1.193(16). (2.5) Had we neglected at all dicretization effect in the kaon decay contant and had we limited ourelve only to the gauge enemble at β = 1.95 and 2.10, the reult for / f π would have been larger by 2.5% getting very cloe to the finding / f π = 1.224(13) obtained in [13]. 3. Calculation of f D, f D and f D / f D The phyical value of f D, f D and f D / f D have been determined by analyzing Φ D = f D MD and the double ratio ( f D / f D )/( / f π ) a function of m c, m, and a 2. We firt performed an interpolation of lattice data for f D, f D and f D / f D to the phyical trange and charm quark mae determined in [5]. Then the dependence of Φ D on the light-quark ma and on the lattice pacing turned out to be decribed well by the imple expreion Φ D = P 1 (1 + P 2 + P 3 m 2 l + P 4 a 2 ). (3.1) Uing the experimental value M D = 1.969GeV allowed u to determine the phyical value of f D. A in the cae of the kaon, the lattice data for Φ D are converted in unit of either r 0 or the ma M c of a fictitiou PS meon, made by one trange-like and one charm-like valence quark 4

5 Φ D (r 0 ) 3/ Cont. Limit Linear Fit Cont. Limit Quadratic Fit Φ D /(M c ) 3/ Cont. Limit Linear Fit Cont. Limit Quadratic Fit 0.13 Figure 2: Chiral and continuum extrapolation of r 3/2 0 Φ D (left) and Φ D /M 3/2 c (right) baed on Eq. (3.1). with mae fixed at r 0 m = 0.22 and r 0 m c = 2.4, repectively. Such a reference ma M c i expected to have dicretization effect cloer to the one of M D. The chiral and continuum extrapolation of Φ D r 3/2 0 and Φ D /M 3/2 c are hown in Fig. 2. The ytematic uncertainty aociated with the chiral extrapolation ha been tudied uing both a linear or a quadratic fit in. We have taken into account alo the error induced by the uncertaintie on the phyical trange and charm quark mae, and we have included it in the tat+fit error. A for FSE, the comparion of our lattice data at the ame quark ma and different lattice volume how that it effect i well within the tatitical uncertainty. All the reult, including the one obtained uing the Z P value from the M1 or M2 method [9], have been combined to get our final reult for f D, namely f D = 242.1(7.6) tat+ f it (1.4) Chiral (2.9) Dic. (0.3) ZP MeV = 242.1(8.3)MeV. (3.2) Our reult obtained at N f = can be compared with the FLAG average f D = 248(6)MeV at N f = 2 and f D = 248.6(2.7)MeV at N f = [12]. The ratio f D / f D can be calculated by analyzing the lattice data of the ratio Φ D /Φ D. However, analyzing intead the double ratio ( f D / f D )/( / f π ) increae the preciion on the determination of f D / f D, becaue of the very mild dependence of the double ratio on [14]. To fit the double ratio we combined the ChPT prediction for f π and with the HMChPT prediction for Φ D /Φ D, obtaining the following formula [ ( f D / f D 9 = P P 2 + 4ĝ2 1 ] K FSE f )ξ l logξ π l / f π 2 K FSE, (3.3) where for the coupling contant ĝ we have conidered the value ĝ = 0.61(7) [15], which among the preently available determination of ĝ maximize the impact of the chiral log in Eq. (3.3). Notice that in Eq. (3.3) any dependence of the double ratio on the lattice pacing i neglected, ince uch a dependence i not viible in the lattice data (ee Fig. 3). We checked that by performing the ame fit in without the data correponding to the coarer lattice pacing (which correpond roughly to keep half of the data) the ame reult for the double ratio i obtained. 5

6 Notice alo in Eq. (3.3) the preence of the FSE correction for both f π and. We decribed them uing the CWW [17] and CDH [10] formulae, repectively, which reproduce the FSE oberved in the data at the ame light-quark ma and lattice pacing, but different lattice volume. An alternative fit with no chiral log wa performed in order to evaluate the ytematic error aociated to chiral extrapolation, namely f D / f D / f π = P 1 (1 + P 2 ) KFSE f π K FSE. (3.4) The chiral extrapolation for the double ratio ( f D / f D )/( / f π ), uing both the ChPT fit (3.3) and the linear one (3.4), i hown in Fig. 3. (f D /f D )/( /f π ) Cont. Limit Linear Fit Cont. Limit Chiral Fit 0.8 Figure 3: Chiral and continuum extrapolation of the double ratio ( f D / f D )/( / f π ) uing both the prediction (3.3) from ChPT and the polynomial expanion (3.4) in the light quark ma. It can be een that the two anätze choen for the chiral extrapolation provide compatible reult for all pion mae within the tatitical uncertaintie. Our data have been extrapolated to m ud and therefore our reult for f D / f D correpond to the QCD iopin ymmetric limit. Iopin breaking effect are etimated to be below the percent level and therefore mall compared to other uncertaintie. The relevant ource of ytematic error for the double ratio ( f D / f D )/( / f π ) i the chiral extrapolation, while for f D / f D alo the uncertainty on / f π ha to be conidered. On the contrary, the error on the phyical trange and charm quark mae a well a the dicretization error and the uncertainty on the RC Z P contribute negligibly. Our final reult for ( f D / f D )/( / f π ) and f D / f D are f D / f D / f π = 1.005(14) tat+ f it (6) Chiral = 1.005(15), (3.5) f D / f D = 1.199(17) tat+ f it (7) Chiral (16) fk / f π = 1.199(25). (3.6) The latter one can be compared with the FLAG average f D / f D = 1.17(5) at N f = 2 and f D / f D = 1.187(12) at N f = 2+1. Notice the remarkable reult for the double ratio given in Eq. (3.5), which 6

7 mean that SU(3) breaking effect in the ratio of PS meon decay contant are the ame in the light and charm ector within a percent accuracy. Finally we combined our reult for f D and f D / f D to obtain for f D the value f D = 201.9(8.0)MeV. (3.7) The FLAG average [12] are f D = 212(8)MeV at N f = 2 and f D = 209.2(3.3)MeV at N f = Acknowledgement We acknowledge the CPU time provided by the PRACE Reearch Infratructure under the project PRA027 QCD Simulation for Flavor Phyic in the Standard Model and Beyond" on the JUGENE BG/P ytem at the Jülich SuperComputing Center (Germany), and by the agreement between INFN and CINECA under the pecific initiative INFN-RM123. Reference [1] Y. Iwaaki, Nucl. Phy. B 258 (1985) 141. [2] R. Frezzotti and G.C. Roi, Nucl. Phy. Proc. Suppl. 128 (2004) 193 [hep-lat/ ]. [3] R. Frezzotti and G.C. Roi, JHEP 0408 (2004) 007 [hep-lat/ ]. [4] K. Oterwalder and E. Seiler, Annal Phy. 110 (1978) 440. [5] N. Carraco et al. [ETM Coll.], PoS LATTICE 2013 (2013) 312 [6] R. Sommer, Nucl. Phy. B 411 (1994) 839 [hep-lat/ ]. [7] R. Baron et al., JHEP 1006 (2010) 111 [arxiv: [hep-lat]]. [8] R. Baron et al. [ETM Coll.], PoS LATTICE 2010 (2010) 123 [arxiv: [hep-lat]]. [9] B. Bloier et al. [ETM Coll.], PoS LATTICE 2011 (2011) 233 [arxiv: [hep-lat]]. [10] G. Colangelo, S. Durr and C. Haefeli, Nucl. Phy. B 721 (2005) 136 [arxiv:hep-lat/ ]. [11] G. M. de Divitii et al., Phy. Rev. D 87 (2013) [arxiv: [hep-lat]]. [12] S. Aoki et al., arxiv: [hep-lat]. For the publihed verion of the FLAG-1 review ee: G. Colangelo et al., Eur. Phy. J. C 71 (2011) 1695 [arxiv: [hep-lat]]. [13] F. Farchioni et al., PoS LATTICE 2010 (2010) 128 [arxiv: [hep-lat]]. [14] D. Becirevic et al., Phy. Lett. B 563 (2003) 150 [hep-ph/ ]. [15] K. Nakamura et al. [Particle Data Group Coll.], J. Phy. G 37 (2010) [16] O. Bar, Phy. Rev. D 82 (2010) [arxiv: [hep-lat]]. [17] G. Colangelo, U. Wenger and J.M.S. Wu, Phy. Rev. D 82 (2010) [arxiv: [hep-lat]]. 7