Electromagnetic mass splittings of the low lying hadrons and quark masses from 2+1 flavor lattice QCD+QED

Results computed in lattice QCD+QED are presented for the electromagnetic mass splittings of the low-lying hadrons. These are used to determine the renormalized, nondegenerate, light quark masses. It is found that m_u^MS¯=2.24(10)(34), m_d^MS¯=4.65(15)(32), and m_s^MS¯=97.6(2.9)(5.5)MeV at the renormalization scale 2 GeV, where the first error is statistical and the second systematic. We find the lowest-order electromagnetic splitting (m_π⁺-m_π⁰)_QED=3.38(23)MeV, the splittings including next-to-leading order, (m_π⁺-m_π⁰)_QED=4.50(23)MeV, (m_K⁺-m_K⁰)_QED=1.87(10)MeV, and the m_u≠m_d contribution to the kaon mass difference, (m_K⁺-m_K⁰)_(mu-md)=-5.840(96)MeV. All errors are statistical only, and the next-to-leading-order pion splitting is only approximate in that it does not contain all next-to-leading-order contributions. We also computed the proton-neutron mass difference, including for the first time, QED interactions in a realistic 2+1 flavor calculation. We find (m_p-m_n)_QED=0.383(68)MeV, (m_p-m_n)_(mu-md)=-2.51(14)MeV (statistical errors only), and the total m_p-m_n=-2.13(16)(70)MeV, where the first error is statistical, and the second, part of the systematic error. The calculations are carried out on QCD ensembles generated by the RBC and UKQCD collaborations, using domain wall fermions and the Iwasaki gauge action (gauge coupling β=2.13 and lattice cutoff a^-1≈1.78GeV). We use two lattice sizes, 16³ and 24³ ((1.8fm)³ and (2.7fm)³), to address finite-volume effects. Noncompact QED is treated in the quenched approximation. The valence pseudoscalar meson masses in our study cover a range of about 250 to 700 MeV, though we use only those up to about 400 MeV to quote final results. We present new results for the electromagnetic low-energy constants in SU(3) and SU(2) partially quenched chiral perturbation theory to the next-to-leading order, obtained from fits to our data. Detailed analysis of systematic errors in our results and methods for improving them are discussed. Finally, new analytic results for SU(2)_L×SU(2)_R-plus-kaon chiral perturbation theory, including the one-loop logs proportional to α_emm, are given.