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] abqs_v1_0.gz(152 Kbytes)
Manuscript Title: NUDENS: a Nilsson-Bardeen-Cooper-Schrieffer code at finite nuclear temperature.
Authors: G. Maino, M. Vaccari, A. Ventura
Program title: NUDENS
Catalogue identifier: ABQS_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 29(1983)375
Programming language: Fortran.
Computer: IBM 370/168.
Operating system: MVS/SP1.
RAM: 26K words
Word size: 32
Peripherals: disc.
Keywords: Nuclear physics, Grand canonical Formalism, Nuclear shell model, Pairing interaction, Nuclear temperature, Entropy, Level densities, Collective model.
Classification: 17.20.

Nature of problem:
Microscopic formalism (NBCS approximation) for the calculation of nuclear thermodynamic functions, especially level densities versus excitation energy and angular momentum. A blocking effect due to unpaired nucleons in odd-A, or odd-odd nuclei is explicitly worked out.

Solution method:
A system of nonlinear equations is solved at zero and finite nuclear temperature in order to determine pairing correlation functions and chemical potentials, to be introduced in grand canonical traces, from which thermodynamic functions are derived. Excitation energy, total angular momentum and ground state correlation functions are taken as input parameters.

Restrictions:
Neutron-proton pairing correlations are neglected: therefore, the formalism should be applied to medium and heavy nuclei ( A>=50). Collective enhancements of nuclear level densities are introduced in a simplified phenomenological form. The equlibrium shape of the nucleus is not allowed to change with excitation energy.

Unusual features:
NUDENS uses the mathematical subroutine NS01AD from the Harwell Subroutine Library, in order to solve a system of n nonlinear equations in n unknowns. This subroutine is included in the deck.

Running time:
The execution time is strongly problem dependent. The test run requires 1 min 20 s on OH-5480 and 1 min 38 s on IBM 370/168: in this case nuclear level densities are computed at three different excitation energies.