Feasibility of the J-PET to monitor the range of therapeutic proton beams.

PURPOSE: The aim of this work is to investigate the feasibility of the Jagiellonian Positron Emission Tomography (J-PET) scanner for intra-treatment proton beam range monitoring. METHODS: The Monte Carlo simulation studies with GATE and PET image reconstruction with CASToR were performed in order to compare six J-PET scanner geometries. We simulated proton irradiation of a PMMA phantom with a Single Pencil Beam (SPB) and Spread-Out Bragg Peak (SOBP) of various ranges. The sensitivity and precision of each scanner were calculated, and considering the setup's cost-effectiveness, we indicated potentially optimal geometries for the J-PET scanner prototype dedicated to the proton beam range assessment. RESULTS: The investigations indicate that the double-layer cylindrical and triple-layer double-head configurations are the most promising for clinical application. We found that the scanner sensitivity is of the order of 10-5 coincidences per primary proton, while the precision of the range assessment for both SPB and SOBP irradiation plans was found below 1 mm. Among the scanners with the same number of detector modules, the best results are found for the triple-layer dual-head geometry. The results indicate that the double-layer cylindrical and triple-layer double-head configurations are the most promising for the clinical application, CONCLUSIONS:: We performed simulation studies demonstrating that the feasibility of the J-PET detector for PET-based proton beam therapy range monitoring is possible with reasonable sensitivity and precision enabling its pre-clinical tests in the clinical proton therapy environment. Considering the sensitivity, precision and cost-effectiveness, the double-layer cylindrical and triple-layer dual-head J-PET geometry configurations seem promising for future clinical application.

Single proton LET characterization with the Timepix detector and artificial intelligence for advanced proton therapy treatment planning.

Objective.Protons have advantageous dose distributions and are increasingly used in cancer therapy. At the depth of the Bragg peak range, protons produce a mixed radiation field consisting of low- and high-linear energy transfer (LET) components, the latter of which is characterized by an increased ionization density on the microscopic scale associated with increased biological effectiveness. Prediction of the yield and LET of primary and secondary charged particles at a certain depth in the patient is performed by Monte Carlo simulations but is difficult to verify experimentally.Approach.Here, the results of measurements performed with Timepix detector in the mixed radiation field produced by a therapeutic proton beam in water are presented and compared to Monte Carlo simulations. The unique capability of the detector to perform high-resolution single particle tracking and identification enhanced by artificial intelligence allowed to resolve the particle type and measure the deposited energy of each particle comprising the mixed radiation field. Based on the collected data, biologically important physics parameters, the LET of single protons and dose-averaged LET, were computed.Main results.An accuracy over 95% was achieved for proton recognition with a developed neural network model. For recognized protons, the measured LET spectra generally agree with the results of Monte Carlo simulations. The mean difference between dose-averaged LET values obtained from measurements and simulations is 17%. We observed a broad spectrum of LET values ranging from a fraction of keVµm-1to about 10 keVµm-1for most of the measurements performed in the mixed radiation fields.Significance.It has been demonstrated that the introduced measurement method provides experimental data for validation of LETDor LET spectra in any treatment planning system. The simplicity and accessibility of the presented methodology make it easy to be translated into a clinical routine in any proton therapy facility.

Optically Stimulated Luminescent Response of the LiMgPO4 Silicone Foils to Protons and Its Dependence on Proton Energy.

Modern radiotherapy (RT) techniques, such as proton therapy, require more and more sophisticated dosimetry methods and materials. One of the newly developed technologies is based on flexible sheets made of a polymer, with the embedded optically stimulated luminescence (OSL) material in the form of powder (LiMgPO4, LMP) and a self-developed optical imaging setup. The detector properties were evaluated to study its potential application in the proton treatment plan verification for eyeball cancer. The data showed a well-known effect of lower luminescent efficiency of the LMP material response to proton energy. The efficiency parameter depends on a given material and radiation quality parameters. Therefore, the detailed knowledge of material efficiency is crucial in establishing a calibration method for detectors exposed to mixed radiation fields. Thus, in the present study, the prototype of the LMP-based silicone foil material was tested with monoenergetic uniform proton beams of various initial kinetic energies constituting the so-called spread-out Bragg peak (SOBP). The irradiation geometry was also modelled using the Monte Carlo particle transport codes. Several beam quality parameters, including dose and the kinetic energy spectrum, were scored. Finally, the obtained results were used to correct the relative luminescence efficiency response of the LMP foils for monoenergetic and spread-out proton beams.

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.

Quantification of biological range uncertainties in patients treated at the Krakow proton therapy centre.

BACKGROUND: Variable relative biological effectiveness (vRBE) in proton therapy might significantly modify the prediction of RBE-weighted dose delivered to a patient during proton therapy. In this study we will present a method to quantify the biological range extension of the proton beam, which results from the application of vRBE approach in RBE-weighted dose calculation. METHODS AND MATERIALS: The treatment plans of 95 patients (brain and skull base patients) were used for RBE-weighted dose calculation with constant and the McNamara RBE model. For this purpose the Monte Carlo tool FRED was used. The RBE-weighted dose distributions were analysed using indices from dose-volume histograms. We used the volumes receiving at least 95% of the prescribed dose (V95) to estimate the biological range extension resulting from vRBE approach. RESULTS: The vRBE model shows higher median value of relative deposited dose and D95 in the planning target volume by around 1% for brain patients and 4% for skull base patients. The maximum doses in organs at risk calculated with vRBE was up to 14 Gy above dose limit. The mean biological range extension was greater than 0.4 cm. DISCUSSION: Our method of estimation of biological range extension is insensitive for dose inhomogeneities and can be easily used for different proton plans with intensity-modulated proton therapy (IMPT) optimization. Using volumes instead of dose profiles, which is the common method, is more universal. However it was tested only for IMPT plans on fields arranged around the tumor area. CONCLUSIONS: Adopting a vRBE model results in an increase in dose and an extension of the beam range, which is especially disadvantageous in cancers close to organs at risk. Our results support the need to re-optimization of proton treatment plans when considering vRBE.

Positronium imaging with the novel multiphoton PET scanner.

In vivo assessment of cancer and precise location of altered tissues at initial stages of molecular disorders are important diagnostic challenges. Positronium is copiously formed in the free molecular spaces in the patient's body during positron emission tomography (PET). The positronium properties vary according to the size of inter- and intramolecular voids and the concentration of molecules in them such as, e.g., molecular oxygen, O2; therefore, positronium imaging may provide information about disease progression during the initial stages of molecular alterations. Current PET systems do not allow acquisition of positronium images. This study presents a new method that enables positronium imaging by simultaneous registration of annihilation photons and deexcitation photons from pharmaceuticals labeled with radionuclides. The first positronium imaging of a phantom built from cardiac myxoma and adipose tissue is demonstrated. It is anticipated that positronium imaging will substantially enhance the specificity of PET diagnostics.

3D Dosimetry Based on LiMgPO4 OSL Silicone Foils: Facilitating the Verification of Eye-Ball Cancer Proton Radiotherapy.

A direct verification of the three-dimensional (3D) proton clinical treatment plan prepared for tumor in the eyeball, using the Eclipse Ocular Proton Planning system (by Varian Medical Systems), has been presented. To achieve this, a prototype of the innovative two-dimensional (2D) circular silicone foils, made of a polymer with the embedded optically stimulated luminescence (OSL) material in powder form (LiMgPO4), and a self-developed optical imaging system, consisting of an illuminating light source and a high-sensitive CCD camera has been applied. A specially designed lifelike eyeball phantom has been used, constructed from 40 flat sheet LMP-based silicone foils stacked and placed together behind a spherical phantom made by polystyrene, all to reflect the curvature of the real eyeball. Two-dimensional OSL signals were captured and further analyzed from each single silicone foil after irradiation using a dedicated patient collimator and a 58.8 MeV modulated proton beam. The reconstructed 3D proton depth dose distribution matches very well with the clinical treatment plan, allowing for the consideration of the new OSL system for further 3D dosimetry applications within the proton radiotherapy area.

Study of relationship between dose, LET and the risk of brain necrosis after proton therapy for skull base tumors.

PURPOSE: We investigated the relationship between RBE-weighted dose (DRBE) calculated with constant (cRBE) and variable RBE (vRBE), dose-averaged linear energy transfer (LETd) and the risk of radiographic changes in skull base patients treated with protons. METHODS: Clinical treatment plans of 45 patients were recalculated with Monte Carlo tool FRED. Radiographic changes (i.e. edema and/or necrosis) were identified by MRI. Dosimetric parameters for cRBE and vRBE were computed. Biological margin extension and voxel-based analysis were employed looking for association of DRBE(vRBE) and LETd with brain edema and/or necrosis. RESULTS: When using vRBE, Dmax in the brain was above the highest dose limits for 38% of patients, while such limit was never exceeded assuming cRBE. Similar values of Dmax were observed in necrotic regions, brain and temporal lobes. Most of the brain necrosis was in proximity to the PTV. The voxel-based analysis did not show evidence of an association with high LETd values. CONCLUSIONS: When looking at standard dosimetric parameters, the higher dose associated with vRBE seems to be responsible for an enhanced risk of radiographic changes. However, as revealed by a voxel-based analysis, the large inter-patient variability hinders the identification of a clear effect for high LETd.