Microscopic study of⁷Li-nucleus potential

The optical potential is an usual and basic tool used in the dynamic analyses of nuclear reactions. At present, most optical potentials are phenomenological. They have several parameters, and they are determined by fitting experimental data. When the experimental data are not sufficient, it is difficult to obtain a reliable phenomenological optical potential. In contrast, the microscopic optical potential (MOP) is theoretically derived from the nucleon-nucleon interaction, it has no free parameters, and does not rely on the experimental data. Therefore, to obtain optical potentials through a microscopic approach represents a goal of nuclear physics. It is of great significance for the analyses of nuclear reactions lacking experimental data.

The studies of nuclear reactions involving a light-particle projectile or ejectile are an important part of nuclear physics and highly useful for practical applications. Thus, we have already obtained the MOPs for the nucleon [1], deuteron [2], triton [3], and^3,4,6He [4-6]. In recent years, the weakly bound⁷Li induced reactions has been a topic of great interest. Breakup, complete and incomplete fusion, and some other reaction mechanisms were addressed by experimental and theoretical nuclear physicists [7,8]. The⁷Li optical potential is required in the theoretical analyses.

To date, there are several⁷Li optical potentials employed to analyze the experimental data. A semi-microscopic optical potential, whose real part is generated by the double folding model and nucleon-nucleon effective interaction and whose imaginary part is in the Woods-Saxon form, is given by Woods et al. [9] and used to analyze elastic scattering data for the¹⁵N and²⁵Mg target. Deshmukh et al. [10] provided a Wood-Saxon form optical potential, while it can only be used for¹¹⁶Sn. An optical potential provided by Camacho et al. [11] meets the dispersion relation of real and imaginary parts; however, it is only suitable for the²⁸Si target. Recently, Xu et al. [12] provided a new global phenomenological optical potential (GOP) based on the present experimental data, which is applicable to a more extensive incident energy and target region.

Since the measured⁷Li scattering data are not sufficient to date, a⁷Li MOP is obtained in the present work by folding the MOPs of its internal nucleons over their density distributions. The isospin-dependent nonrealistic nucleon MOP derived using the Green's function method in our previous work [1,13-15] is adopted to be the MOP for the constituent nucleons. The shell model is applied to construct the internal wave function and generate nucleon density distributions. The⁷Li elastic-scattering angular distributions and reaction cross-sections are calculated by the MOP and compared with the experimental data, and the results are calculated by the GOP [12].

This paper is organized as follows: the theoretical model and formulas of the MOP are presented in Sec. 2; the calculated results and analysis are provided in Sec. 3; the summary and conclusion are given in Sec. 4.

The MOP for⁷Li is generated by the folding model [16] and expressed as

where

$ U_n $ and $ U_p $ represent the MOPs for neutrons and protons, respectively. $ \rho _n $ and $ \rho _p $ are the density distributions of neutrons and protons in the ground state of⁷Li, respectively. $ \vec R $ is the relative coordinate between the centers of mass of the target and⁷Li, and $ \vec r $ is the internal coordinate of⁷Li.

The isospin-dependent nonrealistic nucleon MOP [1,13-15] is adopted to be $ U_n $ and $ U_p $ for neutrons and protons, respectively. A brief introduction for it is given as follows: From the perspective of many-body theory, the nucleon optical potential is equivalent to the mass operator of the single-particle Green's function [17]. Based on the Skyrme nucleon-nucleon effective interaction SKC16 [14], which simultaneously describes nuclear matter properties, ground state properties, and neutron-nucleus scattering, the first- and second-order mass operators of the single-particle Green's function were derived through the nuclear matter and local density approximations. The real part of the nucleon MOP was denoted by the first-order mass operator, and the imaginary part of the nucleon MOP was denoted by the imaginary part of the second-order mass operator. The incident energy of the nucleon is regarded as one seventh of the incident energy of⁷Li.

The shell model is adopted to give an appropriate nucleon density in⁷Li. Because a 1p-shell model space can efficiently describe its structure [18], and as we are here only concerned with the ground-state properties of⁷Li, the harmonic oscillator potential is adopted to describe the mean interaction for the nucleons in⁷Li, and the internal Hamiltonian of⁷Li is expressed as

where $ m $ is the nucleon mass, and $ r_i $ is coordinate of the $ i $ th nucleon in⁷Li relative to the center of mass of⁷Li. $ T_i $ represents the kinetic energy of the $ i $ th nucleon.

As the harmonic oscillator potential is used in the shell model, the ground-state wave function of⁷Li is expressed as

where $ {\cal{A}} $ is the antisymmetrization operator of the nucleons, and $ N $ is the normalization factor. $ \zeta $ represents the spin and isospin part. $ {\Phi _{\rm g.s.}} $ is determined by the parameter $ \beta = m \omega / \hbar $ , under the conditions of meeting antisymmetrization, spin, and parity ( $ I^{\pi} = 3/2^{-} $ ). On basis of the constraint condition,

a set of Jacobi coordinates is used to replace $ r_i $ and expressed as

The value of $ \beta $ is determined by

where $ \sqrt {\left\langle {r_{\rm rms}^2} \right\rangle } $ is the nuclear matter root-mean-square radius of⁷Li and set as 2.50 fm, which was obtained by fitting the reaction cross-section in Ref. [19]. It is convenient to rewrite Eq. (4) as follows:

where

Then, we can obtain a detailed expression for $ \left\langle {r_{\rm rms}^2} \right\rangle $ ,

where

$ {\left\langle {...} \right\rangle _\xi } $ means that the coordinates $ \{ \vec r_i, i = 1-7 \} $ are replaced by the Jacobi coordinates $ \{ \vec \xi _i, i = 1-6 \} $ . The formulas for $ N_i $ are the same as $ A_i $ , while $ {{\hat O}_A} $ is replaced by $ {{\hat O}_N} = 1 $ . Thus we can obtain

and therefore $ \beta $ = 0.2743 fm^-2.

$ \rho _n $ and $ \rho _p $ are defined as

where $ \tau _i $ is the isospin of $ i $ th nucleon. $ \tau _n $ and $ \tau _p $ are the isospins of the neutron and proton, respectively. It would be convenient to calculate $ \rho _i $ first, whose formula is the same as Eq. (10), while only $ {{\hat O}_A} $ is replaced by $ {{\hat O}_{\rho ,i}} = {\delta (\vec r - {{\vec r}_i})} $ . $ \rho _p = \rho _1 + \rho _2 + \rho _5 $ , $ \rho _n = \rho _3 + \rho _4 + \rho _6 + \rho _7 $ , and they have analytical expressions as

where $ a_n $ = 0.0921 fm^-3, $ a_p $ = 0.0621 fm^-3, $ b_n $ = 0.0081 fm^-5and $ b_p $ = 0.0076 fm^-5. The density distributions are plotted inFig. 1.

Figure 1.Neutron ( $\rho _n$ ) and proton ( $\rho _p$ ) density distributions.

The MOP for the⁷Li+⁵⁸Ni collision system at incident⁷Li energies of 10, 100, and 300 MeV is shown inFig. 2as an example. The depth of the real part ( $ V $ ) decreases with the increase of the radius and energy. However, the depth of the imaginary part ( $ W $ ) increases a little first and decreases as the radius increases atEL= 10, 100 MeV, while it decreases monotonously with the increase of the radius at a higher incident energy, 300 MeV. This means the contribution of $ W $ changes from the dominant surface absorption to the volume absorption as the incident energy increases. The real part of the spin-orbit potential $ V_{\rm so}{\vec s} \cdot {\vec l} $ is likewise obtained by the folding model, and $ V_{\rm so} $ is shown inFig. 2, while the imaginary part of the spin-orbit potential is omitted as it is usually very small.

Figure 2.MOP for⁷Li+⁵⁸Ni system (a) real part ( $V$ ), (b) imaginary part ( $W$ ) and (c) real part of spin-orbit potential ( $V_{\rm so}$ ).

The⁷Li elastic-scattering angular distributions and reaction cross-sections are predicted using the MOP. Comparisons are made with experimental data and the results calculated by the GOP [12].

Figure 3shows the elastic-scattering angular distribution for²⁷Al target at incident energies from 6.0 to 24.0 MeV. The result calculated by the MOP is in good agreement with experimental data [20,21] except for the underestimation atEL= 11.0 MeV for large angles. In addition, the MOP result fits the experimental data slightly better than the result calculated by the GOP [12] below 14 MeV at larger angles.

Figure 3.(color online) Calculated elastic-scattering angular distributions in Rutherford ratio for²⁷Al compared with experimental data [20,21]. Solid and dashed lines denote results calculated by MOP and the GOP [12], respectively. Graph results from top to bottom are multiplied 10⁰, 10⁻¹, 10⁻², respectively.

The calculated elastic-scattering angular distribution for⁵⁸Ni target at incident energies from 14.22 MeV to 42.0 MeV is plotted inFig. 4. The MOP efficiently reproduces the experimental data [22-24], except the slight underestimation above 70° at 16.25 MeV and 18.28 MeV, where the GOP performs slightly better.

Figure 4.(color online) Same asFig. 3for⁵⁸Ni. Experimental data are taken from Refs. [22-24]. Results from top to bottom are multiplied by 10⁰, 10⁻¹, 10⁻², respectively.

The elastic-scattering angular distributions for⁶⁵Cu at incident energy 25.0 MeV and⁸⁹Y at incident energy 60.0 MeV are shown inFig. 5. For⁶⁵Cu, the theoretical result from the MOP is lower than the measured values [25] above 70°. Reasonable agreement with the experimental data [26] is obtained for⁸⁹Y.

Figure 5.(color online) Same asFig. 3for⁶⁵Cu and⁸⁹Y. Experimental data are taken from Refs. [25,26]. The data for⁸⁹Y are multiplied by 10⁻².

Figure 6shows the elastic-scattering angular distribution for the¹¹⁶Sn target at incident energies from 18.0 MeV to 35.0 MeV. The calculated result from the MOP is in good agreement with experimental data [10] except at incident energies 22.0, 24.0, and 26.0 MeV at large angles. The MOP reproduces the measurements slightly better than the GOP at relatively lower energies.

Figure 6.(color online) Same asFig. 3for¹¹⁶Sn. Experimental data are taken from Ref. [10]. Results from top to bottom are multiplied by 10⁰, 10⁻¹, 10⁻², respectively.

The calculated elastic-scattering angular distribution for¹³⁸Ba targets is compared with experimental data [27,28] inFig. 7. When the scattering angles are less than 80°, good agreement with experimental data is obtained for the MOP. The GOP works better at larger angles.

Figure 7.(color online) Same asFig. 3for¹³⁸Ba. The experimental data are taken from Refs. [27,28]. The results from top to bottom are multiplied respectively by 10⁰, 10⁻¹, 10⁻².

InFig. 8, the calculated elastic-scattering angular distribution for²⁰⁸Pb is shown from 27.0 MeV to 52.0 MeV. The MOP result is in satisfying agreement with experimental data [29-32] and comparable to the GOP result in fitting the measured data, except for the case at 39 MeV above 70°.

Figure 8.(color online) Same asFig. 3for²⁰⁸Pb. Experimental data are taken from Refs. [29-32]. Results from top to bottom are multiplied by 10⁰, 10⁻¹, 10⁻², respectively

Figure 9shows the elastic-scattering angular distribution at some specific scattering angles for²⁷Al target. The calculated result by the MOP is slightly larger than that by the GOP at incident energies below 15 MeV and has a little better agreement with the measured values [33-36] whenEL $\leqslant $ 9 MeV. However, all calculated results from the MOP and the GOP underestimate the experimental data [37] at incident energies from 9 MeV to 11 MeV for large angles.

Figure 9.(color online) Calculated elastic-scattering angular distributions in Rutherford ratio for²⁷Al at some scattering angles compared with experimental data [33-37]. Solid and dashed lines denote results calculated by MOP and GOP [12], respectively. The results from top to bottom are multiplied by 10⁰, 10⁻¹, 10⁻², respectively.

The reaction cross-sections of⁷Li-induced reactions on¹³C,²⁷Al,²⁸Si,⁶⁴Zn,^natCu,¹¹⁶Sn,¹³⁸Ba, and²⁰⁸Pb are also calculated and shown inFig. 10andFig. 11.Figure 10presents the results for¹³C,²⁸Si,^natCu and²⁰⁸Pb. The theoretical result for¹³C is within the measurement error range [38]. The MOP result for²⁸Si is in good agreement with experimental data [39-43] below 30 MeV; however, it becomes a little larger from 90 MeV to 200 MeV. The reaction cross-section for^natCu is obtained by averaging the reaction cross-sections for⁶³Cu and⁶⁵Cu over the natural abundance. The MOP result is in good agreement with the experimental data [44] except for the energy point of 160 MeV. The MOP result for²⁰⁸Pb reproduces the experimental data [30,45] reasonably well below 70 MeV; however, it yields an underestimation at 300 MeV.Figure 11shows that the MOP reproduces experimental data for²⁷Al [21,46],⁶⁴Zn [47],¹¹⁶Sn [10], and¹³⁸Ba [28] well. The MOP results are comparable to the GOP results in fitting the measured reaction cross-sections except for²⁰⁸Pb.

Figure 10.(color online) Reaction cross-sections calculated by MOP compared with experimental data for¹³C [38],²⁸Si [39-43],^natCu [44] and²⁰⁸Pb [30,45] and results calculated by GOP [12]. Data are shifted upwards by adding 0, 2, 4 and 6 b, respectively. Solid and dashed lines denote results calculated by MOP and GOP, respectively.

Figure 11.(color online) Same asFig. 10for²⁷Al,⁶⁴Zn,¹¹⁶Sn, and¹³⁸Ba. Experimental data are taken from Refs. [10,21,28,46,47].

The application of the MOP to the prediction of⁷Li elastic scattering from light target nuclei was attempted. The elastic-scattering angular distribution for¹⁶O is calculated and compared with experimental data [48-51], as shown inFig. 12. The theoretical result from the MOP is only consistent with the magnitude of measured data in forward angles and gives an overestimation in relatively larger angles. Therefore, the MOP is not suitable for light nuclei. On the one hand, it may be interpreted that the Negele's nuclear density [52] adopted to calculate the nucleon MOP [1,13-15] is not suitable for light nuclei. On the other hand, a light nucleus, such as¹⁶O, has unique structure characteristics and reaction mechanisms [53,54], which may also lead to discrepancy between the MOP results and measured values.

Figure 12.(color online) Same asFig. 3, for¹⁶O. Experimental data are taken from Refs. [48-51]. Results from top to bottom are multiplied by 10⁰, 10^-3, 10^-6, and 10^-9, respectively.

Furthermore, some discrepancies between the calculated and measured elastic-scattering angular distributions appear at relatively larger angles, such as the case for¹³⁸Ba target at 28.0 MeV. To investigate how to improve the MOP to provide a better global agreement with experimental data, the notch perturbation method [55,56] is employed to analyze the sensitivity of the calculated elastic scattering angular distributions to the optical potential. The perturbation is performed by setting $ V $ , $ W $ , or $ V_{\rm so} $ to zero in a region of width 0.5 fm centered at radiusR. The scattering sensitivity is assessed by $ {\chi ^2}/\chi _0^2 $ , where the $ {\chi ^2} $ and $ \chi _0^2 $ are the chi-squares corresponding to the perturbed and original potentials, respectively. $ \chi _0^2 $ is calculated by

where $ {N _{\theta }} $ is the angle numbers of the experimental elastic-scattering angular distributions for⁷Li+¹³⁸Ba atEL= 28.0 MeV. $ \sigma ^{T} _0 (\theta _i) $ , $ \sigma ^{E} (\theta _i) $ and $ \Delta \sigma ^{E} (\theta _i) $ represent the theoretical value without perturbation, experimental value, and experimental error for the $ i $ th measured scattering angle, respectively. The theoretical value with perturbation, $ \sigma ^{T} (\theta _i) $ , is used to calculate $ {\chi ^2} $ in the same method.

Figure 13shows the MOP for the⁷Li+¹³⁸Ba system atEL= 28.0 MeV and $ {\chi ^2}/\chi _0^2 $ . $ {\chi ^2}/\chi _0^2 $ for $ V_{\rm so} $ of the MOP approximately remains at unity, hence it is acceptable to ignore the impact from changing $ V_{\rm so} $ and focus on only $ V $ and $ W $ . The peaks of $ {\chi ^2}/\chi _0^2 $ are located mainly in the surface interaction region 6 <R< 12 fm.

Figure 13.(color online) Notch perturbation analysis of MOP for⁷Li+¹³⁸Ba reaction atEL= 28.0 MeV, (a) $V$ of MOP; (b) $W$ of MOP; (c) radial sensitivity of elastic scattering to MOP. Solid, dashed, and dotted curves in (c) represent results of perturbing $V$ , $W$ , and $V_{\rm so}$ , respectively.

We adjust the $ V $ and $ W $ of the MOP in the sensitive region 6 <R< 12 fm by multiplying $ N_{R} $ and $ N_{I} $ , respectively, and calculate the corresponding $ {\chi ^2}/\chi _0^2 $ . It can be seen inFigure 14shows that a better agreement with experimental data at large angles is obtained when $ N_{R} $ = 0.89 and $ N_{I} $ = 1.00. This implies that a weaker real part in the surface region of the MOP may be more suitable for reproducing the measured data. In contrast, a smaller $ {\chi ^2} $ is obtained when $ N_{R} $ = 1.00 and $ N_{I} $ = 1.89, which means that a stronger imaginary part in the surface region may perform better. A correction of the MOP results is expected from the breakup effect, which is not considered in the folding model, as it onlny provides a repulsive contribution to the real part and an absorptive contribution to the imaginary part in the surface region.

Figure 14.(color online) Calculated elastic-scattering angular distributions in Rutherford ratio for¹³⁸Ba atEL= 28.0 MeV. Solid, dashed, dotted, and dash-dotted lines denote results calculated by MOP, GOP [12], MOP with adjusted $V$ and MOP with adjusted $W$ , respectively. Adjustments of $V$ and $W$ are made only within sensitive region 6 <R< 12 fm.

The⁷Li microscopic optical potential without any free parameters is obtained by the folding model. The internal wave function of⁷Li is obtained by the shell model, and a nucleon MOP base on the Skyrme nucleon-nucleon effective interaction is adopted. The reaction cross-sections and elastic-scattering angular distributions for target from²⁷Al to²⁰⁸Pb at incident energies below 450 MeV are calculated by the⁷Li microscopic optical potential. Generally, reasonable agreement with the experimental data is obtained, and the MOP is comparable to the GOP in reproducing the measurements in numerous cases. However, some discrepancies between the calculated and measured elastic-scattering angular distributions occur at relatively larger angles. The reason is analyzed, and it is found that the MOP can be improved by adding a repulsive contribution to the real part and an absorptive contribution to the imaginary part in the surface region, which may be achieved by considering the breakup effect. This will be the subject of our future study.