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

Dose-response dependencies of OSL dosimeters in conventional linacs and 1.5T MR-linacs: an experimental and Monte Carlo study

Objective. This work focuses on the optically stimulated luminescence dosimetry (OSLD) dose-response characterization, with emphasis on 1.5T MR-Linacs.Approach. Throughout this study, the nanoDots OSLDs (Landauer, USA) were considered. In groups of three, the mean OSLD response was measured in a conventional linac and an MR-Linac under various irradiation conditions to investigate (i) dose-response linearity with and without the 1.5T magnetic field, (ii) signal fading rate and its dependencies, (iii) beam quality, detector orientation and dose rate dependencies in a conventional linac, (iii) potential MR imaging related effects on OSLD response and (iv) detector orientation dependence in an MR-Linac. Monte Carlo calculations were performed to further quantify angular dependence after rotating the detector around its central axis parallel to the magnetic field, and determine the magnetic field correction factors,kB,Q,for all cardinal detector orientations.Main results. OSLD dose-response supralinearity in an MR-Linac setting was found to agree within uncertainties with the corresponding one in a conventional linac, for the axial detector orientation investigated. Signal fading rate does not depend on irradiation conditions for the range of 3-30 d considered. OSLD angular (orientation) dependence is more pronounced under the presence of a magnetic field. OSLDs irradiated with and without real-time T2w MR imaging enabled during irradiation yielded the same response within uncertainties.kB,Qvalues were determined for all three cardinal orientations. Corrections needed reached up to 6.4%. However, if OSLDs are calibrated in the axial orientation and then irradiated in an MR-Linac placed again in the axial orientation (perpendicular to the magnetic field), then simulations suggest thatkB,Qcan be considered unity within uncertainties, irrespective of the incident beam angle.Significance. This work contributes towards OSLD dose-response characterization and relevant correction factors availability. OSLDs are suitable for QA checks in MR-based beam gating applications andin vivodosimetry in MR-Linacs.

Geometric distortion assessment in 3T MR images used for treatment planning in cranial Stereotactic Radiosurgery and Radiotherapy

Magnetic Resonance images (MRIs) are employed in brain Stereotactic Radiosurgery and Radiotherapy (SRS/SRT) for target and/or critical organ localization and delineation. However, MRIs are inherently distorted, which also impacts the accuracy of the Magnetic Resonance Imaging/Computed Tomography (MRI/CT) co-registration process. In this phantom-based study, geometric distortion is assessed in 3T T2-weighted images (T2WIs), while the efficacy of an MRI distortion correction technique is also evaluated. A homogeneous polymer gel-filled phantom was CT-imaged before being irradiated with 26 4-mm Gamma Knife shots at predefined locations (reference control points). The irradiated phantom was MRI-scanned at 3T, implementing a T2-weighted protocol suitable for SRS/SRT treatment planning. The centers of mass of all shots were identified in the 3D image space by implementing an iterative localization algorithm and served as the evaluated control points for MRI distortion detection. MRIs and CT images were spatially co-registered using a mutual information algorithm. The inverse transformation matrix was applied to the reference control points and compared with the corresponding MRI-identified ones to evaluate the overall spatial accuracy of the MRI/CT dataset. The mean image distortion correction technique was implemented, and resulting MRI-corrected control points were compared against the corresponding reference ones. For the scanning parameters used, increased MRI distortion (>1mm) was detected at areas distant from the MRI isocenter (>5cm), while median radial distortion was 0.76mm. Detected offsets were slightly higher for the MRI/CT dataset (0.92mm median distortion). The mean image distortion correction improves geometric accuracy, but residual distortion cannot be considered negligible (0.51mm median distortion). For all three datasets studied, a statistically significant positive correlation between detected spatial offsets and their distance from the MRI isocenter was revealed. This work contributes towards the wider adoption of 3T imaging in SRS/SRT treatment planning. The presented methodology can be employed in commissioning and quality assurance programmes of corresponding treatment workflows.

Dose kernel decomposition for spot-based radiotherapy treatment planning

PURPOSE

Pre-calculation of accurate dose deposition kernels for treatment planning of spot-based radiotherapies, such as Gamma Knife (GK) and Gamma Pod (GP), can be very time-consuming and may require large data storage with an enormous number of possible spots. We proposed a novel kernel decomposition (KD) model to address accurate and fast (real-time) dose calculation with reduced data storage requirements for spot-based treatment planning. The application of the KD model was demonstrated for clinical GK and GP radiotherapy platforms.

METHODS

The dose deposition kernel at each spot (shot position) is modeled as the product of a shift-invariant kernel based on a reference kernel and spatially variant scale factor. The reference kernel, one for each collimator, is defined at the center of the commissioning phantom for GK and at the center of the treatment target for GP and calculated using the Monte Carlo (MC) method. The spatially variant scale factor is defined as the ratio of the mean tissue maximum ratio (TMR) at the candidate shot position to that at the reference kernel position, and the mean TMR map is calculated within the entire volume through parallel beam ray tracing on the density image followed by averaging over all source directions. The proposed KD dose calculations were compared with the MC method and with the GK and GP treatment planning system (TPS) computations for various shot positions and collimator sizes utilizing a phantom and 14 and 12 clinical plans for GK and GP, respectively.

RESULTS

For the phantom study, the KD Gamma index (3%/1 mm) passing rates were greater than 99% (median 100%) relative to the MC doses, except for the shots close to the boundary. The passing rates dropped below 90% for 8 mm (16 mm) shots positioned within ∼1 cm (∼2 cm) of the boundary. For the clinical GK plans, the KD Gamma passing rates were greater than 99% (median 100%) compared to the MC and greater than 92% (median 99%) compared to the TPS. For the clinical GP plans, the KD Gamma passing rates were greater than 95% (median 98%) compared to the MC and greater than 91% (median 97%) compared to the TPS. The scale factors were calculated in sub-seconds with GPU implementation and only need to be calculated once before treatment plan optimization. The calculation of the dose kernel was also within sub-seconds without requiring beam-by-beam calculation commonly done in the TPS.

CONCLUSION

The proposed model can provide an accurate dose and enables real-time dose and derivative calculations by kernel shifting and scaling without pre-calculating or requiring large data storage for GK and GP dose deposition kernels during treatment planning. This model could be useful for spot-based radiotherapy treatment planning by allowing an efficient global fine search for optimal spots.

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.

Dosimetric evaluation of the Leksell GammaPlan™ Convolution dose calculation algorithm

The dosimetric accuracy of the Leksell GammaPlan Convolution calculation algorithm was evaluated through comparison with corresponding Monte Carlo (MC) dosimetric results. MC simulations were based on generated sector phase space files for the 4 mm, 8 mm and 16 mm collimator sizes, using a previous comprehensive Gamma Knife Perfexion^™ source model and validated using film dosimetry. Test cases were designed for the evaluation of the Convolution algorithm involving irradiation of homogeneous and inhomogeneous phantom geometries mimicking clinical cases, with radiation fields created using one sector (single sector), all sectors with the same (single shot) or different (composite shot) collimator sizes. Dose calculations using the Convolution algorithm were found to be in excellent agreement (gamma pass rate greater than 98%, applying 1%/1 mm local dose difference and distance agreement criteria), with corresponding MC calculations, indicating the accuracy of the Convolution algorithm in homogeneous and heterogeneous model geometries. While of minor clinical importance, large deviations were observed for the voxels laying inside air media. The calculated beam on times using the Convolution algorithm were found to increase (up to 7%) relative to the TMR 10 algorithm currently used in clinical practice, especially in a test case mimicking a brain metastasis close to the skull, in excellent agreement with corresponding MC calculations.

Quality assurance for Gamma Knife Perfexion using the Exradin W1 plastic scintillation detector and Lucy phantom

The goal of this work is to validate the use of the Exradin W1 plastic scintillation detector (PSD) to measure profiles and output factors from Gamma Knife Perfexion collimators in a Lucy phantom. The Exradin W1 PSD has a small-volume, near-water-equivalent, energy-independent sensitive element. Output measurements were performed for all 3 collimators (4 mm, 8 mm, and 16 mm) of the Gamma Knife Perfexion system, and these measurements were compared to measurements made with an A16 ion chamber and an EBT3 film and to the nominal values. We showed that a configuration in which the focus or 'shot' moves while the detector remains fixed is essentially equivalent to a configuration in which the focus is fixed while the detector moves. A Lucy phantom containing a PSD was moved in small steps to acquire profiles in all three dimensions. EBT3 film was inserted in the Lucy phantom and exposed to a single shot for each collimator. The relative values for output factors measured with the PSD were 1.000, 0.892, and 0.795, for the 16 mm, 8 mm, and 4 mm collimators, respectively. The values measured with EBT3 film were 1.000, 0.881, and 0.793, and the values measured with the A16 ion chamber were 1.000, 0.883, and 0.727. The nominal output factors for the Gamma Knife Perfexion are 1.000, 0.900, and 0.814, respectively. There was excellent agreement between all profiles measured with the PSD and EBT3 as well as with the treatment planning system data provided by the vendor. In light of our results, the Exradin W1 PSD is well suited for beam quality assurance of a Gamma Knife Perfexion irradiator.

Monte Carlo and experimental determination of correction factors for gamma knife perfexion small field dosimetry measurements

Detector-, field size- and machine-specific correction factors are required for precise dosimetry measurements in small and non-standard photon fields. In this work, Monte Carlo (MC) simulation techniques were used to calculate the [Formula: see text] and [Formula: see text] correction factors for a series of ionization chambers, a synthetic microDiamond and diode dosimeters, used for reference and/or output factor (OF) measurements in the Gamma Knife Perfexion photon fields. Calculations were performed for the solid water (SW) and ABS plastic phantoms, as well as for a water phantom of the same geometry. MC calculations for the [Formula: see text] correction factors in SW were compared against corresponding experimental results for a subset of ionization chambers and diode detectors. Reference experimental OF data were obtained through the weighted average of corresponding measurements using TLDs, EBT-2 films and alanine pellets. [Formula: see text] values close to unity (within 1%) were calculated for most of ionization chambers in water. Greater corrections of up to 6.0% were observed for chambers with relatively large air-cavity dimensions and steel central electrode. A phantom correction of 1.006 and 1.024 (breaking down to 1.014 from the ABS sphere and 1.010 from the accompanying ABS phantom adapter) were calculated for the SW and ABS phantoms, respectively, adding up to [Formula: see text] corrections in water. Both measurements and MC calculations for the diode and microDiamond detectors resulted in lower than unit [Formula: see text] correction factors, due to their denser sensitive volume and encapsulation materials. In comparison, higher than unit [Formula: see text] results for the ionization chambers suggested field size depended dose underestimations (being significant for the 4 mm field), with magnitude depending on the combination of contradicting phenomena associated with volume averaging and electron fluence perturbations. Finally, the presence of 0.5 mm air-gap between the diodes' frontal surface and their phantom-inserts may considerably influence OF measurements, reaching 4.6% for the Razor diode.

Verification of dose profiles generated by the convolution algorithm of the gamma knife® radiosurgery planning system

PURPOSE

A convolution algorithm that takes into account electron-density inhomogeneity was recently introduced to calculate dose distributions for the Gamma Knife (GK) Perfexion™ treatment planning program. The accuracies of the dose distributions computed using the convolution method were assessed using an anthropomorphic phantom and film dosimetry.

METHODS

Absorbed-dose distributions inside a phantom (CIRS Radiosurgery Head Phantom, Model 605) were calculated using the convolution method of the GK treatment-planning software (Leksell Gamma Plan^® version 10.1; LGP) for various combinations of collimator size, location, direction of calculation plane, and number of shots. Computed tomography (CT) images of the phantom and a data set of CT number versus electron density were provided to the LGP. Calculated distributions were exported as digital-image communications in medicine-radiation therapy (DICOM-RT) files. Three types of radiochromic film (GafChromic^® MD-V2-55, MD-V3, and EBT2) were irradiated inside the phantom using GK Perfexion™. Scanned images of the measured films were processed following standard radiochromic film-handling procedures. For a two-dimensional quantitative evaluation, gamma index pass rates (GIPRs) and normalized agreement-test indices (NATIs) were obtained. Image handling and index calculations were performed using a commercial software package (DoseLab Pro version 6.80).

RESULTS

The film-dose calibration data were well fitted with third-order polynomials (R²  ≥ 0.9993). The mean GIPR and NATI of the 93 analyzed films were 99.3 ± 1.1% and 0.8 ± 1.3, respectively, using 3%/1.0 mm criteria. The calculated maximum doses were 4.3 ± 1.7% higher than the measured values for the 4 mm single shots and 1.8 ± 0.7% greater than those for the 8 mm single shots, whereas differences of only 0.3 ± 0.9% were observed for the 16 mm single shots. The accuracy of the calculated distribution was not statistically related to the collimator size, number of shots, or centrality of location (P > 0.05, independent-sample t-test). The plans in the axial planes exhibited poorer agreement with the measured distributions than the plans in the coronal or sagittal planes; however, their GIPR values (≥ 96.9%) were clinically acceptable. The plans for an arbitrary virtual target of volume 1.6 cm³ at an axial plane close to the top of the phantom showed the worst agreement and the greatest fluctuation (GIPR = 96.9 ± 1.2%, NATI = 3.9 ± 1.7).

CONCLUSIONS

The measured accuracies of the dose distributions calculated by the convolution algorithm of the LGP were within the clinically acceptable range (GIPR ≥ 96.9%) for various configurations of collimator size, location, direction of calculation plane, and number of shots. Due to the intrinsic asymmetry in the dose distribution along the z-axis, the treatment plan should also be verified in coronal or sagittal plane.

Stereotactic radiosurgery alone for multiple brain metastases? A review of clinical and technical issues

Over the past three decades several randomized trials have enabled evidence-based practice for patients presenting with limited brain metastases. These trials have focused on the role of surgery or stereotactic radiosurgery (SRS) with or without whole brain radiation therapy (WBRT). As a result, it is clear that local control should be optimized with surgery or SRS in patients with optimal prognostic factors presenting with up to 4 brain metastases. The routine use of adjuvant WBRT remains debatable, as although greater distant brain control rates are observed, there is no impact on survival, and modern outcomes suggest adverse effects from WBRT on patient cognition and quality of life. With dramatic technologic advances in radiation oncology facilitating the adoption of SRS into mainstream practice, the optimal management of patients with multiple brain metastases is now being put forward. Practice is evolving to SRS alone in these patients despite a lack of level 1 evidence to support a clinical departure from WBRT. The purpose of this review is to summarize the current state of the evidence for patients presenting with limited and multiple metastases, and to present an in-depth analysis of the technology and dosimetric issues specific to the treatment of multiple metastases.

Development of a Monte Carlo model for treatment planning dose verification of the Leksell Gamma Knife Perfexion radiosurgery system

Detailed Monte Carlo (MC) modeling of the Leksell Gamma Knife (GK) Perfexion (PFX) collimator system is the only accurate ab initio approach appearing in the literature. As a different approach, in this work, we present a MC model based on film measurement. By adjusting the model parameters and fine-tuning the derived fluence map for each individual source to match the manufacturer's ring output factors, we created a reasonable virtual source model for MC simulations to verify treatment planning dose for the GK PFX radiosurgery system. The MC simulation model was commissioned by simple single shots. Dose profiles and both ring and collimator output factors were compared with the treatment planning system (TPS). Good agreement was achieved for dose profiles especially for the region of plateau (< 2%), while larger difference (< 5%) came from the penumbra region. The maximum difference of the calculated output factor was within 0.7%. The model was further validated by a clinical test case. Good agreement was obtained. The DVHs for brainstem and the skull were almost identical and, for the target, the volume covered by the prescription (12.5 Gy to 50% isodose line) was 95.6% from MC calculation versus 100% from the TPS.