
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 compute the renormalisation factors (Zmatrices) of the $\Delta F=2$
fourquark operators needed for Beyond the Standard Model (BSM) kaon mixing. We
work with nf=2+1 flavours of DomainWall fermions whose chiralflavour
properties are essential to maintain a continuumlike mixing pattern. We
introduce new RISMOM renormalisation schemes, which we argue are better
behaved compared to the commonlyused corresponding RIMOM one. We find that,
once converted to MS, the Zfactors computed through these RISMOM schemes are
in good agreement but differ significantly from the ones computed through the
RIMOM scheme. The RISMOM Zfactors presented here have been used to compute
the BSM neutral kaon mixing matrix elements in the companion paper [1]. We
argue that the renormalisation procedure is responsible for the discrepancies
observed by different collaborations, we will investigate and elucidate the
origin of these differences throughout this work.

We present results for several light hadronic quantities ($f_\pi$, $f_K$,
$B_K$, $m_{ud}$, $m_s$, $t_0^{1/2}$, $w_0$) obtained from simulations of 2+1
flavor domain wall lattice QCD with large physical volumes and nearlyphysical
pion masses at two lattice spacings. We perform a short, O(3)%, extrapolation
in pion mass to the physical values by combining our new data in a simultaneous
chiral/continuum `global fit' with a number of other ensembles with heavier
pion masses. We use the physical values of $m_\pi$, $m_K$ and $m_\Omega$ to
determine the two quark masses and the scale  all other quantities are outputs
from our simulations. We obtain results with subpercent statistical errors and
negligible chiral and finitevolume systematics for these light hadronic
quantities, including: $f_\pi$ = 130.2(9) MeV; $f_K$ = 155.5(8) MeV; the
average up/down quark mass and strange quark mass in the $\bar {\rm MS}$ scheme
at 3 GeV, 2.997(49) and 81.64(1.17) MeV respectively; and the neutral kaon
mixing parameter, $B_K$, in the RGI scheme, 0.750(15) and the $\bar{\rm MS}$
scheme at 3 GeV, 0.530(11).

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 derive the exactly conserved vector, and almost conserved axial currents
for rational approximations to the overlap operator with a general Mobius
kernel. The approach maintains manifest Hermiticity, and allows matrix elements
of the currents to be constructed at no extra cost after solution of the usual
5d system of equations, similar to the original approach of Furman and Shamir
for domain wall Fermions.

The renormalized nexttoleadingorder (NLO) chiral lowenergy constant,
$L_{10}^r$, is determined in a complete nexttonexttoleadingorder (NNLO)
analysis, using a combination of lattice and continuum data for the flavor $ud$
$VA$ correlator and results from a recent chiral sumrule analysis of the
flavorbreaking combination of $ud$ and $us$ $VA$ correlator differences. The
analysis also fixes two combinations of NNLO lowenergy constants, the
determination of which is crucial to the precision achieved for $L_{10}^r$.
Using the results of the flavorbreaking chiral $VA$ sum rule obtained with
current versions of the strange hadronic $\tau$ branching fractions as input,
we find $L_{10}^r(m_\rho )\, =\, 0.00346(32)$. This result represents the
first NNLO determination of $L_{10}^r$ having all inputs under full theoretical
and/or experimental control, and the best current precision for this quantity.

Recent analyses of flavorbreaking hadronic$\tau$decaybased sum rules
produce values of $\vert V_{us}\vert$ $\sim 3\sigma$ low compared to 3family
unitarity expectations. An unresolved systematic issue is the significant
variation in $\vert V_{us}\vert$ produced by different prescriptions for
treating the slowly converging $D=2$ OPE series. We investigate the reliability
of these prescriptions using lattice data for various flavorbreaking
correlators and show the fixedscale prescription is clearly preferred.
Preliminary updates of the conventional $\tau$based, and related mixed
$\tau$electroproductiondatabased, sum rule analyses incorporating Bfactory
results for lowmultiplicity strange $\tau$ decay mode distributions are then
performed. Use of the preferred FOPT $D=2$ OPE prescription is shown to
significantly reduce the discrepancy between 3family unitarity expectations
and the sum rule results.

A combination of lattice and continuum data for the lightquark VA
correlator, supplemented by results from a chiral sumrule analysis of the
flavorbreaking flavor $ud$$us$ VA correlator difference, is shown to make
possible a highprecision NNLO determination of the renormalized NLO chiral
lowenergy constant $L_{10}^r$. Key to this determination is the ability to
simultaneously fix the two combinations of NNLO lowenergy constants also
entering the analysis. With current versions of the strange hadronic $\tau$
branching fractions required as input to the flavorbreaking VA sum rule, we
find $L_{10}^r(m_\rho ) = 0.00346(29)$. This represents both the best current
precision for $L_{10}^r$, and the first NNLO determination having all errors
under full control.

Recent sum rule determinations of V_us, employing flavorbreaking
combinations of hadronic tau decay data, are significantly lower than either
expectations based on 3family unitarity or determinations from K_ell3 and
Gamma[K_mu2]/Gamma[pi_mu2]. We use lattice data to investigate the
accuracy/reliability of the OPE representation of the flavorbreaking
correlator combination entering the tau decay analyses. The behavior of an
alternate correlator combination, constructed to reduce problems associated
with the slow convergence of the D = 2 OPE series, and entering an alternate
sum rule requiring both electroproduction crosssection and hadronic tau decay
data, is also investigated. Preliminary updates of both analyses, with the
lessons learned from the lattice data in mind, are also presented.

We present preliminary results on extractions of the chiral LECs L_10 and
C_87 and constraints on the excited pseudoscalar state pi(1300) and pi(1800)
decay constants obtained from an analysis of lattice data for the flavor ud
light quark VA correlator. A comparison of the results for the correlator to
the corresponding mildlymodeldependent continuum results (based primarily on
experimental hadronic tau decay data) is also given

We describe the computation of the amplitude A_2 for a kaon to decay into two
pions with isospin I=2. The results presented in the letter Phys.Rev.Lett. 108
(2012) 141601 from an analysis of 63 gluon configurations are updated to 146
configurations giving Re$A_2=1.381(46)_{\textrm{stat}}(258)_{\textrm{syst}}
10^{8}$ GeV and Im$A_2=6.54(46)_{\textrm{stat}}(120)_{\textrm{syst}}10^{13}$
GeV. Re$A_2$ is in good agreement with the experimental result, whereas the
value of Im$A_2$ was hitherto unknown. We are also working towards a direct
computation of the $K\to(\pi\pi)_{I=0}$ amplitude $A_0$ but, within the
standard model, our result for Im$A_2$ can be combined with the experimental
results for Re$A_0$, Re$A_2$ and $\epsilon^\prime/\epsilon$ to give
Im$A_0/$Re$A_0= 1.61(28)\times 10^{4}$ . Our result for Im\,$A_2$ implies
that the electroweak penguin (EWP) contribution to $\epsilon^\prime/\epsilon$
is Re$(\epsilon^\prime/\epsilon)_{\mathrm{EWP}} = (6.25 \pm
0.44_{\textrm{stat}} \pm 1.19_{\textrm{syst}}) \times 10^{4}$.

We show that the running of operators which mix under renormalization can be
computed fully nonperturbatively as a product of continuum step scaling
matrices. These step scaling matrices are obtained by taking the "ratio" of Z
matrices computed at different energies in an RIMOM type scheme for which
twisted boundary conditions are an essential ingredient. Our method allows us
to relax the bounds of the RomeSouthampton window. We also explain why such a
method is important in view of the light quark physics program of the RBCUKQCD
collaborations. To illustrate our method, using n_f=2+1 domainwall fermions,
we compute the nonperturbative running matrix of fourquark operators needed
in K>pipi decay and neutral kaon mixing. Our results are then compared to
perturbation theory.

We report a direct lattice calculation of the $K$ to $\pi\pi$ decay matrix
elements for both the $\Delta I=1/2$ and 3/2 amplitudes $A_0$ and $A_2$ on 2+1
flavor, domain wall fermion, $16^3\times32\times16$ lattices. This is a
complete calculation in which all contractions for the required ten, fourquark
operators are evaluated, including the disconnected graphs in which no quark
line connects the initial kaon and final twopion states. These lattice
operators are nonperturbatively renormalized using the RomeSouthampton method
and the quadratic divergences are studied and removed. This is an important but
notoriously difficult calculation, requiring high statistics on a large volume.
In this paper we take a major step towards the computation of the physical
$K\to\pi\pi$ amplitudes by performing a complete calculation at unphysical
kinematics with pions of mass 422\,MeV at rest in the kaon rest frame. With
this simplification we are able to resolve Re$(A_0)$ from zero for the first
time, with a 25% statistical error and can develop and evaluate methods for
computing the complete, complex amplitude $A_0$, a calculation central to
understanding the $\Delta =1/2$ rule and testing the standard model of CP
violation in the kaon system.

We determine the neutral kaon mixing matrix element $B_K$ in the continuum
limit with 2+1 flavors of domain wall fermions, using the Iwasaki gauge action
at two different lattice spacings. These lattice fermions have near exact
chiral symmetry and therefore avoid artificial lattice operator mixing.
We introduce a significant improvement to the conventional NPR method in
which the bare matrix elements are renormalized nonperturbatively in the
RIMOM scheme and are then converted into the MSbar scheme using continuum
perturbation theory. In addition to RIMOM, we introduce and implement four
nonexceptional intermediate momentum schemes that suppress infrared
nonperturbative uncertainties in the renormalization procedure. We compute the
conversion factors relating the matrix elements in this family of RISMOM
schemes and MSbar at oneloop order. Comparison of the results obtained using
these different intermediate schemes allows for a more reliable estimate of the
unknown higherorder contributions and hence for a correspondingly more robust
estimate of the systematic error. We also apply a recently proposed approach in
which twisted boundary conditions are used to control the Symanzik expansion
for offshell vertex functions leading to a better control of the
renormalization in the continuum limit.
We control chiral extrapolation errors by considering both the NLO SU(2)
chiral effective theory, and an analytic mass expansion. We obtain
$B_K^{\msbar}(3 GeV) = 0.529(5)_{stat}(15)_\chi(2)_{FV}(11)_{NPR}$. This
corresponds to $\hat{B}_K = 0.749(7)_{stat}(21)_\chi(3)_{FV}(15)_{NPR}$. Adding
all sources of error in quadrature we obtain $\hat{B}_K =
0.749(27)_{combined}$, with an overall combined error of 3.6%.

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.

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.

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 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).