ORANGE: a Monte Carlo dose engine for radiotherapy

Monte Carlo simulation using PRIMO code as a tool for checking the credibility of commissioning and quality assurance of 6 MV TrueBeam STx varian LINAC

Aim

To validate and implement Monte Carlo simulation using PRIMO code as a tool for checking the credibility of measurements in LINAC initial commissioning and routine Quality Assurance (QA). Relative and absolute doses of 6 MV photon beam from TrueBeam STx Varian Linear Accelerator (LINAC) were simulated and validated with experimental measurement, Analytical Anisotropic Algorithm (AAA) calculation, and golden beam.

Methods and Materials

Varian phase-space files were imported to the PRIMO code and four blocks of jaws were simulated to determine the field size of the photon beam. Water phantom was modeled in the PRIMO code with water equivalent density. Golden beam data, experimental measurement, and AAA calculation results were imported to PRIMO code for gamma comparison.

Results

PRIMO simulations of Percentage Depth Dose (PDD) and in-plane beam profiles had good agreement with experimental measurements, AAA calculations and golden beam. However, PRIMO simulations of cross-plane beam profiles have a better agreement with AAA calculation and golden beam than the experimental measurement. Furthermore, PRIMO simulations of absolute dose agreed well with experimental results with ±0.8% uncertainty.

Conclusion

The PRIMO code has good accuracy and is appropriate for use as a tool to check the credibility of beam scanning and output measurement in initial commissioning and routine QA.

Neutron track length estimator for GATE Monte Carlo dose calculation in radiotherapy

The out-of-field dose in radiation therapy is a growing concern in regards to the late side-effects and secondary cancer induction. In high-energy x-ray therapy, the secondary neutrons generated through photonuclear reactions in the accelerator are part of this secondary dose. The neutron dose is currently not estimated by the treatment planning system while it appears to be preponderant for distances greater than 50 cm from the isocenter. Monte Carlo simulation has become the gold standard for accurately calculating the neutron dose under specific treatment conditions but the method is also known for having a slow statistical convergence, which makes it difficult to be used on a clinical basis. The neutron track length estimator, a neutron variance reduction technique inspired by the track length estimator method has thus been developped for the first time in the Monte Carlo code GATE to allow a fast computation of the neutron dose in radiotherapy. The details of its implementation, as well as the comparison of its performances against the analog MC method, are presented here. A gain of time from 15 to 400 can be obtained by our method, with a mean difference in the dose calculation of about 1% in comparison with the analog MC method.

PRIMO: a graphical environment for the Monte Carlo simulation of Varian and Elekta linacs

BACKGROUND

The accurate Monte Carlo simulation of a linac requires a detailed description of its geometry and the application of elaborate variance-reduction techniques for radiation transport. Both tasks entail a substantial coding effort and demand advanced knowledge of the intricacies of the Monte Carlo system being used.

METHODS

PRIMO, a new Monte Carlo system that allows the effortless simulation of most Varian and Elekta linacs, including their multileaf collimators and electron applicators, is introduced. PRIMO combines (1) accurate physics from the PENELOPE code, (2) dedicated variance-reduction techniques that significantly reduce the computation time, and (3) a user-friendly graphical interface with tools for the analysis of the generated data. PRIMO can tally dose distributions in phantoms and computerized tomographies, handle phase-space files in IAEA format, and import structures (planning target volumes, organs at risk) in the DICOM RT-STRUCT standard.

RESULTS

A prostate treatment, conformed with a high definition Millenium multileaf collimator (MLC 120HD) from a Varian Clinac 2100 C/D, is presented as an example. The computation of the dose distribution in 1.86×3.00×1.86 mm3 voxels with an average 2% standard statistical uncertainty, performed on an eight-core Intel Xeon at 2.67 GHz, took 1.8 h-excluding the patient-independent part of the linac, which required 3.8 h but it is simulated only once.

CONCLUSION

PRIMO is a self-contained user-friendly system that facilitates the Monte Carlo simulation of dose distributions produced by most currently available linacs. This opens the door for routine use of Monte Carlo in clinical research and quality assurance purposes. It is free software that can be downloaded from http://www.primoproject.net.

An efficient numerical tool for dose deposition prediction applied to synchrotron medical imaging and radiation therapy

Medical imaging and radiation therapy are widely used synchrotron-based techniques which have one thing in common: a significant dose delivery to typically biological samples. Among the ways to provide the experimenters with image guidance techniques indicating optimization strategies, Monte Carlo simulation has become the gold standard for accurately predicting radiation dose levels under specific irradiation conditions. A highly important hampering factor of this method is, however, its slow statistical convergence. A track length estimator (TLE) module has been coded and implemented for the first time in the open-source Monte Carlo code GATE/Geant4. Results obtained with the module and the procedures used to validate them are presented. A database of energy-absorption coefficients was also generated, which is used by the TLE calculations and is now also included in GATE/Geant4. The validation was carried out by comparing the TLE-simulated doses with experimental data in a synchrotron radiation computed tomography experiment. The TLE technique shows good agreement versus both experimental measurements and the results of a classical Monte Carlo simulation. Compared with the latter, it is possible to reach a pre-defined statistical uncertainty in about two to three orders of magnitude less time for complex geometries without loss of accuracy.

Verification measurements and clinical evaluation of the iPlan RT Monte Carlo dose algorithm for 6 MV photon energy

This study presents data for verification of the iPlan RT Monte Carlo (MC) dose algorithm (BrainLAB, Feldkirchen, Germany). MC calculations were compared with pencil beam (PB) calculations and verification measurements in phantoms with lung-equivalent material, air cavities or bone-equivalent material to mimic head and neck and thorax and in an Alderson anthropomorphic phantom. Dosimetric accuracy of MC for the micro-multileaf collimator (MLC) simulation was tested in a homogeneous phantom. All measurements were performed using an ionization chamber and Kodak EDR2 films with Novalis 6 MV photon beams. Dose distributions measured with film and calculated with MC in the homogeneous phantom are in excellent agreement for oval, C and squiggle-shaped fields and for a clinical IMRT plan. For a field with completely closed MLC, MC is much closer to the experimental result than the PB calculations. For fields larger than the dimensions of the inhomogeneities the MC calculations show excellent agreement (within 3%/1 mm) with the experimental data. MC calculations in the anthropomorphic phantom show good agreement with measurements for conformal beam plans and reasonable agreement for dynamic conformal arc and IMRT plans. For 6 head and neck and 15 lung patients a comparison of the MC plan with the PB plan was performed. Our results demonstrate that MC is able to accurately predict the dose in the presence of inhomogeneities typical for head and neck and thorax regions with reasonable calculation times (5-20 min). Lateral electron transport was well reproduced in MC calculations. We are planning to implement MC calculations for head and neck and lung cancer patients.

Retrospective monte carlo dose calculations with limited beam weight information

An important unresolved issue in outcomes analysis for lung complications is the effect of poor or completely lacking heterogeneity corrections in previously archived treatment plans. To estimate this effect, we developed a novel method based on Monte Carlo (MC) dose calculations which can be applied retrospectively to RTOG/AAPM-style archived treatment plans (ATP). We applied this method to 218 archived nonsmall cell lung cancer lung treatment plans that were originally calculated either without heterogeneity corrections or with primitive corrections. To retrospectively specify beam weights and wedges, beams were broken into Monte Carlo-generated beamlets, simulated using the VMC++ code, and mathematical optimization was used to match the archived water-based dose distributions. The derived beam weights (and any wedge effects) were then applied to Monte Carlo beamlets regenerated based on the patient computed tomography densities. Validation of the process was performed against five comparable lung treatment plans generated using a commercial convolution/superposition implementation. For the application here (normal lung, esophagus, and planning target volume dose distributions), the agreement was very good. Resulting MC and convolution/superposition values were similar when dose distributions without heterogeneity corrections or dose distributions with corrections were compared. When applied to the archived plans (218), the average absolute percent difference between water-based MC and water-based ATPs, for doses above 2.5% of the maximum dose was 1.8+/-0.6%. The average absolute percent difference between heterogeneity-corrected MC and water-based ATPs increased to 3.1+/-0.9%. The average absolute percent difference between the MC heterogeneity-corrected and the ATP heterogeneity-corrected dose distributions was 3.8+/-1.6% (available in 132/218 archives). The entire dose-volume-histograms for lung, tumor, and esophagus from the different calculation methods, as well as specific dose metrics, were compared. The average difference in maximum lung dose between water-based ATPs and heterogeneity-corrected MC dose distributions was -1.0+/-2.1 Gy. Potential errors in relying on primitive heterogeneity corrections are most evident from a comparison of maximum lung doses, for which the average MC heterogeneity-corrected values were 5.3+/-2.8 Gy less than the ATP heterogeneity-corrected values. We have demonstrated that recalculation of archived dose distributions, without explicit information about beam weights or wedges, is feasible using beamlet-based optimization methods. The method provides heterogeneity-corrected dose data consistent with convolution-superposition calculations and is one feasible approach for improving dosimetric data for outcomes analyses.

Treating voxel geometries in radiation protection dosimetry with a patched version of the Monte Carlo codes MCNP and MCNPX

The question of Monte Carlo simulation of radiation transport in voxel geometries is addressed. Patched versions of the MCNP and MCNPX codes are developed aimed at transporting radiation both in the standard geometry mode and in the voxel geometry treatment. The patched code reads an unformatted FORTRAN file derived from DICOM format data and uses special subroutines to handle voxel-to-voxel radiation transport. The various phases of the development of the methodology are discussed together with the new input options. Examples are given of employment of the code in internal and external dosimetry and comparisons with results from other groups are reported.

Fast dose calculation for stereotactic synchrotron radiotherapy

A hybrid approach is proposed to compute the dose deposited in cancerous and healthy tissues during stereotactic synchrotron radiotherapy treatment. In this approach the computation is divided into two parts: (1) the primary dose is calculated using a deterministic algorithm based on ray casting; (2) the secondary dose (due to scattering and fluorescence) is computed using a hybrid algorithm combining Monte Carlo and a deterministic method. The results obtained for test cases are compared to those obtained with the Monte Carlo method alone (Geant4 code) and found to be in excellent agreement. The proposed simulation scheme makes it possible to simulate dose maps with a single personal computer, with computation time and statistical fluctuations substantially reduced in comparison with conventional Monte Carlo simulations.

Review of electron beam therapy physics

For over 50 years, electron beams have been an important modality for providing an accurate dose of radiation to superficial cancers and disease and for limiting the dose to underlying normal tissues and structures. This review looks at many of the important contributions of physics and dosimetry to the development and utilization of electron beam therapy, including electron treatment machines, dose specification and calibration, dose measurement, electron transport calculations, treatment and treatment-planning tools, and clinical utilization, including special procedures. Also, future changes in the practice of electron therapy resulting from challenges to its utilization and from potential future technology are discussed.