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] aabf_v1_0.gz(12 Kbytes)
Manuscript Title: Monte Carlo simulation of pure SU(2) gauge theory on a body-centered hypercubic lattice.
Authors: W. Celmaster, K.J.M. Moriarty
Program title: BCH
Catalogue identifier: AABF_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 34(1985)415
Programming language: Fortran.
Computer: DEC VAX 11/750.
Operating system: DEC VAX 11/750: VMS 3.6.
RAM: 128K words
Word size: 32
Keywords: Particle physics, Elementary, Quantum chromodynamics, Lattice gauge theory, Su(2) gauge theory, Body-centred hypercubic Lattice, Wilson loops, Monte carlo techniques, Vector processors, Triangular lattice.
Classification: 11.5.

Nature of problem:
The usual hypercubic formulation of lattice gauge theory may be improved by the alternative formulation of the theory of a body- centered hypercubic (BCH) lattice. The method for implementing the Monte Carlo algorithm for calculating Wilson loops and lines in an SU(2) gauge theory on a BCH lattice is described.

Solution method:
SU(2) gauge theory is simulated by importance-sampling Monte Carlo methods on a BCH lattice. The modified heat-bath method of Creutz is used to thermalize the lattice.

Restrictions:
(a) In order to correctly compute the Wilson lines, the lattice must have an even length.
(b) The main memory allocation in the program is for the 12 arrays ALATI,...,ALAT4, BLAT1,...,BLAT4, CLAT1,...,CLAT4 each dimensioned to 12*ISIZE**4 and 3 arrays MSHIFT, MBSHIFT and MSCHIFT dimensioned to 24* ISIZE**4 where ISIZE is the lattice length. Since the VAX is a virtual memory machine the only real limitation on program size is given by the available virtual memory. In practical terms, if the program far exceeds the physical memory (1 mgeword on the test machine) then most of the CPU time is spent doing disc 1/O operations.

Running time:
The test run took 44 mins. on the VAX 11/750 with ISIZE = 6 for 40 sweeps of the lattice and 4 sets of measurements.