%0 Journal Article
%A Aoki, Y.
%A Arthur, R.
%A Blum, T.
%A Boyle, P.
%A Brömmel, D.
%A Christ, N.
%A Dawson, C.
%A Flynn, J.
%A Izubuchi, T.
%A Jin, X.
%A Jung, C.
%A Kelly, C.
%A Li, M.
%A Lichtl, A.
%A Lightman, M.
%A Lin, M.
%A Mawhinney, R.
%A Maynard, C.
%A Ohta, S.
%A Pendleton, B.
%A Sachrajda, C.
%A Scholz, E.
%A Soni, A.
%A Wennekers, J.
%A Zanotti, J.
%A Zhou, R.
%T Continuum Limit Physics from 2+1 Flavor Domain Wall QCD
%J Physical review / D
%V 83
%N 7
%@ 1550-7998
%C [S.l.]
%I Soc.
%M PreJuSER-15439
%P 074508
%D 2011
%Z The calculations reported here were performed on the QCDOC computers [80-82] at Columbia University, Edinburgh University, and at the Brookhaven National Laboratory (BNL). At BNL, the QCDOC computers of the RIKEN-BNL Research Center and the USQCD Collaboration were used. Most important were the computer resources of the Argonne Leadership Class Facility (ALCF) provided under the Incite Program of the U. S. DOE. The very large-scale capability of the ALCF was critical for carrying out the challenging calculations reported here. We also thank the University of Southampton for access to the Iridis computer system used in the calculations of the nonperturbative renormalization factors (with support from UK STFC Grant No. ST/H008888/1). The software used includes: the CPS QCD codes http://qcdoc.phys.columbia.edu/chulwoo/index.html, supported in part by the U.S. DOE SciDAC program; the BAGEL http://www.ph.ed.ac.uk/ paboyle/bagel/Bagel.html assembler kernel generator for many of the high-performance optimized kernels [25]; and the UKHADRON codes. Y. A. is partially supported by JSPS Kakenhi Grant No. 21540289. R. A., P. A. B., B. J. P., and J. M. Z. were partially supported by UK STFC Grant No. ST/G000522/1. T. B. and R. Z. were supported by U.S. DOE Grant No. DE-FG02-92ER40716. D. B., J. M. F., and C. T. S. were partially supported by UK STFC Grant No. ST/G000557/1 and by EU Contract No. MRTN-CT-2006-035482 (Flavianet). N. H. C., M. L., and R. D. M. were supported by U.S. DOE Grant No. DE-FG02-92ER40699. C. J., T. I., and A. S. are partially supported by the U.S. DOE under Contract No. DE-AC02-98CH10886. E. E. S. is partly supported by DFG SFB/TR 55 and by the Research Executive Agency of the European Union under Grant No. PITN-GA-2009-238353 (ITN STRONGnet).
%X We present physical results obtained from simulations using 2 + 1 flavors of domain wall quarks and the Iwasaki gauge action at two values of the lattice spacing a, [a(-1) = 1.73(3) GeV and a(-1) = 2.28(3) GeV]. On the coarser lattice, with 24(3) x 64 x 16 points (where the 16 corresponds to L-s, the extent of the 5th dimension inherent in the domain wall fermion formulation of QCD), the analysis of C. Allton et al. (RBC-UKQCD Collaboration), Phys. Rev. D 78 is extended to approximately twice the number of configurations. The ensembles on the finer 32(3) x 64 x 16 lattice are new. We explain in detail how we use lattice data obtained at several values of the lattice spacing and for a range of quark masses in combined continuum-chiral fits in order to obtain results in the continuum limit and at physical quark masses. We implement this procedure for our data at two lattice spacings and with unitary pion masses in the approximate range 290-420 MeV (225-420 MeV for partially quenched pions). We use the masses of the pi and K mesons and the Omega baryon to determine the physical quark masses and the values of the lattice spacing. While our data in the mass ranges above are consistent with the predictions of next-to-leading order SU(2) chiral perturbation theory, they are also consistent with a simple analytic ansatz leading to an inherent uncertainty in how best to perform the chiral extrapolation that we are reluctant to reduce with model-dependent assumptions about higher order corrections. In some cases, particularly for f(pi), the pion leptonic decay constant, the uncertainty in the chiral extrapolation dominates the systematic error. Our main results include f(pi) = 124(2)(stat)(5)(syst) MeV, f(K)/f(pi) = 1.204(7)(25) where f(K) is the kaon decay constant, m(s)((MS) over bar) (2 GeV) = (96.2 +/- 2.7) MeV and m(s)((MS) over bar) (2 GeV) (3.59 +/- 0.21) MeV (m(s)/m(ud) = 26.8 +/- 1.4) where m(s) and m(ud) are the mass of the strange quark and the average of the up and down quark masses, respectively, [Sigma((MS) over bar) (2 GeV)(1/3) = 256(6) MeV, where Sigma is the chiral condensate, the Sommer scale r(0) = 0.487(9) fm and r(1) = 0.333(9) fm.
%K J (WoSType)
%F PUB:(DE-HGF)16
%9 Journal Article
%U <Go to ISI:>//WOS:000290110100003
%R 10.1103/PhysRevD.83.074508
%U https://juser.fz-juelich.de/record/15439