Clinical validation of a GPU-based Monte Carlo dose engine of a commercial treatment planning system for pencil beam scanning proton therapy

NRG Oncology and PTCOG Patterns of Practice Survey and Consensus Recommendations on Pencil-Beam Scanning Proton Stereotactic Body Radiation Therapy and Hypofractionated Radiation Therapy for Thoracic Malignancies

Stereotactic body radiation therapy (SBRT) and hypofractionation using pencil-beam scanning (PBS) proton therapy (PBSPT) is an attractive option for thoracic malignancies. Combining the advantages of target coverage conformity and critical organ sparing from both PBSPT and SBRT, this new delivery technique has great potential to improve the therapeutic ratio, particularly for tumors near critical organs. Safe and effective implementation of PBSPT SBRT/hypofractionation to treat thoracic malignancies is more challenging than the conventionally fractionated PBSPT because of concerns of amplified uncertainties at the larger dose per fraction. The NRG Oncology and Particle Therapy Cooperative Group Thoracic Subcommittee surveyed proton centers in the United States to identify practice patterns of thoracic PBSPT SBRT/hypofractionation. From these patterns, we present recommendations for future technical development of proton SBRT/hypofractionation for thoracic treatment. Among other points, the recommendations highlight the need for volumetric image guidance and multiple computed tomography-based robust optimization and robustness tools to minimize further the effect of uncertainties associated with respiratory motion. Advances in direct motion analysis techniques are urgently needed to supplement current motion management techniques.

Modelling small block aperture in an in-house developed GPU-accelerated Monte Carlo-based dose engine for pencil beam scanning proton therapy

Purpose. To enhance an in-house graphic-processing-unit accelerated virtual particle (VP)-based Monte Carlo (MC) proton dose engine (VPMC) to model aperture blocks in both dose calculation and optimization for pencil beam scanning proton therapy (PBSPT)-based stereotactic radiosurgery (SRS).Methods and materials. A module to simulate VPs passing through patient-specific aperture blocks was developed and integrated in VPMC based on simulation results of realistic particles (primary protons and their secondaries). To validate the aperture block module, VPMC was first validated by an opensource MC code, MCsquare, in eight water phantom simulations with 3 cm thick brass apertures: four were with aperture openings of 1, 2, 3, and 4 cm without a range shifter, while the other four were with same aperture opening configurations with a range shifter of 45 mm water equivalent thickness. Then, VPMC was benchmarked with MCsquare and RayStation MC for 10 patients with small targets (average volume 8.4 c.c. with range of 0.4-43.3 c.c.). Finally, 3 typical patients were selected for robust optimization with aperture blocks using VPMC.Results. In the water phantoms, 3D gamma passing rate (2%/2 mm/10%) between VPMC and MCsquare was 99.71 ± 0.23%. In the patient geometries, 3D gamma passing rates (3%/2 mm/10%) between VPMC/MCsquare and RayStation MC were 97.79 ± 2.21%/97.78 ± 1.97%, respectively. Meanwhile, the calculation time was drastically decreased from 112.45 ± 114.08 s (MCsquare) to 8.20 ± 6.42 s (VPMC) with the same statistical uncertainties of ~0.5%. The robustly optimized plans met all the dose-volume-constraints (DVCs) for the targets and OARs per our institutional protocols. The mean calculation time for 13 influence matrices in robust optimization by VPMC was 41.6 s and the subsequent on-the-fly 'trial-and-error' optimization procedure took only 71.4 s on average for the selected three patients.Conclusion. VPMC has been successfully enhanced to model aperture blocks in dose calculation and optimization for the PBSPT-based SRS.

A GPU-based fast Monte Carlo code that supports proton transport in magnetic field for radiation therapy

This paper presents the effort to extend a previously reported code ARCHER, a GPU-based Monte Carlo (MC) code for coupled photon and electron transport, into protons including the consideration of magnetic fields. The proton transport is modeled using a Class-II condensed-history algorithm with continuous slowing-down approximation. The model includes ionization, multiple scattering, energy straggling, elastic and inelastic nuclear interactions, as well as deflection due to the Lorentz force in magnetic fields. An additional direction change is added for protons at the end of each step in the presence of the magnetic field. Secondary charge particles, except for protons, are terminated depositing kinetic energies locally, whereas secondary neutral particles are ignored. Each proton is transported step by step until its energy drops to below 0.5 MeV or when the proton leaves the phantom. The code is implemented using the compute unified device architecture (CUDA) platform for optimized GPU thread-level parallelism and efficiency. The code is validated by comparing it against TOPAS. Comparisons of dose distributions between our code and TOPAS for several exposure scenarios, ranging from single square beams in water to patient plan with magnetic fields, show good agreement. The 3D-gamma pass rate with a 2 mm/2% criterion in the region with dose greater than 10% of the maximum dose is computed to be over 99% for all tested cases. Using a single NVIDIA TITAN V GPU card, the computational time of ARCHER is found to range from 0.82 to 4.54 seconds for 1 × 107 proton histories. Compared to a few hours running on TOPAS, this speed improvement is significant. This work presents, for the first time, the performance of a GPU-based MC code to simulate proton transportation magnetic fields, demonstrating the feasibility of accurate and efficient dose calculations in potential magnetic resonance imaging (MRI)-guided proton therapy.

Treatment planning of scanned proton beams in RayStation

The use of scanned proton beams in external beam radiation therapy has seen a rapid development over the past decade. This technique places new demands on treatment planning, as compared to conventional photon-based radiation therapy. In this article, several proton specific functions as implemented in the treatment planning system RayStation are presented. We will cover algorithms for energy layer and spot selection, basic optimization including the handling of spot weight limits, optimization of the linear energy transfer (LET) distribution, robust optimization including the special case of 4D optimization, proton arc planning, and automatic planning using deep learning. We will further present the Monte Carlo (MC) proton dose engine in RayStation to some detail, from the material interpretation of the CT data, through the beam model parameterization, to the actual MC transport mechanism. Useful tools for plan evaluation, including robustness evaluation, and the versatile scripting interface are also described. The overall aim of the paper is to give an overview of some of the key proton planning functions in RayStation, with example usages, and at the same time provide the details about the underlying algorithms that previously have not been fully publicly available.

An optimized framework for cone-beam computed tomography-based online evaluation for proton therapy

BACKGROUND

Clinical evidence has demonstrated that proton therapy can achieve comparable tumor control probabilities compared to conventional photon therapy but with the added benefit of sparing healthy tissues. However, proton therapy is sensitive to inter-fractional anatomy changes. Online pre-fraction evaluation can effectively verify proton dose before delivery to patients, but there is a lack of guidelines for implementing this workflow.

PURPOSE

The purpose of this study is to develop a cone-beam CT-based (CBCT) online evaluation framework for proton therapy that enables knowledge transparency and evaluates the efficiency and accuracy of each essential component.

METHODS

Twenty-three patients with various lesion sites were included to conduct a retrospective study of implementing the proposed CBCT evaluation framework for the clinic. The framework was implemented on the RayStation 11B Research platform. Two synthetic CT (sCT) methods, corrected CBCT (cCBCT), and virtual CT (vCT), were used, and the ground truth images were acquired from the same-day deformed quality assurance CT (dQACT) for the comparisons. The evaluation metrics for the framework include time efficiency, dose-difference distributions (gamma passing rates), and water equivalent thickness (WET) distributions.

RESULTS

The mean online CBCT evaluation times were 1.6 ± 0.3 min and 1.9 ± 0.4 min using cCBCT and vCT, respectively. The dose calculation and deformable image registration dominated the evaluation efficiency, and accounted for 33% and 30% of the total evaluation time, respectively. The sCT generation took another 19% of the total time. Gamma passing rates were greater than 91% and 97% using 1%/1 mm and 2%/2 mm criteria, respectively. When the appropriate sCT was chosen, the target mean WET difference from the reference were less than 0.5 mm. The appropriate sCT method choice determined the uncertainty for the framework, with the cCBCT being superior for head-and-neck patient evaluation and vCT being better for lung patient evaluation.

CONCLUSIONS

An online CBCT evaluation framework was proposed to identify the use of the optimal sCT algorithm regarding efficiency and dosimetry accuracy. The framework is extendable to adopt advanced imaging methods and has the potential to support online adaptive radiotherapy to enhance patient benefits. It could be implemented into clinical use in the future.

Validation of RayStation Monte Carlo dose calculation algorithm for multiple LINACs

PURPOSE

A photon Monte Carlo (MC) model was commissioned for flattened (FF) and flattening filter free (FFF) 6 MV beam energy. The accuracy of this model, as a single model to be used for three beam matched LINACs, was evaluated.

METHODS

Multiple models were created in RayStation v.10A for three linacs equipped with Elekta "Agility" collimator. A clinically commissioned collapsed cone (CC) algorithm (GoldCC), a MC model automatically created from the CC algorithm without further optimization (CCtoMC) and an optimized MC model (GoldMC) were compared with measurements. The validation of the model was performed by following the recommendations of IAEA TRS 430 and comprised of basic validation in a water tank, validation in a heterogeneous phantom and validation of complex IMRT/VMAT paradigms using gamma analysis of calculated and measured dose maps in a 2D-Array.

RESULTS

Dose calculation with the GoldMC model resulted in a confidence level of 3% for point measurements in water tank and heterogeneous phantom for measurements performed in all three linacs. The same confidence level resulted for GoldCC model. Dose maps presented an agreement for all models on par to each other with γ criteria 2%/2mm.

CONCLUSIONS

The GoldMC model showed a good agreement with measured data and is determined to be accurate for clinical use for all three linacs in this study.

Virtual particle Monte Carlo: A new concept to avoid simulating secondary particles in proton therapy dose calculation

BACKGROUND

In proton therapy dose calculation, Monte Carlo (MC) simulations are superior in accuracy but more time consuming, compared to analytical calculations. Graphic processing units (GPUs) are effective in accelerating MC simulations but may suffer thread divergence and racing condition in GPU threads that degrades the computing performance due to the generation of secondary particles during nuclear reactions.

PURPOSE

A novel concept of virtual particle (VP) MC (VPMC) is proposed to avoid simulating secondary particles in GPU-accelerated proton MC dose calculation and take full advantage of the computing power of GPU.

METHODS

Neutrons and gamma rays were ignored as escaping from the human body; doses of electrons, heavy ions, and nuclear fragments were locally deposited; the tracks of deuterons were converted into tracks of protons. These particles, together with primary and secondary protons, are considered to be the realistic particles. Histories of primary and secondary protons were replaced by histories of multiple VPs. Each VP corresponded to one proton (either primary or secondary). A continuous-slowing-down-approximation model, an ionization model, and a large angle scattering event model corresponding to nuclear interactions were developed for VPs by generating probability distribution functions (PDFs) based on simulation results of realistic particles using MCsquare. For efficient calculations, these PDFs were stored in the Compute Unified Device Architecture textures. VPMC was benchmarked with TOPAS and MCsquare in phantoms and with MCsquare in 13 representative patient geometries. Comparisons between the VPMC calculated dose and dose measured in water during patient-specific quality assurance (PSQA) of the selected 13 patients were also carried out. Gamma analysis was used to compare the doses derived from different methods and calculation efficiencies were also compared.

RESULTS

Integrated depth dose and lateral dose profiles in both homogeneous and inhomogeneous phantoms all matched well among VPMC, TOPAS, and MCsquare calculations. The 3D-3D gamma passing rates with a criterion of 2%/2 mm and a threshold of 10% was 98.49% between MCsquare and TOPAS and 98.31% between VPMC and TOPAS in homogeneous phantoms, and 99.18% between MCsquare and TOPAS and 98.49% between VPMC and TOPAS in inhomogeneous phantoms, respectively. In patient geometries, the 3D-3D gamma passing rates with 2%/2 mm/10% between dose distributions from VPMC and MCsquare were 98.56 ± 1.09% in patient geometries. The 2D-3D gamma analysis with 3%/2 mm/10% between the VPMC calculated dose distributions and the 2D measured planar dose distributions during PSQA was 98.91 ± 0.88%. VPMC calculation was highly efficient and took 2.84 ± 2.44 s to finish for the selected 13 patients running on four NVIDIA Ampere GPUs in patient geometries.

CONCLUSION

VPMC was found to achieve high accuracy and efficiency in proton therapy dose calculation.

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.

Fast robust optimization of proton PBS arc therapy plans using early energy layer selection and spot assignment

Objective. Proton pencil-beam scanning arcs (PBS arcs) have gained much attention during the past years, due to its potential for increased clinical benefit compared to conventional proton therapy. Previous studies on PBS arcs have primarily been focused on plan quality, and lately efforts have been made to reduce the delivery time. However, the methods presented so far suffer from slow optimization processes. Approach. We present a new method for fast robust optimization of PBS arc plans. The new method assigns a single energy layer per discretized direction prior to spot weight optimization and reduces the number of initial spots considerably compared to conventional methods. We used the new method for three prostate cancer patients with a prescribed dose to the CTV of 77 GyRBE in 35 fractions. For each of the patients, four plans were created: 2-beam IMPT (2IMPT), 1-beam PBS arc (1Arc), 1-beam PBS arc without focus on reducing upward energy jumps (1Arc_unseq) and two-beam PBS arc (2Arc). Main results. All PBS arc plans show a reduced integral dose compared to their respective 2IMPT plans. In the nominal case, the average CTV D98 and D2 metrics over the three patients were best for the 2Arc, followed by 2IMPT ( D 98 ¯ / D 2 ¯ : 7523/7986 cGyRBE (2IMPT), 7478/7984 cGy (1Arc), 7486/7951 cGy (1Arc_unseq), 7531/7951 cGyRBE (2Arc)). The average robust target coverage in terms of V95 of the voxelwise minimum dose distribution (evaluated over 42 scenarios) was: 98.0% (2IMPT), 88.6% (1Arc), 92.5% (1Arc_unseq), 97.3% (2Arc). The optimization time, including spot selection and spot dose computation, is longest for the 2Arc plan, but is below 6 min for all patients. The maximum estimated delivery time for all types of arc plans is just above 5 min Significance. The ability for efficient treatment planning constitutes an important step towards clinical introduction of proton PBS arcs.