A Monte Carlo based radiation response modelling framework to assess variability of clinical RBE in proton therapy

The clinical implementation of a variable relative biological effectiveness (RBE) in proton therapy is currently controversially discussed. Initial clinical evidence indicates a variable proton RBE, which needs to be verified. In this study, a radiation response modelling framework for assessing clinical RBE variability is established. It was applied to four selected glioma patients (grade III) treated with adjuvant radio(chemo)therapy and who developed late morphological image changes on T1-weighted contrast-enhanced (T1w-CE) magnetic resonance (MR) images within approximately two years of recurrence-free follow-up. The image changes were correlated voxelwise with dose and linear energy transfer (LET) values using univariable and multivariable logistic regression analysis. The regression models were evaluated by the area-under-the-curve (AUC) method performing a leave-one-out cross validation. The tolerance dose TD₅₀ at which 50% of patient voxels experienced toxicity was interpolated from the models. A Monte Carlo (MC) model was developed to simulate dose and LET distributions, which includes variance reduction (VR) techniques to decrease computation time. Its reliability and accuracy were evaluated based on dose calculations of the clinical treatment planning system (TPS) as well as absolute dose measurements performed in the patient specific quality assurance. Morphological image changes were related to a combination of dose and LET. The multivariable models revealed cross-validated AUC values of up to 0.88. The interpolated TD₅₀ curves decreased with increasing LET indicating an increase in biological effectiveness. The MC model reliably predicted average TPS dose within the clinical target volume as well as absolute water phantom dose measurements within 2% accuracy using dedicated VR settings. The observed correlation of dose and LET with late brain tissue damage suggests considering RBE variability for predicting chronic radiation-induced brain toxicities. The MC model simulates radiation fields in patients precisely and time-efficiently. Hence, this study encourages and enables in-depth patient evaluation to assess the variability of clinical proton RBE.

A Bayesian Approach for Measurements of Stray Neutrons at Proton Therapy Facilities: Quantifying Neutron Dose Uncertainty

Bonner sphere measurements are typically analyzed using unfolding codes. It is well known that it is difficult to get reliable estimates of uncertainties for standard unfolding procedures. An alternative approach is to analyze the data using Bayesian parameter estimation. This method provides reliable estimates of the uncertainties of neutron spectra leading to rigorous estimates of uncertainties of the dose. We extend previous Bayesian approaches and apply the method to stray neutrons in proton therapy environments by introducing a new parameterized model which describes the main features of the expected neutron spectra. The parameterization is based on information that is available from measurements and detailed Monte Carlo simulations. The validity of this approach has been validated with results of an experiment using Bonner spheres carried out at the experimental hall of the OncoRay proton therapy facility in Dresden.

Towards clinical application: prompt gamma imaging of passively scattered proton fields with a knife-edge slit camera

Prompt γ-ray imaging with a knife-edge shaped slit camera provides the possibility of verifying proton beam range in tumor therapy. Dedicated experiments regarding the characterization of the camera system have been performed previously. Now, we aim at implementing the prototype into clinical application of monitoring patient treatments. Focused on this goal of translation into clinical operation, we systematically addressed remaining challenges and questions. We developed a robust energy calibration routine and corresponding quality assurance protocols. Furthermore, with dedicated experiments, we determined the positioning precision of the system to 1.1 mm (2σ). For the first time, we demonstrated the application of the slit camera, which was intentionally developed for pencil beam scanning, to double scattered proton beams. Systematic experiments with increasing complexity were performed. It was possible to visualize proton range shifts of 2-5 mm with the camera system in phantom experiments in passive scattered fields. Moreover, prompt γ-ray profiles for single iso-energy layers were acquired by synchronizing time resolved measurements to the rotation of the range modulator wheel of the treatment system. Thus, a mapping of the acquired profiles to different anatomical regions along the beam path is feasible and additional information on the source of potential range shifts can be obtained. With the work presented here, we show that an application of the slit camera in clinical treatments is possible and of potential benefit.

Experimental investigation of irregular motion impact on 4D PET-based particle therapy monitoring

Particle therapy positron emission tomography (PT-PET) is an in vivo and non-invasive imaging technique to monitor treatment delivery in particle therapy. The inevitable patient respiratory motion during irradiation causes artefacts and inaccurate activity distribution in PET images. Four-dimensional (4D) maximum likelihood expectation maximisation (4D MLEM) allows for a compensation of these effects, but has up to now been restricted to regular motion for PT-PET investigations. However, intra-fractional motion during treatment might differ from that during acquisition of the 4D-planning CT (e.g. amplitude variation, baseline drift) and therefore might induce inaccurate 4D PET reconstruction results. This study investigates the impact of different irregular analytical one-dimensional (1D) motion patterns on PT-PET imaging by means of experiments with a radioactive source and irradiated moving phantoms. Three sorting methods, namely phase sorting, equal amplitude sorting and event-based amplitude sorting, were applied to manage the PET list-mode data. The influence of these sorting methods on the motion compensating algorithm has been analysed. The event-based amplitude sorting showed a superior performance and it is applicable for irregular motions with ⩽ 4 mm amplitude elongation and drift. For motion with 10 mm baseline drift, the normalised root mean square error was as high as 10.5% and a 10 mm range deviation was observed.

Efficient laser-ion acceleration from closely stacked ultrathin foils

A new scheme to efficiently accelerate protons by a single linear polarized high-intensity ultrashort laser pulse using multiple ultrathin foils is proposed. The foils are stacked at a spacing comparable to their thickness and subsequently irradiated by the same laser pulse. The foil thicknesses are chosen such that the laser light pressure can displace all electrons out of the foil. The authors present a simple, yet precise dynamical model of the acceleration process from which both optimum foil thickness and spacing can be derived. Extensive two-dimensional (2D) particle-in-cell simulations verify the model predictions and suggest an enhancement of the maximum proton kinetic energy by 30% for the two-foil case compared to a single foil.