Scientific Research

An Academic Publisher

Calculation of Reactions Cross Section for Neutron-Induced Reactions on ^{127}I Isotope

**Author(s)**Leave a comment

^{127}I isotope was calculated using EXIFON code in the energy range of incident particle from 0 MeV to 30 MeV. The code is based on an analytical model for statistical multistep direct and multistep compound reactions (SMD/SMC model). In order to see the effect of nuclear structure on cross sections, the calculation was done using nuclear shell structure effect and without considering shell structure effect. Obtained results and statistical analysis showed that shell structure effect does not give significant changes to the cross-section at considered energy ranges. This shows that EXIFON code is a good tool for investigation of nuclear reaction cross section and is useful in the production of the radioisotopes of Iodine, Antimony and tellurium of high purity and in an efficient manner using cyclotron or nuclear reactors, these isotopes have potential application for field of medical science especially for diagnostics and therapeutic purposes.

KEYWORDS

1. Introduction

The artificially produced radioactive isotopes are important for many different applications [1] . Radioactive isotopes play an important role in the field of medical science in terms of beneficial applications in both diagnosis and therapy purposes [2] .

In radioisotope production programs, nuclear reactions data are mainly needed for optimization of production routes [3] .

At present day, radioisotope production for nuclear medicine is important because of its common use in tomography devices. Both single photon emissions computed tomography (SPECT) and positron emission tomography (PET) are used for diagnosis in nuclear medicine. In particular, the radionuclide of iodine is used for these purposes. Therefore, this radionuclide plays an important role in medical applications and research. For example, Gamma-emitted short-lived ^{123}I and long live ^{124}I [1] isotopes can be used as the diagnostic image in SPECT and PET. Besides, the ^{124}I allows for studying of important organs such as brain and heart [2] . The long-lived ^{125}I isotope is used as a source for internal radiotherapy, bone dosimetry and a biological tracer [3] . Another iodine radionuclide ^{122}I is a very short-lived isotope and used in PET for brain blood-flow studies [4] .

For the last 50 years, the International Atomic Energy Agency (IAEA) Nuclear Data Section (NDS) has been collating, compiling and reviewing nuclear data in a collection of databases and publications, aim of making these data available to a global audience to create an awareness of the wide-ranging data available in support of nuclear-related applications [4] [5] [6] . Nuclear reaction by neutron interactions on bismuth-208 was studied for the production of lead isotopes and thallium isotopes [7] .

Today, the nuclear databases are accessible online through the website provided by the IAEA Nuclear Data Section (NDS). This site offers access to tens of thousands of nuclear data sets that can be used for research, innovation, development and dissemination. By such materials and files, the authors conducted the calculation on the reaction cross-section for neutron-induced reactions on ^{127}I isotope.

2. Theoretical Background

Over the years, the concept of statistical multistep processes has become more and more important for the understanding of nuclear reaction mechanism, especially above 20 MeV. An analytical model for both statistical multistep direct (SMD) and statistical multistep compound (SMD) processes was applied for describing nuclear reactions up to 30 MeV. This can be generalized in several respects:

The extension to higher energies is performed including s-step direct processes for s = 1 up to 5.

The same residual interaction is used for computing both formation and decay of the compound nucleus within SMC as well as SMD processes. Thus, there is no reference to the optical model (OM) reaction cross section. The OM cross section for charged particles was used to simulate coulomb effects in the threshold region only [5] .

α and γ-Processes are included Spin-isospin conservation during the two-body collision is considered.

The calculation of Multiple Particle Emission (MPE) is generalized. Up to three decays of the compound nucleus are considered.

This model is formulated in detail for predicting emission spectra for neutrons, protons, alphas and photons including equilibrium, pre-equilibrium, direct as well as MPE processes in a consistent way.

Calculations are performed with one physical parameter set for several nuclei, several energies, and several reaction types.

2.1. Statistical Multistep Reaction

Statistical multistep models are very successful in describing nuclear reactions at energies up to about 100 MeV. These models enable the description of direct, pre-equilibrium, and equilibrium processes in a consistent way for a wide mass number range and various reaction channels, e.g. neutrons, protons, alpha, and gamma particles.

The application of a statistical multistep model to heavy nuclei requires the consideration of fission as a competing process to particle and “gamma-ray” emissions. Therefore, statistical multistep models should be extended to the fission channel.

In the statistical multistep model, the total emission spectrum of the process $\left(a,xb\right)$ is divided into three main parts [8] [9] ,

$\frac{\text{d}{\sigma}_{a,xb}\left({E}_{a}\right)}{\text{d}{E}_{b}}=\frac{\text{d}{\sigma}_{a,b}^{\text{SMD}}\left({E}_{a}\right)}{\text{d}{E}_{b}}+\frac{\text{d}{\sigma}_{a,b}^{\text{SMC}}\left({E}_{a}\right)}{\text{d}{E}_{b}}+\frac{\text{d}{\sigma}_{a,xb}^{\text{MPE}}\left({E}_{a}\right)}{\text{d}{E}_{b}}$ (2.1)

The first term on the right hand side of Equation (2.1) represents the statistical multistep direct (SMD) part which contains from single-step up to five-step contributions. The second term represents the statistical multistep compound (SMC) emission which is based on a master equation. Both terms together (SMD + SMC) represents the first-chance emission process [5] [7] [10] . The last term of Equation (2.1) represents the multiple particle emission (MPE) reaction which includes the second-chance, third-chance emissions, etc. These terms are summarized below:

$\frac{\text{d}{\sigma}_{a,xb}^{\text{MPE}}\left({E}_{a}\right)}{\text{d}{E}_{b}}=\underset{c}{{\displaystyle \sum}}\frac{\text{d}{\sigma}_{a,cb}\left({E}_{a}\right)}{\text{d}{E}_{b}}+\underset{c,d}{{\displaystyle \sum}}\frac{\text{d}{\sigma}_{a,cdb}\left({E}_{a}\right)}{\text{d}{E}_{b}}+\cdots $ (2.2)

2.2. Activation Cross-Sections

The following relations between the optical model (OM) reaction Cross-section and the energy-integrated partial Cross-sections should be satisfied (at each incident energy $\left({E}_{a}\right)$ )

${\sigma}_{a}^{\text{OM}}=\underset{b}{{\displaystyle \sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\sigma}_{a,b}$ (2.3)

${\sigma}_{a,b}=\underset{c}{{\displaystyle \sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\sigma}_{a,cb}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\sigma}_{a,cb}=\underset{d}{{\displaystyle \sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\sigma}_{a,cdb}$ (2.4)

with ${\sigma}_{a,b}={\sigma}_{a,b}^{\text{SMD}}+{\sigma}_{a,b}^{\text{SMC}}$ the total first-chance emission, in this context, activation Cross-sections are given by

${\sigma}_{a,b\gamma}={\sigma}_{a,b}-\underset{c\ne \gamma}{{\displaystyle \sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\sigma}_{a,cb}$ (2.5)

${\sigma}_{a,cb\gamma}={\sigma}_{a,cb}-\underset{d\ne \gamma}{{\displaystyle \sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\sigma}_{a,cbd}$ (2.6)

where $b,c,d\ne \gamma $

For example, the (n, p)-activation Cross-sections have the form

${\sigma}_{a,p\gamma}={\sigma}_{n,p}-{\sigma}_{n,pn}-{\sigma}_{n,2p}-{\sigma}_{n,p\alpha}$ (2.7)

The SMD Cross-section is a sum over s-step direct processes given by: [11]

$\frac{\text{d}{\sigma}_{a,b}^{SMD}\left({E}_{a}\right)}{\text{d}{E}_{b}}=\underset{s=1}{{\displaystyle \sum}}\frac{\text{d}{\sigma}_{a,b}^{s}\left({E}_{a}\right)}{\text{d}{E}_{b}}$ (2.8)

The SMD Cross-section has the form

$\frac{\text{d}{\sigma}_{a,b}^{\text{SMC}}\left({E}_{a}\right)}{\text{d}{E}_{b}}={\sigma}_{a}^{SMC}\left({E}_{a}\right)\underset{N={N}_{0}}{\overset{{N}^{I}}{{\displaystyle \sum}}}\frac{{\tau}_{N}\left(E\right)}{\hslash}\underset{\left(\Delta V\right)}{{\displaystyle \sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\Gamma}_{N,b}^{\left(\Delta V\right)}\left(E,{E}_{b}\right)\uparrow $ (2.9)

where ${\tau}_{N}$ satisfies the time-integrated master equation

#Math_16# (2.10)

and

${\Gamma}_{N}^{\left(\Delta v\right)}\left(E\right)\text{\hspace{0.05em}}\downarrow \text{\hspace{0.17em}}=2\text{\pi}{I}_{SS}^{2}{\rho}_{N}^{\left(\Delta v\right)}\left(E\right)$ (2.11)

The multiple particle emission is expressed as:

$\frac{\text{d}{\sigma}_{a,xb}^{\text{MPE}}\left({E}_{a}\right)}{\text{d}{E}_{b}}=\underset{c}{{\displaystyle \sum}}\frac{\text{d}{\sigma}_{a,cb}\left({E}_{a}\right)}{\text{d}{E}_{b}}+\underset{cd}{{\displaystyle \sum}}\frac{\text{d}{\sigma}_{a,cdb}\left({E}_{a}\right)}{\text{d}{E}_{b}}+\cdots $ (2.12)

To keep the model tractable, a simple two-body interaction is assumed: [5]

$I\left({r}_{1},{r}_{2}\right)=-4\text{\pi}\frac{{F}_{0}}{A}{\left[{\chi}_{nl}\left(R\right)\right]}^{-4}\delta \left({r}_{1}-{r}_{2}\right)\delta \left({r}_{1}-R\right)$ (2.13)

${F}_{0}=27.5\text{\hspace{0.17em}}\text{MeV}$ taken from nuclear structure considerations [11] .

The factor ${\left[{\chi}_{nl}\left(R\right)\right]}^{-4}$ contains the wave function at the nuclear radius $R={r}_{0}{A}^{1/3}$ .

The single-particle state density of particles $C=n,p,\alpha $ with mass ${\mu}_{c}$ is given by

$\begin{array}{c}\rho \left({E}_{c}\right)=\frac{4\text{\pi}V{\mu}_{c}{\left(2{\mu}_{c}{E}_{c}\right)}^{1/2}}{\left(2\text{\pi}\hslash \right)}\\ =\left(4.48\times {10}^{-3}\text{\hspace{0.17em}}{\text{fm}}^{-3}\cdot {\text{MeV}}^{-3/2}\right){r}_{0}^{3}A{E}_{C}^{1/2}\end{array}$ (2.14)

where $V=\frac{4\text{\pi}{R}^{3}}{3}$ is equal to the nuclear volume [9] .

The single-particle state density of bound particles (at Fermi energy) is then defined by

$g=4\rho \left({E}_{F}\right)$ (2.15)

where factor 4 considers the spin and isospin degeneracy

3. Methodology

EXIFON code is a nuclear reaction software which provides a continuous and smooth description of nuclear reactions over a wide energy and mass range which is based on an analytical model for statistical multistep direct and multistep compound reactions (SMD/SMC model). It predicts emission spectra, angular distributions, and activation Cross-sections for neutrons, protons, alpha particles, and photons. Multiple particle emissions are considered for up to three decays of the compound system. EXIFON is a fast, easy-to-handle code which predicts Cross-sections from one global parameter set. The only adjustable quantity is the pairing shift. The INPEXI code creates input files for EXIFON2.0 from mass and shell-correction tables. The MAKE6 code transforms EXIFON output into an ENDF-6 format file [9] .

The model is based on random matrix physics with the use of the Green’s function formalism [12] [13] . All calculations are performed without any free parameters. Results were presented for bombarding energies below 30 MeV [14] [15] .

3.1. Procedure

3.1.1. Nuclear Model Calculations

Theoretical calculations of Cross-section were performed by nuclear model code EXIFON the program was run and the input and output directory were defined, and then the target nucleus is specified. The incident particle and target nucleus were selected and excitation function in the general options section for this calculation was chooses.

The number of incident energy was specified followed by the first incident energy, and then the incident energy step is also specified. The Cross-section correspond to each particular energy was obtained.

The output data (OUTEXI) for the calculation was then stored in the set output directory. Also, DAT file name is stored in the set output directory.

Secondly, the option without shell effect is also used for each target nucleus, also an output data (OUTEXI) for the calculation is then stored in the set output directory. Also, DAT file is stored in the set output directory.

3.1.2. Shell Structure Effects

The shell structure effects are considered in SMC processes. Under such a situation, the single-particle state density g, in Equation (2.18) is multiplied by the factors

$\left(1+\frac{\delta W}{{E}_{X}}\left[1-\mathrm{exp}\left(-\gamma {E}_{X}\right)\right]\right)$ (2.16)

With $\gamma =0.05\text{\hspace{0.17em}}{\text{MeV}}^{-1}$ and $\delta W$ as the shell correction energy taken from tables [16] . The quantity ${E}_{X}=E\text{\hspace{0.17em}}\text{or}\text{\hspace{0.17em}}U$ denotes the excitation energy of the composite or residual systems respectively.

The calculations in this study were performed with $\left(\delta W\ne 0\right)$ and without $\left(\delta W=0\right)$ shell corrections.

4. Results and Discursions

4.1. Results

The calculated cross-section data for neutron-induced reactions on ^{127}I are given in Table 1 and Table 2.

The calculations in which the shell correction was taken into consideration are denoted by “With” on the graph’s legend, while those without the shell correction effects are denoted by “Without”.

Table 1. Cross section (mb) obtain without shell structure effect of I-127 interactions with neutron particle at different energies (MeV).

Table 2. Cross section (mb) obtain with shell structure effect of I-127 interactions with neutron particle at different energies (MeV).

4.2. Excitation Functions

Excitation function is defined as the graphical plots of cross-section against the energy of the incident particle, it is an important parameter in nuclear data analysis describes the probability that nuclear reactions can occur at particular energy of incident particle. Below are the excitation functions of the reactions.

4.3. Discursion

Neutron on ^{127}I

Iodine-127 is a naturally occurring stable element with 100% abundance isotope of iodine when neutron particle introduced into its nucleus, the nuclear reaction processes occurred and produced other stable and radioactive elements.

Figure 1 is the knockout reaction in which neutron displaced alpha particle from the nucleus of iodine-127. Figure 2 is a compound nucleus formation of iodine-128 and it emits alpha particle to produce Antimony-124.

Figure 1 and Figure 2 shown the excitation function of the interaction between neutron particle and iodine-127 nucleus for the production of Antimony-124 which is one of the radioisotopes of Antimony, it has a half-life of 60 days 21 minute and decays through beta emission and gamma particle.

Figure 3 and Figure 4 shows the excitation function of the production of Antimony-123, is a stable isotope of Antimony.

Figure 5 shown the compound nucleus formation of iodine-128, is a radioactive isotope of iodine with a half-life of twenty-five (25) minutes, it decays either by e^{−} and γs radiations. The shell correction effect was considered in both two graphs and the results show no significant different observed.

Figure 6 is an inelastic collision between the neutron and the nucleus of iodine-127, the graph shows that cross section decreases with increases in the incident energy and Figure 7 is an elastic collision between neutron and nucleus of iodine-127.

Tellurium-127 is an isotope of Tellurium with a half-life of 9.3 h, it decays by beta emission, Figure 8 and Figure 9 shown the compound nucleus formation of this important isotope. The shell correction effect was considered in both two graphs and the results show no significant different observed.

Iodine-126 is a radioactive isotope of iodine with a half-life of thirteen (13) days, it decays either by e^{−}, e^{+} and γs radiations. Figure 10 and Figure 11 shown the compound nucleus formation of this important isotope, the shell correction effect was considered in both two graphs and the results show no significant different observed.

Tellurium-126 is a stable isotope of Tellurium which can be produced through the reaction in Figure 12 and Figure 13. Figure 12 shows that no change observed for both reactions with shell effect correction and without shell effect correction. Figure 13 shows that changes were observed in the cross sections with and without shell correction.

Iodine-125 is a radioisotope of iodine with a half-life of 56 days and it decays through K-capture and gamma particle. Figure 14: Excitation function for production iodine-125.

Figure 15 is an excitation function for the production of Antimony-126 which is a radioactive isotope of Antimony with a half-life of 9 hour and it decays through beta and gamma particle.

Figure 16 is an excitation function for the production of Tellurium-125 which is a radioactive isotope of Tellurium with a half-life of 58 days and it decays

Figure 1. Excitation function of ${}^{127}I\left(n,a\right)$ reaction.

Figure 2. Excitation function of ${}^{127}I\left(n,ag\right)$ reaction.

Figure 3. Excitation function of ${}^{127}I\left(n,na\right)$ reaction.

Figure 4. Excitation function of ${}^{127}I\left(n,an\right)$ reaction.

Figure 5. Excitation function of ${}^{127}I\left(n,g\right)$ reaction.

Figure 6. Excitation function of ${}^{127}I\left(n,ng\right)$ reaction.

Figure 7. Excitation function of ${}^{127}I\left(n,n\right)$ reaction.

Figure 8. Excitation function of ${}^{127}I\left(n,pg\right)$ reaction.

Figure 9. Excitation function of ${}^{127}I\left(n,p\right)$ reaction.

Figure 10. Excitation function of ${}^{127}I\left(n,2ng\right)$ reaction.

Figure 11. Excitation function of ${}^{127}I\left(n,2n\right)$ reaction.

Figure 12. Excitation function of ${}^{127}I\left(n,pn\right)$ reaction.

Figure 13. Excitation function of ${}^{127}I\left(n,np\right)$ reaction.

Figure 14. Excitation function of ${}^{127}I\left(n,3n\right)$ reaction.

Figure 15. Excitation function of ${}^{127}I\left(n,2p\right)$ reaction.

Figure 16. Excitation function of ${}^{127}I\left(n,2np\right)$ reaction.

Table 3. The correlation between different reaction charnels of Cross section without shell structure effect.

through gamma. The shell correction effect was considered and the results show no much difference observed.

Table 3 and Table 4 gives the correlation between exit charnels of the interaction of neutron at different energies with I-127 as a target, both Table 3 and Table 4 shown that the correlation between (n, 2p) charnel and other charnels is not defined, because the value of cross sections of (n, 2p) reaction is zero throughout the energy range considered from 0 MeV to 30 MeV. This shows that this energy range is below the threshold energy for this reaction.

Table 4. The correlation between different reaction charnels of Cross section with shell structure effect.

5. Conclusions

A nuclear reaction in the intermediate-energy region is a matter of interest in some fields of science and technology such as medical radionuclide production. Radioisotope production for nuclear medicine is important because of its common use in tomography devices. Both single photon emissions computed tomography (SPECT) and positron emission tomography (PET) are used for diagnosis in nuclear medicine. In the present work, we have investigated the excitation functions for the formation of ${}^{128,127,126}\text{I}{,}^{123,124,126}\text{Sb}{,}^{127,126,125}\text{Te}$ in the interactions of alpha particle with nucleus of iodine-127 isotope for the production of some medical radioisotopes.

The shell structure effects were considered and the statistical analysis shows that no significant change in the cross-section observed at the energy range from 0 - 30 MeV considered.

This study shows that EXIFON code is a good tool for investigation of nuclear reaction cross section and this research work can be useful in the production of the radioisotopes of Iodine, Antimony and tellurium of high purity and in an efficient manner using cyclotron or nuclear reactors, these isotopes have potential application for field of medical science especially for therapeutic purposes embracing current and possible future needs.

Conflicts of Interest

The authors declare no conflicts of interest.

Cite this paper

^{127}I Isotope.

*International Journal of Medical Physics, Clinical Engineering and Radiation Oncology*,

**6**, 344-359. doi: 10.4236/ijmpcero.2017.63031.

[1] | Art, O. and Aytekin, H. (2015) Calculation of Excitation Functions of Proton, Alpha and Deuteron Induced Reactions for Production of Medical Radioisotopes 122-125I. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 345, 1-8. |

[2] |
Alharbi, A.A. and Azzam, A. (2012) Theoretical Calculations of the Reaction Cross-Sections for Proton-Induced Reactions on Natural Copper Using ALICE-IPPE Code. Applied Radiation and Isotopes, 70, 88-98. https://doi.org/10.1016/j.apradiso.2011.09.009 |

[3] |
Skakun, Y. and Qaim, S.M. (2008) Measurement of Excitation Functions of Helion-Induced Reactions on Enriched Ru Targets for Production of Medically Important 103Pd and 101mRh and Some Other Radionuclides. Applied Radiation and Isotopes, 66, 653-667. https://doi.org/10.1016/j.apradiso.2007.11.013 |

[4] |
Qaim, S.M. (2004) Use of Cyclotrons in Medicine. Radiation Physics and Chemistry, 71, 917-926. https://doi.org/10.1016/j.radphyschem.2004.04.124 |

[5] | Kalka, H. (1992) Hadrons and Nuclei Statistical Multistep Reactions from 1 to 100 MeV. Zeitschrift Physik a Hadron and Nuclei, 299, 289-299. |

[6] | Polster, D. and Kalka, H. (1991) Short Note Fission within a Statistical Multistep Model. Zeitschrift Physik a Hadron and Nuclei, 424, 423-424. |

[7] | Ahmad, I., Ibrahim, Y.Y. and Koki, F.S. (2017) Evaluation of Reaction Cross Section of Radionuclide by Particles Induced Nuclear Reactions Using EXIFON Code. Boson Journal of Modern Physics, 3, 236-244. |

[8] |
Nesaraja, B.C., Linse, K., Spellerberg, S., Sudar, S., Suhaimi, A. and Qaim, S.M. (1999) Excitation Functions of Neutron Induced Reactions on Some Isotopes of Zinc, Gallium, and Germanium in the Energy Range of 6. 2 to 12. 4 MeV. Radiochimica Acta, 9, 1-9. https://doi.org/10.1524/ract.1999.86.12.1 |

[9] | Ebiwonjumi, Y.E.C.B.F. (2014) Determination of Nuclear Reaction Cross-Sections for Neutron-Induced Reactions in Some Odd—A Nuclides. Advance in Physics Theories and Applications, 32, 55-69. |

[10] |
Rao, A.V.M. and Rao, J.R. (1991) Pre-Equilibrium Analysis of the Excitation Functions of (a, xn) Reactions on Silver and Holmium. Il Nuovo Cimento A, 104, 863-874. https://doi.org/10.1007/BF02820560 |

[11] | Muhammed, K., Onimisi, M.Y. and Jonah, S.A. (2011) Investigation of the Shell Effect on Neutron Induced Cross Section of Actinides. Journal of Nuclear and Particle Physics, 1, 6-9. |

[12] | Qaim, S.M. (1986) Recent Development in the Production of 18F, 75, 76, 77Br, and 123I. International Journal of Radiations and Instrumentation Part A. Applied Radiation and Isotopes, 37, 308-310. |

[13] |
Watanabe, Y., Fukahori, T., Kosako, K., Shigyo, N., Murata, T., Yamano, N., Hino, T., Maki, K., Nakashima, H., Odano, N. and Chiba, S. (2005) Nuclear Data Evaluations for JENDL High-Energy File. International Conference on Nuclear Data for Science and Technology, 326-331. https://doi.org/10.1063/1.1945015 |

[14] | Ford, H. (2012) Radiation Safety Information Computational Center Changes to the Rsicc Code and Lawrence Livermore National Laboratory, Liver-More, 1-13. |

[15] | White, J.E., Manneschmidt, J.B., Finch, S.Y. and Dickens, J.K. (1997) Abstracts of Computer Programs and Data Libraries Pertaining to Photon Production Data. Computational Physics and Engineering Division, 32, 101-103. |

[16] | Yamoah, S. and Asamoah, M. (2013) Calculations of Excitation Functions of (n, p), (n, a) and (n, 2n) Reaction Cross-Sections for Stable Isotopes of from Reaction Threshold to 20 MeV. Optics Communication, 3, 100-107. |

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.