Recent developments and comprehensive evaluations of a GPU-based Monte Carlo package for proton therapy

Evaluating the effect of setup uncertainty reduction and adaptation to geometric changes on normal tissue complication probability using online adaptive head and neck intensity modulated proton therapy

Objective. To evaluate the impact of setup uncertainty reduction (SUR) and adaptation to geometrical changes (AGC) on normal tissue complication probability (NTCP) when using online adaptive head and neck intensity modulated proton therapy (IMPT).Approach.A cohort of ten retrospective head and neck cancer patients with daily scatter corrected cone-beam CT (CBCT) was studied. For each patient, two IMPT treatment plans were created: one with a 3 mm setup uncertainty robustness setting and one with no explicit setup robustness. Both plans were recalculated on the daily CBCT considering three scenarios: the robust plan without adaptation, the non-robust plan without adaptation and the non-robust plan with daily online adaptation. Online-adaptation was simulated using an in-house developed workflow based on GPU-accelerated Monte Carlo dose calculation and partial spot-intensity re-optimization. Dose distributions associated with each scenario were accumulated on the planning CT, where NTCP models for six toxicities were applied. NTCP values from each scenario were intercompared to quantify the reduction in toxicity risk induced by SUR alone, AGC alone and SUR and AGC combined. Finally, a decision tree was implemented to assess the clinical significance of the toxicity reduction associated with each mechanism.Main results. For most patients, clinically meaningful NTCP reductions were only achieved when SUR and AGC were performed together. In these conditions, total reductions in NTCP of up to 30.48 pp were obtained, with noticeable NTCP reductions for aspiration, dysphagia and xerostomia (mean reductions of 8.25, 5.42 and 5.12 pp respectively). While SUR had a generally larger impact than AGC on NTCP reductions, SUR alone did not induce clinically meaningful toxicity reductions in any patient, compared to only one for AGC alone.SignificanceOnline adaptive head and neck proton therapy can only yield clinically significant reductions in the risk of long-term side effects when combining the benefits of SUR and AGC.

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 CTDI_vol 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 D₉₈ 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.

MOQUI: an open-source GPU-based Monte Carlo code for proton dose calculation with efficient data structure

Objective.Monte Carlo (MC) codes are increasingly used for accurate radiotherapy dose calculation. In proton therapy, the accuracy of the dose calculation algorithm is expected to have a more significant impact than in photon therapy due to the depth-dose characteristics of proton beams. However, MC simulations come at a considerable computational cost to achieve statistically sufficient accuracy. There have been efforts to improve computational efficiency while maintaining sufficient accuracy. Among those, parallelizing particle transportation using graphic processing units (GPU) achieved significant improvements. Contrary to the central processing unit, a GPU has limited memory capacity and is not expandable. It is therefore challenging to score quantities with large dimensions requiring extensive memory. The objective of this study is to develop an open-source GPU-based MC package capable of scoring those quantities.Approach.We employed a hash-table, one of the key-value pair data structures, to efficiently utilize the limited memory of the GPU and score the quantities requiring a large amount of memory. With the hash table, only voxels interacting with particles will occupy memory, and we can search the data efficiently to determine their address. The hash-table was integrated with a novel GPU-based MC code, moqui.Main results.The developed code was validated against an MC code widely used in proton therapy, TOPAS, with homogeneous and heterogeneous phantoms. We also compared the dose calculation results of clinical treatment plans. The developed code agreed with TOPAS within 2%, except for the fall-off and regions, and the gamma pass rates of the results were >99% for all cases with a 2 mm/2% criteria.Significance.We can score dose-influence matrix and dose-rate on a GPU for a 3-field H&N case with 10 GB of memory using moqui, which would require more than 100 GB of memory with the conventionally used array data structure.

Millisecond speed deep learning based proton dose calculation with Monte Carlo accuracy

Objective.Next generation online and real-time adaptive radiotherapy workflows require precise particle transport simulations in sub-second times, which is unfeasible with current analytical pencil beam algorithms (PBA) or Monte Carlo (MC) methods. We present a deep learning based millisecond speed dose calculation algorithm (DoTA) accurately predicting the dose deposited by mono-energetic proton pencil beams for arbitrary energies and patient geometries.Approach.Given the forward-scattering nature of protons, we frame 3D particle transport as modeling a sequence of 2D geometries in the beam's eye view. DoTA combines convolutional neural networks extracting spatial features (e.g. tissue and density contrasts) with a transformer self-attention backbone that routes information between the sequence of geometry slices and a vector representing the beam's energy, and is trained to predict low noise MC simulations of proton beamlets using 80 000 different head and neck, lung, and prostate geometries.Main results.Predicting beamlet doses in 5 ± 4.9 ms with a very high gamma pass rate of 99.37 ± 1.17% (1%, 3 mm) compared to the ground truth MC calculations, DoTA significantly improves upon analytical pencil beam algorithms both in precision and speed. Offering MC accuracy 100 times faster than PBAs for pencil beams, our model calculates full treatment plan doses in 10-15 s depending on the number of beamlets (800-2200 in our plans), achieving a 99.70 ± 0.14% (2%, 2 mm) gamma pass rate across 9 test patients.Significance.Outperforming all previous analytical pencil beam and deep learning based approaches, DoTA represents a new state of the art in data-driven dose calculation and can directly compete with the speed of even commercial GPU MC approaches. Providing the sub-second speed required for adaptive treatments, straightforward implementations could offer similar benefits to other steps of the radiotherapy workflow or other modalities such as helium or carbon treatments.

CT-on-Rails Versus In-Room CBCT for Online Daily Adaptive Proton Therapy of Head-and-Neck Cancers

PURPOSE

To compare the efficacy of CT-on-rails versus in-room CBCT for daily adaptive proton therapy.

METHODS

We analyzed a cohort of ten head-and-neck patients with daily CBCT and corresponding virtual CT images. The necessity of moving the patient after a CT scan is the most significant difference in the adaptation workflow, leading to an increased treatment execution uncertainty σ. It is a combination of the isocenter-matching σ_i and random patient movements induced by the couch motion σ_m. The former is assumed to never exceed 1 mm. For the latter, we studied three different scenarios with σ_m = 1, 2, and 3 mm. Accordingly, to mimic the adaptation workflow with CT-on-rails, we introduced random offsets after Monte-Carlo-based adaptation but before delivery of the adapted plan.

RESULTS

There were no significant differences in accumulated dose-volume histograms and dose distributions for σ_m = 1 and 2 mm. Offsets with σ_m = 3 mm resulted in underdosage to CTV and hot spots of considerable volume.

CONCLUSION

Since σ_m typically does not exceed 2 mm for in-room CT, there is no clinically significant dosimetric difference between the two modalities for online adaptive therapy of head-and-neck patients. Therefore, in-room CT-on-rails can be considered a good alternative to CBCT for adaptive proton therapy.

Adaptive proton therapy

Radiation therapy treatments are typically planned based on a single image set, assuming that the patient's anatomy and its position relative to the delivery system remains constant during the course of treatment. Similarly, the prescription dose assumes constant biological dose-response over the treatment course. However, variations can and do occur on multiple time scales. For treatment sites with significant intra-fractional motion, geometric changes happen over seconds or minutes, while biological considerations change over days or weeks. At an intermediate timescale, geometric changes occur between daily treatment fractions. Adaptive radiation therapy is applied to consider changes in patient anatomy during the course of fractionated treatment delivery. While traditionally adaptation has been done off-line with replanning based on new CT images, online treatment adaptation based on on-board imaging has gained momentum in recent years due to advanced imaging techniques combined with treatment delivery systems. Adaptation is particularly important in proton therapy where small changes in patient anatomy can lead to significant dose perturbations due to the dose conformality and finite range of proton beams. This review summarizes the current state-of-the-art of on-line adaptive proton therapy and identifies areas requiring further research.

Anatomic changes in head and neck intensity-modulated proton therapy: Comparison between robust optimization and online adaptation

BACKGROUND/PURPOSE

Setup variations and anatomical changes can severely affect the quality of head and neck intensity-modulated proton therapy (IMPT) treatments. The impact of these changes can be alleviated by increasing the plan's robustness a priori, or by adapting the plan online. This work compares these approaches in the context of head and neck IMPT.

MATERIALS/METHODS

A representative cohort of 10 head and neck squamous cell carcinoma (HNSCC) patients with daily cone-beam computed tomography (CBCT) was evaluated. For each patient, three IMPT plans were created: 1- a classical robust optimization (cRO) plan optimized on the planning CT, 2- an anatomical robust optimization (aRO) plan additionally including the two first daily CBCTs and 3- a plan optimized without robustness constraints, but online-adapted (OA) daily, using a constrained spot intensity re-optimization technique only.

RESULTS

The cumulative dose following OA fulfilled the clinical objective of both the high-risk and low-risk clinical target volumes (CTV) coverage in all 10 patients, compared to 8 for aRO and 4 for cRO. aRO did not significantly increase the dose to most organs at risk compared to cRO, although the integral dose was higher. OA significantly reduced the integral dose to healthy tissues compared to both robust methods, while providing equivalent or superior target coverage.

CONCLUSION

Using a simple spot intensity re-optimization, daily OA can achieve superior target coverage and lower dose to organs at risk than robust optimization methods.

Improving Proton Dose Calculation Accuracy by Using Deep Learning

Introduction

Pencil beam (PB) dose calculation is fast but inaccurate due to the approximations when dealing with inhomogeneities. Monte Carlo (MC) dose calculation is the most accurate method but it is time consuming. The aim of this study was to develop a deep learning model that can boost the accuracy of PB dose calculation to the level of MC dose by converting PB dose to MC dose for different tumor sites.

Methods

The proposed model uses the PB dose and CT image as inputs to generate the MC dose. We used 290 patients (90 head and neck, 93 liver, 75 prostate and 32 lung) to train, validate, and test the model. For each tumor site, we performed four numerical experiments to explore various combinations of training datasets.

Results

Training the model on data from all tumor sites together and using the dose distribution of each individual beam as input yielded the best performance for all four tumor sites. The average gamma passing rate (1mm/1%) between the converted and the MC dose was 92.8%, 92.7%, 89.7% and 99.6% for head and neck, liver, lung, and prostate test patients, respectively. The average dose conversion time for a single field was less than 4 seconds. The trained model can be adapted to new datasets through transfer learning.

Conclusions

Our deep learning-based approach can quickly boost the accuracy of PB dose to that of MC dose. The developed model can be added to the clinical workflow of proton treatment planning to improve dose calculation accuracy.

Comparison of weekly and daily online adaptation for head and neck intensity-modulated proton therapy

The high conformality of intensity-modulated proton therapy (IMPT) dose distributions causes treatment plans to be sensitive to geometrical changes during the course of a fractionated treatment. This can be addressed using adaptive proton therapy (APT). One important question in APT is the frequency of adaptations performed during a fractionated treatment, which is related to the question whether plan adaptation has to be done online or offline. The purpose of this work is to investigate the impact of weekly and daily online IMPT plan adaptation on the treatment quality for head and neck patients. A cohort of ten head and neck patients with daily acquired cone-beam CT (CBCT) images was evaluated retrospectively. Dose tracking of the IMPT treatment was performed for three scenarios: base plan with no adaptation (BP), weekly online adaptation (OA_W), and daily online adaptation (OA_D). Both adaptation schemes used an in-house developed online APT workflow, performing Monte Carlo (MC) dose calculations on scatter-corrected CBCTs. IMPT plan adaptation was achieved by only tuning the weights of a subset of beamlets, based on deformable image registration from the planning CT to each CBCT. Although OA_Dmitigated random delivery errors more effectively than OA_Won a fraction per fraction basis, both OA_Wand OA_Dachieved the clinical goals for all ten patients, while BP failed for six cases. In the high-risk CTV, accumulated values of D_98%ranged between 97.15% and 99.73% of the prescription dose for OA_D, with a median of 98.07%. For OA_W, values between 95.02% and 99.26% were obtained, with a median of 97.61% of the prescription dose. Otherwise, the dose to most organs at risk was similar for all three scenarios. Globally, our results suggest that OA_Wcould be used as an alternative approach to OA_Dfor most patients in order to reduce the clinical workload.

Evaluation of an a priori scatter correction algorithm for cone-beam computed tomography based range and dose calculations in proton therapy

BACKGROUND AND PURPOSE

Scatter correction of cone-beam computed tomography (CBCT) projections may enable accurate online dose-delivery estimations in photon and proton-based radiotherapy. This study aimed to evaluate the impact of scatter correction in CBCT-based proton range/dose calculations, in scans acquired in both proton and photon gantries.

MATERIAL AND METHODS

CBCT projections of a Catphan and an Alderson phantom were acquired on both a proton and a photon gantry. The scatter corrected CBCTs (corrCBCTs) and the clinical reconstructions (stdCBCTs) were compared against CTs rigidly registered to the CBCTs (rigidCTs). The CBCTs of the Catphan phantom were segmented by materials for CT number analysis. Water equivalent path length (WEPL) maps were calculated through the Alderson phantom while proton plans optimized on the rigidCT and recalculated on all CBCTs were compared in a gamma analysis.

RESULTS

In medium and high-density materials, the corrCBCT CT numbers were much closer to those of the rigidCT than the stdCBCTs. E.g. in the 50% bone segmentations the differences were reduced from above 300 HU (with stdCBCT) to around 60-70 HU (with corrCBCT). Differences in WEPL from the rigidCT were typically well below 5 mm for the corrCBCTs, compared to well above 10 mm for the stdCBCTs with the largest deviations in the head and thorax regions. Gamma pass rates (2%/2mm) when comparing CBCT-based dose re-calculations to rigidCT calculations were improved from around 80% (with stdCBCT) to mostly above 90% (with corrCBCT).

CONCLUSION

Scatter correction leads to substantial artefact reductions, improving accuracy of CBCT-based proton range/dose calculations.

Linear energy transfer weighted beam orientation optimization for intensity-modulated proton therapy

PURPOSE

In intensity-modulated proton therapy (IMPT), unaccounted-for variation in biological effectiveness contributes to the discrepancy between the constant relative biological effectiveness (RBE) model prediction and experimental observation. It is desirable to incorporate biological doses in treatment planning to improve modeling accuracy and consequently achieve a higher therapeutic ratio. This study addresses this demand by developing a method to incorporate linear energy transfer (LET) into beam orientation optimization (BOO).

METHODS

Instead of RBE-weighted dose, this LET weighted BOO (LETwBOO) framework uses the dose and LET product (LET  ×  D) as the biological surrogate. The problem is formulated with a physical dose fidelity term, a LET  ×  D constraint term, and a group sparsity term. The LET  ×  D of organs at risks is penalized for minimizing the biological effect while maintaining the physical dose objectives. Group sparsity is used to reduce the number of active beams from 600-800 non-coplanar candidate beams to between 2 and 4. This LETwBOO method was tested on three skull base tumor (SBT) patients and three bilateral head-and-neck (H&N) patients. The LETwBOO plans were compared with IMPT plans using manually selected beams with only physical dose constraint (MAN) and the initial MAN plan reoptimized with additional LET  ×  D constraint (LETwMAN).

RESULTS

The LETwBOO plans show superior physical dose and LET  ×  D sparing. On average, the [mean, maximal] doses of organs at risks (OARs) in LETwBOO are reduced by [2.85, 4.6] GyRBE from the MAN plans in the SBT cases and reduced by [0.9, 2.5] GyRBE in the H&N cases, while LETwMAN is comparable to MAN. cLET × Ds of PTVs are comparable in LETwBOO and LETwMAN, where c is a scaling factor of 0.04 μm/keV. On average, in the SBT cases, LETwBOO reduces the OAR [mean, maximal] cLET × D by [1.1, 2.9] Gy from the MAN plans, compared to the reduction by LETwMAN from MAN of [0.7, 1.7] Gy. In the H&N cases, LETwBOO reduces the OAR [mean, maximal] cLET × D by [0.8, 2.6] Gy from the MAN plans, compared to the reduction by LETwMAN from MAN of [0.3, 1.2] Gy.

CONCLUSION

We developed a novel LET weighted BOO method for IMPT to generate plans with improved physical and biological OAR sparing compared with the plans unaccounted for biological effects from BOO.

MR-guided proton therapy: a review and a preview

Background

The targeting accuracy of proton therapy (PT) for moving soft-tissue tumours is expected to greatly improve by real-time magnetic resonance imaging (MRI) guidance. The integration of MRI and PT at the treatment isocenter would offer the opportunity of combining the unparalleled soft-tissue contrast and real-time imaging capabilities of MRI with the most conformal dose distribution and best dose steering capability provided by modern PT. However, hybrid systems for MR-integrated PT (MRiPT) have not been realized so far due to a number of hitherto open technological challenges. In recent years, various research groups have started addressing these challenges and exploring the technical feasibility and clinical potential of MRiPT. The aim of this contribution is to review the different aspects of MRiPT, to report on the status quo and to identify important future research topics.

Methods

Four aspects currently under study and their future directions are discussed: modelling and experimental investigations of electromagnetic interactions between the MRI and PT systems, integration of MRiPT workflows in clinical facilities, proton dose calculation algorithms in magnetic fields, and MRI-only based proton treatment planning approaches.

Conclusions

Although MRiPT is still in its infancy, significant progress on all four aspects has been made, showing promising results that justify further efforts for research and development to be undertaken. First non-clinical research solutions have recently been realized and are being thoroughly characterized. The prospect that first prototype MRiPT systems for clinical use will likely exist within the next 5 to 10 years seems realistic, but requires significant work to be performed by collaborative efforts of research groups and industrial partners.

The TOPAS tool for particle simulation, a Monte Carlo simulation tool for physics, biology and clinical research

PURPOSE

This paper covers recent developments and applications of the TOPAS TOol for PArticle Simulation and presents the approaches used to disseminate TOPAS.

MATERIALS AND METHODS

Fundamental understanding of radiotherapy and imaging is greatly facilitated through accurate and detailed simulation of the passage of ionizing radiation through apparatus and into a patient using Monte Carlo (MC). TOPAS brings Geant4, a reliable, experimentally validated MC tool mainly developed for high energy physics, within easy reach of medical physicists, radiobiologists and clinicians. Requiring no programming knowledge, TOPAS provides all of the flexibility of Geant4.

RESULTS

After 5 years of development followed by its initial release, TOPAS was subsequently expanded from its focus on proton therapy physics to incorporate radiobiology modeling. Next, in 2018, the developers expanded their user support and code maintenance as well as the scope of TOPAS towards supporting X-ray and electron therapy and medical imaging. Improvements have been achieved in user enhancement through software engineering and a graphical user interface, calculational efficiency, validation through experimental benchmarks and QA measurements, and either newly available or recently published applications. A large and rapidly increasing user base demonstrates success in our approach to dissemination of this uniquely accessible and flexible MC research tool.

CONCLUSIONS

The TOPAS developers continue to make strides in addressing the needs of the medical community in applications of ionizing radiation to medicine, creating the only fully integrated platform for four-dimensional simulation of all forms of radiotherapy and imaging with ionizing radiation, with a design that promotes inter-institutional collaboration.

Intensity modulated proton therapy plan generation in under ten seconds

Background: Treatment planning for intensity modulated proton therapy (IMPT) can be significantly improved by reducing the time for plan calculation, facilitating efficient sampling of the large solution space characteristic of IMPT treatments. Additionally, fast plan generation is a key for online adaptive treatments, where the adapted plan needs to be ideally available in a few seconds. However, plan generation is a computationally demanding task and, although dose restoration methods for adaptive therapy have been proposed, computation times remain problematic. Material and methods: IMPT plan generation times were reduced by the development of dedicated graphical processing unit (GPU) kernels for our in-house, clinically validated, dose and optimization algorithms. The kernels were implemented into a coherent system, which performed all steps required for a complete treatment plan generation. Results: Using a single GPU, our fast implementation was able to generate a complete new treatment plan in 5-10 sec for typical IMPT cases, and in under 25 sec for plans to very large volumes such as for cranio-spinal axis irradiations. Although these times did not include the manual input of optimization parameters or a final clinical dose calculation, they included all required computational steps, including reading of CT and beam data. In addition, no compromise was made on plan quality. Target coverage and homogeneity for four patient plans improved (by up to 6%) or remained the same (changes <1%). No worsening of dose-volume parameters of the relevant organs at risk by more than 0.5% was observed. Conclusions: Fast plan generation with a clinically validated dose calculation and optimizer is a promising approach for daily adaptive proton therapy, as well as for automated or highly interactive planning.

Treatment planning for proton therapy: what is needed in the next 10 years?

Treatment planning is the process where the prescription of the radiation oncologist is translated into a deliverable treatment. With the complexity of contemporary radiotherapy, treatment planning cannot be performed without a computerized treatment planning system. Proton therapy (PT) enables highly conformal treatment plans with a minimum of dose to tissues outside the target volume, but to obtain the most optimal plan for the treatment, there are a multitude of parameters that need to be addressed. In this review areas of ongoing improvements and research in the field of PT treatment planning are identified and discussed. The main focus is on issues of immediate clinical and practical relevance to the PT community highlighting the needs for the near future but also in a longer perspective. We anticipate that the manual tasks performed by treatment planners in the future will involve a high degree of computational thinking, as many issues can be solved much better by e.g. scripting. More accurate and faster dose calculation algorithms are needed, automation for contouring and planning is required and practical tools to handle the variable biological efficiency in PT is urgently demanded just to mention a few of the expected improvements over the coming 10 years.

RIDOS: A new system for online computation of the delivered dose distributions in scanning ion beam therapy

PURPOSE

To describe a new system for scanned ion beam therapy, named RIDOS (Real-time Ion DOse planning and delivery System), which performs real time delivered dose verification integrating the information from a clinical beam monitoring system with a Graphic Processing Unit (GPU) based dose calculation in patient Computed Tomography.

METHODS

A benchmarked dose computation algorithm for scanned ion beams has been parallelized and adapted to run on a GPU architecture. A workstation equipped with a NVIDIA GPU has been interfaced through a National Instruments PXI-crate with the dose delivery system of the Italian National Center of Oncological Hadrontherapy (CNAO) to receive in real-time the measured beam parameters. Data from a patient monitoring system are also collected to associate the respiratory phases with each spot during the delivery of the dose. Using both measured and planned spot properties, RIDOS evaluates during the few seconds of inter-spill time the cumulative delivered and prescribed dose distributions and compares them through a fast γ-index algorithm.

RESULTS

The accuracy of the GPU-based algorithms was assessed against the CPU-based ones and the differences were found below 1‰. The cumulative planned and delivered doses are computed at the end of each spill in about 300 ms, while the dose comparison takes approximatively 400 ms. The whole operation provides the results before the next spill starts.

CONCLUSIONS

RIDOS system is able to provide a fast computation of the delivered dose in the inter-spill time of the CNAO facility and allows to monitor online the dose deposition accuracy all along the treatment.

Robust optimization for intensity-modulated proton therapy with soft spot sensitivity regularization

PURPOSE

Proton dose distribution is sensitive to uncertainties in range estimation and patient positioning. Currently, the proton robustness is managed by worst-case scenario optimization methods, which are computationally inefficient. To overcome these challenges, we develop a novel intensity-modulated proton therapy (IMPT) optimization method that integrates dose fidelity with a sensitivity term that describes dose perturbation as the result of range and positioning uncertainties.

METHODS

In the integrated optimization framework, the optimization cost function is formulated to include two terms: a dose fidelity term and a robustness term penalizing the inner product of the scanning spot sensitivity and intensity. The sensitivity of an IMPT scanning spot to perturbations is defined as the dose distribution variation induced by range and positioning errors. To evaluate the sensitivity, the spatial gradient of the dose distribution of a specific spot is first calculated. The spot sensitivity is then determined by the total absolute value of the directional gradients of all affected voxels. The fast iterative shrinkage-thresholding algorithm is used to solve the optimization problem. This method was tested on three skull base tumor (SBT) patients and three bilateral head-and-neck (H&N) patients. The proposed sensitivity-regularized method (SenR) was implemented on both clinic target volume (CTV) and planning target volume (PTV). They were compared with conventional PTV-based optimization method (Conv) and CTV-based voxel-wise worst-case scenario optimization approach (WC).

RESULTS

Under the nominal condition without uncertainties, the three methods achieved similar CTV dose coverage, while the CTV-based SenR approach better spared organs at risks (OARs) compared with the WC approach, with an average reduction of [Dmean, Dmax] of [4.72, 3.38]  GyRBE for the SBT cases and [2.54, 3.33] GyRBE for the H&N cases. The OAR sparing of the PTV-based SenR method was comparable with the WC method. The WC method, and SenR approaches all improved the plan robustness from the conventional PTV-based method. On average, under range uncertainties, the lowest [D95%, V95%, V100%] of CTV were increased from [93.75%, 88.47%, 47.37%] in the Conv method, to [99.28%, 99.51%, 86.64%] in the WC method, [97.71%, 97.85%, 81.65%] in the SenR-CTV method and [98.77%, 99.30%, 85.12%] in the SenR-PTV method, respectively. Under setup uncertainties, the average lowest [D95%, V95%, V100%] of CTV were increased from [95.35%, 94.92%, 65.12%] in the Conv method, to [99.43%, 99.63%, 87.12%] in the WC method, [96.97%, 97.13%, 77.86%] in the SenR-CTV method, and [98.21%, 98.34%, 83.88%] in the SenR-PTV method, respectively. The runtime of the SenR optimization is eight times shorter than that of the voxel-wise worst-case method.

CONCLUSION

We developed a novel computationally efficient robust optimization method for IMPT. The robustness is calculated as the spot sensitivity to both range and shift perturbations. The dose fidelity term is then regularized by the sensitivity term for the flexibility and trade-off between the dosimetry and the robustness. In the stress test, SenR is more resilient to unexpected uncertainties. These advantages in combination with its fast computation time make it a viable candidate for clinical IMPT planning.

A full-scale clinical prototype for proton range verification using prompt gamma-ray spectroscopy

We present a full-scale clinical prototype system for in vivo range verification of proton pencil-beams using the prompt gamma-ray spectroscopy method. The detection system consists of eight LaBr3 scintillators and a tungsten collimator, mounted on a rotating frame. Custom electronics and calibration algorithms have been developed for the measurement of energy- and time-resolved gamma-ray spectra during proton irradiation at a clinical dose rate. Using experimentally determined nuclear reaction cross sections and a GPU-accelerated Monte Carlo simulation, a detailed model of the expected gamma-ray emissions is created for each individual pencil-beam. The absolute range of the proton pencil-beams is determined by minimizing the discrepancy between the measurement and this model, leaving the absolute range of the beam and the elemental concentrations of the irradiated matter as free parameters. The system was characterized in a clinical-like situation by irradiating different phantoms with a scanning pencil-beam. A dose of 0.9 Gy was delivered to a [Formula: see text] cm3 target with a beam current of 2 nA incident on the phantom. Different range shifters and materials were used to test the robustness of the verification method and to calculate the accuracy of the detected range. The absolute proton range was determined for each spot of the distal energy layer with a mean statistical precision of 1.1 mm at a 95% confidence level and a mean systematic deviation of 0.5 mm, when aggregating pencil-beam spots within a cylindrical region of 10 mm radius and 10 mm depth. Small range errors that we introduced were successfully detected and even large differences in the elemental composition do not affect the range verification accuracy. These results show that our system is suitable for range verification during patient treatments in our upcoming clinical study.

Density overwrites of internal tumor volumes in intensity modulated proton therapy plans for mobile lung tumors

The purpose of this study was to investigate internal tumor volume density overwrite strategies to minimize intensity modulated proton therapy (IMPT) plan degradation of mobile lung tumors. Four planning paradigms were compared for nine lung cancer patients. Internal gross tumor volume (IGTV) and internal clinical target volume (ICTV) structures were defined encompassing their respective volumes in every 4DCT phase. The paradigms use different planning CT (pCT) created from the average intensity projection (AIP) of the 4DCT, overwriting the density within the IGTV to account for movement. The density overwrites were: (a) constant filling with 100 HU (C100) or (b) 50 HU (C50), (c) maximum intensity projection (MIP) across phases, and (d) water equivalent path length (WEPL) consideration from beam's-eye-view. Plans were created optimizing dose-influence matrices calculated with fast GPU Monte Carlo (MC) simulations in each pCT. Plans were evaluated with MC on the 4DCTs using a model of the beam delivery time structure. Dose accumulation was performed using deformable image registration. Interplay effect was addressed applying 10 times rescanning. Significantly less DVH metrics degradation occurred when using MIP and WEPL approaches. Target coverage ([Formula: see text] Gy(RBE)) was fulfilled in most cases with MIP and WEPL ([Formula: see text] Gy (RBE)), keeping dose heterogeneity low ([Formula: see text] Gy(RBE)). The mean lung dose was kept lowest by the WEPL strategy, as well as the maximum dose to organs at risk (OARs). The impact on dose levels in the heart, spinal cord and esophagus were patient specific. Overall, the WEPL strategy gives the best performance and should be preferred when using a 3D static geometry for lung cancer IMPT treatment planning. Newly available fast MC methods make it possible to handle long simulations based on 4D data sets to perform studies with high accuracy and efficiency, even prior to individual treatment planning.

Efficiency improvement in proton dose calculations with an equivalent restricted stopping power formalism

The equivalent restricted stopping power formalism is introduced for proton mean energy loss calculations under the continuous slowing down approximation. The objective is the acceleration of Monte Carlo dose calculations by allowing larger steps while preserving accuracy. The fractional energy loss per step length ϵ was obtained with a secant method and a Gauss-Kronrod quadrature estimation of the integral equation relating the mean energy loss to the step length. The midpoint rule of the Newton-Cotes formulae was then used to solve this equation, allowing the creation of a lookup table linking ϵ to the equivalent restricted stopping power L _eq, used here as a key physical quantity. The mean energy loss for any step length was simply defined as the product of the step length with L _eq. Proton inelastic collisions with electrons were added to GPUMCD, a GPU-based Monte Carlo dose calculation code. The proton continuous slowing-down was modelled with the L _eq formalism. GPUMCD was compared to Geant4 in a validation study where ionization processes alone were activated and a voxelized geometry was used. The energy straggling was first switched off to validate the L _eq formalism alone. Dose differences between Geant4 and GPUMCD were smaller than 0.31% for the L _eq formalism. The mean error and the standard deviation were below 0.035% and 0.038% respectively. 99.4 to 100% of GPUMCD dose points were consistent with a 0.3% dose tolerance. GPUMCD 80% falloff positions ([Formula: see text]) matched Geant's [Formula: see text] within 1 μm. With the energy straggling, dose differences were below 2.7% in the Bragg peak falloff and smaller than 0.83% elsewhere. The [Formula: see text] positions matched within 100 μm. The overall computation times to transport one million protons with GPUMCD were 31-173 ms. Under similar conditions, Geant4 computation times were 1.4-20 h. The L _eq formalism led to an intrinsic efficiency gain factor ranging between 30-630, increasing with the prescribed accuracy of simulations. The L _eq formalism allows larger steps leading to a [Formula: see text] algorithmic time complexity. It significantly accelerates Monte Carlo proton transport while preserving accuracy. It therefore constitutes a promising variance reduction technique for computing proton dose distributions in a clinical context.

Experimental assessment of proton dose calculation accuracy in inhomogeneous media

PURPOSE

Proton therapy with Pencil Beam Scanning (PBS) has the potential to improve radiotherapy treatments. Unfortunately, its promises are jeopardized by the sensitivity of the dose distributions to uncertainties, including dose calculation accuracy in inhomogeneous media. Monte Carlo dose engines (MC) are expected to handle heterogeneities better than analytical algorithms like the pencil-beam convolution algorithm (PBA). In this study, an experimental phantom has been devised to maximize the effect of heterogeneities and to quantify the capability of several dose engines (MC and PBA) to handle these.

METHODS

An inhomogeneous phantom made of water surrounding a long insert of bone tissue substitute (1×10×10 cm³) was irradiated with a mono-energetic PBS field (10×10 cm²). A 2D ion chamber array (MatriXX, IBA Dosimetry GmbH) lied right behind the bone. The beam energy was such that the expected range of the protons exceeded the detector position in water and did not attain it in bone. The measurement was compared to the following engines: Geant4.9.5, PENH, MCsquare, as well as the MC and PBA algorithms of RayStation (RaySearch Laboratories AB).

RESULTS

For a γ-index criteria of 2%/2mm, the passing rates are 93.8% for Geant4.9.5, 97.4% for PENH, 93.4% for MCsquare, 95.9% for RayStation MC, and 44.7% for PBA. The differences in γ-index passing rates between MC and RayStation PBA calculations can exceed 50%.

CONCLUSION

The performance of dose calculation algorithms in highly inhomogeneous media was evaluated in a dedicated experiment. MC dose engines performed overall satisfactorily while large deviations were observed with PBA as expected.

Initial development of goCMC: a GPU-oriented fast cross-platform Monte Carlo engine for carbon ion therapy

Monte Carlo (MC) simulation is considered as the most accurate method for calculation of absorbed dose and fundamental physics quantities related to biological effects in carbon ion therapy. To improve its computational efficiency, we have developed a GPU-oriented fast MC package named goCMC, for carbon therapy. goCMC simulates particle transport in voxelized geometry with kinetic energy up to 450 MeV u-1. Class II condensed history simulation scheme with a continuous slowing down approximation was employed. Energy straggling and multiple scattering were modeled. δ-electrons were terminated with their energy locally deposited. Four types of nuclear interactions were implemented in goCMC, i.e. carbon-hydrogen, carbon-carbon, carbon-oxygen and carbon-calcium inelastic collisions. Total cross section data from Geant4 were used. Secondary particles produced in these interactions were sampled according to particle yield with energy and directional distribution data derived from Geant4 simulation results. Secondary charged particles were transported following the condensed history scheme, whereas secondary neutral particles were ignored. goCMC was developed under OpenCL framework and is executable on different platforms, e.g. GPU and multi-core CPU. We have validated goCMC with Geant4 in cases with different beam energy and phantoms including four homogeneous phantoms, one heterogeneous half-slab phantom, and one patient case. For each case [Formula: see text] carbon ions were simulated, such that in the region with dose greater than 10% of maximum dose, the mean relative statistical uncertainty was less than 1%. Good agreements for dose distributions and range estimations between goCMC and Geant4 were observed. 3D gamma passing rates with 1%/1 mm criterion were over 90% within 10% isodose line except in two extreme cases, and those with 2%/1 mm criterion were all over 96%. Efficiency and code portability were tested with different GPUs and CPUs. Depending on the beam energy and voxel size, the computation time to simulate [Formula: see text] carbons was 9.9-125 s, 2.5-50 s and 60-612 s on an AMD Radeon GPU card, an NVidia GeForce GTX 1080 GPU card and an Intel Xeon E5-2640 CPU, respectively. The combined accuracy, efficiency and portability make goCMC attractive for research and clinical applications in carbon ion therapy.