Towards offline PET monitoring of proton therapy at MedAustron

GPU accelerated Monte Carlo scoring of positron emitting isotopes produced during proton therapy for PET verification

Objective.Verification of delivered proton therapy treatments is essential for reaping the many benefits of the modality, with the most widely proposedin vivoverification technique being the imaging of positron emitting isotopes generated in the patient during treatment using positron emission tomography (PET). The purpose of this work is to reduce the computational resources and time required for simulation of patient activation during proton therapy using the GPU accelerated Monte Carlo code FRED, and to validate the predicted activity against the widely used Monte Carlo code GATE.Approach.We implement a continuous scoring approach for the production of positron emitting isotopes within FRED version 5.59.9. We simulate treatment plans delivered to 95 head and neck patients at Centrum Cyklotronowe Bronowice using this GPU implementation, and verify the accuracy using the Monte Carlo toolkit GATE version 9.0.Main results.We report an average reduction in computational time by a factor of 50 when using a local system with 2 GPUs as opposed to a large compute cluster utilising between 200 to 700 CPU threads, enabling simulation of patient activity within an average of 2.9 min as opposed to 146 min. All simulated plans are in good agreement across the two Monte Carlo codes. The two codes agree within a maximum of 0.95σon a voxel-by-voxel basis for the prediction of 7 different isotopes across 472 simulated fields delivered to 95 patients, with the average deviation over all fields being 6.4 × 10^-3σ.Significance.The implementation of activation calculations in the GPU accelerated Monte Carlo code FRED provides fast and reliable simulation of patient activation following proton therapy, allowing for research and development of clinical applications of range verification for this treatment modality using PET to proceed at a rapid pace.

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.

Quantitative assessment of radionuclide production yields in in-beam and offline PET measurements at different proton irradiation facilities

Objective. Reliable radionuclide production yield data are a prerequisite for positron-emission-tomography (PET) based in vivo proton treatment verification. In this context, activation data acquired at two different treatment facilities with different imaging systems were analyzed to provide experimentally determined radionuclide yields in thick targets and were compared with each other to investigate the impact of the respective imaging technique. Approach. Homogeneous thick targets (PMMA, gelatine, and graphite) were irradiated with mono-energetic proton pencil-beams at two distinct energies. Material activation was measured (i) in-beam during and after beam delivery with a double-head prototype PET camera and (ii) offline shortly after beam delivery with a commercial full-ring PET/CT scanner. Integral as well as depth-resolved β +-emitter yields were determined for the dominant positron-emitting radionuclides 11C, 15O, 13N and (in-beam only) 10C. In-beam data were used to investigate the qualitative impact of different monitoring time schemes on activity depth profiles and their quantitative impact on count rates and total activity. Main results. Production yields measured with the in-beam camera were comparable to or higher compared to respective offline results. Depth profiles of radionuclide-specific yields obtained from the double-head camera showed qualitative differences to data acquired with the full-ring camera with a more convex profile shape. Considerable impact of the imaging timing scheme on the activity profile was observed for gelatine only with a range variation of up to 3.5 mm. Evaluation of the coincidence rate and the total number of observed events in the considered workflows confirmed a strongly decreasing rate in targets with a large oxygen fraction. Significance. The observed quantitative and qualitative differences between the datasets underline the importance of a thorough system commissioning. Due to the lack of reliable cross-section data, in-house phantom measurements are still considered a gold standard for careful characterization of the system response and to ensure a reliable beam range verification.

Measurement of biological washout rates depending on tumor vascular status in 15O in-beam rat-PET

Objective. The biological washout of positron emitters should be modeled and corrected in order to achieve quantitative dose range verification in charged particle therapy based on positron emission tomography (PET). This biological washout effect is affected by physiological environmental conditions such as blood perfusion and metabolism, but the correlation to tumour pathology has not been studied yet. Approach. The aim of this study was to investigate the dependence of the biological washout rate on tumour vascular status in rat irradiation. Two types of tumour vascularity conditions, perfused and hypoxic, were modelled with nude rats. The rats were irradiated by a radioactive 15O ion beam and time activity curves were acquired by dynamic in-beam PET measurement. Tumour tissue sections were obtained to observe the histology as well. The biological washout rate was derived using a single-compartment model with two decay components (medium decay, k 2m and slow decay, k 2s ). Main results. All k 2m values in the vascular perfused tumour tissue were higher than the values of the normal tissue. All k 2m values in the hypoxic tumour tissue were much lower than the values of the vascular perfused tumour tissue and slightly lower than the values of the normal tissue. Significance. The dependency of the biological washout on the tumour vasculature conditions was experimentally shown.

The effects of Compton camera data acquisition and readout timing on PG imaging for proton range verification

The purpose of this study was to determine how the characteristics of the data acquisition (DAQ) electronics of a Compton camera (CC) affect the quality of the recorded prompt gamma (PG) interaction data and the reconstructed images, during clinical proton beam delivery. We used the Monte-Carlo-plus-Detector-Effect (MCDE) model to simulate the delivery of a 150 MeV clinical proton pencil beam to a tissue-equivalent plastic phantom. With the MCDE model we analyzed how the recorded PG interaction data changed as two characteristics of the DAQ electronics of a CC were changed: (1) the number of data readout channels; and (2) the active charge collection, readout, and reset time. As the proton beam dose rate increased, the number of recorded PG single-, double-, and triple-scatter events decreased by a factor of 60× for the current DAQ configuration of the CC. However, as the DAQ readout channels were increased and the readout/reset timing decreased, the number of recorded events decreased by <5× at the highest clinical dose rate. The increased number of readout channels and reduced readout/reset timing also resulted in higher quality recorded data. That is, a higher percentage of the recorded double- and triple-scatters were "true" events (caused by a single incident gamma) and not "false" events (caused by multiple incident gammas). The increase in the number and the quality of recorded data allowed higher quality PG images to be reconstructed even at the highest clinical dose rates.

Clinical commissioning of intensity-modulated proton therapy systems: Report of AAPM Task Group 185

Proton therapy is an expanding radiotherapy modality in the United States and worldwide. With the number of proton therapy centers treating patients increasing, so does the need for consistent, high-quality clinical commissioning practices. Clinical commissioning encompasses the entire proton therapy system's multiple components, including the treatment delivery system, the patient positioning system, and the image-guided radiotherapy components. Also included in the commissioning process are the x-ray computed tomography scanner calibration for proton stopping power, the radiotherapy treatment planning system, and corresponding portions of the treatment management system. This commissioning report focuses exclusively on intensity-modulated scanning systems, presenting details of how to perform the commissioning of the proton therapy and ancillary systems, including the required proton beam measurements, treatment planning system dose modeling, and the equipment needed.

Secondary Particle Interactions in a Compton Camera Designed for in vivo Range Verification of Proton Therapy

The purpose of this study was to determine the types, proportions, and energies of secondary particle interactions in a Compton camera (CC) during the delivery of clinical proton beams. The delivery of clinical proton pencil beams ranging from 70 to 200 MeV incident on a water phantom was simulated using Geant4 software (version 10.4). The simulation included a CC similar to the configuration of a Polaris J3 CC designed to image prompt gammas (PGs) emitted during proton beam irradiation for the purpose of in vivo range verification. The interaction positions and energies of secondary particles in each CC detector module were scored. For a 150-MeV proton beam, a total of 156,688(575) secondary particles per 10⁸ protons, primarily composed of gamma rays (46.31%), neutrons (41.37%), and electrons (8.88%), were found to reach the camera modules, and 79.37% of these particles interacted with the modules. Strategies for using CCs for proton range verification should include methods of reducing the large neutron backgrounds and low-energy non-PG radiation. The proportions of interaction types by module from this study may provide information useful for background suppression.