A C-code for the double folding interaction potential of two spherical nuclei

Abstract

We present a C-code designed to obtain the nucleus–nucleus potential by using the double folding model (DFM) and in particular to find the Coulomb barrier. The program calculates the nucleus–nucleus potential as a function of the distance between the centers of mass of colliding nuclei. The most important output parameters are the Coulomb barrier energy and the radius. Since many researchers use a Woods–Saxon profile for the nuclear term of the potential we provide an option in our code for fitting the DFM potential by such a profile. Program summary Program title: DFMSPH Catalogue identifier: AEFH_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEFH_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 5929 No. of bytes in distributed program, including test data, etc.: 115 740 Distribution format: tar.gz Programming language: C Computer: PC Operating system: Windows XP (with the GCC-compiler version 2) RAM: Below 10 Mbyte Classification: 17.9 Nature of problem: The code calculates in a semimicroscopic way the bare interaction potential between two colliding spherical nuclei as a function of the center of mass distance. The height and the position of the Coulomb barrier are found. The calculated potential is approximated by a conventional Woods–Saxon profile near the barrier. Dependence of the barrier parameters upon the characteristics of the effective NN forces (like, e.g. the range of the exchange part of the nuclear term) can be investigated. Solution method: The nucleus–nucleus potential is calculated using the double folding model with the Coulomb and the effective M3Y NN interactions. For the direct parts of the Coulomb and the nuclear terms, the Fourier transform method is used. In order to calculate the exchange parts the density matrix expansion method is applied. Running time: Less than 1 minute using a PC with a 1.60 GHz processor.