• 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.
  • In this document we address an error discovered in the ensemble generation for our calculation of the $I=0$ $K\to\pi\pi$ amplitude (Phys. Rev. Lett. 115, 212001 (2015), arXiv:1505.07863) whereby the same random numbers were used for the two independent quark flavors, resulting in small but measurable correlations between gauge observables separated by 12 units in the y-direction. We conclude that the effects of this error are negligible compared to the overall errors on our calculation.
  • We report the first lattice QCD calculation of the complex kaon decay amplitude $A_0$ with physical kinematics, using a $32^3\times 64$ lattice volume and a single lattice spacing $a$, with $1/a= 1.3784(68)$ GeV. We find Re$(A_0) = 4.66(1.00)(1.26) \times 10^{-7}$ GeV and Im$(A_0) = -1.90(1.23)(1.08) \times 10^{-11}$ GeV, where the first error is statistical and the second systematic. The first value is in approximate agreement with the experimental result: Re$(A_0) = 3.3201(18) \times 10^{-7}$ GeV while the second can be used to compute the direct CP violating ratio Re$(\varepsilon'/\varepsilon)=1.38(5.15)(4.59)\times 10^{-4}$, which is $2.1\sigma$ below the experimental value $16.6(2.3)\times 10^{-4}$. The real part of $A_0$ is CP conserving and serves as a test of our method while the result for Re$(\varepsilon'/\varepsilon)$ provides a new test of the standard-model theory of CP violation, one which can be made more accurate with increasing computer capability.
  • We present new results for the amplitude $A_2$ for a kaon to decay into two pions with isospin $I=2$: Re$A_2 = 1.50(4)_\mathrm{stat}(14)_\mathrm{syst}\times 10^{-8}$ GeV; Im$A_2 = -6.99(20)_\mathrm{stat}(84)_\mathrm{syst}\times 10^{-13}$ GeV. These results were obtained from two ensembles generated at physical quark masses (in the isospin limit) with inverse lattice spacings $a^{-1}=1.728(4)$ GeV and $2.358(7)$ GeV. We are therefore able to perform a continuum extrapolation and hence largely to remove the dominant systematic uncertainty from our earlier results, that due to lattice artefacts. The only previous lattice computation of $K\to\pi\pi$ decays at physical kinematics was performed using an ensemble at a single, rather coarse, value of the lattice spacing ($a^{-1}\simeq 1.37(1)$ GeV). We confirm the observation that there is a significant cancellation between the two dominant contributions to Re$A_2$ which we suggest is an important ingredient in understanding the $\Delta I=1/2$ rule, Re$A_0$/Re$A_2\simeq 22.5$, where the subscript denotes the total isospin of the two-pion final state. Our result for $A_2$ implies that the electroweak penguin contribution to $\epsilon^\prime/\epsilon$ is Re($\epsilon^\prime/\epsilon)_\textrm{EWP}=-(6.6\pm 1.0)\times 10^{-4}$.
  • We present physical results for a variety of light hadronic quantities obtained via a combined analysis of three 2+1 flavour domain wall fermion ensemble sets. For two of our ensemble sets we used the Iwasaki gauge action with beta=2.13 (a^-1=1.75(4) GeV) and beta=2.25 (a^-1=2.31(4) GeV) and lattice sizes of 24^3 x 64 and 32^3 x 64 respectively, with unitary pion masses in the range 293(5)-417(10) MeV. The extent L_s for the 5^th dimension of the domain wall fermion formulation is L_s=16 in these ensembles. In this analysis we include a third ensemble set that makes use of the novel Iwasaki+DSDR (Dislocation Suppressing Determinant Ratio) gauge action at beta = 1.75 (a^-1=1.37(1) GeV) with a lattice size of 32^3 x 64 and L_s=32 to reach down to partially-quenched pion masses as low as 143(1) MeV and a unitary pion mass of 171(1) MeV, while retaining good chiral symmetry and topological tunneling. We demonstrate a significant improvement in our control over the chiral extrapolation, resulting in much improved continuum predictions for the above quantities. The main results of this analysis include the pion and kaon decay constants, f_\pi=127(3)_{stat}(3)_{sys} MeV and f_K = 152(3)_{stat}(2)_{sys} MeV respectively (f_K/f_\pi = 1.199(12)_{stat}(14)_{sys}); the average up/down quark mass and the strange-quark mass in the MSbar-scheme at 3 GeV, m_{ud}(MSbar, 3 GeV) = 3.05(8)_{stat}(6)_{sys} MeV and m_s(MSbar, 3 GeV) = 83.5(1.7)_{stat}(1.1)_{sys}; the neutral kaon mixing parameter in the MSbar-scheme at 3 GeV, B_K(MSbar,3 GeV) = 0.535(8)_{stat}(13)_{sys}, and in the RGI scheme, \hat B_K = 0.758(11)_{stat}(19)_{sys}; and the Sommer scales r_1 = 0.323(8)_{stat}(4)_{sys} fm and r_0 = 0.480(10)_{stat}(4)_{sys} (r_1/r_0 = 0.673(11)_{stat}(3)_{sys}). We also obtain values for the SU(2) ChPT effective couplings, \bar{l_3} = 2.91(23)_{stat}(7)_{sys}$ and \bar{l_4} = 3.99(16)_{stat}(9)_{sys}.
  • We develop and demonstrate techniques needed to compute the long distance contribution to the $K_{L}$-$K_{S}$ mass difference, $\Delta M_K$, in lattice QCD and carry out a first, exploratory calculation of this fundamental quantity. The calculation is performed on 2+1 flavor, domain wall fermion, $16^3\times32$ configurations with a 421 MeV pion mass and an inverse lattice spacing $1/a=1.73$ GeV. We include only current-current operators and drop all disconnected and double penguin diagrams. The short distance part of the mass difference in a 2+1 flavor calculation contains a quadratic divergence cut off by the lattice spacing. Here, this quadratic divergence is eliminated through the GIM mechanism by introducing a valence charm quark. The inclusion of the charm quark makes the complete calculation accessible to lattice methods provided the discretization errors associated with the charm quark can be controlled. The long distance effects are discussed for each parity channel separately. While we can see a clear signal in the parity odd channel, the signal to noise ratio in the parity even channel is exponentially decreasing as the separation between the two weak operators increases. We obtain a mass difference $\Delta M_K$ which ranges from $6.58(30)\times 10^{-12}$ MeV to $11.89(81)\times 10^{-12}$ MeV for kaon masses varying from 563 MeV to 839 MeV. Extensions of these methods are proposed which promise accurate results for both $\Delta M_K$ and $\epsilon_K$, including long distance effects.
  • We review the status of recent calculations by the RBC-UKQCD collaboration of the complex amplitude $A_2$, corresponding to the decay of a kaon to a two pion state with total isospin 2. In particular, we present preliminary results from two new ensembles: $48^3 \times 96$ with $a^{-1}=1.73$ GeV and $64^3 \times 128$ with $a^{-1}=2.3$ GeV, both at physical kinematics. Both ensembles were generated Iwasaki gauge action and domain wall fermion action with 2+1 flavours. These results, in comparison to our earlier ones on a $32^3$ DSDR lattice with $a^{-1}=1.36$ GeV, enable us to significantly reduce the discretization errors. The partial cancellation between the two dominant contractions contributing to Re($A_2$) has been confirmed and we believe that this cancellation is a major contribution to the $\Delta I=1/2$ rule.
  • We report on an ongoing calculation of hadronic matrix elements needed to parameterize K-Kbar mixing in generic BSM scenarios, using domain wall fermions (DWF) at two lattice spacings. Recent work by the SWME collaboration shows a significant disagreement with our previous results for two of these quantities. Since the origin of this disagreement is unknown, it is important to reduce the various uncertainties. In this work, we are using N_f=2+1 DWF with Iwasaki gauge action at inverse lattice spacings of 2.31 and 1.75 GeV, with multiple unitary pions on each ensemble, the lightest being 290 and 330 MeV on the finer and coarser of the two ensembles respectively. This extends previous work by the addition of a second lattice spacing (a^{-1}\approx 1.75 GeV). Renormalization is carried out non-perturbatively in the RI/MOM scheme and converted perturbatively to MSbar.
  • There has been much speculation as to the origin of the \Delta I = 1/2 rule (Re A_0/Re A_2 \simeq 22.5). We find that the two dominant contributions to the \Delta I=3/2, K \to \pi \pi{} correlation functions have opposite signs leading to a significant cancellation. This partial cancellation occurs in our computation of Re A_2 with physical quark masses and kinematics (where we reproduce the experimental value of A_2) and also for heavier pions at threshold. For Re A_0, although we do not have results at physical kinematics, we do have results for pions at zero-momentum with m_\pi{} \simeq 420 MeV (Re A_0/Re A_2=9.1(2.1)) and m_\pi{} \simeq 330 MeV (Re A_0/Re A_2=12.0(1.7)). The contributions which partially cancel in Re A_2 are also the largest ones in Re A_0, but now they have the same sign and so enhance this amplitude. The emerging explanation of the \Delta I=1/2 rule is a combination of the perturbative running to scales of O(2 GeV), a relative suppression of Re A_2 through the cancellation of the two dominant contributions and the corresponding enhancement of Re A_0. QCD and EWP penguin operators make only very small contributions at such scales.
  • We report on the first realistic \emph{ab initio} calculation of a hadronic weak decay, that of the amplitude $A_2$ for a kaon to decay into two \pi-mesons with isospin 2. We find Re$A_2=(1.436\pm 0.063_{\textrm{stat}}\pm 0.258_{\textrm{syst}})\,10^{-8}\,\textrm{GeV}$ in good agreement with the experimental result and for the hitherto unknown imaginary part we find {Im}$\,A_2=-(6.83 \pm 0.51_{\textrm{stat}} \pm 1.30_{\textrm{syst}})\,10^{-13}\,{\rm GeV}$. Moreover combining our result for Im\,$A_2$ with experimental values of Re\,$A_2$, Re\,$A_0$ and $\epsilon^\prime/\epsilon$, we obtain the following value for the unknown ratio Im\,$A_0$/Re\,$A_0$ within the Standard Model: $\mathrm{Im}\,A_0/\mathrm{Re}\,A_0=-1.63(19)_{\mathrm{stat}}(20)_{\mathrm{syst}}\times10^{-4}$. One consequence of these results is that the contribution from Im\,$A_2$ to the direct CP violation parameter $\epsilon^{\prime}$ (the so-called Electroweak Penguin, EWP, contribution) is Re$(\epsilon^\prime/\epsilon)_{\mathrm{EWP}} = -(6.52 \pm 0.49_{\textrm{stat}} \pm 1.24_{\textrm{syst}}) \times 10^{-4}$. We explain why this calculation of $A_2$ represents a major milestone for lattice QCD and discuss the exciting prospects for a full quantitative understanding of CP-violation in kaon decays.
  • 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 use one-loop $\SU(2)_L\times \SU(2)_R$ chiral perturbation theory ($\SU(2)$ ChPT) to study the behaviour of the form-factors for semileptonic $K\to\pi$ decays with the pion mass at $q^2=0$ and at $q^2_{\textrm{max}}=(m_K-m_\pi)^2$, where $q$ is the momentum transfer. At $q^2=0$, the final-state pion has an energy of approximately $m_K/2$ (for $m_K\gg m_\pi$) and so is not soft, nevertheless it is possible to compute the chiral logarithms, i.e. the corrections of $O(m_\pi^2\log(m_\pi^2))$. We envisage that our results at $q^2=0$ will be useful in extrapolating lattice QCD results to physical masses. A consequence of the Callan-Treiman relation is that in the $\SU(2)$ chiral limit ($m_u=m_d=0$), the scalar form factor $f^0$ at $\qsqmax$ is equal to $f^{(K)}/f$, the ratio of the kaon and pion leptonic decay constants in the chiral limit. Lattice results for the scalar form factor at $\qsqmax$ are obtained with excellent precision, but at the masses at which the simulations are performed the results are about 25% below $f^{(K)}/f$ and are increasing only very slowly. We investigate the chiral behaviour of $f^0(\qsqmax)$ and find large corrections which provide a semi-quantitative explanation of the difference between the lattice results and $f^{(K)}/f$. We stress the generality of the relation $f^0_{P\to\pi}(\qsqmax)=f^{(P)}/f$ in the $\SU(2)$ chiral limit, where $P=K,D$ or $B$ and briefly comment on the potential value of using this theorem in obtaining physical results from lattice simulations.
  • By imposing twisted boundary conditions on quark fields it is possible to access components of momenta other than integer multiples of 2pi/L on a lattice with spatial volume L^3. We use Chiral Perturbation Theory to study finite-volume effects with twisted boundary conditions for quantities without final-state interactions, such as meson masses, decay constants and semileptonic form factors, and confirm that they remain exponentially small with the volume. We show that this is also the case for "partially twisted" boundary conditions, in which (some of) the valence quarks satisfy twisted boundary conditions but the sea quarks satisfy periodic boundary conditions. This observation implies that it is not necessary to generate new gluon configurations for every choice of the twist angle, making the method much more practicable. For K->pipi decays we show that the breaking of isospin symmetry by the twisted boundary conditions implies that the amplitudes cannot be determined in general (on this point we disagree with a recent claim).
  • We compute the matrix elements of the operators which contribute to spectator effects in inclusive decays of $B$-mesons. The results agree well with estimates based on the vacuum saturation (factorization) hypothesis. For the ratio of lifetimes of charged and neutral mesons we find $\tau(B^-)/\tau(B_d)=1.03\pm 0.02\pm 0.03$, where the first error represents the uncertainty in our evaluation of the matrix elements, and the second is an estimate of the uncertainty due to the fact that the Wilson coefficient functions have only been evaluated at tree-level in perturbation theory. This result is in agreement with the experimental measurement. We also discuss the implications of our results for the semileptonic branching ratio and the charm yield.
  • We review the application of lattice QCD to the phenomenology of b- and c-quarks. After a short discussion of the lattice techniques used to evaluate hadronic matrix elements and the corresponding systematic uncertainties, we summarise results for leptonic decay constants, B--Bbar mixing, semileptonic and rare radiative decays. A discussion of the determination of heavy quark effective theory parameters is followed by an explanation of the difficulty in applying lattice methods to exclusive nonleptonic decays.
  • We present a high statistics lattice calculation of the B--meson binding energy $\overline{\Lambda}$ of the heavy--quark inside the pseudoscalar B--meson. Our numerical results have been obtained from several independent numerical simulations at $\beta=6.0$, $6.2$ and $6.4$, and using, for the meson correlators, the results obtained by the APE group at the same values of $\beta$. Our best estimate, obtained by combining results at different values of $\beta$, is $\overline{\Lambda}=180^{+30}_{-20}$ MeV. For the $\overline{MS}$ running mass, we obtain $\overline{m}_{b}(\overline{m}_{b})=4.15 \pm 0.05 \pm 0.20$ GeV, in reasonable agreement with previous determinations. The systematic error is the truncation of the perturbative series in the matching condition of the relevant operator of the Heavy Quark Effective Theory.
  • We present the first lattice calculation of the B-meson binding energy $\labar$ and of the kinetic energy $-\lambda_1/2 m_Q$ of the heavy-quark inside the pseudoscalar B-meson. This calculation has required the non-perturbative subtraction of the power divergences present in matrix elements of the Lagrangian operator $\bar h D_4 h$ and of the kinetic energy operator $\bar h \vec D^2 h$. The non-perturbative renormalisation of the relevant operators has been implemented by imposing suitable renormalisation conditions on quark matrix elements, in the Landau gauge. Our numerical results have been obtained from several independent numerical simulations at $\beta=6.0$ and $6.2$, and using, for the meson correlators, the results obtained by the APE group at the same values of $\beta$. Our best estimate, obtained by combining results at different values of $\beta$, is $\labar =190 \err{50}{30}$ MeV. For the $\overline{MS}$ running mass, we obtain $\overline {m}_b(\overline {m}_b) =4.17 \pm 0.06$ GeV, in reasonable agreement with previous determinations. From a subset of 36 configurations, we were only able to establish a loose upper bound on the $b$-quark kinetic energy in a $B$-meson, $\lambda_1=\langle B \vert \bar h \vec{D}^{2} h \vert B \rangle /(2 M_B )<$~1\, GeV$^2$. This shows that a much larger statistical sample is needed to determine this important parameter.