• 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.
  • 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 a computation of B-meson decay constants from lattice QCD simulations within the framework of Heavy Quark Effective Theory for the b-quark. The next-to-leading order corrections in the HQET expansion are included non-perturbatively. Based on Nf=2 gauge field ensembles, covering three lattice spacings a (0.08-0.05)fm and pion masses down to 190MeV, a variational method for extracting hadronic matrix elements is used to keep systematic errors under control. In addition we perform a careful autocorrelation analysis in the extrapolation to the continuum and to the physical pion mass limits. Our final results read fB=186(13)MeV, fBs=224(14)MeV and fBs/fB=1.203(65). A comparison with other results in the literature does not reveal a dependence on the number of dynamical quarks, and effects from truncating HQET appear to be negligible.
  • 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 our final estimate of the b-quark mass from $N_f=2$ lattice QCD simulations using Heavy Quark Effective Theory non-perturbatively matched to QCD at $O(1/m_h)$. Treating systematic and statistical errors in a conservative manner, we obtain $\overline{m}_{\rm b}^{\overline{\rm MS}}(2 {\rm GeV})=4.88(15)$ GeV after an extrapolation to the physical point.
  • We report the final results of the ALPHA collaboration for some B-physics observables: $f_B$, $f_{B_s}$ and $m_b$. We employ CLS configurations with 2 flavors of $O(a)$ improved Wilson fermions in the sea and pion masses ranging down to 190 MeV. The b-quark is treated in HQET to order $1/m_b$. The renormalization, the matching and the improvement were performed non-perturbatively, and three lattice spacings reaching $a=0.048$ fm are used in the continuum extrapolation.
  • The 2012 PDG reports a tension at the level of $3 \sigma$ between two exclusive determinations of $|V_{ub}|$. They are obtained by combining the experimental branching ratios of $B \to \tau \nu$ and $B \to \pi l \nu$ (respectively) with a theoretical computation of the hadronic matrix elements $\fB$ and the $B \to \pi$ form factor $f_+(q^2)$. To understand the tension, improved precision and a careful analysis of the systematics involved are necessary. We report the results of the ALPHA collaboration for $\fB$ from the lattice with 2 flavors of $O(a)$ improved Wilson fermions. We employ HQET, including $1/m_b$ corrections, with pion masses ranging down to $\approx$ 190 MeV. Renormalization and matching were performed non-perturbatively, and three lattice spacings reaching $a^{-1}\approx 4.1$ GeV are used in the continuum extrapolation. We also present progress towards a computation of $f_+(q^2)$, to directly compare two independent exclusive determinations of $|V_{ub}|$ with each other and with inclusive determinations. Additionally, we report on preliminary results for $\fBq{s}$, needed for the analysis of $B_s \to \mu^+\mu^-$.}
  • We present our analysis of B physics quantities using non-perturbatively matched Heavy Quark Effective Theory (HQET) in Nf= 2 lattice QCD on the CLS ensembles. Using all-to-all propagators, HYP-smeared static quarks, and the Generalized Eigenvalue Problem (GEVP) approach with a conservative plateau selection procedure, we are able to systematically control all sources of error. With significantly increased statistics compared to last year, our preliminary results are mb (mb) = 4.22(10)(4)z GeV for the MS b-quark mass, and fB = 193(9)stat (4)_\chi MeV and fBs = 219(12)stat MeV for the B-meson decay constants.
  • 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.
  • 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 analyze the systematic errors made when using the generalized eigenvalue problem to extract energies and matrix elements in lattice gauge theory. Effective theories such as HQET are also discussed. Numerical results are shown for the extraction of ground-state and excited B-meson masses and the ground-state decay constant in the static approximation.
  • We discuss some large effects of dynamical fermions. One is a cutoff effect, others concern the contribution of multi-pion states to correlation functions and are expected to survive the continuum limit. We then turn to the preparation for simulations at small lattice spacings which we are planning down to around a=0.04fm in order to understand the size of O(a^2)-effects of the standard O(a)-improved theory. The dependence of the lattice spacing on the bare coupling is determined through the Schr"odinger functional renormalized coupling.
  • We show that, in four dimensions, the quark propagator is affected by round-off errors for large values of the quark mass am and the time extent T/a even when double precision arithmetics is used to compute it. We introduce a definition of the solver residual which is sensitive to the problem and apply a Schwarz alternating procedure to compute the propagator in a number of sub-domains in the time direction. The effectiveness of the method is demonstrated in a numerical computation of the free one-dimensional Dirac propagator.
  • We introduce an alternative discretization for static quarks on the lattice retaining the O(a) improvement properties of the Eichten-Hill action. In this formulation, statistical fluctuations are reduced by a factor which grows exponentially with Euclidean time, x_0. For the first time, B-meson correlation functions are computed with good statistical precision in the static approximation for x_0>1 fm. At lattice spacings a \approx 0.1 fm, a \approx 0.08 fm and a \approx 0.07 fm the B_s-meson decay constant is determined in static and quenched approximations. A correction due to the finite mass of the b-quark is estimated by combining these static results with a recent determination of F_Ds.
  • We investigate the low-lying eigenvalues of the improved Wilson-Dirac operator in the Schroedinger functional with two dynamical quark flavors. At a lattice spacing of approximately 0.1 fm we find more very small eigenvalues than in the quenched case. These cause problems with HMC-type algorithms and in the evaluation of fermionic correlation functions. Through a simulation at a finer lattice spacing we are able to establish their nature as cutoff-effects
  • We address the locality problem arising in simulations, which take the square root of the staggered fermion determinant as a Boltzmann weight to reduce the number of dynamical quark tastes from four to two. We study analytically and numerically the square root of the staggered fermion operator as a candidate to define a two taste theory from first principles. Although it has the correct weight, this operator is non-local in the continuum limit. Our work serves as a warning that fundamental properties of field theories might be violated when employing blindly the square root trick. The question, whether a local operator reproducing the square root of the staggered fermion determinant exists, is left open.
  • We address the locality problem arising in simulations, which take the square root of the staggered fermion determinant as a Boltzmann weight to reduce the number of dynamical quark tastes. A definition of such a theory necessitates an underlying local fermion operator with the same determinant and the corresponding Green's functions to establish causality and unitarity. We illustrate this point by studying analytically and numerically the square root of the staggered fermion operator. Although it has the correct weight, this operator is non-local in the continuum limit. Our work serves as a warning that fundamental properties of field theories might be violated when employing blindly the square root trick. The question, whether a local operator reproducing the square root of the staggered fermion determinant exists, is left open.
  • We present and discuss results for cutoff effects in the PCAC masses and the mass dependence of r_0 for full QCD and various fermion actions. Our discussion of how one computes mass dependences - here of r_0 - is also relevant for comparisons with chiral perturbation theory.
  • We present two determinations of the coefficient $c_A$ in quenched QCD, needed to build the $O(a)$ improved axial current. The first condition used is the requirement that the PCAC quark mass, as a function of $x_0$, stays flat for a non-trivial spatial phase for the fermions in the Schr\"odinger functional. The second condition is that the PCAC relation for the ground-state and the first excited state at finite L give the same quark mass. Our results confirm previous findings that in the quenched theory the intrinsic $O(a)$ ambiguity of $c_A$ gets relevant around $\beta\simeq6.0$.
  • We present a computation of the decay constant f_Bs in quenched QCD. Our strategy is to combine new precise data from the static approximation with an interpolation of the decay constant around the charm quark mass region. This computation is the first step in demonstrating the feasability of a strategy for f_B in full QCD. The continuum limits in the static theory and at finite mass are taken separately and will be further improved.
  • We present an algorithmic study for the simulation of two massless flavors of O(a) improved Wilson quarks with Schroedinger functional boundary conditions. The algorithm used is Hybrid Monte Carlo with two pseudo-fermion fields as proposed by M. Hasenbusch. A gain in CPU cost of a factor two is reached when compared to one pseudo-fermion field due to the larger possible step-size. This study is integrated in the ALPHA project for the computation of the running of the renormalized quark mass. We include an update on these physics results.