An analysis of the pion mass and pion decay constant is performed using mixed-action Lattice QCD calculations with domain-wall valence quarks on ensembles of rooted, staggered n_f = 2+1 MILC configurations. Calculations were performed at two lattice
spacings of b~0.125 fm and b~0.09 fm, at two strange quark masses, multiple light quark masses, and a number of lattice volumes. The ratios of light quark to strange quark masses are in the range 0.1 <= m_l / m_s <= 0.6, while pion masses are in the range 235 < m_pi < 680 MeV. A two-flavor chiral perturbation theory analysis of the Lattice QCD calculations constrains the Gasser-Leutwyler coefficients bar{l}_3 and bar{l}_4 to be bar{l}_3 = 4.04(40)(+73-55) and bar{l}_4 = 4.30(51)(+84-60). All systematic effects in the calculations are explored, including those from the finite lattice space-time volume, the finite lattice spacing, and the finite fifth dimension in the domain-wall quark action. A consistency is demonstrated between a chiral perturbation theory analysis at fixed lattice spacing combined with a leading order continuum extrapolation, and the mixed-action chiral perturbation theory analysis which explicitly includes the leading order discretization effects. Chiral corrections to the pion decay constant are found to give f_pi / f = 1.062(26)(+42-40) where f is the decay constant in the chiral limit. The most recent scale setting by the MILC Collaboration yields a postdiction of f_pi = 128.2(3.6)(+4.4-6.0)(+1.2-3.3) MeV at the physical pion mass.
We calculate the light hadron spectrum in full QCD using two plus one flavor Asqtad sea quarks and domain wall valence quarks. Meson and baryon masses are calculated on a lattice of spatial size $L approx 2.5$texttt{fm}, and a lattice spacing of $a a
pprox 0.124$texttt{fm}, for pion masses as light as $m_pi approx 300$texttt{MeV}, and compared with the results by the MILC collaboration with Asqtad valence quarks at the same lattice spacing. Two- and three-flavor chiral extrapolations of the baryon masses are performed using both continuum and mixed-action heavy baryon chiral perturbation theory. Both the three-flavor and two-flavor functional forms describe our lattice results, although the low-energy constants from the next-to-leading order SU(3) fits are inconsistent with their phenomenological values. Next-to-next-to-leading order SU(2) continuum formulae provide a good fit to the data and yield and extrapolated nucleon mass consistent with experiment, but the convergence pattern indicates that even our lightest pion mass may be at the upper end of the chiral regime. Surprisingly, our nucleon masses are essentially lineaer in $m_pi$ over our full range of pion masses, and we show this feature is common to all recent dynamical calculations of the nucleon mass. The origin of this linearity is not presently understood, and lighter pion masses and increased control of systematic errors will be needed to resolve this puzzling behavior.
