• 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.
  • Within the HMC algorithm, we discuss how, by using the shadow Hamiltonian and the Poisson brackets, one can achieve a simple factorization in the dependence of the Hamiltonian violations upon either the algorithmic parameters or the parameters specifying the integrator. We consider the simplest case of a second order (nested) Omelyan intergrator and one level of Hasenbusch splitting of the determinant for the simulations of a QCD-like theory (with gauge group SU(2)). Given the specific choice of the integrator, the Poisson brackets reduce to the variances of the molecular dynamics forces. We show how the factorization can be used to optimize in a very economical and simple way both the algorithmic and the intergrator parameters with good accuracy.
  • We compute the leading QED corrections to the hadronic vacuum polarization (HVP) of the photon, relevant for the determination of leptonic anomalous magnetic moments, $a_\ell$. We work in the electroquenched approximation and use dynamical QCD configurations generated by the CLS initiative with two degenerate flavors of non-perturbatively O($a$)-improved Wilson fermions. We consider QED$_{\rm L}$ and QED$_{\rm M}$ to deal with the finite-volume zero modes. We compare results for the Wilson loops with exact analytical determinations. In addition we make sure that the volumes and photon masses used in QED$_{\rm M}$ are such that the correct dispersion relation is reproduced by the energy levels extracted from the charged pions two-point functions. Finally we compare results for pion masses and the HVP between QED$_{\rm L}$ and QED$_{\rm M}$. For the vacuum polarization, corrections with respect to the pure QCD case, at fixed pion masses, turn out to be at the percent level.
  • We consider the possibility of using reweighting techniques in order to correct for the breaking of unitarity when twisted boundary conditions are imposed on valence fermions in simulations of lattice gauge theories. We start by studying the properties of reweighting factors and their variances at tree-level. That leads us to the introduction of a factorization for the fermionic reweighting determinant. In the numerical, stochastic, implementation of the method, we find that the effect of reweighting is negligible in the case of large volumes but it is sizeable when the volumes are small and the twisting angles are large. More importantly, we find that for un-improved Wilson fermions, and in small volumes, the dependence of the critical quark mass on the twisting angle is quite pronounced and results in large violations of the continuum dispersion relation.
  • We present a practical strategy to optimize a set of Hybrid Monte Carlo parameters in simulations of QCD and QCD-like theories. We specialize to the case of mass-preconditioning, with multiple time-step Omelyan integrators. Starting from properties of the shadow Hamiltonian we show how the optimal setup for the integrator can be chosen once the forces and their variances are measured, assuming that those only depend on the mass-preconditioning parameter.
  • We non-perturbatively determine the renormalization factor of the axial vector current in lattice QCD with $N_f=3$ flavors of Wilson-clover fermions and the tree-level Symanzik-improved gauge action. The (by now standard) renormalization condition is derived from the massive axial Ward identity and it is imposed among Schr\"{o}dinger functional states with large overlap on the lowest lying hadronic state in the pseudoscalar channel, in order to reduce kinematically enhanced cutoff effects. We explore a range of couplings relevant for simulations at lattice spacings of $\approx 0.09$ fm and below. An interpolation formula for $Z_A(g_0^2)$, smoothly connecting the non-perturbative values to the 1-loop expression, is provided together with our final results.
  • 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.
  • Imposing twisted boundary conditions on the fermionic fields is a procedure extensively used when evaluating, for example, form factors on the lattice. Twisting is usually performed for one flavour and only in the valence, and this causes a breaking of unitarity. In this work we explore the possibility of restoring unitarity through the reweighting method. We first study some properties of the approach at tree level and then we stochastically evaluate ratios of fermionic determinants for different boundary conditions in order to include them in the gauge averages, avoiding in this way the expensive generation of new configurations for each choice of the twisting angle, $\theta$. As expected the effect of reweighting is negligible in the case of large volumes but it is important when the volumes are small and the twisting angles are large. In particular we find a measurable effect for the plaquette and the pion correlation function in the case of $\theta=\pi/2$ in a volume $16\times 8^3$, and we observe a systematic upward shift in the pion dispersion relation.
  • We present a study of the B spectrum performed in the framework of Heavy Quark Effective Theory expanded to next-to-leading order in 1/m and non-perturbative in the strong coupling. Our analyses have been performed on Nf=2 lattice gauge field ensembles corresponding to three different lattice spacings and a wide range of pion masses. We obtain the Bs-meson mass and hyperfine splittings of the B- and Bs-mesons that are in good agreement with the experimental values and examine the mass difference m_{Bs}-m_B as a further cross-check of our previous estimate of the b-quark mass. We also report on the mass splitting between the first excited state and the ground state in the B and Bs systems.
  • We review recent lattice results for quark masses and low-energy hadronic parameters relevant for flavor physics. We do that by describing the FLAG initiative, with emphasis on its scope and rating criteria. The emerging picture is that while for light quantities a large number of computations using different approaches exist, and this increases the overall confidence on the final averages/estimates, in the heavy-light case the field is less advanced and, with the exception of decay constants, only a few computations are available. The precision reached for the light quantities is such that electromagnetic (EM) corrections, beyond the point-like approximation, are becoming relevant. We discuss recent computations of the spectrum based on direct simulations of QED+QCD. We also present theoretical developments for including EM effects in leptonic decays. We conclude describing recent results for the $K \to \pi \pi$ transition amplitudes and prospects for tackling hadronic decays on the lattice.
  • We report on a non-perturbative computation of the renormalization factor Z_A of the axial vector current in three-flavour O(a) improved lattice QCD with Wilson quarks and tree-level Symanzik improved gauge action and also recall our recent determination of the improvement coefficient c_A. Our normalization and improvement conditions are formulated at constant physics in a Schr\"odinger functional setup. The normalization condition exploits the full, massive axial Ward identity to reduce finite quark mass effects in the evaluation of Z_A and correlators with boundary wave functions to suppress excited state contributions in the pseudoscalar channel.
  • We review a lattice strategy how to non-perturbatively determine the coefficients in the HQET expansion of all components of the heavy-light axial and vector currents, including 1/m_h-corrections. We also discuss recent preliminary results on the form factors parameterizing semi-leptonic B-decays at the leading order in 1/m_h.
  • 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 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 present a strategy how to match the full set of components of the heavy-light axial and vector currents in Heavy Quark Effective Theory (HQET), up to and including 1/m-corrections, to QCD. While the ultimate goal is to apply these matching conditions non-perturbatively, in this study we first have implemented them at tree-level, in order to find good choices of the matching observables with small O(1/m^2) contributions. They can later be employed in the non-perturbative matching procedure which is a crucial part of precision HQET computations of semileptonic decay form factors in lattice QCD.
  • We report on an ongoing non-perturbative determination of the improvement coefficient of the axial current, $c_\mathrm A$, with three flavours of dynamical $\mathrm O(a)$ improved Wilson quarks and tree-level Symanzik improved gauge action. Our computations are based on simulations with the openQCD code. The improvement condition for a range of couplings is formulated with Schr\"odinger functional boundary conditions and imposed along a line of constant physics in parameter space. Our analysis involves correlation functions with boundary wave functions such that a large sensitivity to $c_\mathrm A$ can be reached by exploiting the PCAC relation with two different pseudoscalar states.
  • We use a lattice determination of the hadronic vacuum polarization tensor to study the associated Ward identities and compute the Adler function. The vacuum polarization tensor is computed from a combination of point-split and local vector currents, using two flavours of O($a$)-improved Wilson fermions. Partially twisted boundary conditions are employed to obtain a fine momentum resolution. The modifications of the Ward identities by lattice artifacts and by the use of twisted boundary conditions are monitored. We determine the Adler function from the derivative of the vacuum polarization function over a large region of momentum transfer $q^2$. As a first account of systematic effects, a continuum limit scaling analysis is performed in the large $q^2$ regime.
  • We consider a non-perturbative formulation of an SU(2) massive gauge theory on a space-time lattice, which is also a discretised gauged non-linear chiral model. The lattice model is shown to have an exactly conserved global SU(2) symmetry. If a scaling region for the lattice model exists and the lightest degrees of freedom are spin one vector particles with the same quantum numbers as the conserved current, we argue that the most general effective theory describing their low-energy dynamics must be a massive gauge theory. We present results of a exploratory numerical simulation of the model and find indications for the presence of a scaling region where both a triplet vector and a scalar remain light.
  • We report on our ongoing project to calculate the leading hadronic contribution to the anomalous magnetic moment of the muon a_mu^HLO using two dynamical flavours of non-perturbatively O(a) improved Wilson fermions. In this study, we changed the vacuum polarisation tensor to a combination of local and point-split currents which significantly reduces the numerical effort. Partially twisted boundary conditions allow us to improve the momentum resolution of the vacuum polarisation tensor and therefore the determination of the leading hadronic contribution to (g-2)_mu. We also extended the range of ensembles to include a pion mass below 200 MeV which allows us to check the non-trivial chiral behaviour of a_mu^HLO.
  • We present updated preliminary results for the nucleon electromagnetic form factors for non-perturbatively $\mathcal{O}(a)$ improved Wilson fermions in $N_f=2$ QCD measured on the CLS ensembles. The use of the summed operator insertion method allows us to suppress the influence of excited states in our measurements. A study of the effect that excited state contaminations have on the $Q^2$ dependence of the extracted nucleon form factors may then be made through comparisons of the summation method to standard plateau fits, as well as to excited state fits.
  • We demonstrate that a reduction in the noise-to-signal ratio may be obtained for hadrons at non-zero momenta whilst maintaining a good overlap with the ground state through a generalisation of Gaussian/Wuppertal smearing. The use of an anisotropic smearing wavefunction is motivated by the physical picture of a boosted hadron.
  • We generalize Gaussian/Wuppertal smearing in order to produce non-spherical wave functions. We show that we can achieve a reduction in the noise-to-signal ratio for correlation functions of certain hadrons at non-zero momentum, while at the same time preserving a good projection on the ground state.
  • We report on a non-perturbative determination of the parameters of the lattice Heavy Quark Effective Theory (HQET) Lagrangian and of the time component of the heavy-light axial-vector current with Nf=2 flavors of massless dynamical quarks. The effective theory is considered at the 1/mh order, and the heavy mass mh covers a range from slightly above the charm to beyond the beauty region. These HQET parameters are needed to compute, for example, the b-quark mass, the heavy-light spectrum and decay constants in the static approximation and to order 1/mh in HQET. The determination of the parameters is done non-perturbatively. The computation reported in this paper uses the plaquette gauge action and two different static actions for the heavy quark described by HQET. For the light-quark action we choose non-perturbatively O(a)-improved Wilson fermions.
  • We report on our computation of the leading hadronic contribution to the anomalous magnetic moment of the muon using two dynamical flavours of non-perturbatively O(a) improved Wilson fermions. The strange quark is introduced in the quenched approximation. Partially twisted boundary conditions are applied to improve the momentum resolution in the relevant integral. Our results, obtained at three different values of the lattice spacing, allow for a preliminary study of discretization effects. We explore a wide range of lattice volumes, namely 2 fm < L < 3 fm, with pion masses from 600 to 280 MeV and discuss different chiral extrapolations to the physical point. We observe a non-trivial dependence of a_mu(HLO) on m_pi especially for small pion masses. The final result, a_mu(HLO)=618(64)*10^(-10), is obtained by considering only the quark connected contribution to the vacuum polarization. We present a detailed analysis of systematic errors and discuss how they can be reduced in future simulations.
  • We present an updated analysis of the non-perturbatively renormalized b-quark mass and B meson decay constant based on CLS lattices with two dynamical non-perturbatively improved Wilson quarks. This update incorporates additional light quark masses and lattice spacings in large physical volume to improve chiral extrapolations and to reach the continuum limit. We use Heavy Quark Effective Theory (HQET) including 1/m_b terms with non-perturbative coefficients based on the matching of QCD and HQET developed by the ALPHA collaboration during the past years.