Transport model predictions of 225Ac production cross sections via energetic p, d and α irradiation of 232Th targets
Introduction
There is growing interest in the medical community regarding the use of α-emitting radioisotopes for the treatment of cancer and other diseases. Targeted Alpha Therapy (TAT), one of the techniques used for cancer therapy, takes advantage of the short-range and high linear energy transfer (LET) of particles emitted from radioisotopes. With a LET on the order 100 keV/μm, the emitted α particles can produce a lethal dose to the tumor cell in human tissue and effectively destroy it while minimizing the damage to the healthy cells surrounding the tumor. Actinium-225 (t1/2 = 9.92 d) and its decay chain product 213Bi (t1/2 = 45.59 m) are strong therapeutic candidates for TAT. With the increasing demand for 225Ac for clinical trials and the limited supply of this isotope, new production methods for 225Ac need to be explored.
Currently, the most common method to produce 225Ac is by extracting it from 229Th (t1/2 = 7,880 y), where the 229Th is produced by chemical separation from a supply of 233U (t1/2 = 1.6 × 105 y). Uranium-233 was produced in 1960's and 1970's as part of the molten salt breeder reactor program in the United States. Currently, the 233U is in long-term storage with no immediate or anticipated need. The three main world suppliers of 225Ac are Oak Ridge National Laboratory (ORNL), United States, Institute of Physics and Power Engineering (IPPE), Russia, and the Institute of Transuranium Elements (ITE), Germany, with a combined annual supply of 1–2 Ci of 225Ac (IAEA, 2013). In addition, the Canadian Nuclear Laboratory in Chalk River, Canada, has recently constructed a 10 mCi 229Th/225Ac generator (IAEA, 2018). The 225Ac produced from the decay of 229Th is considered as a “carrier-free” nuclide and supports a limited number of clinical studies. With the increase in clinical studies and potential widespread use of this isotope, the need for additional production channels has become apparent. An alternate approach to produce 225Ac is by high-energy proton spallation of 232Th targets (Duchemin et al., 2015; Ermolaev et al., 2012; Griswold et al., 2018; Robertson et al., 2019; Weidner et al., 2012a, 2012b). After irradiation, 225Ac is chemically separated from the thorium target and other product nuclei like radium, protactinium, as well as other spallation, and fission products. In this production route, the challenge is the presence of 226Ac (t1/2 = 29.37 h) and 227Ac (t1/2 = 21.8 y) contaminants because they cannot be chemically separated from 225Ac. Given its short 29-h half-life, 226Ac has little consequence on the use of 225Ac in TAT.
Because of its much longer half-life 227Ac (and its decay products) is a critical contaminate in 225Ac preparation. Consequently, any methodology, such as those evaluated in our study, which results in co-production of 227Ac with 225Ac must be evaluated for the amount of 227Ac that is produced in order to determine the viability of the method as a potential production mode.
In this paper, we describe the prediction of reaction cross sections for the interaction of energetic protons, deuterons, and α-particles with a232Th target resulting in formation of 225,227Ac, 227,229Th, 223,225Ra, and 229,230,231Pa. We used Monte Carlo particle transport codes PHITS and MCNP6 to calculate the cross sections at incident projectile energies between 10 and 800 MeV/nucleon. The predicted production cross sections were compared with earlier measured data and predictions when available.
Section snippets
PHITS Monte Carlo simulations
The Particle and Heavy Ion Transport code System (PHITS) (Sato et al., 2018) is a Monte Carlo particle transport simulation code developed by the Japan Atomic Energy Agency. PHITS transports all particles over a wide range of energies, and users have the option to select the physics models to be used. Starting with version 3.00 and up, the KUROTAMA model (Iida et al., 2007; Sihver et al., 2014) became the default model for nucleus-nucleus reaction cross sections instead of the NASA model (
MCNP6 Monte Carlo simulations
The Monte Carlo N-Particle (MCNP) transport code, version 6 (MCNP6) (Goorley et al., 2012), was used to perform the same set of simulations that were performed in the previous section, with the same input and with the recommended defaults contained in the Cascade-Excitation Model (CEM03.03). The interaction is first simulated with the Intra Nuclear Cascade (INC) model, and depending on what reaction products are created, the reaction progresses using a choice of several models, including
Results
The yield of a particular product nuclei can occur through several direct and indirect (decay of short-lived precursors) reaction channels. For example, the production of 229Th via proton bombardment of a232Th target can be the result of direct channels like the 232Th(p,p3n)229Th, 232Th(p,d2n)229Th and 232Th(p,tn)229Th reactions. It can also be the result of indirect channels such as the 232Th(p,4n)229Pa(t1/2 = 1.5 d, EC)229Th and 232Th(p,α)229Ac(t1/2 = 62.7 m, β−)229Th reactions. Although
Discussion
In this work, we are investigating the most effective charged particle beam for the production of the reported isotopes. Proton, deuteron, and α projectiles are used for a selected number of isotopes. Fig. 12, Fig. 13 demonstrate the independent cross sections of 225Ac and 229Th from all three projectile species. The PHITS-K model was employed for comparison of these calculated cross section results. Furthermore, with the change of the default nucleus-nucleus reaction cross section model in the
Conclusion
Overall, PHITS and MCNP6 calculations show reasonable agreement with measured independent and cumulative cross sections for the product nuclei reported in this work. The results presented in this work indicate that these transport model calculations are useful tools for predicting and optimizing accelerator production of 229Th and 225Ac. The degree of the codes' abilities to match the measured data provides a qualitative assessment of the codes’ abilities to predict data from similar, but
Author statement
Naser Burahmah: Writing- Original draft, Investigation, Conceptualization, Software, preparation, Visualization Justin Griswold: Investigation, Software, Validation, Writing- Reviewing and Editing. Lawrence Heilbronn: Supervision, Writing- Reviewing and Editing, Investigation, Validation. Saed Mirzadeh: Validation, Writing- Reviewing and Editing.
Declaration of competing interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgments
This research is supported by the U.S. Department of Energy Isotope Program, managed by the Office of Science for isotope R&D and production. This manuscript has been authored by UT-Battelle, LLC, under contract DE-AC05-00OR22725 and DE-SC0020140 with the US Department of Energy (DOE). The US government retains and the publisher, by accepting the article for publication, acknowledges that the US government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce
References (28)
- et al.
Large scale accelerator production of (225)Ac: effective cross sections for 78-192MeV protons incident on (232)Th targets
Appl. Radiat. Isot.
(2016) - et al.
Modern nuclear data evaluation with the TALYS code system
Nucl. Data Sheets
(2012) - et al.
Re'actions nucle'aires de spallation induites sur le thorium par des protons de 150 et 82 MeV
Nucl. Phys.
(1961) - et al.
Cross-sections of the reaction 232Th(p,3n) 230Pa for production of 230U for targeted alpha therapy
Appl. Radiat. Isot.
(2008) - et al.
Comparative study of d- and 6Li-induced reactions on 232Th in terms of breakup and preequilibrium processes
Nucl. Phys. Sect. A
(1986) - et al.
Current status of the “Hybrid Kurotama model” for total reaction cross sections
Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms
(2014) - et al.
Accurate universal parameterization of absorption cross sections
Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms
(1996) - et al.
Accurate universal parameterization of absorption cross sections III – light systems
Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms
(1999) - et al.
Accurate universal parameterization of absorption cross sections II — neutron absorption cross sections
Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms
(1997) - et al.
225Ac and 223Ra production via 800MeV proton irradiation of natural thorium targets
Appl. Radiat. Isot.
(2012)
Proton-induced cross sections relevant to production of 225Ac and 223Ra in natural thorium targets below 200MeV
Appl. Radiat. Isot.
New potentialities of the Li\ège intranuclear cascade model for reactions induced by nucleons and light charged particles
Phys. Rev. C
Cross sections of232Th (p,xn+yp) reactions at Ep= 6.8 to 20.2 MeV
Phys. Scripta
Production of medical isotopes from a thorium target irradiated by light charged particles up to 70 MeV
Phys. Med. Biol.
Cited by (3)
Actinium-225 Targeted Agents: Where Are We Now?
2023, BrachytherapyPa 229 cross section measurements via deuteron irradiation of Th 232
2023, Physical Review C