Effects of spatial size, lattice doubling and source operator on the hadron spectrum with dynamical staggered quarks

We have extended our previous study of the lattice QCD spectrum with 2 flavors of staggered dynamical quarks at $6/g^2=5.6$ and $am_q=0.025$ and 0.01 to larger lattices, with better statistics and with additional sources for the propagators. The additional sources allowed us to estimate the $\Delta$ mass and to measure the masses of all mesons whose operators are local in time. These mesons show good evidence for flavor symmetry restoration, except for the masses of the Goldstone and non-Goldstone pions. PCAC is observed in that $m_\pi^2 \propto m_q$, and $f_\pi$ is estimated. Use of undoubled lattices removes problems with the pion propagator found in our earlier work. Previously we found a large change in the nucleon mass at a quark mass of $am_q=0.01$ when we increased the spatial size from 12 to 16. No such effect is observed at the larger quark mass, $am_q=0.025$. Two kinds of wall source were used, and we have found difficulties in getting consistent results for the nucleon mass between the two sources.Comment: 30 pages PostScript fil

Hadron Spectrum in QCD with Valence Wilson Fermions and Dynamical Staggered Fermions at $6/g^2=5.6

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.Comment: 32 page

QCD thermodynamics with two flavors of Wilson quarks at N_t=6

We report on a study of hadron thermodynamics with two flavors of Wilson quarks on 12^3x6 lattices. We have studied the crossover between the high and low temperature regimes for three values of the hopping parameter, kappa=0.16, 0.17, and 0.18. At each of these values of kappa we have carried out spectrum calculations on 12^3x24 lattices for two values of the gauge coupling in the vicinity of the crossover in order to set an energy scale for our thermodynamics calculations and to determine the critical value of the gauge coupling for which the pion and quark masses vanish. For kappa=0.17 and 0.18 we find coexistence between the high and low temperature regimes over 1,000 simulation time units indicating either that the equilibration time is extremely long or that there is a possibility of a first order phase transition. The pion mass is large at the crossover values of the gauge coupling, but the crossover curve has moved closer to the critical curve along which the pion and quark masses vanish, than it was on lattices with four time slices. In addition, values of the dimensionless quantity T_c/m_rho are in closer agreement with those for staggered quarks than was the case at N_t=4. (A POSTSCRIPT VERSION OF THIS PAPER IS AVAILABLE BY ANONYMOUS FTP FROM sarek.physics.ucsb.edu (128.111.8.250) IN THE FILE pub/wilson_thermo.ps)Comment: 24 page

Properties of the a1 Meson from Lattice QCD

We determine the mass and decay constant of the $a_1$ meson using Monte Carlo simulation of lattice QCD. We find $M_{a_1} = 1250 \pm 80$ MeV and $f_{a_1} = 0.30 \pm 0.03 ~({\rm GeV})^2$, in good agreement with experiment.Comment: 9 page uu-encoded compressed postscript file. version appearing in Phys. Rev. Lett. 74 (1995) 459

Light Quark Masses from Lattice QCD

We present estimates of the masses of light quarks using lattice data. Our main results are based on a global analysis of all the published data for Wilson, Sheikholeslami-Wohlert (clover), and staggered fermions, both in the quenched approximation and with $n_f=2$ dynamical flavors. We find that the values of masses with the various formulations agree after extrapolation to the continuum limit for the $n_f=0$ theory. Our best estimates, in the MSbar scheme at $\mu=2 GeV$, are \mbar=3.4 +- 0.4 +- 0.3 MeV and $m_s = 100 +- 21 +- 10 MeV$ in the quenched approximation. The $n_f=2$ results, \mbar = 2.7 +- 0.3 +- 0.3 MeV and $m_s = 68 +- 12 +- 7 MeV$, are preliminary. (A linear extrapolation in $n_f$ would further reduce these estimates for the physical case of three dynamical flavors.) These estimates are smaller than phenomenological estimates based on sum rules, but maintain the ratios predicted by chiral perturbation theory. The new results have a significant impact on the extraction of $\epsilon'/\epsilon$ from the Standard Model. Using the same lattice data we estimate the quark condensate using the Gell-Mann-Oakes-Renner relation. Again the three formulations give consistent results after extrapolation to $a=0$, and the value turns out to be correspondingly larger, roughly preserving m_s \vev{\bar \psi \psi}.Comment: 32 pages. Package submitted in uufiles format: unpack and tex paper.tex. Modified "axis" source for figures also included. Latex2e document. Uncomment hyperref if available. This is the final published versio

Critical Dynamics of the Hybrid Monte Carlo Algorithm

We investigate the critical dynamics of the Hybrid Monte Carlo algorithm approaching the chiral limit of standard Wilson fermions. Our observations are based on time series of lengths O(5000) for a variety of observables. The lattice sizes are 16^3 x 32 and 24^3 x 40. We work at beta=5.6, and kappa=0.156, 0.157, 0.1575, 0.158, with 0.83 > m_pi/m_rho > 0.55. We find surprisingly small integrated autocorrelation times for local and extended observables. The dynamical critical exponent $z$ of the exponential autocorrelation time is compatible with 2. We estimate the total computational effort to scale between V^2 and V^2.25 towards the chiral limit.Comment: 3 pages, Latex with espcrc2.sty and postscript figures, Talk given at Lattice 9

SESAM and TXL Results for Wilson Action--A Status Report

Results from two studies of full QCD with two flavours of dynamical Wilson fermions are presented. At beta=5.6, the region 0.83 > m_pi/m_rho > 0.56 at m_pia > 0.23 L^{-1} is explored. The SESAM collaboration has generated ensembles of about 200 statistically independent configurations on a 16^3 x 32-lattice at three different kappa-values and is entering the final phase of data analysis. The TXL simulation on a 24^3 x 40-lattice at two kappa-values has reached half statistics and data analysis has started recently, hence most results presented here are preliminary. The focus of this report is fourfold: we demonstrate that algorithmic improvements like fast Krylov solvers and parallel preconditioning recently introduced can be put into practise in full QCD simulations, we present encouraging observations as to the critical dynamics of the Hybrid Monte Carlo algorithm in the approach to the chiral limit, we mention signal improvements of noisy estimator techniques for disconnected diagrams to the pi-N sigma term, and we report on SESAM's results for light hadron spectrum, light quark masses, and heavy quarkonia.Comment: 24 pages, tex + postscript figures, to appear in Proceedings of Int. Workshop "Lattice QCD on Parallel Computers", University of Tsukuba, Japa

Baryon number non-conservation and phase transitions at preheating

Certain inflation models undergo pre-heating, in which inflaton oscillations can drive parametric resonance instabilities. We discuss several phenomena stemming from such instabilities, especially in weak-scale models; generically, these involve energizing a resonant system so that it can evade tunneling by crossing barriers classically. One possibility is a spontaneous change of phase from a lower-energy vacuum state to one of higher energy, as exemplified by an asymmetric double-well potential with different masses in each well. If the lower well is in resonance with oscillations of the potential, a system can be driven resonantly to the upper well and stay there (except for tunneling) if the upper well is not resonant. Another example occurs in hybrid inflation models where the Higgs field is resonant; the Higgs oscillations can be transferred to electroweak (EW) gauge potentials, leading to rapid transitions over sphaleron barriers and consequent B+L violation. Given an appropriate CP-violating seed, we find that preheating can drive a time-varying condensate of Chern-Simons number over large spatial scales; this condensate evolves by oscillation as well as decay into modes with shorter spatial gradients, eventually ending up as a condensate of sphalerons. We study these examples numerically and to some extent analytically. The emphasis in the present paper is on the generic mechanisms, and not on specific preheating models; these will be discussed in a later paper.Comment: 10 pages, 7 figures included, revtex, epsf, references adde

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

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