BiFold: A Python code for the calculation of double-folded (bifold) potentials with density-in/dependent nucleon-nucleon interactions

Abstract

BiFold calculates the density-dependent (DDM3Yn, BDM3Yn, CDM3Yn) or independent double-folded potentials between two colliding spherical nuclei. It is written in a Python package form to give the ability to use the potentials directly in a nuclear reaction/structure code. In addition to using Woods-Saxon/Fermi or Gaussian functions, the code also allows the definition of nuclear matter densities using pre-calculated densities in a data file. The manuscript provides an overview of the double folding model and the use of the code. Program summaryProgram Title: BiFoldCPC Library link to program files:https://doi.org/10.17632/rmdx86p9dg.1Developer's repository link:https://github.com/mkarakoc/BiFoldCode Ocean capsule:https://codeocean.com/capsule/5286231Licensing provisions: GPLv3Programming language: Python 3.xNature of problem: BiFold calculates the real part of the nuclear potential between two colliding spherical nuclei by integrating a density-independent/dependent nucleon-nucleon (NN) interaction [1,2,3] over the nuclear matter densities of the two nuclei. The code is based on M3Y Reid/Paris NN interactions [1,2,3] by default, but it is possible to define custom NN interactions when necessary.Solution method: The code uses the Fourier transform method in spherical coordinates to calculate the potential. The method simplifies the sixfold integration [1] and makes the calculation significantly faster. The integration is done by default using Simpson's integration method, but Filon's integration method is also available.