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

We report initial nucleon structure results computed on lattices with 2+1
dynamical M\"obius domain wall fermions at the physical point generated by the
RBC and UKQCD collaborations. At this stage, we evaluate only connected quark
contributions. In particular, we discuss the nucleon vector and axialvector
form factors, nucleon axial charge and the isovector quark momentum fraction.
From currently available statistics, we estimate the stochastic accuracy of the
determination of $g_A$ and $<x>_{ud}$ to be around 10%, and we expect to
reduce that to 5% within the next year. To reduce the computational cost of our
calculations, we extensively use acceleration techniques such as loweigenmode
deflation and allmodeaveraging (AMA). We present a method for choosing
optimal AMA parameters.

We present lattice QCD calculations of nucleon electromagnetic form factors
using pion masses $m_\pi$ = 149, 202, and 254 MeV and an action with
cloverimproved Wilson quarks coupled to smeared gauge fields, as used by the
BudapestMarseilleWuppertal collaboration. Particular attention is given to
removal of the effects of excited state contamination by calculation at three
sourcesink separations and use of the summation and generalized
penciloffunction methods. The combination of calculation at the nearly
physical mass $m_\pi$ = 149 MeV in a large spatial volume ($m_\pi L_s$ = 4.2)
and removal of excited state effects yields agreement with experiment for the
electric and magnetic form factors $G_E(Q^2)$ and $G_M(Q^2)$ up to $Q^2$ = 0.5
GeV$^2$.

We report the first Lattice QCD calculation using the almost physical pion
mass mpi=149 MeV that agrees with experiment for four fundamental isovector
observables characterizing the gross structure of the nucleon: the Dirac and
Pauli radii, the magnetic moment, and the quark momentum fraction. The key to
this success is the combination of using a nearly physical pion mass and
excluding the contributions of excited states. An analogous calculation of the
nucleon axial charge governing beta decay has inconsistencies indicating a
source of bias at low pion masses not present for the other observables and
yields a result that disagrees with experiment.

Recent results for moments of multiplicity distributions of netprotons and
netelectric charge from the STAR collaboration are compared to lattice QCD
results for higher order fluctuations of baryon number and electric charge by
the WuppertalBudapest collaboration, with the purpose of extracting the
freezeout temperature and chemical potential. All lattice simulations are
performed for a system of 2+1 dynamical quark flavors, at the physical mass for
light and strange quarks; all results are continuum extrapolated. We show that
it is possible to extract an upper value for the freezeout temperature, as
well as precise baryochemical potential values corresponding to the four
highest collision energies of the experimental beam energy scan. Consistency
between the freezeout parameters obtained from baryon number and electric
charge fluctuations is found. The freezeout chemical potentials are now in
agreement with the statistical hadronization model.

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.

We present continuum extrapolated lattice results for the higher order
fluctuations of conserved charges in high temperature Quantum Chromodynamics.
Through the matching of the grand canonical ensemble on the lattice to the net
charge and net baryon distribution realized in heavy ion experiments the
temperature and the chemical potential may be estimated at the time of chemical
freezeout

In lattice QCD computations a substantial amount of work is spent in solving
linear systems arising in Wilson's discretization of the Dirac equations. We
show first numerical results of the extension of the twolevel DD\alpha AMG
method to a true multilevel method based on our parallel MPIC implementation.
Using additional levels pays off, allowing to cut down the core minutes spent
on one system solve by a factor of approximately 700 compared to standard
Krylov subspace methods and yielding another speedup of a factor of 1.7 over
the twolevel approach.

We present our results for ratios of higher order fluctuations of electric
charge as functions of the temperature. These results are obtained in a system
of 2+1 quark flavors at physical quark masses and continuum extrapolated. We
compare them to preliminary data on higher order moments of the net electric
charge distribution from the STAR collaboration. This allows us to determine
the freezeout temperature and chemical potential from first principles. We
also show continuumextrapolated results for ratios of higher order
fluctuations of baryon number. These will allow to test the consistency of the
approach, by comparing them to the corresponding experimental data (once they
become available) and thus extracting the freezeout parameters in an
independent way.

We consider the effect of strong external electromagnetic fields on
thermodynamic observables in QCD, through lattice simulations with 1+1+1
flavors of staggered quarks at physical quark masses. Continuum extrapolated
results are presented for the light quark condensates and for their tensor
polarizations, as functions of the temperature and the magnetic field. We find
the light condensates to undergo inverse magnetic catalysis in the transition
region, in a manner that the transition temperature decreases with growing
magnetic field. We also compare the results to other approaches and lattice
simulations. Furthermore, we relate the tensor polarization to the spin part of
the magnetic susceptibility of the QCD vacuum, and show that this contribution
is diamagnetic.

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.

We present 2+1 flavor Lattice QCD calculations of the nucleon scalar and
tensor charges. Using the BMW cloverimproved Wilson action with pion masses
between 150 and 350 MeV and three sourcesink separations between 0.9 and 1.4
fm, we achieve good control over excitedstate contamination and extrapolation
to the physical pion mass. As a consistency check, we also present results from
calculations using unitary domain wall fermions with pion masses between 300
and 400 MeV, and using domain wall valence quarks and staggered sea quarks with
pion masses between 300 and 600 MeV.

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.

We determine the equation of state of QCD for nonzero chemical potentials via
a Taylor expansion of the pressure. The results are obtained for N_f=2+1
flavors of quarks with physical masses, on various lattice spacings. We present
results for the pressure, interaction measure, energy density, entropy density,
and the speed of sound for small chemical potentials. At low temperatures we
compare our results with the Hadron Resonance Gas model. We also express our
observables along trajectories of constant entropy over particle number. A
simple parameterization is given (the Matlab/Octave script parameterization.m,
submitted to the arXiv along with the paper), which can be used to reconstruct
the observables as functions of T and mu, or as functions of T and S/N.

We study QCD thermodynamics using two flavors of dynamical overlap fermions
with quark masses corresponding to a pion mass of 350 MeV. We determine several
observables on N_t=6 and 8 lattices. All our runs are performed with fixed
global topology. Our results are compared with staggered ones and a nice
agreement is found.

We use the Wilson flow to define the gauge anisotropy at a given physical
scale. We demonstrate the use of the anisotropic flow by performing the tuning
of the bare gauge anisotropy in the treelevel Symanzik action for several
lattice spacings and target anisotropies. We use this method to tune the
anisotropy parameters in full QCD, where we also exploit the diminishing effect
of a well chosen smearing on the renormalization of the fermion anisotropy.

The effect of an external (electro)magnetic field on the finite temperature
transition of QCD is studied. We generate configurations at various values of
the quantized magnetic flux with $N_f=2+1$ flavors of stout smeared staggered
quarks, with physical masses. Thermodynamic observables including the chiral
condensate and susceptibility, and the strange quark number susceptibility are
measured as functions of the field strength. We perform the renormalization of
the studied observables and extrapolate the results to the continuum limit
using $N_t=6,8$ and 10 lattices. We also check for finite volume effects using
various lattice volumes. We find from all of our observables that the
transition temperature $T_c$ significantly decreases with increasing magnetic
field. This is in conflict with various model calculations that predict an
increasing $T_c(B)$. From a finite volume scaling analysis we find that the
analytic crossover that is present at B=0 persists up to our largest magnetic
fields $eB \approx 1 \textmd{GeV}^2$, and that the transition strength
increases mildly up to this $eB\approx1 \textmd{GeV}^2$.

The effect of an external magnetic field on the finite temperature transition
of QCD is studied. We measure thermodynamic observables including the quark
condensates and susceptibilities and the strange quark number susceptibility.
We generate configurations at various values of the quantized magnetic flux
with N_f=2+1 flavors of stout smeared staggered quarks at physical quark
masses. We perform the renormalization of our observables and approach the
continuum limit with N_t=6,8 and 10 lattices. We also check finite volume
effects using various lattice volumes. Our main result is that the transition
temperature significantly decreases with growing magnetic field, and that the
transition remains an analytic crossover up to our largest external field eB ~
1 GeV^2.

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.