Absorbed dose evaluation of Auger electron-emitting radionuclides: impact of input decay spectra on dose point kernels and S-values.

The aim of this study was to investigate the impact of decay data provided by the newly developed stochastic atomic relaxation model BrIccEmis on dose point kernels (DPKs - radial dose distribution around a unit point source) and S-values (absorbed dose per unit cumulated activity) of 14 Auger electron (AE) emitting radionuclides, namely 67Ga, 80mBr, 89Zr, 90Nb, 99mTc, 111In, 117mSn, 119Sb, 123I, 124I, 125I, 135La, 195mPt and 201Tl. Radiation spectra were based on the nuclear decay data from the medical internal radiation dose (MIRD) RADTABS program and the BrIccEmis code, assuming both an isolated-atom and condensed-phase approach. DPKs were simulated with the PENELOPE Monte Carlo (MC) code using event-by-event electron and photon transport. S-values for concentric spherical cells of various sizes were derived from these DPKs using appropriate geometric reduction factors. The number of Auger and Coster-Kronig (CK) electrons and x-ray photons released per nuclear decay (yield) from MIRD-RADTABS were consistently higher than those calculated using BrIccEmis. DPKs for the electron spectra from BrIccEmis were considerably different from MIRD-RADTABS in the first few hundred nanometres from a point source where most of the Auger electrons are stopped. S-values were, however, not significantly impacted as the differences in DPKs in the sub-micrometre dimension were quickly diminished in larger dimensions. Overestimation in the total AE energy output by MIRD-RADTABS leads to higher predicted energy deposition by AE emitting radionuclides, especially in the immediate vicinity of the decaying radionuclides. This should be taken into account when MIRD-RADTABS data are used to simulate biological damage at nanoscale dimensions.

A stochastic cascade model for Auger-electron emitting radionuclides.

To benchmark a Monte Carlo model of the Auger cascade that has been developed at the Australian National University (ANU) against the literature data. The model is applicable to any Auger-electron emitting radionuclide with nuclear structure data in the format of the Evaluated Nuclear Structure Data File (ENSDF). Schönfeld's algorithms and the BrIcc code were incorporated to obtain initial vacancy distributions due to electron capture (EC) and internal conversion (IC), respectively. Atomic transition probabilities were adopted from the Evaluated Atomic Data Library (EADL) for elements with atomic number, Z = 1-100. Atomic transition energies were evaluated using a relativistic Dirac-Fock method. An energy-restriction protocol was implemented to eliminate energetically forbidden transitions from the simulations. Calculated initial vacancy distributions and average energy spectra of 123I, 124I, and 125I were compared with the literature data. In addition, simulated kinetic energy spectra and frequency distributions of the number of emitted electrons and photons of the three iodine radionuclides are presented. Some examples of radiation spectra of individual decays are also given. Good agreement with the published data was achieved except for the outer-shell Auger and Coster-Kronig transitions. Nevertheless, the model needs to be compared with experimental data in a future study.