
We present a lattice QCD spectroscopy study in the isospin singlet,
strangeness $2$ sectors relevant for the conjectured $H$ dibaryon. We employ
both hexaquark and twobaryon interpolating operators to isolate the ground
state in the rest frame and in moving frames. Calculations are performed using
two flavors of O($a$)improved Wilson fermions and a quenched strange quark.
Our initial pointsource method for constructing correlators does not allow for
twobaryon operators at the source; nevertheless, results from using these
operators at the sink indicate that they provide an improved overlap onto the
ground state in comparison with the hexaquark operators. We also present
results, in the rest frame, using a second method based on distillation to
compute a hermitian matrix of correlators with twobaryon operators at both the
source and the sink. This method yields a much more precise and reliable
determination of the groundstate energy. In the flavorSU(3) symmetric case,
we apply L\"uscher's finitevolume quantization condition to the restframe and
movingframe energy levels to determine the $S$wave scattering phase shift,
near and below the twoparticle threshold. For a pion mass of 960 MeV, we find
that there exists a bound $H$ dibaryon with binding energy
${\Delta}E=(19\pm10)$ MeV. In the 27plet (dineutron) sector, the finitevolume
analysis suggests that the existence of a bound state is unlikely.

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.

A new QCD sum rule determination of the leading order hadronic vacuum
polarization contribution to the anomalous magnetic moment of the muon,
$a_{\mu}^{\rm hvp}$, is proposed. This approach combines data on $e^{+}e^{}$
annihilation into hadrons, perturbative QCD and lattice QCD results for the
first derivative of the electromagnetic current correlator at zero momentum
transfer, $\Pi_{\rm EM}^\prime(0)$. The idea is based on the observation that,
in the relevant kinematic domain, the integration kernel $K(s)$, entering the
formula relating $a_{\mu}^{\rm hvp}$ to $e^{+}e^{}$ annihilation data, behaves
like $1/s$ times a very smooth function of $s$, the squared energy. We find an
expression for $a_{\mu}$ in terms of $\Pi_{\rm EM}^\prime(0)$, which can be
calculated in lattice QCD. Using recent lattice results we find a good
approximation for $a_{\mu}^{\rm hvp}$, but the precision is not yet sufficient
to resolve the discrepancy between the $R(s)$ databased results and the
experimentally measured value.

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 preliminary results on the axial form factor $G_A(Q^2)$ and the
induced pseudoscalar form factor $G_P(Q^2)$ of the nucleon. A systematic
analysis of the excitedstate contributions to form factors is performed on the
CLS ensemble `N6' with $m_\pi = 340 \ \text{MeV}$ and lattice spacing $a \sim
0.05 \ \text{fm}$. The relevant threepoint functions were computed with
sourcesink separations ranging from $t_s \sim 0.6 \ \text{fm}$ to $t_s \sim \
1.4 \ \text{fm}$. We observe that the form factors suffer from nontrivial
excitedstate contributions at the sourcesink separations available to us. It
is noted that naive plateau fits underestimate the excitedstate contributions
and that the method of summed operator insertions correctly accounts for these
effects.

We study the ability of a variety of fitting techniques to extract the ground
state matrix elements of the vector current from ratios of nucleon three and
twopoint functions that contain contaminations from excited states. Extending
our highstatistics study of nucleon form factors, we are able to demonstrate
that the treatment of excitedstate contributions in conjunction with
approaching the physical pion mass has a significant impact on the
$Q^2$dependence of the form factors.

We highlight the progress, current status, and open challenges of QCDdriven
physics, in theory and in experiment. We discuss how the strong interaction is
intimately connected to a broad sweep of physical problems, in settings ranging
from astrophysics and cosmology to stronglycoupled, complex systems in
particle and condensedmatter physics, as well as to searches for physics
beyond the Standard Model. We also discuss how success in describing the strong
interaction impacts other fields, and, in turn, how such subjects can impact
studies of the strong interaction. In the course of the work we offer a
perspective on the many research streams which flow into and out of QCD, as
well as a vision for future developments.

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

We investigate the efficiency of single timeslice stochastic sources for the
calculation of light meson masses on the lattice as one varies the quark mass.
Simulations are carried out with Nf = 2 flavours of nonperturbatively O(a)
improved Wilson fermions for pion masses in the range of 450  760 MeV. Results
for pseudoscalar and vector meson twopoint correlation functions computed
using stochastic as well as point sources are presented and compared. At fixed
computational cost the stochastic approach reduces the variance considerably in
the pseudoscalar channel for all simulated quark masses. The vector channel is
more affected by the intrinsic stochastic noise. In order to obtain stable
estimates of the statistical errors and a more pronounced plateau for the
effective vector meson mass, a relatively large number of stochastic sources
must be used.

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 investigate twopoint correlation functions of lefthanded currents
computed in quenched lattice QCD with the NeubergerDirac operator. We consider
two lattice spacings a~0.09,0.12 fm and two different lattice extents L~ 1.5,
2.0 fm; quark masses span both the p and the epsilonregimes. We compare the
results with the predictions of quenched chiral perturbation theory, with the
purpose of testing to what extent the effective theory reproduces quenched QCD
at low energy. In the pregime we test volume and quark mass dependence of the
pseudoscalar decay constant and mass; in the epsilonregime, we investigate
volume and topology dependence of the correlators. While the leading order
behaviour predicted by the effective theory is very well reproduced by the
lattice data in the range of parameters that we explored, our numerical data
are not precise enough to test nexttoleading order effects.

A new method to determine the lowenergy couplings of the $\Delta S=1$ weak
Hamiltonian is presented. It relies on a matching of the topological poles in
$1/m^2$ of threepoint correlators of two pseudoscalar densities and a
fourfermion operator, measured in lattice QCD, to the same observables
computed in the $\epsilon$regime of chiral perturbation theory. We test this
method in a theory with a light charm quark, i.e. with an SU(4) flavour
symmetry. Quenched numerical measurements are performed in a 2 fm box, and
chiral perturbation theory predictions are worked out up to nexttoleading
order. The matching of the two sides allows to determine the weak lowenergy
couplings in the SU(4) limit. We compare the results with a previous
determination, based on threepoint correlators containing two lefthanded
currents, and discuss the merits and drawbacks of the two procedures.

We discuss a new method to determine the lowenergy couplings of the $\Delta
S=1$ weak Hamiltonian in the $\epsilon$regime. It relies on a matching of the
topological poles in $1/m^2$ of threepoint functions of two pseudoscalar
densities and a fourfermion operator computed in lattice QCD, to the same
observables in the Chiral Effective Theory. We present the results of a NLO
computation in chiral perturbation theory of these correlation functions
together with some preliminary numerical results.

We present results for the nonperturbative renormalisation of fourfermion
operators with two flavours of dynamical quarks. We consider both fully
relativistic left currentleft current operators, and a full basis for $\Delta
B=2$ operators with static heavy quarks. The renormalisation group running of
the operators to high energy scales is computed in the continuum limit for a
family of Schroedinger Functional renormalisation schemes, via standard finite
size scaling techniques. The total renormalisation factors relating
renormalisation group invariant to bare operators are computed for a choice of
lattice regularisations.

We compute the leadingorder lowenergy constants of the DeltaS=1 effective
weak Hamiltonian in the quenched approximation of QCD with up, down, strange,
and charm quarks degenerate and light. They are extracted by comparing the
predictions of finite volume chiral perturbation theory with lattice QCD
computations of suitable correlation functions carried out with quark masses
ranging from a few MeV up to half of the physical strange mass. We observe a
large DeltaI=1/2 enhancement in this corner of the parameter space of the
theory. Although matching with the experimental result is not observed for the
DeltaI=1/2 amplitude, our computation suggests large QCD contributions to the
physical DeltaI=1/2 rule in the GIM limit, and represents the first step to
quantify the role of the charm quarkmass in K>pipi amplitudes.

We outline a general strategy for the nonperturbative renormalisation of
composite operators in discretisations based on Neuberger fermions, via a
matching to results obtained with Wilsontype fermions. As an application, we
consider the renormalisation of the fourquark operators entering the Delta S=1
and Delta S=2 effective Hamiltonians. Our results are an essential ingredient
for the determination of the lowenergy constants governing nonleptonic kaon
decays.

We discuss the renormalisation properties of the complete set of $\Delta B =
2$ fourquark operators with the heavy quark treated in the static
approximation. We elucidate the role of heavy quark symmetry and other symmetry
transformations in constraining their mixing under renormalisation. By
employing the Schroedinger functional, a set of nonperturbative
renormalisation conditions can be defined in terms of suitable correlation
functions. As a first step in a fully nonperturbative determination of the
scaledependent renormalisation factors, we evaluate these conditions in
lattice perturbation theory at one loop. Thereby we verify the expected mixing
patterns and determine the anomalous dimensions of the operators at NLO in the
Schroedinger functional scheme. Finally, by employing twistedmass QCD it is
shown how finite subtractions arising from explicit chiral symmetry breaking
can be avoided completely.

We present our preliminary results for threepoint correlation functions
involving the operators entering the $\Delta{S}=1$ effective Hamiltonian with
an active charm quark, obtained using overlap fermions in the quenched
approximation. This is the first computation carried out for valence quark
masses small enough so as to permit a matching to Quenched Chiral Perturbation
Theory in the $\epsilon$regime. The commonly observed large statistical
fluctuations are tamed by means of lowmode averaging techniques, combined with
restrictions to individual topological sectors. We also discuss the matching of
the resulting hadronic matrix elements to the effective lowenergy constants
for $\Delta{S}=1$ transitions. This involves (a) finitevolume corrections
which can be evaluated at NLO in Quenched Chiral Perturbation Theory, and (b)
the shortdistance renormalization of the relevant fourquark operators in
discretizations based on the overlap operator. We discuss perturbative
estimates for the renormalization factors and possible strategies for their
nonperturbative evaluation. Our results can be used to isolate the
longdistance contributions to the $\Delta I=1/2$ rule, coming from physics
effects around the intrinsic QCD scale.

We present a strategy designed to separate several possible origins of the
wellknown enhancement of the Delta{I}=1/2 amplitude in nonleptonic kaon
decays. In particular, we seek to disentangle the contribution of physics at
the typical QCD scale (softgluon exchange) from the effects at the scale of
the charm quark mass. This is achieved by considering QCD with an unphysically
light charm quark, so that the theory possesses an approximate SU(4)_L x
SU(4)_R chiral symmetry. By computing the relevant operator matrix elements and
monitoring their values as the charm quark mass departs from the
SU(4)symmetric situation, the role of the charm quark can be assessed. We
study the influence of the charm quark mass in Chiral Perturbation Theory.
First results from lattice simulations in the SU(4)symmetric limit are also
discussed.

We report on recent work on the determination of lowenergy constants
describing Delta{S}=1 weak transitions, in order to investigate the origins of
the Delta{I}=1/2 rule. We focus on numerical techniques designed to enhance the
statistical signal in threepoint correlation functions computed with overlap
fermions near the chiral limit.