Tanya Kairn

Evaluation of a Gafchromic EBT2 film dosimetry system for radiotherapy quality assurance

This study examines the dosimetric accuracy of Gafchromic EBT2 model radiochromic film for use in radiotherapy quality assurance. In this study, film was scanned using an Epson Perfection V700 flatbed scanner in transmission mode at 75 DPI with the subsequent analysis performed using the red and blue colour channels and ImageJ software. Results of this study suggest that the conversion of film optical density to measured dose should, at present, utilise red channel data only, without application of a blue channel correction to the data. For the batch of film examined here, film uniformity and reproducibility appear to have improved compared with published results using older batches. The orientation of the film on the scanner and the side of the film facing the light source were found to have substantial effects on results. Based on the results of this study, it is possible to recommend the use of EBT2 film in routine quality assurance testing for radiotherapy, in situations where a dose uncertainty of up to 2.8% is acceptable.

Local heterogeneities in early batches of EBT2 film: a suggested solution

To enhance the utility of radiochromic films for high-resolution dosimetry of small and modulated radiotherapy fields, we propose a means to negate the effects of heterogeneities in EBT2 (and other) films. The results of using our simple procedure for evaluating radiation dose in EBT2 film are compared with the results of using the manufacturers recommended procedure as well as a procedure previously established for evaluating dose in older EBT film. It is shown that Newtons ring-like scanning artefacts can be avoided through the use of a plastic frame, to elevate the film above the scanners surface. The effects of film heterogeneity can be minimized by evaluating net optical density, pixelwise, as the logarithm of the ratio of the red-channel pixel value in each pixel of each irradiated film to the red-channel pixel value in the same pixel in the same film prior to irradiation. The application of a blue-channel correction was found to result in increased noise. It is recommended that, when using EBT2 film for radiotherapy quality assurance, the films should be scanned before and after irradiation and analysed using the method proposed herein, without the use of the blue-channel correction, in order to produce dose images with minimal film heterogeneity effects.

Molecular-dynamics simulation of model polymer nanocomposite rheology and comparison with experiment

The shear-rate dependence of viscosity is studied for model polymer melts containing various concentrations of spherical filler particles by molecular-dynamics simulations, and the results are compared with the experimental results for calcium-carbonate-filled polypropylene. Although there are some significant differences in scale between the simulated model polymer composite and the system used in the experiments, some important qualitative similarities in shear behavior are observed. The trends in the steady-state shear viscosities of the simulated polymer-filler system agree with those seen in the experimental results; shear viscosities, zero-shear viscosities, and the rate of shear thinning are all seen to increase with filler content in both the experimental and simulated systems. We observe a significant difference between the filler volume fraction dependence of the zero-shear viscosity of the simulated system and that of the experimental system that can be attributed to a large difference in the ratio of the filler particle radius to the radius of gyration of the polymer molecules. In the simulated system, the filler particles are so small that they only have a weak effect on the viscosity of the composite at low filler volume fraction, but in the experimental system, the viscosity of the composite increases rapidly with increasing filler volume fraction. Our results indicate that there exists a value of the ratio of the filler particle radius to the polymer radius of gyration such that the zero-shear-rate viscosity of the composite becomes approximately independent of the filler particle volume fraction.

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.

Radiotherapy treatment verification using radiological thickness measured with an amorphous silicon electronic portal imaging device: Monte Carlo simulation and experiment

This work validates the use of an amorphous-silicon, flat-panel electronic portal imaging device (a-Si EPID) for use as a gauge of patient or phantom radiological thickness, as an alternative to dosimetry. The response of the a-Si EPID is calibrated by adapting a technique previously applied to scanning liquid ion chamber EPIDs, and the stability, accuracy and reliability of this calibration are explored in detail. We find that the stability of this calibration, between different linacs at the same centre, is sufficient to justify calibrating only one of the EPIDs every month and using the calibration data thus obtained to perform measurements on all of the other linacs. Radiological thickness is shown to provide a reliable means of relating experimental measurements to the results of BEAMnrc Monte Carlo simulations of the linac-phantom-EPID system. For these reasons we suggest that radiological thickness can be used to verify radiotherapy treatment delivery and identify changes in the treatment field, patient position and target location, as well as patient physical thickness.

Adapting a generic BEAMnrc model of the BrainLAB m3 micro-multileaf collimator to simulate a local collimation device

This work is focussed on developing a commissioning procedure so that a Monte Carlo model, which uses BEAMnrcs standard VARMLC component module, can be adapted to match a specific BrainLAB m3 micro-multileaf collimator (microMLC). A set of measurements are recommended, for use as a reference against which the model can be tested and optimized. These include radiochromic film measurements of dose from small and offset fields, as well as measurements of microMLC transmission and interleaf leakage. Simulations and measurements to obtain microMLC scatter factors are shown to be insensitive to relevant model parameters and are therefore not recommended, unless the output of the linear accelerator model is in doubt. Ultimately, this note provides detailed instructions for those intending to optimize a VARMLC model to match the dose delivered by their local BrainLAB m3 microMLC device.

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.

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).