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

Following previous works modeling sources of particles with GAN, we extend the proof of concept for generating back-to-back pairs of gammas with timing information, typically for Monte Carlo simulation of PET imaging. A conditional GAN is trained once from a low statistic simulation in a given attenuation phantom and allows generating various activity source distributions. A new parameterization that improves the training is also proposed. An ideal PET reconstruction algorithm is used to evaluate the quality of the GAN. The proposed method is evaluated on NEMA phantom and CT patient image showing good agreement with reference simulations. 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.&#xD.

Modeling families of particle distributions with conditional GAN for Monte Carlo SPECT simulations

Objective. We propose a method to model families of distributions of particles exiting a phantom with a conditional Generative Adversarial Network (condGAN) during Monte Carlo simulation of SPECT imaging devices. Approach. The proposed condGAN is trained on a low statistics dataset containing the energy, the time, the position and the direction of exiting particles. In addition, it also contains a vector of conditions composed of four dimensions: the initial energy and the position of emitted particles within the phantom (a total of 12 dimensions). The information related to the gammas absorbed within the phantom is also added in the dataset. At the end of the training process, one component of the condGAN, the generator (G), is obtained. Main results. Particles with specific energies and positions of emission within the phantom can then be generated with G to replace the tracking of particle within the phantom, allowing reduced computation time compared to conventional Monte Carlo simulation. Significance. The condGAN generator is trained only once for a given phantom but can generate particles from various activity source distributions.

The OpenGATE ecosystem for Monte Carlo simulation in medical physics

This paper reviews the ecosystem of GATE, an open-source Monte Carlo toolkit for medical physics. Based on the shoulders of Geant4, the principal modules (geometry, physics, scorers) are described with brief descriptions of some key concepts (Volume, Actors, Digitizer). The main source code repositories are detailed together with the automated compilation and tests processes (Continuous Integration). We then described how the OpenGATE collaboration managed the collaborative development of about one hundred developers during almost 20 years. The impact of GATE on medical physics and cancer research is then summarized, and examples of a few key applications are given. Finally, future development perspectives are indicated.

Imaging of polychromatic sources through Compton spectral reconstruction

Objective: Study the performance of a spectral reconstruction method for Compton imaging of polychromatic sources and compare it to standard Compton reconstruction based on the selection of photopeak events. Approach: The proposed spectral and the standard photopeak reconstruction methods are used to reconstruct images from simulated sources emitting simultaneously photons of 140, 245, 364 and 511~keV. Data are simulated with perfect and realistic energy resolutions and including Doppler broadening. We compare photopeak and spectral reconstructed images both qualitatively and quantitatively by means of activity recovery coefficient and spatial resolution. Main results: The presented method allows improving the images of polychromatic sources with respect to standard reconstruction methods. The main reasons for this improvement are the increase of available statistics and the reduction of contamination from higher initial photon energies. The reconstructed images present lower noise, higher activity recovery coefficient and better spatial resolution. The improvements become more sensible as the energy resolution of the detectors decreases. Significance: Compton cameras have been studied for their capability of imaging polychromatic sources, thus allowing simultaneous imaging of multiple radiotracers. In such scenarios, Compton images are conventionally reconstructed for each emission energy independently, selecting only those measured events depositing a total energy within a fixed window around the known emission lines. We propose to employ a spectral image reconstruction method for polychromatic sources, which allows increasing the available statistics by using the information from events with partial energy deposition. The detector energy resolution influences the energy window used to select photopeak events and therefore the level of contamination by higher energies. The spectral method is expected to have a more important impact as the detector resolution worsens. In this paper we focus on energy ranges from nuclear medical imaging and we consider realistic energy resolutions.

Influence of sub-nanosecond time of flight resolution for online range verification in proton therapy using the line-cone reconstruction in Compton imaging

Online ion range monitoring in hadron therapy can be performed via detection of secondary radiation, such as promptγ-rays, emitted during treatment. The promptγemission profile is correlated with the ion depth-dose profile and can be reconstructed via Compton imaging. The line-cone reconstruction, using the intersection between the primary beam trajectory and the cone reconstructed via a Compton camera, requires negligible computation time compared to iterative algorithms. A recent report hypothesised that time of flight (TOF) based discrimination could improve the precision of theγfall-off position (FOP) measured via line-cone reconstruction, where TOF comprises both the proton transit time from the phantom entrance untilγemission, and the flight time of theγ-ray to the detector. The aim of this study was to implement such a method and investigate the influence of temporal resolution on the precision of the FOP. Monte Carlo simulations of a 160 MeV proton beam incident on a homogeneous PMMA phantom were performed using GATE. The Compton camera consisted of a silicon-based scatterer and CeBr₃scintillator absorber. The temporal resolution of the detection system (absorber + beam trigger) was varied between 0.1 and 1.3 ns rms and a TOF-based discrimination method applied to eliminate unlikely solution(s) from the line-cone reconstruction. The FOP was obtained for varying temporal resolutions and its precision obtained from its shift across 100 independentγemission profiles compared to a high statistics reference profile. The optimal temporal resolution for the given camera geometry and 10⁸primary protons was 0.2 ns where a precision of 2.30 ± 0.15 mm (1σ) on the FOP was found. This precision is comparable to current state-of-the-art Compton imaging using iterative reconstruction methods or 1D imaging with mechanically collimated devices, and satisfies the requirement of being smaller than the clinical safety margins.

Advanced Monte Carlo simulations of emission tomography imaging systems with GATE

Built on top of the Geant4 toolkit, GATE is collaboratively developed for more than 15 years to design Monte Carlo simulations of nuclear-based imaging systems. It is, in particular, used by researchers and industrials to design, optimize, understand and create innovative emission tomography systems. In this paper, we reviewed the recent developments that have been proposed to simulate modern detectors and provide a comprehensive report on imaging systems that have been simulated and evaluated in GATE. Additionally, some methodological developments that are not specific for imaging but that can improve detector modeling and provide computation time gains, such as Variance Reduction Techniques and Artificial Intelligence integration, are described and discussed.

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.

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.

Study of sensitivity and resolution for full ring PET prototypes based on continuous crystals and analytical modeling of the light distribution

Sensitivity and spatial resolution are the main parameters to maximize in the performance of a PET scanner. For this purpose, detectors consisting of a combination of continuous crystals optically coupled to segmented photodetectors have been employed. With the use of continuous crystals the sensitivity is increased with respect to the pixelated crystals. In addition, spatial resolution is no longer limited to the crystal size. The main drawback is the difficulty in determining the interaction position. In this work, we present the characterization of the performance of a full ring based on cuboid continuous crystals coupled to SiPMs. To this end, we have employed the simulations developed in a previous work for our experimental detector head. Sensitivity could be further enhanced by using tapered crystals. This enhancement is obtained by increasing the solid angle coverage, reducing the wedge-shaped gaps between contiguous detectors. The performance of the scanners based on both crystal geometries was characterized following NEMA NU 4-2008 standardized protocol in order to compare them. An average sensitivity gain over the entire axial field of view of 13.63% has been obtained with tapered geometry while similar performance of the spatial resolution has been proven with both scanners. The activity at which NECR and true peak occur is smaller and the peak value is greater for tapered crystals than for cuboid crystals. Moreover, a higher degree of homogeneity was obtained in the sensitivity map due to the tighter packing of the crystals, which reduces the gaps and results in a better recovery of homogeneous regions than for the cuboid configuration. Some of the results obtained, such as spatial resolution, depend on the interaction position estimation and may vary if other method is employed.

3D position determination in monolithic crystals coupled to SiPMs for PET

The interest in using continuous monolithic crystals in positron emission tomography (PET) has grown in the last years. Coupled to silicon photomultipliers (SiPMs), the detector can combine high sensitivity and high resolution, the two main factors to be maximized in a positron emission tomograph. In this work, the position determination capability of a detector comprised of a [Formula: see text] mm(3) LYSO crystal coupled to an [Formula: see text]-pixel array of SiPMs is evaluated. The 3D interaction position of γ-rays is estimated using an analytical model of the light distribution including reflections on the facets of the crystal. Monte Carlo simulations have been performed to evaluate different crystal reflectors and geometries. The method has been characterized and applied to different cases. Intrinsic resolution obtained with the position estimation method used in this work, applied to experimental data, achieves sub-millimetre resolution values. Average resolution over the detector surface for 5 mm thick crystal is  ∼0.9 mm FWHM and  ∼1.2 mm FWHM for 10 mm thick crystal. Depth of interaction resolution is close to 2 mm FWHM in both cases, while the FWTM is  ∼5.3 mm for 5 mm thick crystal and  ∼9.6 mm for 10 mm thick crystal.

First Images of a Three-Layer Compton Telescope Prototype for Treatment Monitoring in Hadron Therapy

A Compton telescope for dose monitoring in hadron therapy is under development at IFIC. The system consists of three layers of LaBr3 crystals coupled to silicon photomultiplier arrays. (22)Na sources have been successfully imaged reconstructing the data with an ML-EM code. Calibration and temperature stabilization are necessary for the prototype operation at low coincidence rates. A spatial resolution of 7.8 mm FWHM has been obtained in the first imaging tests.