
The leadingorder electromagnetic and strong isospinbreaking corrections to
the ratio of $K_{\mu 2}$ and $\pi_{\mu 2}$ decay rates are evaluated for the
first time on the lattice, following a method recently proposed. The lattice
results are obtained using the gauge ensembles produced by the European Twisted
Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical quarks. Systematics effects
are evaluated and the impact of the quenched QED approximation is estimated.
Our result for the correction to the treelevel $K_{\mu 2} / \pi_{\mu 2}$ decay
ratio is $1.22\,(16) \%$ to be compared to the estimate $1.12\,(21) \%$ based
on Chiral Perturbation Theory and adopted by the Particle Data Group.

In Carrasco et al. we have recently proposed a method to calculate $O(e^2)$
electromagnetic corrections to leptonic decay widths of pseudoscalar mesons.
The method is based on the observation that the infrared divergent
contributions (that appear at intermediate stages of the calculation and that
cancel in physical quantities thanks to the BlochNordsieck mechanism) are
universal, i.e. depend on the charge and the mass of the meson but not on its
internal structure. In this talk we perform a detailed analysis of the
finitevolume effects associated with our method. In particular we show that
also the leading $1/L$ finitevolume effects are universal and perform an
analytical calculation of the finitevolume leptonic decay rate for a
pointlike meson.

We demonstrate that the leading and nexttoleading finitevolume effects in
the evaluation of leptonic decay widths of pseudoscalar mesons at $O(\alpha)$
are universal, i.e. they are independent of the structure of the meson. This is
analogous to a similar result for the spectrum but with some fundamental
differences, most notably the presence of infrared divergences in decay
amplitudes. The leading nonuniversal, structuredependent terms are of
$O(1/L^2)$ (compared to the $O(1/L^3)$ leading nonuniversal corrections in the
spectrum). We calculate the universal finitevolume effects, which requires an
extension of previously developed techniques to include a dependence on an
external threemomentum (in our case, the momentum of the final state lepton).
The result can be included in the strategy proposed in
Ref.\,\cite{Carrasco:2015xwa} for using lattice simulations to compute the
decay widths at $O(\alpha)$, with the remaining finitevolume effects starting
at order $O(1/L^2)$. The methods developed in this paper can be generalised to
other decay processes, most notably to semileptonic decays, and hence open the
possibility of a new era in precision flavour physics.

Electromagnetic effects in the leptonic decay rates $\pi^+ \to \mu^+ \nu$ and
$K^+ \to \mu^+ \nu$ are evaluated for the first time on the lattice. Following
a method recently proposed in Ref. [1] the emission of virtual photons at
leading order in the electromagnetic coupling is evaluated on the lattice and
the infrared divergence computed for a pointlike meson at finite lattice
volume is subtracted. The physical decay rate is then obtained by adding the
emission of real and virtual photons regularised with a photon mass. Using the
gauge ensembles produced by the European Twisted Mass Collaboration with $N_f =
2 + 1 + 1$ dynamical quarks the feasibility of our approach is demonstrated.
Preliminary results for the electromagnetic corrections to charged (neutral)
pion and kaon masses as well as to the leptonic decay rates of charged pions
and kaons are presented.

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 report on a calculation of the B*Bpi coupling in lattice QCD. The strong
matrix element for a B* to Bpi transition is directly related to the leading
order lowenergy constant in heavy meson chiral perturbation theory (HMChPT)
for B mesons. We carry out our calculation directly at the bquark mass using a
nonperturbatively tuned clover action that controls discretization effects of
order pa and (ma)^n for all n. Our analysis is performed on RBC/UKQCD gauge
configurations using domainwall fermions and the Iwasaki gauge action at two
lattice spacings of ainverse = 1.729(25) GeV, ainverse = 2.281(28) GeV, and
unitary pion masses down to 290 MeV. We achieve good statistical precision and
control all systematic uncertainties, giving a final result for the HMChPT
coupling g_b = 0.56(3)stat(7)sys in the continuum and at the physical
lightquark masses. This is the first calculation performed directly at the
physical bquark mass and lies in the region one would expect from carrying out
an interpolation between previous results at the charm mass and at the static
point.

The rare kaon decays $K\to\pi\ell^+\ell^$ and $K\to\pi\nu\bar{\nu}$ are
flavor changing neutral current (FCNC) processes and hence promising channels
with which to probe the limits of the standard model and to look for signs of
new physics. In this paper we demonstrate the feasibility of lattice
calculations of $K\to\pi\ell^+\ell^$ decay amplitudes for which longdistance
contributions are very significant. We show that the dominant finitevolume
corrections (those decreasing as powers of the volume) are negligibly small and
that, in the fourflavor theory, no new ultraviolet divergences appear as the
electromagnetic current $J$ and the effective weak Hamiltonian $H_W$ approach
each other. In addition, we demonstrate that one can remove the unphysical
terms which grow exponentially with the range of the integration over the time
separation between $J$ and $H_W$. We will now proceed to exploratory numerical
studies with the aim of motivating further experimental measurements of these
decays. Our work extends the earlier study by Isidori, Turchetti and Martinelli
which focussed largely on the renormalization of ultraviolet divergences. In a
companion paper we discuss the evaluation of the longdistance contributions to
$K\to\pi\nu\bar{\nu}$ decays; these contributions are expected to be at the
level of a few percent for $K^+$ decays.

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.

In this paper, for the first time a method is proposed to compute
electromagnetic effects in hadronic processes using lattice simulations. The
method can be applied, for example, to the leptonic and semileptonic decays of
light or heavy pseudoscalar mesons. For these quantities the presence of
infrared divergences in intermediate stages of the calculation makes the
procedure much more complicated than is the case for the hadronic spectrum, for
which calculations already exist. In order to compute the physical widths,
diagrams with virtual photons must be combined with those corresponding to the
emission of real photons. Only in this way do the infrared divergences cancel
as first understood by Bloch and Nordsieck in 1937. We present a detailed
analysis of the method for the leptonic decays of a pseudoscalar meson. The
implementation of our method, although challenging, is within reach of the
present lattice technology.

We report on the first complete calculation of the $K_LK_S$ mass difference,
$\Delta M_K$, using lattice QCD. The calculation is performed on a 2+1 flavor,
domain wall fermion ensemble with a 330MeV pion mass and a 575 MeV kaon mass.
We use a quenched charm quark with a 949 MeV mass to implement
GlashowIliopoulosMaiani cancellation. For these heavierthanphysical
particle masses, we obtain $\Delta M_K =3.19(41)(96)\times 10^{12}$ MeV, quite
similar to the experimental value. Here the first error is statistical and the
second is an estimate of the systematic discretization error. An interesting
aspect of this calculation is the importance of the disconnected diagrams, a
dramatic failure of the OZI rule.

The RBC and UKQCD collaborations have recently proposed a procedure for
computing the K_LK_S mass difference. A necessary ingredient of this procedure
is the calculation of the (nonexponential) finitevolume corrections relating
the results obtained on a finite lattice to the physical values. This requires
a significant extension of the techniques which were used to obtain the
LellouchLuscher factor, which contains the finitevolume corrections in the
evaluation of nonleptonic kaon decay amplitudes. We review the status of our
study of this issue and, although a complete proof is still being developed,
suggest the form of these corrections for general volumes and a strategy for
taking the infinitevolume limit. The general result reduces to the known
corrections in the special case when the volume is tuned so that there is a
twopion state degenerate with the kaon.

We report on a calculation of the $B^*B\pi$ coupling in lattice QCD. The
strong matrix element $\langle B \pi  B^*\rangle$ is directly related to the
leading order lowenergy constant in heavy meson chiral perturbation theory
(HM$\chi$PT) for $B$mesons. We carry out our calculation directly at the
$b$quark mass using a nonperturbatively tuned clover action that controls
discretisation effects of order $\vec{p}a$ and $(ma)^n$ for all $n$. Our
analysis is performed on RBC/UKQCD gauge configurations using domain wall
fermions and the Iwasaki gauge action at two lattice spacings of
$a^{1}=1.73(3)$ GeV, $a^{1}=2.28(3)$ GeV, and unitary pion masses down to 290
MeV. We achieve good statistical precision and control all systematic
uncertainties, giving a final result for the HM$\chi$PT coupling $g_b =
0.569(48)_{stat}(59)_{sys}$ in the continuum and at the physical lightquark
masses. This is the first calculation performed directly at the physical
$b$quark mass and lies in the region one would expect from carrying out an
interpolation between previous results at the charm mass and at the static
point.

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

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.

The large mass of the ninth pseudoscalar meson, the $\eta^\prime$, is
believed to arise from the combined effects of the axial anomaly and the gauge
field topology present in QCD. We report a realistic, 2+1 flavor, lattice QCD
calculation of the $\eta$ and $\eta^\prime$ masses and mixing which confirms
this picture. The physical eigenstates show small octetsinglet mixing with a
mixing angle of $\theta = 14.1(2.8)^\circ$. Extrapolation to physical light
quark mass gives, with statistical errors only, $m_\eta=573(6)$ MeV and
$m_{\eta^\prime}=947(142)$ MeV, consistent with the experimental values of 548
MeV and 958 MeV.

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.

The integrated branching fraction of the process $B\to X_s l^+l^$ is
dominated by resonance background from narrow charmonium states, such as $B\to
X_s\psi\to X_s l^+l^$, which exceeds the nonresonant charmloop contribution
by two orders of magnitude. The origin of this fact is discussed in view of the
general expectation of quarkhadron duality. The situation in $B\to X_s l^+l^$
is contrasted with charmpenguin amplitudes in twobody hadronic B decays of
the type $B\to\pi\pi$, for which it is demonstrated that resonance effects and
the potentially nonperturbative $c\bar c$ threshold region do not invalidate
the standard picture of QCD factorization. This holds irrespective of whether
the charm quark is treated as a light or a heavy quark.

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