Elastic Cross Sections for 3He + 58Ni above the Coulomb Barrier

R. Arceo^{1},
Omar Pedraza^{2},
Luis M. Sandoval^{1},
L. A. López^{2},
C. Álvarez^{1},
F. Hueyotl-Zahuantitla^{1,3},
Ariel Flores-Rosas^{1},
G. Luis Raya^{4},
Jesús Martínez-Castro^{5}

Show more

1. Introduction

We have studied the scattering of nuclei by helium and nickel atoms using the Schrödinger equation at energies up to 35 MeV for the reaction ^{3}He + ^{58}Ni using a radial Woods-Saxon potential. We treated the Schrödinger equation numerically, for the case of the Woods-Saxon potential [1] ; the parameters for this potential were adjusted to coincide with known experimental data.

The value of chi-squared was minimized by using a theoretical model and the experimental data from Fujisawa et al. [2]. The parameters thus obtained are used in the Woods-Saxon potential and we compare the results with known experimental data.

Recently we have results at low energies for the reaction ^{3,4,6He} + ^{58}Ni [2] [3] [4] [5]. We compare our results for this reaction and show that the Woods-Saxon potential agrees with known experimental data.

This paper is divided into four sections as follows. In Section 2, we briefly describe the setup; Section 2.1 is dedicated to obtaining the Woods-Saxon potential. In Section 2.2, we discuss the elastic cross section for the scattering of helium by nickel atoms. In Section 3, the obtained results are shown. Finally, in Section 4 we focus on the discussion of our results.

2. Theory

In this section we describe the procedure used to compute the Woods-Saxon potential produced by a point particle. We then calculate the cross section for the collision of two particles of mass
${m}_{\mathrm{1,2}}$ and atomic number
${Z}_{\mathrm{1,2}}$. Our approach to this problem is numerical, and we make the assumption that the interaction of the incident particle with the rest atom can be accounted for by the effective Woods-Saxon potential which we calculate below. We minimize the value of chi squared from the experimental ^{3}He + ^{58}Ni data [2] and the parameters we obtain are shown in Table 1 and Table 2.

2.1. The Woods-Saxon Potential

The Woods-Saxon potential is a mean field potential for the nucleons (protons and neutrons) inside the atomic nucleus, which is used to describe approximately the forces applied on each nucleon, in the nuclear shell model for the structure of the nucleus.

The standard Woods-Saxon potential [1], as a function of the distance ${r}^{\mathrm{*}}$ from the nuclear center, is defined by:

${V}^{\prime}\left({r}^{*}\right)=-\frac{{V}_{0}}{1+\mathrm{exp}\left(\frac{{r}^{*}-R}{a}\right)},a\ll R$ (1)

Table 1. Parameters obtained for the reaction ^{3}He + ^{58}Ni.

Table 2. Parameters obtained with the derivative in the complex term of the Woods-Saxon potential for the reaction ^{3}He + ^{58}Ni.

where ${V}_{0}$ (with dimensions of energy, MeV) represents the potential well depth, a is a length representing the “surface thickness” of the nucleus, and $R={r}_{0}{A}^{1/3}$ is the nuclear radius where ${r}_{0}=1.25\text{\hspace{0.17em}}\text{fm}$ and A is the atomic mass number.

It is interesting to examine the consequences of the radial effective Woods-Saxon potential, ${V}_{WS}\left({r}^{\mathrm{*}}\right)$, by using both real and imaginary terms in experiments such as scattering events. We do so in the following section, where we include the Coulomb interaction potential ${V}_{C}\left(r\right)$. The total radial effective potential used is

$V\left(r\right)={V}_{C}\left(r\right)+{V}_{WS}\left({r}^{*}\right),$ (2)

$V\left(r\right)=\frac{1.44{Z}_{1,2}}{r}-\frac{{V}_{0}}{1+\mathrm{exp}\left(\frac{{r}^{*}-{R}_{1}}{{a}_{1}}\right)}-\frac{i{W}_{0}}{1+\mathrm{exp}\left(\frac{{r}^{*}-{R}_{2}}{{a}_{2}}\right)}.$ (3)

In Figure 1, we show the real and imaginary parts of the Woods-Saxon potential.

2.2. The Schrödinger Equation with the Woods-Saxon Potential

In this section we solve the radial Schrödinger equation using the radial potential (Equation (3)). The Schrödinger equation is,

$\left[-\frac{{\hslash}^{2}}{2\mu}{\nabla}^{2}+V\left(r\right)\right]\Psi \left(r\right)=E\Psi \left(r\right)\mathrm{,}$ (4)

where $\mu =\frac{{m}_{1}{m}_{2}}{{m}_{1}+{m}_{2}}$ is the reduced mass for a two-particle system, E is the

energy and $V\left(r\right)$ is the radial effective potential calculated in the previous section.

We introduce $U\left(r\right)$, where

$\Psi \left(r\right)=\Psi \left(r\right)=\frac{U\left(r\right)}{r},$ (5)

and the Schrödinger Equation (4) is solved by the method of separation of variables. For the radial component we obtain

Figure 1. The solid line is for the Woods-Saxon potential using the parameters obtained from the Colorado group (Table I-a from Ref. [2] ) and the dashed line is for the Coulomb potential.

${{U}^{\u2033}}_{l}\left(r\right)+\frac{2\mu}{{\hslash}^{2}}\left[E-V\left(r\right)\right]{U}_{l}\left(r\right)-\frac{l\left(l+1\right)}{{r}^{2}}{U}_{l}\left(r\right)=0.$ (6)

The radial equation takes the final form,

$\begin{array}{l}{{U}^{\u2033}}_{l}\left(r\right)+\frac{2\mu}{{\hslash}^{2}}\left[E-\frac{1.44{Z}_{\mathrm{1,2}}}{r}+\frac{{V}_{0}}{1+\mathrm{exp}\left(\frac{{r}^{\mathrm{*}}-{R}_{1}}{{a}_{1}}\right)}+\frac{i{W}_{0}}{1+\mathrm{exp}\left(\frac{{r}^{\mathrm{*}}-{R}_{2}}{{a}_{2}}\right)}\right]{U}_{l}\left(r\right)\\ -\frac{l\left(l+1\right)}{{r}^{2}}{U}_{l}\left(r\right)=0.\end{array}$ (7)

The next step is to determine the set of the parameters for the Woods-Saxon potential. In the Table 1 and Table 2 we show the parameters obtained by minimizing the chi squared value,

${\chi}^{2}={\displaystyle \underset{i=1}{\overset{N}{\sum}}}{\left[\frac{{\sigma}_{th}\left({\theta}_{i}\right)-{\sigma}_{exp}\left({\theta}_{i}\right)}{\Delta {\sigma}_{exp}\left({\theta}_{i}\right)}\right]}^{2}\mathrm{.}$ (8)

The calculations for this analysis were done using the experimental data from Fujisawa et al. [2].

In Table 1 and Table 2 we show the parameters obtained by minimizing chi squared and using the experimental data from Fujisawa et al. [2].

The numerical techniques necessary to solve the Schrödinger equation with a radial potential are explained in chapter 3, Equation (3.28) of Ref. [6]. The solutions of ${U}_{l}$ from Equation (7) are introduced in the S matrix (Eq. 10.58 of Ref. [6] ), which is,

${S}_{l}=\frac{{U}_{l}\left({r}_{n-1}\right){r}_{n}{h}_{l}^{-}\left(k{r}_{n}\right)-{U}_{l}\left({r}_{n}\right){r}_{n-1}{h}_{l}^{-}\left(k{r}_{n-1}\right)}{{U}_{l}\left({r}_{n}\right){r}_{n-1}{h}_{l}^{+}\left(k{r}_{n-1}\right)-{U}_{l}\left({r}_{n-1}\right){r}_{n}{h}_{l}^{+}\left(k{r}_{n}\right)}\mathrm{,}$ (9)

where the S matrix is evaluated in the last two points on a mesh of size $\delta $ ( $r=0,\delta ,2\delta ,\cdots ,n\delta $ ). ${U}_{l}$ are the solutions to the Schrödinger equation with the potential previously calculated and ${h}_{l}$ are the spherical Hankel functions defined in Eq. 10.52 of Ref. [6]. The scattering amplitude for a partial wave decomposition in terms of the S matrix is,

$f\left(\theta \right)=\frac{1}{2ik}{\displaystyle \underset{l=0}{\overset{\infty}{\sum}}}\left(2l+1\right){P}_{l}\left(\mathrm{cos}\theta \right)\left({S}_{l}-1\right).$ (10)

For states with well defined spin and isospin the elastic and total cross section of nucleon-nucleon scattering into a solid angle element $\text{d}\Omega $ is given by the scattering amplitude $f\left(\theta \right)$ of the reaction

$\frac{\text{d}\sigma}{\text{d}\Omega}={\left|f\left(\theta \right)\right|}^{2}\mathrm{,}$ (11)

${\sigma}_{T}=\frac{4\pi}{k}Im\left[f\left({0}^{\circ}\right)\right]\mathrm{,}$ (12)

where k is the center-of-mass momentum and $f\left({0}^{\circ}\right)$ is the forward amplitude.

The reaction cross section is defined as the subtraction from the integral of the elastic cross section from the total cross section,

${\sigma}_{R}\equiv \frac{2\pi}{{k}^{2}}{\displaystyle \underset{\mathcal{l}=0}{\overset{\infty}{\sum}}}\left(2\mathcal{l}+1\right)Re\left(1-{S}_{\mathcal{l}}\right)-{\displaystyle \int \sigma \left(\theta \right)\text{d}\Omega}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\mathrm{.}$ (13)

Doing the integration gives

${\sigma}_{R}=\frac{2\pi}{{k}^{2}}{\displaystyle \underset{\mathcal{l}=0}{\overset{\infty}{\sum}}}\left(2\mathcal{l}+1\right)Re\left(1-{S}_{\mathcal{l}}\right)-\frac{\pi}{{k}^{2}}{\displaystyle \underset{\mathcal{l}=0}{\overset{\infty}{\sum}}}\left(2\mathcal{l}+1\right)\left[{\left|{S}_{\mathcal{l}}\right|}^{2}+1-2Re{S}_{\mathcal{l}}\right]\mathrm{,}$ (14)

${\sigma}_{R}=\frac{\pi}{{k}^{2}}{\displaystyle \underset{\mathcal{l}=0}{\overset{\infty}{\sum}}}\left(2\mathcal{l}+1\right)\left(1-{\left|{S}_{\mathcal{l}}\right|}^{2}\right)\mathrm{.}$ (15)

The results from the calculations are shown in the next section.

3. Elastic Cross Section

With the analysis performed, we proceed to evaluate numerically the Equations (11)-(12) and (15). We compare the theoretical results with experimental data for elastic cross sections for elastic scattering of helium by nickel atoms [2]. This comparison is made explicitly in Figures 2-5.

In Figures 2-4 the elastic cross section is analyzed for helium by nickel atoms. We evaluate the elastic cross section at energies from T_{Lab} = 24.15, 27.64 and 34.14 MeV considering the radial effective Woods-Saxon potential and setting the parameters to adjust the experimental points (see Table 1 and Table 2).

In Figure 5 we show the differential cross section at energies above the Coulomb barrier.

Figure 6 shows the total and reaction cross section for the interaction of helium by nickel atoms. We evaluate the cross sections at energies up to ${T}_{Lab}=35\text{\hspace{0.17em}}\text{MeV}$ considering the radial effective Woods-Saxon potential and setting the parameters to adjust the experimental points (see Table 1).

Figure 2. The differential cross section for^{ 3}He + ^{58}Ni is plotted as a function of the angle at the energy of 24.15 MeV. The solid line is for the Woods-Saxon potential, the dashed line is with the derivative in the complex term for the Woods-Saxon potential, and the dashed-dot line is for the Coulomb potential. The experimental points come from [2].

Figure 3. The differential cross section for ^{3}He + ^{58}Ni is plotted as a function of the angle at the energy of 27.64 MeV. The lines and the experimental points have the same meaning as in Figure 2.

Figure 4. The differential cross section for ^{3}He + ^{58}Ni is plotted as a function of the angle at the energy of 34.14 MeV. The lines and the experimental points have the same meaning as in Figure 2.

Figure 5. ^{3}He + ^{58}Ni elastic scattering for energies at the barrier Coulomb plotted as a function of the angular distribution. The lines are for the case of a Woods-Saxon potential.

Figure 6. The total and integrated elastic scattering cross section for ^{3}He + ^{58}Ni are plotted as a function of the incident energy. The experimental points come from [2].

4. Conclusions

In this work, we present a numerical solution of the radial Schrödinger equation using a Woods-Saxon potential. We have examined the scattering of helium atoms via nickel. The scattered from alpha particles via nickel atoms was performed and the use of an imaginary term in the Woods-Saxon potential gives a better fit to the experimental data. The parameters for the Woods-Saxon potential were varied until ${\chi}^{2}$ was minimized and they are shown in Table 1 and Table 2. The values obtained are better in comparison to those of the Colorado group (Table I-a from Ref. [2] ) at the energies of 24.15, 27.64 and 34.14 MeV.

Finally, the total cross section and integrated elastic scattering cross section are calculated and compared with experimental data. We obtain excellent agreement with the experimental data of Fujisawa et al. (Ref. [2] ) for the set of parameters obtained in Table 1 and Table 2.

Acknowledgements

This work was partially supported by FCFM-UNACH. R. Arceo acknowledges the support of the PFCE 2018.

References

[1] Woods, R.D. and Saxon, D.S. (1954) Diffuse Surface Optical Model for Nucleon-Nuclei Scattering. Physical Review, 95, 577-578.

https://doi.org/10.1103/PhysRev.95.577

[2] Fujisawa, T., Kamitsubo, H., Wada, T. and Igarashi, M. (1969) Optical-Model Analysis of Elastic Scattering of ^{3}He Particles from ^{58}Ni at 24.15, 27.64 and 34.14 MeV. Journal of the Physical Society of Japan, 27, 278-290.

https://doi.org/10.1143/JPSJ.27.278

[3] Jansen, R.H.J., de Heer, F.J., Luyken, H.J., van Wingerden, B. and Blaauw, H.J. (1976) Absolute Differential Cross Sections for Elastic Scattering of Electrons by Helium, Neon, Argon and Molecular Nitrogen. Journal of Physics B, 9, 185-212.

https://doi.org/10.1088/0022-3700/9/2/009

[4] Baron, N., Leonard, R.F. and Stewart, W.M. (1971) Alpha-Gamma Angular Correlations in ^{12}C, ^{24}Mg, ^{58}Ni, and ^{120}Sn. Physical Review C, 4, 1159-1173.

https://doi.org/10.1103/PhysRevC.4.1159

[5] Albinski, J., Budzanowski, A., Dabrowski, H., Rogalska, Z., Wiktor, S., Rebel, H., Srivastava, D.K., Alderliesten, C., Bojowald, J., Oelert, W., Mayer-B öricke, C. and Turek, P. (1985) α-Particle Scattering from Ni Isotopes at E_{α} = 172.5 MeV. Nuclear Physics A, 445, 477-494.

https://doi.org/10.1016/0375-9474(85)90453-1

[6] Gibbs, W.R. (2006) Computation in Modern Physics. 3rd Edition, World Scientific Publishing, Singapore.