Computer Physics Communications Program LibraryPrograms 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_0Distribution 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. |

Disclaimer | ScienceDirect | CPC Journal | CPC | QUB |