Lattice QCD Production on Commodity Clusters at Fermilab

We describe the construction and results to date of Fermilab's three Myrinet-networked lattice QCD production clusters (an 80-node dual Pentium III cluster, a 48-node dual Xeon cluster, and a 128-node dual Xeon cluster). We examine a number of aspects of performance of the MILC lattice QCD code running on these clusters.Comment: Talk from the 2003 Computing in High Energy and Nuclear Physics (CHEP03), La Jolla, Ca, USA, March 2003, 6 pages, LaTeX, 8 eps figures. PSN TUIT00

Scaling of Pseudo-Critical Couplings in Two-Flavour QCD

We study the scaling behaviour of the pseudo-critical couplings for the chiral phase transition in two-flavour QCD. We show that all existing results from lattice simulations on lattices with temporal extent $N_\tau = 4$, 6 and 8 can be mapped onto a universal scaling curve. The relevant combination of critical exponents, $\beta\delta$, is consistent with the scaling behaviour expected for a second order phase transition with $O(4)$ exponents. At present, scaling according to the $O(2)$ symmetry group can, however, not be ruled out.Comment: 8 pages, NSF-ITP 93-12

An adaptive pseudo-spectral method for reaction diffusion problems

The spectral interpolation error was considered for both the Chebyshev pseudo-spectral and Galerkin approximations. A family of functionals I sub r (u), with the property that the maximum norm of the error is bounded by I sub r (u)/J sub r, where r is an integer and J is the degree of the polynomial approximation, was developed. These functionals are used in the adaptive procedure whereby the problem is dynamically transformed to minimize I sub r (u). The number of collocation points is then chosen to maintain a prescribed error bound. The method is illustrated by various examples from combustion problems in one and two dimensions

Heavy Dynamical Fermions in Lattice QCD

It is expected that the only effect of heavy dynamical fermions in QCD is to renormalize the gauge coupling. We derive a simple expression for the shift in the gauge coupling induced by $N_f$ flavors of heavy fermions. We compare this formula to the shift in the gauge coupling at which the confinement-deconfinement phase transition occurs (at fixed lattice size) from numerical simulations as a function of quark mass and $N_f$. We find remarkable agreement with our expression down to a fairly light quark mass. However, simulations with eight heavy flavors and two light flavors show that the eight flavors do more than just shift the gauge coupling. We observe confinement-deconfinement transitions at $\beta=0$ induced by a large number of heavy quarks. We comment on the relevance of our results to contemporary simulations of QCD which include dynamical fermions.Comment: COLO-HEP-311, 26 pages and 6 postscript figures; file is a shar file and all macros are (hopefully) include

Computer programs for reduction of microphotometer data

Five computer programs for analyzing magnetic tape recordings of digital data from microphotomete

The beta function and equation of state for QCD with two flavors of quarks

We measure the pressure and energy density of two flavor QCD in a wide range of quark masses and temperatures. The pressure is obtained from an integral over the average plaquette or psi-bar-psi. We measure the QCD beta function, including the anomalous dimension of the quark mass, in new Monte Carlo simulations and from results in the literature. We use it to find the interaction measure, E-3p, yielding non-perturbative values for both the energy density E and the pressure p. uuencoded compressed PostScript file Revised version should work on more PostScript printers.Comment: 24 page

Ideal evolution of MHD turbulence when imposing Taylor-Green symmetries

We investigate the ideal and incompressible magnetohydrodynamic (MHD) equations in three space dimensions for the development of potentially singular structures. The methodology consists in implementing the four-fold symmetries of the Taylor-Green vortex generalized to MHD, leading to substantial computer time and memory savings at a given resolution; we also use a re-gridding method that allows for lower-resolution runs at early times, with no loss of spectral accuracy. One magnetic configuration is examined at an equivalent resolution of $6144^3$ points, and three different configurations on grids of $4096^3$ points. At the highest resolution, two different current and vorticity sheet systems are found to collide, producing two successive accelerations in the development of small scales. At the latest time, a convergence of magnetic field lines to the location of maximum current is probably leading locally to a strong bending and directional variability of such lines. A novel analytical method, based on sharp analysis inequalities, is used to assess the validity of the finite-time singularity scenario. This method allows one to rule out spurious singularities by evaluating the rate at which the logarithmic decrement of the analyticity-strip method goes to zero. The result is that the finite-time singularity scenario cannot be ruled out, and the singularity time could be somewhere between $t=2.33$ and $t=2.70.$ More robust conclusions will require higher resolution runs and grid-point interpolation measurements of maximum current and vorticity.Comment: 18 pages, 13 figures, 2 tables; submitted to Physical Review

New measure of electron correlation

We propose to quantify the "correlation" inherent in a many-electron (or many-fermion) wavefunction by comparing it to the unique uncorrelated state that has the same single-particle density operator as it does.Comment: Final version to appear in PR