% IMPORTANT: The following is UTF-8 encoded.  This means that in the presence
% of non-ASCII characters, it will not work with BibTeX 0.99 or older.
% Instead, you should use an up-to-date BibTeX implementation like “bibtex8” or
% “biber”.

@ARTICLE{Aoki:15439,
      author       = {Aoki, Y. and Arthur, R. and Blum, T. and Boyle, P. and
                      Brömmel, D. and Christ, N. and Dawson, C. and Flynn, J. and
                      Izubuchi, T. and Jin, X. and Jung, C. and Kelly, C. and Li,
                      M. and Lichtl, A. and Lightman, M. and Lin, M. and
                      Mawhinney, R. and Maynard, C. and Ohta, S. and Pendleton, B.
                      and Sachrajda, C. and Scholz, E. and Soni, A. and Wennekers,
                      J. and Zanotti, J. and Zhou, R.},
      title        = {{C}ontinuum {L}imit {P}hysics from 2+1 {F}lavor {D}omain
                      {W}all {QCD}},
      journal      = {Physical review / D},
      volume       = {83},
      number       = {7},
      issn         = {1550-7998},
      address      = {[S.l.]},
      publisher    = {Soc.},
      reportid     = {PreJuSER-15439},
      pages        = {074508},
      year         = {2011},
      note         = {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).},
      abstract     = {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.},
      keywords     = {J (WoSType)},
      cin          = {JSC},
      ddc          = {530},
      cid          = {I:(DE-Juel1)JSC-20090406},
      pnm          = {Scientific Computing (FUEK411) / 411 - Computational
                      Science and Mathematical Methods (POF2-411)},
      pid          = {G:(DE-Juel1)FUEK411 / G:(DE-HGF)POF2-411},
      shelfmark    = {Astronomy $\&$ Astrophysics / Physics, Particles $\&$
                      Fields},
      typ          = {PUB:(DE-HGF)16},
      UT           = {WOS:000290110100003},
      doi          = {10.1103/PhysRevD.83.074508},
      url          = {https://juser.fz-juelich.de/record/15439},
}