
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
electroquenched 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 masssplittings.

We review lattice results related to pion, kaon, D and Bmeson physics with
the aim of making them easily accessible to the particle physics community.
More specifically, we report on the determination of the lightquark 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 lowenergy 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 heavyquark sector, we provide results for mc and mb (also new
compared to the previous review), as well as those for D and Bmeson decay
constants, form factors, and mixing parameters. These are the heavyquark
quantities most relevant for the determination of CKM matrix elements and the
global CKM unitaritytriangle 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
ydirection. 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
standardmodel 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 twopion 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
partiallyquenched 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 strangequark mass in the MSbarscheme 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
MSbarscheme 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 currentcurrent 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 RBCUKQCD 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 KKbar 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 nonperturbatively 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 zeromomentum 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 \pimesons
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 socalled 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
CPviolation in kaon decays.

We present results for the first two moments of the lightcone distribution
amplitudes of the pion and kaon pseudoscalar 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
nonperturbative renormalisation which we perform using the RomeSouthampton
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 oneloop $\SU(2)_L\times \SU(2)_R$ chiral perturbation theory
($\SU(2)$ ChPT) to study the behaviour of the formfactors for semileptonic
$K\to\pi$ decays with the pion mass at $q^2=0$ and at
$q^2_{\textrm{max}}=(m_Km_\pi)^2$, where $q$ is the momentum transfer. At
$q^2=0$, the finalstate 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 CallanTreiman 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 semiquantitative
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
finitevolume effects with twisted boundary conditions for quantities without
finalstate 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 treelevel 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
cquarks. 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, BBbar 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 Bmeson binding
energy $\overline{\Lambda}$ of the heavyquark inside the pseudoscalar
Bmeson. 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 Bmeson binding energy
$\labar$ and of the kinetic energy $\lambda_1/2 m_Q$ of the heavyquark inside
the pseudoscalar Bmeson. This calculation has required the nonperturbative
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 nonperturbative 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.