• We present a lattice calculation of the nucleon iso-vector axial and induced pseudoscalar form factors on the CLS ensembles using $N_{\rm f}=2$ dynamical flavours of non-perturbatively $\mathcal{O}(a)$-improved Wilson fermions and an $\mathcal{O}(a)$-improved axial current together with the pseudoscalar density. Excited-state 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 excited-state 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 hard-scattering measurements, often using a variety of hadronic observables. The second approach, based on first-principle operator definitions of PDFs, uses lattice QCD to compute directly some PDF-related quantities, such as their moments. Motivated by recent progress in both approaches, in this document we present an overview of lattice-QCD and global-analysis techniques used to determine unpolarized and polarized proton PDFs and their moments. We provide benchmark numbers to validate present and future lattice-QCD 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 light-by-light scattering (HLbL) using lattice QCD. First we present a position-space approach to the HLbL contribution in the muon g-2, where we focus on exploratory studies of the pion-pole 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 double-virtual 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 pion-pole 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 light-by-light forward scattering amplitudes using vector four-point 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 meson-pole contributions to the hadronic light-by-light contribution to the $(g-2)$ 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 source-sink separations ($\sim 1.0\,\mathrm{fm}$ to $\sim 1.4\,\mathrm{fm}$) are used to assess the excited-state contamination. To gain further insight, the generalized pencil-of-functions approach has been implemented to reduce the excited-state contamination in the relevant two- and three-point functions. Preliminary results are shown for the isovector nucleon charges from vector, axial vector and tensor derivative (twist-2) 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 2-flavour 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 centre-of-mass 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 $g-2$.
  • 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 truncated-solver 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 zero-temperature 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 axial-vector 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 centre-of-mass 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 2-flavour 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 $g-2$ 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 all-mode-averaging (AMA) when used in conjunction with a locally deflated SAP-preconditioned 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 non-perturbatively O(a)-improved Wilson quarks, using O(10,000) measurements to pursue the signal out to source-sink separations of $t_s\sim 1.5$ fm. Our results suggest that the axial charge is suffering from a significant amount (5-10%) of excited-state contamination at source-sink separations of up to $t_s\sim 1.2$ fm, whereas the excited-state 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 non-perturbatively $\mathrm{O}(a)$-improved Wilson fermions and open temporal boundary conditions. Systematic effects due to excited-state contamination are accounted for using both two-state 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 non-perturbative 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 high-statistics 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 H-dibaryon using two flavors of $\mathcal{O}(a)$ improved Wilson fermions. We employ local six-quark interpolating operators at the source with a combination of local six-quark and two-baryon operators at the sink with the appropriate quantum numbers of the H-dibaryon and its coupling to the two-baryon channels. We find that the two-baryon operators provide an improved overlap onto the ground state in comparison to the local six-quark operators. We also apply L\"uscher's finite volume formalism to obtain information on the nature of the infinite-volume interaction of two particles. Further, the momentum projection to three moving frames enables the isolation of the pole in the infinite-volume 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 finite-volume analysis fails to conclusively show the existence of an infinite-volume bound state.
  • We extend our study of the pion scalar radius in two-flavour 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 low-energy constant $\overline{\ell}_4$ of NLO chiral perturbation theory.
  • The anomalous magnetic moment of the muon, $g_\mu-2$, 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_\mu-2$. 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 leading-order hadronic contribution to $g_\mu-2$ 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 time-momentum 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 non-perturbatively $O(a)$-improved Wilson fermions. All-to-all 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 $(g-2)_\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 two-point functions using operators consisting of the two independent local products of six positive-parity-projected quarks with the appropriate quantum numbers, which belong to the singlet and 27-plet irreducible representations of flavor SU(3). To expand this basis, we also independently vary the quark-field smearing, and apply a new scheme to reduce the noise caused by smearing. We then find the ground-state 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 local-conserved two-point lattice vector correlation functions, computed on a subset of light two-flavor 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 four-momentum method. Both methods are systematically compared as the calculations approach the physical point.
  • We review lattice results related to pion, kaon, D- and B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark 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 low-energy 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 heavy-quark quantities significantly expands the FLAG scope with respect to the previous review. Therefore, for this review, we focus on D- and B-meson decay constants, form factors, and mixing parameters, since these are most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle 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 non-perturbatively $\mathcal{O}(a)$-improved Wilson fermions, including both the connected and the disconnected contribution to the relevant correlation functions. We employ the calculation of all-to-all propagators using stochastic sources and a generalized hopping parameter expansion. From the form factor data at vanishing momentum transfer, $Q^2=0$, and two non-vanishing $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 next-to-leading 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 Karsten-Wilczek class in the quenched approximation. This requires fixing the counterterms, which arise due to hypercubic symmetry breaking induced by the Karsten-Wilczek term. Non-perturbative 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 non-perturbative 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 Karsten-Wilczek 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 Boson-like behaviour.