
We present a firstprinciples lattice QCD+QED calculation at physical pion
mass of the leadingorder 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 Rratio 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, Rratio statistical,
and Rratio systematic errors given separately. This is the currently most
precise determination of the leadingorder hadronic vacuum polarization
contribution to the muon anomalous magnetic moment. In addition, we present the
first lattice calculation of the lightquark 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 setup and present a status update on the decay
constants $f_D$, $f_{D_s}$, the charm quark mass, heavylight and heavystrange
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
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.

We present results for the leading hadronic contribution to the muon
anomalous magnetic moment due to strange quarkconnected vacuum polarisation
effects. Simulations were performed using RBCUKQCD'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 noisereduction 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$ physicalpionmass 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 RBCUKQCD domain wall fermion ensembles, to $SU(2)$ partially
quenched chiral perturbation theory at nextto leading order (NLO) and
nexttonextto 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 Swave $\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 percentscale 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
threepoint 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 nonperturbatively O(a)
improved Wilson fermions for pion masses in the range of 450  760 MeV. Results
for pseudoscalar and vector meson twopoint 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 lowestlying 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 domainwall fermion
phenomenology programme, we calculate the first two moments of the lightcone
distribution amplitudes of the pseudoscalar mesons pion and kaon and the
(longitudinallypolarised) vector mesons rho, Kstar and phi. We obtain the
desired quantities with good precision and are able to discern the expected
quarkmass 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 lightcone 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$ domainwall 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 nonperturbatively 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 ongoing 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
semileptonic form factors in the range of momentum transfers 0 <~ q^2 <=
q^2_{max}=(mKmpi)^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.

One of the major challenges of particle physics has been to gain an indepth
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 CPviolating 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 ongoing project to compute mesonic and baryonic two and
threepoint 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 DDHMC 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 2point 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 gaugefixed 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 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 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 timeslice 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 nexttoleading 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
nonperturbative renormalization to relate lattice regularized quark masses to
RIMOM 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 KKbar mixing parameter B_K. The coefficients
relating the bare lattice operators to those in the RI/MOM scheme are computed
nonperturbatively 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 lightcone 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.