Scott Crowe

Monte Carlo-based diode design for correction-less small field dosimetry

Due to their small collecting volume, diodes are commonly used in small field dosimetry. However, the relative sensitivity of a diode increases with decreasing small field size. Conversely, small air gaps have been shown to cause a significant decrease in the sensitivity of a detector as the field size is decreased. Therefore, this study uses Monte Carlo simulations to look at introducing air upstream to diodes such that they measure with a constant sensitivity across all field sizes in small field dosimetry. Varying thicknesses of air were introduced onto the upstream end of two commercial diodes (PTW 60016 photon diode and PTW 60017 electron diode), as well as a theoretical unenclosed silicon chip using field sizes as small as 5 mm × 5 mm. The metric D(w,Q)/D(Det,Q) used in this study represents the ratio of the dose to a point of water to the dose to the diode active volume, for a particular field size and location. The optimal thickness of air required to provide a constant sensitivity across all small field sizes was found by plotting D(w,Q)/D(Det,Q) as a function of introduced air gap size for various field sizes, and finding the intersection point of these plots. That is, the point at which D(w,Q)/D(Det,Q) was constant for all field sizes was found. The optimal thickness of air was calculated to be 3.3, 1.15 and 0.10 mm for the photon diode, electron diode and unenclosed silicon chip, respectively. The variation in these results was due to the different design of each detector. When calculated with the new diode design incorporating the upstream air gap, k(f(clin),f(msr))(Q(clin),Q(msr)) was equal to unity to within statistical uncertainty (0.5%) for all three diodes. Cross-axis profile measurements were also improved with the new detector design. The upstream air gap could be implanted on the commercial diodes via a cap consisting of the air cavity surrounded by water equivalent material. The results for the unclosed silicon chip show that an ideal small field dosimetry diode could be created by using a silicon chip with a small amount of air above it.

Technical Note: Modeling a complex micro‐multileaf collimator using the standard BEAMnrc distribution

PURPOSE The component modules in the standard BEAMnrc istribution may appear to be insufficient to model micro-multileaf collimators that have trifaceted leaf ends and complex leaf profiles. This note indicates, however, that accurate Monte Carlo simulations of radiotherapy beams defined by a complex collimation device can be completed using BEAMnrcs standard VARMLC component module. METHODS That this simple collimator model can produce spatially and dosimetrically accurate microcollimated fields is illustrated using comparisons with ion chamber and film measurements of the dose deposited by square and irregular fields incident on planar, homogeneous water phantoms. RESULTS Monte Carlo dose calculations for on-axis and off-axis fields are shown to produce good agreement with experimental values, even on close examination of the penumbrae. CONCLUSIONS The use of a VARMLC model of the micro-multileaf collimator, along with a commissioned model of the associated linear accelerator, is therefore recommended as an alternative to the development or use of in-house or third-party component modules for simulating stereotactic radiotherapy and radiosurgery treatments. Simulation parameters for the VARMLC model are provided which should allow other researchers to adapt and use this model to study clinical stereotactic radiotherapy treatments.

THE DEVELOPMENT OF A MONTE CARLO SYSTEM TO VERIFY RADIOTHERAPY TREATMENT DOSE CALCULATIONS

Recent advances in the planning and delivery of radiotherapy treatments have resulted in improvements in the accuracy and precision with which therapeutic radiation can be administered. As the complexity of the treatments increases it becomes more difficult to predict the dose distribution in the patient accurately. Monte Carlo methods have the potential to improve the accuracy of the dose calculations and are increasingly being recognised as the “gold standard” for predicting dose deposition in the patient. In this study, software has been developed that enables the transfer of treatment plan information from the treatment planning system to a Monte Carlo dose calculation engine. A database of commissioned linear accelerator models (Elekta Precise and Varian 2100CD at various energies) has been developed using the EGSnrc/BEAMnrc Monte Carlo suite. Planned beam descriptions and CT images can be exported from the treatment planning system using the DICOM framework. The information in these files is combined with an appropriate linear accelerator model to allow the accurate calculation of the radiation field incident on a modelled patient geometry. The Monte Carlo dose calculation results are combined according to the monitor units specified in the exported plan. The result is a 3D dose distribution that could be used to verify treatment planning system calculations. The software, MCDTK (Monte Carlo Dicom ToolKit), has been developed in the Java programming language and produces BEAMnrc and DOSXYZnrc input files, ready for submission on a high-performance computing cluster. The code has been tested with the Eclipse (Varian Medical Systems), Oncentra MasterPlan (Nucletron B.V.) and Pinnacle3 (Philips Medical Systems) planning systems. In this study the software was validated against measurements in homogenous and heterogeneous phantoms. Monte Carlo models are commissioned through comparison with quality assurance measurements made using a large square field incident on a homogenous volume of water. This study aims to provide a valuable confirmation that Monte Carlo calculations match experimental measurements for complex fields and heterogeneous media.

An experimental extrapolation technique using the Gafchromic EBT3 film for relative output factor measurements in small x-ray fields

PURPOSE An experimental extrapolation technique is presented, which can be used to determine the relative output factors for very small x-ray fields using the Gafchromic EBT3 film. METHODS Relative output factors were measured for the Brainlab SRS cones ranging in diameters from 4 to 30 mm(2) on a Novalis Trilogy linear accelerator with 6 MV SRS x-rays. The relative output factor was determined from an experimental reducing circular region of interest (ROI) extrapolation technique developed to remove the effects of volume averaging. This was achieved by scanning the EBT3 film measurements with a high scanning resolution of 1200 dpi. From the high resolution scans, the size of the circular regions of interest was varied to produce a plot of relative output factors versus area of analysis. The plot was then extrapolated to zero to determine the relative output factor corresponding to zero volume. RESULTS Results have shown that for a 4 mm field size, the extrapolated relative output factor was measured as a value of 0.651 ± 0.018 as compared to 0.639 ± 0.019 and 0.633 ± 0.021 for 0.5 and 1.0 mm diameter of analysis values, respectively. This showed a change in the relative output factors of 1.8% and 2.8% at these comparative regions of interest sizes. In comparison, the 25 mm cone had negligible differences in the measured output factor between zero extrapolation, 0.5 and 1.0 mm diameter ROIs, respectively. CONCLUSIONS This work shows that for very small fields such as 4.0 mm cone sizes, a measureable difference can be seen in the relative output factor based on the circular ROI and the size of the area of analysis using radiochromic film dosimetry. The authors recommend to scan the Gafchromic EBT3 film at a resolution of 1200 dpi for cone sizes less than 7.5 mm and to utilize an extrapolation technique for the output factor measurements of very small field dosimetry.

Treatment plan complexity metrics for predicting IMRT pre-treatment quality assurance results.

The planning of IMRT treatments requires a compromise between dose conformity (complexity) and deliverability. This study investigates established and novel treatment complexity metrics for 122 IMRT beams from prostate treatment plans. The Treatment and Dose Assessor software was used to extract the necessary data from exported treatment plan files and calculate the metrics. For most of the metrics, there was strong overlap between the calculated values for plans that passed and failed their quality assurance (QA) tests. However, statistically significant variation between plans that passed and failed QA measurements was found for the established modulation index and for a novel metric describing the proportion of small apertures in each beam. The ‘small aperture score’ provided threshold values which successfully distinguished deliverable treatment plans from plans that did not pass QA, with a low false negative rate.

Retrospective evaluation of dosimetric quality for prostate carcinomas treated with 3D conformal, intensity modulated and volumetric modulated arc radiotherapy

This study examines and compares the dosimetric quality of radiotherapy treatment plans for prostate carcinoma across a cohort of 163 patients treated across five centres: 83 treated with three‐dimensional conformal radiotherapy (3DCRT), 33 treated with intensity modulated radiotherapy (IMRT) and 47 treated with volumetric modulated arc therapy (VMAT).

Technical Note: Relationships between gamma criteria and action levels: Results of a multicenter audit of gamma agreement index results.

PURPOSE The aim of this work was to use a multicenter audit of modulated radiotherapy quality assurance (QA) data to provide a practical examination of gamma evaluation criteria and action level selection. The use of the gamma evaluation method for patient-specific pretreatment QA is widespread, with most commercial solutions implementing the method. METHODS Gamma agreement indices were calculated using the criteria 1%/1 mm, 2%/2 mm, 2%/3 mm, 3%/2 mm, 3%/3 mm, and 5%/3 mm for 1265 pretreatment QA measurements, planned at seven treatment centers, using four different treatment planning systems, delivered using three different delivery systems (intensity-modulated radiation therapy, volumetric-modulated arc therapy, and helical tomotherapy) and measured using three different dose measurement systems. The sensitivity of each pair of gamma criteria was evaluated relative to the gamma agreement indices calculated using 3%/3 mm. RESULTS A linear relationship was observed for 2%/2 mm, 2%/3 mm, and 3%/2 mm. This result implies that most beams failing at 3%/3 mm would also fail for those criteria, if the action level was adjusted appropriately. Some borderline plans might be passed or failed depending on the relative priority (tighter tolerance) used for dose difference or distance to agreement evaluation. Dosimeter resolution and treatment modality were found to have a smaller effect on the results of QA measurements than the number of dimensions (2D or 3D) over which the gamma evaluation was calculated. CONCLUSIONS This work provides a method (and a large sample of results) for calculating equivalent action levels for different gamma evaluation criteria. This work constitutes a valuable guide for clinical decision making and a means to compare published gamma evaluation results from studies using different evaluation criteria. More generally, the data provided by this work support the recommendation that gamma criteria that specifically prioritize the property of greatest clinical importance for each treatment modality of anatomical site should be selected when using gamma evaluations for modulated radiotherapy QA. It is therefore suggested that departments using the gamma evaluation as a QA analysis tool should consider the relative importance of dose difference and distance to agreement, when selecting gamma evaluation criteria.

Monte Carlo verification of gel dosimetry measurements for stereotactic radiotherapy.

The quality assurance of stereotactic radiotherapy and radiosurgery treatments requires the use of small-field dose measurements that can be experimentally challenging. This study used Monte Carlo simulations to establish that PAGAT dosimetry gel can be used to provide accurate, high-resolution, three-dimensional dose measurements of stereotactic radiotherapy fields. A small cylindrical container (4 cm height, 4.2 cm diameter) was filled with PAGAT gel, placed in the parietal region inside a CIRS head phantom and irradiated with a 12-field stereotactic radiotherapy plan. The resulting three-dimensional dose measurement was read out using an optical CT scanner and compared with the treatment planning prediction of the dose delivered to the gel during the treatment. A BEAMnrc/DOSXYZnrc simulation of this treatment was completed, to provide a standard against which the accuracy of the gel measurement could be gauged. The three-dimensional dose distributions obtained from Monte Carlo and from the gel measurement were found to be in better agreement with each other than with the dose distribution provided by the treatment planning systems pencil beam calculation. Both sets of data showed close agreement with the treatment planning systems dose distribution through the centre of the irradiated volume and substantial disagreement with the treatment planning system at the penumbrae. The Monte Carlo calculations and gel measurements both indicated that the treated volume was up to 3 mm narrower, with steeper penumbrae and more variable out-of-field dose, than predicted by the treatment planning system. The Monte Carlo simulations allowed the accuracy of the PAGAT gel dosimeter to be verified in this case, allowing PAGAT gel to be utilized in the measurement of dose from stereotactic and other radiotherapy treatments, with greater confidence in the future.

Use of 3D printed materials as tissue-equivalent phantoms

This study used the specific example of 3D printing with acrylonitrile butadiene styrene (ABS) as a means to investigate the potential usefulness of benchtop rapid prototyping as a technique for producing patient specific phantoms for radiotherapy dosimetry. Three small cylinders and one model of a human lung were produced via in-house 3D printing with ABS, using 90%, 50%, 30% and 10% ABS infill densities. These phantom samples were evaluated in terms of their geometric accuracy, tissue equivalence and radiation hardness, when irradiated using a range of clinical radiotherapy beams. The measured dimensions of the small cylindrical phantoms all matched their planned dimensions, within 1mm. The lung phantom was less accurately matched to the lung geometry on which it was based, due to simplifications introduced during the phantom design process. The mass densities, electron densities and linear attenuation coefficients identified using CT data, as well as the results of film measurements made using megavoltage photon and electron beams, indicated that phantoms printed with ABS, using infill densities of 30% or more, are potentially useful as lung- and tissue-equivalent phantoms for patient-specific radiotherapy dosimetry. All cylindrical 3D printed phantom samples were found to be unaffected by prolonged radiation and to accurately match their design specifications. However, care should be taken to avoid oversimplifying anatomical structures when printing more complex phantoms.

Technical Note: Preliminary investigations into the use of a functionalised polymer to reduce diffusion in Fricke gel dosimeters

PURPOSE A modification of the existing PVA-FX hydrogel has been made to investigate the use of a functionalised polymer in a Fricke gel dosimetry system to decrease Fe(3+) diffusion. METHODS The chelating agent, xylenol orange, was chemically bonded to the gelling agent, polyvinyl alcohol (PVA) to create xylenol orange functionalised PVA (XO-PVA). A gel was created from the XO-PVA (20% w/v) with ferrous sulfate (0.4 mM) and sulfuric acid (50 mM). RESULTS This resulted in an optical density dose sensitivity of 0.014 Gy(-1), an auto-oxidation rate of 0.0005 h(-1), and a diffusion rate of 0.129 mm(2) h(-1); an 8% reduction compared to the original PVA-FX gel, which in practical terms adds approximately 1 h to the time span between irradiation and accurate read-out. CONCLUSIONS Because this initial method of chemically bonding xylenol orange to polyvinyl alcohol has inherently low conversion, the improvement on existing gel systems is minimal when compared to the drawbacks. More efficient methods of functionalising polyvinyl alcohol with xylenol orange must be developed for this system to gain clinical relevance.