1.College of Physics and Electronic Science, Shanxi Datong University, Datong 037009, China 2.Key Laboratory of Nuclear Data, China Institute of Atomic Energy, P.O.BOX(275-41), Beijing 102413, China 3.College of Physics, Guangxi Normal University, Guilin 541004,China 4.Institute of Applied Physics and Computational Mathematics, Beijing 100094, China 5.Department of Physics, Nankai University, Tianjin 300071, China Received Date:2020-07-13 Available Online:2020-12-01 Abstract:We construct a new global optical model potential to describe the elastic scattering of 12C. The experimental data of elastic-scattering angular distributions and total reaction cross sections for targets from 24Mg to 209Bi are considered below 200 MeV within the framework of the optical model. The results calculated using the derived global optical potential are then compared with the existing experimental data. The reliability of the global optical potential is further tested by predicting the elastic scattering data out of the mass and energy ranges, within which the global potential parameters are determined, and reasonable results are also obtained.
HTML
--> --> --> -->
2.1.Form of the optical model potential
The optical model potential is generally composed of a real part, an imaginary part, and a Coulomb potential:
where $ V_{\rm R} $ is the real part potential, and $ W_{\rm S} $ and $ W_{\rm V} $ are the surface and volume absorption imaginary part potentials, respectively, and $ V_{\rm C}(r) $ is the Coulomb potential. As discussed in a previous work [9], only the central components have been included in the optical potential because spin-orbit and other high-order spin-dependent couplings have been shown to have little or no influence on the shape of the elastic-scattering angular distributions and reaction cross sections analyzed in the present study. The Woods-Saxon form factor is adopted in our analyses for both the real and imaginary parts of the global OMP. The real part and the surface absorption and volume absorption imaginary parts of the OMP are respectively expressed as
where $ R_{\rm i} $ and $ a_{\rm i} $ , with $ {\rm i = R, S, V} $ , are the radius parameters and diffuseness parameters corresponding to the real and surface and volume absorption imaginary parts of the potential, respectively. The radius parameter is expressed as $ R_{\rm i} = r_{\rm i}A^{1/3} $ , and A is the target mass number. The radius parameter of the real part potential can be further given by $ r_{\rm R} = r_{{\rm R}_{0}}+r_{{\rm R}_{1}}A^{1/3} $. The energy-dependent potential depths $ V_{\rm R}(E) $, $ W_{\rm S}(E) $, and $ W_{\rm V}(E) $ are parameterized as
where Z and z are the charge of the target and projectile, respectively. 22.2.Parametrization of the optical model potential -->
2.2.Parametrization of the optical model potential
Experimental data for 12C elastic scattering for nuclei from 24Mg to 209Bi below 200 MeV are analyzed in the present work. All the data were obtained from the nuclear database EXFOR at www-nds.iaea.org/exfor/exfor.htm. Details of these data are listed in Table 1.
Table1.The elastic-scattering angular distributions database for 12C projectiles. E is the incident energy for different targets in the laboratory system.
As noted in our previous works [9-11], all these elastic scattering data are fitted simultaneously using the improved APMN code [12]. The 12C OMP parameters are optimized by minimizing $ \chi^{2} $, which is defined as follows:
where N is the number of data points. $ \sigma_{i}^{\rm th} $, $ \sigma_{i}^{\rm ex} $, and $ \Delta\sigma $ are the theoretical and experimental differential cross sections and the corresponding experimental error, respectively. In searching for the parameters, we first give all the potential parameter reasonable boundaries in the varied region using some physical limitation. The uncertainty in the experimental data was also taken from the EXFOR nuclear database. Details of applying the method used here to obtain the best OMP parameters can be found in Refs. [9-11]. The final fifteen optical potential parameters, including the strength, radius, and diffuseness of the real and imaginary potentials, are presented in Table 2.
parameter
value
unit
$ V_{0} $
277.772
MeV
$ V_{1} $
?0.278
$ V_{2} $
?0.0001
MeV?1
$ W_{0} $
56.059
MeV
$ W_{1} $
?0.0546
$ U_{0} $
5.0
MeV
$ U_{1} $
0.279
$ r_{\rm R_{0}} $
1.158
fm
$ r_{\rm R_{1}} $
0.0273
fm
$ r_{\rm S} $
1.161
fm
$ r_{\rm V} $
1.627
fm
$ r_{\rm C} $
1.1
fm
$ a_{\rm R} $
0.770
fm
$ a_{\rm S} $
0.851
fm
$ a_{\rm V} $
0.545
fm
Table2.Global OMP parameters for 12C projectiles
-->
3.1.Elastic-scattering angular distributions
The elastic-scattering angular distributions of the reactions induced by 12C are calculated using the present global OMP and are systematically analyzed by comparing the theoretical results with the corresponding experimental data. For light heavy ion systems, the total mass number of the projectile and target is usually less than fifty ($ A_{ p} $ + $ A_{ T} $$ < $ 50). The behavior of the elastic-scattering angular distribution is strikingly different from that of the heavier targets at energies near the Coulomb barrier. It presents a clear oscillatory pattern because of its sensitivity to the nuclear interior. The elastic-scattering angular distributions for 24Mg targets are compared with measurements [13, 14] from the same experiment at incident energies from 16.0 to 40.0 MeV. The results fit the experimental data well below $ 100^{\circ} $ and produce the average at the backward-angle, where the data present a strong oscillatory pattern. We also predict the angular distribution of 12C + 27Al and compare it with the existing experimental data [40-43]. Similarly, agreement is found for the forward-angles, where there is no oscillatory pattern, while the results produce an average description for backward-angles. The comparison between the calculations and the corresponding experimental data for 24Mg is presented in Fig. 1. Figure1. Comparison of the calculated elastic-scattering angular distributions with the experimental data [13, 14] for 24Mg.
For the 12C + 28Si system, the OMP calculations of the elastic-scattering angular distribution are also compared with different experimental data [15-20], as shown in Fig. 2. It is clear that the experimental angular distributions are smooth below 30.0 MeV. Moreover, no oscillation is observed at the incident energies of 21.0 and 23.0 MeV, especially for backward-angles [19]. The results calculated by the global OMP thus reproduce the experimental data well. However, between 32.0 and 40.0 MeV, the experimental data for the angular distributions exhibit an upward trend above $ 60^{\circ} $, for which the calculations of the optical model do not provide a good fit [15]. In the energy range from 48.0 to 186.4 MeV, there is good agreement between the calculations and the experiment values except for individual energies, for which occasional irregular structures are apparent at some scattering angles. The elastic angular distributions between 23.0 and 31.92 MeV were also measured [21]; from 25.0 to 31.92 MeV, they exhibit strong diffraction oscillations above $ 100^{\circ} $ . It is known that results calculated using the global OMP can describe the overall pattern of the elastic cross section well, but they have no oscillating behavior in the backward-angle range, where the oscillating behavior may be explained by introducing Regge-pole form terms in the optical model S-matrix [21, 44]. Comparisons between the calculations and the experimental data are displayed in Fig. 3. It can be seen that there are some discrepancies between the calculations and the measurements. In particular, strong diffraction oscillations in the backward-angle will be studied in our next work to obtain the local OMPs of these light heavy ion systems. Figure2. Same as Fig. 1, but for 28Si [15-20].
The results of our global OMP analysis of the elastic-scattering angular distribution for 32S are shown in Fig. 4. The data in Ref. [22] are only for the angular distributions in the range from $ 72^{\circ} $ to $ 170^{\circ} $. Compared with the experimental data [22], the calculations agree well with the measurements below 20.87 MeV. At energies of 21.87 and 22.87 MeV, the experimental elastic scattering cross sections exhibit smooth oscillations, and the calculations are slightly larger than the data above $ 120^{\circ} $. Moreover, below 120°, the global OMP also gives a good description of the data [23] at 35.78 MeV. Above $ 120^{\circ} $, it only produces an average of the cross sections where the data exhibit a strong oscillation pattern. In general, a similarity among the experimental measurements is that weak and strong oscillations and occasional irregular structures are clearly exhibited for the light targets 24Mg, 27Al, 28Si, and 32S at different incident energies. The global OMP cannot reasonably describe oscillating behaviors in backward-angles since these light nuclei are considered to exhibit significant $ \alpha $ clustering. In particular, for the 12C + 24Mg system, oscillating behavior in the backward-angle is evident even at energies under the Coulomb barrier, while for the 12C + 28Si system at the same energy, the data are smooth and non-oscillating. The reason may be related to the fact that 24Mg and 28Si have quite different cluster properties [19]. For these targets, the global OMP can give a good description of the elastic-scattering angular distribution in forward-angles (below approximately $ 100^{\circ} $), although there are discrepancies between the calculations and the experimental data for backward-angles. In future work, we will focus on the oscillating behavior of these heavy-ion systems by considering other reaction mechanisms to improve the present calculations. Figure4. Same as Fig. 1, but for 32S [22, 23].
Figure 5 presents the fit of the calculations to the experimental data [25, 26] for 40Ca at incident energies of 50.96 and 180.0 MeV and compares the predictions and experimental data [26] at 300.0 and 420.0 MeV. Finally, it provides comparisons for isotopes, $ ^{42,48} $Ca, at incident energies of approximately 50 MeV. The calculations fit the data well for 40Ca at incident energies of 180.0, 300.0, and 420.0 MeV and for 48Ca at 47.13 MeV. However, a pronounced oscillation superimposed on this overall shape is observed for 40Ca at 50.96 MeV in the backward-angle range. For 42Ca, a weak oscillation can also be seen at 49.89 MeV. Nevertheless, the results calculated using the global OMP can give a reasonable fit for these targets. Figure5. Same as Fig. 1, but for $ ^{40,42,48} $Ca [25, 26].
The calculations of the elastic angular distributions for 39K, 50Cr, and 56Fe are shown in Fig. 6. This comparison indicates that a reasonable fit with the experimental data [24, 27-29] is obtained using the global OMP except for 56Fe at the incident energy of 60.0 MeV, for which the calculations are slightly larger than the measurements [28]. In Fig. 7, we compare the elastic angular distribution calculated using the obtained global OMP with the corresponding experimental data [29-31] at different incident energies for the reaction of 12C + 58Ni. A good description of the elastic scattering cross section data is obtained at the incident energies of 26.0, 27.0, 27.5, 28.0, 28.5, and 124.5 MeV, while the calculation at 60.0 MeV is smaller than the data in Ref. [31]. Figure6. Same as Fig. 1, but for 39K, 50Cr, and 56Fe [24, 27-29].
The elastic-scattering angular distributions are also calculated using the obtained global OMP for some medium mass targets. For the 12C + 90Zr system, the calculations are compared with the corresponding experimental data [26, 33, 34], and they are in good agreement with the existing experimental data at all incident energies. These results are presented in Fig. 8. Figure8. Same as Fig. 1, but for 90Zr [26, 33, 34].
For 64Ni targets, the elastic-scattering angular distributions were only measured at one energy point. Figure 9 presents the comparison between the calculations and the measurements [32] at 48.0 MeV. Good agreement on the elastic angular distribution is obtained between them. We compare the elastic angular distributions of 92Mo with the data [34] in Fig. 9. The calculations are very consistent with the corresponding experimental data at 60.0 and 90.0 MeV. Figure9. Same as Fig. 1, but for 64Ni and 92Mo [32, 34].
For heavier targets, the angular distributions are also calculated from the global OMP. Figure 10 provides comparisons of the calculations with the corresponding experimental data for 208Pb [26, 30, 37, 38] at incident energies from 54.5 to 420 MeV. Across this energy range, the global OMP gives a good approximation to these data. Similarly, for 209Bi, the agreement between the calculations and the experimental data [38, 39] is very good from 58.9 to 118.0 MeV. This result is displayed in Fig. 11. Figure10. Same as Fig. 1, but for 208Pb [26, 30, 37, 38].
The angular distributions of the elastic scattering for different targets are also measured at the same incident energy. The full curves in Fig. 12 illustrate the comparison between the global OMP calculations and the data [33, 35] at an incident energy of 66.0 MeV for the isotopic chains of $ {\rm Sn} $ and $ {\rm Zr} $. The overall agreement is also very good for these isotopes. In Fig. 13, the elastic-scattering angular distributions at incident energies of 73.5 MeV are compared with the experimental data [36] for $ ^{194,198} $Pt targets. Good agreement is again achieved. Figure12. Calculated elastic-scattering angular distributions compared with the experimental data [33, 35] at incident 12C energies of 66.0 MeV.
The elastic scattering of 12C ions by $ ^{\rm nat.} $Al, $ ^{\rm nat.} $Si, and $ ^{\rm nat.} $Ti is also investigated at a backscattering angle of $ 170^{\circ} $ below 30 MeV [45]. The experimental data are not included in the fitting data. We compare the results predicted by the global OMP with the data and obtain perfect agreement between them, as shown in Fig. 14. Figure14. Calculated elastic-scattering angular distributions at the same incident angle compared with the experimental data [38] for different targets.
For the 12C + 197Au system, the quasi-elastic scattering angular distributions are also measured in the energy range from 56 to 82 MeV [46]. A comparison between the predictions and the existing experimental data is presented in Fig. 15. It is found that the calculations agree well with the experimental values except for those at 82 MeV. Moreover, the experimental data of quasi-elastic scattering angular distributions include elastic scattering and contributions from inelastic scattering and neutron transfer. Figure15. Comparison between the predictions of the optical model and the experimental data for the quasi-elastic scattering angular distributions of 197Au [46].
The angular distributions for actinides are also predicted and compared with the corresponding experimental data. The results for 236U are shown in Fig. 16. As can be seen from the figure, a good theoretical fit to the experimental data [47] is obtained. Figure16. Comparison between the predictions of the optical model and the experimental data for the elastic scattering angular distributions of 236U [47].
23.2.Reaction cross sections -->
3.2.Reaction cross sections
The total reaction cross section is also an important observable in the optical model and is further investigated using the global OMP. The total reaction cross sections for 27Al and 28Si are predicted and compared with the existing experimental data; the predicted results are consistent with almost all of the data [48-52]. The total reaction cross sections for 40Ca and 56Fe are also calculated using the global OMP. Compared with the corresponding experimental data [26, 48], the calculations fit the data well in the error range for these targets. The results for 28Si and 40Ca are shown in Fig. 17. Figure17. Comparison between the optical model calculations and the experimental data [26, 49-52] of 12C reaction cross sections for 28Si and 40Ca.
For copper targets, the total reaction cross sections are measured for $ ^{\rm nat.} $Cu above 200 MeV [53]. We compare the theoretical results of $ ^{63,65} $Cu with the existing experimental data [53] and find that reasonable agreement is achieved between them. Similarly, the calculations for 64Zn are in good agreement with the experimental data [48]. For 90Zr and 208Pb, there are also data for the total reaction cross sections [26], which are extracted from the measured angular distributions of the elastic scattering process. The model calculations are in good agreement with these data. The results are shown in Fig. 18. Figure18. Same as Fig. 17, but for 90Zr and 208Pb [26].
To further investigate the structure of 12C, we compare the total reaction cross sections calculated for different systems. Among these systems are combinations of stable and unstable weakly bound nuclei 9Be [9] and 6Li [10], halo nuclei 6,8He [54], and $ \alpha $ cluster structure nuclei 12C and 4He [6] on the same medium mass target 90Zr. As in Ref. [55], suitable scaling of the results is achieved by dividing the cross sections by the factor $ (A_{p}^{1/3}+ A_{T}^{1/3})^{2} $ and the energy by the factor $ Z_{p}Z_{T}/(A_{p}^{1/3}+A_{T}^{1/3}) $, where $ Z_{p}(Z_{T}) $ and $ A_{p}(A_{T}) $ are the charge and mass number of the projectile (target), respectively. The results for the total reaction cross sections for the different systems are shown in Fig. 19. It is evident that the total reaction cross sections for the systems involving halo nuclei $ ^{6,8} $He are highest among the investigated systems. The reason is that the Coulomb polarization favors the neutrons in the halo residing in the region between the core and the target, which then enhances the reaction probabilities [56]. The total reaction cross sections involving the weakly bound nuclei 9Be and 6Li are higher compared to those involving the cluster projectiles 12C and 4He. The extra contribution to the total reaction cross sections must come from the weakly bound nature of the nuclei as the breakup channel is one of the dominant reaction channels [57]. For the $ \alpha $ cluster structure nuclei 12C and 4He, the total reaction cross sections are the smallest among the investigated systems. Figure19. (color online) Comparison of the total reaction cross sections for different projectiles on the medium mass target 90Zr.