Enhancement of radiation exposure risk from β-emitter radionuclides due to Internal Bremsstrahlung effect: A Monte Carlo study of 90Y case

Monte Carlo Simulations Corroborate PET-Measured Discrepancies in Activity Assessments of Commercial 90Y Vials

In a recent multicenter study, discrepancies between PET/CT-measured activity and vendor-calibrated activity for 90Y glass and resin microspheres were found. In the present work, the origin of these discrepancies was investigated by Monte Carlo (MC) simulations. Methods: Three vial configurations, containing 90Y-chloride, 90Y-labeled glass microspheres, and 90Y-labeled resin microspheres, were modeled with GAMOS, and the electric signal generated in an activity meter was simulated. Energy deposition was scored in the activity meter-active regions and converted into electric current per unit activity. Internal bremsstrahlung (IB) photons, always accompanying β-decay, were simulated in addition to 90Y decays. The electric current per source activity obtained for 90Y glass and resin microspheres, Iglass and Iresin, was compared in terms of relative percentage difference with that of 90Y-chloride ([Formula: see text] and [Formula: see text]) and each other (δ). The findings of this work were compared with the ones obtained through PET measurements in the multicenter study. Results: With the inclusion of IB photons as primary particles in MC simulations, the [Formula: see text] and [Formula: see text] results were 24.6% ± 3.9% and -15.0% ± 2.2%, respectively, whereas δ was 46.5% ± 1.9%, in very good agreement with the values reported in the multicenter study. Conclusion: The MC simulations performed in this study indicate that the discrepancies recently found between PET/CT-measured activity and vendor-calibrated activity for 90Y glass and resin microspheres can be attributed to differences in the geometry of the respective commercial vials and to the metrologic approach adopted for activity meter calibration with a 90Y-chloride liquid source. Furthermore, IB photons were shown to play a relevant role in determining the electric current in the activity meter.

Updating 90Y Voxel S-Values including internal Bremsstrahlung: Monte Carlo study and development of an analytical model

PURPOSE

Internal Bremsstrahlung (IB) is a process accompanying β-decay but neglected in Voxel S-Values (VSVs) calculation. Aims of this work were to calculate, through Monte Carlo (MC) simulation, updated 90Y-VSVs including IB, and to develop an analytical model to evaluate 90Y-VSVs for any voxel size of practical interest.

METHODS

GATE (Geant4 Application for Tomographic Emission) was employed for simulating voxelized geometries of soft tissue, with voxels sides l ranging from 2 to 6 mm, in steps of 0.5 mm. The central voxel was set as a homogeneous source of 90Y when IB photons are not modelled. For each l, the VSVs were computed for 90Y decays alone and for 90Y + IB. The analytical model was then built through fitting procedures of the VSVs including IB contribution.

RESULTS

Comparing GATE-VSVs with and without IB, differences between + 25% and + 30% were found for distances from the central voxel larger than the maximum β-range. The analytical model showed an agreement with MC simulations within ± 5% in the central voxel and in the Bremsstrahlung tails, for any l value examined, and relative differences lower than ± 40%, for other distances from the source.

CONCLUSIONS

The presented 90Y-VSVs include for the first time the contribution due to IB, thus providing a more accurate set of dosimetric factors for three-dimensional internal dosimetry of 90Y-labelled radiopharmaceuticals and medical devices. Furthermore, the analytical model constitutes an easy and fast alternative approach for 90Y-VSVs estimation for non-standard voxel dimensions.

Internal Bremsstrahlung, the missing process in beta decay Monte Carlo simulation: The relevance in 32P Dose-Point-Kernel estimation

PURPOSE

In nuclear medicine, Dose Point Kernels (DPKs), representing the energy deposited all around a point isotropic source, are extensively used for dosimetry and are usually obtained by Monte Carlo (MC) simulations. For beta-decaying nuclides, DPK is usually estimated neglecting Internal Bremsstrahlung (IB) emission, a process always accompanying the beta decay and consisting in the emission of photons having a continuous spectral distribution. This work aims to study the significance of IB emission for DPK estimation in the case of ³²P and provide DPK values corrected for the IB photon contribution.

METHODS

DPK, in terms of the scaled absorbed dose fraction, F(R/X₉₀), was first estimated by GAMOS MC simulation using the standard beta decay spectrum of ³²P, F_β(R/X₉₀). Subsequently, an additional source term accounting for IB photons and their spectral distribution was defined and used for a further MC simulation, thus evaluating the contribution of IB emission to DPK values, F_β+IB(R/X₉₀). The relative percent difference, δ, between the DPKs obtained by the two approaches, F_β+IB vs. F_β, was studied as a function of the radial distance, R.

RESULTS

As far as the energy deposition is mainly due to the beta particles, IB photons does not significantly contribute to DPK; conversely, for larger R, F_β+IB values are higher by 30-40% than F_β.

CONCLUSIONS

The inclusion of IB emission in the MC simulations for DPK estimations is recommended, as well as the use of the DPK values corrected for IB photons, here provided.

Technical note: The contribution of internal bremsstrahlung to the 90 Y dose point kernel

BACKGROUND

Internal dosimetry has an increasing role in the planning and verification of nuclear medicine therapies with radiopharmaceuticals. Dose Point Kernels (DPKs), quantifying the energy deposition all around a point source, in a homogenous medium, are extensively used for 3D dosimetry and nowadays are mostly evaluated by Monte Carlo (MC) simulation. To our knowledge, DPK for beta emitters is estimated neglecting the continuous photon emission due to the Internal Bremsstrahlung (IB), whose contribution to the absorbed dose can be relevant beyond the maximum range of betas, as evidenced in recent works.

PURPOSE

Aim of this study was to investigate and quantify, by means of MC simulations, the contribution of IB photons to DPK calculated for ⁹⁰ Y and provide the updated ⁹⁰ Y DPK.

METHODS

The overall radiation due to the decay of a ⁹⁰ Y point source, placed at the centre of concentric water shells of increasing radii from 0.02 cm to 20 cm, was simulated with GAMOS, including the IB source term whose spectral distribution was described by an analytical model. Energy deposition was scored in the shells as a function of the distance from the source, R, and DPK was estimated in terms of the scaled absorbed dose fraction, F(R/X₉₀ ), where X₉₀ is the range within which the beta particles deposit 90% of their energy.

RESULTS

A comparison between the two simulated absorbed dose distributions, calculated with or without IB, clearly shows that the latter (incomplete) choice is consistent with the findings of other Authors and systematically underestimates the absorbed dose imparted to the tissue. ⁹⁰ Y DPK values currently used are underestimated by 20%-34% for R>2X₉₀ .

CONCLUSIONS

The revised values provided in this work suggest that the inclusion of IB emission in DPK evaluations is advisable for pure beta emitters.

Relevance of Internal Bremsstrahlung photons from 90Y decay: an experimental and Monte Carlo study

Internal Bremsstrahlung (IB) is a continuous electromagnetic radiation accompanying beta decay; however, this process is not considered in radiation protection studies, particularly when estimating exposure from beta-decaying radionuclides. The aims of the present work are: i) to show that neglecting the IB process in Monte Carlo (MC) simulation leads to an underestimation of the energy deposited in a ionization chamber, in the case of a high-energy pure beta emitter such as Yttrium-90 (⁹⁰Y), and ii) to determine the most reliable choice of source term for ⁹⁰Y IB to be used in MC simulations. For this radionuclide, commonly employed in nuclear medicine and radiochemistry applications, experimental data acquired with a well ionization chamber have been compared with Monte Carlo (MC) calculations carried out in the GAMOS framework. Simulations that do not include the effect of the IB process, are found to give results underestimating the experimental values by 12-14%. Consequently, two models for the IB energy spectra, previously described by Italiano et al. [1], have been implemented using MC simulation and a good agreement has been achieved with one of them. We therefore conclude that inclusion of IB process in Monte Carlo simulation packages is advisable for a more accurate and complete treatment of electromagnetic interactions.