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] adik_v1_0.tar.gz(265 Kbytes)
Manuscript Title: Partial level densities for nuclear data calculations.
Authors: M. Avrigeanu, V. Avrigeanu
Program title: PLD
Catalogue identifier: ADIK_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 112(1998)191
Programming language: Fortran.
Computer: PC 486.
Operating system: DOS.
RAM: 491K words
Keywords: Nuclear physics, Spectr, Level scheme, Partial nuclear level Density, Nuclear level density, Single-particle level Density, Equidistant-spacing Model, Preequilibrium emission, Reactions nuclear.
Classification: 17.15.

Nature of problem:
This Fortran code is a collection of subroutines for calculation of the partial nuclear level densities (PLD) mainly used in preequilibrium nuclear reaction models, by using 14 formalisms for the related partial state densities (PSD).

Solution method:
The main approaches to the calculation of the partial state density, based on the equidistant spacing model (ESM), are used. Composite (recommended) formulas including optionally various corrections, i.e. the advanced pairing and shell correction in addition to the Pauli effect, and average energy-dependent single-particle level (s.p.l.) densities for the excited particles and holes, are also involved. The density of the particle-hole bound states is moreover comprised, and the effects of an exact correction for the Pauli-exclusion principle are considered.

Although the quantum-mechanical s.p.l. density and the continuum effect can also be reproduced by a corresponding Fermi-gas formula, to be used accordingly within the average energy-dependent PSD in multistep reaction models, this effect is not included. The calculation of PLD with linear momentum, of first interest for modelling preequilibrium- emission angular distributions, is not available either.

Unusual features:
The PSD functions have been optimized for their independent use, in order to provide tools for PSD/PLD users (see [5]). The related drawback is the increase in execution time, while a proper use would involve the calculation of some coefficients only once in the main program. Second, the PLD.FOR has been organized so that various formulas and versions may be tried as well as the comparison between their predictions.

Running time:
The execution time is strongly problem-dependent: it is roughly proportional to both the number of the excitation energies and the exciton configurations considered in the calculation, while consistent differences arise when various PSD formalisms are used. Twenty-eight sample cases with 73 versions which require from 0.1 to 1661 s on a PC Pentium/166MHz are provided.

[1] M. Avrigeanu and V. Avrigeanu, J. Phys. G 20 (1994) 613.
[2] M. Avrigeanu, A. Harangozo, and V. Avrigeanu, Rom. J. Phys. 41 (1996) 77.
[3] V. Avrigeanu, A. Florescu, A. Sandulescu, and W. Greiner, Phys. Rev. C 52 (1995) 1765.
[4] M. Avrigeanu, A. Harangozo, V. Avrigeanu, and A.N. Antonov, Phys. Rev. C 54 (1996) 2538; ibid. 56 (1997) 1633.
[5] M. Avrigeanu, I. Stetcu, and V. Avrigeanu, Development of data file with partial level densities for nuclear data calculations. Final report for IAEA Research Contract No. 8886, 15.12.1995-15.12.1997, http://tndln1.ifa.ro/~vavrig/absf8886.html.