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] acfg_v1_0.gz(9 Kbytes)
Manuscript Title: Monte Carlo simulation of pure U(N) and SU(N) gauge theories on a simplicial lattice.
Authors: J.-M. Drouffe, K.J.M. Moriarty, C.N. Mouhas
Program title: LATTICE
Catalogue identifier: ACFG_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 30(1983)249
Programming language: Fortran.
Computer: VAX 11/780.
Operating system: VMS 3.2.
RAM: 110K words
Word size: 32
Keywords: Elementary, Particle physics, Qcd, Lattice gauge theory, Yang-mills theory, U(n) and su(n) gauge Theories, Simplicial lattice, Wilson loops, Monte carlo methods.
Classification: 11.5.

Nature of problem:
Lattice gauge theory is usually formulated on a hypercubical lattice. However, in order to check universality it is necessary to formulate the theory on another lattice, e.g. the simplicial lattice. The method for mplementing the Monte Carlo algorithm for pure U(N) and SU(N) gauge theories on the recently introduced simplicial lattice is presented. The calculational technique for calculating all Wilson loops containing up to 13 triangles or 9 squares is explained.

Solution method:
Pure U(N) and SU(N) gauge theories are simulated by Monte Carlo methods on a four-dimensional simplicial space-time lattice. The method of Metropolis et al. is used to equilibriate our lattice.

Restrictions:
The only restriction on the program is the size of the links array which is dimensioned to
10(ISIZE)**4*2N**2,
where ISIZE is the number of lattices sites in any space-time dimension and N**2 is the number of elements in an U(N) or SU(N) matrix. This effectively sets the limit to the size of lattice which will reside in real memory. Of course, the VAX 11/780 is a virtual memory machine and we could use this facility to examine larger lattices.

Running time:
The test run took 6 h and 17 min on the VAX 11/780 with ISIZE=6 and N=2 for 42 sweeps through the lattice.