No Arabic abstract
We have performed fits of the pseudoscalar masses and decay constants, from a variety of RBC-UKQCD domain wall fermion ensembles, to $SU(2)$ partially quenched chiral perturbation theory at next-to leading order (NLO) and next-to-next-to leading order (NNLO). We report values for 9 NLO and 8 linearly independent combinations of NNLO partially quenched low energy constants, which we compare to other lattice and phenomenological determinations. We discuss the size of successive terms in the chiral expansion and use our large set of low energy constants to make predictions for mass splittings due to QCD isospin breaking effects and the S-wave $pi pi$ scattering lengths. We conclude that, for the range of pseudoscalar masses explored in this work, $115~mathrm{MeV} lesssim m_{rm PS} lesssim 430~mathrm{MeV}$, the NNLO $SU(2)$ expansion is quite robust and can fit lattice data with percent-scale accuracy.
We measure the pion mass and decay constant on ensembles generated by the Wuppertal-Budapest Collaboration, and extract the NLO low-energy constants l_3 and l_4 of SU(2) chiral perturbation theory. The data are generated in 2+1 flavor simulations with Symanzik glue and 2-fold stout-smeared staggered fermions, with pion masses varying from 135 MeV to 400 MeV, lattice scales between 0.7 GeV and 2.0 GeV, and m_s kept at its physical value. Furthermore, by excluding the lightest mass points, we are able to test the reliability of SU(2) chPT as a tool to extrapolate towards the physical point from higher pion masses.
We have simulated QCD using 2+1 flavors of domain wall quarks on a $(2.74 {rm fm})^3$ volume with an inverse lattice scale of $a^{-1} = 1.729(28)$ GeV. The up and down (light) quarks are degenerate in our calculations and we have used four values for the ratio of light quark masses to the strange (heavy) quark mass in our simulations: 0.217, 0.350, 0.617 and 0.884. We have measured pseudoscalar meson masses and decay constants, the kaon bag parameter $B_K$ and vector meson couplings. We have used SU(2) chiral perturbation theory, which assumes only the up and down quark masses are small, and SU(3) chiral perturbation theory to extrapolate to the physical values for the light quark masses. While next-to-leading order formulae from both approaches fit our data for light quarks, we find the higher order corrections for SU(3) very large, making such fits unreliable. We also find that SU(3) does not fit our data when the quark masses are near the physical strange quark mass. Thus, we rely on SU(2) chiral perturbation theory for accurate results. We use the masses of the $Omega$ baryon, and the $pi$ and $K$ mesons to set the lattice scale and determine the quark masses. We then find $f_pi = 124.1(3.6)_{rm stat}(6.9)_{rm syst} {rm MeV}$, $f_K = 149.6(3.6)_{rm stat}(6.3)_{rm syst} {rm MeV}$ and $f_K/f_pi = 1.205(0.018)_{rm stat}(0.062)_{rm syst}$. Using non-perturbative renormalization to relate lattice regularized quark masses to RI-MOM masses, and perturbation theory to relate these to $bar{rm MS}$ we find $ m_{ud}^{bar{rm MS}}(2 {rm GeV}) = 3.72(0.16)_{rm stat}(0.33)_{rm ren}(0.18)_{rm syst} {rm MeV}$ and $m_{s}^{bar{rm MS}}(2 {rm GeV}) = 107.3(4.4)_{rm stat}(9.7)_{rm ren}(4.9)_{rm syst} {rm MeV}$.
We have performed global fits of $f_{pi}$ and $m_{pi}$, from a variety of RBC-UKQCD domain wall fermion ensembles, to $SU(2)$ partially quenched chiral perturbation theory at NNLO. We report values for 9 NLO and 8 linearly independent combinations of NNLO partially quenched low energy constants, which we compare to other lattice and phenomenological determinations. We discuss the convergence of the expansion and use our large set of low energy constants to make predictions for the pion mass splitting due to QCD isospin breaking effects and the s-wave $pi pi$ scattering lengths. We conclude that, for the range of pseudoscalar masses explored in this work, $115~mathrm{MeV} lesssim m_{rm PS} lesssim 430~mathrm{MeV}$, the NNLO $SU(2)$ expansion is quite robust and can fit lattice data with percent-scale accuracy.
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 results for K to pi and K to 0 matrix elements to next-to-leading order in 2+1 flavor partially quenched chiral perturbation theory. Results are presented for both the Delta I=1/2 and 3/2 channels, for chiral operators corresponding to current-current, gluonic penguin, and electroweak penguin 4-quark operators. These formulas are useful for studying the chiral behavior of currently available 2+1 flavor lattice QCD results, from which the low energy constants of the chiral effective theory can be determined. The low energy constants of these matrix elements are necessary for an understanding of the Delta I=1/2 rule, and for calculations of epsilon/epsilon using current lattice QCD simulations.