Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] abdh_v1_0.gz(28 Kbytes)
Manuscript Title: A vectorized code for the pseudofermion simulation of QCD with dynamical quarks.
Authors: M. Campostrini, K.J.M. Moriarty, J. Potvin, C. Rebbi
Program title: PSEUDOF
Catalogue identifier: ABDH_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 50(1988)395
Programming language: Fortran.
Computer: CDC CYBER 205, MODEL 682.
Operating system: CDC CYBER 200 VSOS 2.1.6.
RAM: 3500K words
Word size: 64
Keywords: Particle physics, Elementary, Quantum chromodynamics, Hadron masses, Pseudofermion, Metropolis and heat Bath algorithms.
Classification: 11.5.

Nature of problem:
Monte Carlo simulations of Quantum Chromodynamics (QCD) often neglect the effects of virtual quark loops. These effects can be included using the pseudofermion algorithm, giving a complete solution of QCD and a more accurate description of hadronic phenomena.

Solution method:
The exact simulation of dynamical quark effects is impractical even on lattices of moderate size. The pseudofermion algorithm is an efficient approximation scheme that can be tuned to any desired accuracy. For every Monte Carlo sweep, the contribution of virtual qq- pairs is computed with a Monte Carlo simulation over auxiliary pseudofermionic fields, and added to the standard bosonic action to upgrade the gauge fields. The two Monte Carlo processes and the evaluation of the Wilson loops and of the quark propagator are optimized using techniques similar to those of Comp. Phys. Commun. 32(1984)1, 36(1985)1 and 42(1986)175.

(a) The space and time lattice size NS and NT must be even.
(b) With half-word precision arithmetic, 110 * NS**3 * NT words of fast memory are required.