Validation of complex radiotherapy techniques using polymer gel dosimetry

Modern radiotherapy delivers highly conformal dose distributions to irregularly shaped target volumes while sparing the surrounding normal tissue. Due to the complex planning and delivery techniques, dose verification and validation of the whole treatment workflow by end-to-end tests became much more important and polymer gel dosimeters are one of the few possibilities to capture the delivered dose distribution in 3D. The basic principles and formulations of gel dosimetry and its evaluation methods are described and the available studies validating device-specific geometrical parameters as well as the dose delivery by advanced radiotherapy techniques, such as 3D-CRT/IMRT and stereotactic radiosurgery treatments, the treatment of moving targets, online-adaptive magnetic resonance-guided radiotherapy as well as proton and ion beam treatments, are reviewed. The present status and limitations as well as future challenges of polymer gel dosimetry for the validation of complex radiotherapy techniques are discussed.

Radiation Dosimetry by Use of Radiosensitive Hydrogels and Polymers: Mechanisms, State-of-the-Art and Perspective from 3D to 4D

Gel dosimetry was developed in the 1990s in response to a growing need for methods to validate the radiation dose distribution delivered to cancer patients receiving high-precision radiotherapy. Three different classes of gel dosimeters were developed and extensively studied. The first class of gel dosimeters is the Fricke gel dosimeters, which consist of a hydrogel with dissolved ferrous ions that oxidize upon exposure to ionizing radiation. The oxidation results in a change in the nuclear magnetic resonance (NMR) relaxation, which makes it possible to read out Fricke gel dosimeters by use of quantitative magnetic resonance imaging (MRI). The radiation-induced oxidation in Fricke gel dosimeters can also be visualized by adding an indicator such as xylenol orange. The second class of gel dosimeters is the radiochromic gel dosimeters, which also exhibit a color change upon irradiation but do not use a metal ion. These radiochromic gel dosimeters do not demonstrate a significant radiation-induced change in NMR properties. The third class is the polymer gel dosimeters, which contain vinyl monomers that polymerize upon irradiation. Polymer gel dosimeters are predominantly read out by quantitative MRI or X-ray CT. The accuracy of the dosimeters depends on both the physico-chemical properties of the gel dosimeters and on the readout technique. Many different gel formulations have been proposed and discussed in the scientific literature in the last three decades, and scanning methods have been optimized to achieve an acceptable accuracy for clinical dosimetry. More recently, with the introduction of the MR-Linac, which combines an MRI-scanner and a clinical linear accelerator in one, it was shown possible to acquire dose maps during radiation, but new challenges arise.

Target localization accuracy in frame‐based stereotactic radiosurgery: Comparison between MR‐only and MR/CT co‐registration approaches

Purpose

In frame‐based Gamma Knife (GK) stereotactic radiosurgery two treatment planning workflows are commonly employed; one based solely on magnetic resonance (MR) images and the other based on magnetic resonance/computed tomography (MR/CT) co‐registered images. In both workflows, target localization accuracy (TLA) can be deteriorated due to MR‐related geometric distortions and/or MR/CT co‐registration uncertainties. In this study, the overall TLA following both clinical workflows is evaluated for cases of multiple brain metastases.

Methods

A polymer gel‐filled head phantom, having the Leksell stereotactic headframe attached, was CT‐imaged and irradiated by a GK Perfexion unit. A total of 26 4‐mm shots were delivered at 26 locations directly defined in the Leksell stereotactic space (LSS), inducing adequate contrast in corresponding T2‐weighted (T2w) MR images. Prescribed shot coordinates served as reference locations. An additional MR scan was acquired to implement the “mean image” distortion correction technique. The TLA for each workflow was assessed by comparing the radiation‐induced target locations, identified in MR images, with corresponding reference locations. Using T1w MR and CT images of 15 patients (totaling 81 lesions), TLA in clinical cases was similarly assessed, considering MR‐corrected data as reference. For the MR/CT workflow, both global and region of interest (ROI)‐based MR/CT registration approaches were studied.

Results

In phantom measurements, the MR‐corrected workflow demonstrated unsurpassed TLA (median offset of 0.2 mm) which deteriorated for MR‐only and MR/CT workflows (median offsets of 0.8 and 0.6 mm, respectively). In real‐patient cases, the MR‐only workflow resulted in offsets that exhibit a significant positive correlation with the distance from the MR isocenter, reaching 1.1 mm (median 0.6 mm). Comparable results were obtained for the MR/CT‐global workflow, although a maximum offset of 1.4 mm was detected. TLA was improved with the MR/CT‐ROI workflow resulting in median/maximum offsets of 0.4 mm/1.1 mm.

Conclusions

Subpixel TLA is achievable in all workflows. For the MR/CT workflow, a ROI‐based MR/CT co‐registration approach could considerably increase TLA and should be preferred instead of a global registration.

Simplified method for determining dose to a non-water phantom through the use of ND,w and IAEA TRS 483 for the GammaPod

PURPOSE

GammaPod, a breast stereotactic radiosurgery device, utilizes 25 rotating Co-60 sources to deliver highly conformal dose distributions. The GammaPod system requires that reference dosimetry be performed in a specific vendor-supplied poly-methylmethacrylate (PMMA) phantom. The nonstandard nature of GammaPod dosimetry, in both the phantom material and machine-specific reference (msr), prohibits use of the American Association of Physicists in Medicine Task Group 51 (TG-51) protocol. This study proposes a practical method using TRS 483 to make the reference dosimetry procedure simpler and to reduce overall uncertainties.

METHODS

The dose to PMMA (D_PMMA) is determined under msr conditions using TRS 483 with an Exradin A1SL chamber placed in a PMMA phantom. The conversion factor, which converts from the dose-to-water (D_w) in broad-beam Co-60 reference geometry to D_PMMA in the msr small field Co-60 (Q_msr) geometry, is derived using the Monte Carlo simulations and procedure described in TRS 483.

RESULTS

The new conversion factor value for an Exradin A1SL chamber is 0.974. When combined with N_D,w, D_PMMA differs by 0.5% from the TG-21/N_x method and 0.2% from the IROC values. Uncertainty decreased from 2.2% to 1.6%.

CONCLUSION

We successfully implemented TRS 483 reference dosimetry protocols utilizing N_D,w for the GammaPod in the PMMA phantom. These results show not only agreement between measurements performed with the previously published method and independent thermoluminescent dosimetry measurements but also reductions in uncertainty. This also provides readers with a pathway to develop their own IAEA TRS 483 factor for any new small field machine that may be developed.

A novel and effective method for validation and measurement of output factors for Leksell Gamma Knife® Icon™ using TRS 483 protocol

The objective of this work was to identify the exact location of the effective point of measurement (EPM) of four different active detectors to compare the relative collimator output factors (ROF) of Leksell Gamma Knife (LGK) according to IAEA TRS-483 recommendations. ROF was measured at the center of the spherical LGK-Solid Water (LGK-SW) Phantom for three (4-, 8-, and 16-mm in diameter) collimators using four (PTW-TN60008, PTW-TN60016, PTW-TN60017, and PTW-60019 diode/diamond) detectors. Since diode detectors have a much smaller sensitive volume than the PTW-31010 ion chamber used for reference dosimetry, its EPM might not be at the center of the phantom, or (100, 100, 100) of the Leksell Coordinate System, particularly in the z-direction. Hence for each diode detector, a CBCT image was acquired after it was inserted into the phantom, from which the z-Leksell coordinate of EPM was determined. Relative collimator output factors was then measured by focusing GK beams on the determined EPM of each diode. Measured ROFs were compared with the vendor-provided values in GK treatment planning system. For validation, a plan was generated to measure the output of 4-mm collimator for PTW-TN60017 at various couch locations along the z-axis. For PTW-TN60008, the percentage variations were 0.6 ± 0.4%, and -0.8 ± 0.2% for 4 and 8-mm collimators, respectively. For PTW-TN60016, the percentage variations were 0.8 ± 0.0%, and 0.2 ± 0.1%, respectively. The percentage variations were -3.3 ± 0.0% and -0.9 ± 0.1%, respectively, for PTW-TN60017, and -0.5 ± 0.0% and -0.8 ± 0.2%, respectively, for PTW-TN60019. Center of the measured profile for PTW-TN60017 was only 0.1 mm different from that identified using the CBCT. In conclusion, we have developed a simple and effective method to determine the EPMs of diode detectors when inserted into the existing LGK-SW phantom. With the acquired positional information and using TRS-483 protocol, good agreements were obtained between the measured ROFs and manufacturer recommended values.

A systematic review of clinical applications of polymer gel dosimeters in radiotherapy

Radiotherapy has rapidly improved because of the use of new equipment and techniques. Hence, the appeal for a feasible and accurate three-dimensional (3D) dosimetry system has increased. In this regard, gel dosimetry systems are accurate 3D dosimeters with high resolution. This systematic review evaluates the clinical applications of polymer gel dosimeters in radiotherapy. To find the clinical applications of polymer gel dosimeters in radiotherapy, a full systematic literature search was performed on the basis of the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines in electronic databases up to January 31, 2017, with use of search-related terms in the titles and abstracts of articles. A total of 765 articles were screened in accordance with our inclusion and exclusion criteria. Eventually, 53 articles were included in the study. The findings show that most clinical applications of polymer gel dosimeters relate to external radiotherapy. Most of the gel dosimeters studied have acceptable dose accuracy as a 3D dosimeter with high resolution. It is difficult to judge which is the best polymer gel dosimeter to use in a clinical setting, because each gel dosimeter has advantages and limitations. For example, methacrylic acid-based gel dosimeters have high dose sensitivity and low toxicity, while their dose response is beam energy dependent; in contrast, N-isopropylacrylamide gel dosimeters have low dose resolution, but their sensitivity is lower and they are relatively toxic.

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.

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.

PRESAGE 3D dosimetry accurately measures Gamma Knife output factors

Small-field output factor measurements are traditionally very difficult because of steep dose gradients, loss of lateral electronic equilibrium, and dose volume averaging in finitely sized detectors. Three-dimensional (3D) dosimetry is ideal for measuring small output factors and avoids many of these potential challenges of point and 2D detectors. PRESAGE 3D polymer dosimeters were used to measure the output factors for the 4 mm and 8 mm collimators of the Leksell Perfexion Gamma Knife radiosurgery treatment system. Discrepancies between the planned and measured distance between shot centers were also investigated. A Gamma Knife head frame was mounted onto an anthropomorphic head phantom. Special inserts were machined to hold 60 mm diameter, 70 mm tall cylindrical PRESAGE dosimeters. The phantom was irradiated with one 16 mm shot and either one 4 mm or one 8 mm shot, to a prescribed dose of either 3 Gy or 4 Gy to the 50% isodose line. The two shots were spaced between 30 mm and 60 mm apart and aligned along the central axis of the cylinder. The Presage dosimeters were measured using the DMOS-RPC optical CT scanning system. Five independent 4 mm output factor measurements fell within 2% of the manufacturer's Monte Carlo simulation-derived nominal value, as did two independent 8 mm output factor measurements. The measured distances between shot centers varied by ± 0.8 mm with respect to the planned shot displacements. On the basis of these results, we conclude that PRESAGE dosimetry is excellently suited to quantify the difficult-to-measure Gamma Knife output factors.

Assessment and characterization of the total geometric uncertainty in Gamma Knife radiosurgery using polymer gels

PURPOSE

This work proposes and implements an experimental methodology, based on polymer gels, for assessing the total geometric uncertainty and characterizing its contributors in Gamma Knife (GK) radiosurgery.

METHODS

A treatment plan consisting of 26, 4-mm GK single shot dose distributions, covering an extended region of the Leksell stereotactic space, was prepared and delivered to a polymer gel filled polymethyl methacrylate (PMMA) head phantom (16 cm diameter) used to accurately reproduce every link in the GK treatment chain. The center of each shot served as a "control point" in the assessment of the GK total geometric uncertainty, which depends on (a) the spatial dose delivery uncertainty of the PERFEXION GK unit used in this work, (b) the spatial distortions inherent in MR images commonly used for target delineation, and (c) the geometric uncertainty contributor associated with the image registration procedure performed by the Leksell GammaPlan (LGP) treatment planning system (TPS), in the case that registration is directly based on the apparent fiducial locations depicted in each MR image by the N-shaped rods on the Leksell localization box. The irradiated phantom was MR imaged at 1.5 T employing a T2-weighted pulse sequence. Four image series were acquired by alternating the frequency encoding axis and reversing the read gradient polarity, thus allowing the characterization of the MR-related spatial distortions.

RESULTS

MR spatial distortions stemming from main field (B0) inhomogeneity as well as from susceptibility and chemical shift phenomena (also known as sequence dependent distortions) were found to be of the order of 0.5 mm, while those owing to gradient nonlinearities (also known as sequence independent distortions) were found to increase with distance from the MR scanner isocenter extending up to 0.47 mm at an Euclidean distance of 69.6 mm. Regarding the LGP image registration procedure, the corresponding average contribution to the total geometric uncertainty ranged from 0.34 to 0.80 mm. The average total geometric uncertainty, which also includes the GK spatial dose delivery uncertainty, was found equal to (0.88 ± 0.16), (0.88 ± 0.26), (1.02 ± 0.09), and (1.15 ± 0.24) mm for the MR image series acquired with the read gradient polarity (direction) set toward right, left, posterior, and anterior, respectively.

CONCLUSIONS

The implemented methodology seems capable of assessing the total geometric uncertainty, as well as of characterizing its contributors, ascribed to the entire GK treatment delivery (i.e., from MR imaging to GK dose delivery) for an extended region of the Leksell stereotactic space. Results obtained indicate that the selection of both the frequency encoding axis and the read gradient polarity during MRI acquisition may affect the magnitude as well as the spatial components of the total geometric uncertainty.

MRI-based polymer gel dosimetry for validating plans with multiple matrices in Gamma Knife stereotactic radiosurgery

One of treatment planning techniques with Leksell GammaPlan (LGP) for Gamma Knife stereotactic radiosurgery (GKSRS) uses multiple matrices with multiple dose prescriptions. Computational complexity increases when shots are placed in multiple matrices with different grid sizes. Hence, the experimental validation of LGP calculated dose distributions is needed for those cases. For the current study, we used BANG3 polymer gel contained in a head-sized glass bottle to simulate the entire treatment process of GKSRS. A treatment plan with three 18 mm shots and one 8 mm shot in separate matrices was created with LGP. The prescribed maximum dose was 8 Gy to three shots and 16 Gy to one of the 18 mm shots. The 3D dose distribution recorded in the gel dosimeter was read using a Siemens 3T MRI scanner. The scanning parameters of a CPMG pulse sequence with 32 equidistant echoes were as follows: TR = 7 s, echo step = 13.6 ms, field-of-view = 256 mm × 256 mm, and pixel size = 1 mm × 1 mm. Interleaved acquisition mode was used to obtain 15 to 45 2-mm-thick slices. Using a calibration relationship between absorbed dose and the spin-spin relaxation rate (R2), we converted R2 images to dose images. MATLAB-based in-house programs were used for R2 estimation and dose comparison. Gamma-index analysis for the 3D data showed gamma values less than unity for 86% of the voxels. Through this study we accomplished the first application of polymer gel dosimetry for a true comparison between measured 3D dose distributions and LGP calculations for plans using multiple matrices for multiple targets.