
We present a lattice calculation of the nucleon isovector axial and induced
pseudoscalar form factors on the CLS ensembles using $N_{\rm f}=2$ dynamical
flavours of nonperturbatively $\mathcal{O}(a)$improved Wilson fermions and an
$\mathcal{O}(a)$improved axial current together with the pseudoscalar density.
Excitedstate effects in the extraction of the form factors are treated using a
variety of methods, with a detailed discussion of their respective merits. The
chiral and continuum extrapolation of the results is performed both using
formulae inspired by Heavy Baryon Chiral Perturbation Theory (HBChPT) and a
global approach to the form factors based on a chiral effective theory (EFT)
including axial vector mesons. Our results indicate that careful treatment of
excitedstate effects is important in order to obtain reliable results for the
axial form factors of the nucleon, and that the main remaining error stems from
the systematic uncertainties of the chiral extrapolation. As final results, we
quote $g_{\rm A} = 1.278 \pm 0.068\genfrac{}{}{0pt}{1}{+0.000}{0.087}$,
$\langle r_{\rm A}^2\rangle = 0.360 \pm
0.036\genfrac{}{}{0pt}{1}{+0.080}{0.088}~\mathrm{fm}^2$, and $g_{\rm P} = 7.7
\pm 1.8 \genfrac{}{}{0pt}{1}{+0.8}{2.0}$ for the axial charge, axial charge
radius and induced pseudoscalar charge, respectively, where the first error is
statistical and the second is systematic.

In the framework of quantum chromodynamics (QCD), parton distribution
functions (PDFs) quantify how the momentum and spin of a hadron are divided
among its quark and gluon constituents. Two main approaches exist to determine
PDFs. The first approach, based on QCD factorization theorems, realizes a QCD
analysis of a suitable set of hardscattering measurements, often using a
variety of hadronic observables. The second approach, based on firstprinciple
operator definitions of PDFs, uses lattice QCD to compute directly some
PDFrelated quantities, such as their moments. Motivated by recent progress in
both approaches, in this document we present an overview of latticeQCD and
globalanalysis techniques used to determine unpolarized and polarized proton
PDFs and their moments. We provide benchmark numbers to validate present and
future latticeQCD calculations and we illustrate how they could be used to
reduce the PDF uncertainties in current unpolarized and polarized global
analyses. This document represents a first step towards establishing a common
language between the two communities, to foster dialogue and to further improve
our knowledge of PDFs.

We briefly review several activities at Mainz related to hadronic
lightbylight scattering (HLbL) using lattice QCD. First we present a
positionspace approach to the HLbL contribution in the muon g2, where we
focus on exploratory studies of the pionpole contribution in a simple model
and the lepton loop in QED in the continuum and in infinite volume. The second
part describes a lattice calculation of the doublevirtual pion transition form
factor F_{pi^0 gamma^* gamma^*}(q_1^2, q_2^2) in the spacelike region with
photon virtualities up to 1.5 GeV^2 which paves the way for a lattice
calculation of the pionpole contribution to HLbL. The third topic involves
HLbL forward scattering amplitudes calculated in lattice QCD which can be
described, using dispersion relations (HLbL sum rules), by gamma^* gamma^* >
hadrons fusion cross sections and then compared with phenomenological models.

We present our preliminary results on the calculation of hadronic
lightbylight forward scattering amplitudes using vector fourpoint
correlation functions computed on the lattice. Using a dispersive approach,
forward scattering amplitudes can be described by $\gamma^* \gamma^* \to$
hadrons fusion cross sections and then compared with phenomenology. We show
that only a few states are needed to reproduce our data. In particular, the sum
rules considered in this study imply relations between meson$\gamma\gamma$
couplings and provide valuable information about individual form factors which
are often used to estimate the mesonpole contributions to the hadronic
lightbylight contribution to the $(g2)$ of the muon.

We report on an analysis of the average quark momentum fraction of the
nucleon and related quantities using $N_\mathrm{f}=2+1$ Wilson fermions.
Computations are performed on four CLS ensembles covering three values of the
lattice spacing at pion masses down to $M_\pi \approx 200\,\mathrm{MeV}$.
Several sourcesink separations ($\sim 1.0\,\mathrm{fm}$ to $\sim
1.4\,\mathrm{fm}$) are used to assess the excitedstate contamination. To gain
further insight, the generalized penciloffunctions approach has been
implemented to reduce the excitedstate contamination in the relevant two and
threepoint functions. Preliminary results are shown for the isovector nucleon
charges from vector, axial vector and tensor derivative (twist2) operators.

We present an exploratory study of leading isospin breaking effects in
mesonic masses using $O(a)$ improved Wilson fermions. Isospin symmetry is
explicitly broken by distinct masses and electric charges of the up and down
quarks. In order to be able to make use of existing isosymmetric QCD gauge
ensembles we apply reweighting techniques. The path integral describing QCD+QED
is expanded perturbatively in powers of the light quarks' mass deviations and
the electromagnetic coupling. We employ QED$_{\mathrm{L}}$ as a finite volume
formulation of QED.

Results are presented from an ongoing study of the $\rho$ resonance. The
focus is on CLS 2flavour ensembles generated using $\mathcal{O}(a)$ improved
Wilson fermions with pion masses ranging from $265$ to $437$ $\mathrm{MeV}$.
The energy levels are extracted by solving the GEVP of correlator matrices,
created with the distillation approach involving $\rho$ and $\pi\pi$
interpolators. The study is done in the centreofmass frame and several moving
frames. One aim of this work is to extract the timelike pion form factor after
applying the L\"uscher formalism. We therefore plan to integrate this study
with the existing Mainz programme for the calculation of the hadronic vacuum
polarization contribution to the muon $g2$.

We present updated results on the nucleon electromagnetic form factors and
axial coupling calculated using CLS ensembles with $N_\mathrm{f}=2+1$ dynamical
flavours of Wilson fermions. The measurements are performed on large, fine
lattices with a pseudoscalar mass reaching down to 200 MeV. The
truncatedsolver method is employed to reduce the variance of the measurements.
Estimation of the matrix elements is challenging due to large contamination
from excited states and further investigation is necessary to bring these
effects under control.

We study the thermal transition of QCD with two degenerate light flavours by
lattice simulations using $O(a)$improved Wilson quarks. Temperature scans are
performed at a fixed value of $N_t = (aT)^{1}=16$, where $a$ is the lattice
spacing and $T$ the temperature, at three fixed zerotemperature pion masses
between 200 MeV and 540 MeV. In this range we find that the transition is
consistent with a broad crossover. As a probe of the restoration of chiral
symmetry, we study the static screening spectrum. We observe a degeneracy
between the transverse isovector vector and axialvector channels starting from
the transition temperature. Particularly striking is the strong reduction of
the splitting between isovector scalar and pseudoscalar screening masses around
the chiral phase transition by at least a factor of three compared to its value
at zero temperature. In fact, the splitting is consistent with zero within our
uncertainties. This disfavours a chiral phase transition in the $O(4)$
universality class.

Results are presented from an ongoing study of the $\rho$ resonance. We use
the distillation approach in order to create correlator matrices involving
$\rho$ and $\pi\pi$ interpolators. The study is done in a centreofmass frame
and several moving frames. We are able to extract energy levels by solving the
GEVP of those correlator matrices. The initial exploratory study is being done
on a CLS 2flavour lattice with a pion mass of $451$ $\mathrm{MeV}$ using
$\mathcal{O}(a)$ improved Wilson fermions. One aim of this work is to extract
the timelike pion form factor after applying the L\"uscher formalism. We also
plan to integrate this study with the existing Mainz programme for the
calculation of the hadronic vacuum polarization contribution to the muon $g2$
and will extend our study to lower pion masses and larger lattices in the
future, including ensembles with $2+1$ flavours.

We study the performance of allmodeaveraging (AMA) when used in conjunction
with a locally deflated SAPpreconditioned solver, determining how to optimize
the local block sizes and number of deflation fields in order to minimize the
computational cost for a given level of overall statistical accuracy. We find
that AMA enables a reduction of the statistical error on nucleon charges by a
factor of around two at the same cost when compared to the standard method. As
a demonstration, we compute the axial, scalar and tensor charges of the nucleon
in $N_f=2$ lattice QCD with nonperturbatively O(a)improved Wilson quarks,
using O(10,000) measurements to pursue the signal out to sourcesink
separations of $t_s\sim 1.5$ fm. Our results suggest that the axial charge is
suffering from a significant amount (510%) of excitedstate contamination at
sourcesink separations of up to $t_s\sim 1.2$ fm, whereas the excitedstate
contamination in the scalar and tensor charges seems to be small.

We present an update of our lattice QCD study of the vacuum polarization
function using O$(a)$improved $N_ {\rm f} =2$ Wilson fermions with increased
statistics and a large set of momenta. The resulting points are highly
correlated and thus require a correlated fitting procedure. We employ an
extended frequentist method to estimate the systematic uncertainties due to the
momentum dependence and to the continuum and chiral extrapolations. We present
preliminary results of the leading order hadronic contribution of the anomalous
magnetic moment of the muon $\left(a_\mu^{\mathrm{HLO}}\right)$ at the physical
point for $u,d,s$ and $c$ valence quarks.

We present preliminary results on the electromagnetic form factors and axial
charge of the nucleon from ensembles generated by the CLS effort with
$N_\mathrm{f}=2+1$ flavours of nonperturbatively $\mathrm{O}(a)$improved
Wilson fermions and open temporal boundary conditions. Systematic effects due
to excitedstate contamination are accounted for using both twostate fits and
the method of summed operator insertions. This exploratory analysis
demonstrates the viability of obtaining precision baryon observables with
$N_\mathrm{f}=2+1$ flavours of Wilson fermions on fine lattices, aiming towards
controlled chiral and continuum limits in the future.

The electromagnetic coupling receives significant contributions to its
running from nonperturbative QCD effects. We present an update of a lattice
QCD study of the Adler function and of its application to the determination of
the running of the QED coupling. We perform a highstatistics computation with
two flavours of O$(a)$ improved Wilson fermions in a large range of momentum
transfer $Q^2$. The running of the electromagnetic coupling, including
contributions from $u$, $d$, $s$ and $c$ valence quarks, is compared to
phenomenological determinations at intermediate $Q^2$ values. An extension of
this study to the determination of the hadronic contributions to the running of
the weak mixing angle is also described.

We present preliminary results from a lattice QCD calculation of the
Hdibaryon using two flavors of $\mathcal{O}(a)$ improved Wilson fermions. We
employ local sixquark interpolating operators at the source with a combination
of local sixquark and twobaryon operators at the sink with the appropriate
quantum numbers of the Hdibaryon and its coupling to the twobaryon channels.
We find that the twobaryon operators provide an improved overlap onto the
ground state in comparison to the local sixquark operators. We also apply
L\"uscher's finite volume formalism to obtain information on the nature of the
infinitevolume interaction of two particles. Further, the momentum projection
to three moving frames enables the isolation of the pole in the infinitevolume
scattering amplitude. Preliminary results at pion masses of 450 MeV and 1 GeV
clearly indicate the presence of states below the $\Lambda \Lambda$ threshold
while a finitevolume analysis fails to conclusively show the existence of an
infinitevolume bound state.

We extend our study of the pion scalar radius in twoflavour lattice QCD to
include two additional lattice spacings as well as lighter pion masses,
enabling us to perform a combined chiral and continuum extrapolation. We find
discretisation artefacts to be small for the radius, and confirm the importance
of the disconnected diagrams in reproducing the correct chiral behaviour. Our
final result for the scalar radius of the pion at the physical point is
$\left\langle r^2\right\rangle^\pi_{\rm S}=0.600\pm0.052$ fm$^2$, corresponding
to a value of $\overline{\ell}_4=4.54\pm0.30$ for the lowenergy constant
$\overline{\ell}_4$ of NLO chiral perturbation theory.

The anomalous magnetic moment of the muon, $g_\mu2$, is one of the most
promising observables to identify hints for physics beyond the Standard Model.
QCD contributions are currently responsible for the largest fraction of the
overall theoretical uncertainty in $g_\mu2$. The possibility to determine
these hadronic contributions from first principles through lattice QCD
calculations has triggered a number of recent studies. Recent proposals to
improve the accuracy of lattice determinations are reported. We present an
update of our studies of the leadingorder hadronic contribution to $g_\mu2$
with improved Wilson fermions.

The contribution from hadronic vacuum polarisation effects is responsible for
a large fraction of the theoretical uncertainty in the running of the QED
coupling. The current level of uncertainty has become a limitation for
electroweak precision tests. We use lattice QCD simulations with two flavours
of O$(a)$ improved Wilson fermions to determine the Adler function in a broad
range of the momentum transfer $Q^2$. The running of the QED coupling,
including valence contributions from $u$, $d$, $s$ and $c$ quarks, is compared
to phenomenological results at intermediate $Q^2$ values. In the large $Q^2$
regime, the lattice determination of the Adler function is fitted to
perturbation theory in order to examine the feasibility of a determination of
the strong coupling constant.

The hadronic vacuum polarization can be determined from the vector correlator
in a mixed timemomentum representation. We explicitly calculate the
disconnected contribution to the vector correlator, both in the $N_f = 2$
theory and with an additional quenched strange quark, using nonperturbatively
$O(a)$improved Wilson fermions. Alltoall propagators are computed using
stochastic sources and a generalized hopping parameter expansion. Combining the
result with the dominant connected contribution, we are able to estimate an
upper bound for the systematic error that arises from neglecting the
disconnected contribution in the determination of $(g2)_\mu$.

We present early results from a lattice QCD study seeking a bound
$H$dibaryon using $N_f=2$ flavors of $O(a)$ improved Wilson fermions and a
quenched strange quark. We compute a matrix of twopoint functions using
operators consisting of the two independent local products of six
positiveparityprojected quarks with the appropriate quantum numbers, which
belong to the singlet and 27plet irreducible representations of flavor SU(3).
To expand this basis, we also independently vary the quarkfield smearing, and
apply a new scheme to reduce the noise caused by smearing. We then find the
groundstate mass by solving the generalized eigenvalue problem. We show
results from ensembles with pion masses 451 MeV and 1 GeV, and compare with
other lattice calculations.

We compute the Adler function on the lattice from vacuum polarization data
with twisted boundary conditions using numerical derivatives. The study is
based on CLS ensembles with two flavours of $O(a)$ improved Wilson fermions. We
extrapolate the lattice data for the Adler function to the continuum limit and
to the physical pion mass and analyze its dependence on the momentum transfer.
We discuss the application of this method to the extraction of the $u,d$
contribution to $a_\mu^{\mathrm{HLO}}$.

We extend our analysis of the leading hadronic contribution to the anomalous
magnetic moment of the muon using the mixed representation method to study its
chiral behavior. We present results derived from localconserved twopoint
lattice vector correlation functions, computed on a subset of light twoflavor
ensembles made available to us through the CLS effort with pion masses as low
as 190 MeV. The data is analyzed also using the more standard fourmomentum
method. Both methods are systematically compared as the calculations approach
the physical point.

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 semileptonic K > pi transition at zero
momentum transfer, as well as the decay constant ratio fK/fpi of decay
constants 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 and review the determination of the BK parameter of neutral kaon mixing.
The inclusion of heavyquark quantities significantly expands the FLAG scope
with respect to the previous review. Therefore, for this review, we focus on D
and Bmeson decay constants, form factors, and mixing parameters, since these
are most relevant for the determination of CKM matrix elements and the global
CKM unitaritytriangle fit. In addition we review the status of lattice
determinations of the strong coupling constant alpha_s.

We calculate the scalar form factor of the pion using two dynamical flavors
of nonperturbatively $\mathcal{O}(a)$improved Wilson fermions, including both
the connected and the disconnected contribution to the relevant correlation
functions. We employ the calculation of alltoall propagators using stochastic
sources and a generalized hopping parameter expansion. From the form factor
data at vanishing momentum transfer, $Q^2=0$, and two nonvanishing $Q^2$ we
obtain an estimate for the scalar radius $\left<r^2\right>^\pi_{_{\rm S}}$ of
the pion at one value of the lattice spacing and for five different pion
masses. Using Chiral Perturbation Theory at nexttoleading order, we find
$\left<r^2\right>^\pi_{_{\rm S}}=0.635\pm0.016$ fm$^2$ at the physical pion
mass (statistical error only). This is in good agreement with the
phenomenological estimate from $\pi\pi$scattering. The inclusion of the
disconnected contribution is essential for achieving this level of agreement.

We have performed the first numerical study of minimally doubled fermions of
the KarstenWilczek class in the quenched approximation. This requires fixing
the counterterms, which arise due to hypercubic symmetry breaking induced by
the KarstenWilczek term. Nonperturbative renormalisation criteria are
formulated after a detailed study of the parameter dependence of mesonic
observables. Minimisation of the mass anisotropy of the pseudoscalar ground
state fixes nonperturbative renormalisation conditions for the counterterm
coefficients. These anisotropies are mapped out by probing different euclidean
components of the transfer matrix through calculations of the pseudoscalar
ground state mass in different directions. The chiral behaviour of the
pseudoscalar ground state is studied with the tuned KarstenWilczek action for
multiple lattice spacings. Light pseudoscalar masses ($ M_{PS} \lesssim
250\,MeV $) were achieved in the quenched approximation without encountering
exceptional configurations. The presence of quenched chiral logarithms is
studied under the tentative assumption of Goldstone Bosonlike behaviour.