Abstract:Background: Recent developments in alpha and beta emitting radionuclide therapy highlight the importance of developing efficient methods for patient-specific dosimetry. Traditional tabulated methods such as Medical Internal Radiation Dose (MIRD) estimate the dose at the organ level while more recent numerical methods based on Monte Carlo (MC) simulations are able to calculate dose at the voxel level. A precalculated MC (PMC) approach was developed in this work as an alternative to time-consuming fully simulated MC. Once the spatial distribution of alpha and beta emitters is determined using imaging and/or numerical methods, the PMC code can be used to achieve an accurate voxelized 3D distribution of the deposited energy without relying on full MC calculations. Purpose: To implement the PMC method to calculate energy deposited by alpha and beta particles emitted from the Ra-224 decay chain. Methods: The GEANT4 (version 10.7) MC toolkit was used to generate databases of precalculated tracks to be integrated in the PMC code as well as to benchmark its output. In this regard, energy spectra of alpha and beta particles emitted by the Ra-224 decay chain were generated using GAMOS (version 6.2.0) and imported into GEANT4 macro files. Either alpha or beta emitting sources were defined at the center of a homogeneous phantom filled with various materials such as soft tissue, bone, and lung where particles were emitted either mono-directionally (for database generation) or isotropically (for benchmarking). Two heterogeneous phantoms were used to demonstrate PMC code compatibility with boundary crossing events. Each precalculated database was generated step-by-step by storing particle track information from GEANT4 simulations followed by its integration in a PMC code developed in MATLAB. For a user-defined number of histories, one of the tracks in a given database was selected randomly and rotated randomly to reflect an isotropic emission. Afterward, deposited energy was divided between voxels based on step length in each voxel using a ray-tracing approach. The radial distribution of deposited energy was benchmarked against fully simulated MC calculations using GEANT4. The effect of the GEANT4 parameter StepMax on the accuracy and speed of the code was also investigated. Results: In the case of alpha decay, primary alpha particles show the highest contribution (>99%) in deposited energy compared to their secondary particles. In most cases, protons act as the main secondary particles in the deposition of energy. However, for a lung phantom, using a range cutoff parameter of 10 µm on primary alpha particles yields a higher contribution of secondary electrons than protons. Differences between deposited energy calculated by PMC and fully simulated MC are within 2% for all alpha and beta emitters in homogeneous and heterogeneous phantoms. Additionally, statistical uncertainties are less than 1% for voxels with doses higher than 5% of the maximum dose. Moreover, optimization of the parameter StepMax is necessary to achieve the best tradeoff between code accuracy and speed. Conclusions: The PMC code shows good performance for dose calculations deposited by alpha and beta emitters. As a stand-alone algorithm, it is suitable to be integrated into clinical treatment planning systems.

Calculation of the absorbed dose to 99mTc-HMPAO, 99mTc-DMSA, and 99mTc-sulfur colloid injected patients by using Monte Carlo and MIRD methods

Monte Carlo-based estimation of patient absorbed dose in 99mTc-DMSA, -MAG3, and -DTPA SPECT imaging using the University of Florida (UF) phantoms

Pre-treatment dosimetry in90Y-SIRT: Is it possible to optimise SPECT reconstruction parameters and calculation methods for accurate dosimetry?

Model-based versus specific dosimetry in diagnostic context: comparison of three dosimetric approaches

Investigation of radiological properties of imaging agents used in nuclear medicine with different methods and GATE/Geant4 simulation program

Comparison of Measured and Calculated Dose Rates Near Nuclear Medicine Patients

Development of a stand-alone precalculated Monte Carlo code to calculate the dose by alpha and beta emitters from the Ra-224 decay chain

Relevance of artefacts in99mTc-MAA SPECT scans on pre-therapy patient-specific90Y TARE internal dosimetry: a GATE Monte Carlo study

Establishing Updated Safety Standards for Independent 99mTc-MAA SPECT/CT Treatment Planning in Radioembolization

Assessment of four dose calculation algorithms using IAEA-TECDOC-1583 with medium dependency correction factor (K med ) application

An international phantom study of inter-site variability in Technetium-99m image quantification: analyses from the TARGET radioembolization study

Assessment of Radiation Dose in Nuclear Cardiovascular Imaging Using Realistic Computational Models

Patient-specific Effective Dose Estimation from Dose-Length Product in Lung Computed Tomography Using Monte Carlo Simulation

Monte Carlo CT Dose Calculation: A Comparison Between Experiment and Simulation Using ARCHER-CT

Calculations of beta radiation doses from multiwell Phosphorus-32 and Yttrium-90 applicators designed to be used in the treatment of superficial skin tumors: comparison of Monte Carlo and analytical methods

Dose Calculation Deviations Induced by Fractional Image-Guided-couch-shifts for Varian Halcyon MV Cone Beam CT

A Monte Carlo study for optimizing the detector of SPECT imaging using a XCAT human phantom

Internal absorbed dose calculation in body organs due to injection of Rhenium-188 labeled to Mu-9 antibody

Pre-treatment radiotherapy dose verification using Monte Carlo doselet modulation in a spherical phantom

SU‐FF‐T‐422: Dependency of Proton Dose on Tissue Composition Investigated Using Monte Carlo Models of Hounsfield Number Conversion and Cadaver‐Based Anatomical Data

Phantom study for 90 Y liver radioembolization dosimetry with a long axial field-of-view PET/CT