• 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 two-baryon 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 point-source method for constructing correlators does not allow for two-baryon 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 two-baryon operators at both the source and the sink. This method yields a much more precise and reliable determination of the ground-state energy. In the flavor-SU(3) symmetric case, we apply L\"uscher's finite-volume quantization condition to the rest-frame and moving-frame energy levels to determine the $S$-wave scattering phase shift, near and below the two-particle 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 27-plet (dineutron) sector, the finite-volume 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 finite-size effects by combining the Gounaris-Sakurai parameterisation of the timelike pion form factor with the L\"uscher formalism. The impact of quark-disconnected diagrams and the precision of the scale determination is discussed and included in our final result in two-flavour 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 long-distance contribution can be accurately constrained by a dedicated spectrum calculation in the iso-vector 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 time-momentum representation. We correct our results for finite-volume effects by combining the Gounaris-Sakurai 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 quark-disconnected 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 sub-percent 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)$ data-based results and the experimentally measured value.
  • 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 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 low-energy 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 heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle 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 three-point 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 excited-state 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 three-point functions were computed with source-sink 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 non-trivial excited-state contributions at the source-sink separations available to us. It is noted that naive plateau fits underestimate the excited-state 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 two-point functions that contain contaminations from excited states. Extending our high-statistics study of nucleon form factors, we are able to demonstrate that the treatment of excited-state 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 QCD-driven 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 strongly-coupled, complex systems in particle and condensed-matter 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 excited-state 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 excited-state contamination). We confirm our earlier finding that a reliable extraction of the axial charge must be based on a method which eliminates excited-state contaminations. Similar conclusions apply to our EM form factor calculations . The measurements are calculated using the CLS ensembles with non-perturbatively 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 three-point 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 non-perturbatively O(a) improved Wilson fermions for pion masses in the range of 450 - 760 MeV. Results for pseudoscalar and vector meson two-point 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 lowest-lying 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 on-going 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 on-going project to compute mesonic and baryonic two- and three-point 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 DD-HMC 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 two-point correlation functions of left-handed currents computed in quenched lattice QCD with the Neuberger-Dirac 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 epsilon-regimes. 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 p-regime we test volume and quark mass dependence of the pseudoscalar decay constant and mass; in the epsilon-regime, 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 next-to-leading order effects.
  • A new method to determine the low-energy couplings of the $\Delta S=1$ weak Hamiltonian is presented. It relies on a matching of the topological poles in $1/m^2$ of three-point correlators of two pseudoscalar densities and a four-fermion 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 next-to-leading order. The matching of the two sides allows to determine the weak low-energy couplings in the SU(4) limit. We compare the results with a previous determination, based on three-point correlators containing two left-handed currents, and discuss the merits and drawbacks of the two procedures.
  • We discuss a new method to determine the low-energy 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 three-point functions of two pseudoscalar densities and a four-fermion 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 non-perturbative renormalisation of four-fermion operators with two flavours of dynamical quarks. We consider both fully relativistic left current-left 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 leading-order low-energy 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 quark-mass in K-->pipi amplitudes.
  • We outline a general strategy for the non-perturbative renormalisation of composite operators in discretisations based on Neuberger fermions, via a matching to results obtained with Wilson-type fermions. As an application, we consider the renormalisation of the four-quark operators entering the Delta S=1 and Delta S=2 effective Hamiltonians. Our results are an essential ingredient for the determination of the low-energy constants governing non-leptonic kaon decays.
  • We discuss the renormalisation properties of the complete set of $\Delta B = 2$ four-quark 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 non-perturbative renormalisation conditions can be defined in terms of suitable correlation functions. As a first step in a fully non-perturbative determination of the scale-dependent 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 twisted-mass QCD it is shown how finite subtractions arising from explicit chiral symmetry breaking can be avoided completely.
  • We present our preliminary results for three-point 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 low-mode averaging techniques, combined with restrictions to individual topological sectors. We also discuss the matching of the resulting hadronic matrix elements to the effective low-energy constants for $\Delta{S}=1$ transitions. This involves (a) finite-volume corrections which can be evaluated at NLO in Quenched Chiral Perturbation Theory, and (b) the short-distance renormalization of the relevant four-quark operators in discretizations based on the overlap operator. We discuss perturbative estimates for the renormalization factors and possible strategies for their non-perturbative evaluation. Our results can be used to isolate the long-distance 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 well-known enhancement of the Delta{I}=1/2 amplitude in non-leptonic kaon decays. In particular, we seek to disentangle the contribution of physics at the typical QCD scale (soft-gluon 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 low-energy 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 three-point correlation functions computed with overlap fermions near the chiral limit.