arxiv: v1 [hep-lat] 21 Nov 2011

Transcription

1 Deta I=3/2 K to pi-pi decays with neary physica kinematics arxiv: v1 [hep-at] 21 Nov 2011 University of Southampton, Schoo of Physics and Astronomy, Highfied, Southampton, SO17 1BJ, United Kingdom E-mai: ejg4g08@soton.ac.uk Matthew Lightman Department of Physics, Coumbia University, New York, NY USA Department of Physics, Washington University, St. Louis, MO USA E-mai: mightman@physics.wust.edu The I = 3/2 K ππ decay ampitude is cacuated on RBC/UKQCD , L s = 32 dynamica attices with favors of domain wa fermions using the Disocation Suppressing Determinant Ratio and Iwasaki gauge action. The cacuation is performed cose to the physica pion mass (m π = 142.9(1.1) MeV) and with a singe attice spacing (a 1 = 1.375(9) GeV.) We find Re(A 2 ) = (1.436 ± stat ± syst ) 10 8 GeV and Im(A 2 ) = ( 6.29 ± 0.46 stat ± 1.20 syst ) GeV. These resuts are combined with the experimenta resut for ε /ε to predict Im(A 0 ) = 5.32(64) stat (71) syst GeV within the Standard Mode. We aso perform a reweighting anaysis to invesitgate the effects of partia quenching in the ight-quark sector of our cacuation. Foowing reweighting we find Re(A 2 ) = (1.52 ± 0.14 stat ) 10 8 GeV and Im(A 2 ) = ( 6.47 ± 0.55 stat ) GeV, which are consistent with our main resuts. XXIX Internationa Symposium on Lattice Fied Theory Juy Squaw Vaey, Lake Tahoe, Caifornia Speaker. c Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercia-ShareAike Licence.

2 1. Introduction The cacuation of K ππ decay ampitudes is motivated by a desire to understand the I = 1/2 rue and CP vioation in kaon decays. Such a cacuation is a non-perturbative probem requiring attice techniques to make progress. Previous attice cacuations have reied on the quenched approximation and uncertain chira extrapoations [1, 2, 3, 4]. In this tak we present the resuts from the first reaistic attice cacuation of a K ππ decay ampitude, where we simuate the two-body decay directy on the attice at neary-physica kinematics. We proceed by evauating matrix eements of the S = 1 effective Hamitonian H eff = VusV ud C i Q i (1.1) i where C i are Wison coefficients and Q i are four-quark operators. In this tak we consider ony the I = 3/2 transition, in which case ony three operators contribute in equation (1.1). We find it convenient to evauate unphysica K + π + π + matrix eements of the foowing operators: Q (27,1) = (s i d i ) L (u j d j ) L, Q (8,8) = (s i d i ) L (u j d j ) R and Q (8,8)mx = (s i d j ) L (u j d i ) R, (1.2) where the operators are abeed according to their transformation under SU(3) L SU(3) R, and the abes i and j on the quark fieds abe coour. These are reated to the physica K + π + π 0 matrix eements via the Wigner-Eckart theorem. 2. Detais of the Simuation The anaysis is performed on a singe ensembe of favour domain wa fermions (DWF) with Disocation Suppressing Determinant Ratio (DSDR)+Iwasaki gauge action at β = The attice size is and the extent of the fifth dimension is L s = 32. The inverse attice spacing is a 1 = 1.375(9) GeV. The ensembe is genereated with sea-quark masses am = and am h = 0.045, corresponding to a unitary pion mass of approximatey 170 MeV. We find the residua mass to be am res = (1). The correation functions are cacuated with vaence quark masses of m = and m s = 0.049, corresponding to a vaence pion mass of m π = 142.9(1.1) MeV and kaon mass of m K = 511.3(3.9) MeV. A tota of 63 gauge configurations, each separated by 8 moecuar dynamics time units, are incuded in the anaysis. Quark propagators with periodic and antiperiodic boundary conditions in the time direction were computed on each configuration with a source at t = 0. They were then combined so as to effectivey doube the time extent of the attice. Meson correation functions formed using the average of the propagators with periodic and antiperiodic boundary conditions can be interpreted as containing forward propagating mesons originating at time t = 0, whereas those cacuated with the antisymmetric combination can be interpreted as containing backward propagating mesons originating from a source at t = 64. Strange-quark propagators, aso with Periodic + Antiperiodic combinations, were generated with sources at t K = 20, 24, 28, 32, 36, 40 and 44 in order to cacuate K ππ correation functions with kaon sources at these times, whie the two-pion sources remained at either t = 0 or t = 64. Thus we coud achieve time separations between the kaon and two pions of 20, 24, 28 and 32 attice time units in two different ways which increased the statistics. 2

3 These separations were chosen so that the signas from the kaon and two pions did not decay into noise before reaching the four-quark operator Q i. For physica K ππ decays in the CM frame, the fina state pions have equa and opposite non-zero momentum. We achieve this for the π + π + fina state by giving momentum to the d-quark. The u- and s-quark propagators are generated with Couomb gauge fixed wa sources and periodic spatia boundary conditions. Simiary, the d-quark propagators used in the zero-momentum pion and two-pion correation functions are computed with Couomb gauge fixed wa sources and periodic spatia boundary conditions. However, the d-quark propagators used in the two-pion and K ππ correation functions with non-zero momentum pions are generated with antiperiodic spatia boundary conditions and cosine sources. By imposing antiperiodic spatia boundary conditions the aowed quark momenta are p n = (π + 2πn)/L, where L is the spatia extent of the attice, corresponding to a ground-state momentum of π/l in the direction in which antiperiodic boundary conditions have been used. In practice we impose antiperiodic boundary conditions in two spatia directions, which aows us to simuate pions with ground state momentum ± 2π/L. This decision is motivated by the expectation that pions with p = 2π/L wi correspond to a two-pion fina state with energy cose to m K. The use of cosine sources in K + π + π + decays is described in [5]. 3. Anaysis We extract the K ππ matrix eement M by fitting a constant to the eft hand side of (3.1) CKππ i (t) C K (t K t)c ππ (t) = M i. (3.1) Z K Z ππ CKππ i is the K ππ correator with a kaon source at t K, i abes the four-quark operator Q i which is inserted at time t, and Z K and Z ππ are cacuated from the kaon and two-pion correators respectivey, whose sources are at t = 0. The eft hand side of equation (3.1) is potted in Figure 1 for each operator. The figure demonstrates that sufficienty far from the kaon and two-pion sources we are justified in fitting to a constant. The fit resuts for M i /(Z K Z ππ ) are indicated on the pot. We aso use a quotient method to extract the two-pion energy, as we find this improves the statistica precision of the fits. We fit the quotient of correators C ππ /(C π ) 2 Ae E t to extract E = E ππ 2E π. We then get the two-pion energy by cacuating E + 2E π, where in the case of p = 0, E π = m π whie for p = 2π/L, E π is found from a 2 parameter fit to the pion correation function which aso has p = 2π/L. The numerica resuts for a the meson masses and energies which we extract from the correation functions are given in Tabe 1. From Tabe 1 we see that the kaon mass is not exacty equa to the two-pion energy, so our cacuation is not quite on-she. This is taken into account when estimating the systematic error on the fina resuts. The finite voume matrix eements are reated to the infinite voume ampitudes A i (where i abes the 4-quark operator) using the Leouch-Lüscher factor [6, 7]. In particuar we have [ ] A i = 1 1 φ + δ L 3/2 m K E ππ M i (3.2) 2 πq π q π q π where the quantity in square brackets contains the effects of the Leouch-Lüscher factor beyond the free fied normaization, δ is the s-wave phase shift, q π is a dimensioness quantity reated to 3

4 units m π m K E π,2 E ππ,0 E ππ,2 m K E ππ,2 attice (32) (91) (14) (63) (35) (35) MeV 142.9(1.1) 511.3(3.9) 238.8(2.4) 288.0(2.2) 492.6(5.5) 18.7(4.8) Tabe 1: Resuts for meson masses and energies. The subscripts 0, 2 denote p = 0, p = 2π/L respectivey. (a) (27, 1) operator (b) (8, 8) operator (c) (8, 8)mix operator Figure 1: K ππ quotient pots for p = 2π/L. The two pion source is at t = 0 whie the kaon source is at t = 24. The dashed ine shows the error on the fit. the individua pion momentum k π via q π = k π L/2π and φ is a kinematic function defined in [6]. The pion momentum k π is cacuated using the dispersion reation E ππ = 2 k 2 π + m 2 π, and differs from 2π/L due to interactions between the two pions. Once q π is known, δ can be cacuated using the Lüscher quantisation condition [8], nπ = δ(k π ) + φ(q π ). As discussed in [5], the phase shift derivative is cacuated using the phenomenoogica curve of [9]. This is necessary because we ony have two vaues of the two-pion energy from which to extract the phase shift. The ampitudes A i are reated to the physica decay ampitude A 2 via 3 A 2 = a 3 2 G FV ud Vus C i (µ)z i j (µ)a j, (3.3) i, j where C i are the Wison coefficients and Z i j are the renormaization constants. The Wison coefficients must be evauated at the same scae and scheme as the renormaization constants. The renormaization constants are first evauated in the RI-SMOM(/q,/q) scheme [10]. In order to minimize discretization effects, this procedure takes pace at a reativey ow energy µ 0 = GeV. A non-perturbative step-scaing function is then used to convert these resuts to a scae of 3 GeV, at which point a perturbative matching to the MS-NDR scheme is possibe. The Wison coefficients are known in the NDR scheme at the W-mass scae, and can be perturbativey run to the desired matching point of 3 GeV [11]. 4. Resuts We cacuate A 2 for the four different separations between the kaon source and two-pion source. Our fina resut, presented in equation (4.1), is an error weighted average over these four 4

5 resuts. The first error in equation (4.1) is a statistica error, where the statistica uncertainties in the ampitude A i and attice spacing (4% in tota) are combined in quadrature with the statistica error on the renormaization constants (0.8% for ReA 2 and 6% for Im(A 2 )). Re(A 2 ) = (1.436 ± stat ± syst ) 10 8 GeV, Im(A 2 ) = ( 6.29 ± 0.46 stat ± 1.20 syst ) GeV. (4.1) The second error in equation (4.1) is systematic. The systematic errors in our cacuation, for (Re(A 2 ), Im(A 2 )) respectivey, are from attice artifacts (15%, 15%), uncertainty in the phase shift derivative (0.32% 0.32%), finite voume (6.2% and 6.8%), partia quenching (3.5% and 1.7%), uncertainties in the renormaization procedure (1.7% and 4.7%), unphysica kinematics (3.0% and 0.22%), and perturbative truncation in the evauation of the Wison coefficients (7.1% 8.1%). Combining in quadrature, we find the systematic errors to be 18% for Re(A 2 ) and 19% for Im(A 2 ). Further detais on how these errors are estimanted can be found in [12, 13]. 5. Reweighting We use the technique of reweighting to test the consequences of the partia quenching in the ight-quark sector of our cacuation. The reweighting is performed in 30 increments from the simuated mass m sea = down to a vaue of m sea = , corresponding to the vaence ightquark mass. The resuts are shown in Figure 2. The rightmost point in Figure 2(a) shows the resut for Re(A 2 ) before reweighting, whie the remaining points show the resuts after reweighting to the mass indicated on the x-axis, ending with m sea = for the eftmost point. Simiary Figure 2(b) shows the effects of reweighting on Im(A 2 ). Examining the figures, it can be seen that the errors on Re(A 2 ) and Im(A 2 ) grow, but the centra vaues remain unchanged within the errors. Since reweighting effectivey reduces the number of configurations contributing to the observabes [14] it is natura that the statistica error shoud increase. However, the observation that the centra vaues are unchanged confirms that partia quenching in the ight quark does not introduce a significant source of systematic error. The fina resuts after reweighting are shown in Tabe 2 where they are compared with the resuts before reweighting. m = m = (reweighted) Re(A 2 ) 1.436(63) 10 8 GeV 1.52(14) 10 8 GeV Im(A 2 ) 6.29(46) GeV 6.47(55) GeV Tabe 2: A 2 before and after reweighting. 6. Prediction for Im(A 0 ) Assuming isospin symmetry, the CP-vioating parameter ε /ε can be expressed in terms of Re(A 0 ), Im(A 0 ), Re(A 2 ) and Im(A 2 ) according to equation (6.1): ( ) ε Re = ω [ Im(A2 ) ε 2 ε Re(A 2 ) Im(A ] 0), (6.1) Re(A 0 ) 5

6 (a) Reweighting Re(A 2 ) (b) Reweighting Im(A 2 ) Figure 2: Reweighting A 2 from m sea = to m sea = ( ) where ω = Re(A 2 )/Re(A 0 ). Re ε ε, ε, ω and Re(A 0 ) are known experimentay and presented in Tabe 3. Combining these known factors with our attice resut for Im(A 2 )/Re(A 2 ) we can determine the the unknown quantity Im(A 0 ) within the Standard Mode, finding Im(A 0 ) = 5.32(64) stat (71) syst GeV. (6.2) The error on ImA 0 is obtained by combining the errors on the quantities in Tabe 3 in quadrature. In equation (6.3) beow we compare the reative contribution to Im(A 0 )/Re(A 0 ) from Im(A 2 )/Re(A 2 ) and the term containing the experimentay known contributions: Im(A 0 ) Re(A 0 ) = Im(A 2 ) Re(A 2 ) 2 ε ω ε ε (6.3) 1.60(19) stat (21) syst 10 4 = 4.38(34) stat (95) syst (18) Thus we see that whie the error on the determination of Im(A 0 ) is dominated by the uncertainty in the experimenta vaue of ε /ε, the contribution of Im(A 2 )/Re(A 2 ) to Im(A 0 ) is significant (about 25% in the determination of Im(A 0 )/Re(A 0 )). Re(ε /ε) (1.65 ± 0.26) 10 3 ω (12) ε (2.228 ± 0.011) 10 3 Re(A 0 ) (18) 10 7 GeV Im(A 2 )/Re(A 2 ) (attice) 4.38(34) stat (95) syst 10 5 Tabe 3: Experimenta vaues of the components of equation (6.1) used in the determination of Im(A 0 ), together with the resuts for Im(A 2 )/Re(A 2 ) from these proceedings. 6

7 7. Concusions We have presented preiminary resuts for the I = 3/2 K ππ decay ampitude on 32 3 attices with favours of DWF and the Iwasaki-DSDR gauge action. We find Re(A 2 ) = (1.436 ± stat ± syst ) 10 8 GeV, Im(A 2 ) = ( 6.29 ± 0.46 stat ± 1.20 syst ) GeV. Our resut for Re(A 2 ) is in good agreement with the experimenta resut of 1.479(4) 10 8 GeV obtained from K + decays. In the future we pan to undertake simiar cacuations of Re(A 0 ) and Im(A 0 ) [15], aowing ε /ε to be cacuated from first principes for the first time. In the mean time we make the prediction, based on the Standard Mode, that Im(A 0 ) = 5.32(64) stat (71) sys GeV. We thank a of our coeagues in the RBC and UKQCD coaborations for their contribution to these resuts, for hepfu discussions and the deveopment and support of the QCDOC hardware and software infrastructure which was essentia to this work. In addition we acknowedge Coumbia University, RIKEN, BNL, ANL, and the U.S. DOE for providing the faciities on which this work was performed. This research used resources of the Argonne Leadership Computing Faciity at Argonne Nationa Laboratory, which is supported by the Office of Science of the U.S. DOE under contract DE-AC02-06CH This work was supported in part by U.S. DOE grant number DE- FG02-92ER E.G. is supported by an STFC studentship and grant ST/G000557/1 and by EU contract MRTN-CT (Favianet). References [1] T. Bum et a. [RBC Coaboration], Phys. Rev. D 68 (2003) [arxiv:hep-at/ ]. [2] C.H. Kim and N.H. Christ, Nuc. Phys. Proc. Supp. 119 (2003) 365 [arxiv:hep-at/ ] [3] J. Noaki et a., Phys. Rev. D 68 (2003) [arxiv:hep-at/ ] [4] N.H. Christ and Li, S., PoS(Lattice 2008)272, [arxiv:hep-at/ ] [5] E. Goode and M. Lightman, PoS(Lattice 2010)313, 2010 [arxiv: v1 [hep-at]] [6] Laurent. Leouch and Martin Lüscher, Commun. Math. Phys 219 (2001) 31 [7] C.-J.D. Lin et a., Nuc. Phys. B619 (2001) 467 [arxiv:hep-at/ ] [8] M. Lüscher, Nuc. Phys. B 354, 531 (1991) [9] A. Schenk, Nuc. Phys. B 363 (1991) 97. [10] N. Garron and A. T. Lyte, PoS(Lattice 2011)335, 2011 [11] G. Buchaa, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68 (1996) 1125 [arxiv:hep-ph/ ]. [12] T. Bum et a. [RBC Coaboration and UKQCD Coaboration], manuscript in preparation [13] M. Lightman, Doctora Thesis (2011) [14] Y. Aoki et a. [RBC Coaboration and UKQCD Coaboration], Phys. Rev. D 83 (2011) [arxiv: [hep-at]]. [15] Q. Liu, PoS(Lattice 2011)287, [arxiv: [hep-at]] 7