
We review lattice results related to pion, kaon, D and Bmeson physics with
the aim of making them easily accessible to the particle physics community.
More specifically, we report on the determination of the lightquark masses,
the form factor f+(0), arising in the semileptonic K > pi transition at zero
momentum transfer, as well as the decay constant ratio fK/fpi and its
consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe
the results obtained on the lattice for some of the lowenergy constants of
SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the
determination of the BK parameter of neutral kaon mixing as well as the
additional four B parameters that arise in theories of physics beyond the
Standard Model. The latter quantities are an addition compared to the previous
review. For the heavyquark sector, we provide results for mc and mb (also new
compared to the previous review), as well as those for D and Bmeson decay
constants, form factors, and mixing parameters. These are the heavyquark
quantities most relevant for the determination of CKM matrix elements and the
global CKM unitaritytriangle fit. Finally, we review the status of lattice
determinations of the strong coupling constant alpha_s.

We present results from the first largescale study of two flavor QCD using
domain wall fermions (DWF), a chirally symmetric fermion formulation which has
proven to be very effective in the quenched approximation. We work on lattices
of size 16^3x32, with a lattice cutoff of a^{1}\approx 1.7 GeV, and dynamical
(or sea) quark masses in the range m_{strange}/2 \simle m_{sea} \simle
m_{strange}. After discussing the algorithmic and implementation issues
involved in simulating dynamical DWF, we report on the lowlying hadron
spectrum, decay constants, static quark potential, and the important kaon weak
matrix element describing indirect CP violation in the Standard Model, B_K. In
the latter case we include the effect of nondegenerate quark masses (m_s \neq
m_u = m_d), finding B_K(MSbar, 2 GeV) = 0.495(18).

The DWF formulation becomes increasingly problematic at gauge couplings for
which $a^{1}<2$ GeV, where the roughness of the gauge field leads to increased
explicit chiral symmetry breaking ($\mres$). This problem becomes especially
severe for sufficiently strong coupling where the underlying 4dimensional
Wilson theory is in the Aoki phase. We review our attempts to find a suitable
modification of the gauge and/or the fermion action which would allow the DWF
method to work reliably at stronger coupling.

We report on studies of the chiral properties of dynamical domain wall
fermions combined with the DBW2 gauge action for different gauge couplings and
fermion masses. For quenched theories, the DBW2 action gives a residual chiral
symmetry breaking much smaller than what was found with more traditional
choices for the gauge action. Our goal is to investigate the possibilities
which this and further improvements provide for the study of QCD thermodynamics
and other simulations at stronger couplings.

We study the chiral properties of quenched domain wall fermions with several
gauge actions. We demonstrate that the residual chiral symmetry breaking, which
is present for a finite number of lattice sites in the fifth dimension ($L_s$),
can be substantially suppressed using improved gauge actions. In particular the
Symanzik action, the Iwasaki action, and a renormalization group improved gauge
action, called doubly blocked Wilson (DBW2), are studied and compared to the
Wilson action. All improved gauge actions studied show a reduction in the
additive residual quark mass, $\mres$. Remarkably, in the DBW2 case $\mres$ is
roughly two orders of magnitude smaller than the Wilson gauge action at
$a^{1}=2$ GeV and $L_s=16$. Significant reduction in $\mres$ is also realized
at stronger gauge coupling corresponding to $a^{1}=1.3$ GeV. As our numerical
investigation indicates, this reduction is achieved by reducing the number of
topological lattice dislocations present in the gauge field configurations. We
also present detailed results for the quenched light hadron spectrum and the
pion decay constant using the DBW2 gauge action.

We report the results of a calculation of the K > pi pi matrix elements
relevant for the $\DIhalf$ rule and $\epe$ in quenched lattice QCD using domain
wall fermions at a fixed lattice spacing $a^{1} \sim 2$ GeV. Working in the
threequark effective theory, where only the u, d and s quarks enter and which
is known perturbatively to nexttoleading order, we calculate the lattice K
> pi and K > 0> matrix elements of dimension six, fourfermion operators.
Through lowest order chiral perturbation theory these yield K > pi pi matrix
elements, which we then normalize to continuum values through a
nonperturbative renormalization technique. For the ratio of isospin amplitudes
A_0/A_2 we find a value of $25.3 \pm 1.8$ (statistical error only) compared
to the experimental value of 22.2, with individual isospin amplitudes 1020%
below the experimental values. For $\epe$, using known central values for
standard model parameters, we calculate $(4.0 \pm 2.3) \times 10^{4}$
(statistical error only) compared to the current experimental average of $(17.2
\pm 1.8) \times 10^{4}$. Because we find a large cancellation between the I =
0 and I = 2 contributions to $\epe$, the result may be very sensitive to the
approximations employed. Among these are the use of: quenched QCD, lowest order
chiral perturbation theory and continuum perturbation theory below 1.3 GeV. We
have also calculated the kaon B parameter, B_K and find $B_{K,\bar{MS}}(2 {\rm
GeV}) = 0.532(11)$. Although currently unable to give a reliable systematic
error, we have control over statistical errors and more simulations will yield
information about the effects of the approximations on this firstprinciples
determination of these important quantities.

In the dilute instanton gas model of the QCD vacuum, one expects a strong
spatial correlation between chirality and the maxima of the Dirac eigenvectors
with small eigenvalues. Following Horvath, {\it et al.} we examine this
question using lattice gauge theory within the quenched approximation. We
extend the work of those authors by using weaker coupling, $\beta=6.0$, larger
lattices, $16^4$, and an improved fermion formulation, domain wall fermions. In
contrast with this earlier work, we find a striking correlation between the
magnitude of the chirality density, $\psi^\dagger(x)\gamma^5\psi(x)$, and the
normal density, $\psi^\dagger(x)\psi(x)$, for the lowlying Dirac eigenvectors.

We find the renormalisation coefficients of the quark field and the flavour
nonsinglet fermion bilinear operators for the domain wall fermion action, in
the regularisation independent (RI) renormalisation scheme. Our results are
from a quenched simulation, on a 16^3x32 lattice, with beta=6.0 and an extent
in the fifth dimension of 16. We also discuss the expected effects of the
residual chiral symmetry breaking inherent in a domain wall fermion simulation
with a finite fifth dimension, and study the evidence for both explicit and
spontaneous chiral symmetry breaking effects in our numerical results. We find
that the relations between different renormalisation factors predicted by
chiral symmetry are, to a good approximation, satisfied by our results and that
systematic effects due to the (low energy) spontaneous chiral symmetry breaking
and zeromodes can be controlled. Our results are compared against the
perturbative predictions for both their absolute value and renormalisation
scale dependence.

Quenched QCD simulations on three volumes, $8^3 \times$, $12^3 \times$ and
$16^3 \times 32$ and three couplings, $\beta=5.7$, 5.85 and 6.0 using domain
wall fermions provide a consistent picture of quenched QCD. We demonstrate that
the small induced effects of chiral symmetry breaking inherent in this
formulation can be described by a residual mass ($\mres$) whose size decreases
as the separation between the domain walls ($L_s$) is increased. However, at
stronger couplings much larger values of $L_s$ are required to achieve a given
physical value of $\mres$. For $\beta=6.0$ and $L_s=16$, we find
$\mres/m_s=0.033(3)$, while for $\beta=5.7$, and $L_s=48$,
$\mres/m_s=0.074(5)$, where $m_s$ is the strange quark mass. These values are
significantly smaller than those obtained from a more naive determination in
our earlier studies. Important effects of topological near zero modes which
should afflict an accurate quenched calculation are easily visible in both the
chiral condensate and the pion propagator. These effects can be controlled by
working at an appropriately large volume. A nonlinear behavior of $m_\pi^2$ in
the limit of small quark mass suggests the presence of additional infrared
subtlety in the quenched approximation. Good scaling is seen both in masses and
in $f_\pi$ over our entire range, with inverse lattice spacing varying between
1 and 2 GeV.

The domain wall formulation of lattice fermions is expected to support
accurate chiral symmetry, even at finite lattice spacing. Here we attempt to
use this new fermion formulation to simulate twoflavor, finite temperature QCD
near the chiral phase transition. In this initial study, a variety of quark
masses, domain wall heights and domain wall separations are explored using an
8^3 x 4 lattice. Both the expectation value of the Wilson line and the chiral
condensate show the temperature dependence expected for the QCD phase
transition. Further, the desired chiral properties are seen for the chiral
condensate, suggesting that the domain wall fermion formulation may be an
effective approach for the numerical study of QCD at finite temperature.

This paper has been withdrawn.

We report on simulations of quenched QCD using domain wall fermions, where we
focus on basic questions about the formalism and its ability to produce
expected low energy hadronic physics for light quarks. The work reported here
is on quenched $8^3 \times 32$ lattices at $\beta = 5.7$ and 5.85, using values
for the length of the fifth dimension between 10 and 48. We report results for
parameter choices which lead to the desired number of flavors, a study of
undamped modes in the extra dimension and hadron masses.

The domain wall approach to lattice fermions employs an additional dimension,
in which gauge fields are merely replicated, to separate the chiral components
of a Dirac fermion. It is known that in the limit of infinite separation in
this new dimension, domain wall fermions have exact zero modes, even for gauge
fields which are not smooth. We explore the effects of finite extent in the
fifth dimension on the zero modes for both smooth and nonsmooth topological
configurations and find that a fifth dimension of around ten sites is
sufficient to clearly show zero mode effects. This small value for the extent
of the fifth dimension indicates the practical utility of this technique for
numerical simulations of QCD.

We present results from numerical simulations of full, two flavor QCD
thermodynamics at N_t=4 with domain wall fermions. For the first time a
numerical simulation of the full QCD phase transition displays a low
temperature phase with spontaneous chiral symmetry breaking but intact flavor
symmetry and a high temperature phase with the full SU(2) x SU(2) chiral flavor
symmetry.

The first stage in the construction of the 0.8 Teraflops Supercomputer at
Columbia, a working, two node parallel computer, has been successfully
completed. The next stage, a 512 node, 26 Gigaflops prototype, is in its final
construction phase. A general description and current status of the hardware
and software is presented.