Experimental verification of a Monte Carlo-based MLC simulation model for IMRT dose calculation

MLC parameters from static fields to VMAT plans: an evaluation in a RT-dedicated MC environment (PRIMO)

Background

PRIMO is a graphical environment based on PENELOPE Monte Carlo (MC) simulation of radiotherapy beams able to compute dose distribution in patients, from plans with different techniques. The dosimetric characteristics of an HD-120 MLC (Varian), simulated using PRIMO, were here compared with measurements, and also with Acuros calculations (in the Eclipse treatment planning system, Varian).

Materials and methods

A 10 MV FFF beam from a Varian EDGE linac equipped with the HD-120 MLC was used for this work. Initially, the linac head was simulated inside PRIMO, and validated against measurements in a water phantom. Then, a series of different MLC patterns were established to assess the MLC dosimetric characteristics. Those tests included: i) static fields: output factors from MLC shaped fields (2 × 2 to 10 × 10 cm²), alternate open and closed leaf pattern, MLC transmitted dose; ii) dynamic fields: dosimetric leaf gap (DLG) evaluated with sweeping gaps, tongue and groove (TG) effect assessed with profiles across alternate open and closed leaves moving across the field. The doses in the different tests were simulated in PRIMO and then compared with EBT3 film measurements in solid water phantom, as well as with Acuros calculations. Finally, MC in PRIMO and Acuros were compared in some clinical cases, summarizing the clinical complexity in view of a possible use of PRIMO as an independent dose calculation check.

Results

Static output factor MLC tests showed an agreement between MC calculated and measured OF of 0.5%. The dynamic tests presented DLG values of 0.033 ± 0.003 cm and 0.032 ± 0.006 cm for MC and measurements, respectively. Regarding the TG tests, a general agreement between the dose distributions of 1–2% was achieved, except for the extreme patterns (very small gaps/field sizes and high TG effect) were the agreement was about 4–5%. The analysis of the clinical cases, the Gamma agreement between MC in PRIMO and Acuros dose calculation in Eclipse was of 99.5 ± 0.2% for 3%/2 mm criteria of dose difference/distance to agreement.

Conclusions

MC simulations in the PRIMO environment were in agreement with measurements for the HD-120 MLC in a 10 MV FFF beam from a Varian EDGE linac. This result allowed to consistently compare clinical cases, showing the possible use of PRIMO as an independent dose calculation check tool.

Beam modeling and beam model commissioning for Monte Carlo dose calculation-based radiation therapy treatment planning: Report of AAPM Task Group 157

Dose calculation plays an important role in the accuracy of radiotherapy treatment planning and beam delivery. The Monte Carlo (MC) method is capable of achieving the highest accuracy in radiotherapy dose calculation and has been implemented in many commercial systems for radiotherapy treatment planning. The objective of this task group was to assist clinical physicists with the potentially complex task of acceptance testing and commissioning MC-based treatment planning systems (TPS) for photon and electron beam dose calculations. This report provides an overview on the general approach of clinical implementation and testing of MC-based TPS with a specific focus on models of clinical photon and electron beams. Different types of beam models are described including those that utilize MC simulation of the treatment head and those that rely on analytical methods and measurements. The trade-off between accuracy and efficiency in the various source-modeling approaches is discussed together with guidelines for acceptance testing of MC-based TPS from the clinical standpoint. Specific recommendations are given on methods and practical procedures to commission clinical beam models for MC-based TPS.

Development of a Geant4-based independent patient dose validation system with an elaborate multileaf collimator simulation model

Despite the improvements in the dose calculation models of the commercial treatment planning systems (TPS), their ability to accurately predict patient dose is still limited. One of the limitations is caused by the simplified model of the multileaf collimator (MLC). The aim of this study was to develop a Monte Carlo (MC) method‐based independent patient dose validation system with an elaborate MLC model for more accurate dose evaluation. Varian Clinac 2300 IX was simulated using Geant4 toolkits, after which MC commissioning with measurements was performed to validate the simulation model. A DICOM‐RT interface was developed to obtain the beam delivery conditions including the hundreds of MLC motions. Finally, the TPS dose distributions were compared with the MC dose distributions for water phantom cases and a patient case. Our results show that the TPS overestimated the absolute abutting leakage dose in the closed MLC field, with about 20% more of the maximum dose than that of the MC calculation. For water phantom cases, the dose distributions inside the target region were almost identical with the dose difference of less than 2%, while the dose near the edge of the target shows difference about 10% between Geant4 and TPS due to geometrical differences in MLC model. For the patient analysis, the Geant4 and TPS doses of all organs were matched well within 1.4% of the prescribed dose. However, for organs located in areas with high ratio of leaf pairs with distances less than 10 mm leaf pair (LP_(<10mm)), the maximum dose of TPS was overestimated by about 3% of the prescribed dose. These dose comparison results demonstrate that our system for calculating the patient dose is quite accurate. Furthermore, if the MLC sequences in treatment plan have a large ratio of LP_(short), more than 3% dose difference in normal tissue could be seen.

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.

Linking computer-aided design (CAD) to Geant4-based Monte Carlo simulations for precise implementation of complex treatment head geometries

Most of the treatment head components of medical linear accelerators used in radiation therapy have complex geometrical shapes. They are typically designed using computer-aided design (CAD) applications. In Monte Carlo simulations of radiotherapy beam transport through the treatment head components, the relevant beam-generating and beam-modifying devices are inserted in the simulation toolkit using geometrical approximations of these components. Depending on their complexity, such approximations may introduce errors that can be propagated throughout the simulation. This drawback can be minimized by exporting a more precise geometry of the linac components from CAD and importing it into the Monte Carlo simulation environment. We present a technique that links three-dimensional CAD drawings of the treatment head components to Geant4 Monte Carlo simulations of dose deposition.

An analytical approach for optimizing the leaf design of a multi-leaf collimator in a linear accelerator

In this study, we present an analytical approach for optimizing the leaf design of a multi-leaf collimator (MLC) in a linear accelerator. Because leaf designs vary between vendors, our goal is to characterize and quantify the effects of different compromises which have to be made between performance parameters. Subsequently, an optimal leaf design for an earlier proposed six-bank MLC which combines a high-resolution field-shaping ability with a large field size is determined. To this end a model of the linac is created that includes the following parameters: the source size, the maximum field size, the distance between source and isocenter, and the leaf's design parameters. First, the optimal radius of the leaf tip was found. This optimum was defined by the requirement that the fluence intensity should fall from 80% of the maximum value to 20% in a minimal distance, defining the width of the fluence penumbra. A second requirement was that this penumbra width should be constant when a leaf moves from one side of the field to the other. The geometric, transmission and total penumbra width (80-20%) were calculated depending on the design parameters. The analytical model is in agreement with Elekta, Varian and Siemens collimator designs. For leaves thinner than 4 cm, the transmission penumbra becomes dominant, and for leaves close to the source the geometric penumbra plays a role. Finally, by choosing the leaf thickness of 3.5 cm, 4 cm and 5 cm from the lowest to the highest bank, respectively, an optimal leaf design for a six-bank MLC is achieved.

Management of three-dimensional intrafraction motion through real-time DMLC tracking

Tumor tracking using a dynamic multileaf collimator (DMLC) represents a promising approach for intrafraction motion management in thoracic and abdominal cancer radiotherapy. In this work, we develop, empirically demonstrate, and characterize a novel 3D tracking algorithm for real-time, conformal, intensity modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT)-based radiation delivery to targets moving in three dimensions. The algorithm obtains real-time information of target location from an independent position monitoring system and dynamically calculates MLC leaf positions to account for changes in target position. Initial studies were performed to evaluate the geometric accuracy of DMLC tracking of 3D target motion. In addition, dosimetric studies were performed on a clinical linac to evaluate the impact of real-time DMLC tracking for conformal, step-and-shoot (S-IMRT), dynamic (D-IMRT), and VMAT deliveries to a moving target. The efficiency of conformal and IMRT delivery in the presence of tracking was determined. Results show that submillimeter geometric accuracy in all three dimensions is achievable with DMLC tracking. Significant dosimetric improvements were observed in the presence of tracking for conformal and IMRT deliveries to moving targets. A gamma index evaluation with a 3%-3 mm criterion showed that deliveries without DMLC tracking exhibit between 1.7 (S-IMRT) and 4.8 (D-IMRT) times more dose points that fail the evaluation compared to corresponding deliveries with tracking. The efficiency of IMRT delivery, as measured in the lab, was observed to be significantly lower in case of tracking target motion perpendicular to MLC leaf travel compared to motion parallel to leaf travel. Nevertheless, these early results indicate that accurate, real-time DMLC tracking of 3D tumor motion is feasible and can potentially result in significant geometric and dosimetric advantages leading to more effective management of intrafraction motion.

Spatial dependence of MLC transmission in IMRT delivery

In complex intensity-modulated radiation therapy cases, a considerable amount of the total dose may be delivered through closed leaves. In such cases an accurate knowledge of spatial characteristics of multileaf collimator (MLC) transmission is crucial, especially for the treatment of large targets with split fields. Measurements with an ionization chamber, radiographic films (EDR2, EBT) and EPID are taken to characterize all relevant effects related to MLC transmission for various field sizes and depths. Here we present a phenomenological model to describe MLC transmission, whereby the main focus is the off-axis decrease of transmission for symmetric and asymmetric fields as well as on effects due to the tongue and groove design of the leaves, such as interleaf transmission and the tongue and groove effect. Data obtained with the four different methods are presented, and the utility of each measurement method to determine the necessary model parameters is discussed. With the developed model, it is possible to predict the relevant MLC effects at any point in the phantom for arbitrary jaw settings and depths.