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.

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.

Mid-position treatment strategy for locally advanced lung cancer: a dosimetric study

OBJECTIVE

The internal target volume (ITV) strategy generates larger planning target volumes (PTVs) in locally advanced non-small cell lung cancer (LA-NSCLC) than the Mid-position (Mid-p) strategy. We investigated the benefit of the Mid-p strategy regarding PTV reduction and dose to the organs at risk (OARs).

METHODS

44 patients with LA-NSCLC were included in a randomized clinical study to compare ITV and Mid-p strategies. GTV were delineated by a physician on maximum intensity projection images and on Mid-p images from four-dimensional CTs. CTVs were obtained by adding 6 mm uniform margin for microscopic extension. CTV to PTV margins were calculated using the van Herk's recipe for setup and delineation errors. For the Mid-p strategy, the mean target motion amplitude was added as a random error. For both strategies, three-dimensional conformal plans delivering 60-66 Gy to PTV were performed. PTVs, dose-volume parameters for OARs (lung, esophagus, heart, spinal cord) were reported and compared.

RESULTS

With the Mid-p strategy, the median of volume reduction was 23.5 cm3 (p = 0.012) and 8.8 cm3 (p = 0.0083) for PTVT and PTVN respectively; the median mean lung dose reduction was 0.51 Gy (p = 0.0057). For 37.1% of the patients, delineation errors led to smaller PTV with the ITV strategy than with the Mid-p strategy.

CONCLUSION

PTV and mean lung dose were significantly reduced using the Mid-p strategy. Delineation uncertainty can unfavorably impact the advantage.

ADVANCES IN KNOWLEDGE

To the best of our knowledge, this is the first dosimetric comparison study between ITV and Mid-p strategies for LA-NSCLC.

Concepts and terms for dose/volume parameters in carbon-ion radiotherapy: Conclusions of the ULICE taskforce

PURPOSE

The Union of Light Ion Centers in Europe (ULICE) program addressed the need for uniting scientific results for carbon-ion radiation therapy obtained by several institutions worldwide in different fields of excellence, and translating them into a real benefit to the community. Particularly, the concepts for dose/volume parameters developed in photon radiotherapy cannot be extrapolated to high linear energy transfer particles.

METHODS AND MATERIALS

The ULICE-WP2 taskforce included radiation oncologists involved in carbon-ion radiation therapy and International Commission on Radiation Units and Measurements, radiation biologists, expert physicists in the fields of carbon-ion radiation therapy, microdosimetry, biological modeling and image-guided radiotherapy. Consensual reports emerged from multiple discussions within both the restricted group and the wider ULICE community. Public deliverables were produced and disseminated to the European Commission.

RESULTS

Here we highlight the disparity in practices between treating centers, then address the main topics to finally elaborate specific recommendations. Although it appears relatively simple to add geometrical margins around the clinical target volume to obtain the planning target volume as performed in photon radiotherapy, this procedure is not appropriate for carbon-ion radiation therapy. Due to the variation of the radiation quality in depth, there is no generic relative biological effectiveness value for carbon-ions outside of an isolated point, for a given fractionation and specific experimental conditions. Absorbed dose and "equieffective dose" for specified conditions must always be reported.

CONCLUSIONS

This work contributed to the development of standard operating procedures for carbon-ion radiation therapy clinical trials. These procedures are now being applied, particularly in the first phase III international, multicenter trial (PHRC Étoile).

An empirical model for calculation of the collimator contamination dose in therapeutic proton beams

Collimators are used as lateral beam shaping devices in proton therapy with passive scattering beam lines. The dose contamination due to collimator scattering can be as high as 10% of the maximum dose and influences calculation of the output factor or monitor units (MU). To date, commercial treatment planning systems generally use a zero-thickness collimator approximation ignoring edge scattering in the aperture collimator and few analytical models have been proposed to take scattering effects into account, mainly limited to the inner collimator face component. The aim of this study was to characterize and model aperture contamination by means of a fast and accurate analytical model. The entrance face collimator scatter distribution was modeled as a 3D secondary dose source. Predicted dose contaminations were compared to measurements and Monte Carlo simulations. Measurements were performed on two different proton beam lines (a fixed horizontal beam line and a gantry beam line) with divergent apertures and for several field sizes and energies. Discrepancies between analytical algorithm dose prediction and measurements were decreased from 10% to 2% using the proposed model. Gamma-index (2%/1 mm) was respected for more than 90% of pixels. The proposed analytical algorithm increases the accuracy of analytical dose calculations with reasonable computation times.

Monte Carlo simulation of prompt γ-ray emission in proton therapy using a specific track length estimator

A Monte Carlo (MC) variance reduction technique is developed for prompt-γ emitters calculations in proton therapy. Prompt-γ emitted through nuclear fragmentation reactions and exiting the patient during proton therapy could play an important role to help monitoring the treatment. However, the estimation of the number and the energy of emitted prompt-γ per primary proton with MC simulations is a slow process. In order to estimate the local distribution of prompt-γ emission in a volume of interest for a given proton beam of the treatment plan, a MC variance reduction technique based on a specific track length estimator (TLE) has been developed. First an elemental database of prompt-γ emission spectra is established in the clinical energy range of incident protons for all elements in the composition of human tissues. This database of the prompt-γ spectra is built offline with high statistics. Regarding the implementation of the prompt-γ TLE MC tally, each proton deposits along its track the expectation of the prompt-γ spectra from the database according to the proton kinetic energy and the local material composition. A detailed statistical study shows that the relative efficiency mainly depends on the geometrical distribution of the track length. Benchmarking of the proposed prompt-γ TLE MC technique with respect to an analogous MC technique is carried out. A large relative efficiency gain is reported, ca. 10(5).

A Monte Carlo pencil beam scanning model for proton treatment plan simulation using GATE/GEANT4

This work proposes a generic method for modeling scanned ion beam delivery systems, without simulation of the treatment nozzle and based exclusively on beam data library (BDL) measurements required for treatment planning systems (TPS). To this aim, new tools dedicated to treatment plan simulation were implemented in the Gate Monte Carlo platform. The method was applied to a dedicated nozzle from IBA for proton pencil beam scanning delivery. Optical and energy parameters of the system were modeled using a set of proton depth-dose profiles and spot sizes measured at 27 therapeutic energies. For further validation of the beam model, specific 2D and 3D plans were produced and then measured with appropriate dosimetric tools. Dose contributions from secondary particles produced by nuclear interactions were also investigated using field size factor experiments. Pristine Bragg peaks were reproduced with 0.7 mm range and 0.2 mm spot size accuracy. A 32 cm range spread-out Bragg peak with 10 cm modulation was reproduced with 0.8 mm range accuracy and a maximum point-to-point dose difference of less than 2%. A 2D test pattern consisting of a combination of homogeneous and high-gradient dose regions passed a 2%/2 mm gamma index comparison for 97% of the points. In conclusion, the generic modeling method proposed for scanned ion beam delivery systems was applicable to an IBA proton therapy system. The key advantage of the method is that it only requires BDL measurements of the system. The validation tests performed so far demonstrated that the beam model achieves clinical performance, paving the way for further studies toward TPS benchmarking. The method involves new sources that are available in the new Gate release V6.1 and could be further applied to other particle therapy systems delivering protons or other types of ions like carbon.