
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 compute the slope and curvature, at vanishing fourmomentum transfer
squared, of the leading order hadron vacuum polarization function, using
lattice QCD. Calculations are performed with 2+1+1 flavors of staggered
fermions directly at the physical values of the quark masses and in volumes of
linear extent larger than 6fm. The continuum limit is carried out using six
different lattice spacings. All connected and disconnected contributions are
calculated, up to and including those of the charm.

We review lattice results related to pion, kaon, D and Bmeson physics with
the aim of making them easily accessible to the particle physics community.
More specifically, we report on the determination of the lightquark masses,
the form factor f+(0), arising in the semileptonic K > pi transition at zero
momentum transfer, as well as the decay constant ratio fK/fpi and its
consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe
the results obtained on the lattice for some of the lowenergy constants of
SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the
determination of the BK parameter of neutral kaon mixing as well as the
additional four B parameters that arise in theories of physics beyond the
Standard Model. The latter quantities are an addition compared to the previous
review. For the heavyquark sector, we provide results for mc and mb (also new
compared to the previous review), as well as those for D and Bmeson decay
constants, form factors, and mixing parameters. These are the heavyquark
quantities most relevant for the determination of CKM matrix elements and the
global CKM unitaritytriangle fit. Finally, we review the status of lattice
determinations of the strong coupling constant alpha_s.

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).

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 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.

Scale setting is of central importance in lattice QCD. It is required to
predict dimensional quantities in physical units. Moreover, it determines the
relative lattice spacings of computations performed at different values of the
bare coupling, and this is needed for extrapolating results into the continuum.
Thus, we calculate a new quantity, $w_0$, for setting the scale in lattice QCD,
which is based on the Wilson flow like the scale $t_0$ (M. Luscher, JHEP 1008
(2010) 071). It is cheap and straightforward to implement and compute. In
particular, it does not involve the delicate fitting of correlation functions
at asymptotic times. It typically can be determined on the few permil level.
We compute its continuum extrapolated value in 2+1flavor QCD for physical and
nonphysical pion and kaon masses, to allow for massindependent scale setting
even away from the physical mass point. We demonstrate its robustness by
computing it with two very different actions (one of them with staggered, the
other with Wilson fermions) and by showing that the results agree for physical
quark masses in the continuum limit.

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.

We study the spectra of heavylight and heavyheavy mesons containing charm
quarks, including higher spin states. We use two sets of $N_f = 2 + 1$ gauge
configurations, one set from QCDSF using the SLiNC action, and the other
configurations from the BudapestMarseilleWuppertal collaboration, using the
HEX smeared clover action. To extract information about the excited states, we
choose a suitable basis of operators to implement the variational method.

Ordinary matter is described by six fundamental parameters: three couplings
(gravitational, electromagnetic and strong) and three masses: the electron's
(m_e) and those of the up (m_u) and down (m_d) quarks. An additional mass
enters through quantum fluctuations: the strange quark mass (m_s). The three
couplings and m_e are known with an accuracy of better than a few per mil.
Despite their importance, $m_u$, $m_d$ (their average m_{ud}) and m_s are
relatively poorly known: e.g. the Particle Data Group quotes them with
conservative errors close to 25%. Here we determine these quantities with a
precision below 2% by performing ab initio lattice quantum chromodynamics (QCD)
calculations, in which all systematics are controlled. We use pion and quark
masses down to (and even below) their physical values, lattice sizes of up to 6
fm, and five lattice spacings to extrapolate to continuum spacetime. All
necessary renormalizations are performed nonperturbatively.

Based on a series of lattice calculations we determine the ratio FK/Fpi in
QCD. With experimental data from kaon decay and nuclear double beta decay, we
obtain a precise determination of Vus. Our simulation includes 2+1 flavours
of sea quarks, with three lattice spacings, large volumes and a simulated pion
mass reaching down to about 190 MeV for a full control over the systematic
uncertainties.

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.

We give details of our precise determination of the light quark masses
m_{ud}=(m_u+m_d)/2 and m_s in 2+1 flavor QCD, with simulated pion masses down
to 120 MeV, at five lattice spacings, and in large volumes. The details concern
the action and algorithm employed, the HMC force with HEX smeared clover
fermions, the choice of the scale setting procedure and of the input masses.
After an overview of the simulation parameters, extensive checks of algorithmic
stability, autocorrelation and (practical) ergodicity are reported. To
corroborate the good scaling properties of our action, explicit tests of the
scaling of hadron masses in N_f=3 QCD are carried out. Details of how we
control finite volume effects through dedicated finite volume scaling runs are
reported. To check consistency with SU(2) Chiral Perturbation Theory the
behavior of M_\pi^2/m_{ud} and F_\pi as a function of m_{ud} is investigated.
Details of how we use the RI/MOM procedure with a separate continuum limit of
the running of the scalar density R_S(\mu,\mu') are given. This procedure is
shown to reproduce the known value of r_0m_s in quenched QCD. Input from
dispersion theory is used to split our value of m_{ud} into separate values of
m_u and m_d. Finally, our procedure to quantify both systematic and statistical
uncertainties is discussed.

We determine the ratio FK/Fpi in QCD with Nf=2+1 flavors of sea quarks, based
on a series of lattice calculations with three different lattice spacings,
large volumes and a simulated pion mass reaching down to about 190 MeV. We
obtain FK/Fpi=1.192 +/ 0.007(stat) +/ 0.006(syst). This result is then used
to give an updated value of the CKM matrix element Vus. The unitarity
relation for the first row of this matrix is found to be well observed.

More than 99% of the mass of the visible universe is made up of protons and
neutrons. Both particles are much heavier than their quark and gluon
constituents, and the Standard Model of particle physics should explain this
difference. We present a full abinitio calculation of the masses of protons,
neutrons and other light hadrons, using lattice quantum chromodynamics. Pion
masses down to 190 mega electronvolts are used to extrapolate to the physical
point with lattice sizes of approximately four times the inverse pion mass.
Three lattice spacings are used for a continuum extrapolation. Our results
completely agree with experimental observations and represent a quantitative
confirmation of this aspect of the Standard Model with fully controlled
uncertainties.

We present a framework for phenomenological lattice QCD calculations which
makes use of a tree level Symanzink improved action for gluons and stoutlink
Wilson fermions. We give details of our efficient HMC/RHMC algorithm and
present a scaling study of the lowlying N_f=3 baryon spectrum. We find a
scaling region that extends to a~<0.16fm and conclude that our action and
algorithm are suitable for large scale phenomenological investigations of
N_f=2+1 QCD. We expect this conclusion to hold for other comparable actions.

Some algorithmic details of our $N_f=2+1$ QCD mixed action simulations with
overlap valence and improved Wilson sea quarks are presented.

Using quenched and unquenched chiral perturbation theory we compute vector
and axial current twopoint functions at finite volume and fixed gauge field
topology, in the socalled \epsilonregime of QCD. A comparison of these
results with finite volume lattice calculations allows to determine the
parameters of the corresponding chiral Lagrangians.

Using quenched chiral perturbation theory, we compute the longdistance
behaviour of twopoint functions of flavour nonsinglet axial and vector
currents in a finite volume, for small quark masses, and at a fixed gaugefield
topology. We also present the corresponding predictions for the unquenched
theory at fixed topology. These results can in principle be used to measure the
lowenergy constants of the chiral Lagrangian, from lattice simulations in
volumes much smaller than one pion Compton wavelength. We show that quenching
has a dramatic effect on the vector correlator, which is argued to vanish to
all orders, while the axial correlator appears to be a robust observable only
moderately sensitive to quenching.

We present the results of a lattice computation of the form factors for B^0
>pi^ l^+nu_l decays near zerorecoil. These results will allow a
determination of the CKM matrix element Vub when measurements of the
differential decay rate become available. We also provide models for
extrapolation of the form factors and rate to the full recoil range. Our
computation is performed in the quenched approximation to QCD on a 24^3x48
lattice at beta=6.2, using a nonperturbatively O(a)improved action. The
masses of all light valence quarks involved are extrapolated to their physical
values.