Quantification of deformable image registration uncertainties for dose accumulation on head and neck cancer proton treatments.

PURPOSE: Head and neck cancer (HNC) patients in radiotherapy require adaptive treatment plans due to anatomical changes. Deformable image registration (DIR) is used in adaptive radiotherapy, e.g. for deformable dose accumulation (DDA). However, DIR's ill-posedness necessitates addressing uncertainties, often overlooked in clinical implementations. DIR's further clinical implementation is hindered by missing quantitative commissioning and quality assurance tools. This study evaluates one pathway for more quantitative DDA uncertainties. METHODS: For five HNC patients, each with multiple repeated CTs acquired during treatment, a simultaneous-integrated boost (SIB) plan was optimized. Recalculated doses were warped individually using multiple DIRs from repeated to reference CTs, and voxel-by-voxel dose ranges determined an error-bar for DDA. Followed by evaluating, a previously proposed early-stage DDA uncertainty estimation method tested for lung cancer, which combines geometric DIR uncertainties, dose gradients and their directional dependence, in the context of HNC. RESULTS: Applying multiple DIRs show dose differences, pronounced in high dose gradient regions. The patient with largest anatomical changes (-13.1 % in ROI body volume), exhibited 33 % maximum uncertainty in contralateral parotid, with 54 % of voxels presenting an uncertainty >5 %. Accumulation over multiple CTs partially mitigated uncertainties. The estimation approach predicted 92.6 % of voxels within ±5 % to the reference dose uncertainty across all patients. CONCLUSIONS: DIR variations impact accumulated doses, emphasizing DDA uncertainty quantification's importance for HNC patients. Multiple DIR dose warping aids in quantifying DDA uncertainties. An estimation approach previously described for lung cancer was successfully validated for HNC, for SIB plans, presenting different dose gradients, and for accumulated treatments.

The impact of motion on onboard MRI-guided pencil beam scanned proton therapy treatments.

Objective.Online magnetic resonance imaging (MRI) guidance could be especially beneficial for pencil beam scanned (PBS) proton therapy of tumours affected by respiratory motion. For the first time to our knowledge, we investigate the dosimetric impact of respiratory motion on MRI-guided proton therapy compared to the scenario without magnetic field.Approach.A previously developed analytical proton dose calculation algorithm accounting for perpendicular magnetic fields was extended to enable 4D dose calculations. For two geometrical phantoms and three liver and two lung patient cases, static treatment plans were optimised with and without magnetic field (0, 0.5 and 1.5 T). Furthermore, plans were optimised using gantry angle corrections (0.5 T +5° and 1.5 T +15°) to reproduce similar beam trajectories compared to the 0 T reference plans. The effect of motion was then considered using 4D dose calculations without any motion mitigation and simulating 8-times volumetric rescanning, with motion for the patient cases provided by 4DCT(MRI) data sets. Each 4D dose calculation was performed for different starting phases and the CTV dose coverageV95%and homogeneityD5%-D95%were analysed.Main results.For the geometrical phantoms with rigid motion perpendicular to the beam and parallel to the magnetic field, a comparable dosimetric effect was observed independent of the magnetic field. Also for the five 4DCT(MRI) cases, the influence of motion was comparable for all magnetic field strengths with and without gantry angle correction. On average, the motion-induced decrease in CTVV95%from the static plan was 17.0% and 18.9% for 1.5 T and 0.5 T, respectively, and 19.9% without magnetic field.Significance.For the first time, this study investigates the combined impact of magnetic fields and respiratory motion on MR-guided proton therapy. The comparable dosimetric effects irrespective of magnetic field strength indicate that the effects of motion for future MR-guided proton therapy may not be worse than for conventional PBS proton therapy.

Automatic lung segmentation of magnetic resonance images: A new approach applied to healthy volunteers undergoing enhanced Deep-Inspiration-Breath-Hold for motion-mitigated 4D proton therapy of lung tumors.

Background and purpose: Respiratory suppression techniques represent an effective motion mitigation strategy for 4D-irradiation of lung tumors with protons. A magnetic resonance imaging (MRI)-based study applied and analyzed methods for this purpose, including enhanced Deep-Inspiration-Breath-Hold (eDIBH). Twenty-one healthy volunteers (41-58 years) underwent thoracic MR scans in four imaging sessions containing two eDIBH-guided MRIs per session to simulate motion-dependent irradiation conditions. The automated MRI segmentation algorithm presented here was critical in determining the lung volumes (LVs) achieved during eDIBH. Materials and methods: The study included 168 MRIs acquired under eDIBH conditions. The lung segmentation algorithm consisted of four analysis steps: (i) image preprocessing, (ii) MRI histogram analysis with thresholding, (iii) automatic segmentation, (iv) 3D-clustering. To validate the algorithm, 46 eDIBH-MRIs were manually contoured. Sørensen-Dice similarity coefficients (DSCs) and relative deviations of LVs were determined as similarity measures. Assessment of intrasessional and intersessional LV variations and their differences provided estimates of statistical and systematic errors. Results: Lung segmentation time for 100 2D-MRI planes was âˆ¼ 10 s. Compared to manual lung contouring, the median DSC was 0.94 with a lower 95 % confidence level (CL) of 0.92. The relative volume deviations yielded a median value of 0.059 and 95 % CLs of -0.013 and 0.13. Artifact-based volume errors, mainly of the trachea, were estimated. Estimated statistical and systematic errors ranged between 6 and 8 %. Conclusions: The presented analytical algorithm is fast, precise, and readily available. The results are comparable to time-consuming, manual segmentations and other automatic segmentation approaches. Post-processing to remove image artifacts is under development.

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.

Technical note: Towards more realistic 4DCT(MRI) numerical lung phantoms.

BACKGROUND: Numerical 4D phantoms, together with associated ground truth motion, offer a flexible and comprehensive data set for realistic simulations in radiotherapy and radiology in target sites affected by respiratory motion. PURPOSE: We present an openly available upgrade to previously reported methods for generating realistic 4DCT lung numerical phantoms, which now incorporate respiratory ribcage motion and improved lung density representation throughout the breathing cycle. METHODS: Density information of reference CTs, toget her with motion from multiple breathing cycle 4DMRIs have been combined to generate synthetic 4DCTs (4DCT(MRI)s). Inter-subject correspondence between the CT and MRI anatomy was first established via deformable image registration (DIR) of binary masks of the lungs and ribcage. Ribcage and lung motions were extracted independently from the 4DMRIs using DIR and applied to the corresponding locations in the CT after post-processing to preserve sliding organ motion. In addition, based on the Jacobian determinant of the resulting deformation vector fields, lung densities were scaled on a voxel-wise basis to more accurately represent changes in local lung density. For validating this process, synthetic 4DCTs, referred to as 4DCT(CT)s, were compared to the originating 4DCTs using motion extracted from the latter, and the dosimetric impact of the new features of ribcage motion and density correction were analyzed using pencil beam scanned proton 4D dose calculations. RESULTS: Lung density scaling led to a reduction of maximum mean lung Hounsfield units (HU) differences from 45 to 12 HU when comparing simulated 4DCT(CT)s to their originating 4DCTs. Comparing 4D dose distributions calculated on the enhanced 4DCT(CT)s to those on the original 4DCTs yielded 2%/2 mm gamma pass rates above 97% with an average improvement of 1.4% compared to previously reported phantoms. CONCLUSIONS: A previously reported 4DCT(MRI) workflow has been successfully improved and the resulting numerical phantoms exhibit more accurate lung density representations and realistic ribcage motion.

Treatment planning comparison for head and neck cancer between photon, proton, and combined proton-photon therapy - From a fixed beam line to an arc.

BACKGROUND AND PURPOSE: This study investigates whether combined proton-photon therapy (CPPT) improves treatment plan quality compared to single-modality intensity-modulated radiation therapy (IMRT) or intensity-modulated proton therapy (IMPT) for head and neck cancer (HNC) patients. Different proton beam arrangements for CPPT and IMPT are compared, which could be of specific interest concerning potential future upright-positioned treatments. Furthermore, it is evaluated if CPPT benefits remain under inter-fractional anatomical changes for HNC treatments. MATERIAL AND METHODS: Five HNC patients with a planning CT and multiple (4-7) repeated CTs were studied. CPPT with simultaneously optimized photon and proton fluence, single-modality IMPT, and IMRT treatment plans were optimized on the planning CT and then recalculated and reoptimized on each repeated CT. For CPPT and IMPT, plans with different degrees of freedom for the proton beams were optimized. Fixed horizontal proton beam line (FHB), gantry-like, and arc-like plans were compared. RESULTS: The target coverage for CPPT without adaptation is insufficient (average V95%=88.4 %), while adapted plans can recover the initial treatment plan quality for target (average V95%=95.5 %) and organs-at-risk. CPPT with increased proton beam flexibility increases plan quality and reduces normal tissue complication probability of Xerostomia and Dysphagia. On average, Xerostomia NTCP reductions compared to IMRT are -2.7 %/-3.4 %/-5.0 % for CPPT FHB/CPPT Gantry/CPPT Arc. The differences for IMPT FHB/IMPT Gantry/IMPT Arc are + 0.8 %/-0.9 %/-4.3 %. CONCLUSION: CPPT for HNC needs adaptive treatments. Increasing proton beam flexibility in CPPT, either by using a gantry or an upright-positioned patient, improves treatment plan quality. However, the photon component is substantially reduced, therefore, the balance between improved plan quality and costs must be further determined.

A fast analytical dose calculation approach for MRI-guided proton therapy.

Objective.Magnetic resonance (MR) is an innovative technology for online image guidance in conventional radiotherapy and is also starting to be considered for proton therapy as well. For MR-guided therapy, particularly for online plan adaptations, fast dose calculation is essential. Monte Carlo (MC) simulations, however, which are considered the gold standard for proton dose calculations, are very time-consuming. To address the need for an efficient dose calculation approach for MRI-guided proton therapy, we have developed a fast GPU-based modification of an analytical dose calculation algorithm incorporating beam deflections caused by magnetic fields.Approach.Proton beams (70-229 MeV) in orthogonal magnetic fields (0.5/1.5 T) were simulated using TOPAS-MC and central beam trajectories were extracted to generate look-up tables (LUTs) of incremental rotation angles as a function of water-equivalent depth. Beam trajectories are then reconstructed using these LUTs for the modified ray casting dose calculation. The algorithm was validated against MC in water, different materials and for four example patient cases, whereby it has also been fully incorporated into a treatment plan optimisation regime.Main results.Excellent agreement between analytical and MC dose distributions could be observed with sub-millimetre range deviations and differences in lateral shifts <2 mm even for high densities (1000 HU). 2%/2 mm gamma pass rates were comparable to the 0 T scenario and above 94.5% apart for the lung case. Further, comparable treatment plan quality could be achieved regardless of magnetic field strength.Significance.A new method for accurate and fast proton dose calculation in magnetic fields has been developed and successfully implemented for treatment plan optimisation.

Patient-specific quality assurance for deformable IMRT/IMPT dose accumulation: Proposition and validation of energy conservation based validation criterion.

BACKGROUND: Deformable image registration (DIR)-based dose accumulation (DDA) is regularly used in adaptive radiotherapy research. However, the applicability and reliability of DDA for direct clinical usage are still being debated. One primary concern is the validity of DDA, particularly for scenarios with substantial anatomical changes, for which energy-conservation problems were observed in conceptual studies. PURPOSE: We present and validate an energy-conservation (EC)-based DDA validation workflow and further investigate its usefulness for actual patient data, specifically for lung cancer cases. METHODS: For five non-small cell lung cancer (NSCLC) patients, DDA based on five selective DIR methods were calculated for five different treatment plans, which include one intensity-modulated photon therapy (IMRT), two intensity-modulated proton therapy (IMPT), and two combined proton-photon therapy (CPPT) plans. All plans were optimized on the planning CT (planCT) acquired in deep inspiration breath-hold (DIBH) and were re-optimized on the repeated DIBH CTs of three later fractions. The resulting fractional doses were warped back to the planCT using each DIR. An EC-based validation of the accumulation process was implemented and applied to all DDA results. Correlations between relative organ mass/volume variations and the extent of EC violation were then studied using Bayesian linear regression (BLR). RESULTS: For most OARs, EC violation within 10% is observed. However, for the PTVs and GTVs with substantial regression, severe overestimation of the fractional energy was found regardless of treatment type and applied DIR method. BLR results show that EC violation is linearly correlated to the relative mass variation (R^2 > 0.95) and volume variation (R^2 > 0.60). CONCLUSION: DDA results should be used with caution in regions with high mass/volume variation for intensity-based DIRs. EC-based validation is a useful approach to provide patient-specific quality assurance of the validity of DDA in radiotherapy.

Large anatomical changes in head-and-neck cancers - A dosimetric comparison of online and offline adaptive proton therapy.

Purpose: This work evaluates an online adaptive (OA) workflow for head-and-neck (H&N) intensity-modulated proton therapy (IMPT) and compares it with full offline replanning (FOR) in patients with large anatomical changes. Methods: IMPT treatment plans are created retrospectively for a cohort of eight H&N cancer patients that previously required replanning during the course of treatment due to large anatomical changes. Daily cone-beam CTs (CBCT) are acquired and corrected for scatter, resulting in 253 analyzed fractions. To simulate the FOR workflow, nominal plans are created on the planning-CT and delivered until a repeated-CT is acquired; at this point, a new plan is created on the repeated-CT. To simulate the OA workflow, nominal plans are created on the planning-CT and adapted at each fraction using a simple beamlet weight-tuning technique. Dose distributions are calculated on the CBCTs with Monte Carlo for both delivery methods. The total treatment dose is accumulated on the planning-CT. Results: Daily OA improved target coverage compared to FOR despite using smaller target margins. In the high-risk CTV, the median D98 degradation was 1.1 % and 2.1 % for OA and FOR, respectively. In the low-risk CTV, the same metrics yield 1.3 % and 5.2 % for OA and FOR, respectively. Smaller setup margins of OA reduced the dose to all OARs, which was most relevant for the parotid glands. Conclusion: Daily OA can maintain prescription doses and constraints over the course of fractionated treatment, even in cases of large anatomical changes, reducing the necessity for manual replanning in H&N IMPT.

Early outcome after craniospinal irradiation with pencil beam scanning proton therapy for children, adolescents and young adults with brain tumors.

Central nervous system (CNS) tumors are the most common solid malignancies in children and adolescents and young adults (C-AYAs). Craniospinal irradiation (CSI) is an essential treatment component for some malignancies, but it can also lead to important toxicity. Pencil beam scanning proton therapy (PBSPT) allows for a minimization of dose delivered to organs at risk and, thus, potentially reduced acute and late toxicity. This study aims to report the clinical outcomes and toxicity rates after CSI for C-AYAs treated with PBSPT. Seventy-one C-AYAs (median age: 7.4 years) with CNS tumors were treated with CSI between 2004 and 2021. Medulloblastoma (n = 42: 59%) and ependymoma (n = 8; 11%) were the most common histologies. Median prescribed total PBSPT dose was 54 GyRBE (range: 18-60.4), and median prescribed craniospinal dose was 24 GyRBE (range: 18-36.8). Acute and late toxicities were coded according to Common Terminology Criteria for Adverse Events. After a median follow-up of 24.5 months, the estimated 2-year local control, distant control, and overall survival were 86.3%, 80.5%, and 84.7%, respectively. Late grade ≥3 toxicity-free rate was 92.6% at 2 years. Recurrent and metastatic tumors were associated with worse outcome. In conclusion, excellent tumor control with low toxicity rates was observed in C-AYAs with brain tumors treated with CSI using PBSPT.

Optically stimulated luminescence dosimeters for simultaneous measurement of point dose and dose-weighted LET in an adaptive proton therapy workflow.

Purpose: To demonstrate the suitability of optically stimulated luminescence detectors (OSLDs) for accurate simultaneous measurement of the absolute point dose and dose-weighted linear energy transfer (LETD) in an anthropomorphic phantom for experimental validation of daily adaptive proton therapy. Methods: A clinically realistic intensity-modulated proton therapy (IMPT) treatment plan was created based on a CT of an anthropomorphic head-and-neck phantom made of tissue-equivalent material. The IMPT plan was optimized with three fields to deliver a uniform dose to the target volume covering the OSLDs. Different scenarios representing inter-fractional anatomical changes were created by modifying the phantom. An online adaptive proton therapy workflow was used to recover the daily dose distribution and account for the applied geometry changes. To validate the adaptive workflow, measurements were performed by irradiating Al2O3:C OSLDs inside the phantom. In addition to the measurements, retrospective Monte Carlo simulations were performed to compare the absolute dose and dose-averaged LET (LETD) delivered to the OSLDs. Results: The online adaptive proton therapy workflow was shown to recover significant degradation in dose conformity resulting from large anatomical and positioning deviations from the reference plan. The Monte Carlo simulations were in close agreement with the OSLD measurements, with an average relative error of 1.4% for doses and 3.2% for LETD. The use of OSLDs for LET determination allowed for a correction for the ionization quenched response. Conclusion: The OSLDs appear to be an excellent detector for simultaneously assessing dose and LET distributions in proton irradiation of an anthropomorphic phantom. The OSLDs can be cut to almost any size and shape, making them ideal for in-phantom measurements to probe the radiation quality and dose in a predefined region of interest. Although we have presented the results obtained in the experimental validation of an adaptive proton therapy workflow, the same approach can be generalized and used for a variety of clinical innovations and workflow developments that require accurate assessment of point dose and/or average LET.

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.

Low-Dose Computed Tomography Scanning Protocols for Online Adaptive Proton Therapy of Head-and-Neck Cancers.

PURPOSE: To evaluate the suitability of low-dose CT protocols for online plan adaptation of head-and-neck patients. METHODS: We acquired CT scans of a head phantom with protocols corresponding to CT dose index volume CTDIvol in the range of 4.2-165.9 mGy. The highest value corresponds to the standard protocol used for CT simulations of 10 head-and-neck patients included in the study. The minimum value corresponds to the lowest achievable tube current of the GE Discovery RT scanner used for the study. For each patient and each low-dose protocol, the noise relative to the standard protocol, derived from phantom images, was applied to a virtual CT (vCT). The vCT was obtained from a daily CBCT scan corresponding to the fraction with the largest anatomical changes. We ran an established adaptive workflow twice for each low-dose protocol using a high-quality daily vCT and the corresponding low-dose synthetic vCT. For a relative comparison of the adaptation efficacy, two adapted plans were recalculated in the high-quality vCT and evaluated with the contours obtained through deformable registration of the planning CT. We also evaluated the accuracy of dose calculation in low-dose CT volumes using the standard CT protocol as reference. RESULTS: The maximum differences in D98 between low-dose protocols and the standard protocol for the high-risk and low-risk CTV were found to be 0.6% and 0.3%, respectively. The difference in OAR sparing was up to 3%. The Dice similarity coefficient between propagated contours obtained with low-dose and standard protocols was above 0.982. The mean 2%/2 mm gamma pass rate for the lowest-dose image, using the standard protocol as reference, was found to be 99.99%. CONCLUSION: The differences between low-dose protocols and the standard scanning protocol were marginal. Thus, low-dose CT protocols are suitable for online adaptive proton therapy of head-and-neck cancers. As such, considering scanning protocols used in our clinic, the imaging dose associated with online adaption of head-and-neck cancers treated with protons can be reduced by a factor of 40.

Impact of spot reduction on the effectiveness of rescanning in pencil beam scanned proton therapy for mobile tumours.

Objective. In pencil beam scanning proton therapy, individually calculated and positioned proton pencil beams, also referred to as 'spots', are used to achieve a highly conformal dose distributions to the target. Recent work has shown that this number of spots can be substantially reduced, resulting in shorter delivery times without compromising dosimetric plan quality. However, the sensitivity of spot-reduced plans to tumour motion is unclear. Although previous work has shown that spot-reduced plans are slightly more sensitive to small positioning inaccuracies of the individual pencil beams, the resulting shorter delivery times may allow for more rescanning. The aim of this study was to assess the impact of tumour motion and the effectiveness of 3D volumetric rescanning for spot-reduced treatment plans.Approach.Three liver and two lung cancer patients with non-negligible motion amplitudes were analysed. Conventional and probabilistic internal target volume definitions were used for planning considering single or multiple breathing cycles respectively. For each patient, one clinical and two spot-reduced treatment plans were created using identical field geometries. 4D dynamic dose calculations were then performed and resulting target coverage (V95%), dose homogeneity (D5%-D95%) and hot spots (D2%) evaluated for 1-25 rescans.Main results. Over all patients investigated, spot reduction reduced the number of spots by 91% in comparison to the clinical plan, reducing field delivery times by approximately 50%. This reduction, together with the substantially increased dose per spot resulting from the spot reduction process, allowed for more rescans in the same amount of time as for clinical plans and typically improved dosimetric parameters, in some cases to values better than the reference static (3D calculated) plans. However, spot-reduced plans had an increased possibility of interference with the breathing cycle, especially for simulations of perfectly repeatable breathing.Significance.For the patients analysed in this study, spot-reduced plans were found to be a valuable option to increase the efficiency of 3D volumetric rescanning for motion mitigation, if attention is paid to possible interference patterns.

Ultra-fast pencil beam scanning proton therapy for locally advanced non-small-cell lung cancers: Field delivery within a single breath-hold.

PURPOSE: The use of motion mitigation techniques such as breath-hold can reduce the dosimetric uncertainty of lung cancer proton therapy. We studied the feasibility of pencil beam scanning (PBS) proton therapy field delivery within a single breath-hold at PSI's Gantry 2. METHODS: In PBS proton therapy, the delivery time for a field is determined by the beam-on time and the dead time between proton spots (the time required to change the energy and/or lateral position). We studied ways to reduce beam-on and lateral scanning time, without sacrificing dosimetric plan quality, aiming at a single field delivery time of 15 seconds at maximum. We tested this approach on 10 lung cases with varying target volumes. To reduce the beam-on time, we increased the beam current at the isocenter by developing new beam optics for PSI's PROSCAN beamline and Gantry 2. To reduce the dead time between the spots, we used spot-reduced plan optimization. RESULTS: We found that it is possible to achieve conventional fractionated (2 Gy(RBE)/fraction) and hypofractionated (6 Gy(RBE)/fraction) field delivery times within a single breath-hold (<15 sec) for a variety non-small-cell lung cancer cases. CONCLUSION: In summary, the combination of spot reduction and improved beam line transmission is a promising approach for the treatment of mobile tumours within clinically achievable breath-hold durations.

Synthetic 4DCT(MRI) lung phantom generation for 4D radiotherapy and image guidance investigations.

PURPOSE: Respiratory motion is one of the major challenges in radiotherapy. In this work, a comprehensive and clinically plausible set of 4D numerical phantoms, together with their corresponding "ground truths," have been developed and validated for 4D radiotherapy applications. METHODS: The phantoms are based on CTs providing density information and motion from multi-breathing-cycle 4D Magnetic Resonance imagings (MRIs). Deformable image registration (DIR) has been utilized to extract motion fields from 4DMRIs and to establish inter-subject correspondence by registering binary lung masks between Computer Tomography (CT) and MRI. The established correspondence is then used to warp the CT according to the 4DMRI motion. The resulting synthetic 4DCTs are called 4DCT(MRI)s. Validation of the 4DCT(MRI) workflow was conducted by directly comparing conventional 4DCTs to derived synthetic 4D images using the motion of the 4DCTs themselves (referred to as 4DCT(CT)s). Digitally reconstructed radiographs (DRRs) as well as 4D pencil beam scanned (PBS) proton dose calculations were used for validation. RESULTS: Based on the CT image appearance of 13 lung cancer patients and deformable motion of five volunteer 4DMRIs, synthetic 4DCT(MRI)s with a total of 871 different breathing cycles have been generated. The 4DCT(MRI)s exhibit an average superior-inferior tumor motion amplitude of 7 ± 5 mm (min: 0.5 mm, max: 22.7 mm). The relative change of the DRR image intensities of the conventional 4DCTs and the corresponding synthetic 4DCT(CT)s inside the body is smaller than 5% for at least 81% of the pixels for all studied cases. Comparison of 4D dose distributions calculated on 4DCTs and the synthetic 4DCT(CT)s using the same motion achieved similar dose distributions with an average 2%/2 mm gamma pass rate of 90.8% (min: 77.8%, max: 97.2%). CONCLUSION: We developed a series of numerical 4D lung phantoms based on real imaging and motion data, which give realistic representations of both anatomy and motion scenarios and the accessible "ground truth" deformation vector fields of each 4DCT(MRI). The open-source code and motion data allow foreseen users to generate further 4D data by themselves. These numeric 4D phantoms can be used for the development of new 4D treatment strategies, 4D dose calculations, DIR algorithm validations, as well as simulations of motion mitigation and different online image guidance techniques for both proton and photon radiation therapy.

NTCP Modeling for High-Grade Temporal Radionecroses in a Large Cohort of Patients Receiving Pencil Beam Scanning Proton Therapy for Skull Base and Head and Neck Tumors.

PURPOSE: To develop a normal tissue complication probability model including clinical and dosimetric parameters for high-grade temporal lobe radionecroses (TRN) after pencil beam scanning proton therapy. METHODS AND MATERIALS: We included data on 299 patients with skull base and head and neck tumors treated with pencil-beam scan proton therapy, with a total dose of ≥60 GyRBE (relative biological effectiveness) from May 2004 to November 2018. We considered 9 clinical and 27 dosimetric parameters for the structure-wise modeling of high-grade (grade ≥2) TRN. After eliminating strongly cross-correlated variables, we generated logistic regression models using least absolute shrinkage and selection operator regression. We performed bootstrapping to assess parameter selection robustness and evaluated model performance via cross-correlation by assessing the area under the curve of receiver operating characteristic curves and calibration with a Hosmer-Lemeshow test statistic. RESULTS: After a median radiologic follow-up of 51.5 months (range, 4-190), 27 patients (9%) developed grade ≥2 TRN. Eleven patients had bitemporal necrosis, resulting in 38 events in 598 temporal lobes for structure-wise analysis. During our bootstrapping analysis, we found that the highest selection frequency was for prescription dose, followed by age, V40Gy (%), hypertension, and dose to at least 1 cc (D1cc) (Gy) in the temporal lobe. During our cross-validation, we found that age*prescription-dose*D1cc (Gy)*hypertension was superior in all described test statistics. We built full cohort structure-wise and patient-wise models with maximum area under the curve of receiver operating characteristic curves of 0.79 (structure-wise) and 0.76 (patient-wise). CONCLUSIONS: While developing a logistic regression normal tissue complication probability model to predict grade ≥2 TRN, the best fit was found for the model containing age, prescription dose, D1cc (Gy), and hypertensive blood pressure as risk factors. External validation will be the next step to improve generalizability and potential introduction into clinical routine.

Clinical necessity of multi-image based (4DMIB) optimization for targets affected by respiratory motion and treated with scanned particle therapy - A comprehensive review.

4D multi-image-based (4DMIB) optimization is a form of robust optimization where different uncertainty scenarios, due to anatomy variations, are considered via multiple image sets (e.g., 4DCT). In this review, we focused on providing an overview of different 4DMIB optimization implementations, introduced various frameworks to evaluate the robustness of scanned particle therapy affected by breathing motion and summarized the existing evidence on the necessity of using 4DMIB optimization clinically. Expected potential benefits of 4DMIB optimization include more robust and/or interplay-effect-resistant doses for the target volume and organs-at-risk for indications affected by anatomical variations (e.g., breathing, peristalsis, etc.). Although considerable literature is available on the research and technical aspects of 4DMIB, clinical studies are rare and often contain methodological limitations, such as, limited patient number, motion amplitude, motion and delivery time structure considerations, number of repeat CTs, etc. Therefore, the data are not conclusive. In addition, multiple studies have found that robust 3D optimized plans result in dose distributions within the set clinical tolerances and, therefore, are suitable for a treatment of moving targets with scanned particle therapy. We, therefore, consider the clinical necessity of 4DMIB optimization, when treating moving targets with scanned particle therapy, as still to be demonstrated.

A quantitative FLASH effectiveness model to reveal potentials and pitfalls of high dose rate proton therapy.

PURPOSE: In ultrahigh dose rate radiotherapy, the FLASH effect can lead to substantially reduced healthy tissue damage without affecting tumor control. Although many studies show promising results, the underlying biological mechanisms and the relevant delivery parameters are still largely unknown. It is unclear, particularly for scanned proton therapy, how treatment plans could be optimized to maximally exploit this protective FLASH effect. MATERIALS AND METHODS: To investigate the potential of pencil beam scanned proton therapy for FLASH treatments, we present a phenomenological model, which is purely based on experimentally observed phenomena such as potential dose rate and dose thresholds, and which estimates the biologically effective dose during FLASH radiotherapy based on several parameters. We applied this model to a wide variety of patient geometries and proton treatment planning scenarios, including transmission and Bragg peak plans as well as single- and multifield plans. Moreover, we performed a sensitivity analysis to estimate the importance of each model parameter. RESULTS: Our results showed an increased plan-specific FLASH effect for transmission compared with Bragg peak plans (19.7% vs. 4.0%) and for single-field compared with multifield plans (14.7% vs. 3.7%), typically at the cost of increased integral dose compared to the clinical reference plan. Similar FLASH magnitudes were found across the different treatment sites, whereas the clinical benefits with respect to the clinical reference plan varied strongly. The sensitivity analysis revealed that the threshold dose as well as the dose per fraction strongly impacted the FLASH effect, whereas the persistence time only marginally affected FLASH. An intermediate dependence of the FLASH effect on the dose rate threshold was found. CONCLUSIONS: Our model provided a quantitative measure of the FLASH effect for various delivery and patient scenarios, supporting previous assumptions about potentially promising planning approaches for FLASH proton therapy. Positive clinical benefits compared to clinical plans were achieved using hypofractionated, single-field transmission plans. The dose threshold was found to be an important factor, which may require more investigation.

Commissioning and quality assurance of a novel solution for respiratory-gated PBS proton therapy based on optical tracking of surface markers.

We present the commissioning and quality assurance of our clinical protocol for respiratory gating in pencil beam scanning proton therapy for cancer patients with moving targets. In a novel approach, optical tracking has been integrated in the therapy workflow and used to monitor respiratory motion from multiple surrogates, applied on the patients' chest. The gating system was tested under a variety of experimental conditions, specific to proton therapy, to evaluate reaction time and reproducibility of dose delivery control. The system proved to be precise in the application of beam gating and allowed the mitigation of dose distortions even for large (1.4cm) motion amplitudes, provided that adequate treatment windows were selected. The total delivered dose was not affected by the use of gating, with measured integral error within 0.15cGy. Analysing high-resolution images of proton transmission, we observed negligible discrepancies in the geometric location of the dose as a function of the treatment window, with gamma pass rate greater than 95% (2%/2mm) compared to stationary conditions. Similarly, pass rate for the latter metric at the 3%/3mm level was observed above 97% for clinical treatment fields, limiting residual movement to 3mm at end-exhale. These results were confirmed in realistic clinical conditions using an anthropomorphic breathing phantom, reporting a similarly high 3%/3mm pass rate, above 98% and 94%, for regular and irregular breathing, respectively. Finally, early results from periodic QA tests of the optical tracker have shown a reliable system, with small variance observed in static and dynamic measurements.