Transport model predictions of 225Ac production cross sections via energetic p, d and α irradiation of 232Th targets

https://doi.org/10.1016/j.apradiso.2021.109676Get rights and content

Highlights

  • Cross section measurements for proton, deuteron, and alphainduced reactions on Natural Thorium.

  • Comparison between Monte Carlo transport codes (PHITS and MCNP6) and experimental measurement have been reported.

  • Cross section measurements for 225,227Ac, 227,229Th, 223,225Ra, and 229,230,231Pa.

  • Reasonable agreement with experimental data on the systems considered.

Abstract

Monte Carlo transport codes PHITS and MCNP6 were used to calculate the production cross sections of 225,227Ac, 227,229Th, 223,225Ra, and 229,230,231Pa via the bombardment of a232Th target with energetic protons, deuterons, and α-particles. The incident projectile energies ranged between 10 and 800 MeV/nucleon. When possible, the predicted production cross sections were compared with the available experimental data and other predictions. 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 unmeasured, projectile/target systems. In addition, a comparison between calculated cross sections and data may provide insight into possible improvements in the physics models employed by those transport codes.

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)

  • J.W. Weidner et al.

    Proton-induced cross sections relevant to production of 225Ac and 223Ra in natural thorium targets below 200MeV

    Appl. Radiat. Isot.

    (2012)
  • A. Boudard et al.

    New potentialities of the Li\ège intranuclear cascade model for reactions induced by nucleons and light charged particles

    Phys. Rev. C

    (2013)
  • A. Celler et al.

    Cross sections of232Th (p,xn+yp) reactions at Ep= 6.8 to 20.2 MeV

    Phys. Scripta

    (1981)
  • C. Duchemin et al.

    Production of medical isotopes from a thorium target irradiated by light charged particles up to 70 MeV

    Phys. Med. Biol.

    (2015)
  • Cited by (3)

    View full text