• We present a first-principles lattice QCD+QED calculation at physical pion mass of the leading-order hadronic vacuum polarization contribution to the muon anomalous magnetic moment. The total contribution of up, down, strange, and charm quarks including QED and strong isospin breaking effects is found to be $a_\mu^{\rm HVP~LO}=715.4(16.3)(9.2) \times 10^{-10}$, where the first error is statistical and the second is systematic. By supplementing lattice data for very short and long distances with experimental R-ratio data using the compilation of Ref. [1], we significantly improve the precision of our calculation and find $a_\mu^{\rm HVP~LO} = 692.5(1.4)(0.5)(0.7)(2.1) \times 10^{-10}$ with lattice statistical, lattice systematic, R-ratio statistical, and R-ratio systematic errors given separately. This is the currently most precise determination of the leading-order hadronic vacuum polarization contribution to the muon anomalous magnetic moment. In addition, we present the first lattice calculation of the light-quark QED correction at physical pion mass.
  • We calculate the strong isospin breaking and QED corrections to meson masses and the hadronic vacuum polarization in an exploratory study on a $64\times24^3$ lattice with an inverse lattice spacing of $a^{-1}=1.78$ GeV and an isospin symmetric pion mass of $m_\pi=340$ MeV. We include QED in an electro-quenched setup using two different methods, a stochastic and a perturbative approach. We find that the electromagnetic correction to the leading hadronic contribution to the anomalous magnetic moment of the muon is smaller than $1\%$ for the up quark and $0.1\%$ for the strange quark, although it should be noted that this is obtained using unphysical light quark masses. In addition to the results themselves, we compare the precision which can be reached for the same computational cost using each method. Such a comparison is also made for the meson electromagnetic mass-splittings.
  • We present a determination of the corrections to Dashen's theorem and of the individual up and down quark masses from a lattice calculation based on quenched QED and $N_f=2+1$ QCD simulations with 5 lattice spacings down to 0.054 fm. The simulations feature lattice sizes up to 6 fm and average up-down quark masses all the way down to their physical value. For the parameter which quantifies violations to Dashens's theorem we obtain $\epsilon=0.73(2)(5)(17)$, where the first error is statistical, the second is systematic, and the third is an estimate of the QED quenching error. For the light quark masses we obtain, $m_u=2.27(6)(5)(4) \, MeV$ and $m_d=4.67(6)(5)(4) \, MeV$ in the $\overline{MS}$ scheme at $2 \, GeV$ and the isospin breaking ratios $m_u/m_d=0.485(11)(8)(14)$, $R=38.2(1.1)(0.8)(1.4)$ and $Q=23.4(0.4)(0.3)(0.4)$. Our results exclude the $m_u=0$ solution to the strong CP problem by more than 24 standard deviations.
  • We present results for several light hadronic quantities ($f_\pi$, $f_K$, $B_K$, $m_{ud}$, $m_s$, $t_0^{1/2}$, $w_0$) obtained from simulations of 2+1 flavor domain wall lattice QCD with large physical volumes and nearly-physical pion masses at two lattice spacings. We perform a short, O(3)%, extrapolation in pion mass to the physical values by combining our new data in a simultaneous chiral/continuum `global fit' with a number of other ensembles with heavier pion masses. We use the physical values of $m_\pi$, $m_K$ and $m_\Omega$ to determine the two quark masses and the scale - all other quantities are outputs from our simulations. We obtain results with sub-percent statistical errors and negligible chiral and finite-volume systematics for these light hadronic quantities, including: $f_\pi$ = 130.2(9) MeV; $f_K$ = 155.5(8) MeV; the average up/down quark mass and strange quark mass in the $\bar {\rm MS}$ scheme at 3 GeV, 2.997(49) and 81.64(1.17) MeV respectively; and the neutral kaon mixing parameter, $B_K$, in the RGI scheme, 0.750(15) and the $\bar{\rm MS}$ scheme at 3 GeV, 0.530(11).
  • We present a QCD calculation of the $u$, $d$ and $s$ scalar quark contents of nucleons based on $47$ lattice ensembles with $N_f = 2+1$ dynamical sea quarks, $5$ lattice spacings down to $0.054\,\text{fm}$, lattice sizes up to $6\,\text{fm}$ and pion masses down to $120\,\text{MeV}$. Using the Feynman-Hellmann theorem, we obtain $f^N_{ud} = 0.0405(40)(35)$ and $f^N_s = 0.113(45)(40)$, which translates into $\sigma_{\pi N}=38(3)(3)\,\text{MeV}$, $\sigma_{sN}=105(41)(37)\,\text{MeV}$ and $y_N=0.20(8)(8)$ for the sigma terms and the related ratio, where the first errors are statistical and the second are systematic. Using isospin relations, we also compute the individual up and down quark contents of the proton and neutron (results in the main text).
  • We present results for the leading hadronic contribution to the muon anomalous magnetic moment due to strange quark-connected vacuum polarisation effects. Simulations were performed using RBC--UKQCD's $N_f=2+1$ domain wall fermion ensembles with physical light sea quark masses at two lattice spacings. We consider a large number of analysis scenarios in order to obtain solid estimates for residual systematic effects. Our final result in the continuum limit is $a_\mu^{(2)\,{\rm had},\,s}=53.1(9)\left(^{+1}_{-3}\right)\times10^{-10}$.
  • In a previous letter (arXiv:1306.2287) we determined the isospin mass splittings of the baryon octet from a lattice calculation based on quenched QED and $N_f{=}2{+}1$ QCD simulations with 5 lattice spacings down to $0.054~\mathrm{fm}$, lattice sizes up to $6~\mathrm{fm}$ and average up-down quark masses all the way down to their physical value. Using the same data we determine here the corrections to Dashen's theorem and the individual up and down quark masses. For the parameter which quantifies violations to Dashens's theorem, we obtain $\epsilon=0.73(2)(5)(17)$, where the first error is statistical, the second is systematic, and the third is an estimate of the QED quenching error. For the light quark masses we obtain, $m_u=2.27(6)(5)(4)~\mathrm{MeV}$ and $m_d=4.67(6)(5)(4)~\mathrm{MeV}$ in the $\bar{\mathrm{MS}}$ scheme at $2~\mathrm{GeV}$ and the isospin breaking ratios $m_u/m_d=0.485(11)(8)(14)$, $R=38.2(1.1)(0.8)(1.4)$ and $Q=23.4(0.4)(0.3)(0.4)$. Our results exclude the $m_u=0$ solution to the strong CP problem by more than $24$ standard deviations.
  • By using lattice QCD computations we determine the sigma terms and strangeness content of all octet baryons by means of an application of the Hellmann-Feynman theorem. In addition to polynomial and rational expressions for the quark mass dependence of octet members, we use SU(3) covariant baryon chiral perturbation theory to perform the extrapolation to the physical up and down quark masses. Our N_f=2+1 lattice ensembles include pion masses down to about 190 MeV in large volumes (M_\pi L > 4), and three values of the lattice spacing. Our main results are the nucleon sigma term \sigma_{\pi N} = 39(4)(^{+18}_{-7}) and the strangeness content y_{N} = 0.20(7)(^{+13}_{-17}). Under the assumption of validity of covariant baryon \chi PT in our range of masses one finds y_{N} = 0.276(77)(^{+90}_{-62}).
  • The rare kaon decays $K\to\pi\ell^+\ell^-$ and $K\to\pi\nu\bar{\nu}$ are flavor changing neutral current (FCNC) processes and hence promising channels with which to probe the limits of the standard model and to look for signs of new physics. In this paper we demonstrate the feasibility of lattice calculations of $K\to\pi\ell^+\ell^-$ decay amplitudes for which long-distance contributions are very significant. We show that the dominant finite-volume corrections (those decreasing as powers of the volume) are negligibly small and that, in the four-flavor theory, no new ultraviolet divergences appear as the electromagnetic current $J$ and the effective weak Hamiltonian $H_W$ approach each other. In addition, we demonstrate that one can remove the unphysical terms which grow exponentially with the range of the integration over the time separation between $J$ and $H_W$. We will now proceed to exploratory numerical studies with the aim of motivating further experimental measurements of these decays. Our work extends the earlier study by Isidori, Turchetti and Martinelli which focussed largely on the renormalization of ultraviolet divergences. In a companion paper we discuss the evaluation of the long-distance contributions to $K\to\pi\nu\bar{\nu}$ decays; these contributions are expected to be at the level of a few percent for $K^+$ decays.
  • We report the first lattice QCD calculation of the hadronic vacuum polarization disconnected contribution to the muon anomalous magnetic moment at physical pion mass. The calculation uses a refined noise-reduction technique which enabled the control of statistical uncertainties at the desired level with modest computational effort. Measurements were performed on the $48^3 \times 96$ physical-pion-mass lattice generated by the RBC and UKQCD collaborations. We find $a_\mu^{\rm HVP~(LO)~DISC} = -9.6(3.3)(2.3)\times 10^{-10}$, where the first error is statistical and the second systematic.
  • 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.
  • The existence and stability of atoms rely on the fact that neutrons are more massive than protons. The measured mass difference is only 0.14\% of the average of the two masses. A slightly smaller or larger value would have led to a dramatically different universe. Here, we show that this difference results from the competition between electromagnetic and mass isospin breaking effects. We performed lattice quantum-chromodynamics and quantum-electrodynamics computations with four nondegenerate Wilson fermion flavors and computed the neutron-proton mass-splitting with an accuracy of $300$ kilo-electron volts, which is greater than $0$ by $5$ standard deviations. We also determine the splittings in the $\Sigma$, $\Xi$, $D$ and $\Xi_{cc}$ isospin multiplets, exceeding in some cases the precision of experimental measurements.
  • Electromagnetic effects are increasingly being accounted for in lattice quantum chromodynamics computations. Because of their long-range nature, they lead to large finite-size effects over which it is important to gain analytical control. Nonrelativistic effective field theories provide an efficient tool to describe these effects. Here we argue that some care has to be taken when applying these methods to quantum electrodynamics in a finite volume.
  • While electromagnetic and up-down quark mass difference effects on octet baryon masses are very small, they have important consequences. The stability of the hydrogen atom against beta decay is a prominent example. Here we include these effects by adding them to valence quarks in a lattice QCD calculation based on $N_f=2+1$ simulations with 5 lattice spacings down to 0.054 fm, lattice sizes up to 6 fm and average up-down quark masses all the way down to their physical value. This allows us to gain control over all systematic errors, except for the one associated with neglecting electromagnetism in the sea. We compute the octet baryon isomultiplet mass splittings, as well as the individual contributions from electromagnetism and the up-down quark mass difference. Our results for the total splittings are in good agreement with experiment.
  • Indirect CP violation in K \rightarrow {\pi}{\pi} decays plays a central role in constraining the flavor structure of the Standard Model (SM) and in the search for new physics. For many years the leading uncertainty in the SM prediction of this phenomenon was the one associated with the nonperturbative strong interaction dynamics in this process. Here we present a fully controlled lattice QCD calculation of these effects, which are described by the neutral kaon mixing parameter B_K . We use a two step HEX smeared clover-improved Wilson action, with four lattice spacings from a\approx0.054 fm to a\approx0.093 fm and pion masses at and even below the physical value. Nonperturbative renormalization is performed in the RI-MOM scheme, where we find that operator mixing induced by chiral symmetry breaking is very small. Using fully nonperturbative continuum running, we obtain our main result B_K^{RI}(3.5GeV)=0.531(6)_{stat}(2)_{sys}. A perturbative 2-loop conversion yields B_K^{MSbar-NDR}(2GeV)=0.564(6)_{stat}(3)_{sys}(6)_{PT}, which is in good agreement with current results from fits to experimental data.
  • At the precision reached in current lattice QCD calculations, electromagnetic effects are becoming numerically relevant. We will present preliminary results for electromagnetic corrections to light hadron masses, based on simulations in which a $\mathrm{U}(1)$ degree of freedom is superimposed on $N_f=2+1$ QCD configurations from the BMW collaboration.
  • A status report is given for a joint project of the Budapest-Marseille-Wuppertal collaboration and the Regensburg group to study the quark mass-dependence of octet baryons in SU(3) Baryon XPT. This formulation is expected to extend to larger masses than Heavy-Baryon XPT. Its applicability is tested with 2+1 flavor data which cover three lattice spacings and pion masses down to about 190 MeV, in large volumes. Also polynomial and rational interpolations in M_\pi^2 and M_K^2 are used to assess the uncertainty due to the ansatz. Both frameworks are combined to explore the precision to be expected in a controlled determination of the nucleon sigma term and strangeness content.
  • While the masses of light hadrons have been extensively studied in lattice QCD simulations, there exist only a few exploratory calculations of the strong decay widths of hadronic resonances. We will present preliminary results of a computation of the rho meson width obtained using $N_f=2+1$ flavor simulations. The work is based on L\"uscher's formalism and its extension to moving frames.