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] adgz_v1_0.gz(39 Kbytes)
Manuscript Title: PRESCOLD: calculation of the half-life for alpha decay, cluster radioactivity and cold fission processes.
Authors: M. Goncalves, S.B. Duarte, F. Garcia, O. Rodriguez
Program title: PRESCOLD
Catalogue identifier: ADGZ_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 107(1997)246
Programming language: Fortran.
Computer: Intel 80386+80387.
Operating system: MS-DOS 6.00, Windows 95, Windows NT 4.0.
RAM: 9M words
Word size: 16
Keywords: Nuclear physics, Heavy ion, Inertia coefficient, Gamow penetrability, Half-life calculation, Alpha decay, Cluster radioactivity, Cold fission.
Classification: 17.7.

Nature of problem:
Spontaneous cold fission processes have appeared as a new source of information about the nuclear properties in the ground state [1,2]. Theoretical attempts in respect to the unification of alpha decay, cluster radioactivity and cold fission phenomena have appeared in the literature recently [3-5]. The inertia coefficient as well as the shape parametrization used to describe the dynamical evolution of the separating dinuclear system are of key importance for the half-life calculation. The present numerical code calculates the half-lives of alpha decay, cluster radioactivity and cold fission processes by using four alternative inertial coefficients.

Solution method:
Geometric and incompressibility constraint relations are used explicitly in reducing the number of collective variables in the molecular phase of the fissioning system to calculate the Gamow penetrability factor. Varying Mass Asymmetry Shape (VMAS) and Constant Mass Asymmetry Shape (CMAS) modes are used to calculate the inertial coefficient in the framework of "effective mass coefficient" [6] and Werner-Wheeler approximation [7]. The fundamental Q-value energy calculated in terms of the mass excesses taken from the last version of the Nuclear Mass Data Table [8] is regarded in the calculation.

All possible structural information of the fragmenting nucleus in the calculation are present only in the Q-values. The multidimensional potential surface is reduced to an onedimensional problem. The effective liquid drop model description for potential barrier is used.

Running time:
Depends on the choice for calculations. To compute one particular decay of one parent nucleus the running time is approximately 2 seconds in a typical Pentium 150MHz PC.

[1] A. Sandulescu et al., Phys. Rev. C54 (1996) 258.
[2] S. Singh et al., J. Phys. G18 (1992) 1243.
[3] S. Sham et al., Phys. Lett. B399 (1997) 8.
[4] W. Greiner, Lectures on the VI Andre Swieca Brazilian Summer School, 1997, Sao Paulo, Brazil; W. Greiner, Heavy Ion Physics 2 (1995) 23.
[5] E. Stefanescu et al., Phys. Rev. C53 (1996) 3014.
[6] S.B. Duarte and M. Goncalves, Phys. Rev. C53 (1996) 2309.
[7] D.N. Poenaru, J.A. Maruhn, W. Greiner, M. Ivascu, D. Mazilu, I. Ivascu, Z. Phys. A-Atomic Nuclei 333 (1989) 291.
[8] G. Audi and A.H. Wapstra, Nucl, Phys. A595 (1995) 409.