Position-probability-sampled Monte Carlo calculation of VMAT, 3DCRT, step-shoot IMRT, and helical tomotherapy dose distributions using BEAMnrc/DOSXYZnrc

A fast jaw-tracking model for VMAT and IMRT Monte Carlo simulations

Modern radiotherapy techniques involve routine use of volumetric arc therapy (VMAT) and intensity modulated radiotherapy (IMRT) with jaw‐tracking – dynamic motion of the secondary collimators (jaws) in tandem with multi‐leaf collimators (MLCs). These modalities require accurate dose calculations for the purposes of treatment planning and dose verification. Monte Carlo (MC) methods for radiotherapy dose calculation are widely accepted as capable of achieving high accuracy. This paper presents an efficiency‐enhancement method for secondary collimator modeling, presented in the context of a tool for MC‐based dose second checks. The model constitutes an accuracy trade‐off in the source model for the sake of efficiency enhancement, but maintains the advantages of MC transport in patient heterogeneities. The secondary collimator model is called Flat‐Absorbing‐Jaw‐Tracking (FAJT). Transmission through and scatter from the secondary collimators is neglected, and jaws are modeled as perfectly absorbing planes. To couple the motion of secondary collimators with MLCs for jaw‐tracking, the FAJT model was built into the VCU‐MLC model. Gamma‐index analysis of the dose distributions from FAJT against the full BEAMnrc MC simulations showed over 99% pass rate for a range of open fields, two clinical IMRT, and one VMAT treatment plan, for 2%/2 mm criteria above 10%. Using FAJT, the simulation speed of the secondary collimators for open fields increased by a factor of 237, 1489, and 1395 for 4 × 4, 10 × 10, and 30 × 30 cm², respectively. In general, clinically oriented simulation times are reduced from “hours” to “minutes” on identical hardware. Results for nine representative clinical cases (seven with jaw‐tracking) are presented. The average 2%/2 mm γ‐test success rate above the 80% isodose was 96.8% when tested against the EPIDose electronic portal image‐based dose reconstruction method and 97.3% against the Eclipse analytical anisotropic algorithm.

Experimental Validation of Monte Carlo Simulations Based on a Virtual Source Model for TomoTherapy in a RANDO Phantom

A virtual source model for Monte Carlo simulations of helical TomoTherapy has been developed previously by the authors. The purpose of this work is to perform experiments in an anthropomorphic (RANDO) phantom with the same order of complexity as in clinical treatments to validate the virtual source model to be used for quality assurance secondary check on TomoTherapy patient planning dose. Helical TomoTherapy involves complex delivery pattern with irregular beam apertures and couch movement during irradiation. Monte Carlo simulation, as the most accurate dose algorithm, is desirable in radiation dosimetry. Current Monte Carlo simulations for helical TomoTherapy adopt the full Monte Carlo model, which includes detailed modeling of individual machine component, and thus, large phase space files are required at different scoring planes. As an alternative approach, we developed a virtual source model without using the large phase space files for the patient dose calculations previously. In this work, we apply the simulation system to recompute the patient doses, which were generated by the treatment planning system in an anthropomorphic phantom to mimic the real patient treatments. We performed thermoluminescence dosimeter point dose and film measurements to compare with Monte Carlo results. Thermoluminescence dosimeter measurements show that the relative difference in both Monte Carlo and treatment planning system is within 3%, with the largest difference less than 5% for both the test plans. The film measurements demonstrated 85.7% and 98.4% passing rate using the 3 mm/3% acceptance criterion for the head and neck and lung cases, respectively. Over 95% passing rate is achieved if 4 mm/4% criterion is applied. For the dose–volume histograms, very good agreement is obtained between the Monte Carlo and treatment planning system method for both cases. The experimental results demonstrate that the virtual source model Monte Carlo system can be a viable option for the accurate dose calculation of helical TomoTherapy.

A review of the use and potential of the GATE Monte Carlo simulation code for radiation therapy and dosimetry applications

In this paper, the authors' review the applicability of the open-source GATE Monte Carlo simulation platform based on the GEANT4 toolkit for radiation therapy and dosimetry applications. The many applications of GATE for state-of-the-art radiotherapy simulations are described including external beam radiotherapy, brachytherapy, intraoperative radiotherapy, hadrontherapy, molecular radiotherapy, and in vivo dose monitoring. Investigations that have been performed using GEANT4 only are also mentioned to illustrate the potential of GATE. The very practical feature of GATE making it easy to model both a treatment and an imaging acquisition within the same framework is emphasized. The computational times associated with several applications are provided to illustrate the practical feasibility of the simulations using current computing facilities.

A virtual source model for Monte Carlo simulation of helical tomotherapy

The purpose of this study was to present a Monte Carlo (MC) simulation method based on a virtual source, jaw, and MLC model to calculate dose in patient for helical tomotherapy without the need of calculating phase‐space files (PSFs). Current studies on the tomotherapy MC simulation adopt a full MC model, which includes extensive modeling of radiation source, primary and secondary jaws, and multileaf collimator (MLC). In the full MC model, PSFs need to be created at different scoring planes to facilitate the patient dose calculations. In the present work, the virtual source model (VSM) we established was based on the gold standard beam data of a tomotherapy unit, which can be exported from the treatment planning station (TPS). The TPS‐generated sinograms were extracted from the archived patient XML (eXtensible Markup Language) files. The fluence map for the MC sampling was created by incorporating the percentage leaf open time (LOT) with leaf filter, jaw penumbra, and leaf latency contained from sinogram files. The VSM was validated for various geometry setups and clinical situations involving heterogeneous media and delivery quality assurance (DQA) cases. An agreement of <1% was obtained between the measured and simulated results for percent depth doses (PDDs) and open beam profiles for all three jaw settings in the VSM commissioning. The accuracy of the VSM leaf filter model was verified in comparing the measured and simulated results for a Picket Fence pattern. An agreement of <2% was achieved between the presented VSM and a published full MC model for heterogeneous phantoms. For complex clinical head and neck (HN) cases, the VSM‐based MC simulation of DQA plans agreed with the film measurement with 98% of planar dose pixels passing on the 2%/2 mm gamma criteria. For patient treatment plans, results showed comparable dose‐volume histograms (DVHs) for planning target volumes (PTVs) and organs at risk (OARs). Deviations observed in this study were consistent with literature. The VSM‐based MC simulation approach can be feasibly built from the gold standard beam model of a tomotherapy unit. The accuracy of the VSM was validated against measurements in homogeneous media, as well as published full MC model in heterogeneous media.

PACS numbers: 87.53.‐j, 87.55.K‐

Monte Carlo modeling of HD120 multileaf collimator on Varian TrueBeam linear accelerator for verification of 6X and 6X FFF VMAT SABR treatment plans

A Monte Carlo (MC) validation of the vendor-supplied Varian TrueBeam 6 MV flattened (6X) phase-space file and the first implementation of the Siebers-Keall MC MLC model as applied to the HD120 MLC (for 6X flat and 6X flattening filter-free (6X FFF) beams) are described. The MC model is validated in the context of VMAT patient-specific quality assurance. The Monte Carlo commissioning process involves: 1) validating the calculated open-field percentage depth doses (PDDs), profiles, and output factors (OF), 2) adapting the Siebers-Keall MLC model to match the new HD120-MLC geometry and material composition, 3) determining the absolute dose conversion factor for the MC calculation, and 4) validating this entire linac/MLC in the context of dose calculation verification for clinical VMAT plans. MC PDDs for the 6X beams agree with the measured data to within 2.0% for field sizes ranging from 2 × 2 to 40 × 40 cm2. Measured and MC profiles show agreement in the 50% field width and the 80%-20% penumbra region to within 1.3 mm for all square field sizes. MC OFs for the 2 to 40 cm2 square fields agree with measurement to within 1.6%. Verification of VMAT SABR lung, liver, and vertebra plans demonstrate that measured and MC ion chamber doses agree within 0.6% for the 6X beam and within 2.0% for the 6X FFF beam. A 3D gamma factor analysis demonstrates that for the 6X beam, &gt; 99% of voxels meet the pass criteria (3%/3 mm). For the 6X FFF beam, &gt; 94% of voxels meet this criteria. The TrueBeam accelerator delivering 6X and 6X FFF beams with the HD120 MLC can be modeled in Monte Carlo to provide an independent 3D dose calculation for clinical VMAT plans. This quality assurance tool has been used clinically to verify over 140 6X and 16 6X FFF TrueBeam treatment plans.

A bias-free, automated planning tool for technique comparison in radiotherapy - application to nasopharyngeal carcinoma treatments

In this study a novel, user‐independent automated planning technique was developed to objectively compare volumetric‐modulated arc therapy (VMAT) and intensity‐modulated radiotherapy (IMRT) for nasopharyngeal carcinoma planning, and to determine which technique offers a greater benefit for parotid‐sparing and dose escalation strategies. Ten patients were investigated, with a standard prescription of three dose levels to the target volumes (70, 63, and 56 Gy), using a simultaneous integrated boost in 33 fractions. The automated tool was used to investigate three planning strategies with both IMRT and VMAT: clinically acceptable plan creation, parotid dose sparing, and dose escalation. Clinically acceptable plans were achieved for all patients using both techniques. For parotid‐sparing, automated planning reduced the mean dose to a greater extent using VMAT rather than IMRT (17.0 Gy and 19.6 Gy, respectively, p<0.01). For dose escalation to the mean of the main clinical target volume, neither VMAT nor IMRT offered a significant benefit over the other. The OAR‐limiting prescriptions for VMAT ranged from 84‐98 Gy, compared to 76‐110 Gy for IMRT. Employing a user‐independent planning technique, it was possible to objectively compare VMAT and IMRT for nasopharyngeal carcinoma treatment strategies. VMAT offers a parotid‐sparing improvement, but no significant benefit was observed for dose escalation to the primary target.

PACS numbers: 87.55.D‐, 87.55.kd

Impact of the number of discrete angles used during dose computation for TomoTherapy treatments

PURPOSE

To quantify systematically the effect on accuracy of discretizing gantry rotation during the dose calculation process of TomoTherapy treatments.

METHODS

Up to version 4.0.x included, TomoTherapy treatment planning system (TPS) approximates gantry rotation by computing dose from 51 discrete angles corresponding to the center of the projections used to control the binary multileaf collimator. Potential effects on dose computation accuracy for off-axis targets and low modulation factors have been shown previously for a few treatment configurations. In versions 4.1.x and later, TomoTherapy oversamples the projections to better account for gantry rotation, but only during full scatter optimization and final calculation (i.e., not during optimization in "beamlet" mode). The effect on accuracy of changing the number of angles was quantified with the following framework: (1) predict the impact of the discretization of gantry rotation for various modulation factors, target sizes, and off-axis positions using a simplified analytical algorithm; (2) perform regular quality assurance using measurements with EDR2 radiographic films; (3) isolating the effect of changing the number of discretized angles only (51, 153, and 459) using a previously validated Monte Carlo model (TomoPen). The diameters of the targets were 2, 3, and 5 cm; off-axis central positions of target volumes were 5, 10 and 15, and 17 cm (when accepted by the treatment unit); planned modulation factors were 1.3 and 2.0.

RESULTS

For extreme configurations (3 cm tumor, 1.3 modulation factor, 15 cm off-axis position), effects on dose distributions were significant with 89.3% and 95.4% of the points passing gamma tests with 2%∕2 mm and 3%∕3 mm criteria, respectively, for TPS software version 4.0.x (51 gantry angles). The passing rate was 100% for both gamma criteria for the 4.1.x version (153 gantry angles). Those differences could be attributed almost completely to gantry motion discretization using TomoPen. Using 51 gantry angles for dose computation, TomoPen reproduced within statistical uncertainties (<1% standard deviation) dose distributions computed with version 4.0.x. Using 153 and 459 gantry angles, TomoPen reproduced within statistical uncertainties measurements and dose distributions computed with version 4.1.x.

CONCLUSIONS

When low modulation factors and significant off-axis positions are used, accounting for gantry rotation during dose computation using at least 153 gantry angles is required to ensure optimal accuracy.

A modular method to handle multiple time-dependent quantities in Monte Carlo simulations

A general method for handling time-dependent quantities in Monte Carlo simulations was developed to make such simulations more accessible to the medical community for a wide range of applications in radiotherapy, including fluence and dose calculation. To describe time-dependent changes in the most general way, we developed a grammar of functions that we call 'Time Features'. When a simulation quantity, such as the position of a geometrical object, an angle, a magnetic field, a current, etc, takes its value from a Time Feature, that quantity varies over time. The operation of time-dependent simulation was separated into distinct parts: the Sequence samples time values either sequentially at equal increments or randomly from a uniform distribution (allowing quantities to vary continuously in time), and then each time-dependent quantity is calculated according to its Time Feature. Due to this modular structure, time-dependent simulations, even in the presence of multiple time-dependent quantities, can be efficiently performed in a single simulation with any given time resolution. This approach has been implemented in TOPAS (TOol for PArticle Simulation), designed to make Monte Carlo simulations with Geant4 more accessible to both clinical and research physicists. To demonstrate the method, three clinical situations were simulated: a variable water column used to verify constancy of the Bragg peak of the Crocker Lab eye treatment facility of the University of California, the double-scattering treatment mode of the passive beam scattering system at Massachusetts General Hospital (MGH), where a spinning range modulator wheel accompanied by beam current modulation produces a spread-out Bragg peak, and the scanning mode at MGH, where time-dependent pulse shape, energy distribution and magnetic fields control Bragg peak positions. Results confirm the clinical applicability of the method.