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] abqg_v1_0.gz(6 Kbytes)
Manuscript Title: Liquid drop model deformation energies of nuclei with axial symmetry and reflection asymmetry.
Authors: D.N. Poenaru, M. Ivascu
Catalogue identifier: ABQG_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 16(1978)85
Programming language: Fortran.
Computer: IBM 370/135.
Operating system: DOS/VS-310.
RAM: 33K words
Word size: 8
Keywords: Low energy, Nuclear physics, Structure, Liquid drop model, Finite range, Collective model, Collision, Krappe-nix integral, Coulomb energy, Surface energy.
Classification: 17.20.

Nature of problem:
The deformation dependent terms of the nuclear potential energy in the Myers-Swiatecki's liquid drop model are: Coulomb and surface energies. For nuclear shapes with axial symmetry, they are expressed by two fold and simple integrals respectively. The surface energy was redefined by Krappe and Nix in order to take into account the finite range of nuclear forces. Calculation of this energy requires a three fold integration.

Solution method:
The multiple integrals involved are computed numerically by Gauss- Legendre quadrature. Because the nuclear shape is usually described by a different equation in the neck region, the accuracy and speed is improved by a suitable choice of four subintervals of integration.

The number of meshpoints for numerical quadrature are chosen to assure the required accuracy for a given shape. Only the following values are allowed: n=4, 8, 12 and 16 for each of the four subintervals.

Unusual features:
The program as it stands is used for nuclear energies. But with few changes it may be used for calculating the electrostatic energy and the surface area of any body with axial symmetry and uniform distribution of matter and electric charge, assuming that the charge and matter density is a single-valued function of Z.

Running time:
Typical running time on the IBM 370/135 for n=8 and one deformation set is 12 s.