MCNP-based computational model for the Leksell Gamma Knife

Monte-Carlo techniques for radiotherapy applications I: introduction and overview of the different Monte-Carlo codes

Introduction:

The dose calculation plays a crucial role in many aspects of contemporary clinical radiotherapy treatment planning process. It therefore goes without saying that the accuracy of the dose calculation is of very high importance. The gold standard for absorbed dose calculation is the Monte-Carlo algorithm.

Methods:

This first of two papers gives an overview of the main openly available and supported codes that have been widely used for radiotherapy simulations.

Results:

The paper aims to provide an overview of Monte-Carlo in the field of radiotherapy and point the reader in the right direction of work that could help them get started or develop their existing understanding and use of Monte-Carlo algorithms in their practice.

Conclusions:

It also serves as a useful companion to a curated collection of papers on Monte-Carlo that have been published in this journal.

Monte Carlo methods for device simulations in radiation therapy

Monte Carlo (MC) simulations play an important role in radiotherapy, especially as a method to evaluate physical properties that are either impossible or difficult to measure. For example, MC simulations (MCSs) are used to aid in the design of radiotherapy devices or to understand their properties. The aim of this article is to review the MC method for device simulations in radiation therapy. After a brief history of the MC method and popular codes in medical physics, we review applications of the MC method to model treatment heads for neutral and charged particle radiation therapy as well as specific in-room devices for imaging and therapy purposes. We conclude by discussing the impact that MCSs had in this field and the role of MC in future device design.

Comparison of penumbra regions produced by ancient Gamma knife model C and Gamma ART 6000 using Monte Carlo MCNP6 simulation

The accuracy of penumbral measurements in radiotherapy is pivotal because dose planning computers require accurate data to adequately modeling the beams, which in turn are used to calculate patient dose distributions. Gamma knife is a non-invasive intracranial technique based on principles of the Leksell stereotactic system for open deep brain surgeries, invented and developed by Professor Lars Leksell. The aim of this study is to compare the penumbra widths of Leksell Gamma Knife model C and Gamma ART 6000. Initially, the structure of both systems were simulated by using Monte Carlo MCNP6 code and after validating the accuracy of simulation, beam profiles of different collimators were plotted. MCNP6 beam profile calculations showed that the penumbra values of Leksell Gamma knife model C and Gamma ART 6000 for 18, 14, 8 and 4 mm collimators are 9.7, 7.9, 4.3, 2.6 and 8.2, 6.9, 3.6, 2.4, respectively. The results of this study showed that since Gamma ART 6000 has larger solid angle in comparison with Gamma Knife model C, it produces better beam profile penumbras than Gamma Knife model C in the direct plane.

On the development of a comprehensive MC simulation model for the Gamma Knife Perfexion radiosurgery unit

This work presents a comprehensive Monte Carlo (MC) simulation model for the Gamma Knife Perfexion (PFX) radiosurgery unit. Model-based dosimetry calculations were benchmarked in terms of relative dose profiles (RDPs) and output factors (OFs), against corresponding EBT2 measurements. To reduce the rather prolonged computational time associated with the comprehensive PFX model MC simulations, two approximations were explored and evaluated on the grounds of dosimetric accuracy. The first consists in directional biasing of the (60)Co photon emission while the second refers to the implementation of simplified source geometric models. The effect of the dose scoring volume dimensions in OF calculations accuracy was also explored. RDP calculations for the comprehensive PFX model were found to be in agreement with corresponding EBT2 measurements. Output factors of 0.819  ±  0.004 and 0.8941  ±  0.0013 were calculated for the 4 mm and 8 mm collimator, respectively, which agree, within uncertainties, with corresponding EBT2 measurements and published experimental data. Volume averaging was found to affect OF results by more than 0.3% for scoring volume radii greater than 0.5 mm and 1.4 mm for the 4 mm and 8 mm collimators, respectively. Directional biasing of photon emission resulted in a time efficiency gain factor of up to 210 with respect to the isotropic photon emission. Although no considerable effect on relative dose profiles was detected, directional biasing led to OF overestimations which were more pronounced for the 4 mm collimator and increased with decreasing emission cone half-angle, reaching up to 6% for a 5° angle. Implementation of simplified source models revealed that omitting the sources' stainless steel capsule significantly affects both OF results and relative dose profiles, while the aluminum-based bushing did not exhibit considerable dosimetric effect. In conclusion, the results of this work suggest that any PFX simulation model should be benchmarked in terms of both RDP and OF results.

A simple method to prolong the service life of radioactive sources for external radiotherapy

A radioactive source is usually replaced and disposed after being used for a certain amount of time (usually a half‐life). In this study, a simple method is proposed to prolong its service life. Instead of replacing the used source with a new source of full activity, a new source of less activity is added in the source holder in front of the used one, so that the total activity of two sources is equal to the initial activity of the used source or even higher. Similarly, more sources can be added to the previous ones. Attenuation of front source(s) to the back source(s) was evaluated with exponential attenuation equation, and variation of source‐focus distance (SFD) with inverse square law for Leksell 4C Gamma Knife, which served as an example of external radiotherapy units. When the number of front sources increased from 1 to 3, the relative air kerma decreased from 36.5% to 5.0%. Both the attenuation effect and SFD variation contributed to the decrease in air kerma, with the former being the major factor. If the height of the source can be decreased in some way, such as increasing the specific activity of sources, the sources can be used more efficiently. The method prolongs the service life of sources by several factors, and reduces the expense of source exchange and reclamation.

PACS number: 87.56.bg

MAGAT gel and EBT2 film-based dosimetry for evaluating source plugging-based treatment plan in Gamma Knife stereotactic radiosurgery

This work illustrates a procedure to assess the overall accuracy associated with Gamma Knife treatment planning using plugging. The main role of source plugging or blocking is to create dose falloff in the junction between a target and a critical structure. We report the use of MAGAT gel dosimeter for verification of an experimental treatment plan based on plugging. The polymer gel contained in a head‐sized glass container simulated all major aspects of the treatment process of Gamma Knife radiosurgery. The 3D dose distribution recorded in the gel dosimeter was read using a 1.5T MRI scanner. Scanning protocol was: CPMG pulse sequence with 8 equidistant echoes, TR=7 s, echo step=14 ms, pixel size=0.5 mm x 0.5 mm, and slice thickness of 2 mm. Using a calibration relationship between absorbed dose and spin‐spin relaxation rate (R2), we converted R2 images to dose images. Volumetric dose comparison between treatment planning system (TPS) and gel measurement was accomplished using an in‐house MATLAB‐based program. The isodose overlay of the measured and computed dose distribution on axial planes was in close agreement. Gamma index analysis of 3D data showed more than 94% voxel pass rate for different tolerance criteria of 3%/2 mm, 3%/1 mm and 2%/2 mm. Film dosimetry with GAFCHROMIC EBT 2 film was also performed to compare the results with the calculated TPS dose. Gamma index analysis of film measurement for the same tolerance criteria used for gel measurement evaluation showed more than 95% voxel pass rate. Verification of gamma plan calculated dose on account of shield is not part of acceptance testing of Leksell Gamma Knife (LGK). Through this study we accomplished a volumetric comparison of dose distributions measured with a polymer gel dosimeter and Leksell GammaPlan (LGP) calculations for plans using plugging. We propose gel dosimeter as a quality assurance (QA) tool for verification of plug‐based planning.

PACS number: 87.53.Ly, 87.55.‐x, 87.56.N‐

FLUKA Monte Carlo simulation for the Leksell Gamma Knife Perfexion radiosurgery system: homogeneous media

The purpose of this work is to investigate the capability of the FLUKA Monte Carlo (MC) code to simulate the Elekta Leksell Gamma Knife Perfexion (LGK-PFX) and reproduce the Treatment Planning System (TPS) Leksell GammaPlan version 8.2 (LGP) dose calculations for the case of a water equivalent phantom target. Thanks to the collaboration with Elekta Instruments AB, the collimation system geometry, the source positions and all the involved material have been simulated in detail. The relative linear dose distribution along the three coordinate axes, for each collimator size, and the Relative Output Factors (ROF) have been investigated. The simulation has been validated comparing simulated linear dose profiles with measurements performed with EBT radiochromic films. The acceptance criterion between experimental data and FLUKA results is based on the gamma index (GI) method. The FLUKA MC calculation for the ROF provided the values of 0.920 for the 8 mm collimators and 0.800 for the 4 mm collimators. These values are in good agreement with the Elekta reference data of 0.924 and 0.805 respectively. The percentage difference between calculated and reference values for the ROF is under 1% and within the FLUKA uncertainty. Also the simulated relative dose profiles show a good agreement with the LGP calculation expressed by means of the gamma index method. This established accuracy proves that FLUKA is a suitable and powerful tool in order to reproduce successfully the LGP calculations for the homogeneous media.

Study of a spherical phantom for Gamma knife dosimetry

Four 16 cm diameter spherical phantoms were modeled in this study: a homogenous water phantom, and three water phantoms with 1 cm thick shell each made of different materials (PMMA, Plastic Water™ and polystyrene). The PENELOPE Monte Carlo code was utilized to simulate photon beams from the Leksell Gamma Knife (LGK) unit and to determine absorbed dose to water (Dw) from a single 18 mm beam delivered to each phantom. A score spherical volume of 0.007 cm³ was used to simulate the dimensions of the sensitive volume of the Exradin A‐16 ionization chamber, in the center of the phantom. In conclusion, the PMMA shell filled with water required a small correction for the determination of the absorbed dose, while remaining within the statistical uncertainty of the calculations (±0.71). Plastic Water™ and polystyrene shells can be used without correction. There is a potential advantage to measuring the 4 mm helmet output using these spherical water phantoms.

PACS numbers: 87.10.Rt, 87.50.cm

Measurement of relative output factors for the 8 and 4 mm collimators of Leksell Gamma Knife Perfexion by film dosimetry

Three types of films, Kodak EDR2, Gafchromic EBT, and Gafchromic MD-V2-55, were used to measure relative output factors of 4 and 8 mm collimators of the Leksell Gamma Knife Perfexion. The optical density to dose calibration curve for each of the film types was obtained by exposing the films to a range of known doses. Ten data points were acquired for each of the calibration curves in the dose ranges from 0 to 4 Gy, 0 to 8 Gy, and 0 to 80 Gy for Kodak EDR2, Gafchromic EBT, and Gafchromic MD-V2-55 films, respectively. For the measurement of relative output factors, five films of each film type were exposed to a known dose. All films were scanned using EPSON EXPRESSION 10000 XL scanner with 200 dpi resolution in 16 bit gray scale for EDR2 film and 48 bit color scale for Gafchromic films. The scanned images were imported in the red channel for both Gafchromic films. The background corrections from an unexposed film were applied to all films. The output factors obtained from film measurements were in a close agreement both with the Monte Carlo calculated values of 0.924 and 0.805 for 8 and 4 mm collimators, respectively. These values are provided by the vendor and used as default values in the vendor's treatment planning system. The largest differences were noted for the Kodak EDR 2 films (-2.1% and -4.5% for 8 and 4 mm collimators, respectively). The best agreement observed was for EBT Gafchromic film (-0.8% and +0.6% differences for 8 and 4 mm collimators, respectively). Based on the present values, no changes in the default relative output factor values were made in the treatment planning system.