A photon source model for alpha-emitter radionuclides.

Objective. A Monte Carlo virtual source model named PHID (photon from Ion decay) that generates photons emitted in the complex decay chain process of alpha-emitter radionuclides is proposed, typically for use during the simulation of SPECT image acquisition. Approach. Given an alpha-emitter radionuclide, the PHID model extracts from Geant4 databases the photon emission lines from all decaying daughters for both isometric transition and atomic relaxation processes. According to a given time range, abundances and activities in the decay chain are considered thanks to the Bateman equations, taking into account the decay rates and the initial abundances. Main results. PHID is evaluated by comparison with analog Monte Carlo simulation. It generates photons with the correct energy and temporal distribution, avoiding the costly simulation of the complete decay chain thus decreasing the computation time. The exact time gain depends on the simulation setup. As an example, it is 30× faster for simulating 1 MBq of 225 Ac in water for 1 section Moreover, for 225 Ac, PHID was also compared to a simplified source model with the two main photon emission lines (218 and 440 keV). PHID shows that 2 times more particles are simulated and 60% more counts are detected in the images. Significance. PHID can simulate any alpha-emitter radionuclide available in the Geant4 database. As a limitation, photons emitted from Bremsstrahlung are ignored, but they represent only 0.7% of the photons above 30 keV and are not significant for SPECT imaging. PHID is open-source, available in GATE 10, and eases the investigation of imaging photon emission from alpha emitters.