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

Abstract

This is a new version of the DFMSPH code designed to obtain the nucleus–nucleus potential by using the double folding model (DFM) and in particular to find the Coulomb barrier. The new version uses the charge, proton, and neutron density distributions provided by the user. Also we added an option for fitting the DFM potential by the Gross–Kalinowski profile. The main functionalities of the original code (e.g. the nucleus–nucleus potential as a function of the distance between the centers of mass of colliding nuclei, the Coulomb barrier characteristics, etc.) have not been modified. New version program summaryProgram title: DFMSPH14Catalog identifier: AEFH_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFH_v2_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland.Licensing provisions: GNU General Public License, version 3No. of lines in distributed program, including test data, etc.: 7211No. of bytes in distributed program, including test data, etc.: 114404Distribution format: tar.gzProgramming language: CComputer: PC and MacOperation system: Windows XP and higher, MacOS, Unix/LinuxMemory required to execute with typical data: below 10 MbyteClassification: 17.9Catalog identifier of previous version: AEFH_v1_0Journal reference of previous version: Comp. Phys. Comm. 181 (2010) 168Does the new version supersede the previous version?: YesNature of physical 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 an analytical profile (Woods–Saxon or Gross–Kalinowski) 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.Method of solution: 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.Typical running time: less than 1 minute.Reason for new version: Many users asked us how to implement their own density distributions in the DFMSPH. Now this option has been added. Also we found that the calculated Double-Folding Potential (DFP) is approximated more accurately by the Gross–Kalinowski (GK) profile. This option has been also added.