
Observation of neutrinoless double beta decay, a lepton number violating
process that has been proposed to clarify the nature of neutrino masses, has
spawned an enormous worldwide experimental effort. Relating nuclear decay
rates to highenergy, beyond the Standard Model (BSM) physics requires detailed
knowledge of nonperturbative QCD effects. Using lattice QCD and taking
advantage of effective field theory methods, we compute the modelindependent
leadingorder matrix elements of shortrange operators, which arise due to
heavy BSM mediators, that contribute to this decay. Contributions from
shortrange operators may prove to be equally important to or even more
important than those from longrange Majorana neutrino exchange.

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 compute the hadronic matrix elements of the fourquark operators needed
for the study of $K^0{\bar K^0}$ mixing beyond the Standard Model. We have
used $n_f=2+1$ flavours of domain wall fermion (DWF) at two values of the
lattice spacing ($a\approx0.08$ and $a\approx0.11\text{fm}$) and with lightest
unitary pion mass of $\approx300\,\text{MeV}$. Renormalisation is performed
nonperturbatively and the impact of different intermediate momentum schemes is
investigated.

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 study the heavy quark mass dependence of selected observables constructed
from heavylight meson correlation functions in smallvolume twoflavour
lattice QCD after taking the continuum limit. The light quark mass is tuned to
zero, whereas the range of available heavy quark masses $m_h$ covers a region
extending from around the charm to beyond the bottom quark mass scale. This
allows entering the asymptotic massscaling regime as $1/m_h \to 0$ and
performing wellcontrolled extrapolations to the infinitemass limit. Our
results are then compared to predictions obtained in the static limit of
continuum Heavy Quark Effective Theory (HQET), in order to verify
nonperturbatively that HQET is an effective theory of QCD. While in general we
observe a nice agreement at the few% level, we find it to be less convincing
for the smallvolume pseudoscalar decay constant when perturbative matching is
involved.

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

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.

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 present a computation of Bmeson decay constants from lattice QCD
simulations within the framework of Heavy Quark Effective Theory for the
bquark. The nexttoleading order corrections in the HQET expansion are
included nonperturbatively. Based on Nf=2 gauge field ensembles, covering
three lattice spacings a (0.080.05)fm and pion masses down to 190MeV, a
variational method for extracting hadronic matrix elements is used to keep
systematic errors under control. In addition we perform a careful
autocorrelation analysis in the extrapolation to the continuum and to the
physical pion mass limits. Our final results read fB=186(13)MeV, fBs=224(14)MeV
and fBs/fB=1.203(65). A comparison with other results in the literature does
not reveal a dependence on the number of dynamical quarks, and effects from
truncating HQET appear to be negligible.

We report our final estimate of the bquark mass from $N_f=2$ lattice QCD
simulations using Heavy Quark Effective Theory nonperturbatively matched to
QCD at $O(1/m_h)$. Treating systematic and statistical errors in a conservative
manner, we obtain $\overline{m}_{\rm b}^{\overline{\rm MS}}(2 {\rm
GeV})=4.88(15)$ GeV after an extrapolation to the physical point.

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.

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.

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.

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.

We report the final results of the ALPHA collaboration for some Bphysics
observables: $f_B$, $f_{B_s}$ and $m_b$. We employ CLS configurations with 2
flavors of $O(a)$ improved Wilson fermions in the sea and pion masses ranging
down to 190 MeV. The bquark is treated in HQET to order $1/m_b$. The
renormalization, the matching and the improvement were performed
nonperturbatively, and three lattice spacings reaching $a=0.048$ fm are used
in the continuum extrapolation.

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.

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

The 2012 PDG reports a tension at the level of $3 \sigma$ between two
exclusive determinations of $V_{ub}$. They are obtained by combining the
experimental branching ratios of $B \to \tau \nu$ and $B \to \pi l \nu$
(respectively) with a theoretical computation of the hadronic matrix elements
$\fB$ and the $B \to \pi$ form factor $f_+(q^2)$. To understand the tension,
improved precision and a careful analysis of the systematics involved are
necessary. We report the results of the ALPHA collaboration for $\fB$ from the
lattice with 2 flavors of $O(a)$ improved Wilson fermions. We employ HQET,
including $1/m_b$ corrections, with pion masses ranging down to $\approx$ 190
MeV. Renormalization and matching were performed nonperturbatively, and three
lattice spacings reaching $a^{1}\approx 4.1$ GeV are used in the continuum
extrapolation. We also present progress towards a computation of $f_+(q^2)$, to
directly compare two independent exclusive determinations of $V_{ub}$ with
each other and with inclusive determinations. Additionally, we report on
preliminary results for $\fBq{s}$, needed for the analysis of $B_s \to
\mu^+\mu^$.}

We present our analysis of B physics quantities using nonperturbatively
matched Heavy Quark Effective Theory (HQET) in Nf= 2 lattice QCD on the CLS
ensembles. Using alltoall propagators, HYPsmeared static quarks, and the
Generalized Eigenvalue Problem (GEVP) approach with a conservative plateau
selection procedure, we are able to systematically control all sources of
error. With significantly increased statistics compared to last year, our
preliminary results are mb (mb) = 4.22(10)(4)z GeV for the MS bquark mass, and
fB = 193(9)stat (4)_\chi MeV and fBs = 219(12)stat MeV for the Bmeson decay
constants.

We compute the hadronic matrix elements of the fourquark operators needed
for the study of neutral kaon mixing beyond the Standard Model (SM). We use
nf=2+1 flavours of domainwall fermions (DWF) which exhibit good chiralflavour
symmetry. The renormalization is performed nonperturbatively through the
RIMOM scheme and our results are converted perturbatively to MSbar. The
computation is performed on a single lattice spacing a=0.086 fm with a lightest
unitary pion mass of 290 MeV. The various systematic errors, including the
discretisation effects, are estimated and discussed. Our results confirm a
previous quenched study, where large ratios of nonSM to SM matrix elements
were obtained.

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 compute the renormalization factors of fourquark operators needed for the
study of $K\to\pi\pi$ decay in the $\Delta I=3/2$ channel. We evaluate the
Zfactors at a low energy scale ($\mu_0=1.145 \GeV$) using four different
nonexceptional RISMOM schemes on a large, coarse lattice ($a\sim 0.14\fm$) on
which the bare matrix elements are also computed. Then we compute the
universal, nonperturbative, scale evolution matrix of these renormalization
factors between $\mu_0$ and $3\GeV$. We give the numerical results for the
different steps of the computation in two different nonexceptional lattice
schemes, and the connection to $\msbar$ at $3\GeV$ is made using oneloop
perturbation theory.