• 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 finite-precision 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 Wilson--Dirac 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 time-reversal 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 non-singlet quark masses) occuring for non-chirally invariant lattice fermions. This procedure enables a wide range of quark masses to be probed, including the case with a heavy up-down 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 next-to-leading 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 non-perturbative renormalization to relate lattice regularized quark masses to RI-MOM 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 state-of-the-art 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.8-4.9$, leading to a $\Delta \beta =0.6-0.7$ which is higher than the perturbative 2-loop 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 higher-order 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 pseudo-scalar 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 unit-length Hybrid Monte Carlo (HMC) trajectories (about a Gigaflop/sec-year). 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 non-dynamical "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$.