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

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

Applications of spectral methods to turbulent magnetofluids in space and fusion research

Recent and potential applications of spectral method computation to incompressible, dissipative magnetohydrodynamics are surveyed. Linear stability problems for one dimensional, quasi-equilibria are approachable through a close analogue of the Orr-Sommerfeld equation. It is likely that for Reynolds-like numbers above certain as-yet-undetermined thresholds, all magnetofluids are turbulent. Four recent effects in MHD turbulence are remarked upon, as they have displayed themselves in spectral method computations: (1) inverse cascades; (2) small-scale intermittent dissipative structures; (3) selective decays of ideal global invariants relative to each other; and (4) anisotropy induced by a mean dc magnetic field. Two more conjectured applications are suggested. All the turbulent processes discussed are sometimes involved in current carrying confined fusion magnetoplasmas and in space plasmas

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

Two-Flavor Staggered Fermion Thermodynamics at N_t = 12

We present results of an ongoing study of the nature of the high temperature crossover in QCD with two light fermion flavors. These results are obtained with the conventional staggered fermion action at the smallest lattice spacing to date---approximately 0.1 fm. Of particular interest are a study of the temperature of the crossover a determination of the induced baryon charge and baryon susceptibility, the scalar susceptibility, and the chiral order parameter, used to test models of critical behavior associated with chiral symmetry restoration. From our new data and published results for N_t = 4, 6, and 8, we determine the QCD magnetic equation of state from the chiral order parameter using O(4) and mean field critical exponents and compare it with the corresponding equation of state obtained from an O(4) spin model and mean field theory. We also present a scaling analysis of the Polyakov loop, suggesting a temperature dependent ``constituent quark free energy.''Comment: LaTeX 25 pages, 15 Postscript figure

Thermodynamics of lattice QCD with two light quark flavours on a 16^3 x 8 lattice II

We have extended our earlier simulations of the high temperature behaviour of lattice QCD with two light flavours of staggered quarks on a $16^3 \times 8$ lattice to lower quark mass (m_q=0.00625). The transition from hadronic matter to a quark-gluon plasma is observed at $6/g^2=5.49(2)$ corresponding to a temperature of $T_c \approx 140$MeV. We present measurements of observables which probe the nature of the quark-gluon plasma and serve to distinguish it from hadronic matter. Although the transition is quite abrupt, we have seen no indications that it is first order.Comment: 23 pages, RevteX, 6 encapsulated postscript figure

A precise determination of T_c in QCD from scaling

Existing lattice data on the QCD phase transition are analyzed in renormalized perturbation theory. In quenched QCD it is found that T_c scales for lattices with only 3 time slices, and that T_c/Lambda_msbar=1.15 \pm 0.05. A preliminary estimate in QCD with two flavours of dynamical quarks shows that this ratio depends on the quark mass. For realistic quark masses we estimate T_c/Lambda_msbar=0.49 \pm 0.02. We also investigate the equation of state in quenched QCD at 1-loop order in renormalised perturbation theory.Comment: 7 pages, 5 eps figures; improved error analysis yields smaller errors on T_

Testing Algorithms for Finite Temperature Lattice QCD

We discuss recent algorithmic improvements in simulating finite temperature QCD on a lattice. In particular, the Rational Hybrid Monte Carlo(RHMC) algorithm is employed to generate lattice configurations for 2+1 flavor QCD. Unlike the Hybrid R Algorithm, RHMC is reversible, admitting a Metropolis accept/reject step that eliminates the $\mathcal{O}(\delta t^2)$ errors inherent in the R Algorithm. We also employ several algorithmic speed-ups, including multiple time scales, the use of a more efficient numerical integrator, and Hasenbusch pre-conditioning of the fermion force.Comment: 4 pages, 2 figures, poster presented at International Conference on Strong and Electroweak Matter 2006 (SEWM2006), BNL, May 10-13, 200