We review lattice results related to pion, kaon, D- and B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark 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 low-energy 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 heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. Finally, we review the status of lattice determinations of the strong coupling constant alpha_s.

We present results from the first large-scale 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 cut-off 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 low-lying 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 non-degenerate quark masses (m_s \neq m_u = m_d), finding B_K(MS-bar, 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 4-dimensional 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 three-quark effective theory, where only the u, d and s quarks enter and which is known perturbatively to next-to-leading order, we calculate the lattice K --> pi and K --> |0> matrix elements of dimension six, four-fermion operators. Through lowest order chiral perturbation theory these yield K --> pi pi matrix elements, which we then normalize to continuum values through a non-perturbative 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 10-20% 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 first-principles 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 low-lying Dirac eigenvectors.

We find the renormalisation coefficients of the quark field and the flavour non-singlet 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 zero-modes 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 non-linear 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 two-flavor, 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 non-smooth 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.