We report on our calculation of K \to \pi vector form factor by numerical
simulations of two-flavor QCD on a 16^3x32x12 lattice at a \simeq 0.12 fm using
domain-wall quarks and DBW2 glue. Our preliminary result at a single sea quark
mass correponding to m_PS/m_V \simeq 0.53 shows a good agreement with previous
estimate in quenched QCD and that from a phenomenological model.Comment: 6 pages, 5 figures, poster presented at Lattice2005 (Weak matrix
  elements); v2: a reference adde

Dawson, C.

Izubuchi, T.

Kaneko, T.

Sasaki, S.

Soni, A.

English

arXiv

We report on our calculation of K {yields} {pi} vector form factor by numerical simulations of two-flavor QCD on a 16{sup 3} x 32 x 12 lattice at a {approx_equal} 0.12 fm using domain-wall quarks and DBW2 glue. Our preliminary result at a single sea quark mass corresponding to m{sub PS}/m{sub V} {approx_equal} 0.53 shows a good agreement with previous estimate in quenched QCD and that from a phenomenological model

Dawson, T.

UNT (University of North Texas) Digital Library

N AT 1 0 N A t LAB 0 R A T 0  RU BNL-75638-2006-CP K(13) Form Factor With Two Flavors of Dynamical Domain Wall Quarks A. Sonil, C. Dawson2, T. I z ~ b u c h i ~ ~  3, T. Kaneko4? ’,S. Sasaki2y Physics Department, Brookhaven National Laboratory, Upton, NY 11973 USA RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, NY 1 1973 ’Institute of Theoretical Physics, Kanazawa University, Ishikawa 920- 1 192, Japan ‘High Energy Accelerator Research Organization (KEK), Ibaraki 305-080 1, Japan 5The Graduate University for Advanced Studies, Ibaraki 305-080 1, Japan 6Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 1 2 Preserited at the 23rd International Symposium on Lattice Field Theoiy Dublin, Ireland 25-30 July 2005 March 2006 Physics Department Brook h aven National La bo rat o ry P.O. Box 5000 Upton, NY 1 1973-5000 www.bnl.gov Notice: This manuscript has been authored by employees of Brookhaven Science Associates, LLC under Contract No. DE-AC02-98CH10886 with the U.S. Department of Energy. The publisher by accepting the manuscript for publication acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. This preprint is intended for publication in a journal or proceedings. Since changes may be made before publication, it may not be cited or reproduced without the author’s permission. DISCLAIMER This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, nor any of their contractors, subcontractors, or their employees,. makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or any third party's use or the results of such use of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof or its contractors or subcontractors. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. P KEK-CP- 174 BNL-HET-05/23 KL3 form factor with two-flavors of dynamical domain-wall quarks RBC Collaboration: C. Dawson RIKEN-BNL Research Centel; Brookhaven National Laboratoiy, Uptoit, NY 11973, USA T. lzubuchi Institute of Theoretical Physics, Kanazawa UniversiQ, Ishikawa 920-1 192, Japan RIKEN-BNL Research Centei; Brooklaaven National Laboratoiy Upton, NY 11973, USA T. Kaneko* High Energy Accelerator Research Organization (KEK), Ibaraki 305-0801, Japan The Graduate Universit), for Advanced Studies, Ibaraki 305-0801, Japan S. Sasaki Departtileiit of Physics, University of Tokyo, Tokyo 113-0033, Japan RIKEN-BNL Research Center; Brookhaven National Laboratory9 Upton, NY 11973, USA A. Soni Playsics Department, Brooklzaven National Laboratory, Upton, NY II973, USA We report on our calculation of K --+ n vector form factor by numerical simulations of two-flavor QCD on a 163 x 32 x 12 lattice at a N 0.12 fm using domain-wall quarks and DBW2 glue. Our preliminary result at a single sea quark mass correponding to rnpslinv N 0.53 shows a good agreement with previous estimate in quenched QCD and that from a phenomenological model. XXIIIrd Interriatiorial Syiiposiuin on Lattice Field Theory 25-30 JLL~JJ  2005 Trinity College. Dubliri, Ireland Speaker. @ Copyright owned hy the author(s) under the terms irf the Creative Cammons Attribution-NonCominercial-ShareAlikeLicence. http://pos.sissa.it/ Kl3 foriiz factor with twoJEavois of dynamical doiizairi-wall quarks T. Kaneko 1. Introduction Currently. K i nlvl (Kl3) decay provides the most precise determination of the Cabibbo- Kobayashi-Maskawa (CKM) matrix element /VJl,l through I?= ~~,,~2~f+(0)~2, where r is the decay rate and f+(q2) is the form factor defined from the K + K matrix element of the weak vector current (44 I v, I K ( P ) )  = (Py + PLj f+ (2)  + ( P y  - P L )  f- (2) 7 q2 = ( P  - P Y  ( 1.1 1 While the PDG values for the CKM matrix elements [l] show a 2 0  deviation from the CKM unit arity IV,rr,12 + (V,rs12 + IV,,/,1' = 1 - 6, 6 = 0.0033( 15), (1.2) recent experiments for r [2, 3, 4, 51 prefer a slightly larger value of IV,,,I, which is consistent with the unitarity (6 =O).  However, in order to make a definite conclusion on this issue, f+ (0) has to be determined theoretically with an accuracy of about 1 %. In chiral perturbation theory (ChPT), f+(O) is expanded in powers of meson masses M K .  M,. and Mq f+(o) = l + f i + f 4 + ' * ' >  ( h = O ( M g , ~ , q ) ) *  (1.3) Thanks to the Ademollo-Gatto theorem [6], which states that SU(3) breaking effects in f+(O) start at O((m,  - muc/)2j, poorly-known low-energy constants of the chiral Lagrangian do not enter the ChPT formula for f2 [7], and hence it is precisely determined as -0.023. However, this is not the case for the higher order coirections f i l  (~124)  [8]. Therefore, a phenomenological estimate of fd based on the quark model [9] has been used in previous analyses of IV,,,I. In order to estimate f+ (0) without relying on any phenomenological model, lattice calculations have been carried out first in the quenched approximation [lo], and later in unquenched QCD [ 11, 121. In these calculations, however, conventional Wilson- or Kogut-Susskind-type quark actions are employed, and hence chiral properties of f+(O) may be significantly affected by the explicit breaking of chiral or flavor symmetry. Since the ChPT formula for f2 plays a crucial role in the chiral extrapolation of lattice data, it is advantageous to use a quark action, which posses chiral symmetry even at finite lattice spacing. In this work, we calculate f+ (0) in two-flavor dynamical QCD using the DBW2 gauge [13] and the domain-wall quark actions [14], with which the hadron spectrum and the kaon B parameter show good chiral properties [ 151. 2. Simulation method Our calculations are carried out on a 163 x 32 lattice with statistics of 4750 HMC trajectories. While we simulate a single value for the lattice spacing =1.69(5) GeV, we expect that the scal- ing violation in f+ is not large, since the employed lattice action is (automatically) O(n)-improved. The size of the fifth dimension is set to Ls = 12, which leads to the residual mass of a few MeV. We refer to Ref. [ 151 for details on the gauge ensembles used in this study. At the moment. we complete our calculation at a single sea quark mass nz,,,/= 0.02, which is roughly half the physical strange quark mass. Three heavier masses in, = 0.03, 0.04 and 0.05 are employed for the valence strange quarks. 337 I 2  Kl3 f o m  factol- with two flavors of dyriamical domain-wall quarks T. Kaneko 3. Extraction of form factor The so-called scalar form factor at qiax = ( M K  - Mn)2 can be.extracted from'the double ratio proposed in Ref. [16], If0 (Cl iax  1 I 3 (3.2) c4K" (t ) t'; 0,O) C,.K (t , t' ; 0,O) (MY + WJ2 R(t,t ') = C,""(t, t';O, 0) CY(t, P ; O ,  0) t:(t'--t)-i- 4MKMn CLQ(t,t';p,p') = ( o e ( ~ ' , t ' ) V ~ ( x , t ) O ~ ( O , O ) )  e-ip'(x'-x)e-ipx, (P,Q = x or K ) ,  (3.3) where  on(^) represents the interpolating operator for pion (kaon). As shown in Fig. 1, fo(qkax) is determined with an accuracy of 5 0.1 %, since various uncertainties of the three-point function. such as the statistical fluctuation, are canceled at least partially in the ratio. x. X' To study the q2 dependence of the form factor, we calculate from a ratio (3 .5)  EK(P) +EdP') F ( p  /) ) P  > ___$ C y y t ,  t'; p, p') cyt; 0) C y  - t ;  0) cy (t , t' ; 0,O) CK (t ; p) Cn (t' - t ; p') t:(l'-c)--- M K  + Mn where C"(K)(t;p) represents the pion (kaon) propagator with the spatial momentum p. Figure 1 shows a plot of F ( y , p ' ) .  The accuracy of F ( p , p ' )  is typically 5- 10% for the spatial momenta In order to convert F ( p , p ' )  to f 0 ( q 2 ) ,  5 (q2)  is estimated from the double ratio proposed in IPI, IP'I 5 Jz. Ref. [ 1 01 (3.6) 337 I3 K L ~  form factor with two flavors o f  dynainical domain-wall quarks T. Kaneko We observe that <(q2) has a mild dependence on the valence quark mass and its magnitude is typically e -0.01 with 50- 100% error. 4. Interpolation to q2 =O Here, we test two methods to determine f+(O) (=fo(O)). 0 method-1: As in Ref.[lO], we first calculate fo (q2)  at each simulated q' from fo(qiax). F ( p , i )  and <(q2),  and then interpolate f o (q2 )  to q2=0. e method-2: As in Ref.[l 11, F ( p , p ' )  (5 (q2)) is interpolated (extrapolated) to q2 =0, and then fo (0)  is calculated at q2=0. The interpolation of F(p ,p ' )  is carried out using data with fixed IpI (or lp'l) so that we can unambiguously identify lp'l ([pi) corresponding to q2 =0, which is needed to convert F(y,p')l+, to fo(0) .  We repeat this analysis for two'data sets with p = 0 and p' = 0, and take the average of results for fo (0). We test quadratic and polar fits @(O) 8 ( q 2 )  = ~ 1 - c1 q 2  (4.1 ) (4.2) for the interpolation of fo in method-1 and F for method-2, while only the quadratic fit is tested for 5. We observe that four possible ways (quadratic or polar fit for method-1 or 2) give mutually consistent results with an accuracy of 5 1 %. This is because we have very accurate data of fo (&,) near q2 = 0, and hence the uncertainty due to the choice of the interpolation method is not large. However, we observe that method-2 leads to a slightly smaller error of fo (0) than method-1. since the kinematical factor (EK(~)-E,(~'))/(EK(~)+E~(~')) in Eq. (3.4) is not large at q2 =0, and hence the uncertainty of 5 (0) has small influence to fo(0) in method-2. Therefore, we employ f+(O) obtained from method-2 with the quadratic fit in the following. -0.25 -0.20 -0.15 -0.10 -0.05 0.W 0.03 0.01 ' 1 ' I " ' I ' I " ' I Y- Y- Figure 2: Interpolation of fo(q2)  (left figure) and F ( p ,  p') with p'= O (right figure) to q2 = 0. Both figures show data at in,  = 0.04. 337 I 4  K13 forin factor with twojc~vors  of dynumical dornain-wall quarks T. Kaneko 5. Chiral extrapolation From the Ademollo-Gatto theorem, the higher order correction is proportional to ( 1 7 7 ~ - n z ~ ~ ~ , ) ~ .  Therefore, as in Ref.[10], it is convenient to consider a ratio for the chiral extrapolation of f+(O). In this analysis, f2 at the simulated quark masses is cal by using the ChPT formula in unquenched QCD [7]’. (5.1 j (5.2) dated Our results for R A ~  are plotted in Fig. 3. Since the single sea quark mass is simulated so far, we extrapolate R A ~  only in terms of the valence quark mass by naively assuming that its sea quark mass dependence is not large below the simulated sea quark mass. From a simple linear fit shown in Fig. 3, we obtain f+(O) = 0.955( 12), which is consistent with both of the previous phenomenological estimate 0.961(8) [9] and the quenched result 0.960(9) [ 101. We note that recent unquenched calculations in Refs.[ 11, 121 also obtained similar values for f+(O). 6. Conclusions We have calculated f +  (0) in two-flavor QCD using the domain-wall quarks. Our preliminary result is consistent with the previous phenomenological estimate with an accuracy of 1% level. this work lattice, N,=2, domain-wall Q Leutwyler-Roos: quark model Becirevic et al.: lattice, N,=O. clover 0.0 Figure 3: Quark mass dependence of R A ~ .  Open circle and square represent values corresponding to the phenomenological estimate [9] and the quenched result [ 101. ‘Very recently, f2 is calculated in partially quenched ChPT (PQChPT) [17]. We confirm that f+ (0) obtained from the chiral extrapolation with the PQChPT formula for f2 is consistent with the result presented in the text within the statistical errors. 337 I 5  K L ~  .foriii.factor with two flavors o f  dytiamical domain-wall quarks T. Kaneko Our estimate combined with recent experimental results of r leads to lvcsl which is consistent with the CKM unitarity. For instance, by using I? from the E865 Collaboration 1121, we obtain lv,J I = 0.229(4) which leads to IV,</I2 + IVr,J2 + IV,/,I2 = 1 - 6 ,  6 = -0.001(2). (6.1 ) We note that, however, our result has an additional uncertainty arising from the fact that we have not taken the limit of the physical sea quark mass. To remove this uncertainty, our simulations at two different sea quark masses are in progress. The nice consistency between our and previous quenched estimates of f+(O) may suggest that the systematic error due to the quenched approximation for strange quarks is not large.'This point, however, has to be confirmed by extending our calculation to three-flavor QCD. We thank RIKEN, BNL and the US .  Department of Energy for providing the facilities essen- tial for this work. The work of TK is supported in part by the Grant-in-Aid of the Japanese Ministry of Education (Nos.17740171). The work of AS was supported in part by US DOE Contract No. DE-AC02-98CH 10886. References [ 1 J S. Eidelman et al., Phys. Lett. B 592, 1 (2004). [2] A. Slier et al. (E865 Collaboration), Phys. Rev. Lett. 91, 261802 (2003). [3] T. Alexopoulos et al. (KTeV Collaboration), Phys. Rev. Lett. 93, 181802 (2004). [4] A. Lai et a/. (NA 48 Collaboration), Phys. Lett. B 602.41 (2004). [5] P. Franzini, econf C0406271 (2004); G. Lanfranchi (KLOE Collaboration), Int. J. Mod. Phys A 20, 493 (2005). [6] R.E. Behrends and A. Sirlin, Phys. Rev. Lett. 4, 186 (1960); M. Ademollo and R. Gatto, Phys. Rev. [7] J. Gasser and H. Leutwyler, Nucl. Phys. B 250,517 (1985). [8] J. Bijnens, G. Colangelo and G. Ecker, Phys. Lett. B 441,437 (1998); P. Post and K. Schilcher, Eur. Lett. 13, 264 (1964). Phys. J. C 25,427 (2002); J. Bijnens and P. Talavera, Nucl. Phys. B 669,341 (2003). [9] H. Leutwyler and M.Roos, Z. Phys. C 25,91 (1984). [lo] D. BeCireviC et al., Nucl. Phys. B 705,339 (2005); Eur. Phys. J. A24S1,69 (2005). [l 11 N. Tsutsui et al. (JLQCD Collaboration), in these proceedings. [ 121 M. Okamoto (Fermilab-MILC-HPQCD Collaborations), hep-lat/04 12044. [13] T. Takaishi, Phys. Rev. D 54, 1050 (1996); Ph. de Forcrand et al. (QCD-TARO Collaboration), Nucl. Phys. B 577,263 (2000). [14] D.B. Kaplan, Phys. Lett. B 288,342 (1992); Y. Shamir, Nucl. Phys. B 406.90 (1993). [ 151 Y. Aoki et al. (RBC Collaboration), hep-lat/O411006. [16] S. Hashimoto et a/.,  Phys. Rev. D 61,014502 (2000) [ 171 D. Bedirevie, G. Martinelli and G. Villadoro, hep-lat/0508013. 337 I 6  

K_l3 form factor with two-flavors of dynamical domain-wall quarks

Abstract

Similar works

Full text

Available Versions

UNT (University of North Texas) Digital Library