Generalized eMC implementation for Monte Carlo dose calculation of electron beams from different machine types

Robust optimization and assessment of dynamic trajectory and mixed-beam arc radiotherapy: a preliminary study

Objective.Dynamic trajectory radiotherapy (DTRT) and dynamic mixed-beam arc therapy (DYMBARC) exploit non-coplanarity and, for DYMBARC, simultaneously optimized photon and electron beams. Margin concepts to account for set-up uncertainties during delivery are ill-defined for electron fields. We develop robust optimization for DTRT&DYMBARC and compare dosimetric plan quality and robustness for both techniques and both optimization strategies for four cases.Approach.Cases for different treatment sites and clinical target volume (CTV) to planning target volume (PTV) margins,m, were investigated. Dynamic gantry-table-collimator photon paths were optimized to minimize PTV/organ-at-risk (OAR) overlap in beam's-eye-view and minimize potential photon multileaf collimator (MLC) travel. For DYMBARC plans, non-isocentric partial electron arcs or static fields with shortened source-to-surface distance (80 cm) were added. Direct aperture optimization (DAO) was used to simultaneously optimize MLC-based intensity modulation for both photon and electron beams yielding deliverable PTV-based DTRT&DYMBARC plans. Robust-optimized plans used the same paths/arcs/fields. DAO with stochastic programming was used for set-up uncertainties with equal weights in all translational directions and magnitudeδsuch thatm= 0.7δ. Robust analysis considered random errors in all directions with or without an additional systematic error in the worst 3D direction for the adjacent OARs.Main results.Electron contribution was 7%-41% of target dose depending on the case and optimization strategy for DYMBARC. All techniques achieved similar CTV coverage in the nominal (no error) scenario. OAR sparing was overall better in the DYMBARC plans than in the DTRT plans and DYMBARC plans were generally more robust to the considered uncertainties. OAR sparing was better in the PTV-based than in robust-optimized plans for OARs abutting or overlapping with the target volume, but more affected by uncertainties.Significance.Better plan robustness can be achieved with robust optimization than with margins. Combining electron arcs/fields with non-coplanar photon trajectories further improves robustness and OAR sparing.

Commissioning and initial validation of Eclipse eMC algorithm for the electron FLASH research extension (FLEX) system for pre-clinical studies

PURPOSE

To investigate the feasibility of commissioning the 16 MeV electron FLASH Extension (FLEX) in the commercial treatment planning system (TPS) for biomedical research with cell and mouse models, and in silico treatment planning studies.

METHODS

To commission the FLEX system with the electron Monte Carlo (eMC) algorithm in the commercial TPS, radiochromic film was used to measure the vendor-recommended beam data. Once the beam model was generated for the eMC algorithm, supplemental measurements were collected for validation purposes and compared against the TPS-calculated results. Additionally, the newly commissioned 16 MeV FLASH beam was compared to the corresponding 16 MeV conventional electron beam.

RESULTS

The eMC algorithm effectively modeled the FLEX system. The eMC-calculated PDDs and profiles for the 16 MeV electron FLASH beam agreed with measured values within 1%, on average, for 6 × 6 cm2 and 10 × 10 cm2 applicators. Flatness and symmetry deviated by less than 1%, while FWHM and penumbra agreed within 1 mm for both eMC-calculated and measured profiles. Additionally, the small field (i.e., 2-cm diameter cutout) that was measured for validation purposes agreed with TPS-calculated results within 1%, on average, for both the PDD and profiles. The FLASH and conventional dose rate 16 MeV electron beam were in agreement in regard to energy, but the profiles for larger field sizes began to deviate (>10 × 10 cm2) due to the forward-peaked nature of the FLASH beam. For cell irradiation experiments, the measured and eMC-calculated in-plane and cross-plane absolute dose profiles agreed within 1%, on average.

CONCLUSIONS

The FLEX system was successfully commissioned in the commercial TPS using the eMC algorithm, which accurately modeled the forward-peaked nature of the FLASH beam. A commissioned TPS for FLASH will be useful for pre-clinical cell and animal studies, as well as in silico FLASH treatment planning studies for future clinical implementation.

Delivery time reduction for mixed photon-electron radiotherapy by using photon MLC collimated electron arcs

Objective. Electron arcs in mixed-beam radiotherapy (Arc-MBRT) consisting of intensity-modulated electron arcs with dynamic gantry rotation potentially reduce the delivery time compared to mixed-beam radiotherapy containing electron beams with static gantry angle (Static-MBRT). This study aims to develop and investigate a treatment planning process (TPP) for photon multileaf collimator (pMLC) based Arc-MBRT.Approach. An existing TPP for Static-MBRT plans is extended to integrate electron arcs with a dynamic gantry rotation and intensity modulation using a sliding window technique. The TPP consists of a manual setup of electron arcs, and either static photon beams or photon arcs, shortening of the source-to-surface distance for the electron arcs, initial intensity modulation optimization, selection of a user-defined number of electron beam energies based on dose contribution to the target volume and finally, simultaneous photon and electron intensity modulation optimization followed by full Monte Carlo dose calculation. Arc-MBRT plans, Static-MBRT plans, and photon-only plans were created and compared for four breast cases. Dosimetric validation of two Arc-MBRT plans was performed using film measurements.Main results. The generated Arc-MBRT plans are dosimetrically similar to the Static-MBRT plans while outperforming the photon-only plans. The mean heart dose is reduced by 32% on average in the MBRT plans compared to the photon-only plans. The estimated delivery times of the Arc-MBRT plans are similar to the photon-only plans but less than half the time of the Static-MBRT plans. Measured and calculated dose distributions agree with a gamma passing rate of over 98% (3% global, 2 mm) for both delivered Arc-MBRT plans.Significance. A TPP for Arc-MBRT is successfully developed and Arc-MBRT plans showed the potential to improve the dosimetric plan quality similar as Static-MBRT while maintaining short delivery times of photon-only treatments. This further facilitates integration of pMLC-based MBRT into clinical practice.

Auto-commissioning of a Monte Carlo electron beam model with application to photon MLC shaped electron fields

Objective.Presently electron beam treatments are delivered using dedicated applicators. An alternative is the usage of the already installed photon multileaf collimator (pMLC) enabling efficient electron treatments. Currently, the commissioning of beam models is a manual and time-consuming process. In this work an auto-commissioning procedure for the Monte Carlo (MC) beam model part representing the beam above the pMLC is developed for TrueBeam systems with electron energies from 6 to 22 MeV.Approach.The analytical part of the electron beam model includes a main source representing the primary beam and a jaw source representing the head scatter contribution each consisting of an electron and a photon component, while MC radiation transport is performed for the pMLC. The auto-commissioning of this analytical part relies on information pre-determined from MC simulations, in-air dose profiles and absolute dose measurements in water for different field sizes and source to surface distances (SSDs). For validation calculated and measured dose distributions in water were compared for different field sizes, SSDs and beam energies for eight TrueBeam systems. Furthermore, a sternum case in an anthropomorphic phantom was considered and calculated and measured dose distributions were compared at different SSDs.Main results.Instead of the manual commissioning taking up to several days of calculation time and several hours of user time, the auto-commissioning is carried out in a few minutes. Measured and calculated dose distributions agree generally within 3% of maximum dose or 2 mm. The gamma passing rates for the sternum case ranged from 96% to 99% (3% (global)/2 mm criteria, 10% threshold).Significance.The auto-commissioning procedure was successfully implemented and applied to eight TrueBeam systems. The newly developed user-friendly auto-commissioning procedure allows an efficient commissioning of an MC electron beam model and eases the usage of advanced electron radiotherapy utilizing the pMLC for beam shaping.

Efficiency enhancements of a Monte Carlo beamlet based treatment planning process: implementation and parameter study

Objective.The computational effort to perform beamlet calculation, plan optimization and final dose calculation of a treatment planning process (TPP) generating intensity modulated treatment plans is enormous, especially if Monte Carlo (MC) simulations are used for dose calculation. The goal of this work is to improve the computational efficiency of a fully MC based TPP for static and dynamic photon, electron and mixed photon-electron treatment techniques by implementing multiple methods and studying the influence of their parameters.Approach.A framework is implemented calculating MC beamlets efficiently in parallel on each available CPU core. The user can specify the desired statistical uncertainty of the beamlets, a fractional sparse dose threshold to save beamlets in a sparse format and minimal distances to the PTV surface from which 2 × 2 × 2 = 8 (medium) or even 4 × 4 × 4 = 64 (large) voxels are merged. The compromise between final plan quality and computational efficiency of beamlet calculation and optimization is studied for several parameter values to find a reasonable trade-off. For this purpose, four clinical and one academic case are considered with different treatment techniques.Main results.Setting the statistical uncertainty to 5% (photon beamlets) and 15% (electron beamlets), the fractional sparse dose threshold relative to the maximal beamlet dose to 0.1% and minimal distances for medium and large voxels to the PTV to 1 cm and 2 cm, respectively, does not lead to substantial degradation in final plan quality compared to using 2.5% (photon beamlets) and 5% (electron beamlets) statistical uncertainty and no sparse format nor voxel merging. Only OAR sparing is slightly degraded. Furthermore, computation times are reduced by about 58% (photon beamlets), 88% (electron beamlets) and 96% (optimization).Significance.Several methods are implemented improving computational efficiency of beamlet calculation and plan optimization of a fully MC based TPP without substantial degradation in final plan quality.

A hybrid column generation and simulated annealing algorithm for direct aperture optimization

The purpose of this work was to develop a hybrid column generation (CG) and simulated annealing (SA) algorithm for direct aperture optimization (H-DAO) and to show its effectiveness in generating high quality treatment plans for intensity modulated radiation therapy (IMRT) and mixed photon-electron beam radiotherapy (MBRT). The H-DAO overcomes limitations of the CG-DAO with two features improving aperture selection (branch-feature) and enabling aperture shape changes during optimization (SA-feature). The H-DAO algorithm iteratively adds apertures to the plan. At each iteration, a branch is created for each field provided. First, each branch determines the most promising aperture of its assigned field and adds it to a copy of the current apertures. Afterwards, the apertures of each branch undergo an MU-weight optimization followed by an SA-based simultaneous shape and MU-weight optimization and a second MU-weight optimization. The next H-DAO iteration continues the branch with the lowest objective function value. IMRT and MBRT treatment plans for an academic, a brain and a head and neck case generated using the CG-DAO and H-DAO were compared. For every investigated case and both IMRT and MBRT, the H-DAO leads to a faster convergence of the objective function value with number of apertures compared to the CG-DAO. In particular, the H-DAO needs about half the apertures to reach the same objective function value as the CG-DAO. The average aperture areas are 27% smaller for H-DAO than for CG-DAO leading to a slightly larger discrepancy between optimized and final dose. However, a dosimetric benefit remains. The H-DAO was successfully developed and applied to IMRT and MBRT. The faster convergence with number of apertures of the H-DAO compared to the CG-DAO allows to select a better compromise between plan quality and number of apertures.

Quantitative evaluation of dosimetric uncertainties in electron therapy by measurement and calculation using the electron Monte Carlo dose algorithm in the Eclipse treatment planning system

In the electron beam radiation therapy, customized blocks are mostly used to shape treatment fields to generate conformal doses. The goal of this study is to investigate quantitatively dosimetric uncertainties associated with heterogeneities, detectors used in the measurement of the beam data commissioning, and modeling of the interactions of high energy electrons with tissue. These uncertainties were investigated both by measurements with different detectors and calculations using electron Monte Carlo algorithm (eMC) in the Eclipse treatment planning system. Dose distributions for different field sizes were calculated using eMC and measured with a multiple-diode-array detector (MapCheck2) for cone sizes ranging from 6 to 25 cm. The dose distributions were calculated using the CT images of the MapCheck2 and water-equivalent phantoms. In the umbra region (<20% isodose line), the eMC underestimated dose by a factor of 3 for high energy electron beams due to lack of consideration of bremsstrahlung emitted laterally that was not accounted by eMC in the low dose region outside the field. In the penumbra (20%-80% isodose line), the eMC overestimated dose (40%) for high energy 20 MeV electrons compared to the measured dose with small diodes in the high gradient dose region. This was mainly due to lack of consideration of volume averaging of the ion chamber used in beam data commissioning which was input to the eMC dose calculation algorithm. Large uncertainties in the CT numbers (25%) resulted from the image artifacts in the CT images of the MapCheck2 phantom due to metal artifacts. The eMC algorithm used the electron and material densities extracted from the CT numbers which resulted large dosimetric uncertainties (10%) in the material densities and corresponding stopping power ratios. The dose calculations with eMC are associated with large uncertainties particularly in penumbra and umbra regions and around heterogeneities which affect the low dose level that cover nearby normal tissue or critical structures.

TriB-RT: Simultaneous optimization of photon, electron and proton beams

PURPOSE

To develop a novel treatment planning process (TPP) with simultaneous optimization of modulated photon, electron and proton beams for improved treatment plan quality in radiotherapy.

METHODS

A framework for fluence map optimization of Monte Carlo (MC) calculated beamlet dose distributions is developed to generate treatment plans consisting of photon, electron and spot scanning proton fields. Initially, in-house intensity modulated proton therapy (IMPT) plans are compared to proton plans created by a commercial treatment planning system (TPS). A triple beam radiotherapy (TriB-RT) plan is generated for an exemplary academic case and the dose contributions of the three particle types are investigated. To investigate the dosimetric potential, a TriB-RT plan is compared to an in-house IMPT plan for two clinically motivated cases. Benefits of TriB-RT for a fixed proton beam line with a single proton field are investigated.

RESULTS

In-house optimized IMPT are of at least equal or better quality than TPS-generated proton plans, and MC-based optimization shows dosimetric advantages for inhomogeneous situations. Concerning TriB-RT, for the academic case, the resulting plan shows substantial contribution of all particle types. For the clinically motivated case, improved sparing of organs at risk close to the target volume is achieved compared to IMPT (e.g. myelon and brainstem [Formula: see text] -37%) at cost of an increased low dose bath (healthy tissue V _10% +22%). In the scenario of a fixed proton beam line, TriB-RT plans are able to compensate the loss in degrees of freedom to substantially improve plan quality compared to a single field proton plan.

CONCLUSION

A novel TPP which simultaneously optimizes photon, electron and proton beams was successfully developed. TriB-RT shows the potential for improved treatment plan quality and is especially promising for cost-effective single-room proton solutions with a fixed beamline in combination with a conventional linac delivering photon and electron fields.

Development of an extended Macro Monte Carlo method for efficient and accurate dose calculation in magnetic fields

MOTIVATION

Progress in the field of magnetic resonance (MR)-guided radiotherapy has triggered the need for fast and accurate dose calculation in presence of magnetic fields. The aim of this work is to satisfy this need by extending the macro Monte Carlo (MMC) method to enable dose calculation for photon, electron, and proton beams in a magnetic field.

METHODS

The MMC method is based on the transport of particles in macroscopic steps through an absorber by sampling the relevant physical quantities from a precalculated database containing probability distribution functions. To enable MMC particle transport in a magnetic field, a transformation accounting for the Lorentz force is applied for each macro step by rotating the sampled position and direction around the magnetic field vector. The transformed position and direction distributions on local geometries are validated against full MC for electron and proton pencil beams. To enable photon dose calculation, an in-house MC algorithm is used for photon transport and interaction. Emerging secondary charged particles are passed to MMC for transport and energy deposition. The extended MMC dose calculation accuracy and efficiency is assessed by comparison with EGSnrc (photon and electron beams) and Geant4 (proton beam) calculated dose distributions of different energies and homogeneous magnetic fields for broad beams impinging on water phantoms with bone and lung inhomogeneities.

RESULTS

The geometric transformation on the local geometries is able to reproduce the results of full MC for all investigated settings (difference in mean value and standard deviation <1%). Macro Monte Carlo calculated dose distributions in a homogeneous magnetic field are in agreement with EGSnrc and Geant4, respectively, with gamma passing rates >99.6% (global 2%, 2 mm and 10% threshold criteria) for all situations. MMC achieves a substantial efficiency gain of up to a factor of 21 (photon beam), 66 (electron beam), and 356 (proton beam) compared to EGSnrc or Geant4.

CONCLUSION

Efficient and accurate dose calculation in magnetic fields was successfully enabled by utilizing the developed extended MMC transport method for photon, electron, and proton beams.

Part 2: Dynamic mixed beam radiotherapy (DYMBER): Photon dynamic trajectories combined with modulated electron beams

PURPOSE

The purpose of this study was to develop a treatment technique for dynamic mixed beam radiotherapy (DYMBER) utilizing increased degrees of freedom (DoF) of a conventional treatment unit including different particle types (photons and electrons), intensity and energy modulation and dynamic gantry, table, and collimator rotations.

METHODS

A treatment planning process has been developed to create DYMBER plans combining photon dynamic trajectories (DTs) and step and shoot electron apertures collimated with the photon multileaf collimator (pMLC). A gantry-table path is determined for the photon DTs with minimized overlap of the organs at risk (OARs) with the target. In addition, an associated dynamic collimator rotation is established with minimized area between the pMLC leaves and the target contour. pMLC sequences of photon DTs and electron pMLC apertures are then simultaneously optimized using direct aperture optimization (DAO). Subsequently, the final dose distribution of the electron pMLC apertures is calculated using the Swiss Monte Carlo Plan (SMCP). The pMLC sequences of the photon DTs are then re-optimized with a finer control point resolution and with the final electron dose distribution taken into account. Afterwards, the final photon dose distribution is calculated also using the SMCP and summed together with the one of the electrons. This process is applied for a brain and two head and neck cases. The resulting DYMBER dose distributions are compared to those of dynamic trajectory radiotherapy (DTRT) plans consisting only of photon DTs and clinically applied VMAT plans. Furthermore, the deliverability of the DYMBER plans is verified in terms of dosimetric accuracy, delivery time and collision avoidance. For this purpose, The DYMBER plans are delivered to Gafchromic EBT3 films placed in an anthropomorphic head phantom on a Varian TrueBeam linear accelerator.

RESULTS

For each case, the dose homogeneity in the target is similar or better for DYMBER compared to DTRT and VMAT. Averaged over all three cases, the mean dose to the parallel OARs is 16% and 28% lower, D2% to the serial OARs is 17% and 37% lower and V10% to normal tissue is 12% and 4% lower for the DYMBER plans compared to the DTRT and VMAT plans, respectively. The DYMBER plans are delivered without collision and with a 4-5 min longer delivery time than the VMAT plans. The absolute dose measurements are compared to calculation by gamma analysis using 2% (global)/2 mm criteria with passing rates of at least 99%.

CONCLUSIONS

A treatment technique for DYMBER has been successfully developed and verified for its deliverability. The dosimetric superiority of DYMBER over DTRT and VMAT indicates utilizing increased DoF to be the key to improve brain and head and neck radiation treatments in future.

Monte Carlo systems used for treatment planning and dose verification

General-purpose radiation transport Monte Carlo codes have been used for estimation of the absorbed dose distribution in external photon and electron beam radiotherapy patients since several decades. Results obtained with these codes are usually more accurate than those provided by treatment planning systems based on non-stochastic methods. Traditionally, absorbed dose computations based on general-purpose Monte Carlo codes have been used only for research, owing to the difficulties associated with setting up a simulation and the long computation time required. To take advantage of radiation transport Monte Carlo codes applied to routine clinical practice, researchers and private companies have developed treatment planning and dose verification systems that are partly or fully based on fast Monte Carlo algorithms. This review presents a comprehensive list of the currently existing Monte Carlo systems that can be used to calculate or verify an external photon and electron beam radiotherapy treatment plan. Particular attention is given to those systems that are distributed, either freely or commercially, and that do not require programming tasks from the end user. These systems are compared in terms of features and the simulation time required to compute a set of benchmark calculations.

Full Monte Carlo and measurement-based overall performance assessment of improved clinical implementation of eMC algorithm with emphasis on lower energy range

New version 13.6.23 of the electron Monte Carlo (eMC) algorithm in Varian Eclipse™ treatment planning system has a model for 4MeV electron beam and some general improvements for dose calculation. This study provides the first overall accuracy assessment of this algorithm against full Monte Carlo (MC) simulations for electron beams from 4MeV to 16MeV with most emphasis on the lower energy range. Beams in a homogeneous water phantom and clinical treatment plans were investigated including measurements in the water phantom. Two different material sets were used with full MC: (1) the one applied in the eMC algorithm and (2) the one included in the Eclipse™ for other algorithms. The results of clinical treatment plans were also compared to those of the older eMC version 11.0.31. In the water phantom the dose differences against the full MC were mostly less than 3% with distance-to-agreement (DTA) values within 2mm. Larger discrepancies were obtained in build-up regions, at depths near the maximum electron ranges and with small apertures. For the clinical treatment plans the overall dose differences were mostly within 3% or 2mm with the first material set. Larger differences were observed for a large 4MeV beam entering curved patient surface with extended SSD and also in regions of large dose gradients. Still the DTA values were within 3mm. The discrepancies between the eMC and the full MC were generally larger for the second material set. The version 11.0.31 performed always inferiorly, when compared to the 13.6.23.

Deterministic model for the transport of energetic particles: Application in the electron radiotherapy

A new deterministic method for calculating the dose distribution in the electron radiotherapy field is presented. The aim of this work was to validate our model by comparing it with the Monte Carlo simulation toolkit, GEANT4. A comparison of the longitudinal and transverse dose deposition profiles and electron distributions in homogeneous water phantoms showed a good accuracy of our model for electron transport, while reducing the calculation time by a factor of 50. Although the Bremsstrahlung effect is not yet implemented in our model, we propose here a method that solves the Boltzmann kinetic equation and provides a viable and efficient alternative to the expensive Monte Carlo modeling.