About CCFULLYPE
CCFULLYPE is a modified version of the code CCFULL for deep sub-barrier fusion reactions based on an adiabatic approach. The original CCFULL code is written by K. Hagino. In comparison to the original CCFULL, the Yukawa-plus-exponential (YPE) model is adopted as a basic heavy ion-ion potential and the damping factor which describes a smooth transition from sudden to adiabatic processes is implemented in the code. In this version, the adiabatic one-body potential is approximated with a third-order polynomial function. You can easily handle the code only by adding a few lines in the input file of the original CCFULL.
Note that CCFULLYPE, different from the original CCFULL, does not take into account the full order of the couplings. All definitions follow Esbensen's one to compare the sudden model, namely, the coupling potential is expanded up to the second order of amplitude.
You can download the CCFULYPE code from here.
Examples of input files for 64Ni + 64Ni and 16O + 208Pb (ccfull.inp and input) are also contained.
Input Parameters
Lines from 1 to 8 are the same as those of the original CCFULL.
Line 9: RMAX, DR, ISW
RMAX and DR are the same as those of the original CCFULL.
If ISW = 0, a calculation is identical to that of the original CCFULL.
If ISW = 1 (recommended), in a calculation, the position of the incoming wave boundary condition (IWBC) is searched at each angular momentum. Thus, the speed of a calculation becomes slow, but its calculated result is stabilized. If you see strange fluctuations in calculated fusion cross section, use ISW = 1.
Line 10: IYPE, R0YPE, A0YPE, EGS, HOMG
If IYPE = 0, the Woods-Saxon (WS) potential is used and the code is identical to the original CCFULL. The parameters of the WS potential are taken from the line 7.
If IYPE = 1, the YPE potential is used. Then, EGS is ignored and is taken from the liquid-drop model.
If IYPE = 2 (recommended), the YPE potential is used and the value of EGS has to be inputted by a user.
R0YPE and A0YPE are the radius and diffuseness parameters of the YPE potential.
EGS is the energy at the ground state. If IYPE = 1, EGS is calculated with the liquid-drop model. If IYPE = 2, EGS is taken from your input value. EGS should be estimated by the reaction Q value calculated using experimental mass tables (e.g. Audi's table).
HOMG is the curvature of the potential energy at the ground state. If HOMG = 0.0 (recommended), the value of HOMG is automatically determined by a systematic curve.
Line 11: IDAMP, R0DAMP, A0DAMP
if IDMAP = 0, the damping factor is ignored and the code is identical to the original CCFULL.
if IDAMP = 1, the damping factor is switched on.
R0DAMP and A0DAMP are the damping radius and diffuseness parameters.
The details of the YPE potential and the damping factor are described in [T. Ichikawa, Phy. Rev. C92, 064604 (2015)]
Hints for Fitting Data
First, you should start to fit data with IDAMP = 0 by adjusting R0YPE and coupling strengths. In many cases, you do not need to vary A0YPE. After you obtain enough enhancements of fusion cross section at extremely low incident energies, the damping factor should be switched on by IDAMP = 1. After that, you should adjust R0DAMP and A0DAMP so as to reproduce the fusion hindrance.