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

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 summaryProgram title: DFMSPHCatalogue identifier: AEFH_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFH_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 5929No. of bytes in distributed program, including test data, etc.: 115 740Distribution format: tar.gzProgramming language: CComputer: PCOperating system: Windows XP (with the GCC-compiler version 2)RAM: Below 10 MbyteClassification: 17.9Nature 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.