
We introduce a new algorithm for finding the eigenvalues and eigenvectors of
Hermitian matrices within a specified region, based upon the LANSO algorithm of
Parlett and Scott. It uses selective reorthogonalization to avoid the
duplication of eigenpairs in finiteprecision arithmetic, but uses a new bound
to decide when such reorthogonalization is required, and only reorthogonalizes
with respect to eigenpairs within the region of interest. We investigate its
performance for the Hermitian WilsonDirac operator (\gamma_5D) in lattice
quantum chromodynamics, and compare it with previous methods.

It is usually assumed that the infrared divergences that occur for virtual
photons in Feynman diagrams cancel against the corresponding divergences in
phase space integrals for the emission of real soft photons. If this were the
case then why do the corresponding phase space integrals for the absorption of
real soft photons not diverge too, as one might expect on grounds of
timereversal symmetry, and if so what cancels these divergences? We show that
this apparent time asymmetry is an artifact of the way that the diagrams are
calculated, and we introduce a procedure for such calculations that is
manifestly time symmetric. Moreover, we show that it is not mathematically
meaningful to distinguish arbitrarily soft massless particles of positive and
negative energy.

QCD results are presented for a 2+1 flavour fermion clover action (which we
call the SLiNC action). A method of tuning the quark masses to their physical
values is discussed. In this method the singlet quark mass is kept fixed, which
solves the problem of different renormalisations (for singlet and nonsinglet
quark masses) occuring for nonchirally invariant lattice fermions. This
procedure enables a wide range of quark masses to be probed, including the case
with a heavy updown quark mass and light strange quark mass. Preliminary
results show the correct splittings for the baryon (octet and) decuplet
spectrum.

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

We present results for light meson masses and psedoscalar meson decay
constants in 2+1 flavour domain wall QCD with the DBW2 and Iwasaki gauge
actions, using lattices with linear sizes in the range 1.6 to 2.2fm and $u$ and
$d$ quark masses as low as one quarter of the strange quark mass. All data were
generated on the QCDOC machines at the University of Edinburgh and Brookhaven
National Laboratory. Despite large residual masses and a limited number of sea
quark mass values with which to perform chiral extrapolations, our results
agree with experiment and scale within errors.

We present results for some of the light baryon masses and their excited
states in 2+1 flavour domain wall QCD. We considered several lattice spacings,
with the DBW2 and Iwasaki gauge actions and different sea quark masses on a
volume of $16^3\times32$ and a fifth dimension of size 8. All data were
generated on the QCDOC machines. Despite large residual massses and a limited
number of sea quark mass values with which to perform chiral extrapolations,
our results are in reasonable agreement with experiment and scale within
errors. Finite size effects on most ensembles appear to be small.

We present first results from a simulation of quenched overlap fermions with
improved gauge field action. Among the quantities we study are the spectral
properties of the overlap operator, the chiral condensate and topological
charge, quark and hadron masses, and selected nucleon matrix elements. To make
contact with continuum physics, we compute the renormalization constants of
quark bilinear operators in perturbation theory and beyond.

We study QCD with two flavors of dynamical Wilson fermions at $\beta = 5.5$
and three values of $\kappa$. The corresponding pion masses are 0.375, 0.324
and 0.262 in lattice units, with pion to rho mass ratios of 0.76, 0.71 and
0.62, respectively. We use the configurations to compute the heavy quark
potential, leading to lattice spacings of 0.110, 0.105 and 0.099 fm, and to
compute spectroscopy for several different valence quark $\kappa$'s.

The theoretical justification of the Hybrid Monte Carlo algorithm depends
upon the molecular dynamics trajectories within it being exactly reversible. If
computations were carried out with exact arithmetic then it would be easy to
ensure such reversibility, but the use of approximate floating point arithmetic
inevitably introduces violations of reversibility. In the absence of evidence
to the contrary, we are usually prepared to accept that such rounding errors
can be made small enough to be innocuous, but in certain circumstances they are
exponentially amplified and lead to blatantly erroneous results. We show that
there are two types of instability of the molecular dynamics trajectories which
lead to this behavior, instabilities due to insufficiently accurate numerical
integration of Hamilton's equations, and intrinsic chaos in the underlying
continuous fictitious time equations of motion themselves. We analyze the
former for free field theory, and show that it is essentially a finite volume
effect. For the latter we propose a hypothesis as to how the Liapunov exponent
describing the chaotic behavior of the fictitious time equations of motion for
an asymptotically free quantum field theory behaves as the system is taken to
its continuum limit, and explain why this means that instabilities in molecular
dynamics trajectories are not a significant problem for Hybrid Monte Carlo
computations. We present data for pure $SU(3)$ gauge theory and for QCD with
dynamical fermions on small lattices to illustrate and confirm some of our
results.

We test the asymptotic scaling behavior of stateoftheart simulations of
QCD with two flavors of light Wilson fermions. This is done by matching $\pi$
and $\rho$ masses on lattices of size $16^3\times32$ and $8^3\times16$. We find
that at $\beta=6/g^2=5.3$ matching is not possible over a range extending down
to $\beta=3.5$. The large lattice data at $\beta=5.5$ matches the small lattice
values at $\beta=4.9(1)$ leading to a shift $\Delta\beta=0.6(1)$, considerably
larger than the perturbative prediction of $0.45$. In both cases we conclude
that the simulations are very far from the asymptotic scaling region.

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.

We test the cooling algorithm with gluonic and staggered hadronic
spectroscopy on $SU(3)$ gauge field configurations generated with two flavors
of staggered dynamical fermions. We find cooling is not reliable as the basis
for improved hadronic operators. We also find that performing cooling sweeps to
reveal more clearly the topological properties of the gauge fields eliminates
the spin structure of the hadron spectrum.

We present data for the scaling behavior of lattice QCD with two flavors of
light Wilson fermions. This is done by matching $\pi$ n and $\rho$ masses at
the two lattice sizes of $16^3\times32 $ and $8^3\times16$. We find that at
$\beta=6/g^2\equiv5.3$ the matching does not occur over a range extending down
to $\beta=3.5$. For $\beta=5.5$ matching may be achieved at $\beta=4.84.9$,
leading to a $\Delta \beta =0.60.7$ which is higher than the perturbative
2loop value of 0.45. In both cases we conclude that the simulations are very
far from the perturbative scaling region.

We compute the heavy quark potential on configurations generated by the
HEMCGC collaboration with dynamical staggered fermions at $6/g^2 = 5.6$ and
with dynamical Wilson fermions at $6/g^2 = 5.3$. The computations are done on
$16^3 \times 32$ lattices, corresponding to physical sizes of about 1.6 and 2.3
fm, respectively. Up to the distances probed no sign of string breaking is
detectable. We also compute the recently proposed scale $r_0$ defined by $r_0^2
F(r_0) = 1.65$.

We have measured some simple matrix elements for pseudoscalar and vector
mesons made of Wilson valence quarks and staggered sea quarks at $\beta=5.6$ at
sea quark masses $am_q=0.01$ and 0.025. Our measurements include the decay
constants of pseudoscalars (including $f_D$), the wave function at the origin
(or decay constant) of vector mesons, and the calculation of quark masses from
current algebra. The effects of sea quarks on the simulations are small. We
make comparisons to quenched simulations at similar values of the lattice
spacing ($1/a \simeq 2$ GeV).

We describe recent results concerning the behavior of lattice QCD with light
dynamical Wilson and Staggered quarks. We show that it is possible to reach
regions of parameter space with light pions $m_\pi\approx 0.2/a$ using Wilson
fermions. If the Hybrid Molecular Dynamics (HMD) algorithm is used with the
same parameters it gives incorrect results. We also present preliminary results
using a higherorder integration scheme.

We present an analysis of hadronic spectroscopy for Wilson valence quarks
with dynamical staggered fermions at lattice coupling $6/g^2 = \beta=5.6$ at
sea quark mass $am_q=0.01$ and 0.025, and of Wilson valence quarks in quenched
approximation at $\beta=5.85$ and 5.95, both on $16^3 \times 32$ lattices. We
make comparisons with our previous results with dynamical staggered fermions at
the same parameter values but on $16^4$ lattices doubled in the temporal
direction.

We present preliminary results from the 1991 HEMCGC simulations with
staggered dynamical fermions on a $16^3 \times 32$ lattice at $\beta = 5.6$
with sea quark masses $am_q = 0.025$ and 0.01. The spectroscopy was done both
for staggered valence quarks with mass equal to the sea quark masses and for
Wilson valence quarks at six different values for $\kappa$, 0.1320, 0.1410,
0.1525, 0.1565, 0.1585, and 0.1600. In addition to the measurements performed
in our earlier work, we also measured the $\Delta$ and other `extended' hadrons
for staggered valence quarks and pseudoscalar decay constants and vector meson
matrix elements, the wave function at the origin, for Wilson valence quarks.

We describe recent results obtained as part of the High Energy Monte Carlo
Grand Challenge (HEMCGC) project concerning the behaviour of lattice QCD with
light dynamical Wilson quarks. We show that it is possible to reach regions of
parameter space with light pions $m_\pi\ll0.2/a$, but that the equilibration
time for such a system is at least of the order of 1,000 unitlength Hybrid
Monte Carlo (HMC) trajectories (about a Gigaflop/secyear). If the Hybrid
Molecular Dynamics (HMD) algorithm is used with the same parameters it gives
incorrect results.

Lattice QCD with 2 light staggered quark flavours is being simulated on a
$16^3\times8$ lattice to study the transition from hadronic matter to a quark
gluon plasma. We have completed runs at $m_q=0.0125$ and are extending this to
$m_q=0.00625$. We also examine the addition of a nondynamical "strange" quark.
Thermodynamic order parameters are being measured across the transition and
further into the plasma phase, as are various screening lengths. No evidence
for a first order transition is seen, and we estimate the transition
temperature to be $TY_c=143(7) MeV$.