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 of225Ac in water for 1 section Moreover, for225Ac, 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.

Annihilation photon GAN source model for PET Monte Carlo simulation.

Objective.Following previous works on virtual sources model with Generative Adversarial Network (GAN), we extend the proof of concept for generating back-to-back pairs of gammas with timing information, typically for Monte Carlo simulation of Positron Emission Tomography(PET) imaging.Approach.A conditional GAN is trained once from a low statistic simulation in a given attenuation phantom and enables the generation of various activity source distributions. GAN training input is a set of gammas exiting a phantom, tracked from a source of positron emitters, described by position, direction and energy. A new parameterization that improves the training is also proposed. An ideal PET reconstruction algorithm is used to evaluate the quality of the GAN.Main results.The proposed method is evaluated on National Electrical Manufacturers Association (NEMA) International Electrotechnical Commission (IEC) phantoms and with CT patient image showing good agreement with reference simulations. The proportions of 2-gammas, 1-gammas and absorbed-gammas are respected to within one percent, image profiles matched and recovery coefficients were close with less than 5% difference. GAN tends to blur gamma energy peak, e.g. 511 keV.Significance.Once trained, the GAN generator can be used as input source for Monte Carlo simulations of PET imaging systems, decreasing the computational time with speedups up to ×400 according to the configurations.

Modeling complex particles phase space with GAN for Monte Carlo SPECT simulations: a proof of concept.

A method is proposed to model by a generative adversarial network the distribution of particles exiting a patient during Monte Carlo simulation of emission tomography imaging devices. The resulting compact neural network is then able to generate particles exiting the patient, going towards the detectors, avoiding costly particle tracking within the patient. As a proof of concept, the method is evaluated for single photon emission computed tomography (SPECT) imaging and combined with another neural network modeling the detector response function (ARF-nn). A complete rotating SPECT acquisition can be simulated with reduced computation time compared to conventional Monte Carlo simulation. It also allows the user to perform simulations with several imaging systems or parameters, which is useful for imaging system design.

MACACO II test-beam with high energy photons.

The IRIS group at IFIC Valencia is developing a three-layer Compton camera for treatment monitoring in proton therapy. The system is composed of three detector planes, each made of a [Formula: see text] monolithic crystal coupled to a SiPM array. Having obtained successful results with the first prototype (MACACO) that demonstrated the feasibility of the proposed technology, a second prototype (MACACO II) with improved performance has been developed, and is the subject of this work. The new system has an enhanced detector energy resolution which translates into a higher spatial resolution of the telescope. The image reconstruction method has also been improved with an accurate model of the sensitivity matrix. The device has been tested with high energy photons at the National Accelerator Centre (CNA, Seville). The tests involved a proton beam of 18 MeV impinging on a graphite target, to produce 4.4 MeV photons. Data were taken at different system positions of the telescope with the first detector at 65 and 160 mm from the target, and at different beam intensities. The measurements allowed successful reconstruction of the photon emission distribution at two target positions separated by 5 mm in different telescope configurations. This result was obtained both with data recorded in the first and second telescope planes (two interaction events) and, for the first time in beam experiments, with data recorded in the three planes (three interaction events).

Image reconstruction for a multi-layer Compton telescope: an analytical model for three interaction events.

Compton Cameras are electronically collimated photon imagers suitable for sub-MeV to few MeV gamma-ray detection. Such features are desirable to enable in vivo range verification in hadron therapy, through the detection of secondary Prompt Gammas. A major concern with this technique is the poor image quality obtained when the incoming gamma-ray energy is unknown. Compton Cameras with more than two detector planes (multi-layer Compton Cameras) have been proposed as a solution, given that these devices incorporate more signal sequences of interactions to the conventional two interaction events. In particular, three interaction events convey more spectral information as they allow inferring directly the incident gamma-ray energy. A three-layer Compton Telescope based on continuous Lanthanum (III) Bromide crystals coupled to Silicon Photomultipliers is being developed at the IRIS group of IFIC-Valencia. In a previous work we proposed a spectral reconstruction algorithm for two interaction events based on an analytical model for the formation of the signal. To fully exploit the capabilities of our prototype, we present here an extension of the model for three interaction events. Analytical expressions of the sensitivity and the System Matrix are derived and validated against Monte Carlo simulations. Implemented in a List Mode Maximum Likelihood Expectation Maximization algorithm, the proposed model allows us to obtain four-dimensional (energy and position) images by using exclusively three interaction events. We are able to recover the correct spectrum and spatial distribution of gamma-ray sources when ideal data are employed. However, the uncertainties associated to experimental measurements result in a degradation when real data from complex structures are employed. Incorrect estimation of the incident gamma-ray interaction positions, and missing deposited energy associated with escaping secondaries, have been identified as the causes of such degradation by means of a detailed Monte Carlo study. As expected, our current experimental resolution and efficiency to three interaction events prevents us from correctly recovering complex structures of radioactive sources. However, given the better spectral information conveyed by three interaction events, we expect an improvement of the image quality of conventional Compton imaging when including such events. In this regard, future development includes the incorporation of the model assessed in this work to the two interaction events model in order to allow using simultaneously two and three interaction events in the image reconstruction.

CCMod: a GATE module for Compton camera imaging simulation.

Compton cameras are gamma-ray imaging systems which have been proposed for a wide variety of applications such as medical imaging, nuclear decommissioning or homeland security. In the design and optimization of such a system Monte Carlo simulations play an essential role. In this work, we propose a generic module to perform Monte Carlo simulations and analyses of Compton Camera imaging which is included in the open-source GATE/Geant4 platform. Several digitization stages have been implemented within the module to mimic the performance of the most commonly employed detectors (e.g. monolithic blocks, pixelated scintillator crystals, strip detectors...). Time coincidence sorter and sequence coincidence reconstruction are also available in order to aim at providing modules to facilitate the comparison and reproduction of the data taken with different prototypes. All processing steps may be performed during the simulation (on-the-fly mode) or as a post-process of the output files (offline mode). The predictions of the module have been compared with experimental data in terms of energy spectra, angular resolution, efficiency and back-projection image reconstruction. Consistent results within a 3-sigma interval were obtained for the energy spectra except for low energies where small differences arise. The angular resolution measure for incident photons of 1275 keV was also in good agreement between both data sets with a value close to 13°. Moreover, with the aim of demonstrating the versatility of such a tool the performance of two different Compton camera designs was evaluated and compared.