
Lattice QCD calculations including the effects of one or more nondegenerate
sea quark flavors are conventionally performed using the Rational Hybrid Monte
Carlo (RHMC) algorithm, which computes the square root of the determinant of
$\mathscr{D}^{\dagger} \mathscr{D}$, where $\mathscr{D}$ is the Dirac operator.
The special case of two degenerate quark flavors with the same mass is
described directly by the determinant of $\mathscr{D}^{\dagger} \mathscr{D}$
 in particular, no square root is necessary  enabling a variety of
algorithmic developments, which have driven down the cost of simulating the
light (up and down) quarks in the isospinsymmetric limit of equal masses. As a
result, the relative cost of single quark flavors  such as the strange or
charm  computed with RHMC has become more expensive. This problem is even
more severe in the context of our measurements of the $\Delta I = 1/2$ $K
\rightarrow \pi \pi$ matrix elements on lattice ensembles with $G$parity
boundary conditions, since $G$parity is associated with a doubling of the
number of quark flavors described by $\mathscr{D}$, and thus RHMC is needed for
the isospinsymmetric light quarks as well. In this paper we report on our
implementation of the exact one flavor algorithm (EOFA) introduced by the TWQCD
collaboration for simulations including single flavors of domain wall quarks.
We have developed a new preconditioner for the EOFA Dirac equation, which both
reduces the cost of solving the Dirac equation and allows us to reuse the bulk
of our existing highperformance code. Coupling these improvements with careful
tuning of our integrator, the time per accepted trajectory in the production of
our 2+1 flavor $G$parity ensembles with physical pion and kaon masses has been
decreased by a factor of 4.2.

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 have performed global fits of $f_{\pi}$ and $m_{\pi}$, from a variety of
RBCUKQCD domain wall fermion ensembles, to $SU(2)$ partially quenched chiral
perturbation theory at 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
convergence of the expansion and use our large set of low energy constants to
make predictions for the pion mass splitting 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 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 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%.

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.

We calculate the equation of state in 2+1 flavor QCD at finite temperature
with physical strange quark mass and almost physical light quark masses using
lattices with temporal extent Nt=8. Calculations have been performed with two
different improved staggered fermion actions, the asqtad and p4 actions.
Overall, we find good agreement between results obtained with these two O(a^2)
improved staggered fermion discretization schemes. A comparison with earlier
calculations on coarser lattices is performed to quantify systematic errors in
current studies of the equation of state. We also present results for
observables that are sensitive to deconfining and chiral aspects of the QCD
transition on Nt=6 and 8 lattices. We find that deconfinement and chiral
symmetry restoration happen in the same narrow temperature interval. In an
Appendix we present a simple parametrization of the equation of state that can
easily be used in hydrodynamic model calculations. In this parametrization we
also incorporated an estimate of current uncertainties in the lattice
calculations which arise from cutoff and quark mass effects. We estimate these
systematic effects to be about 10 MeV

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.

QCDOC is a massively parallel supercomputer whose processing nodes are based
on an applicationspecific integrated circuit (ASIC). This ASIC was
customdesigned so that crucial lattice QCD kernels achieve an overall
sustained performance of 50% on machines with several 10,000 nodes. This strong
scalability, together with low power consumption and a price/performance ratio
of $1 per sustained MFlops, enable QCDOC to attack the most demanding lattice
QCD problems. The first ASICs became available in June of 2003, and the testing
performed so far has shown all systems functioning according to specification.
We review the hardware and software status of QCDOC and present performance
figures obtained in real hardware as well as in simulation.

A status report is given of the QCDOC project, a massively parallel computer
optimized for lattice QCD using systemonachip technology. We describe
several of the hardware and software features unique to the QCDOC architecture
and present performance figures obtained from simulating the current VHDL
design of the QCDOC chip with singlecycle accuracy.