skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: TU-AB-BRC-04: Commissioning of a New MLC Model for the GEPTS Monte Carlo System: A Model Based On the Leaf and Interleaf Effective Density

Abstract

Purpose: To commission a new MLC model for the GEPTS Monte Carlo system. The model is based on the concept of leaves and interleaves effective densities Methods: GEPTS is a Monte Carlo system to be used for external beam planning verification. GEPTS incorporates detailed photon and electron transport algorithms (Med.Phys. 29, 2002, 835). A new GEPTS model for the Varian Millennium MLC is presented. The model accounts for: 1) thick (1 cm) and thin (0.5 cm) leaves, 2) tongue-and-groove design, 3) High-Transmission (HT) and Low-Transmission (LT) interleaves, and 4) rounded leaf end. Leaf (and interleaf) height is set equal to 6 cm. Instead of modeling air gaps, screw holes, and complex leaf heads, “effective densities” are assigned to: 1) thin leaves, 2) thick leaves, 3) HT-, and 4) LT-interleaves. Results: The new MLC model is used to calculate dose profiles for Closed-MLC and Tongue-and-Groove fields at 5 cm depth for 6, 10 and 15 MV Varian beams. Calculations are compared with 1) Pin-point ionization chamber transmission ratios and 2) EBT3 Radiochromic films. Pinpoint readings were acquired beneath thick and thin leaves, and HT and LT interleaves. The best fit of measured dose profiles was obtained for the following parameters: Thick-leafmore » density = 16.1 g/cc, Thin-leaf density = 17.2 g/cc; HT Interleaf density = 12.4 g/cc, LT Interleaf density = 14.3 g/cc; Interleaf thickness = 1.1 mm. Attached figures show comparison of calculated and measured transmission ratios for the 3 energies. Note this is the only study where transmission profiles are compared with measurements for 3 different energies. Conclusion: The new MLC model reproduces transmission measurements within 0.1%. The next step is to implement the MLC model for real plans and quantify the improvement in dose calculation accuracy gained using this model for IMRT plans with high modulation factors.« less

Authors:
; ;  [1]
  1. Fox Chase Cancer Center, Philadelphia, PA (United States)
Publication Date:
OSTI Identifier:
22653933
Resource Type:
Journal Article
Resource Relation:
Journal Name: Medical Physics; Journal Volume: 43; Journal Issue: 6; Other Information: (c) 2016 American Association of Physicists in Medicine; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
60 APPLIED LIFE SCIENCES; 61 RADIATION PROTECTION AND DOSIMETRY; COMMISSIONING; HIGH-TC SUPERCONDUCTORS; IONIZATION CHAMBERS; LEAVES; MONTE CARLO METHOD; PLANNING; RADIATION DOSES; SIMULATION; THICKNESS; TRANSMISSION

Citation Formats

Chibani, O, Tahanout, F, and Ma, C. TU-AB-BRC-04: Commissioning of a New MLC Model for the GEPTS Monte Carlo System: A Model Based On the Leaf and Interleaf Effective Density. United States: N. p., 2016. Web. doi:10.1118/1.4957398.
Chibani, O, Tahanout, F, & Ma, C. TU-AB-BRC-04: Commissioning of a New MLC Model for the GEPTS Monte Carlo System: A Model Based On the Leaf and Interleaf Effective Density. United States. doi:10.1118/1.4957398.
Chibani, O, Tahanout, F, and Ma, C. Wed . "TU-AB-BRC-04: Commissioning of a New MLC Model for the GEPTS Monte Carlo System: A Model Based On the Leaf and Interleaf Effective Density". United States. doi:10.1118/1.4957398.
@article{osti_22653933,
title = {TU-AB-BRC-04: Commissioning of a New MLC Model for the GEPTS Monte Carlo System: A Model Based On the Leaf and Interleaf Effective Density},
author = {Chibani, O and Tahanout, F and Ma, C},
abstractNote = {Purpose: To commission a new MLC model for the GEPTS Monte Carlo system. The model is based on the concept of leaves and interleaves effective densities Methods: GEPTS is a Monte Carlo system to be used for external beam planning verification. GEPTS incorporates detailed photon and electron transport algorithms (Med.Phys. 29, 2002, 835). A new GEPTS model for the Varian Millennium MLC is presented. The model accounts for: 1) thick (1 cm) and thin (0.5 cm) leaves, 2) tongue-and-groove design, 3) High-Transmission (HT) and Low-Transmission (LT) interleaves, and 4) rounded leaf end. Leaf (and interleaf) height is set equal to 6 cm. Instead of modeling air gaps, screw holes, and complex leaf heads, “effective densities” are assigned to: 1) thin leaves, 2) thick leaves, 3) HT-, and 4) LT-interleaves. Results: The new MLC model is used to calculate dose profiles for Closed-MLC and Tongue-and-Groove fields at 5 cm depth for 6, 10 and 15 MV Varian beams. Calculations are compared with 1) Pin-point ionization chamber transmission ratios and 2) EBT3 Radiochromic films. Pinpoint readings were acquired beneath thick and thin leaves, and HT and LT interleaves. The best fit of measured dose profiles was obtained for the following parameters: Thick-leaf density = 16.1 g/cc, Thin-leaf density = 17.2 g/cc; HT Interleaf density = 12.4 g/cc, LT Interleaf density = 14.3 g/cc; Interleaf thickness = 1.1 mm. Attached figures show comparison of calculated and measured transmission ratios for the 3 energies. Note this is the only study where transmission profiles are compared with measurements for 3 different energies. Conclusion: The new MLC model reproduces transmission measurements within 0.1%. The next step is to implement the MLC model for real plans and quantify the improvement in dose calculation accuracy gained using this model for IMRT plans with high modulation factors.},
doi = {10.1118/1.4957398},
journal = {Medical Physics},
number = 6,
volume = 43,
place = {United States},
year = {Wed Jun 15 00:00:00 EDT 2016},
month = {Wed Jun 15 00:00:00 EDT 2016}
}
  • Purpose: We have previously developed a GPU-OpenCL-based MC dose engine named goMC with built-in analytical linac beam model. To move goMC towards routine clinical use, we have developed an automatic beam-commissioning method, and an efficient source sampling strategy to facilitate dose calculations for real treatment plans. Methods: Our commissioning method is to automatically adjust the relative weights among the sub-sources, through an optimization process minimizing the discrepancies between calculated dose and measurements. Six models built for Varian Truebeam linac photon beams (6MV, 10MV, 15MV, 18MV, 6MVFFF, 10MVFFF) were commissioned using measurement data acquired at our institution. To facilitate dose calculationsmore » for real treatment plans, we employed inverse sampling method to efficiently incorporate MLC leaf-sequencing into source sampling. Specifically, instead of sampling source particles control-point by control-point and rejecting the particles blocked by MLC, we assigned a control-point index to each sampled source particle, according to MLC leaf-open duration of each control-point at the pixel where the particle intersects the iso-center plane. Results: Our auto-commissioning method decreased distance-to-agreement (DTA) of depth dose at build-up regions by 36.2% averagely, making it within 1mm. Lateral profiles were better matched for all beams, with biggest improvement found at 15MV for which root-mean-square difference was reduced from 1.44% to 0.50%. Maximum differences of output factors were reduced to less than 0.7% for all beams, with largest decrease being from1.70% to 0.37% found at 10FFF. Our new sampling strategy was tested on a Head&Neck VMAT patient case. Achieving clinically acceptable accuracy, the new strategy could reduce the required history number by a factor of ∼2.8 given a statistical uncertainty level and hence achieve a similar speed-up factor. Conclusion: Our studies have demonstrated the feasibility and effectiveness of our auto-commissioning approach and new efficient source sampling strategy, implying the potential of our GPU-based MC dose engine goMC for routine clinical use.« less
  • Purpose: One of the most accurate methods for radiation transport is Monte Carlo (MC) simulation. Long computation time prevents its wide applications in clinic. We have recently developed a fast MC code for carbon ion therapy called GPU-based OpenCL Carbon Monte Carlo (goCMC) and its accuracy in physical dose has been established. Since radiobiology is an indispensible aspect of carbon ion therapy, this study evaluates accuracy of goCMC in biological dose and microdosimetry by benchmarking it with FLUKA. Methods: We performed simulations of a carbon pencil beam with 150, 300 and 450 MeV/u in a homogeneous water phantom using goCMCmore » and FLUKA. Dose and energy spectra for primary and secondary ions on the central beam axis were recorded. Repair-misrepair-fixation model was employed to calculate Relative Biological Effectiveness (RBE). Monte Carlo Damage Simulation (MCDS) tool was used to calculate microdosimetry parameters. Results: Physical dose differences on the central axis were <1.6% of the maximum value. Before the Bragg peak, differences in RBE and RBE-weighted dose were <2% and <1%. At the Bragg peak, the differences were 12.5% caused by small range discrepancy and sensitivity of RBE to beam spectra. Consequently, RBE-weighted dose difference was 11%. Beyond the peak, RBE differences were <20% and primarily caused by differences in the Helium-4 spectrum. However, the RBE-weighted dose agreed within 1% due to the low physical dose. Differences in microdosimetric quantities were small except at the Bragg peak. The simulation time per source particle with FLUKA was 0.08 sec, while goCMC was approximately 1000 times faster. Conclusion: Physical doses computed by FLUKA and goCMC were in good agreement. Although relatively large RBE differences were observed at and beyond the Bragg peak, the RBE-weighted dose differences were considered to be acceptable.« less
  • Purpose: To create a Varian TrueBeam 6 MV FFF Monte Carlo model using BEAMnrc/EGSnrc that accurately reproduces the Varian representative dataset, followed by tuning the model’s source parameters to accurately reproduce in-house measurements. Methods: A BEAMnrc TrueBeam model for 6 MV FFF has been created by modifying a validated 6 MV Varian CL21EX model. Geometric dimensions and materials were adjusted in a trial and error approach to match the fluence and spectra of TrueBeam phase spaces output by the Varian VirtuaLinac. Once the model’s phase space matched Varian’s counterpart using the default source parameters, it was validated to match 10more » × 10 cm{sup 2} Varian representative data obtained with the IBA CC13. The source parameters were then tuned to match in-house 5 × 5 cm{sup 2} PTW microDiamond measurements. All dose to water simulations included detector models to include the effects of volume averaging and the non-water equivalence of the chamber materials, allowing for more accurate source parameter selection. Results: The Varian phase space spectra and fluence were matched with excellent agreement. The in-house model’s PDD agreement with CC13 TrueBeam representative data was within 0.9% local percent difference beyond the first 3 mm. Profile agreement at 10 cm depth was within 0.9% local percent difference and 1.3 mm distance-to-agreement in the central axis and penumbra regions, respectively. Once the source parameters were tuned, PDD agreement with microDiamond measurements was within 0.9% local percent difference beyond 2 mm. The microDiamond profile agreement at 10 cm depth was within 0.6% local percent difference and 0.4 mm distance-to-agreement in the central axis and penumbra regions, respectively. Conclusion: An accurate in-house Monte Carlo model of the Varian TrueBeam was achieved independently of the Varian phase space solution and was tuned to in-house measurements. KO acknowledges partial support by the CREATE Medical Physics Research Training Network grant of the Natural Sciences and Engineering Research Council (Grant number: 432290).« less
  • Purpose: To develop a general method for human tissue characterization with dual-and multi-energy CT and evaluate its performance in determining elemental compositions and the associated proton stopping power relative to water (SPR) and photon mass absorption coefficients (EAC). Methods: Principal component analysis is used to extract an optimal basis of virtual materials from a reference dataset of tissues. These principal components (PC) are used to perform two-material decomposition using simulated DECT data. The elemental mass fraction and the electron density in each tissue is retrieved by measuring the fraction of each PC. A stoichiometric calibration method is adapted to themore » technique to make it suitable for clinical use. The present approach is compared with two others: parametrization and three-material decomposition using the water-lipid-protein (WLP) triplet. Results: Monte Carlo simulations using TOPAS for four reference tissues shows that characterizing them with only two PC is enough to get a submillimetric precision on proton range prediction. Based on the simulated DECT data of 43 references tissues, the proposed method is in agreement with theoretical values of protons SPR and low-kV EAC with a RMS error of 0.11% and 0.35%, respectively. In comparison, parametrization and WLP respectively yield RMS errors of 0.13% and 0.29% on SPR, and 2.72% and 2.19% on EAC. Furthermore, the proposed approach shows potential applications for spectral CT. Using five PC and five energy bins reduces the SPR RMS error to 0.03%. Conclusion: The proposed method shows good performance in determining elemental compositions from DECT data and physical quantities relevant to radiotherapy dose calculation and generally shows better accuracy and unbiased results compared to reference methods. The proposed method is particularly suitable for Monte Carlo calculations and shows promise in using more than two energies to characterize human tissue with CT.« less
  • Purpose: To introduce egs-brachy, a new, fast, and versatile Monte Carlo code for brachytherapy applications. Methods: egs-brachy is an EGSnrc user-code based on the EGSnrc C++ class library (egs++). Complex phantom, applicator, and source model geometries are built using the egs++ geometry module. egs-brachy uses a tracklength estimator to score collision kerma in voxels. Interaction, spectrum, energy fluence, and phase space scoring are also implemented. Phase space sources and particle recycling may be used to improve simulation efficiency. HDR treatments (e.g. stepping source through dwell positions) can be simulated. Standard brachytherapy seeds, as well as electron and miniature x-ray tubemore » sources are fully modelled. Variance reduction techniques for electron source simulations are implemented (Bremsstrahlung cross section enhancement, uniform Bremsstrahlung splitting, and Russian Roulette). TG-43 parameters of seeds are computed and compared to published values. Example simulations of various treatments are carried out on a single 2.5 GHz Intel Xeon E5-2680 v3 processor core. Results: TG-43 parameters calculated with egs-brachy show excellent agreement with published values. Using a phase space source, 2% average statistical uncertainty in the PTV ((2mm){sup 3} voxels) can be achieved in 10 s for 100 {sup 125}I or {sup 103}Pd seeds in a 36.2 cm{sup 3} prostate PTV, 31 s for 64 {sup 103}Pd seeds in a 64 cm{sup 3} breast PTV, and 56 s for a miniature x-ray tube in a 27 cm{sup 3} breast PTV. Comparable uncertainty is reached in 12 s in a (1 mm){sup 3} water voxel 5 mm away from a COMS 16mm eye plaque with 13 {sup 103}Pd seeds. Conclusion: The accuracy of egs-brachy has been demonstrated through benchmarking calculations. Calculation times are sufficiently fast to allow full MC simulations for routine treatment planning for diverse brachytherapy treatments (LDR, HDR, miniature x-ray tube). egs-brachy will be available as free and open-source software to the medical physics research community. This work is partially funded by the Canada Research Chairs program, the Natural Sciences and Engineering Research Council of Canada, and the Ontario Ministry of Research and Innovation (Ontario Early Researcher Award).« less