Programs in Physics & Physical Chemistry
|[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.
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.
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.
With few changes one can compute the gravitational energy instead of the Coulomb energy.
Typical running time for n = m = 12 and only one deformation set is 50 s.
|Disclaimer | ScienceDirect | CPC Journal | CPC | QUB|