B3Y-FETAL effective interaction in the folding analysis of elastic scattering of ¹⁶O + ¹⁶O

B3Y-FETAL effective interaction in the folding analysis of elastic scattering of 16 O + 16 O 1 corresp first-author I. Ochala, ojoagbamiisaiah1234@gmail.com I. Ochala, ojoagbamiisaiah1234@gmail.com ojoagbamiisaiah1234@gmail.com 2 Department of Physics, Kogi State University, P.M.B 1008, Anyigba, Nigeria Department of Physics, Kogi State University, P.M.B 1008, Anyigba, Nigeria Department of Physics, Benue State University, Makurdi, Nigeria Department of Physics, Benue State University, Makurdi, Nigeria In this paper, a new M3Y-type effective nucleon-nucleon interaction,derived based on the lowest order constrained variational approach (LOCV) and termed B3Y-Fetal, has been used in DDM3Y1, BDM3Y1, BDM3Y2, and BDM3Y3 density-dependent versions in a heavy ion (HI) optical potential based on four types of a real folded potential and a phenomenological Woods-Saxon imaginary potential to study the elastic scattering of the 16 O+ 16 O nuclear system within the framework of the optical model (OM) by computing the associated differential cross-sections at various incident energies. The results of the folding analyses have shown the DDB3Y1-Fetal and BDB3Y1-Fetal, out of the four folded potentials, give a reasonably better description of the elastic data of the nuclear system. These best-fit folded potentials are followed, in performance, by the BDB3Y2-Fetal, with the BDB3Y3-Fetal potential coming last. This performance trend was also demonstrated by the optical potentials based on the M3Y-Reid interaction. Furthermore, the best-fit folded potentials, renormalized by a factor N R of approximately 0.9, have been shown to reproduce the energy dependence of the real optical potential for 16 O scattering found in previous optical model analyses creditably well. In excellent agreement with previous works, they have also been identified in this work to belong to the family of deep refractive potentials because they have been able to reproduce and consistently describe the evolution of Airy-like structures, at large scattering angles, observed in the 16 O scattering data at different energies. Finally, a comparison of the performances of B3Y-Fetal and M3Y-Reid effective interactions undertaken in this work has shown impressive agreement between them. B3Y-fetal Folding analysis Elastic scattering Optical model Incident energies 1 Introduction 1 The folding model is a microscopic optical potential that has been developed to meet the need for an optical model (OM) to give nucleon-nucleus and nucleus-nucleus scattering [ 1 1 ]. For years, it has been used to calculate the elastic potential for heavy ions and the inelastic form factor [ 2 2 ]. It is believed to be the first-order term of the optical potential derived from Feshbach’s theory of nuclear reactions [ 3 3 , 4 4 ], and it generates the real parts of both nucleon-nucleus and nucleus-nucleus potentials [ 5 5 , 6 6 ]. Its remarkable success in describing the observed elastic scattering in many systems shows that it produces the dominant part of the real optical potential. In this work, we use a newer approach in OM analysis, in which the optical potential is derived in the framework of macroscopic approximation using nuclear matter properties, nuclear surface, density distributions of nuclei, and an effective interaction. Here, the nuclear matter properties are incorporated in the OM analysis through the use of the double folding model [ 3 3 ]. In this approach, the optical potential is generated by folding an effective nucleon-nucleon (NN) interaction over the ground-state density distributions of the two nuclei [ 1 1 , 7 7 ]. Of all effective interactions, the M3Y-type interaction, which is based on a realistic NN force, is a popular choice for use in the double folding model [ 8 8 , 9 9 ]. Consequently, different density-dependent versions of the designed M3Y-Reid and M3Y-Paris interactions within the Hartree-Fock (HF) formalism that generate different values of the incompressibility, K, of nuclear matter, have been severally used either in the single folding [ 4 4 ] or double folding models [ 2 2 , 7 7 ]. In the same vein, the new M3Y-type effective interaction, derived by Fiase and co-researchers [ 10 10 ] and based on the lowest order constrained variational approach (LOCV), called B3Y-Fetal [ 11 11 ] for the purpose of identification, has been successfully applied to nuclear matter calculations with values of incompressibility K that have been found to demonstrate excellent agreement with previous work with the M3Y-Reid and M3Y-Paris interactions. The new effective interaction is used herein in its DDM3Y and BDM3Y density-dependent versions to perform a detailed analysis of the refractive nucleus-nucleus scattering of a nuclear system. The B3Y-Fetal is used alongside the M3Y-Reid interaction in the DDM3Y1, BDM3Y1, BDM3Y2, and BDM3Y3 versions for the folding analysis following the approach of Khoa et al. [ 7 7 , 8 8 ] with the intent to determine and compare the performance of the first with the second interaction. To reproduce the empirical energy dependence of the nucleon-nucleus optical potential, a weak intrinsic energy dependence is introduced into the new effective interaction. Because the energy dependence of the real optical potential comes mostly from the exchange term [ 7 7 ], the study of the elastic scattering of chosen light heavy ions is carried out in this work at various collision energies of interest using the exact finite-range evaluation of the exchange potential by means of local approximation within the context of the double folding model. This approach has been used with immense success by many researchers [ 12 12 ]. Because the nuclear system considered here is spin-saturated, only the isoscalar component of the new M3Y-type effective interaction is chiefly employed in all optical model analyses. The folded potential arising from the double folding model is useful only for the real part of the optical potential, and it is supplemented by a phenomenological imaginary part [ 13 13 ]. The choice of the nuclear system for this study is based on reasons related directly to its peculiarities, one of which is that it exhibits weak absorption in elastic scattering so that refractive scattering occurs in it. The refractive scattering in this system is known to provide important information on the equation of state (EOS) of nuclear matter [ 2 2 , 7 7 ] because the double folding model, with an effective NN interaction appropriately adjusted to nuclear matter properties, enters the optical model analyses as a mean-field effect that produces a very deep real component of the optical potential for the elastic channel. In this way, analysis of the chosen system provides a unique opportunity to gain clear insight into the nature of the EOS of normal nuclear matter. Furthermore, refractive scattering in this system plays a great role in determination of the real component of the heavy-ion (HI) optical potential down to small inter-nuclear distances, giving the possibility of testing various theoretical models [ 2 2 , 14 14 ] for HI optical potential. Because B3Y-Fetal is an effective interaction based on a theoretical model, different from the theoretical models for other effective interactions that are popularly used for constructing HI optical potentials, it is strongly believed that the study of the elastic scattering of the chosen light heavy ions, using a HI optical potential derived from the new effective interaction, affords ample opportunity to test and affirm the theoretical model used for developing the effective interaction. In view of the foregoing, we determined to apply the newly developed B3Y-Fetal interaction to the folding analysis of the chosen light heavy ions within the framework of the OM in the present work, in anticipation that the expected results will provide insight into those qualities that will help determine the true place of the new effective interaction among the numerous theoretical models for HI optical potential. It is also hoped that the data arising from the nuclear reactions will give a clear picture of the EOS of nuclear matter, the nuclear reaction dynamics of the nuclear systems involved, and the relative performance of the effective interactions used. In this light, the remainder of this paper is organized as follows. Section 2 provides an overview of the density- and energy-dependent B3Y-Fetal interactions. In Sect. 3, explanation of the basic features of the double folding model employed for the study of nuclear reactions is provided. In Sect. 4, we present the results of our calculated cross-sections and compare them with experimental data, while Sect. 5 focusses on summary and discussion. 2 Density- and Energy-Dependent B3Y-Fetal Interaction 1 The B3Y-Fetal interaction is a variational effective interaction derived from the lowest-order constrained variational approach [ 10 10 ]. The matrix elements were calculated using a harmonic oscillator basis, and the details of the calculation were reported in [ 10 10 ] and in our earlier report, in which we applied it to symmetric nuclear matter (SNM) calculations [ 11 11 ]. In this study, we use the isoscalar part of the B3Y-Fetal effective interaction [ 10 10 ], whose Yukawa strengths (in MeV) are given in Table 1 Table 1 along with those of the M3Y-Reid [ 15 15 ]: 1 1 Interaction 1 1 i 1 1 1/ μ i 1 1 Y D 1 1 Y EX 1 1   1 1   1 1 (fm) 1 1 (MeV) 1 1 (MeV) 1 1 B3Y-Fetal 1 1 1 1 1 0.25 1 1 10472.13 1 1 499.63 1 1   1 1 2 1 1 0.40 1 1 -2203.11 1 1 -1347.77 1 1   1 1 3 1 1 1.414 1 1 0.0 1 1 -7.8474 1 1 M3Y-Reid 1 1 1 1 1 0.25 1 1 7999.00 1 1 4631.375 1 1   1 1 2 1 1 0.40 1 1 -2134.25 1 1 -1787.125 1 1   1 1 3 1 1 1.414 1 1 0.0 1 1 -7.8474 Here, the radial shape of the isoscalar M3Y-type effective interaction is expressed in terms of three Yukawas as [ 15 15 ]: <math display="block" id="M1"><mrow><msup><mi>v</mi><mrow><mtext>D(EX)</mtext></mrow></msup><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover></mstyle><msup><mi>Y</mi><mrow><mtext>D(EX)</mtext></mrow></msup><mo stretchy="false">(</mo><mi>i</mi><mo stretchy="false">)</mo><mfrac><mrow><mi>exp</mi><mo stretchy="false">(</mo><mo>−</mo><msub><mi>μ</mi><mi>i</mi></msub><mi>r</mi><mo stretchy="false">)</mo></mrow><mrow><msub><mi>μ</mi><mi>i</mi></msub><mi>r</mi></mrow></mfrac><mn>,</mn></mrow></math> where the functions Y D(EX) are represented in terms of the singlet even (SE), singlet odd (SO), triplet even (TE), and triplet odd (TO) angular momenta channels [ 15 15 ]. To be able to reproduce the saturation properties of cold symmetric nuclear matter and satisfy the requirement for folding analysis, the B3Y-Fetal interaction is used in the energy- and density-dependent form [ 12 12 , 16 16 ]: <math display="block" id="M2"><mrow><msup><mi>v</mi><mrow><mtext>D(EX)</mtext></mrow></msup><mo stretchy="false">(</mo><mi>E</mi><mn>,</mn><mi>ρ</mi><mn>,</mn><mi>r</mi><mo stretchy="false">)</mo><mo>=</mo><mi>g</mi><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mi>h</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><msup><mi>v</mi><mrow><mtext>D(EX)</mtext></mrow></msup><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo></mrow></math> where h ( ρ ) is the density-dependent factor and g ( E ) <math id="M3"><mo>≃</mo></math> 1-0.002E for M3Y-Reid [ 7 7 , 8 8 ] as well as B3Y-Fetal; E is the incident nucleon energy. This isoscalar energy- and density-dependent interaction has been tested and found to work well in folding model analyses of refractive α -nucleus [ 8 8 ] and nucleus-nucleus [ 8 8 ] elastic scattering. The explicit forms of h ( ρ ) used in this study are [ 7 7 , 8 8 ]: <math display="block" id="M4"><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><mi>C</mi><mo stretchy="false">(</mo><mn>1</mn><mo>+</mo><mi>y</mi><msup><mi>e</mi><mrow><mo>−</mo><mi>z</mi><mi>ρ</mi></mrow></msup><mo stretchy="false">)</mo></mrow></math> for the DDM3Yn ( n = 1) interaction and <math display="block" id="M5"><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><mi>C</mi><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>y</mi><msup><mi>ρ</mi><mi>z</mi></msup><mo stretchy="false">)</mo></mrow></math> for the BDM3Y n ( n = 0, 1, 2, 3) interactions. The parameters C, y , and z of the density dependences are such that they reproduce the saturation properties of nuclear matter at density ρ 0 = 0.17 fm -3 , with a binding energy of E/A = 16 MeV within HF calculations. To ascertain the functional viability of the B3Y-Fetal interaction, we applied it to symmetric nuclear matter studies in [ 11 11 ] in its DDM3Y n and BDM3Y n density-dependent versions and found it to demonstrate excellent agreement with both the M3Y-Reid and M3Y-Paris interactions. The curves of the EOS of cold symmetric nuclear matter obtained with the DDB3Y1-, BDB3Y0-, BDB3Y1-, BDB3Y2- and BDB3Y3-Fetal interactions, found to have incompressibility values K 0 = 176, 196, 235, 351, and 467 MeV, respectively, were well illustrated in [ 11 11 ]. Each of these interactions has been found to reproduce the saturation of nuclear matter at density ρ = 0.17 fm -3 and binding energy per nucleon <math id="M6"><mrow><mfrac><mi>E</mi><mi>A</mi></mfrac></mrow></math> = 16 MeV correctly, demonstrating good agreement with the work of [ 12 12 ], which is a good basis for certainty and confidence that B3Y-Fetal will perform well in folding analyses. Importantly, DDB3Y1-, BDB3Y0-, BDB3Y1-, BDB3Y2- and BDB3Y3-Fetal interactions and their corresponding K-values are known to distinguish different nuclear EOSs in nuclear reactions and provide vital information on the nuclear matter EOS. In particular, the DDM3Y1 and BDM3Y1 versions of the B3Y-Fetal used in this work, the most realistic interactions, have given for SNM at equilibrium incompressibility K 0 <math id="M7"><mo>≃</mo></math> 176 – 235 MeV, and this compares favorably well with K 0 <math id="M8"><mo>≃</mo></math> 171 – 232 MeV and K 0 <math id="M9"><mo>≃</mo></math> 176 – 270 MeV predicted by the same versions of the M3Y-Reid and M3Y-Paris effective interactions, respectively, in [ 12 12 ]. When this is compared with the experimental estimate of K 0 = 200–350 MeV based on studies on giant monopole resonances [ 17 17 ] and the theoretical estimates of K 0 = 220±50 MeV and K 0 <math id="M10"><mo>≃</mo></math> 250–270 MeV [ 12 12 ], based on non-relativistic and relativistic mean-field models, respectively [ 18 18 ], for SNM, it is evident that the prediction of the B3Y-Fetal interaction is in good agreement. It is with this strong and positive expectation that we use it in the optical model analyses in this paper. For the folding calculation in this paper, the isoscalar component of the density- and energy-dependent B3Y-Fetal interaction is chiefly employed in all optical model analyses, based on the fact that the system considered involves spin-saturated nuclei. Specifically, the B3Y-Fetal is used in its DDM3Y1, BDM3Y1, BDM3Y2 and BDM3Y3 versions, along with the M3Y-Reid interaction, for folding analysis following the approach of Khoa et al. [ 7 7 , 8 8 ], with the intent of determining and comparing the performance of the first with the second interaction. 3 B3Y-Fetal in Folding Analysis 1 The double folding model is a well-established theory and computational procedure that provides an adequate description of the average properties of nucleus-nucleus scattering. In this model, the direct and exchange parts of the nucleus-nucleus potential are obtained by folding the densities of the projectile and target nuclei with the direct and exchange components of the effective NN interaction in Eq. (2). Thus, the resulting optical potential is <math display="block" id="M11"><mrow><msub><mi>V</mi><mtext>F</mtext></msub><mo stretchy="false">(</mo><mi>E</mi><mn>,</mn><mi>R</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>V</mi><mtext>D</mtext></msub><mo stretchy="false">(</mo><mi>E</mi><mn>,</mn><mi>R</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>V</mi><mrow><mtext>EX</mtext></mrow></msub><mo stretchy="false">(</mo><mi>E</mi><mn>,</mn><mi>R</mi><mo stretchy="false">)</mo></mrow></math> where V D ( E, R ) and V EX ( E, R ) are the direct and exchange components of the real folded potential, respectively, shown explicitly in [ 7 7 , 12 12 ]. The direct component of the folded potential is local in coordinate space, but the exchange is nonlocal because of the effect of anti-symmetrization occasioned by single-nucleon knock-on exchange, making its evaluation difficult; hence, a much simpler approximation, called the zero-range potential, has been used in many nuclear-nucleus calculations [ 19 19 , 20 20 ]. This approach has been used with success in the nucleus-nucleus calculations of heavy-ion optical potential at low energies, but has been found inadequate in refractive scattering, where the data are sensitive to the real optical potential over a wider radial domain. As a consequence, in this work, the study of the HI optical potential for the elastic scattering of the chosen nuclear system at various collision energies of interest is carried out using the finite-range evaluation of the exchange potential by means of the local momentum approximation proposed by Brieva and Rook [ 21 21 ]. This approach has been used with immense success by many researchers previously [ 7 7 , 9 9 ]. Therefore, to solve the non-locality problem, the exact, consistent microscopic approximation of the exchange potential developed by Khoa et al. [ 7 7 ] is employed in this work. Their approach produces an accurate local approximation, resulting from treating the relative motion locally as a plane wave [ 8 8 ]. The exchange potential is computed using an iterative method in this work, ensuring a careful and accurate evaluation of the knock-on exchange effects to obtain a realistic energy dependence of the real folded potential V ( E,R ), shown from previous calculations [ 1 1 , 7 7 ] to have a much stronger dependence on the exchange term than the intrinsic energy dependence g ( E ) in Eq. (2). Thus, the volume integral of the real folded potential per interacting nucleon pair, which gives a useful insight into the energy dependence of the real optical potential, is expressed as [ 7 7 ]: <math display="block" id="M12"><mrow><msub><mi>J</mi><mtext>R</mtext></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn>4</mn><mi>π</mi><msub><mi>N</mi><mtext>R</mtext></msub></mrow><mrow><msub><mi>A</mi><mn>1</mn></msub><msub><mi>A</mi><mn>2</mn></msub></mrow></mfrac><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mn>0</mn><mi>∞</mi></msubsup></mrow></mstyle><mo stretchy="false">[</mo><msub><mi>V</mi><mtext>D</mtext></msub><mo stretchy="false">(</mo><mi>E</mi><mn>,</mn><mi>r</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>V</mi><mrow><mtext>EX</mtext></mrow></msub><mo stretchy="false">(</mo><mi>E</mi><mn>,</mn><mi>r</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo><msup><mi>r</mi><mn>2</mn></msup><mtext>d</mtext><mi>r</mi></mrow></math> where A 1 and A 2 represent the two colliding nuclei. In our case, in which a conscious effort is made to assess the performance of the real folded potential derived from the B3Y-Fetal in relation to the M3Y-Reid-based folded potential, this is a very important distinguishing parameter. Now, considering the fact that the B3Y-Fetal and the M3Y-Reid effective interactions are real, the potential resulting from the double folding model is real and constitutes the real part of the optical potential in the OM analyses. This real potential, accounting for elastic scattering, is supplemented by an imaginary phenomenological potential that accounts for absorption into non-elastic channels. Therefore, as it is usually practiced in this hybrid approach, the imaginary phenomenological part has parameters that are adjusted to give the best fit to the elastic scattering data. Here, we chose the Woods-Saxon shape as the imaginary potential, so that the optical potential in the OM calculations is [ 7 7 ]: <math display="block" id="M13"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>U</mi><mo stretchy="false">(</mo><mi>E</mi><mn>,</mn><mi>R</mi><mo stretchy="false">)</mo><mo>=</mo></mrow></mtd><mtd columnalign="left"><mrow><msub><mi>N</mi><mtext>R</mtext></msub><mo stretchy="false">[</mo><msub><mi>V</mi><mtext>D</mtext></msub><mo stretchy="false">(</mo><mi>E</mi><mn>,</mn><mi>R</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>V</mi><mrow><mtext>EX</mtext></mrow></msub><mo stretchy="false">(</mo><mi>E</mi><mn>,</mn><mi>R</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>−</mo><mi>i</mi><msub><mi>W</mi><mtext>V</mtext></msub><msup><mrow><mo stretchy="false">[</mo><mn>1</mn><mo>+</mo><mi>exp</mi><mo stretchy="false">(</mo><mrow><mo>(</mo><mrow><mfrac><mrow><mi>R</mi><mo>−</mo><msub><mi>R</mi><mtext>V</mtext></msub></mrow><mrow><msub><mi>a</mi><mtext>V</mtext></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>+</mo><mn>4</mn><mi>i</mi><msub><mi>W</mi><mtext>D</mtext></msub><msub><mi>a</mi><mtext>D</mtext></msub><mfrac><mtext>d</mtext><mrow><mtext>d</mtext><mi>R</mi></mrow></mfrac><msup><mrow><mo stretchy="false">[</mo><mn>1</mn><mo>+</mo><mi>exp</mi><mo stretchy="false">(</mo><mrow><mo>(</mo><mrow><mfrac><mrow><mi>R</mi><mo>−</mo><msub><mi>R</mi><mtext>D</mtext></msub></mrow><mrow><msub><mi>a</mi><mtext>D</mtext></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo stretchy="false">)</mo><mo stretchy="false">]</mo></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mn>,</mn></mrow></mtd></mtr></mtable></mrow></math> where W V and W D , a V and a D , and R V and R D are the respective strengths, diffuseness, and radii of the volume (WSV) and surface (WSD) terms of the Woods-Saxon function, respectively [ 7 7 ]. These parameters, alongside the normalization factor N R , were adjusted to obtain the best fits to the elastic scattering angular distributions. The surface term of Eq. (7) was added for possible improvement of the agreement with the data [ 7 7 ] at certain incident angles. All nuclear densities in the folding model were taken as a two-parameter Fermi distribution with parameters that can reproduce the Shell-Model densities for the nuclei under consideration. All computational analyses of the optical model were carried out using the ECIS-97 code. These OM analyses are used to study the elastic scattering of 16 O+ 16 O at incident energies of 145, 250, 350, and 480 MeV. 4 Results and Discussion 1 Generally, the results obtained with the heavy-ion nucleus-nucleus optical potential used for the study of the elastic scattering of the nuclear system chosen in the present study demonstrate good agreement with results of previous works [ 7 7 , 8 8 ]. The results are presented in both tabular and graphical forms herein and discussed accordingly. For economy of space, we have combined some graphs that have similar features into one figure, and salient points relating to them are appropriately discussed. In discussing the results, attention is given primarily to the results obtained with the optical potential derived from the B3Y-Fetal effective interaction because it is the principal interaction used for folding analyses in this work. The performance of the optical potential based on the M3Y-Reid effective interaction is discussed only when it becomes necessary to compare it with that of the B3Y-Fetal effective interaction. The different versions of the M3Y-Reid and B3Y-Fetal interactions in Tables 2 Tables 2 and 3 3 are based on different versions of DDM3Yn and BDM3Yn density dependences and are grouped together in different groups for the purpose of comparison. Therefore, the DDM3Y1-Reid and DDB3Y1- Fetal interactions are derived from the DDM3Y1 density-dependent version; BDM3Y1-Reid and BDB3Y1-Fetal interactions are from the BDM3Y1 density dependence, BDM3Y2-Reid and BDB3Y2-Fetal are based on the BDM3Y2 density dependence, and BDM3Y3-Reid and BDB3Y3- Fetal interactions are based on the BDM3Y3 density dependence. These density dependencies form the basis for comparison of the B3Y-Fetal with the M3Y-Reid effective interaction in the nuclear reaction considered in this study. 1 1 Potential 1 1 N R 1 1 - J R 1 1 <math id="M14"><mrow><msup><mrow><mo>〈</mo><msubsup><mi>r</mi><mtext>R</mtext><mn>2</mn></msubsup><mo>〉</mo></mrow><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></mrow></msup></mrow></math> 1 1 W V 1 1 R V 1 1 a V 1 1 σ R 1 1 χ 2 1 1   1 1   1 1 (MeVfm 3 ) 1 1 (fm) 1 1 (MeV) 1 1 (fm) 1 1 (fm) 1 1 (mb) 1 1   2 1 [ 16 O+ 16 O, E Lab =145 MeV] 1 1 DDM3Y1-REID 1 1 0.8526 1 1 336.3 1 1 4.2130 1 1 17.533 1 1 6.0028 1 1 0.6905 1 1 1703 1 1 26.2 1 1   1 1 (0.9445) 1 1 (339.0) 1 1 (4.1574) 1 1 (17.133) 1 1 (6.0028) 1 1 (0.6905) 1 1 (1693) 1 1 (53.4) 1 1 DDB3Y1-FETAL 1 1 0.8548 1 1 332.1 1 1 4.1813 1 1 17.638 1 1 6.0028 1 1 0.6905 1 1 1700 1 1 18.3 1 1 BDM3Y1-REID 1 1 0.8782 1 1 340.5 1 1 4.2135 1 1 18.767 1 1 5.9726 1 1 0.7001 1 1 1725 1 1 20.4 1 1   1 1 (0.9653) 1 1 (340.3) 1 1 (4.1597) 1 1 (17.490) 1 1 (5.9726) 1 1 (0.7011) 1 1 (1702) 1 1 (55.4) 1 1 BDB3Y1-FETAL 1 1 0.8688 1 1 332.1 1 1 4.1812 1 1 17.929 1 1 5.9726 1 1 0.7001 1 1 1706 1 1 19.0 1 1 BDM3Y2-REID 1 1 0.9210 1 1 344.2 1 1 4.2205 1 1 19.829 1 1 5.9176 1 1 0.7128 1 1 1737 1 1 22.4 1 1   1 1 (1.0154) 1 1 (344.4) 1 1 (4.1706) 1 1 (18.531) 1 1 (5.9176) 1 1 (0.7128) 1 1 (1714) 1 1 (59.8) 1 1 BDB3Y2-FETAL 1 1 0.9108 1 1 336.4 1 1 4.1867 1 1 18.823 1 1 5.9176 1 1 0.7128 1 1 1716 1 1 20.5 1 1 BDM3Y3-REID 1 1 0.9345 1 1 352.1 1 1 4.2286 1 1 22.926 1 1 5.8798 1 1 0.6552 1 1 1679 1 1 27.4 1 1   1 1 (1.0588) 1 1 (341.2) 1 1 (4.2075) 1 1 (22.601) 1 1 (5.8798) 1 1 (0.6552) 1 1 (1673) 1 1 (70.9) 1 1 BDB3Y3-FETAL 1 1 0.9605 1 1 309.7 1 1 4.3863 1 1 22.335 1 1 5.8798 1 1 0.6552 1 1 1669 1 1 29.1 2 1 [ 16 O+ 16 O, E Lab =250 MeV] 1 1 Potential 1 1 N R 1 1 - J R 1 1 <math id="M15"><mrow><msup><mrow><mo>〈</mo><msubsup><mi>r</mi><mtext>R</mtext><mn>2</mn></msubsup><mo>〉</mo></mrow><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></mrow></msup></mrow></math> 1 1 W V 1 1 R V 1 1 a V 1 1 σ R 1 1 χ 2 1 1   1 1   1 1 (MeVfm 3 ) 1 1 (fm) 1 1 (MeV) 1 1 (fm) 1 1 (fm) 1 1 (mb) 1 1 DDM3Y1-REID 1 1 0.7637 1 1 283.1 1 1 4.2134 1 1 29.427 1 1 5.6184 1 1 0.7368 1 1 1772 1 1 14.6 1 1   1 1 (0.8905) 1 1 (303.7) 1 1 (4.1650) 1 1 (28.498) 1 1 (5.6184) 1 1 (0.7368) 1 1 (1762) 1 1 (46.0) 1 1 DDB3Y1-FETAL 1 1 0.8209 1 1 298.6 1 1 4.1817 1 1 29.281 1 1 5.6184 1 1 0.7368 1 1 1771 1 1 10.1 1 1 BDM3Y1-REID 1 1 0.7835 1 1 285.4 1 1 4.2142 1 1 29.188 1 1 5.6418 1 1 0.7349 1 1 1777 1 1 15.6 1 1   1 1 (0.9099) 1 1 (304.7) 1 1 (4.1675) 1 1 (28.429) 1 1 (5.6318) 1 1 (0.7349) 1 1 (1763) 1 1 (48.4) 1 1 BDB3Y1-FETAL 1 1 0.8366 1 1 299.4 1 1 4.1819 1 1 28.933 1 1 5.6418 1 1 0.7349 1 1 1774 1 1 10.4 1 1 BDM3Y2-REID 1 1 0.8135 1 1 285.7 1 1 4.2219 1 1 29.177 1 1 5.6600 1 1 0.7299 1 1 1783 1 1 17.5 1 1   1 1 (0.9510) 1 1 (306.3) 1 1 (4.1788) 1 1 (28.360) 1 1 (5.6600) 1 1 (0.7299) 1 1 (1767) 1 1 (57.0) 1 1 BDB3Y2-FETAL 1 1 0.8727 1 1 301.7 1 1 4.1879 1 1 28.745 1 1 5.6600 1 1 0.7339 1 1 1778 1 1 11.8 1 1 BDM3Y3-REID 1 1 0.8233 1 1 291.9 1 1 4.2293 1 1 29.238 1 1 5.7002 1 1 0.7131 1 1 1767 1 1 18.0 1 1   1 1 (0.9866) 1 1 (301.8) 1 1 (4.2160) 1 1 (28.270) 1 1 (5.7002) 1 1 (0.7002) 1 1 (1762) 1 1 (74.6) 1 1 BDB3Y3-FETAL 1 1 0.8270 1 1 247.8 1 1 4.3999 1 1 28.459 1 1 5.7002 1 1 0.7131 1 1 1755 1 1 16.6 2 1 [ 16 O+ 16 O, E Lab =350 MeV] 1 1 Potential 1 1 N R 1 1 - J R 1 1 <math id="M16"><mrow><msup><mrow><mo>〈</mo><msubsup><mi>r</mi><mtext>R</mtext><mn>2</mn></msubsup><mo>〉</mo></mrow><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></mrow></msup></mrow></math> 1 1 W V 1 1 R V 1 1 a V 1 1 σ R 1 1 χ 2 1 1   1 1   1 1 (MeVfm 3 ) 1 1 (fm) 1 1 (MeV) 1 1 (fm) 1 1 (fm) 1 1 (mb) 1 1 DDM3Y1-REID 1 1 0.8450 1 1 295.4 1 1 4.2149 1 1 33.274 1 1 5.6099 1 1 0.6534 1 1 1660 1 1 9.0 1 1   1 1 (0.9218) 1 1 (299.4) 1 1 (4.1728) 1 1 (31.628) 1 1 (5.6099) 1 1 (0.6534) 1 1 (1642) 1 1 (40.2) 1 1 DDB3Y1-FETAL 1 1 0.8566 1 1 292.7 1 1 4.1833 1 1 32.669 1 1 5.6099 1 1 0.6534 1 1 1652 1 1 7.9 1 1 BDM3Y1-REID 1 1 0.8724 1 1 299.6 1 1 4.2160 1 1 32.176 1 1 5.6766 1 1 0.6470 1 1 1669 1 1 7.9 1 1   1 1 (0.9435) 1 1 (300.9) 1 1 (4.1754) 1 1 (31.430) 1 1 (5.6366) 1 1 (0.6470) 1 1 (1642) 1 1 (41.2) 1 1 BDB3Y1-FETAL 1 1 0.8735 1 1 293.6 1 1 4.1837 1 1 31.516 1 1 5.6766 1 1 0.6470 1 1 1659 1 1 7.9 1 1 BDM3Y2-REID 1 1 0.9150 1 1 302.6 1 1 4.2242 1 1 32.5583 1 1 5.7042 1 1 0.6330 1 1 1663 1 1 8.8 1 1   1 1 (0.9942) 1 1 (304.8) 1 1 (4.1870) 1 1 (30.980) 1 1 (5.6999) 1 1 (0.6330) 1 1 (1645) 1 1 (43.4) 1 1 BDB3Y2-FETAL 1 1 0.9164 1 1 297.4 1 1 4.1903 1 1 31.585 1 1 5.7042 1 1 0.6330 1 1 1651 1 1 8.0 1 1 BDM3Y3-REID 1 1 0.9141 1 1 306.1 1 1 4.2309 1 1 33.458 1 1 5.7148 1 1 0.6301 1 1 1671 1 1 9.4 1 1   1 1 (1.0369) 1 1 (301.7) 1 1 (4.2248) 1 1 (31.914) 1 1 (5.7048) 1 1 (0.6301) 1 1 (1652) 1 1 (48.6) 1 1 BDB3Y3-FETAL 1 1 1.0556 1 1 294.8 1 1 4.4149 1 1 36.870 1 1 5.7148 1 1 0.6301 1 1 1701 1 1 14.6 2 1 [ 16 O+ 16 O, E Lab =480 MeV] 1 1 DDM3Y1-REID 1 1 0.8018 1 1 259.9 1 1 4.2184 1 1 33.505 1 1 5.3784 1 1 0.7434 1 1 1662 1 1 8.7 1 1   1 1 (0.8520) 1 1 (259.8) 1 1 (4.1838) 1 1 (32.675) 1 1 (5.3984) 1 1 (0.7434) 1 1 (1657) 1 1 (6.7) 1 1 DDB3Y1-FETAL 1 1 0.7987 1 1 251.6 1 1 4.1871 1 1 32.790 1 1 5.3784 1 1 0.7434 1 1 1652 1 1 8.5 1 1 BDM3Y1-REID 1 1 0.8119 1 1 258.5 1 1 4.2198 1 1 32.531 1 1 5.42774 1 1 0.7377 1 1 1664 1 1 8.9 1 1   1 1 (0.8741) 1 1 (261.6) 1 1 (4.1866) 1 1 (32.517) 1 1 (5.4274) 1 1 (0.7377) 1 1 (1660) 1 1 (6.9) 1 1 BDB3Y1-FETAL 1 1 0.8044 1 1 249.2 1 1 4.1878 1 1 31.715 1 1 5.4274 1 1 0.7377 1 1 1653 1 1 9.2 1 1 BDM3Y2-REID 1 1 0.9765 1 1 299.3 1 1 4.2288 1 1 33.183 1 1 5.7948 1 1 0.6030 1 1 1638 1 1 19.8 1 1   1 1 (0.9116) 1 1 (262.2) 1 1 (4.1986) 1 1 (32.988) 1 1 (5.6948) 1 1 (0.6030) 1 1 (1588) 1 1 (10.5) 1 1 BDB3Y2-FETAL 1 1 0.9783 1 1 292.5 1 1 4.1950 1 1 33.999 1 1 5.7008 1 1 0.6030 1 1 1602 1 1 13.0 1 1 BDM3Y3-REID 1 1 0.8272 1 1 257.4 1 1 4.2342 1 1 33.700 1 1 5.4213 1 1 0.7503 1 1 1692 1 1 11.0 1 1   1 1 (1.0224) 1 1 (278.9) 1 1 (4.2370) 1 1 (33.521) 1 1 (5.5213) 1 1 (0.7503) 1 1 (1735) 1 1 (10.6) 1 1 BDB3Y3-FETAL 1 1 0.9445 1 1 240.6 1 1 4.4370 1 1 35.121 1 1 5.4213 1 1 0.7503 1 1 1709 1 1 11.0 The results obtained in Ref. [ 7 7 ] using the M3Y-Reid interaction are in brackets. 1 1 Potential 1 1 N R 1 1 W V 1 1 R V 1 1 a V 1 1 W D 1 1 R D 1 1 a D 1 1 σ R 1 1 χ 2 1 1   1 1   1 1 (MeV) 1 1 (fm) 1 1 (fm) 1 1 (MeV) 1 1 (fm) 1 1 (fm) 1 1 (mb) 2 1 [ 16 O+ 16 O, E Lab =250 MeV (WSV+WSD)] 1 1 DDM3Y1-REID 1 1 0.7764 1 1 37.433 1 1 4.7915 1 1 0.9062 1 1 3.5230 1 1 5.8309 1 1 0.4268 1 1 1830 1 1 11.2 1 1   1 1 (0.8905) 1 1 (36.641) 1 1 (4.7915) 1 1 (0.9062) 1 1 (3.5230) 1 1 (5.8309) 1 1 (0.4268) 1 1 (1822) 1 1 (37.4) 1 1 DDB3Y1-FETAL 1 1 0.8336 1 1 38.086 1 1 4.7915 1 1 0.9062 1 1 3.5230 1 1 5.8309 1 1 0.4268 1 1 1836 1 1 9.8 1 1 BDM3Y1-REID 1 1 0.7945 1 1 38.620 1 1 4.7039 1 1 0.9212 1 1 3.9878 1 1 5.8064 1 1 0.4267 1 1 1835 1 1 10.7 1 1   1 1 (0.9099) 1 1 (37.876) 1 1 (4.7083) 1 1 (0.9219) 1 1 (3.9878) 1 1 (5.8064) 1 1 (0.4267) 1 1 (1829) 1 1 (37.2) 1 1 BDB3Y1-FETAL 1 1 0.8499 1 1 39.2118 1 1 4.7038 1 1 0.9212 1 1 3.9878 1 1 5.8064 1 1 0.4267 1 1 1841 1 1 9.4 1 1 BDM3Y2-REID 1 1 0.8259 1 1 38.056 1 1 4.6464 1 1 0.9313 1 1 5.1053 1 1 5.7267 1 1 0.4231 1 1 1827 1 1 10.0 1 1   1 1 (0.9510) 1 1 (38.528) 1 1 (4.6464) 1 1 (0.9313) 1 1 (5.1053) 1 1 (5.7267) 1 1 0.4231 1 1 (1831) 1 1 (37.6) 1 1 BDB3Y2-FETAL 1 1 0.8830 1 1 38.804 1 1 4.6464 1 1 0.9313 1 1 5.1013 1 1 5.7267 1 1 0.4231 1 1 1833 1 1 9.0 1 1 BDM3Y3-REID 1 1 0.8269 1 1 39.258 1 1 4.6109 1 1 0.9355 1 1 5.6530 1 1 5.7518 1 1 0.4012 1 1 1830 1 1 9.5 1 1   1 1 (0.9866) 1 1 (39.836) 1 1 (4.6109) 1 1 (0.9355) 1 1 (5.6530) 1 1 (5.7518) 1 1 0.4012 1 1 (1835) 1 1 (37.6) 1 1 BDB3Y3-FETAL 1 1 0.8836 1 1 39.745 1 1 4.6109 1 1 0.9355 1 1 5.6530 1 1 5.7518 1 1 0.4012 1 1 1834 1 1 10.8 2 1 [ 16 O+ 16 O, E Lab =350 MeV (WSV+WSD)] 1 1 DDM3Y1-REID 1 1 0.8416 1 1 29.9001 1 1 5.7554 1 1 0.6488 1 1 7.9655 1 1 4.7371 1 1 0.3637 1 1 1683 1 1 8.1 1 1   1 1 (0.9218) 1 1 (25.265) 1 1 (5.7554) 1 1 (0.6488) 1 1 (7.9655) 1 1 (4.7371) 1 1 (0.3637) 1 1 (1634) 1 1 (39.8) 1 1 DDB3Y1-FETAL 1 1 0.8452 1 1 29.676 1 1 5.7554 1 1 0.6488 1 1 7.9655 1 1 4.7371 1 1 0.3637 1 1 1678 1 1 10.9 1 1 BDM3Y1-REID 1 1 0.8560 1 1 30.079 1 1 5.7578 1 1 0.6481 1 1 7.4643 1 1 4.7473 1 1 0.3823 1 1 1685 1 1 8.2 1 1   1 1 (0.9435) 1 1 (25.429) 1 1 (5.7578) 1 1 (0.6481) 1 1 (7.4643) 1 1 (4.7473) 1 1 (0.3823) 1 1 (1637) 1 1 (40.6) 1 1 BDB3Y1-FETAL 1 1 0.8624 1 1 29.814 1 1 5.7578 1 1 0.6481 1 1 7.4643 1 1 4.7473 1 1 0.3823 1 1 1679 1 1 10.8 1 1 BDM3Y2-REID 1 1 0.9057 1 1 30.704 1 1 5.7612 1 1 0.6458 1 1 6.2634 1 1 4.8143 1 1 0.4240 1 1 1690 1 1 8.4 1 1   1 1 (0.9942) 1 1 (25.929) 1 1 (5.7612) 1 1 (0.6458) 1 1 (6.2634) 1 1 (4.8143) 1 1 (0.4240) 1 1 (1645) 1 1 (42.3) 1 1 BDB3Y2-FETAL 1 1 0.9037 1 1 30.226 1 1 5.7612 1 1 0.6458 1 1 6.2634 1 1 4.8143 1 1 0.4240 1 1 1682 1 1 10.3 1 1 BDM3Y3-REID 1 1 0.9050 1 1 31.690 1 1 5.7692 1 1 0.6309 1 1 4.5614 1 1 4.9566 1 1 0.5031 1 1 1680 1 1 8.8 1 1   1 1 (1.0369) 1 1 (26.304) 1 1 (5.7692) 1 1 (0.6309) 1 1 (4.5614) 1 1 (4.9566) 1 1 (0.5031) 1 1 (1644) 1 1 (51.5) 1 1 BD3Y3-FETAL 1 1 0.8363 1 1 28.753 1 1 5.6792 1 1 0.6309 1 1 4.5614 1 1 4.9566 1 1 0.5031 1 1 1646 1 1 18.4 4.1 Total Real Folded Potential 2 The optical model analyses performed in this work have proven beyond reasonable doubt that nuclear matter calculations [ 11 11 ] alone do not suffice for an explicit description of the form and character of the B3Y-Fetal effective interaction. Therefore, the optical potential results shown in Tables 2 Tables 2 and 3 3 , as well as Figs. 1 Figs. 1 , 2 2 , 3 3 , 4 4 , 5 5 , 6 6 and 7 7 , are more revealing. The foremost revelation emanating from the folding analyses is that of form. Figures 1 Figures 1 and 2 2 unambiguously show that one of the underlying reasons for the performance difference is closely associated with the fact that different effective interactions have different forms. Figure 1 Figure 1 reveals that the optical potential based on the B3Y-Fetal interaction has a repulsive direct term (Eq.(5)), which is combined with an attractive exchange term (Eq. (5)) of large magnitude to produce a total real folded potential that is attractive. In contrast, Figure 2 Figure 2 shows that the optical potential involving the M3Y-Reid interaction has an attractive direct term and an attractive exchange term of approximately the same magnitude that combine to form an attractive total real folded potential. This basic difference between the two interactions was observed during folding analyses to have been translated into the response of each effective interaction to the elastic scattering data in various ways. Figure 1 Figure 1 shows the repulsive direct V D (E,R) and attractive exchange V EX (E,R) contributions to the total real folded potential for the 16 O+ 16 O system, computed with the DDB3Y1-FETAL interaction at incident energies from 9 to 30 MeV/nucleon. The plots clearly show that most of the energy dependence of the HI optical potential comes from the exchange term, which dominates at smaller inter-nuclear distances, especially at low energies; this suggests that the density-dependent contribution from V EX is also much stronger than that from V D . However, the contribution of V EX becomes comparable to and even less than that of V D as the overlap density decreases exponentially in the surface region. Thus, the total HI potential is dominated by the direct component at distances R ≥ 7 fm, and this explains why the simple folding model [ 1 1 ] has been very successful in cases where the scattering data is only sensitive to the tail of the optical potential. It is also evident from Fig. 1 Fig. 1 that the exchange potential becomes less dominant with increasing incident energy, so that both the exchange and direct potentials have approximately the same strength even at small distances. The result is that the total real HI potential becomes increasingly repulsive with incident energy; this is in good agreement with the work of [ 7 7 ]. θ r, the fact that the absolute strength of the total real folded potential is governed to some extent by the shape of the chosen effective NN forces makes it quite clear that the real folded potentials based on the B3Y-Fetal become more repulsive at higher incident energies than their counterparts based on M3Y-Reid, whose direct and exchange components are attractive, as shown in Fig. 2 Fig. 2 . Therefore, when compared with the total HI potential of DDM3Y1-Reid, the total HI potential of DDB3Y1-Fetal becomes more repulsive with increasing bombarding energy because of its peculiar nature. 4.2 Fits to 16 O+ 16 O Elastic Data 2 The elastic scattering of this system was studied at incident energies of 145, 250, 350, and 480 MeV. All experimental data were obtained from Khoa et al. [ 7 7 , 8 8 ] through helpful communication. The data at E Lab = 250, 350, and 480 MeV were originally measured using the Q3D magnetic spectrometer at the cyclotron of the Hahn-Meitner Institute (HMI) in Berlin, Germany [ 22 22 - 24 24 22-24 ], while the data at E Lab = 145 MeV were measured by Sugiyama et al. at the Japan Atomic Energy Research Institute (JAERI) in Tokai [ 25 25 ]; all these data were obtained by Khoa et al. [ 7 7 , 8 8 ] through collaboration. The results produced by various real folded potentials with the volume WS imaginary potential are presented in Table 2 Table 2 and plotted in Fig. 3 Fig. 3 . In each case, a reasonable agreement with the experimental data throughout the entire angular range can be observed. The deep DDB3Y1-Fetal and BDB3Y1-Fetal folded potentials proved to be the best, reproducing the data over the entire angular range. The BDB3Y3-Fetal potential, followed by the BDB3Y2-Fetal potential, had the poorest performance, providing an incorrect description of the data at large angles. On average, the renormalization of the real folded potentials in these cases is approximately 0.85, giving J R values of the real optical potential that are very close to those obtained by Khoa and co-researchers [ 7 7 ]. To improve the agreement of the calculated cross-sections with the data in the large-angle (refractive) region at 250 and 350 MeV, a surface WS term, WSD, was added to the imaginary potential, with the real folded potentials kept fixed to determine the effects of the imaginary potential. Consequently, the surface WSD potential did improve fits to the data at 250 and 350 MeV, respectively. The improved fits are respectively recorded in Table 3 Table 3 and plotted in Figs. 4 Figs. 4 and 5 5 . As recorded in Tables 3a Tables 3a and 3b 3b , the improved fits are more pronounced at 250 MeV than at 350 MeV. Clearly, a comparison of Tables 3a Tables 3a and 2b 2b , representing E Lab = 250 MeV, shows that the values of χ 2 in Table 3a Table 3a are generally lower; this indicates improved fits. Quantitatively, the values of χ 2 in Table 3a Table 3a obtained in this work with the DDM3Y1-REID, BDM3Y1-REID, BDM3Y2-REID, and BDM3Y3-REID are 9.0, 7.9, 8.8 and 9.4 (with WSD) versus 14.6, 15.6, 17.5 and 18.0, respectively, in Table 1b Table 1b (without WSD). Comparing these two sets of values of χ 2 , the differences (reductions) in the values of χ 2 obtained with the DDM3Y1-REID, BDM3Y1-REID, BDM3Y2-REID, and BDM3Y3-REID potentials with the inclusion of WSD are 5.6, 7.7, 8.7, and 8.6, respectively; these represent improved fits. Similarly, the B3Y-Fetal-based potentials DDB3Y1-FETAL, BDB3Y1-FETAL, BDB3Y2-FETAL, and BDB3Y3-FETAL are shown in Table 3a Table 3a to produce 9.8, 9.4, 9.0, and 10.8, respectively, as values of χ 2 at 250 MeV (with WSD). When these values are compared with 10.1, 10.4, 11.8, and 16.6, obtained with the same set of folded potentials (arranged in the same order) as the values of χ 2 in Table 2b Table 2b (without WSD), the differences (reductions) in the values of χ 2 with the inclusion of WSD are 0.3, 1.0, 2.8, and 5.8, respectively, indicating some improvement of the fits at 250 MeV. Here, without WSD, the fits are worse for the M3Y-REID than for the B3Y-FETAL-based potentials in Table 2b Table 2b ; therefore, with the inclusion of WSD in Table 3a Table 3a , the improvement of the fits is more pronounced for the former than for the latter. The fact is that the inclusion of WSD in our optical model analysis has actually improved the fits, but the improvement cannot be clearly seen with the eye in Figs. 4 Figs. 4 and 5 5 . This is in agreement with the results of Ref. [ 7 7 ], which suggest that a combination of the imaginary potentials makes the absorption at the surface strong enough to reproduce the diffractive part of the scattering, while the absorption at small distances becomes weaker (upper part of Figure 5 Figure 5 ), giving rise to a broad maximum in the elastic cross-section. (the lower part of Figure 5 Figure 5 ). With the M3Y-Reid-based potentials, the improvement of fits at 350 MeV is marginal. In Table 3b Table 3b , the DDM3Y1-REID, BDM3Y1-REID, BDM3Y2-REID, and BDM3Y3-REID produced 8.1, 8.2, 8.4, and 8.8 respectively as values of χ 2 at 350 MeV (with WSD added), compared with 9.0, 7.9, 8.8, and 9.4, respectively, in Table 2c Table 2c (without WSD), and the corresponding differences (reductions) are 0.8, 0.3, 0.4, and 0.6, respectively. In this case, a difference of only 0.3 does not indicate an improved fit, whereas the others are marginal differences that represent improved fits. When the differences obtained with these potentials at 250 MeV are compared with the differences observed at 350 MeV, it is evident that the improved fits were more pronounced at 250 MeV than at 350 MeV. Figure 6 Figure 6 shows the improved fit obtained with the M3Y-Reid at 350 MeV to show that the folded potentials derived from the B3Y-Fetal are in good agreement by comparison. Even though the fits have been improved at these energies, the relative differences in the shapes of the different folded potentials remain practically the same, but it is noteworthy that the best fit to the large-angle data is given by the DDM3Y1 and BDM3Y1 versions of the folded potentials, which are the most realistic folded potentials. From past work [ 5 5 , 7 7 ], several OM analyses of the data at 350 MeV have revealed that the broad maximum near 50, identified as a remnant of the primary rainbow, (with the first Airy minimum around 44), can only be reproduced by real optical potentials that are deep enough to belong to the group of refractive potentials. Accordingly, the results of the present work have shown the best-fit folded potentials that consistently produce the evolution of Airy oscillation with energy to be deeply refractive. To show and confirm this, we have followed the procedure of Khoa and collaborators [ 7 7 ] to plot the elastic scattering cross-sections arising from the unsymmetrized calculations of the 16 O+ 16 O system at 145, 250, and 350 MeV obtained with the best-fit folded potential, DDB3Y1-Fetal, in Figure 7 Figure 7 for the sole purpose of studying the evolution of the rainbow pattern, in which Airy minima were observed to shift smoothly toward small scattering angles with increasing incident energy. The plots show the existence of Airy minima at the aforementioned energies. Thus, at an energy of 145 MeV, the first Airy minimum occurs at 100, followed by the second at 75 and the third at 48. What is interesting here is the continuous shifting of the Airy minima toward the diffractive region as the incident energy increases. Accordingly, it can be seen from Figure 10 that the first Airy minimum shifts from approximately 100 at the energy of 145 MeV to approximately 66 at 250 MeV (solid curve), and at 350 MeV, it is seen to have shifted to approximately 44. This shows that our version of the folding model, derived from the B3Y-Fetal, has proven to be a reliable HI optical potential at different energies with the DDB3Y1- or BDB3Y1-Fetal potential as a trustworthy tool to predict scattering cross-sections at large angles. 4.3 Energy Dependence of the Real HI Potential 2 Many folding analysis results have shown that the renormalization factor N R of the real folded potential should be very close to unity for the model used for the nucleus-nucleus potential to be adjudged satisfactorily realistic [ 7 7 , 8 8 ]. However, higher-order terms such as dynamic polarization correction, which is mostly repulsive and surface-peaked, contribute to both the real and imaginary parts of the HI potential, where the folded potential is a first-order term. N R smaller than one might be considered reasonable [ 7 7 ]. Consequently, the N R values ranging from 0.8 to 0.9 obtained in the present folding analyses can be considered reasonable. Furthermore, all N R values obtained using the best-fit real folded potential (DDB3Y1-Fetal and DDM3Y1-Reid) in the analyses of 16 O and 12 C scattering data are plotted against the incident energy in the upper regions of Figs. 8 Figs. 8 and 9 9 . The plots show that N R depends weakly on energy and varies around a value of 0.9. The BDM3Y1 exhibits the same behavior. Therefore, it is evident that the intrinsic energy dependence arising from the exchange effects and the weak explicit energy dependence g ( E ) factor introduced into the effective NN interaction in Eq. (6) are actually the main sources of the total energy dependence of the real HI optical potential. In the lower regions of Figs. 8 Figs. 8 and 9 9 , the volume integral of the real folded potential per interacting nucleon pair J R is seen to decrease smoothly from approximately 340 MeV to 150 MeV fm 3 as the incident energy increases from approximately 10 to 120 MeV/nucleon. This trend agrees well with previous studies. [ 7 7 ], which shows that prediction of the energy dependence of the real HI potential using the folding model in this work is reliably satisfactory. Regarding the overall performance of the two effective interactions in the OM analyses of elastic scattering data, it is worth noting that the B3Y-Fetal has so far demonstrated acceptable agreement with the M3Y-Reid. Tables 2 Tables 2 and 3 3 show that the values of the renormalization factor N R obtained with the B3Y-Fetal are slightly higher than those obtained with the M3Y-Reid, so that the corresponding reaction cross-sections σ R for the former are lower than those for the latter in almost all cases. This level of performance by B3Y-Fetal is quite impressive. A well-established trend that both interactions portray is one in which the optical potential becomes progressively repulsive as the incident energy increases, which is in agreement with the work of Nuclear Physics scholars [ 1 1 , 7 7 ]. 4.4 EOS of Normal Nuclear Matter 2 The folding analyses in this study have shown that the DDM3Y1-Fetal ( K = 176 MeV) and BDM3Y1-Fetal ( K = 235 MeV) are the best-fit and most realistic optical potentials. This is also true for the folded potentials derived from M3Y-Reid (DDM3Y1- and BDM3Y1-Reid) effective interactions. When these results bear on the results of nuclear matter calculations reproduced and discussed in Section 2 (page 6), a strong factual link for establishing the nature of the EOS of symmetric nuclear matter (SNM) using incompressibility values clearly evolves. It is known that the incompressibilities of 351 MeV and 467 MeV produced by the BDB3Y2-Fetal and BDB3Y3-Fetal represent stiff EOS of SNMs [ 3 3 , 12 12 ]. These interactions are shown in Figs. 1 Figs. 1 and 2 2 to be more repulsive and in Figs. 3 Figs. 3 , 4 4 , 5 5 to provide a very bad description of the elastic data of the 16 O+ 16 O system, whereas the DDB3Y1- and BDB3Y1-Fetal, with incompresibilities of 176 MeV and 235 MeV, respectively, are shown to be better and in good agreement. Because the well-known standard is that the theoretical EOSs that predict higher K 0 values of approximately 300 MeV are often called “stiff”, whereas those that predict smaller K 0 values of approximately 200 MeV are said to be “soft” [ 12 12 ], the agreement demonstrated by the DDB3Y1- and BDB3Y1-Fetal herein shows that cold nuclear matter possibly has an underlying soft EOS. An additional important conclusion may be drawn from the present study of refractive HI scattering: the corresponding nuclear incompressibility K = 170 - 240 MeV [ 7 7 ] is the most realistic EOS for cold nuclear matter. It is hoped that this conclusion will be of great interest for some studies [ 26 26 , 27 27 ] in Astrophysics where the nuclear EOS [ 28 28 ] is a key component of the calculation. 5 Conclusion 1 The focus of this paper is the study of the elastic scattering of 16 O+ 16 O nuclei using a heavy-ion nucleus-nucleus optical potential based on four types of real folded potential derived from a new M3Y-type effective interaction called the B3Y-Fetal interaction within the framework of an optical model. The major findings of this study are summarized in the following concluding statements. • Folding analyses are observed in this work to be a computational necessity for a complete description of the form and character of the B3Y-Fetal interaction. Figure 1 Figure 1 provides convincing information concerning the form and character of this new effective interaction, which was not revealed in nuclear matter calculations in our earlier paper [ 11 11 ]. • The basic difference between the optical potentials derived from the M3Y-Reid and B3Y-Fetal interactions, in terms of form and character, is closely tied to their direct components. The optical potential based on the B3Y-Fetal has a repusive direct component, whereas the one based on M3Y-Reid has an attractive direct component. • The heavy-ion optical potential derived from the B3Y-Fetal has been found to be deep and attractive at small inter-nuclear distances, especially at low energies, as shown in Fig. 1 Fig. 1 . Furthermore, it has demonstrated weak absorption in the chosen nuclear system. • Results have shown the DDB3Y1-Fetal and BDB3Y1-Fetal to be the best-fit folded potentials, followed by the BDB3Y2-Fetal, with the BDB3Y3-Fetal as the optical potential with the poorest fit to the elastic data of the nuclear system at all energies considered in this work. This trend of performance of the optical potentials, which is also true for the optical potentials based on the M3Y-Reid interaction, is again in excellent agreement with the findings of Khoa et al. [ 7 7 ]. • The best-fit folded potentials, DDB3Y1- and BDB3Y1-Fetal, were found to reproduce reasonably well the energy dependence of the real optical potential for 16 O scattering found in previous OM analyses. The results obtained herein also show them to be of the family of deep refractive potentials, as they consistently describe Airy-like structures at large scattering angles observed in the 16 O scattering data at different energies. In addition, the folding analyses in this study of this family of folded potentials, known to have microscopically generated a nuclear matter equation of state with incompressibility K = 176–235 MeV [ 11 11 ], suggest strongly that cold nuclear matter is governed by a soft EOS. On a final note, the performance of the B3Y-Fetal in optical model analyses has been compared with that of the M3Y-Reid interaction, and the agreement between them has been found to be quite impressive, as they were seen to overlap and overtake each other with a marginal performance gap at various points in our OM analyses. References G.R. Satchler and W. G. Love . Folding Model Potentials from Realistic Interactions for Heavy-ion Scattering . Physics Reports (Review Section of Physics Letters). 55 , 90081-4 ( 1979 ). doi: 10.1016/0370-1573(79)90081-4 http://doi.org/10.1016/0370-1573(79)90081-4 . Folding Model Potentials from Realistic Interactions for Heavy-ion Scattering D. T. Khoa , V. W. Oertzen , H.G. Bohlen et al ., Study of Diffractive and Refractive Structure in the Elastic 16O+16O Scattering at Incident Energies Ranging from 124 to 1120 MeV . Nuclear Physics A. 672 , 00856-8 ( 2000 ). doi: 10.1016/S0375-9474(99)00856-8 http://doi.org/10.1016/S0375-9474(99)00856-8 . Study of Diffractive and Refractive Structure in the Elastic 16O+16O Scattering at Incident Energies Ranging from 124 to 1120 MeV M.E. Brandan and G.R. Satchler . The Interaction between Light Heavy Ions and What It Tells Us . Physics Reports. 285 , 00048-8 ( 1997 ). doi: 10.1016/S0370-1573(96)00048-8 http://doi.org/10.1016/S0370-1573(96)00048-8 . The Interaction between Light Heavy Ions and What It Tells Us D.T. Khoa , E. Khan , G. Colo et al ., Folding Model Analysis of Elastic and Inelastic Proton Scattering on Sulphur Isotopes . Nuclear Physics A. 706 , 00866-7 ( 2002 ). doi: 10.1016/S0375-9474(02)00866-7 http://doi.org/10.1016/S0375-9474(02)00866-7 . Folding Model Analysis of Elastic and Inelastic Proton Scattering on Sulphur Isotopes M. Katsuma , Y. Sakuragi , S. Okabe et al ., Coupled-Channels Study of the Nuclear Rainbow Phenomenon for the 12C+ 16O System . Progress of Theoretical Physics. 107 , 377 ( 2002 ). doi: 10.1143/PTP.107.377 http://doi.org/10.1143/PTP.107.377 . Coupled-Channels Study of the Nuclear Rainbow Phenomenon for the 12C+ 16O System S. M. Wong . Introductory Nuclear Physics . 2 ed. Toronto: Prentice-Hall International, Inc . ( 2004 ). doi: 10.1002/9783527617906 http://doi.org/10.1002/9783527617906 . Introductory Nuclear Physics D. T. Khoa , V. W. Oertzen and H. G. Bohlen . Double-Folding Model for Heavy-ion Optical Potential: Revised and Applied to Study 12C and 16O Elastic Scattering . Physical Review C. 49 , 3 ( 1994 ). doi: 1103/PhysRevC.49 http://doi.org/1103/PhysRevC.49 . Double-Folding Model for Heavy-ion Optical Potential: Revised and Applied to Study 12C and 16O Elastic Scattering D.T. Khoa , G.R. Satchler and W. V. Oertzen . Nuclear Incompressibility and Density Dependent NN Interactions in the Folding Model for Nucleus Potentials . Physical Review C. 56 , ( 1997 ). doi: 10.1103/PhysRevC.56 http://doi.org/10.1103/PhysRevC.56 . Nuclear Incompressibility and Density Dependent NN Interactions in the Folding Model for Nucleus Potentials G. R. Satchler . Introduction to Nuclear Reactions 2 ed. Macmillan Education Ltd . London . 20532-8 ( 1990 ). doi: 10.1007/978-1-349-20532-8 http://doi.org/10.1007/978-1-349-20532-8 . Introduction to Nuclear Reactions J. O. Fiase , K.R.S. Devan and A. Hosaka . Mass Dependence of M3Y-Type Interactions and the Effects of Tensor Correlations . Physical Review C. 66 , 014004 ( 2002 ). doi: 10.1103/PhysRevC.66.014004 http://doi.org/10.1103/PhysRevC.66.014004 . Mass Dependence of M3Y-Type Interactions and the Effects of Tensor Correlations I. Ochala and J. O. Fiase . Symmetric Nuclear Matter Calculations - A Variational Approach . Physical Review C. 98 , 064001 ( 2018 ). doi: 10.1103/PhysRevC.98.064001 http://doi.org/10.1103/PhysRevC.98.064001 . Symmetric Nuclear Matter Calculations - A Variational Approach D. T. Khoa and V. W. Oertzen . Refractive Alpha-Nucleus Scattering; A Probe for the Incompressibility of Cold Nuclear Matter . Physics Letters B. 342 , 01393-Q ( 1995 ). doi: 10.1016/0370-2693(94)01393-Q http://doi.org/10.1016/0370-2693(94)01393-Q . Refractive Alpha-Nucleus Scattering; A Probe for the Incompressibility of Cold Nuclear Matter P.E. Hodgson . The Nucleon Optical Model . World Scientific Publishing Co. Ltd. London . 2317 ( 1994 ). doi: 10.1142/2317 http://doi.org/10.1142/2317 . The Nucleon Optical Model S. Hamada , N. Burtebayev , N. Amangeldi et al ., Phenomenological and Semi-microscopic Analysis for 16O and 12C Elastically Scattering on the Nucleus of 16O and 12C at Energies Near Coulomb Barrier . Journal of Physics: Conference Series. 381 , 012130 ( 2012 ). doi: 10.1088/1742-6596/381/1/012130 http://doi.org/10.1088/1742-6596/381/1/012130 . Phenomenological and Semi-microscopic Analysis for 16O and 12C Elastically Scattering on the Nucleus of 16O and 12C at Energies Near Coulomb Barrier G. Bertsch , J. Borsowicz , H. McManus et al ., Interactions for Inelastic Scattering Derived from Realistic Potentials . Nuclear Physics A. 284 , 90392-X ( 1977 ). doi: 10.1016/0375-9474(77)90392-X http://doi.org/10.1016/0375-9474(77)90392-X . Interactions for Inelastic Scattering Derived from Realistic Potentials I. Gontcharn , D. J. Hinde , M. Dasgupta et al ., Double Folding Nucleus-Nucleus Potential Applied to Heavy-Ion Fusion Reactions . Physical Review C. 69 , 02461 ( 2004 ). doi: 10.1103/PhysRevC.69.02461 http://doi.org/10.1103/PhysRevC.69.02461 . Double Folding Nucleus-Nucleus Potential Applied to Heavy-Ion Fusion Reactions S. Shlomo and D. H. Youngblood . Nuclear Matter Incompressibility and Giant Monopole Resonance . Nuclear Physics A. 569 , 90121-X ( 1994 ). doi: 10.1016/0375-9474(94)90121-X http://doi.org/10.1016/0375-9474(94)90121-X . Nuclear Matter Incompressibility and Giant Monopole Resonance G. Colo and N. V. Giai . Theoretical Understanding of the Nuclear Matter Incompressibility: Where Do We Stand? Nuclear Physics A. 731 , 014 ( 2004 ). doi: 10.1016/j.nuclphysa.2003.11.014 http://doi.org/10.1016/j.nuclphysa.2003.11.014 . Theoretical Understanding of the Nuclear Matter Incompressibility: Where Do We Stand? K.C. Panda , B.C. Sahu and J. Bhoi . Accuracy of Simple Folding Model in the Calculation of the Direct Part of Real α-α Interaction Potential . PRAMANA Journal of Physics. 82 , 0737-2 ( 2014 ). doi: 10.1007/s12043-014-0737-2 http://doi.org/10.1007/s12043-014-0737-2 . Accuracy of Simple Folding Model in the Calculation of the Direct Part of Real α-α Interaction Potential L. Trache , A. Azhari , H. L. Clark et al ., Optical Model Potential Involving Loosely Bound p-Shell Nuclei Around 10 MeV/ nucleon . Physical Review C. 61 , 024612 ( 2000 ). doi: 10.1103/PhysRevC.61.024612 http://doi.org/10.1103/PhysRevC.61.024612 . Optical Model Potential Involving Loosely Bound p-Shell Nuclei Around 10 MeV/ nucleon F. A. Brieva and J. R. Rook . Nucleon-Nucleus Optical Model (1). Nuclear Matter Approach . Nuclear Physics A. 291 , 90322-0 ( 1977 ). doi: 10.1016/0375-9474(77)90322-0 http://doi.org/10.1016/0375-9474(77)90322-0 . Nucleon-Nucleus Optical Model (1). Nuclear Matter Approach E. Stiliaris , H. G. Bohlen , P. Frobrich et al ., Nuclear Rainbow Structures in the Elastic Scattering of 16O on 16O at ELab = 350 MeV . Phys. Lett. B. 223 , 91604-3 ( 1989 ). doi: 10.1016/0370-2693(89)91604-3 http://doi.org/10.1016/0370-2693(89)91604-3 . Nuclear Rainbow Structures in the Elastic Scattering of 16O on 16O at ELab = 350 MeV H. G. Bohlen , E. Stiliaris , B. Gebauer et al ., Refractive Scattering and Reactions, Comparison of Two Systems:16O + 16O and 20Ne + 12C Z . Phys. A. 346 , BF01306079 ( 1993 ). doi: 10.1007/BF01306079 http://doi.org/10.1007/BF01306079 . Refractive Scattering and Reactions, Comparison of Two Systems:16O + 16O and 20Ne + 12C Z G. Bartnitzky , A. Blazevic , H. G. Bohlen et al ., Model-Unrestricted Nucleus-Nucleus Scattering Potentials from Measurement and Analysis of 16O+16O Scattering Phys. Lett. B. 365 , 01292-3 ( 1996 ). doi: 10.1016/0370-2693(95)01292-3 http://doi.org/10.1016/0370-2693(95)01292-3 . Model-Unrestricted Nucleus-Nucleus Scattering Potentials from Measurement and Analysis of 16O+16O Scattering Y. Sugiyama , Y. Tomita , H. Ikezoe et al ., Observation of Airy Oscillation for the 16O+16O System at ELab = 145 MeV Phys. Lett. B. 312 , 90482-W ( 1993 ). doi: 10.1016/0370-2693(93)90482-W http://doi.org/10.1016/0370-2693(93)90482-W . Observation of Airy Oscillation for the 16O+16O System at ELab = 145 MeV C.A. Bertulani and A. Gade . Nuclear Astrophysics with Radioactive Beams . Physics Reports. 485 , 09.002(2010) . doi: 10.1016/j.Physrep.2009.09.002 http://doi.org/10.1016/j.Physrep.2009.09.002 . Nuclear Astrophysics with Radioactive Beams L. H. Chien , D. T. Khoa , D. C. Cuong et al ., Consistent Mean-Field Description of the 12C+12C Optical Potential at Low Energies and the Astrophysical S-Factor . Physical Review C. 98 , 064001 ( 2018 ). doi: 10.1103/PhysRevC.98.064001 http://doi.org/10.1103/PhysRevC.98.064001 . Consistent Mean-Field Description of the 12C+12C Optical Potential at Low Energies and the Astrophysical S-Factor N. L. Anh , N. H. Phuc , D. T. Khoa , et al ., Folding Model Approach to the Elastic p+12, 13C Scattering at Low Energies and Radiative Capture 12, 13C(p,γ) Reactions . Nuclear Physics A. 1006 , 122078 ( 2021 ). doi: 10.1016/j.nuclphysa.2020.122078 http://doi.org/10.1016/j.nuclphysa.2020.122078 . Folding Model Approach to the Elastic p+12, 13C Scattering at Low Energies and Radiative Capture 12, 13C(p,γ) Reactions Both authors contributed to the study conception and design. Material preparation, data gathering and analysis were performed by Isaiah Ochala, supervised by Joseph O. Fiase. The first draft of the manuscript was written by Isaiah Ochala and both authors commented on previous versions of the manuscript. The authors jointly read and approved final manuscript. Prof. Dao T. Khoa of the Institute for Nuclear Science and Technology (INST), Vietnam, is specifically and gratefully acknowledged for providing helpful academic materials and excellent guidance in this work. 10.1007/s41365-021-00920-z 2020-12-16 2021-05-25 2021-05-26 2021-08-05 2021-08-01 2021-09-07 © 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. 2021 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-021-00920-z http://dx.doi.org/10.1007/s41365-021-00920-z . 2021 Nuclear Science and Techniques 08 32