Assessment of a new multileaf collimator concept usingGEANT4Monte Carlo simulations

Experimental validation of a linac head Geant4 model under a grid computing environment

Background and purpose: This work aims to present the strategy to simulate a clinical linear accelerator based on the geometry provided by the manufacturer and summarize the corresponding experimental validation. Simulations were performed with the Geant4 Monte Carlo code under a grid computing environment. The objective of this contribution is reproducing therapeutic dose distributions in a water phantom with an accuracy less than 2%. Materials and methods: A Geant4 Monte Carlo model of an Elekta Synergy linear accelerator has been established, the simulations were launched in a large grid computing platform. Dose distributions are calculated for a 6 MV photon beam with treatment fields ranging from 5 × 5 cm2 to 20 × 20 cm 2 at a source - surface distance of 100 cm. Results: A high degree of agreement is achieved between the simulation results and the measured data, with dose differences of about 1.03% and 1.96% for the percentage depth dose curves and lateral dose profiles, respectively. This agreement is evaluated by the gamma index comparisons. Over 98% of the points for all simulations meet the restrictive acceptability criteria of 2%/2 mm. Conclusion: We have demonstrated the possibility to establish an accurate linac head Monte Carlo model for dose distribution simulations and quality assurance. Percentage depth dose curves and beam quality indices are in perfect agreement with the measured data with an accuracy of better than 2%.

CPU time optimization and precise adjustment of the Geant4 physics parameters for a VARIAN 2100 C/D gamma radiotherapy linear accelerator simulation using GAMOS

We have verified the GAMOS/Geant4 simulation model of a 6 MV VARIAN Clinac 2100 C/D linear accelerator by the procedure of adjusting the initial beam parameters to fit the percentage depth dose and cross-profile dose experimental data at different depths in a water phantom. Thanks to the use of a wide range of field sizes, from 2  ×  2 cm² to 40  ×  40 cm², a small phantom voxel size and high statistics, fine precision in the determination of the beam parameters has been achieved. This precision has allowed us to make a thorough study of the different physics models and parameters that Geant4 offers. The three Geant4 electromagnetic physics sets of models, i.e. Standard, Livermore and Penelope, have been compared to the experiment, testing the four different models of angular bremsstrahlung distributions as well as the three available multiple-scattering models, and optimizing the most relevant Geant4 electromagnetic physics parameters. Before the fitting, a comprehensive CPU time optimization has been done, using several of the Geant4 efficiency improvement techniques plus a few more developed in GAMOS.

Dosimetric impact assessment using a general algorithm in geant4 simulations for a complex-shaped multileaf collimator

PURPOSE

We have developed an inhouse algorithm for the multileaf collimator (MLC) geometry model construction with an appropriate accuracy for dosimetric tests. Our purpose is to build a complex type of MLC and analyze the influence of the modeling parameters on the dose calculation.

METHODS

Using radiochromic films as detector the following tests were done: (I) Density test field: to compare measured and calculated dose distributions in order to determine the tungsten alloy physical density value. (II) Leaf ends test field: to verify the penumbra shape sensitivity against the discretization level set to simulate the curved leaf ends. (III) MLC-closed field: to obtain the value of the air gap between opposite leaves for a closed configuration which completes the modeling of the MLC leakage radiation. (IV) Picket-fence field: to fit the leaf tilt angle with respect of the divergent ray emerging from the source.

RESULTS

For a 18.5g/cm³ density value we have obtained a maximum, minimum and mean leakage values of 0.43%, 0.36% and 0.38%, similar to the experimental ones. The best discretization level in the leaf ends field shows a 5.51mm FWHM, very close to the measured value (5.49mm). An air gap of 370μm has been used in the simulation for the separation between opposite leaves. Using a 0.44° tilt angle, we found the same pattern as the experimental values.

CONCLUSIONS

Our code can reproduce complex MLC designs with a submilimetric dosimetric accuracy which implies the necessary background for dose calculation of high clinical interest small fields.

Monte Carlo implementation, validation, and characterization of a 120 leaf MLC

PURPOSE

Recently, the new high definition multileaf collimator (HD120 MLC) was commercialized by Varian Medical Systems providing high resolution in the center section of the treatment field. The aim of this work is to investigate the characteristics of the HD120 MLC using Monte Carlo (MC) methods.

METHODS

Based on the information of the manufacturer, the HD120 MLC was implemented into the already existing Swiss MC Plan (SMCP). The implementation has been configured by adjusting the physical density and the air gap between adjacent leaves in order to match transmission profile measurements for 6 and 15 MV beams of a Novalis TX. These measurements have been performed in water using gafchromic films and an ionization chamber at an SSD of 95 cm and a depth of 5 cm. The implementation was validated by comparing diamond measured and calculated penumbra values (80%-20%) for different field sizes and water depths. Additionally, measured and calculated dose distributions for a head and neck IMRT case using the DELTA(4) phantom have been compared. The validated HD120 MLC implementation has been used for its physical characterization. For this purpose, phase space (PS) files have been generated below the fully closed multileaf collimator (MLC) of a 40 × 22 cm(2) field size for 6 and 15 MV. The PS files have been analyzed in terms of energy spectra, mean energy, fluence, and energy fluence in the direction perpendicular to the MLC leaves and have been compared with the corresponding data using the well established Varian 80 leaf (MLC80) and Millennium M120 (M120 MLC) MLCs. Additionally, the impact of the tongue and groove design of the MLCs on dose has been characterized.

RESULTS

Calculated transmission values for the HD120 MLC are 1.25% and 1.34% in the central part of the field for the 6 and 15 MV beam, respectively. The corresponding ionization chamber measurements result in a transmission of 1.20% and 1.35%. Good agreement has been found for the comparison between transmission profiles resulting from MC simulations and film measurements. The simulated and measured values for the penumbra agreed within <0.5 mm for all field sizes, depths, and beam energies, and a good agreement has been found between the measured and the calculated dose distributions for the IMRT case. The total energy spectra are almost identical for the three MLCs. However, the mean energy, fluence and energy fluence are significantly different. Due to the different leaf widths of the MLCs, the shape of these distributions is different, each representing its leave structure. Due to the increase in width from the inner to the outer HD120 MLC leaves, the fluence and energy fluence clearly decrease below the outer leaves. The MLC80 and the M120 MLC resulted in an increase of the fluence and energy fluence compared with those resulted for the HD120 MLC. The dose reduction can exceed 20% compared with the dose of the open field due to the tongue and groove design of the HD120 MLC.

CONCLUSIONS

The HD120 MLC has been successfully implemented into the SMCP. Comparisons between MC calculations and measurements show very good agreement. The SMCP is now able to calculate accurate dose distributions for treatment plans using the HD120 MLC.

Leakage of the Siemens 160 MLC multileaf collimator on a dual energy linear accelerator

Multileaf collimators (MLCs) have been in clinical use for many years and meanwhile are commonly used to deliver intensity-modulated radiotherapy (IMRT) beams. For this purpose it is important to know their dosimetric properties precisely, one of them being inter- and intraleaf leakage. The Siemens 160 MLC features a single focus design with flat-sided and tilted leaves instead of tongue-and-groove. The leakage performance of the 160 MLC was investigated on a dual energy linear accelerator Siemens ARTISTE with 6 MV and 18 MV photon energies. While the intraleaf leakage amounted to nearly the same dose for 6 and for 18 MV, a much higher interleaf leakage for 6 MV was measured. It could be reduced by simply rotating the collimator, and also by changing the voltage applied to the beam steering coils. The leakage of the 160 MLC is shown to be sensitive to beam alignment. This is of special interest for dual energy accelerators, as the two focal spots of both energies, neither in position nor in shape, do not necessarily always coincide. As a consequence of that, a higher leakage can be expected for one out of two energies for the 160 MLC.

Suggesting a new design for multileaf collimator leaves based on Monte Carlo simulation of two commercial systems

Due to intensive use of multileaf collimators (MLCs) in clinics, finding an optimum design for the leaves becomes essential. There are several studies which deal with comparison of MLC systems, but there is no article with a focus on offering an optimum design using accurate methods like Monte Carlo. In this study, we describe some characteristics of MLC systems including the leaf tip transmission, beam hardening, leakage radiation and penumbra width for Varian and Elekta 80-leaf MLCs using MCNP4C code. The complex geometry of leaves in these two common MLC systems was simulated. It was assumed that all of the MLC systems were mounted on a Varian accelerator and with a similar thickness as Varian's and the same distance from the source. Considering the obtained results from Varian and Elekta leaf designs, an optimum design was suggested combining the advantages of three common MLC systems and the simulation results of this proposed one were compared with the Varian and the Elekta. The leakage from suggested design is 29.7% and 31.5% of the Varian and Elekta MLCs. In addition, other calculated parameters of the proposed MLC leaf design were better than those two commercial ones. Although it shows a wider penumbra in comparison with Varian and Elekta MLCs, taking into account the curved motion path of the leaves, providing a double focusing design will solve the problem. The suggested leaf design is a combination of advantages from three common vendors (Varian, Elekta and Siemens) which can show better results than each one. Using the results of this theoretical study may bring about superior practical outcomes.

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.

6 MV dosimetric characterization of the 160 MLC, the new Siemens multileaf collimator

New technical developments constantly aim at improving the outcome of radiation therapy. With the use of a computer-controlled multileaf collimator (MLC), the quality of the treatment and the efficiency in patient throughput is significantly increased. New MLC designs aim to further enhance the advantages. In this article, we present the first detailed experimental investigation of the new 160 MLC, Siemens Medical Solutions. The assessment included the experimental investigation of typical MLC characteristics such as leakage, tongue-and-groove effect, penumbra, leaf speed, and leaf positioning accuracy with a 6 MV treatment beam. The leakage is remarkably low with an average of 0.37% due to a new design principle of slightly tilted leaves instead of the common tongue-and-groove design. But due to the tilt, the triangular tongue-and-groove effect occurs. Its magnitude of approximately 19% is similar to the dose defect measured for MLCs with the common tongue-and-groove design. The average longitudinal penumbra measured at depth d(max) = 15 mm with standard 100 x 100 mm2 fields is 4.1 +/- 0.5 mm for the central range and increases to 4.9 +/- 1.3 mm for the entire field range of 400 x 400 mm2. The increase is partly due to the single-focusing design and the large distance between the MLC and the isocenter enabling a large patient clearance. Regarding the leaf speed, different velocity tests were performed. The positions of the moving leaves were continuously recorded with the kilovoltage-imaging panel. The maximum leaf velocities measured were 42.9 +/- 0.6 mm/s. In addition, several typical intensity-modulated radiation therapy treatments were performed and the delivery times compared to the Siemens OPTIFOCUS MLC. An average decrease of 11% in delivery time was observed. The experimental results presented in this article indicate that the dosimetric characteristics of the 160 MLC are capable of improving the quality of dose delivery with respect to precision and dose conformity.

Siemens multileaf collimator characterization and quality assurance approaches for intensity-modulated radiotherapy

Application of the multileaf collimator (MLC) has evolved from replacing blocks to create treatment fields to creating photon fluence modulation for intensity-modulated radiotherapy (IMRT). Multileaf collimator system performance requirements are far more stringent for such applications and will require increased performance for future applications, such as motion tracking. This article reviews Siemens MLC systems, including a technical description and dosimetric characteristics of 56-, 82-, and 160-leaf designs. Routine quality assurance of MLC for IMRT necessitates frequent and critical assessment of MLC leaf position calibration errors that can present in many different ways (e.g., accuracy, reproducibility, longevity, hysteresis, and collimator/gantry angle dependencies). Several techniques for measuring these errors are presented, along with qualitative and quantitative techniques for analyzing results. In particular, increased accuracy of leaf position measurement at variable gantry angles is enabled by spatial transformations to electronic portal imaging device position quantified by calibration protocols introduced with megavoltage cone beam. Measured values of X-ray transmission (intra-leaf, inter-leaf, and through abutting leaf pairs) and penumbra (leaf end, leaf tongue, leaf groove) are presented with an evaluation of their characterization by a treatment-planning system. The dosimetric impact of planning system model inadequacies is demonstrated for collimator scatter, dose profile values within 30 mm of the field edge, and the resultant effect demonstrated on clinical cases. Finally, a description of automated quality assurance delivery, analysis, and calibration protocols applicable for the specific vendor's system is provided.

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.

Reducing stray radiation dose to patients receiving passively scattered proton radiotherapy for prostate cancer

Proton beam radiotherapy exposes healthy tissue to stray radiation emanating from the treatment unit and secondary radiation produced within the patient. These exposures provide no known benefit and may increase a patient's risk of developing a radiogenic second cancer. The aim of this study was to explore strategies to reduce stray radiation dose to a patient receiving a 76 Gy proton beam treatment for cancer of the prostate. The whole-body effective dose from stray radiation, E, was estimated using detailed Monte Carlo simulations of a passively scattered proton treatment unit and an anthropomorphic phantom. The predicted value of E was 567 mSv, of which 320 mSv was attributed to leakage from the treatment unit; the remainder arose from scattered radiation that originated within the patient. Modest modifications of the treatment unit reduced E by 212 mSv. Surprisingly, E from a modified passive-scattering device was only slightly higher (109 mSv) than from a nozzle with no leakage, e.g., that which may be approached with a spot-scanning technique. These results add to the body of evidence supporting the suitability of passively scattered proton beams for the treatment of prostate cancer, confirm that the effective dose from stray radiation was not excessive, and, importantly, show that it can be substantially reduced by modest enhancements to the treatment unit.