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] acva_v1_0.gz(13 Kbytes)
Manuscript Title: Self-consistent solution to a fermion determinant with space dependent fields.
Authors: R. Alkofer, H. Weigel
Program title: chisol
Catalogue identifier: ACVA_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 82(1994)30
Programming language: Fortran.
Computer: SUN-SPARC.
Operating system: UNIX, MS-DOS.
RAM: 572K words
Peripherals: disc.
Keywords: Particle physics, Elementary, Empirical model, Nambu-jona-lasinio model, Chiral soliton, Hedgehog ansatz, Functional determinant, Static energy functional.
Classification: 11.6.

Nature of problem:
Quark fields are integrated out from the Nambu-Jona-Lasinio model yielding a determinant which involves composite meson fields. The fermion determinant is evaluated for static but space-dependent meson field configurations. The configuration which minimizes the corresponding energy functional is constructed self-consistently.

Solution method:
A basis for the quark spinors is provided suited for spherically symmetric configurations. The one-particle Dirac Hamiltonian involving the space-dependent meson field is diagonalized using this basis. The fermion determinant is calculated by summing over the resulting eigenvalues. The corresponding eigenfunctions serve to set up the equation of motion for the meson field. The solution to this equation is obtained by iteration.

Only one dynamical field is considered. This is generalized in an accompanying paper [1].

Running time:
The running time strongly depends on the size of the momentum space basis as well as the number of lattice points. Typical parameters yield for a single iteration 240s (SUN), 275s (DEC), 95s (HP), 1h (PC; 386/387 20 MHz) cpu-time.

[1] (Same as ref 10 in long write-up) (If needed, storage may be reduced by choosing a smaller value for the parameter NMAX e.g. for NMAX=80 (the minimal value needed for the test run) the binary size is 480KB.)