• We present a first-principles lattice QCD+QED calculation at physical pion mass of the leading-order hadronic vacuum polarization contribution to the muon anomalous magnetic moment. The total contribution of up, down, strange, and charm quarks including QED and strong isospin breaking effects is found to be $a_\mu^{\rm HVP~LO}=715.4(16.3)(9.2) \times 10^{-10}$, where the first error is statistical and the second is systematic. By supplementing lattice data for very short and long distances with experimental R-ratio data using the compilation of Ref. [1], we significantly improve the precision of our calculation and find $a_\mu^{\rm HVP~LO} = 692.5(1.4)(0.5)(0.7)(2.1) \times 10^{-10}$ with lattice statistical, lattice systematic, R-ratio statistical, and R-ratio systematic errors given separately. This is the currently most precise determination of the leading-order hadronic vacuum polarization contribution to the muon anomalous magnetic moment. In addition, we present the first lattice calculation of the light-quark QED correction at physical pion mass.
  • We review the domain wall charm physics program of the RBC and UKQCD collaborations based on simulations including ensembles with physical pion mass. We summarise our current set-up and present a status update on the decay constants $f_D$, $f_{D_s}$, the charm quark mass, heavy-light and heavy-strange bag parameters and the ratio $\xi$.
  • We calculate the strong isospin breaking and QED corrections to meson masses and the hadronic vacuum polarization in an exploratory study on a $64\times24^3$ lattice with an inverse lattice spacing of $a^{-1}=1.78$ GeV and an isospin symmetric pion mass of $m_\pi=340$ MeV. We include QED in an electro-quenched setup using two different methods, a stochastic and a perturbative approach. We find that the electromagnetic correction to the leading hadronic contribution to the anomalous magnetic moment of the muon is smaller than $1\%$ for the up quark and $0.1\%$ for the strange quark, although it should be noted that this is obtained using unphysical light quark masses. In addition to the results themselves, we compare the precision which can be reached for the same computational cost using each method. Such a comparison is also made for the meson electromagnetic mass-splittings.
  • 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 leading hadronic contribution to the muon anomalous magnetic moment due to strange quark-connected vacuum polarisation effects. Simulations were performed using RBC--UKQCD's $N_f=2+1$ domain wall fermion ensembles with physical light sea quark masses at two lattice spacings. We consider a large number of analysis scenarios in order to obtain solid estimates for residual systematic effects. Our final result in the continuum limit is $a_\mu^{(2)\,{\rm had},\,s}=53.1(9)\left(^{+1}_{-3}\right)\times10^{-10}$.
  • We report the first lattice QCD calculation of the hadronic vacuum polarization disconnected contribution to the muon anomalous magnetic moment at physical pion mass. The calculation uses a refined noise-reduction technique which enabled the control of statistical uncertainties at the desired level with modest computational effort. Measurements were performed on the $48^3 \times 96$ physical-pion-mass lattice generated by the RBC and UKQCD collaborations. We find $a_\mu^{\rm HVP~(LO)~DISC} = -9.6(3.3)(2.3)\times 10^{-10}$, where the first error is statistical and the second systematic.
  • We have performed fits of the pseudoscalar masses and decay constants, from a variety of RBC-UKQCD domain wall fermion ensembles, to $SU(2)$ partially quenched chiral perturbation theory at next-to leading order (NLO) and next-to-next-to leading order (NNLO). We report values for 9 NLO and 8 linearly independent combinations of NNLO partially quenched low energy constants, which we compare to other lattice and phenomenological determinations. We discuss the size of successive terms in the chiral expansion and use our large set of low energy constants to make predictions for mass splittings due to QCD isospin breaking effects and the S-wave $\pi \pi$ scattering lengths. We conclude that, for the range of pseudoscalar masses explored in this work, $115~\mathrm{MeV} \lesssim m_{\rm PS} \lesssim 430~\mathrm{MeV}$, the NNLO $SU(2)$ expansion is quite robust and can fit lattice data with percent-scale accuracy.
  • 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.
  • As part of the UKQCD and RBC collaborations' N_f=2+1 domain-wall fermion phenomenology programme, we calculate the first two moments of the light-cone distribution amplitudes of the pseudoscalar mesons pion and kaon and the (longitudinally-polarised) vector mesons rho, K-star and phi. We obtain the desired quantities with good precision and are able to discern the expected quark-mass dependence of SU(3)-flavour breaking effects. An important ingredient of the calculation is the nonperturbative renormalisation of lattice operators using the RI'/MOM technique.
  • We calculated the first two moments of the light-cone distribution amplitudes for the pseudoscalar mesons ($\pi$ and $K$) and the longitudinally polarised vector mesons ($\rho$, $K^*$ and $\phi$) as part of the UKQCD and RBC collaborations' $N_f=2+1$ domain-wall fermion phenomenology programme. These quantities were obtained with a good precision and, in particular, the expected effects of $SU(3)$-flavour symmetry breaking were observed. Operators were renormalised non-perturbatively and extrapolations to the physical point were made, guided by leading order chiral perturbation theory. The main results presented are for two volumes, $16^3\times 32$ and $24^3\times 64$, with a common lattice spacing. Preliminary results for a lattice with a finer lattice spacing, $32^3\times64$, are discussed and a first look is taken at the use of twisted boundary conditions to extract distribution amplitudes.
  • 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.
  • Using partially twisted boundary conditions we compute the K->pi semi-leptonic form factors in the range of momentum transfers 0 <~ q^2 <= q^2_{max}=(mK-mpi)^2 in lattice QCD with N_f=2+1 dynamical flavours. In this way we are able to determine f+(0) without any interpolation in the momentum transfer, thus eliminating one source of systematic error. This study confirms our earlier phenomenological ansatz for the strange quark mass dependence of the scalar form factor. We identify and estimate potentially significant NNLO effects in the chiral expansion that guides the extrapolation of the data to the physical point. Our main result is f+(0) = 0.9599(34)(^{+31}_{-43})(14)$, where the first error is statistical, the second error is due to the uncertainties in the chiral extrapolation of the lattice data and the last error is an estimate of potential discretisation effects.
  • M. Antonelli, D.M. Asner, D. Bauer, T. Becher, M. Beneke, A.J. Bevan, M. Blanke, C. Bloise, M. Bona, A. Bondar, C. Bozzi, J. Brod, A.J. Buras, N. Cabibbo, A. Carbone, G. Cavoto, V. Cirigliano, M. Ciuchini, J.P. Coleman, D.P. Cronin-Hennessy, J.P. Dalseno, C.H. Davies, F. DiLodovico, J. Dingfelder, Z. Dolezal, S. Donati, W. Dungel, U. Egede, G. Eigen, R. Faccini, T. Feldmann, F. Ferroni, J.M. Flynn, E. Franco, M. Fujikawa, I.K. Furic, P. Gambino, E. Gardi, T.J. Gershon, S. Giagu, E. Golowich, T. Goto, C. Greub, C. Grojean, D. Guadagnoli, U. A. Haisch, R.F. Harr, A.H. Hoang, T. Hurth, G. Isidori, D.E. Jaffe, A. Jüttner, S. Jäger, A. Khodjamirian, P. Koppenburg, R. V. Kowalewski, P. Krokovny, A. S. Kronfeld, J. Laiho, G. Lanfranchi, T. E. Latham, J. Libby, A. Limosani, D. Lopes Pegna, C. D. Lu, V. Lubicz, E. Lunghi, V. G. Lüth, K. Maltman, W. J. Marciano, E. C. Martin, G. Martinelli, F. Martinez-Vidal, A. Masiero, V. Mateu, F. Mescia, G. Mohanty, M. Moulson, M. Neubert, H. Neufeld, S. Nishida, N. Offen, M. Palutan, P. Paradisi, Z. Parsa, E. Passemar, M. Patel, B. D. Pecjak, A. A. Petrov, A. Pich, M. Pierini, B. Plaster, A. Powell, S. Prell, J. Rademaker, M. Rescigno, S. Ricciardi, P. Robbe, E. Rodrigues, M. Rotondo, R. Sacco, C. J. Schilling, O. Schneider, E. E. Scholz, B. A. Schumm, C. Schwanda, A. J. Schwartz, B. Sciascia, J. Serrano, J. Shigemitsu, I. J. Shipsey, A. Sibidanov, L. Silvestrini, F. Simonetto, S. Simula, C. Smith, A. Soni, L. Sonnenschein, V. Sordini, M. Sozzi, T. Spadaro, P. Spradlin, A. Stocchi, N. Tantalo, C. Tarantino, A. V. Telnov, D. Tonelli, I. S. Towner, K. Trabelsi, P. Urquijo, R. S. Van de Water, R. J. Van Kooten, J. Virto, G. Volpi, R. Wanke, S. Westhoff, G. Wilkinson, M. Wingate, Y. Xie, J. Zupan
    Feb. 19, 2010 hep-ph
    One of the major challenges of particle physics has been to gain an in-depth understanding of the role of quark flavor and measurements and theoretical interpretations of their results have advanced tremendously: apart from masses and quantum numbers of flavor particles, there now exist detailed measurements of the characteristics of their interactions allowing stringent tests of Standard Model predictions. Among the most interesting phenomena of flavor physics is the violation of the CP symmetry that has been subtle and difficult to explore. Till early 1990s observations of CP violation were confined to neutral $K$ mesons, but since then a large number of CP-violating processes have been studied in detail in neutral $B$ mesons. In parallel, measurements of the couplings of the heavy quarks and the dynamics for their decays in large samples of $K, D$, and $B$ mesons have been greatly improved in accuracy and the results are being used as probes in the search for deviations from the Standard Model. In the near future, there will be a transition from the current to a new generation of experiments, thus a review of the status of quark flavor physics is timely. This report summarizes the results of the current generation of experiments that is about to be completed and it confronts these results with the theoretical understanding of the field.
  • 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.
  • In this paper we investigate the benefits of using Z(2)xZ(2) single timeslice stochastic sources for the calculation of light quark physics on the lattice. Meson 2-point correlators measured using sources stochastic in only spin and those stochastic in both spin and colour indices are compared to point source correlators on the unit gauge and on a 16^3 x 32 Domain Wall QCD ensemble. It is found that the use of stochastic sources gives a considerable improvement in statistics for the same computational cost. The neutral kaon mixing matrix element B_K is also calculated on this ensemble with stochastic sources, but we conclude that the stochastic method offers no significant advantage over the traditional gauge-fixed wall source approach which already offers an exact volume average. We also discuss the application to semileptonic form factors in conjunction with partially twisted boundary conditions.
  • We present the first results for the Kl3 form factor from simulations with 2+1 flavours of dynamical domain wall quarks. Combining our result, namely f_+(0)=0.964(5), with the latest experimental results for Kl3 decays leads to |V_{us}|=0.2249(14), reducing the uncertaintity in this important parameter. For the O(p^6) term in the chiral expansion we obtain \Delta f=-0.013(5).
  • We present results for the first two moments of the light-cone distribution amplitudes of the pion and kaon pseudo-scalar mesons and of the rho, K* and phi vector mesons. The calculations are performed on the RBC/UKQCD collaborations' ensembles generated with the Iwasaki gauge action and with 2+1 flavours of domain wall fermions. In addition we also provide some results on the necessary non-perturbative renormalisation which we perform using the Rome-Southampton method. We discuss the benefits of the momentum source approach such as much smaller statistical errors and the possibility to see effects of the discretisation.
  • We compute the electromagnetic form factor of a "pion" with mass m_pi=330MeV at low values of Q^2\equiv -q^2, where q is the momentum transfer. The computations are performed in a lattice simulation using an ensemble of the RBC/UKQCD collaboration's gauge configurations with Domain Wall Fermions and the Iwasaki gauge action with an inverse lattice spacing of 1.73(3)GeV. In order to be able to reach low momentum transfers we use partially twisted boundary conditions using the techniques we have developed and tested earlier. For the pion of mass 330MeV we find a charge radius given by <r_pi^2>_{330MeV}=0.354(31)fm^2 which, using NLO SU(2) chiral perturbation theory, extrapolates to a value of <r_pi^2>=0.418(31)fm^2 for a physical pion, in agreement with the experimentally determined result. We confirm that there is a significant reduction in computational cost when using propagators computed from a single time-slice stochastic source compared to using those with a point source; for m_pi=330MeV and volume (2.74fm)^3 we find the reduction is approximately a factor of 12.
  • We have simulated QCD using 2+1 flavors of domain wall quarks on a $(2.74 {\rm fm})^3$ volume with an inverse lattice scale of $a^{-1} = 1.729(28)$ GeV. The up and down (light) quarks are degenerate in our calculations and we have used four values for the ratio of light quark masses to the strange (heavy) quark mass in our simulations: 0.217, 0.350, 0.617 and 0.884. We have measured pseudoscalar meson masses and decay constants, the kaon bag parameter $B_K$ and vector meson couplings. We have used SU(2) chiral perturbation theory, which assumes only the up and down quark masses are small, and SU(3) chiral perturbation theory to extrapolate to the physical values for the light quark masses. While next-to-leading order formulae from both approaches fit our data for light quarks, we find the higher order corrections for SU(3) very large, making such fits unreliable. We also find that SU(3) does not fit our data when the quark masses are near the physical strange quark mass. Thus, we rely on SU(2) chiral perturbation theory for accurate results. We use the masses of the $\Omega$ baryon, and the $\pi$ and $K$ mesons to set the lattice scale and determine the quark masses. We then find $f_\pi = 124.1(3.6)_{\rm stat}(6.9)_{\rm syst} {\rm MeV}$, $f_K = 149.6(3.6)_{\rm stat}(6.3)_{\rm syst} {\rm MeV}$ and $f_K/f_\pi = 1.205(0.018)_{\rm stat}(0.062)_{\rm syst}$. Using non-perturbative renormalization to relate lattice regularized quark masses to RI-MOM masses, and perturbation theory to relate these to $\bar{\rm MS}$ we find $ m_{ud}^{\bar{\rm MS}}(2 {\rm GeV}) = 3.72(0.16)_{\rm stat}(0.33)_{\rm ren}(0.18)_{\rm syst} {\rm MeV}$ and $m_{s}^{\bar{\rm MS}}(2 {\rm GeV}) = 107.3(4.4)_{\rm stat}(9.7)_{\rm ren}(4.9)_{\rm syst} {\rm MeV}$.
  • We present a calculation of the renormalization coefficients of the quark bilinear operators and the K-Kbar mixing parameter B_K. The coefficients relating the bare lattice operators to those in the RI/MOM scheme are computed non-perturbatively and then matched perturbatively to the MSbar scheme. The coefficients are calculated on the RBC/UKQCD 2+1 flavor dynamical lattice configurations. Specifically we use a 16^3 x 32 lattice volume, the Iwasaki gauge action at beta=2.13 and domain wall fermions with L_s=16.
  • We present results for light meson masses and pseudoscalar decay constants from the first of a series of lattice calculations with 2+1 dynamical flavors of domain wall fermions and the Iwasaki gauge action. The work reported here was done at a fixed lattice spacing of about 0.12 fm on a 16^3\times32 lattice, which amounts to a spatial volume of (2 fm)^3 in physical units. The number of sites in the fifth dimension is 16, which gives m_{res} = 0.00308(4) in these simulations. Three values of input light sea quark masses, m_l^{sea} \approx 0.85 m_s, 0.59 m_s and 0.33 m_s were used to allow for extrapolations to the physical light quark limit, whilst the heavier sea quark mass was fixed to approximately the physical strange quark mass m_s. The exact rational hybrid Monte Carlo algorithm was used to evaluate the fractional powers of the fermion determinants in the ensemble generation. We have found that f_\pi = 127(4) MeV, f_K = 157(5) MeV and f_K/f_\pi = 1.24(2), where the errors are statistical only, which are in good agreement with the experimental values.
  • We present results for the couplings of light vector mesons to vector and tensor currents and on the low moments of the light-cone distribution amplitudes of the pion and kaon. The calculations are performed on the RBC and UKQCD collaborations' ensembles generated with the Iwasaki gauge action and with 2+1 flavours of domain wall fermions. The (preliminary) results for the ratios of the couplings of the vector meson to the vector and tensor currents (f_V and f^T_V respectively) in the MSbar scheme at 2GeV are: f_\rho^T/f_\rho=0.681(20); f_{K^\ast}^T/f_{K^\ast}=0.712(11) and f_\phi^T/f_\phi=0.751(9). For the first moment of the kaon's distribution amplitude we find (in the same scheme and at the same scale) <xi>_K=0.029(2) and for the second moment <xi^2>_\pi=0.28(3) and <xi^2>_K=0.27(2).
  • The introduction of partially twisted boundary conditions allows weak and electromagnetic form factors to be evaluated at specified values of the hadronic momenta (and hence momentum transfers) in lattice simulations. We present and demonstrate this technique for the computation of the K->pi semileptonic form factor at zero momentum transfer and for the electromagnetic form factor of the pion at arbitrarily small momentum transfers. These exploratory computations are carried out in full QCD with 3 flavours of sea quarks, but with only two values of m_u=m_d which limits our ability to perform the chiral extrapolations. The results should therefore be viewed primarily as a demonstration of the feasibility of the method. For the K-> pi form factor we compare the new technique to the conventional approach and for the pion form factor we assess our results for very small momentum transfer with the help of chiral perturbation theory.