Continuum limit physics from 2+1 flavor domain wall QCD

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 243x64x16 points (where the 16 corresponds to Ls, 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 323x64x16 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 π and K mesons and the Ω 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π, the pion leptonic decay constant, the uncertainty in the chiral extrapolation dominates the systematic error. Our main results include fπ=124(2)stat(5)syst MeV, fK/fπ=1.204(7)(25) where fK is the kaon decay constant, msMS(2 GeV)=(96.2±2.7) MeV and mudMS(2 GeV)=(3.59±0.21) MeV (ms/mud=26.8±1.4) where ms and mud are the mass of the strange quark and the average of the up and down quark masses, respectively, [ΣMS(2 GeV)]1/3=256(6) MeV, where Σ is the chiral condensate, the Sommer scale r0=0.487(9) fm and r1=0.333(9) fm.