Investigating the effect of entrance channel mass asymmetry on fusion reactions using the Skyrme energy density formalism

Investigating the effect of entrance channel mass asymmetry on fusion reactions using the Skyrme energy density formalism 1 first-author 1 1 corresp Corresponding author, nbadiger@gmail.com Corresponding author, nbadiger@gmail.com nbadiger@gmail.com https://orcid.org/0000-0002-1633-7455 Department of Studies in Physics, Karnatak University, Dharwad 580003, India Department of Studies in Physics, Karnatak University, Dharwad 580003, India In the present investigations, the fusion cross-sections for the formation of 200 Pb compound nucleus (CN) using 16 O + 184 W, 30 Si + 170 Er, and 40 Ar + 160 Gd nuclear reactions at energies above the Coulomb barrier were calculated to understand the effect of entrance channel mass asymmetry ( α ) on the fusion reactions; the Skyrme energy density formalism (SEDF) was used for this calculation. The SEDF uses the Hartree–Fock–Bogolyubov (HFB) computational program with Skyrme forces such as SkM*, SLy4, and SLy5 to obtain the nucleus-nucleus potential parameters for the above reactions. Using the SEDF model with SkM*, SLy4, and SLy5 interaction forces, the theoretical fusion cross-sections were determined to be above the barrier energy and compared with the available experimental fusion cross-sections. The results show a close agreement between the theoretical and experimental values for all selected systems at energies well above the barrier. However, near the barrier energies, the theoretical values were observed to be higher than the experimental values. Skyrme force Energy density formalism Hartree–Fock–Bogolyubov Thomas–Fermi model Coupled-channel calculation 1. Introduction 1 The inter-nuclear potential is one of the most important factors to describe the heavy-ion nuclear reaction dynamics [ 1 1 ]. It is well known that the study of heavy-ion fusion reactions, in which the interacting nuclei effectively interact with each other, will probe the pre-saddle region as well as the field inside the Coulomb barrier extensively [ 2 2 ]. Further, the study of heavy-ion nuclear reactions can provide major information about the barrier between the interacting nuclei, which hinder s the synthesis of superheavy elements. Conversely, by studying the shape and size of the barrier distribution, the dynamics of the interacting nuclei can be investigated [ 3 3 ]. For a given system, the height and width of the barrier—in accordance with the semi-classical theory—remain same, and it is interesting to evaluate the effect of incident energy on such a barrier. Furthermore, theoretically, the reaction barrier is not very well defined and to some extent, it depends on the way it is calculated [ 4 4 , 5 5 ]. For classical simplification, we are assuming that the interaction barrier is estimated from the interaction energy of the interacting nuclei, as a function of the relative distance (R) between them. However, the dynamical calculations of nuclear reactions were also performed by solving the classical equations of motion and using the liquid drop model [ 6 6 , 7 7 ]. The estimated fusion barrier is also a deciding parameter for evaluating the probability of interaction of the heavy ions. In other words, the fusion barrier height and width are the deciding factors to evaluate the probability of fusion reactions. Therefore, the aim of the present investigations is to calculate the interaction potential using the energy density formalism (EDF) [ 8 8 ] with Skyrme forces [ 9 9 ]. Considering this , the total energy of the fermions of the system is given as a function of the one-body density [ 10 10 ]. Further, different Skyrme force gives different barrier height, position, and curvature [ 11 11 , 12 12 ]. Moreover, this force is a zero-range force and density dependent by nature. This model can be applied for inter-nuclear potential, which consistently considers the saturation property of the n uclear matter and it is known as the EDF. This model assumes that the density of each colliding ion remains the same at all distances during the collision. Using this formalism and various types of Skyrme forces, the earlier studies have shown that this method can account for the elastic scattering and fusion barrier height of 16 O + 16 O reaction [ 13 13 ]. Brink and Stancu [ 14 14 ] have investigated the applicability of this method using the Skyrme energy functional. They have shown that both the EDF and proximity potential are consistent with each other. The same work has also been performed by Dobrowolski et al. [ 15 15 ] using a higher order Thomas–Fermi approximation (TFA) for the kinetic energy and spin-orbit-densities. Recently, the EDF has been applied to simplify the coupled-channel calculations for heavy-ion nuclear fusion reactions [ 16 16 - 18 18 16-18 ]. The Skyrme energy density formalism (SEDF) is given as a sum of the spin-orbit-densities dependent universal functions; the universal function parameters are obtained for different Skyrme forces [ 19 19 ]. This method is advantageous in solving the barrier modification problems by introducing the effects of either modifying the Fermi density parameters (the half-density radii and surface thickness for exact SEDF calculations) or constants of the parameterized universal functions, as mentioned in reference [ 14 14 ]. Further, for nuclear systems, the semi-classical extended Thomas–Fermi (ETF) approach is extended at finite temperature [ 20 20 , 21 21 ] and is limited for free energy and entropy density functions only. In the semi-classical formalism, we have included the temperature effects, in the following spin-orbit density <math id="M1"><mover accent="true"><mi>J</mi><mo>→</mo></mover></math> <math id="M2"><mrow><mrow><mo>(</mo><mi>ρ</mi><mo>)</mo></mrow></mrow></math> via the nuclear density, taken as the two-parameter Fermi-density function [ 22 22 ]. This follows from the fact that the HF density approximately takes a shape of the Fermi type at higher temperatures; the density becomes flat in the interior region of smaller " r " values. Furthermore, the shell model density distributions, the Fermi density distributions, and the total density interaction potential are approximately the same, at least in the surface region of relevance for the heavy-ion collision. In actual case, the barrier height decreases with increase in temperature. However, there is a vital point of difference in temperature dependence of the two approaches, whereas the order of the barrier remains the same in the semi-classical case. The different filling of the shell model states changes the approach of the same microscopic shell models . Therefore, we have applied the EDF with Skyrme forces SkM*, SLy4 and SLy5 to obtain the fusion reaction cross-sections for the systems that form the same compound nucleus (CN) in the mass region ~ 200 amu, without adjusting any parameters in the Skyrme forces. The nuclear reactions 16 O + 184 W, 30 Si + 170 Er, and 40 Ar + 160 Gd, which were studied in the present investigations, form the same CN 200 Pb with different entrance channel mass asymmetry ( α ). The effect of mass asymmetry on the fusion reaction and the role of incompressibility ( K ) in nuclear dynamics are discussed in the following sections. The calculated fusion cross-sections are compared with the others experimental values and with cross-section obtained by the coupled-channel CCFULL code [ 23 23 ]. The CCFULL code calculates the fusion cross-section and mean-angular momentum of the CN under the influence of coupling between the relative motion and several nuclear collective motions. Using this formalism, the fusion cross-sections for the above systems in the given energy range well above the barrier were determined. 2. Theoretical background 1 It is well known that the ion-ion interaction potential can be generated using the SEDF [ 24 24 ]. Using the SEDF, it is possible to understand the fusion of two nuclei at various energies, predict the mass of nuclei as accurate as the semi-empirical mass formula, predict the ground state properties of the nuclei, and predict the fission barrier heights for all nuclei available in the periodic table [ 25 25 ]. Generally, in nuclear reactions the total potential is a combination of the nuclear, Coulomb and spin parts. However, in the SEDF, the total potential V T ( R ) is the combination of the nuclear part V N ( R ), which is the inclusion of spin-orbit interaction, and the Coulomb part V C ( R ). This total potential is given as: <math display="block" id="M3"><mrow><msub><mi>V</mi><mtext>T</mtext></msub><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow><mtext> </mtext><mo>=</mo><msub><mi>V</mi><mtext>N</mtext></msub><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow><mtext> </mtext><mo>+</mo><msub><mi>V</mi><mtext>C</mtext></msub><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow><mo>,</mo><mo> </mo><mo> </mo></mrow></math> <math display="block" id="M4"><mrow><msub><mi>V</mi><mtext>T</mtext></msub><mo stretchy="false">(</mo><mi>R</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>V</mi><mtext>B</mtext></msub><mo>−</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>μ</mi><mn>2</mn></msup><msup><mi>ω</mi><mn>2</mn></msup><msup><mrow><mo stretchy="false">(</mo><mi>R</mi><mo>−</mo><msub><mi>R</mi><mtext>B</mtext></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow></math> where equation (2) is a parabolic approximation of the total potential described by equation (1) around the barrier. The height V B and position R B of the fusion barrier can be obtained by applying the conditions [ 26 26 ]: <math display="block" id="M5"><mrow><msub><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mtext>d</mtext><msub><mi>V</mi><mtext>T</mtext></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><mi>r</mi></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mrow><mi>r</mi><mo>=</mo><msub><mi>R</mi><mtext>B</mtext></msub></mrow></msub><mo>=</mo><mn>0</mn><mtext> and </mtext><msub><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><msup><mtext>d</mtext><mn>2</mn></msup><msub><mi>V</mi><mtext>T</mtext></msub><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow><mrow><mtext>d</mtext><msup><mi>r</mi><mn>2</mn></msup></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mrow><mo stretchy="false">(</mo><mi>r</mi><mo>=</mo><msub><mi>R</mi><mtext>B</mtext></msub><mo stretchy="false">)</mo></mrow></msub><mo>≤</mo><mn>0.</mn></mrow></math> In the SEDF, V N ( R ) between the projectile and target is a function of the distance between the two nuclei. This potential is determined using the Hamiltonian H that depends on the nucleonic density ρ , kinetic energy density τ , and spin density J . It is given by <math display="block" id="M6"><mrow><msub><mi>V</mi><mtext>N</mtext></msub><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow><mtext> </mtext><mo>=</mo><msub><mi>H</mi><mrow><mn>12</mn></mrow></msub><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow><mtext> </mtext><mo>−</mo><msub><mi>H</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow><mtext> </mtext><mo>−</mo><msub><mi>H</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow><mo>,</mo></mrow></math> <math display="block" id="M7"><mrow><mo> </mo><msub><mi>H</mi><mrow><mn>12</mn></mrow></msub><mo stretchy="false">(</mo><mi>R</mi><mo stretchy="false">)</mo><mo>=</mo><mo>∫</mo><mi>H</mi><mo stretchy="false">[</mo><msub><mi>ρ</mi><mn>1</mn></msub><msub><mrow/><mi>p</mi></msub><mo stretchy="false">(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo stretchy="false">)</mo><mo>+</mo><msub><mi>ρ</mi><mn>2</mn></msub><msub><mrow/><mi>p</mi></msub><mo stretchy="false">(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo>−</mo><mover accent="true"><mi>R</mi><mo>→</mo></mover><mo stretchy="false">)</mo><mo>,</mo><msub><mi>ρ</mi><mn>1</mn></msub><msub><mrow/><mi>n</mi></msub><mo stretchy="false">(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo stretchy="false">)</mo><mo>+</mo><msub><mi>ρ</mi><mn>2</mn></msub><msub><mrow/><mi>n</mi></msub><mo stretchy="false">(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo>−</mo><mover accent="true"><mi>R</mi><mo>→</mo></mover><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mtext>d</mtext><mi>r</mi><mo>,</mo></mrow></math> <math display="block" id="M8"><mrow><msub><mi>H</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow><mo>=</mo><mtext> </mtext><mo>∫</mo><mi>H</mi><mrow><mo>[</mo><mrow><msub><mi>ρ</mi><mrow><mn>1</mn><mi>p</mi></mrow></msub><mrow><mo>(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo>)</mo></mrow><mo>,</mo><msub><mi>ρ</mi><mrow><mn>1</mn><mi>n</mi></mrow></msub><mrow><mo>(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo>)</mo></mrow></mrow><mo>]</mo></mrow><mtext>d</mtext><mi>r</mi><mo>,</mo></mrow></math> <math display="block" id="M9"><mrow><msub><mi>H</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow><mo>=</mo><mtext> </mtext><mo>∫</mo><mi>H</mi><mrow><mo>[</mo><mrow><msub><mi>ρ</mi><mrow><mn>2</mn><mi>p</mi></mrow></msub><mrow><mo>(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo>)</mo></mrow><mo>,</mo><msub><mi>ρ</mi><mrow><mn>2</mn><mi>n</mi></mrow></msub><mrow><mo>(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo>)</mo></mrow></mrow><mo>]</mo></mrow><mtext>d</mtext><mi>r</mi><mo>,</mo></mrow></math> In these equations, ρ 1n , ρ 2n are the frozen neutron densities and ρ 1p, ρ 2p are the frozen proton densities of the interacting nuclei. The H ( <math id="M10"><mover accent="true"><mi>r</mi><mo>→</mo></mover></math> <math id="M11"><mrow><mo stretchy="false">)</mo><mo> </mo></mrow></math> is associated with the Skyrme interaction and is given by: <math display="block" id="M12"><mtable columnalign="left"><mtr><mtd><mi>H</mi><mo stretchy="false">(</mo><mover accent="true"><mi>r</mi><mo>→</mo></mover><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msup><mi>ℏ</mi><mn>2</mn></msup><mi>τ</mi></mrow><mrow><mn>2</mn><mi>m</mi></mrow></mfrac><mo>+</mo><mfrac><mrow><msub><mi>t</mi><mn>0</mn></msub></mrow><mn>2</mn></mfrac><mrow><mo>[</mo><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>x</mi><mn>0</mn></msub></mrow><mn>2</mn></mfrac><msup><mi>ρ</mi><mn>2</mn></msup><mo>−</mo><mo stretchy="false">(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msubsup><mi>ρ</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>ρ</mi><mi>p</mi><mn>2</mn></msubsup><mo stretchy="false">)</mo></mrow><mo>]</mo></mrow></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><mo stretchy="false">(</mo><msub><mi>t</mi><mn>1</mn></msub><mo>+</mo><msub><mi>t</mi><mn>2</mn></msub><mo stretchy="false">)</mo><mi>ρ</mi><mi>τ</mi><mo>+</mo><mfrac><mn>1</mn><mn>8</mn></mfrac><mo stretchy="false">(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>+</mo><msub><mi>t</mi><mn>1</mn></msub><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msub><mi>ρ</mi><mi>n</mi></msub><msub><mi>τ</mi><mi>n</mi></msub><mo>+</mo><msub><mi>ρ</mi><mi>p</mi></msub><msub><mi>τ</mi><mi>p</mi></msub><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mn>1</mn><mrow><mn>16</mn></mrow></mfrac><mo stretchy="false">(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>−</mo><mn>3</mn><msub><mi>t</mi><mn>1</mn></msub><mo stretchy="false">)</mo><mi>ρ</mi><msup><mo>∇</mo><mn>2</mn></msup><mi>ρ</mi></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mn>1</mn><mrow><mn>32</mn></mrow></mfrac><mo stretchy="false">(</mo><mn>3</mn><msub><mi>t</mi><mn>1</mn></msub><mo>+</mo><msub><mi>t</mi><mn>2</mn></msub><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msub><mi>ρ</mi><mi>n</mi></msub><msup><mo>∇</mo><mn>2</mn></msup><msub><mi>ρ</mi><mi>n</mi></msub><mo>+</mo><msub><mi>ρ</mi><mi>p</mi></msub><msup><mo>∇</mo><mn>2</mn></msup><msub><mi>ρ</mi><mi>p</mi></msub><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mn>1</mn><mrow><mn>16</mn></mrow></mfrac><mo stretchy="false">(</mo><msub><mi>t</mi><mn>1</mn></msub><mo>−</mo><msub><mi>t</mi><mn>2</mn></msub><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msubsup><mover accent="true"><mi>J</mi><mo>→</mo></mover><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mover accent="true"><mi>J</mi><mo>→</mo></mover><mi>p</mi><mn>2</mn></msubsup><mo stretchy="false">)</mo><mo>+</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><msub><mi>t</mi><mn>3</mn></msub><msub><mi>ρ</mi><mi>n</mi></msub><msub><mi>ρ</mi><mi>p</mi></msub><mi>ρ</mi></mtd></mtr><mtr><mtd><mo>−</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>W</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mover accent="true"><mo>∇</mo><mo>→</mo></mover><mo>⋅</mo><mover accent="true"><mi>J</mi><mo>→</mo></mover><mo>+</mo><msub><mi>ρ</mi><mi>n</mi></msub><mover accent="true"><mo>∇</mo><mo>→</mo></mover><mo>⋅</mo><msub><mover accent="true"><mi>J</mi><mo>→</mo></mover><mi>n</mi></msub><mo>+</mo><msub><mi>ρ</mi><mi>p</mi></msub><mover accent="true"><mo>∇</mo><mo>→</mo></mover><mo>⋅</mo><msub><mover accent="true"><mi>J</mi><mo>→</mo></mover><mi>p</mi></msub><mo stretchy="false">)</mo><mo>,</mo></mtd></mtr></mtable></math> where ρ = ρ n + ρ p, τ = τ n + τ p, and <math id="M13"><mrow><mover accent="true"><mi>J</mi><mo>→</mo></mover><mo>=</mo><msub><mover accent="true"><mi>J</mi><mo>→</mo></mover><mtext>n</mtext></msub><mo>+</mo><msub><mover accent="true"><mi>J</mi><mo>→</mo></mover><mtext>p</mtext></msub></mrow></math> . The parameters ρ , τ , and J are the local densities of the nucleons, kinetic energy, and spin-orbit, respectively. The subscripts n and p represent neutron and proton, respectively. The parameters of different Skyrme forces t 0 , x 0 , t 1 , t 2 , t 3 and W 0 are estimated by fitting the different properties of nuclei [ 27 27 ]. The mass of the nucleons is indicated by m. The parameters τ and J are estimated by the semi-classical ETF model [ 28 28 ]. Further, the fusion barrier is a combination of the nuclear and coulomb potentials, as given in equation (1). Therefore, the coulomb potential between the two-interacting heavy-ions can be expressed as: <math display="block" id="M14"><mtable columnalign="left"><mtr><mtd><msub><mi>V</mi><mtext>c</mtext></msub><mrow><mo>(</mo><mrow><mi>R</mi><mo>,</mo><mi>θ</mi></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>Z</mi><mn>1</mn></msub><msub><mi>Z</mi><mn>2</mn></msub><msup><mi>e</mi><mn>2</mn></msup></mrow><mi>R</mi></mfrac><mo>+</mo><msqrt><mrow><mfrac><mn>9</mn><mrow><mn>20</mn><mi>π</mi></mrow></mfrac></mrow></msqrt><mo> </mo><mfrac><mrow><msub><mi>Z</mi><mn>1</mn></msub><msub><mi>Z</mi><mn>2</mn></msub><msup><mi>e</mi><mn>2</mn></msup></mrow><mrow><msup><mi>R</mi><mn>3</mn></msup></mrow></mfrac><mo> </mo><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><msubsup><mi>R</mi><mi>i</mi><mn>2</mn></msubsup><msub><mi>β</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><msub><mi>P</mi><mn>2</mn></msub><mrow><mo>(</mo><mrow><mi>cos</mi><msub><mi>θ</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>+</mo><mo> </mo><mfrac><mn>3</mn><mrow><mn>7</mn><mi>π</mi></mrow></mfrac><mo> </mo><mfrac><mrow><msub><mi>Z</mi><mn>1</mn></msub><msub><mi>Z</mi><mn>2</mn></msub><msup><mi>e</mi><mn>2</mn></msup></mrow><mrow><msup><mi>R</mi><mn>3</mn></msup></mrow></mfrac><mo> </mo><mo> </mo><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><msubsup><mi>R</mi><mi>i</mi><mn>2</mn></msubsup><msup><mrow><mo>[</mo><mrow><msub><mi>β</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><msub><mi>P</mi><mn>2</mn></msub><mrow><mo>(</mo><mrow><mi>cos</mi><msub><mi>θ</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow><mn>2</mn></msup><mo>,</mo></mtd></mtr></mtable></math> where <math id="M15"><mrow><mo> </mo><msub><mi>θ</mi><mi>i</mi></msub></mrow></math> , R i , and β 2 i are the angle between radius vector and symmetry axis, effective sharp radius, and the quadruple deformation of ith nucleus, respectively. <math id="M16"><mrow><msub><mi>R</mi><mi>i</mi></msub><mo>=</mo><mn>1.28</mn><mo> </mo><msubsup><mi>A</mi><mi>i</mi><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo>−</mo><mn>0.76</mn><mo>+</mo><mn>0.8</mn><mo> </mo><msubsup><mi>A</mi><mi>i</mi><mrow><mo>−</mo><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup></mrow></math> , as given by Blocki [ 29 29 , 30 30 ]. During collisions, the nuclei may have several orientations. Therefore, it is preferred to consider the average over all possible relative orientations from 0 o to 180 o for both deformed nuclei in the study of fusion reactions [ 31 31 ]. Since both the nuclei are in the same plane i.e. <math id="M17"><mrow><msub><mi>P</mi><mn>2</mn></msub><mrow><mo>(</mo><mrow><mi>cos</mi><msub><mi>θ</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><mrow><mn>3</mn><mi>cos</mi><msub><mi>θ</mi><mi>i</mi></msub><msup><mrow/><mn>2</mn></msup><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></math> , the average possible orientation is unity [ 32 32 ]. In SEDF, the nucleon density distribution is determined using the ion-ion interaction potential. The density distribution of the interacting nuclei can be obtained by different approximation methods. However, the self-consistent Hartree–Fock–Bogolyubov (HFB) method can precisely determine the Skyrme effective interaction using a complete quantum approach [ 33 33 ]. The density distribution of the nucleons in the nuclei is calculated by the code HFBRAD(v1.00), which does not include the deformation of the nucleus [ 33 33 ]. <math display="block" id="M18"><mrow><msub><mi>ρ</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>ρ</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub></mrow><mrow><mrow><mo>[</mo><mrow><mn>1</mn><mo>+</mo><mtext>exp</mtext><mrow><mo>(</mo><mrow><mfrac><mrow><mi>r</mi><mo>−</mo><msub><mi>R</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub></mrow><mrow><msub><mi>a</mi><mi>i</mi></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow></mfrac><mo> </mo></mrow></math> The central density, <math id="M19"><mrow><msub><mi>ρ</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub></mrow></math> , is given by <math display="block" id="M20"><mrow><msub><mi>ρ</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn>3</mn><msub><mi>A</mi><mi>i</mi></msub></mrow><mrow><mn>4</mn><mi>π</mi><msubsup><mi>R</mi><mrow><mn>0</mn><mi>i</mi></mrow><mn>3</mn></msubsup></mrow></mfrac><msup><mrow><mrow><mo>[</mo><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><msup><mi>π</mi><mn>2</mn></msup><msubsup><mi>a</mi><mi>i</mi><mn>2</mn></msubsup></mrow><mrow><msubsup><mi>R</mi><mrow><mn>0</mn><mi>i</mi></mrow><mn>2</mn></msubsup></mrow></mfrac></mrow><mo>]</mo></mrow></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow></math> The half density radii <math id="M21"><mrow><msub><mi>R</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub><mtext> </mtext></mrow></math> and the surface thickness <math id="M22"><mrow><msub><mi>a</mi><mi>i</mi></msub></mrow></math> are determined by applying a polynomial fit to the experimentally determined values of the nuclei in the mass range from 4 to 238 amu [ 21 21 , 34 34 ]. 3. Result and discussions 1 Heavy ion-induced fusion reaction is a sensitive tool to probe the shape, structure, and size of an atomic nucleus as well as the reaction dynamics of a di-nuclear system (DNS) [ 34 34 ]. In the fusion of heavy ions, the relative kinetic energy of an entrance channel converts into the internal excitation energy of the CN and this internal excitation energy is distributed among all the nucleons in the nucleus. The fusion reactions can be explained based on the energy density distribution of the interacting nucleons in the nuclei . Therefore, in the present work, calculations of the SEDF were carried out using energy density distributions. Systems, such as 16 O + 184 W, 30 Si + 170 Er, and 40 Ar + 160 Gd were used, which form the same CN, 200 Pb. The ion-ion potential parameters are determined using compatible Skyrme forces like SkM*, SLy4, and SLy5. It is important to notice that in the present investigations, the reactions with different entrance channels and α were selected deliberately to form the same CN. The α are different with respect to the Businaro–Gallone mass asymmetry ( α BG ) [ 35 35 ]. Some fitting parameters are required for the corresponding Skyrme forces to calculate the ion-ion potential parameters using the HFB method. These parameters are tabulated in Table 1 Table 1 with their respective forces. The Woods–Saxon potential parameters, calculated using the Akyuz–Winther (AW) formalism, along with α and the deformation parameters for the corresponding systems are given in Table 2. The local ρ , τ , and J density distributions of the nuclei for all the above reactions were determined using the HFB program. A typical local densities distribution of 16 O + 184 W for both neutron and proton are shown in Figs. 1 Figs. 1 and 2 2 . The Skyrme energy density function H ( <math id="M23"><mover accent="true"><mi>r</mi><mo>→</mo></mover></math> ), as given in Eq. (8), was evaluated using local densities, which are obtained from the Skyrme forces of the HFB program. The interaction nuclear potential V N ( R ), as given in Eq. (4), was determined using this H ( <math id="M24"><mover accent="true"><mtext>r</mtext><mo>→</mo></mover></math> ). The nuclear potential for all the Skyrme forces for the selected reaction systems is shown in Figure 3 Figure 3 . In this figure, we notice that the depth of the potential depends on the K of nuclear matter and the entrance channel α of the reaction; K is a significant component of the equation of state (EOS) of nuclear matter and it is one of the interesting subjects in the study of heavy-ion fusion reactions. For the investigation of the fusion process, different K values were used, which result from the corresponding effective interactions [ 36 36 , 37 37 ]. The results obtained in our calculations show that the theoretical values are sensitive to the values of K . Therefore, it is very interesting to describe the nuclear fusion reaction using different effective interactions with different K and α values. Consequently, by knowing the K value, we can understand the variations in K of nuclear matter for various energies of projectile. 1 1 Parameters 1 1 SkM* 1 1 SLy4 1 1 SLy5 1 1 t 0 1 1 -2645.0 1 1 -2488.913 1 1 -2488.354 1 1 t 1 1 1 410.0 1 1 486.818 1 1 484.230 1 1 t 2 1 1 -135.0 1 1 -546.395 1 1 -556.690 1 1 t 3 1 1 15595.0 1 1 13777.0 1 1 13757.0 1 1 x 0 1 1 0.09 1 1 0.834 1 1 0.776 1 1 x 1 1 1 0.0 1 1 -0.344 1 1 -0.317 1 1 x 2 1 1 0.0 1 1 -1.000 1 1 -1.000 1 1 x 3 1 1 0.0 1 1 1.354 1 1 1.263 1 1 γ 1 1 1/6 1 1 1/6 1 1 1/6 1 1 W 0 1 1 120.0 1 1 123.0 1 1 125.0 For the calculation of EDF based on each Skyrme force, we have estimated the Fermi density distributions using the parameters obtained from the fitting results of the HFB. The calculated density distributions’ diffuseness parameters of the nucleons for the projectiles and targets of the systems 16 O + 184 W, 30 Si + 170 Er, and 40 Ar + 160 Gd are shown in Figs. 1 Figs. 1 and 2 2 using Skyrme forces SkM*, SLy4 and SLy5. The interaction potentials were obtained for above mentioned reactions using the calculated densities with the respective Skyrmes. The properties of a calculated fusion barrier potential, viz. V 0 , a 0 , and R B , are given in Table 3 Table 3 for the corresponding Skyrmes. The results clearly indicate that V 0 increases with the increasing K value. In other words, the increase in the value of K increases the depth of reaction potential pocket. This observed increase in the potential depth may be due to a vital role of the peripheral nucleons in the heavy-ion reaction. Another set of potential parameters for the above systems was calculated using the AW formalism for the CCFULL calculation ( Table 2 Table 2 ) [ 38 38 ]. It is also to be noted that the system 16 O + 184 W has α ˃ α BG, and the systems 30 Si + 170 Er and 40 Ar + 160 Gd have α ˂ α BG . In case of incompressibility, SkM* [ 39 39 ] has lesser incompressibility as compared to SLy4 and SLy5 [ 40 40 ]. The nuclear potential parameters obtained from different Skyrme forces and with different incompressibility of the HFB are given in Table 3 Table 3 . The Coulomb part of the nuclear reaction was also estimated using equation (9) for all the systems. The total potential for the three Skyrme forces were determined for all the selected systems using the nuclear and coulomb potentials. There exists a potential pocket that is responsible for the capture of the nuclei. The deformation of the interacting nuclei is neglected here to avoid the long duration of the calculation [ 41 41 , 42 42 ]. With respect to the α, the system 40 Ar + 160 Gd may have less capture probability; therefore, the fission barrier height appears to be very less compared to the remaining systems, and hence, the contribution of the quasi-fission may lead to less fusion probability. To confirm this, our group is planning to experimentally measure these type of systems in the above mentioned energy range using a recoil mass separator. By incorporating these potential parameters—derived from the different Skyrme forces—in the CCFULL code as an input, we have calculated the fusion cross-section for reactions 16 O + 184 W, 30 Si + 170 Er, and 40 Ar + 160 Gd. The fusion cross-sections of the CN are calculated using the following equation: 1 1 Skyrme forces 1 1 Incompressibility K (MeV) 1 1 Systems 1 1 V 0 (MeV) 1 1 a 0 (fm) 1 1 R B (fm) 1 1 SkM* 1 1 216.65 1 1 16O+184 W 1 1 49.50 1 1 0.41 1 1 6.9 1 1 1 1 1 1 30Si+170Er 1 1 48.00 1 1 0.40 1 1 7.2 1 1 1 1 1 1 40Ar+160Gd 1 1 46.00 1 1 0.41 1 1 7.1 1 1 SLy4 1 1 229.90 1 1 16O+184 W 1 1 69.50 1 1 0.459 1 1 7.3 1 1 1 1 1 1 30Si+170Er 1 1 66.80 1 1 0.45 1 1 7.2 1 1 1 1 1 1 40Ar+160Gd 1 1 64.10 1 1 0.45 1 1 7.19 1 1 SLy5 1 1 229.90 1 1 16O+184 W 1 1 70.20 1 1 0.45 1 1 7.3 1 1 1 1 1 1 30Si+170Er 1 1 67.10 1 1 0.455 1 1 7.2 1 1 1 1 1 1 40Ar+160Gd 1 1 64.80 1 1 0.46 1 1 7.19 1 1 Reaction Systems 1 1 Bass Barrier (MeV) 2 1 Mass asymmetry 2 1 Deformation 3 1 Potential parameters using AW 1 1 1 1 1 1 α 1 1 α BG 1 1 β 2p 1 1 β 2T 1 1 V 0 (MeV) 1 1 R 0 (fm) 1 1 a 0 (fm) 1 1 16 O+ 184 W 1 1 74.19 1 1 0.84 1 1 0.83 1 1 0.364 1 1 0.236 1 1 63.90 1 1 1.20 1 1 0.65 1 1 30 Si+ 170 Er 1 1 113.24 1 1 0.70 1 1 0.83 1 1 0.315 1 1 0.336 1 1 72.93 1 1 1.20 1 1 0.67 1 1 40 Ar+ 160 Gd 1 1 100.43 1 1 0.60 1 1 0.83 1 1 0.251 1 1 0.353 1 1 77.57 1 1 1.20 1 1 0.68 <math display="block" id="M25"><mrow><msub><mtext>σ</mtext><mrow><mtext>fus</mtext></mrow></msub><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mtext>π</mtext><mrow><msup><mi>k</mi><mn>2</mn></msup></mrow></mfrac><munder><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mi>l</mi></munder><mrow><mo>(</mo><mrow><mn>2</mn><mi>l</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo></mrow><msub><mi>P</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow></mrow></math> where <math id="M26"><mrow><msub><mi>P</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow></mrow></math> is the transmission probability or penetrability for the l -wave scattering and k is the local wave number for the l th partial wave [ 23 23 ]. It has been used to explain the fusion cross-section in CCFULL at energies near and above the barrier [ 38 38 ]. The theoretical fusion cross-section at various excitation energies have been compared with the available experimental values [ 43 43 , 44 44 ] and the theoretical fusion cross-sections obtained by the CCFULL code with Woods–Saxon potential parameters of the AW formalism as the inputs. Figure 4 Figure 4 gives the comparison of fusion cross-sections. As shown in this figure, the density distributions of the heavy ions, which were calculated by different methods and using different parameters of Skyrmes, influence the fusion barrier, which are calculated using the EDF. With this motivation, we have carried out a comparative study of the fusion barrier using SkM*, SLy4, and SLy5 forces, in addition to the Fermi density distribution including the parameters determined in the SEDF. The total fusion cross-sections were computed from the estimated potentials and compared with the available experimental data. By varying the K values with different α of the same CN, a good agreement was obtained between the calculated fusion cross-sections derived from different Skyrme forces and the available experimental fusion cross-sections at different bombarding energies . It was also observed that as the K value increases, the depth of the reaction capture pocket increases. In other words, it increases the probability of the fusion and level density . Similarly, with the increase in α, the fusion probability would also increase. It implies that α between the projectile and the target nuclei is responsible for the deformation of the DNS during the dynamic process. If α value is large, the DNS will be less deformed and consequently, the contribution of the fission and fission-like probabilities will be less. Conversely, if α value is small, the DNS will be more deformed, and therefore, the contribution of the fission or fission-like probability will be more. Further, K also depends on the nucleons lying on the surface or on the periphery of the nuclei, which are more actively involved during the interaction, even though the system is less asymmetric. In other words, K would enhance the probability of the capture cross-section even for less asymmetric systems. 4. Conclusion 1 In the present investigations, different Skyrme forces were used to estimate the fusion cross-sections at energies well above the Coulomb barrier for the reactions 16 O + 184 W, 30 Si + 170 Er, and 40 Ar + 160 Gd. In this study, three systems were selected such that the first reaction has α ≥ α BG , the second one has α ˂ α BG , and the third one has α ˂˂ α BG . The results of the SEDF calculations show that the calculated values are in close agreement with the experimental values of 16 O + 184 W system, indicating that whenever α is close to α BG the fusion cross-section values, based on Skyrme forces, closely agree with the experimental values. On the other hand, in the case of 30 Si + 170 Er ( α ˂ α BG ), the calculated fusion cross-section based on Skyrme forces were found to be higher than the experimental values. It is interesting to note that as the energy increases the deviation between the theoretical and experimental values decreases. In the case of 40 Ar + 160 Gd, there are no experimental data on fusion cross-sections. Therefore, we have compared the calculated fusion cross-sections based on Skyrme forces with the values obtained from the CCFULL code. Interestingly, the calculated theoretical values closely agree with the experimental values in the higher energy region, but in the lower energy region, the values predicted by CCFULL are higher than the theoretical values. Therefore, we may conclude that the fusion cross-section values, obtained with different Skyrme forces, can be used not only to estimate the fusion cross-section for any projectile and target combination, but it can also throw light on the effect of entrance channel α on the fusion cross-section . Further, there is no need to change the parameters of Skyrme force to modify the barrier as the theoretical framework based on the Skyrme forces can reproduce the corresponding experiment values. For the fusion cross-sections at sub-barrier energies, different Skyrme forces are required for the different reactions to reproduce the experimental data. A single Skyrme force cannot reproduce the fusion cross-sections over the entire energy range for the selected reactions. References W. Kohn , L. J. Sham , Quantum density oscillations in an inhomogeneous electron gas . Phys. Rev . 137 , A1697 ( 1965 ). https://doi.org/10.1103/PhysRev.137.A1697 https://doi.org/10.1103/PhysRev.137.A1697 Quantum density oscillations in an inhomogeneous electron gas P. Hohenberg , W. Kohn , Inhomogeneous electron gas . Phys. Rev . 136 , B864 ( 1964 ). https://doi.org/10.1103/PhysRev.136.B864 https://doi.org/10.1103/PhysRev.136.B864 Inhomogeneous electron gas R. Kumar , M. K. Sharma , R. K. Gupta , Fusion reaction cross-sections using the Wong model within Skyrme energy density based semiclassical extended Thomas Fermi approach . Nucl. Phys. A 870 - 871 , 42 - 57 ( 2011 ). https://doi.org/10.1016/j.nuclphysa.2011.09.010 https://doi.org/10.1016/j.nuclphysa.2011.09.010 Fusion reaction cross-sections using the Wong model within Skyrme energy density based semiclassical extended Thomas Fermi approach R. K. Gupta , Dalip Singh ,  Raj Kumar et al ., Universal functions of nuclear proximity potential for Skyrme nucleus-nucleus interaction in a semiclassical approach . J. Phys. G: Nucl. Part. Phys . 36 , 075104 ( 2009 ). Universal functions of nuclear proximity potential for Skyrme nucleus-nucleus interaction in a semiclassical approach K. A. Brueckner , J. R. Buchler , R. C. Clark et al ., Phys. Rev . 181 , 1534 ( 1969 ). Phys. Rev H. Holm , W. Greiner , Influence of nuclear forces on the Coulomb barrier in heavy-ion reactions . Phys. Rev. Lett . 24 , 404 ( 1970 ). https://doi.org/10.1103/PhysRevLett.24.404 https://doi.org/10.1103/PhysRevLett.24.404 Influence of nuclear forces on the Coulomb barrier in heavy-ion reactions A. S. Jensen , C. Y. Wong , Distortion of nuclei in heavy ion reactions . Phys. Lett. B 32 , 567 - 570 ( 1970 ). https://doi.org/10.1016/0370-2693(70)90545-9 https://doi.org/10.1016/0370-2693(70)90545-9 Distortion of nuclei in heavy ion reactions H. Holm , D. Scharnweber , W. Scheid et al ., Coulomb distortion in heavy-ion reactions . Z. Phys . 231 , 450 ( 1970 ). https://doi.org/10.1007/BF01642535 https://doi.org/10.1007/BF01642535 Coulomb distortion in heavy-ion reactions P. W. Riesenfeld , T. D. Thomas , Effect of nuclear deformability on reaction cross sections . Phys. Rev. C 2 , 711 ( 1970 ). https://doi.org/10.1103/PhysRevC.2.711 https://doi.org/10.1103/PhysRevC.2.711 Effect of nuclear deformability on reaction cross sections L. Wilets . A. Goldberg , F. H. Lewis , New theoretical approach to nuclear heavy-ion scattering . Phys. Rev . 2 , 1576 ( 1970 ). https://doi.org/10.1103/PhysRevC.2.1576 https://doi.org/10.1103/PhysRevC.2.1576 New theoretical approach to nuclear heavy-ion scattering K. A. Brueckner , J. R. Buchler , M. M. Kelly , New theoretical approach to nuclear heavy-ion scattering . Phys. Rev . 173 , 944 ( 1968 ). https://doi.org/10.1103/PhysRev.173.944 https://doi.org/10.1103/PhysRev.173.944 New theoretical approach to nuclear heavy-ion scattering T. H. R. Skyrme et al ., Phil . Mag . 1 , 1043 ( 1956 ); Phil T. H. R. Skyrme et al ., Phil . Nucl. Phys. A 9 , 615 ( 1959 ). Phil C. Ngow , B. Tamain , J. Galin et al , Calculation of interaction barriers using the energy density formalism . Nucl. Phys. A 240 , 353 ( 1975 ). https://doi.org/10.1016/0375-9474(75)90335-8 https://doi.org/10.1016/0375-9474(75)90335-8 Calculation of interaction barriers using the energy density formalism D. M. Brink , FL. Stancu , Interaction potential between two 16O nuclei derived from the Skyrme interaction . Nucl. Phys. A 243 , 175 - 188 ( 1975 ). https://doi.org/10.1016/0375-9474(75)90027-5 https://doi.org/10.1016/0375-9474(75)90027-5 Interaction potential between two 16O nuclei derived from the Skyrme interaction A. Dobrowolski , K. Pomorski , J. Bartel , Mean-field description of fusion barriers with Skyrme`s interaction . Nucl. Phys. A 729 , 713 ( 2003 ). https://doi.org/10.1016/j.nuclphysa.2003.09.008 https://doi.org/10.1016/j.nuclphysa.2003.09.008 Mean-field description of fusion barriers with Skyrme`s interaction L. Kaur , R. Kumari , On the role of deformed Coulomb potential in fusion using energy density formalism . Pramana-Journal of Physics 85 , 639 ( 2015 ). https://doi.org/10.1007/s12043-014-0898-z https://doi.org/10.1007/s12043-014-0898-z On the role of deformed Coulomb potential in fusion using energy density formalism O. N. Ghodsi , F. Torabi , Comparative study of fusion barriers using Skyrme interactions and the energy density functional . Phys. Rev. C 92 , 064612 ( 2015 ). https://doi.org/10.1103/PhysRevC.92.064612 https://doi.org/10.1103/PhysRevC.92.064612 Comparative study of fusion barriers using Skyrme interactions and the energy density functional O. N. Ghodsi , F. Torabi , Effect of nuclear matter incompressibility on the 16O + 208Pb system . Phys. Rev. C 93 , 064611 ( 2016 ). https://doi.org/10.1103/PhysRevC.93.064611 https://doi.org/10.1103/PhysRevC.93.064611 Effect of nuclear matter incompressibility on the 16O + 208Pb system R. J. Lombard , The energy density formalism in nuclei . Annals of Physics 77 , 380 - 413 ( 1973 ). https://doi.org/10.1016/0003-4916(73)90422-3 https://doi.org/10.1016/0003-4916(73)90422-3 The energy density formalism in nuclei M. Brack , C. Guet , H. B. Hakansson , Selfconsistent semiclassical description of average nuclear properties - A link between microscopic and macroscopic models . Phys. Rep . 123 , 275 ( 1985 ). Selfconsistent semiclassical description of average nuclear properties - A link between microscopic and macroscopic models J. Bartel , M. Brack , M. Durand , Extended Thomas-Fermi theory at finite temperature . Nucl. Phys. A 445 , 263 - 303 ( 1985 ). https://doi.org/10.1016/0375-9474(85)90071-5 https://doi.org/10.1016/0375-9474(85)90071-5 Extended Thomas-Fermi theory at finite temperature D. Singh et al ., Proceedings of DAE-BRNS symposium on Nuclear Physics , Vol B46 , P- 254 ( 2003 ). Proceedings of DAE-BRNS symposium on Nuclear Physics K. Hagino , N. Rowley , A. T. Kruppa , A program for coupled-channel calculations with all order couplings for heavy-ion fusion reactions . Comp. Phys. Comm . 123 , 143 - 152 ( 1999 ). https://doi.org/10.1016/S0010-4655(99)00243-X https://doi.org/10.1016/S0010-4655(99)00243-X A program for coupled-channel calculations with all order couplings for heavy-ion fusion reactions A. S. Umar , V. E. Oberacker , Three-dimensional unrestricted time-dependent Hartree-Fock fusion calculations using the full Skyrme interaction . Phys. Rev. C 73 , 054607 ( 2006 ). https://doi.org/10.1103/PhysRevC.73.054607 https://doi.org/10.1103/PhysRevC.73.054607 Three-dimensional unrestricted time-dependent Hartree-Fock fusion calculations using the full Skyrme interaction J. W. Negele , The mean-field theory of nuclear structure and dynamics . Rev. Mod. Phys . 54 , 913 ( 1982 ). https://doi.org/10.1103/RevModPhys.54.913 https://doi.org/10.1103/RevModPhys.54.913 The mean-field theory of nuclear structure and dynamics K. T. R. Davies et al ., in Treatise on heavy ion science , edited by D. A. Bromley ( Plenum , New York , 1985 ), Vol. 03 P.0 3 . in Treatise on heavy ion science J. P. Svenne , Fourteen years of self-consistent field calculations: What has been learned . Adv. Nucl. Phys . 11 , 179 ( 1979 ). Fourteen years of self-consistent field calculations: What has been learned R. K. Gupta , D. Singh , W. Greiner , Semiclassical and microscopic calculations of the spin-orbit density part of the Skyrme nucleus-nucleus interaction potential with temperature effects included . Phys. Rev. C 75 , 024603 ( 2007 ). https://doi.org/10.1103/PhysRevC.75.024603 https://doi.org/10.1103/PhysRevC.75.024603 Semiclassical and microscopic calculations of the spin-orbit density part of the Skyrme nucleus-nucleus interaction potential with temperature effects included J. Blocki , J. Randrup , W. J. Swiatecki et al ., Proximity forces . Ann. Phys. (NY) 105 , 427 ( 1977 ). https://doi.org/10.1016/0003-4916(77)90249-4 https://doi.org/10.1016/0003-4916(77)90249-4 Proximity forces R. K. Puri , R. K. Gupta , Analytical formulation of the ion-ion interaction potential including spin density term in energy density formalism . Int. J. Mod. Phys. E 1 , 269 ( 1992 ). https://doi.org/10.1142/S0218301392000138 https://doi.org/10.1142/S0218301392000138 Analytical formulation of the ion-ion interaction potential including spin density term in energy density formalism G. B. Arfken et al , Mathematical Methods for Physicists: A comprehensive guide , 7 th edition, Academic Press is an imprint of Elsevier 225 Wyman Street , Waltham, MA 02451, USA . Mathematical Methods for Physicists: A comprehensive guide Ishita Sharma , R. Kumar , M. K. Sharma , Probing the role of Skyrme interactions on the fission dynamics of the 6Li + 238U reaction . Eur. Phys. J. A 53 , 140 ( 2017 ). https://doi.org/10.1140/epja/i2017-12331-5 https://doi.org/10.1140/epja/i2017-12331-5 Probing the role of Skyrme interactions on the fission dynamics of the 6Li + 238U reaction K. Bennaceur , Coordinate-space solution of the Skyrme-Hartree-Fock-Bogolyubov equations within spherical symmetry . The program HFBRAD (v1.00). Comp. Phys. Comm. 168 , 96 ( 2005 ). https://doi.org/10.1016/j.cpc.2005.02.002 https://doi.org/10.1016/j.cpc.2005.02.002 Coordinate-space solution of the Skyrme-Hartree-Fock-Bogolyubov equations within spherical symmetry A. Shamlath , E. Prasad , N. Madhavan et al ., Fusion and quasifission studies in reactions forming Rn via evaporation residue measurements . Phys. Rev. C 95 , 034610 ( 2017 ). https://doi.org/10.1103/PhysRevC.95.034610 https://doi.org/10.1103/PhysRevC.95.034610 Fusion and quasifission studies in reactions forming Rn via evaporation residue measurements B. B. Back , R. R. Betts , J. E. Gindler et al ., Angular distribution in heavy-ion induced fission . Phys. Rev. C 32 , 195 ( 1985 ). https://doi.org/10.1103/PhysRevC.32.195 https://doi.org/10.1103/PhysRevC.32.195 Angular distribution in heavy-ion induced fission S. Misicu , H. Esbensen , Signature of shallow potentials in deep sub-barrier fusion reactions . Phys. Rev. C 75 , 034606 ( 2007 ). https://doi.org/10.1103/PhysRevC.75.034606 https://doi.org/10.1103/PhysRevC.75.034606 Signature of shallow potentials in deep sub-barrier fusion reactions H. Esbensen , S. Misicu , Hindrance of 16O + 208Pb fusion at extreme sub-barrier energies . Phys. Rev. C 76 , 054609 ( 2007 ). https://doi.org/10.1103/PhysRevC.76.054609 https://doi.org/10.1103/PhysRevC.76.054609 Hindrance of 16O + 208Pb fusion at extreme sub-barrier energies K. Hagino , N. Takigawa , Subbarrier fusion reactions and many-particle quantum tunneling . Progress of Theoretical Physics , 128 , 1061 - 1106 ( 2012 ). https://doi.org/10.1143/PTP.128.1061 https://doi.org/10.1143/PTP.128.1061 Subbarrier fusion reactions and many-particle quantum tunneling J. Bartel , P. Quentin , Towards a better parameterization of Skyrme-like effective forces: A critical study of the SkM force . Nucl. Phys. A 386 , 79 ( 1982 ). Towards a better parameterization of Skyrme-like effective forces: A critical study of the SkM force E. Chabanat , P. Bonche , P. Haensel et al ., A Skyrme parameterization from subnuclear to neutron star densities part II . Nuclei far from stability. Nucl. Phys. A 635 , 231 ( 1998 ). https://doi.org/10.1016/S0375-9474(98)00180-8 https://doi.org/10.1016/S0375-9474(98)00180-8 A Skyrme parameterization from subnuclear to neutron star densities part II J. Dobaczewski , J. Dudek , Solution of the Skyrme-Hartree-Fock equations in the Cartesian deformed harmonic oscillator basis I . The method. Comput. Phys. Comm . 102 , 166 ( 1997 ). Solution of the Skyrme-Hartree-Fock equations in the Cartesian deformed harmonic oscillator basis I J. Dobaczewski , J. Dudek , Solution of the Skyrme-Hartree-Fock equations in the Cartesian deformed harmonic oscillator basis I . The program HFODD. Comput. Phys. Comm . 102 , 183 ( 1997 ). Solution of the Skyrme-Hartree-Fock equations in the Cartesian deformed harmonic oscillator basis I P. D. Shidling , N. M. Badiger , S. Nath et al ., Fission hindrance studies in 200Pb: Evaporation residue cross section and spin distribution measurements . Phys. Rev. C 74 , 064603 ( 2006 ). https://doi.org/10.1103/PhysRevC.74.064603 https://doi.org/10.1103/PhysRevC.74.064603 Fission hindrance studies in 200Pb: Evaporation residue cross section and spin distribution measurements G. Mohanto , N. Madhavan , S. Nath et al ., Entrance channel effect on ER spin distribution . Nucl. Phys. A 890 , 62 ( 2012 ). https://doi.org/10.1016/j.nuclphysa.2012.07.004 https://doi.org/10.1016/j.nuclphysa.2012.07.004 Entrance channel effect on ER spin distribution 10.1007/s41365-020-00802-w 2020-03-18 2020-07-22 2020-07-27 2020-08-29 2020-09-01 2021-04-10 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2020 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-020-00802-w http://dx.doi.org/10.1007/s41365-020-00802-w . 2020 Nuclear Science and Techniques 9 31