
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.

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 present results on both the restoration of the spontaneously broken chiral
symmetry and the effective restoration of the anomalously broken U(1)_A
symmetry in finite temperature QCD at zero chemical potential using lattice
QCD. We employ domain wall fermions on lattices with fixed temporal extent
N_\tau = 8 and spatial extent N_\sigma = 16 in a temperature range of T = 139 
195 MeV, corresponding to lattice spacings of a \approx 0.12  0.18 fm. In
these calculations, we include two degenerate light quarks and a strange quark
at fixed pion mass m_\pi = 200 MeV. The strange quark mass is set near its
physical value. We also present results from a second set of finite temperature
gauge configurations at the same volume and temporal extent with slightly
heavier pion mass. To study chiral symmetry restoration, we calculate the
chiral condensate, the disconnected chiral susceptibility, and susceptibilities
in several meson channels of different quantum numbers. To study U(1)_A
restoration, we calculate spatial correlators in the scalar and pseudoscalar
channels, as well as the corresponding susceptibilities. Furthermore, we also
show results for the eigenvalue spectrum of the Dirac operator as a function of
temperature, which can be connected to both U(1)_A and chiral symmetry
restoration via BanksCasher relations.

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 study the region of the QCD phase transition using 2+1 flavors of domain
wall fermions (DWF) and a $16^3 \times 8$ lattice volume with a fifth dimension
of $L_s = 32$. The disconnected light quark chiral susceptibility, quark number
susceptibility and the Polyakov loop suggest a chiral and deconfining crossover
transition lying between 155 and 185 MeV for our choice of quark mass and
lattice spacing. In this region the lattice scale deduced from the Sommer
parameter $r_0$ is $a^{1} \approx 1.3$ GeV, the pion mass is $\approx 300$ MeV
and the kaon mass is approximately physical. The peak in the chiral
susceptibility implies a pseudo critical temperature $T_c = 171(10)(17)$ MeV
where the first error is associated with determining the peak location and the
second with our unphysical light quark mass and nonzero lattice spacing. The
effects of residual chiral symmetry breaking on the chiral condensate and
disconnected chiral susceptibility are studied using several values of the
valence $L_s$.

We present results on the equation of state in QCD with two light quark
flavors and a heavier strange quark. Calculations with improved staggered
fermions have been performed on lattices with temporal extent Nt =4 and 6 on a
line of constant physics with almost physical quark mass values; the pion mass
is about 220 MeV, and the strange quark mass is adjusted to its physical value.
High statistics results on large lattices are obtained for bulk thermodynamic
observables, i.e. pressure, energy and entropy density, at vanishing quark
chemical potential for a wide range of temperatures, 140 MeV < T < 800 MeV. We
present a detailed discussion of finite cutoff effects which become
particularly significant for temperatures larger than about twice the
transition temperature. At these high temperatures we also performed
calculations of the trace anomaly on lattices with temporal extent Nt=8.
Furthermore, we have performed an extensive analysis of zero temperature
observables including the light and strange quark condensates and the static
quark potential at zero temperature. These are used to set the temperature
scale for thermodynamic observables and to calculate renormalized observables
that are sensitive to deconfinement and chiral symmetry restoration and become
order parameters in the infinite and zero quark mass limits, respectively.

We present the first results for neutral kaon mixing using 2+1 flavors of
domain wall fermions. A new approach is used to extrapolate to the physical up
and down quark masses from our numerical studies with pion masses in the range
240  420 MeV; only $SU(2)_L \times SU(2)_R$ chiral symmetry is assumed and
the kaon is not assumed to be light. Our main result is $B_K^{\bar{\rm MS}}(2
\mathrm{GeV}) = 0.524(10)(28)$ where the first error is statistical and the
second incorporates estimates for all systematic errors.

We present a study of the flavor symmetry breaking in the pion spectrum for
various improved staggered fermion actions. To study the effects of link
fattening and tadpole improvement, we use three different variants of the p4
action  p4fat3, p4fat7, and p4fat7tad. These are compared to Asqtad and also
to naive staggered. To study the pattern of symmetry breaking, we measure all
15 meson masses in the 4flavor staggered theory. The measurements are done on
a quenched gauge background, generated using a oneloop improved Symanzik
action with $\beta=10/g^2 = 7.40, 7.75,$ and 8.00, corresponding to lattice
spacings of approximately a = .31 fm., .21 fm., and .14 fm. We also study how
the lattice scale set by the $\rho$ mass on each of these ensembles compares to
one set by the static quark potential.

We present results for the static interquark potential, light meson and
baryon masses, and light pseudoscalar meson decay constants obtained from
simulations of domain wall QCD with one dynamical flavour approximating the $s$
quark, and two degenerate dynamical flavours with input bare masses ranging
from $m_s$ to $m_s/4$ approximating the $u$ and $d$ quarks. We compare these
quantities obtained using the Iwasaki and DBW2 improved gauge actions, and
actions with larger rectangle coefficients, on $16^3\times32$ lattices. We seek
parameter values at which both the chiral symmetry breaking residual mass due
to the finite lattice extent in the fifth dimension and the Monte Carlo time
history for topological charge are acceptable for this set of quark masses at
lattice spacings above 0.1 fm. We find that the Iwasaki gauge action is best,
demonstrating the feasibility of using QCDOC to generate ensembles which are
good representations of the QCD path integral on lattices of up to 3 fm in
spatial extent with lattice spacings in the range 0.090.13 fm. Despite large
residual masses and a limited number of sea quark mass values with which to
perform chiral extrapolations, our results for light hadronic physics scale and
agree with experimental measurements within our statistical uncertainties.

We study the finite temperature transition in QCD with three flavors of equal
masses using the R and RHMC algorithm on lattices with temporal extent
N_{\tau}=4 and 6. For the transition temperature in the continuum limit we find
r_0 T_c=0.429(8) for the light pseudoscalar mass corresponding to the end
point of the 1st order transition region. When comparing the results obtained
with the R and RHMC algorithms for p4fat3 action we see no significant
stepsize errors down to a lightest pseudoscalar mass of m_{ps} r_0=0.4.

We present a detailed calculation of the transition temperature in QCD with
two light and one heavier (strange) quark mass on lattices with temporal extent
N_t =4 and 6. Calculations with improved staggered fermions have been performed
for various light to strange quark mass ratios in the range, 0.05 <= m_l/m_s <=
0.5, and with a strange quark mass fixed close to its physical value. From a
combined extrapolation to the chiral (m_l > 0) and continuum (aT = 1/N_t > 0)
limits we find for the transition temperature at the physical point T_c r_0 =
0.457(7) where the scale is set by the Sommerscale parameter r_0 defined as
the distance in the static quark potential at which the slope takes on the
value, (dV_qq(r)/dr)_r=r_0 = 1.65/r_0^2. Using the currently best known value
for r_0 this translates to a transition temperature T_c = 192(7)(4)MeV. The
transition temperature in the chiral limit is about 3% smaller. We discuss
current ambiguities in the determination of T_c in physical units and also
comment on the universal scaling behavior of thermodynamic quantities in the
chiral limit.

We present numerical results for the kaon Bparameter, B_K, determined in the
quenched approximation of lattice QCD. Our simulations are performed using
domainwall fermions and the renormalization group improved, DBW2 gauge action
which combine to give quarks with good chiral symmetry at finite lattice
spacing. Operators are renormalized nonperturbatively using the RI/MOM scheme.
We study scaling by performing the simulation on two different lattices with
a^{1} = 1.982(30) and 2.914(54) GeV. We combine this quenched scaling study
with an earlier calculation of B_K using two flavors of dynamical, domainwall
quarks at a single lattice spacing to obtain
B_K(MS,NDR,mu=2GeV)=0.563(21)(39)(30), were the first error is statistical, the
second systematic (without quenching errors) and the third estimates the error
due to quenching.

We discuss the implementation and properties of the quenched approximation in
the calculation of the leftright, strong penguin contributions (i.e. Q_6) to
epsilon'/epsilon. The coefficient of the new chiral logarithm, discovered by
Golterman and Pallante, which appears at leading order in quenched chiral
perturbation theory is evaluated using both the method proposed by those
authors and by an improved approach which is free of power divergent
corrections. The result implies a large quenching artifact in the contribution
of Q_6 to epsilon'/epsilon. This failure of the quenched approximation affects
only the strong penguin operators and so does not affect the Q_8 contribution
to epsilon'/epsilon nor Re A_0, Re A_2 and thus the Delta I=1/2 rule at tree
level in chiral perturbation theory.

The architecture of a new class of computers, optimized for lattice QCD
calculations, is described. An individual node is based on a single integrated
circuit containing a PowerPC 32bit integer processor with a 1 Gflops 64bit
IEEE floating point unit, 4 Mbyte of memory, 8 Gbit/sec nearestneighbor
communications and additional control and diagnostic circuitry. The machine's
name, QCDOC, derives from ``QCD On a Chip''.