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

QCD thermodynamics is considered using Wilson fermions in the fixed scale
approach. The temperature dependence of the renormalized chiral condensate,
quark number susceptibility and Polyakov loop is measured at four lattice
spacings allowing for a controlled continuum limit. The light quark masses are
fixed to heavier than physical values in this first study. Finite volume
effects are ensured to be negligible by using approriately large box sizes. The
final continuum results are compared with staggered fermion simulations
performed in the fixed N_t approach. The same continuum renormalization
conditions are used in both approaches and the final results agree perfectly.

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.

We perform an investigation of the static quarkquarkpotential both in the
confined and the deconfined phase. We discuss conceptual and technical problems
and present first results of an exploratory numerical investigation.

In lattice QCD it is possible, in principle, to determine the parameters in
the effective chiral lagrangian (including weak interaction couplings) by
performing numerical simulations in the $\epsilon$regime, i.e. at quark
masses where the physical extent of the lattice is much smaller than the
Compton wave length of the pion. The use of a formulation of the lattice theory
that preserves chiral symmetry is attractive in this context, but the numerical
implementation of any such approach requires special care in this kinematical
situation due to the presence of some very low eigenvalues of the Dirac
operator. We discuss a set of techniques (lowmode preconditioning and
adaptedprecision algorithms in particular) that make such computations
numerically safe and more efficient by a large factor.

We present the results of an exploratory numerical study of two dimensional
QCD with overlap fermions. We have performed extensive simulations for U(N_c)
and SU(N_c) color groups with N_c=2, 3, 4 and coupling constants chosen to
satisfy the 't Hooft condition g^2 N_c =const=4/3. We have computed the meson
spectrum and decay constants, the topological susceptibility and the chiral
condensate. For U(N_c) gauge groups, our results indicate that the
WittenVeneziano relation is satisfied within our statistical errors and that
the chiral condensate for N_f=1 is compatible with a nonzero value. Our
results exhibit universality in N_c and confirm once more the excellent chiral
properties of the overlapDirac operator.

We present results of a quenched QCD simulation with overlap fermions on a
lattice of volume V = 16^3X32 at beta=6.0, which corresponds approximatively to
a lattice cutoff of 2 GeV and an extension of 1.4 fm. From the twopoint
correlation functions of bilinear operators we extract the pseudoscalar meson
masses and the corresponding decay constants. From the GMOR relation we
determine the chiral condensate and, by using the Kmeson mass as experimental
input, we compute the sum of the strange and average updown quark masses (m_s
+ \hat m). The needed logarithmic divergent renormalization constant Z_S is
computed with the RI/MOM nonperturbative renormalization technique. Since the
overlap preserves chiral symmetry at finite cutoff and volume, no divergent
quark mass and chiral condensate additive renormalizations are required and the
results are O(a) improved.