
We present results of calculations of the hadronic vacuum polarisation
contribution to the muon anomalous magnetic moment. Specifically, we focus on
controlling the infrared regime of the vacuum polarisation function. Our
results are corrected for finitesize effects by combining the GounarisSakurai
parameterisation of the timelike pion form factor with the L\"uscher formalism.
The impact of quarkdisconnected diagrams and the precision of the scale
determination is discussed and included in our final result in twoflavour QCD,
which carries an overall uncertainty of 6%. We present preliminary results
computed on ensembles with $N_f=2+1$ dynamical flavours and discuss how the
longdistance contribution can be accurately constrained by a dedicated
spectrum calculation in the isovector channel.

We present a calculation of the hadronic vacuum polarization contribution to
the muon anomalous magnetic moment, $a_\mu^{\mathrm hvp}$, in lattice QCD
employing dynamical up and down quarks. We focus on controlling the infrared
regime of the vacuum polarization function. To this end we employ several
complementary approaches, including Pad\'e fits, time moments and the
timemomentum representation. We correct our results for finitevolume effects
by combining the GounarisSakurai parameterization of the timelike pion form
factor with the L\"uscher formalism. On a subset of our ensembles we have
derived an upper bound on the magnitude of quarkdisconnected diagrams and
found that they decrease the estimate for $a_\mu^{\mathrm hvp}$ by at most 2%.
Our final result is $a_\mu^{\mathrm hvp}=(654\pm32\,{}^{+21}_{23})\cdot
10^{10}$, where the first error is statistical, and the second denotes the
combined systematic uncertainty. Based on our findings we discuss the prospects
for determining $a_\mu^{\mathrm hvp}$ with subpercent precision.

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 results for the nucleon electromagnetic form factors, including
the momentum transfer dependence and derived quantities (charge radii and
magnetic moment). The analysis is performed using O(a) improved Wilson fermions
in Nf=2 QCD measured on the CLS ensembles. Particular focus is placed on a
systematic evaluation of the influence of excited states in threepoint
correlation functions, which lead to a biased evaluation, if not accounted for
correctly. We argue that the use of summed operator insertions and fit
ans\"atze including excited states allow us to suppress and control this
effect. We employ a novel method to perform joint chiral and continuum
extrapolations, by fitting the form factors directly to the expressions of
covariant baryonic chiral effective field theory. The final results for the
charge radii and magnetic moment from our lattice calculations include, for the
first time, a full error budget. We find that our estimates are compatible with
experimental results within their overall uncertainties.

We present a computation of Bmeson decay constants from lattice QCD
simulations within the framework of Heavy Quark Effective Theory for the
bquark. The nexttoleading order corrections in the HQET expansion are
included nonperturbatively. Based on Nf=2 gauge field ensembles, covering
three lattice spacings a (0.080.05)fm and pion masses down to 190MeV, a
variational method for extracting hadronic matrix elements is used to keep
systematic errors under control. In addition we perform a careful
autocorrelation analysis in the extrapolation to the continuum and to the
physical pion mass limits. Our final results read fB=186(13)MeV, fBs=224(14)MeV
and fBs/fB=1.203(65). A comparison with other results in the literature does
not reveal a dependence on the number of dynamical quarks, and effects from
truncating HQET appear to be negligible.

We present updated results for the nucleon axial charge and electromagnetic
(EM) form factors, which include a significant increase in statistics for all
ensembles (up to 4000 measurements), as well as the addition of ensembles with
pion masses down to $m_\pi\sim195$ MeV. We also present results for the average
quark momentum fraction. The new data allows us to perform a thorough study of
the systematic effects encountered in the lattice extraction. We concentrate on
systematic effects due to excitedstate contaminations for each of the
quantities, which we check using several different time separations between the
operators at the source and sink through a comparison of plateau fits and the
summed operator insertion method (which provides a mechanism to suppress the
excitedstate contamination). We confirm our earlier finding that a reliable
extraction of the axial charge must be based on a method which eliminates
excitedstate contaminations. Similar conclusions apply to our EM form factor
calculations . The measurements are calculated using the CLS ensembles with
nonperturbatively O(a) improved Wilson fermions in $N_f=2$ QCD.

We report our final estimate of the bquark mass from $N_f=2$ lattice QCD
simulations using Heavy Quark Effective Theory nonperturbatively matched to
QCD at $O(1/m_h)$. Treating systematic and statistical errors in a conservative
manner, we obtain $\overline{m}_{\rm b}^{\overline{\rm MS}}(2 {\rm
GeV})=4.88(15)$ GeV after an extrapolation to the physical point.

We report the final results of the ALPHA collaboration for some Bphysics
observables: $f_B$, $f_{B_s}$ and $m_b$. We employ CLS configurations with 2
flavors of $O(a)$ improved Wilson fermions in the sea and pion masses ranging
down to 190 MeV. The bquark is treated in HQET to order $1/m_b$. The
renormalization, the matching and the improvement were performed
nonperturbatively, and three lattice spacings reaching $a=0.048$ fm are used
in the continuum extrapolation.

The 2012 PDG reports a tension at the level of $3 \sigma$ between two
exclusive determinations of $V_{ub}$. They are obtained by combining the
experimental branching ratios of $B \to \tau \nu$ and $B \to \pi l \nu$
(respectively) with a theoretical computation of the hadronic matrix elements
$\fB$ and the $B \to \pi$ form factor $f_+(q^2)$. To understand the tension,
improved precision and a careful analysis of the systematics involved are
necessary. We report the results of the ALPHA collaboration for $\fB$ from the
lattice with 2 flavors of $O(a)$ improved Wilson fermions. We employ HQET,
including $1/m_b$ corrections, with pion masses ranging down to $\approx$ 190
MeV. Renormalization and matching were performed nonperturbatively, and three
lattice spacings reaching $a^{1}\approx 4.1$ GeV are used in the continuum
extrapolation. We also present progress towards a computation of $f_+(q^2)$, to
directly compare two independent exclusive determinations of $V_{ub}$ with
each other and with inclusive determinations. Additionally, we report on
preliminary results for $\fBq{s}$, needed for the analysis of $B_s \to
\mu^+\mu^$.}

We present our analysis of B physics quantities using nonperturbatively
matched Heavy Quark Effective Theory (HQET) in Nf= 2 lattice QCD on the CLS
ensembles. Using alltoall propagators, HYPsmeared static quarks, and the
Generalized Eigenvalue Problem (GEVP) approach with a conservative plateau
selection procedure, we are able to systematically control all sources of
error. With significantly increased statistics compared to last year, our
preliminary results are mb (mb) = 4.22(10)(4)z GeV for the MS bquark mass, and
fB = 193(9)stat (4)_\chi MeV and fBs = 219(12)stat MeV for the Bmeson decay
constants.

We report on our calculation of the nucleon axial charge gA in QCD with two
flavours of dynamical quarks. A detailed investigation of systematic errors is
performed, with a particular focus on contributions from excited states to
threepoint correlation functions. The use of summed operator insertions allows
for a much better control over such contamination. After performing a chiral
extrapolation to the physical pion mass, we find gA=1.223 +/ 0.063 (stat)
+0.035 0.060 (syst), in good agreement with the experimental value.

Lattice simulations of QCD have produced precise estimates for the masses of
the lowestlying hadrons which show excellent agreement with experiment. By
contrast, lattice results for the vector and axial vector form factors of the
nucleon show significant deviations from their experimental determination. We
present results from our ongoing project to compute a variety of form factors
with control over all systematic uncertainties. In the case of the pion
electromagnetic form factor we employ partially twisted boundary conditions to
extract the pion charge radius directly from the linear slope of the form
factor near vanishing momentum transfer. In the nucleon sector we focus
specifically on the possible contamination from contributions of higher excited
states. We argue that summed correlation functions offer the possibility of
eliminating this source of systematic error. As an illustration of the method
we discuss our results for the axial charge, gA, of the nucleon.

We present an update on our ongoing project to compute hadronic observables
for Nf=2 flavours of O(a) improved Wilson fermions at small lattice spacings.
The procedure to determine the lattice scale via the mass of the Omega baryon
is described. Furthermore we present preliminary results for the pion form
factor computed using twisted boundary conditions, and report on the
implementation of a novel approach to determine the contribution of the
hadronic vacuum polarisation to the anomalous magnetic moment of the muon.

We report on our ongoing project to compute mesonic and baryonic two and
threepoint correlation functions in simulations using Nf=2 flavours of O(a)
improved Wilson quarks and the Wilson plaquette action. We present performance
figures for the DDHMC algorithm on commodity cluster hardware and discuss the
issue of critical slowing down, which is particularly pronounced for the
topological charge. The effectiveness of stochastic noise sources and Jacobi
smearing are investigated. Our preliminary results obtained at three quark
masses on 96x48^3 at beta=5.5 imply that the lattice spacing is about 0.06 fm,
while the smallest pion mass in the current runs is around 360 MeV, which
corresponds to m_pi*L=5.3.

We analyze the systematic errors made when using the generalized eigenvalue
problem to extract energies and matrix elements in lattice gauge theory.
Effective theories such as HQET are also discussed. Numerical results are shown
for the extraction of groundstate and excited Bmeson masses and the
groundstate decay constant in the static approximation.

We discuss some large effects of dynamical fermions. One is a cutoff effect,
others concern the contribution of multipion states to correlation functions
and are expected to survive the continuum limit. We then turn to the
preparation for simulations at small lattice spacings which we are planning
down to around a=0.04fm in order to understand the size of O(a^2)effects of
the standard O(a)improved theory. The dependence of the lattice spacing on the
bare coupling is determined through the Schr"odinger functional renormalized
coupling.

We show that, in four dimensions, the quark propagator is affected by
roundoff errors for large values of the quark mass am and the time extent T/a
even when double precision arithmetics is used to compute it. We introduce a
definition of the solver residual which is sensitive to the problem and apply a
Schwarz alternating procedure to compute the propagator in a number of
subdomains in the time direction. The effectiveness of the method is
demonstrated in a numerical computation of the free onedimensional Dirac
propagator.

We introduce an alternative discretization for static quarks on the lattice
retaining the O(a) improvement properties of the EichtenHill action. In this
formulation, statistical fluctuations are reduced by a factor which grows
exponentially with Euclidean time, x_0. For the first time, Bmeson correlation
functions are computed with good statistical precision in the static
approximation for x_0>1 fm. At lattice spacings a \approx 0.1 fm, a \approx
0.08 fm and a \approx 0.07 fm the B_smeson decay constant is determined in
static and quenched approximations. A correction due to the finite mass of the
bquark is estimated by combining these static results with a recent
determination of F_Ds.

We investigate the lowlying eigenvalues of the improved WilsonDirac
operator in the Schroedinger functional with two dynamical quark flavors. At a
lattice spacing of approximately 0.1 fm we find more very small eigenvalues
than in the quenched case. These cause problems with HMCtype algorithms and in
the evaluation of fermionic correlation functions. Through a simulation at a
finer lattice spacing we are able to establish their nature as cutoffeffects

We address the locality problem arising in simulations, which take the square
root of the staggered fermion determinant as a Boltzmann weight to reduce the
number of dynamical quark tastes from four to two. We study analytically and
numerically the square root of the staggered fermion operator as a candidate to
define a two taste theory from first principles. Although it has the correct
weight, this operator is nonlocal in the continuum limit. Our work serves as a
warning that fundamental properties of field theories might be violated when
employing blindly the square root trick. The question, whether a local operator
reproducing the square root of the staggered fermion determinant exists, is
left open.

We address the locality problem arising in simulations, which take the square
root of the staggered fermion determinant as a Boltzmann weight to reduce the
number of dynamical quark tastes. A definition of such a theory necessitates an
underlying local fermion operator with the same determinant and the
corresponding Green's functions to establish causality and unitarity. We
illustrate this point by studying analytically and numerically the square root
of the staggered fermion operator. Although it has the correct weight, this
operator is nonlocal in the continuum limit. Our work serves as a warning that
fundamental properties of field theories might be violated when employing
blindly the square root trick. The question, whether a local operator
reproducing the square root of the staggered fermion determinant exists, is
left open.

We present and discuss results for cutoff effects in the PCAC masses and the
mass dependence of r_0 for full QCD and various fermion actions. Our discussion
of how one computes mass dependences  here of r_0  is also relevant for
comparisons with chiral perturbation theory.

We present two determinations of the coefficient $c_A$ in quenched QCD,
needed to build the $O(a)$ improved axial current. The first condition used is
the requirement that the PCAC quark mass, as a function of $x_0$, stays flat
for a nontrivial spatial phase for the fermions in the Schr\"odinger
functional. The second condition is that the PCAC relation for the groundstate
and the first excited state at finite L give the same quark mass. Our results
confirm previous findings that in the quenched theory the intrinsic $O(a)$
ambiguity of $c_A$ gets relevant around $\beta\simeq6.0$.

We present a computation of the decay constant f_Bs in quenched QCD. Our
strategy is to combine new precise data from the static approximation with an
interpolation of the decay constant around the charm quark mass region. This
computation is the first step in demonstrating the feasability of a strategy
for f_B in full QCD. The continuum limits in the static theory and at finite
mass are taken separately and will be further improved.

We present an algorithmic study for the simulation of two massless flavors of
O(a) improved Wilson quarks with Schroedinger functional boundary conditions.
The algorithm used is Hybrid Monte Carlo with two pseudofermion fields as
proposed by M. Hasenbusch. A gain in CPU cost of a factor two is reached when
compared to one pseudofermion field due to the larger possible stepsize. This
study is integrated in the ALPHA project for the computation of the running of
the renormalized quark mass. We include an update on these physics results.