
We present a firstprinciples lattice QCD+QED calculation at physical pion
mass of the leadingorder 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 Rratio 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, Rratio statistical,
and Rratio systematic errors given separately. This is the currently most
precise determination of the leadingorder hadronic vacuum polarization
contribution to the muon anomalous magnetic moment. In addition, we present the
first lattice calculation of the lightquark 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
electroquenched 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 masssplittings.

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 updown
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 nearlyphysical
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 subpercent statistical errors and
negligible chiral and finitevolume 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
FeynmanHellmann 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 quarkconnected vacuum polarisation
effects. Simulations were performed using RBCUKQCD'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 updown
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
HellmannFeynman 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 longdistance
contributions are very significant. We show that the dominant finitevolume
corrections (those decreasing as powers of the volume) are negligibly small and
that, in the fourflavor 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 longdistance 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 noisereduction 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$ physicalpionmass 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 RBCUKQCD domain wall fermion ensembles, to $SU(2)$ partially
quenched chiral perturbation theory at nextto leading order (NLO) and
nexttonextto 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 Swave $\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 percentscale 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 quantumchromodynamics and quantumelectrodynamics
computations with four nondegenerate Wilson fermion flavors and computed the
neutronproton masssplitting with an accuracy of $300$ kiloelectron 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 longrange nature, they
lead to large finitesize 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 updown 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 updown 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 updown 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 cloverimproved
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 RIMOM 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 2loop conversion
yields B_K^{MSbarNDR}(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
BudapestMarseilleWuppertal collaboration and the Regensburg group to study
the quark massdependence of octet baryons in SU(3) Baryon XPT. This
formulation is expected to extend to larger masses than HeavyBaryon 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.