A validation framework to assess performance of commercial deformable image registration in lung radiotherapy.

INTRODUCTION: Deformable image registration (DIR) can play an important role in the context of adaptive radiotherapy. The AAPM Task Group 132 (TG-132) has described several quantitative measures for DIR error assessment but they can only be accurately defined when there is a ground-truth present in high-contrast regions. This work aims to set out a framework to obtain optimal results for CT-CT lung DIR in clinical setting for a commercially available system by quantifying the DIR performance in both low- and high-contrast regions. METHODS: Five publicly available thorax datasets were used to assess the DIR quality. A "Ghost fiducial" method was implemented by windowing the contrast in a new feature provided by Varian Velocity v4.1. Target registration error (TRE) of the landmarks and Dice-similarity coefficient of the tumour were calculated at three different contrast settings to assess the algorithm in high- and low-contrast scenarios. RESULTS: For the original unedited dataset, higher resolution DIR methods showed best performance acceptable within the recommended limit according to TG-132, when actual displacements were less than 10 mm. The relation of the actual displacement of the landmarks and TRE shows the limited capacity of the algorithm to deal with movements larger than 10 mm. CONCLUSION: This work found the performance of DIR methods and settings available in Varian Velocity v4.1 to be a function of contrast level as well as extent of motion. This highlights the need for multiple metrics to assess different aspects of DIR performance for various applications related to low-contrast and/or high-contrast regions.

Out-of-field in vivo dosimetry using TLD in SABR for primary kidney cancer involving mixed photon fields.

PURPOSE: To assess out-of-field dose using three different variants of LiF thermoluminescence dosimeters (TLD) for ten patients who underwent stereotactic ablative body radiotherapy (SABR) for primary renal cell carcinoma (RCC) and compare with treatment planning system (TPS) dose calculations. METHODS AND MATERIALS: Thermoluminescent dosimeter (TLD) measurements were conducted at 20, 30, 40 and 50cm from isocentre on ten patients undergoing SABR for primary RCC. Three types of high-sensitivity LiF:Mg,Cu,P TLD material with different 6Li/7Li isotope ratios were used. Patient plans were calculated using Eclipse Anisotropic Analytical Algorithm (AAA) for clinical evaluation and recalculated using Pencil Beam Convolution (PBC) algorithm for comparison. RESULTS: Both AAA and PBC showed diminished accuracy for photon doses at increasing distance out-of-field. At 50cm, measured photon dose was 0.3cGy normalised to a 10Gy prescription on average with only small variation across all patients. This is likely due to the leakage component of the out-of-field dose. The 6Li-enriched TLD materials showed increased signal attributable to additional neutron contribution. CONCLUSION: LiF:Mg,Cu,P TLD containing 6Li is sensitive enough to measure out-of-field dose 50cm from isocentre however will over-estimate the photon component of out-of-field dose in high energy treatments due to the presence of thermal neutrons. 7Li enriched materials which are insensitive to neutrons are therefore required for accurate photon dosimetry. Neutron signal has been shown here to increase with MUs and is higher for patients treated using certain non coplanar beam arrangements. Further work is required to convert this additional neutron signal to dose.

DOSE AND GAMMA-RAY SPECTRA FROM NEUTRON-INDUCED RADIOACTIVITY IN MEDICAL LINEAR ACCELERATORS FOLLOWING HIGH-ENERGY TOTAL BODY IRRADIATION.

Production of radioisotopes in medical linear accelerators (linacs) is of concern when the beam energy exceeds the threshold for the photonuclear interaction. Staff and patients may receive a radiation dose as a result of the induced radioactivity in the linac. Gamma-ray spectroscopy was used to identify the isotopes produced following the delivery of 18 MV photon beams from a Varian 21EX and an Elekta Synergy. The prominent radioisotopes produced include 187W, 63Zn, 56Mn, 24Na and 28Al in both linac models. The dose rate was measured at the beam exit window (12.6 µSv in the first 10 min) following 18 MV total body irradiation (TBI) beams. For a throughput of 24 TBI patients per year, staff members are estimated to receive an annual dose of up to 750 µSv at the patient location. This can be further reduced to 65 µSv by closing the jaws before re-entering the treatment bunker.

Effect of light source instability on uniformity of 3D reconstructions from a cone beam optical CT scanner.

Temporally varying light intensity during acquisition of projection images in an optical CT scanner can potentially be misinterpreted as physical properties of the sample. This work investigated the impact of LED light source intensity instability on measured attenuation coefficients. Different scenarios were investigated by conducting one or both of the reference and data scans in a 'cold' scanner, where the light source intensity had not yet stabilised. Uniform samples were scanned to assess the impact on measured uniformity. The orange (590 nm) light source decreased in intensity by 29 % over the first 2 h, while the red (633 nm) decreased by 9 %. The rates of change of intensity at 2 h were 0.1 and 0.03 % respectively over a 5 min period-corresponding to the scan duration. The normalisation function of the reconstruction software does not fully account for the intensity differences and discrepancies remain. Attenuation coefficient inaccuracies of up to 8 % were observed for data reconstructed from projection images acquired with a cold scanner. Increased noise was observed for most cases where one or both of the scans was acquired without sufficient warm-up. The decrease in accuracy and increase in noise were most apparent for data reconstructed from reference and data scans acquired with a cold scanner on different days.

A collimated detection system for assessing leakage dose from medical linear accelerators at the patient plane.

Leakage radiation from linear accelerators can make a significant contribution to healthy tissue dose in patients undergoing radiotherapy. In this work thermoluminescent dosimeters (LiF:Mg,Cu,P TLD chips) were used in a focused lead cone loaded with TLD chips for the purpose of evaluating leakage dose at the patient plane. By placing the TLDs at one end of a stereotactic cone, a focused measurement device is created; this was tested both in and out of the primary beam of a Varian 21-iX linac using 6 MV photons. Acrylic build up material of 1.2 cm thickness was used inside the cone and measurements made with either one or three TLD chips at a given distance from the target. Comparing the readings of three dosimeters in one plane inside the cone offered information regarding the orientation of the cone relative to a radiation source. Measurements in the patient plane with the linac gantry at various angles demonstrated that leakage dose was approximately 0.01% of the primary beam out of field when the cone was pointed directly towards the target and 0.0025% elsewhere (due to scatter within the gantry). No specific 'hot spots' (e.g., insufficient shielding or gaps at abutments) were observed. Focused cone measurements facilitate leakage dose measurements from the linac head directly at the patient plane and allow one to infer the fraction of leakage due to 'direct' photons (along the ray-path from the bremsstrahlung target) and that due to scattered photons.

Source position verification and dosimetry in HDR brachytherapy using an EPID.

PURPOSE: Accurate treatment delivery in high dose rate (HDR) brachytherapy requires correct source dwell positions and dwell times to be administered relative to each other and to the surrounding anatomy. Treatment delivery inaccuracies predominantly occur for two reasons: (i) anatomical movement or (ii) as a result of human errors that are usually related to incorrect implementation of the planned treatment. Electronic portal imaging devices (EPIDs) were originally developed for patient position verification in external beam radiotherapy and their application has been extended to provide dosimetric information. The authors have characterized the response of an EPID for use with an (192)Ir brachytherapy source to demonstrate its use as a verification device, providing both source position and dosimetric information. METHODS: Characterization of the EPID response using an (192)Ir brachytherapy source included investigations of reproducibility, linearity with dose rate, photon energy dependence, and charge build-up effects associated with exposure time and image acquisition time. Source position resolution in three dimensions was determined. To illustrate treatment verification, a simple treatment plan was delivered to a phantom and the measured EPID dose distribution compared with the planned dose. RESULTS: The mean absolute source position error in the plane parallel to the EPID, for dwells measured at 50, 100, and 150 mm source to detector distances (SDD), was determined to be 0.26 mm. The resolution of the z coordinate (perpendicular distance from detector plane) is SDD dependent with 95% confidence intervals of ± 0.1, ± 0.5, and ± 2.0 mm at SDDs of 50, 100, and 150 mm, respectively. The response of the EPID is highly linear to dose rate. The EPID exhibits an over-response to low energy incident photons and this nonlinearity is incorporated into the dose calibration procedure. A distance (spectral) dependent dose rate calibration procedure has been developed. The difference between measured and planned dose is less than 2% for 98.0% of pixels in a two-dimensional plane at an SDD of 100 mm. CONCLUSIONS: Our application of EPID dosimetry to HDR brachytherapy provides a quality assurance measure of the geometrical distribution of the delivered dose as well as the source positions, which is not possible with any current HDR brachytherapy verification system.

Performance of 12 DIR algorithms in low-contrast regions for mass and density conserving deformation.

PURPOSE: Deformable image registration (DIR) has become a key tool for adaptive radiotherapy to account for inter- and intrafraction organ deformation. Of contemporary interest, the application to deformable dose accumulation requires accurate deformation even in low contrast regions where dose gradients may exist within near-uniform tissues. One expects high-contrast features to generally be deformed more accurately by DIR algorithms. The authors systematically assess the accuracy of 12 DIR algorithms and quantitatively examine, in particular, low-contrast regions, where accuracy has not previously been established. METHODS: This work investigates DIR algorithms in three dimensions using deformable gel (DEFGEL) [U. J. Yeo, M. L. Taylor, L. Dunn, R. L. Smith, T. Kron, and R. D. Franich, "A novel methodology for 3D deformable dosimetry," Med. Phys. 39, 2203-2213 (2012)], for application to mass- and density-conserving deformations. CT images of DEFGEL phantoms with 16 fiducial markers (FMs) implanted were acquired in deformed and undeformed states for three different representative deformation geometries. Nonrigid image registration was performed using 12 common algorithms in the public domain. The optimum parameter setup was identified for each algorithm and each was tested for deformation accuracy in three scenarios: (I) original images of the DEFGEL with 16 FMs; (II) images with eight of the FMs mathematically erased; and (III) images with all FMs mathematically erased. The deformation vector fields obtained for scenarios II and III were then applied to the original images containing all 16 FMs. The locations of the FMs estimated by the algorithms were compared to actual locations determined by CT imaging. The accuracy of the algorithms was assessed by evaluation of three-dimensional vectors between true marker locations and predicted marker locations. RESULTS: The mean magnitude of 16 error vectors per sample ranged from 0.3 to 3.7, 1.0 to 6.3, and 1.3 to 7.5 mm across algorithms for scenarios I to III, respectively. The greatest accuracy was exhibited by the original Horn and Schunck optical flow algorithm. In this case, for scenario III (erased FMs not contributing to driving the DIR calculation), the mean error was half that of the modified demons algorithm (which exhibited the greatest error), across all deformations. Some algorithms failed to reproduce the geometry at all, while others accurately deformed high contrast features but not low-contrast regions-indicating poor interpolation between landmarks. CONCLUSIONS: The accuracy of DIR algorithms was quantitatively evaluated using a tissue equivalent, mass, and density conserving DEFGEL phantom. For the model studied, optical flow algorithms performed better than demons algorithms, with the original Horn and Schunck performing best. The degree of error is influenced more by the magnitude of displacement than the geometric complexity of the deformation. As might be expected, deformation is estimated less accurately for low-contrast regions than for high-contrast features, and the method presented here allows quantitative analysis of the differences. The evaluation of registration accuracy through observation of the same high contrast features that drive the DIR calculation is shown to be circular and hence misleading.

Direct MC conversion of absorbed dose to graphite to absorbed dose to water for 60Co radiation.

The ARPANSA calibration service for (60)Co gamma rays is based on a primary standard graphite calorimeter that measures absorbed dose to graphite. Measurements with the calorimeter are converted to the absorbed dose to water using the calculation of the ratio of the absorbed dose in the calorimeter to the absorbed dose in a water phantom. ARPANSA has recently changed the basis of this calculation from a photon fluence scaling method to a direct Monte Carlo (MC) calculation. The MC conversion uses an EGSnrc model of the cobalt source that has been validated against water tank and graphite phantom measurements, a step that is required to quantify uncertainties in the underlying interaction coefficients in the MC code. A comparison with the Bureau International des Poids et Mesures (BIPM) as part of the key comparison BIPM.RI(I)-K4 showed an agreement of 0.9973 (53).

Spectral differences in 6 MV beams with matched PDDs and the effect on chamber response.

The Australian Radiation Protection and Nuclear Safety Agency (ARPANSA) has installed an Elekta Synergy platform linac to establish a direct megavoltage primary standard calibration service, instead of relying on calibrations derived from (60)Co. One of the 6 MV beams of the ARPANSA linac has been approximately matched to the Varian high energy platform 6 MV photon beam. The electron beam energy was adjusted to match the percentage depth dose (PDD) curve and TPR(20,10). This work calculates the error introduced when using a calibration factor from this Elekta Synergy Platform linac on a Varian high-energy platform beam at 6 MV. Monte Carlo models of the Varian and matched Elekta accelerator accurately predict the measured PDDs and profiles, but show significantly different energy spectra, resulting mainly from differences in target thickness between the two accelerators. Monte Carlo modelling of the energy correction factor k(Q) of a secondary standard NE2561 chamber shows a difference of 0.4% between the Varian and the Varian-matched Elekta beams. Although small, this is a significant discrepancy for primary standard calibrations. Similar variations are expected for chambers of similar construction, and additional variations may occur with other linac manufacturers. The work has also investigated the design of a custom flattening filter to precisely match the energy spectrum of the Varian beam on the Elekta platform.

A phantom for testing of 4D-CT for radiotherapy of small lesions.

PURPOSE: The use of time-resolved four-dimensional computed tomography (4D-CT) in radiotherapy requires strict quality assurance to ensure the accuracy of motion management protocols. The aim of this work was to design and test a phantom capable of large amplitude motion for use in 4D-CT, with particular interest in small lesions typical for stereotactic body radiotherapy. METHODS: The phantom of "see-saw" design is light weight, capable of including various sample materials and compatible with several surrogate marker signal acquisition systems. It is constructed of polymethylmethacrylate (Perspex) and its movement is controlled via a dc motor and drive wheel. It was tested using two CT scanners with different 4D acquisition methods: the Philips Brilliance Big Bore CT (helical scan, pressure belt) and a General Electric Discovery STE PET∕CT (axial scan, infrared marker). Amplitudes ranging from 1.5 to 6.0 cm and frequencies of up to 40 cycles per minute were used to study the effect of motion on image quality. Maximum intensity projections (MIPs), as well as average intensity projections (AIPs) of moving objects were investigated and their quality dependence on the number of phase reconstruction bins assessed. RESULTS: CT number discrepancies between moving and stationary objects were found to have no systematic dependence on amplitude, frequency, or specific interphase variability. MIP-delineated amplitudes of motion were found to match physical phantom amplitudes to within 2 mm for all motion scenarios tested. Objects undergoing large amplitude motions (>3.0 cm) were shown to cause artefacts in MIP and AIP projections when ten phase bins were assigned. This problem can be mitigated by increasing the number of phase bins in a 4D-CT scan. CONCLUSIONS: The phantom was found to be a suitable tool for evaluating the image quality of 4D-CT motion management technology, as well as providing a quality assurance tool for intercenter∕intervendor testing of commercial 4D-CT systems. When imaging objects with large amplitudes, the completeness criterion described here indicates the number of phase bins required to prevent missing data in MIPs and AIPs. This is most relevant for small lesions undergoing large motions.

Is it sensible to "deform" dose? 3D experimental validation of dose-warping.

PURPOSE: Strategies for dose accumulation in deforming anatomy are of interest in radiotherapy. Algorithms exist for the deformation of dose based on patient image sets, though these are sometimes contentious because not all such image calculations are constrained by physical laws. While tumor and organ motion has been a key area of study for a considerable amount of time, deformation is of increasing interest. In this work, we demonstrate a full 3D experimental validation of results from a range of dose deformation algorithms available in the public domain. METHODS: We recently developed the first tissue-equivalent, full 3D deformable dosimetric phantom-"DEFGEL." To assess the accuracy of dose-warping based on deformable image registration (DIR), we have measured doses in undeformed and deformed states of the DEFGEL dosimeter and compared these to planned doses and warped doses. In this way we have directly evaluated the accuracy of dose-warping calculations for 11 different algorithms. We have done this for a range of stereotactic irradiation schemes and types and magnitudes of deformation. RESULTS: The original Horn and Schunck algorithm is shown to be the best performing of the 11 algorithms trialled. Comparing measured and dose-warped calculations for this method, it is found that for a 10 × 10 mm(2) square field, γ(3%∕3mm) = 99.9%; for a 20 × 20 mm(2) cross-shaped field, γ(3%∕3mm) = 99.1%; and for a multiple dynamic arc (0.413 cm(3) PTV) treatment adapted from a patient treatment plan, γ(3%∕3mm) = 95%. In each case, the agreement is comparable to-but consistently ∼1% less than-comparison between measured and calculated (planned) dose distributions in the absence of deformation. The magnitude of the deformation, as measured by the largest displacement experienced by any voxel in the volume, has the greatest influence on the accuracy of the warped dose distribution. Considering the square field case, the smallest deformation (∼9 mm) yields agreement of γ(3%∕3mm) = 99.9%, while the most significant deformation (∼20 mm) yields agreement of γ(3%∕3mm) = 96.7%. CONCLUSIONS: We have confirmed that, for a range of mass and density conserving deformations representative of those observable in anatomical targets, DIR-based dose-warping can yield accurate predictions of the dose distribution. Substantial differences can be seen between the results of different algorithms indicating that DIR performance should be scrutinized before application todose-warping. We have demonstrated that the DEFGEL deformable dosimeter can be used to evaluate DIR performance and the accuracy of dose-warping results by direct measurement.

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 system's pencil beam calculation. Both sets of data showed close agreement with the treatment planning system's 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.

Robust calculation of effective atomic numbers: the Auto-Z(eff) software.

PURPOSE: The most appropriate method of evaluating the effective atomic number necessitates consideration of energy-dependent behavior. Previously, this required quite laborious calculation, which is why many scientists revert to over-simplistic power-law methods. The purpose of this work is to develop user-friendly software for the robust, energy-dependent computation of effective atomic numbers relevant within the context of medical physics, superseding the commonly employed simplistic power law approaches. METHOD: Visual Basic was used to develop a GUI allowing the straightforward calculation of effective atomic numbers. Photon interaction cross section matrices are constructed for energies spanning 10 keV to 10 GeV and elements Z = 1-100. Coefficients for composite media are constructed via linear additivity of the fractional constituents and contrasted against the precalculated matrices at each energy, thereby associating an effective atomic number through interpolation of adjacent cross section data. Uncertainties are of the order of 1-2%. RESULTS: Auto-Z(eff) allows rapid (∼0.6 s) calculation of effective atomic numbers for a range of predefined or user-specified media, allowing estimation of radiological properties and comparison of different media (for instance assessment of water equivalence). The accuracy of Auto-Z(eff) has been validated against numerous published theoretical and experimental predictions, demonstrating good agreement. The results also show that commonly employed power-law approaches are inaccurate, even in their intended regime of applicability (i.e., photoelectric regime). Furthermore, comparing the effective atomic numbers of composite materials using power-law approaches even in a relative fashion is shown to be inappropriate. CONCLUSION: Auto-Z(eff) facilitates easy computation of effective atomic numbers as a function of energy, as well as average and spectral-weighted means. The results are significantly more accurate than normal power-law predictions. The software is freely available to interested readers, who are encouraged to contact the authors.

A novel methodology for 3D deformable dosimetry.

PURPOSE: Interfraction and intrafraction variation in anatomic structures is a significant challenge in contemporary radiotherapy. The objective of this work is to develop a novel tool for deformable structure dosimetry, using a tissue-equivalent deformable gel dosimeter that can reproducibly simulate targets subject to deformation. This will enable direct measurement of integrated doses delivered in different deformation states, and the verification of dose deforming algorithms. METHODS: A modified version of the nPAG polymer gel has been used as a deformable 3D dosimeter and phantom to investigate doses delivered to deforming tissue-equivalent geometry. The deformable gel (DEFGEL) dosimeter/phantom is comprised of polymer gel in a latex membrane, moulded (in this case) into a cylindrical geometry, and deformed with an acrylic compressor. Fifteen aluminium fiducial markers (FM) were implanted into DEFGEL phantoms and the reproducibility of deformation was determined via multiple computed tomography (CT) scans in deformed and nondeformed states before and after multiple (up to 150) deformations. Dose was delivered to the DEFGEL phantom in three arrangements: (i) without deformation, (ii) with deformation, and (iii) cumulative exposures with and without deformation, i.e., dose integration. Irradiations included both square field and a stereotactic multiple dynamic arc treatment adapted from a patient plan. Doses delivered to the DEFGEL phantom were read out using cone beam optical CT. RESULTS: Reproducibility was verified by observation of interscan shifts of FM locations (as determined via CT), measured from an absolute reference point and in terms of inter-FM distance. The majority (76%) of points exhibited zero shift, with others shifting by one pixel size consistent with setup error as confirmed with a control sample. Comparison of dose profiles and 2D isodose distributions from the three arrangements illustrated complex spatial redistribution of dose in all three dimensions occurring as a result of the change in shape of the target between irradiations, even for a relatively simple deformation. Discrepancies of up to 30% of the maximum dose were evident from dose difference maps for three orthogonal planes taken through the isocenter of a stereotactic field. CONCLUSIONS: This paper describes the first use of a tissue-equivalent, 3D dose-integrating deformable phantom that yields integrated or redistributed dosimetric information. The proposed methodology readily yields three-dimensional (3D) dosimetric data from radiation delivery to the DEFGEL phantom in deformed and undeformed states. The impacts of deformation on dose distributions were readily seen in the isodose contours and line profiles from the three arrangements. It is demonstrated that the system is potentially capable of reproducibly emulating the physical deformation of an organ, and therefore can be used to evaluate absorbed doses to deformable targets and organs at risk in three dimensions and to validate deformation algorithms applied to dose distributions.

Assessment of leakage doses around the treatment heads of different linear accelerators.

Out-of-field doses to untargeted organs may have long-term detrimental health effects for patients treated with radiotherapy. It has been observed that equivalent treatments delivered to patients with different accelerators may result in significant differences in the out-of-field dose. In this work, the points of leakage dose are identified about the gantry of several treatment units. The origin of the observed higher doses is investigated. LiF:Mg,Cu,P thermoluminescent dosimetry has been employed to quantify the dose at a several points around the linac head of various linear accelerators (linacs): a Varian 600C, Varian 21-iX, Siemens Primus and Elekta Synergy-II. Comparisons are also made between different energy modes, collimator rotations and field sizes. Significant differences in leaked photon doses were identified when comparing the various linac models. The isocentric-waveguide 600C generally exhibits the lowest leakage directed towards the patient. The Siemens and Elekta models generally produce a greater leakage than the Varian models. The leakage 'hotspots' are evident on the gantry section housing the waveguide on the 21-iX. For all machines, there are significant differences in the x and y directions. Larger field sizes result in a greater leakage at the interface plate. There is a greater leakage around the waveguide when operating in a low-energy mode, but a greater leakage for the high-energy mode at the linac face. Of the vendors investigated, the Varian 600C showed the lowest average leakage dose. The Varian 21-iX showed double the dose of the 600C. The Elekta Synergy-II had on average four times the dose leakage than the 600C, and the Siemens Primus showed an average of five times that of the 600C. All vendors show strong differences in the x and y directions. The results offer the potential for patient-positioning strategies, linac choice and shielding strategies to reduce the leakage dose to patients.

A programmable motion phantom for quality assurance of motion management in radiotherapy.

A commercially available motion phantom (QUASAR, Modus Medical) was modified for programmable motion control with the aim of reproducing patient respiratory motion in one dimension in both the anterior-posterior and superior-inferior directions, as well as, providing controllable breath-hold and sinusoidal patterns for the testing of radiotherapy gating systems. In order to simulate realistic patient motion, the DC motor was replaced by a stepper motor. A separate 'chest-wall' motion platform was also designed to accommodate a variety of surrogate marker systems. The platform employs a second stepper motor that allows for the decoupling of the chest-wall and insert motion. The platform's accuracy was tested by replicating patient traces recorded with the Varian real-time position management (RPM) system and comparing the motion platform's recorded motion trace with the original patient data. Six lung cancer patient traces recorded with the RPM system were uploaded to the motion platform's in-house control software and subsequently replicated through the phantom motion platform. The phantom's motion profile was recorded with the RPM system and compared to the original patient data. Sinusoidal and breath-hold patterns were simulated with the motion platform and recorded with the RPM system to verify the systems potential for routine quality assurance of commercial radiotherapy gating systems. There was good correlation between replicated and actual patient data (P 0.003). Mean differences between the location of maxima in replicated and patient data-sets for six patients amounted to 0.034 cm with the corresponding minima mean equal to 0.010 cm. The upgraded motion phantom was found to replicate patient motion accurately as well as provide useful test patterns to aid in the quality assurance of motion management methods and technologies.

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 BEAMnrc's 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.

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 BEAMnrc's 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.

Stereotactic fields shaped with a micro-multileaf collimator: systematic characterization of peripheral dose.

Despite the highly localized doses that may be delivered via stereotactic radiotherapy, a small dose is nonetheless delivered to out-of-field regions, which may cause detriment to the patient. In this work, a systematic set of dose measurements have been undertaken up to a distance of 45 cm from the isocentre, for stereotactic fields shaped by a BrainLAB mini-multileaf collimator (MMLC) mounted on a Varian 600C linear accelerator. A range of treatment parameters were varied so as to determine the factors of greatest influence and establish relationships with dose. The commercial treatment planning software (TPS) miscalculates the dose to out-of-field regions. Measured dose decreases consistently out to 45 cm, whereas the TPS decreases out to 10-15 cm, at which point the predicted dose is constant. At the 5-10 cm off-axis distance (OAD), measurements indicate doses of about 5-10% of the dose at the isocentre, 1% at 15 cm OAD and 0.1% at 45 cm OAD. There are several observed trends. Greater MMLC field sizes (with static jaw) result in higher out-of-field dose, as do shallower depths. The source-to-surface distance does not greatly influence peripheral dose. However, the results given in this work do indicate that simple treatment arrangements, such as preferable collimator rotation, would in certain cases reduce out-of-field dose by an order of magnitude. Peripheral dose raises questions of treatment optimization, particularly in cases where patients have a long life expectancy in which secondary effects may become manifest, such as in the treatment of paediatric patients or those with a non-malignant primary. For instance, for a 20 Gy hypo-fractionated treatment, dose to out-of-field regions is of the order of cGy-a substantial dose in radiation protection terms.