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] abqh_v1_0.gz(18 Kbytes)
Manuscript Title: Folded Yukawa-plus-exponential model PES for nuclei with different charge densities.
Authors: D.N. Poenaru, M. Ivascu, D. Mazilu
Program title: FYPEDIFC
Catalogue identifier: ABQH_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 19(1980)205
Programming language: Fortran.
Computer: IBM 370/135.
Operating system: DOS/VS-310.
RAM: 118K words
Word size: 8
Keywords: Nuclear physics, Liquid drop model, Folded yukawa-plus Expotential, Charge density asymmetry, Potential energy surface, Theoretical methods, Optical model.
Classification: 17.9, 17.16.

Nature of problem:
The deformation dependent terms of the potential energy in the folded Yukawa-plus-exponential model for nuclei with different charge densities are: volume, Coulomb and surface energies. They are expressed by an analytical relationship, two- and three-fold integrals, respectively, if one assumes nuclear shapes with axial symmetry. One has to consider separately the self-energies of the two fragments and their interaction energy.

Solution method:
For a given mass and charge asymmetry, the multiple integrals are computed numerically by Gauss-Legendre quadrature with 6 X 6 grid points, which are used to obtain by interpolation an array of 51 X 51 elements and to plot the potential-energy-surface (PES) map with the line printer. An iteration method is used by the saddle-point searching computer code.

Restrictions:
The 4m X 4m and 4m X 4m X n mesh points for numerical quadrature are chosen to assure the required accuracy for a given shape and charge density asymmetry. Only values m,n <= 64 are allowed.

Unusual features:
With few changes one can compute the gravitational energy instead of the Coulomb energy.

Running time:
Typical running time for n = m = 12 and only one deformation set is 50 s.