Evaluation of a novel 4D in vivo dosimetry system.

A prototype of a new 4D in vivo dosimetry system capable of simultaneous real-time position monitoring and dose measurement has been developed. The radiation positioning system (RADPOS) is controlled by a computer and combines two technologies: MOSFET radiation detector coupled with an electromagnetic positioning device. Special software has been developed that allows sampling position and dose either manually or automatically in user-defined time intervals. Preliminary tests of the new device include a dosimetric evaluation of the detector in 60Co, 6 MV, and 18 MV beams and measurements of spatial position stability and accuracy. In addition, the effect of metals and other materials on the performance of the positioning system has been investigated. Results show that the RADPOS system can measure in-air dose profiles that agree, on average, within 3%-5% of diode measurements for the energies tested. The response of the detector is isotropic within 1.6% (1 SD) with a maximum deviation of +/- 4.0% over 360 degrees. The maximum variation in the calibration coefficient over field sizes from 6 x 6 to 25 x 25 cm2 was 2.3% for RADPOS probe with the high sensitivity MOSFET and 4.6% for the probe with the standard sensitivity MOSFET. Of the materials tested, only aluminum, lead, and brass caused shifts in the RADPOS read position. The magnitude of the shift varied between materials and size of the material sample. Nonmagnetic stainless steel (Grade 304) caused a distortion of less than 2 mm when placed within 10 mm of the detector; therefore, it can provide a reasonable alternative to other metals if required. The results of the preliminary tests indicate that the device can be used for in vivo dosimetry in 60Co and high-energy beams from linear accelerators.

Radiotherapy dosimetry using a commercial OSL system.

A commercial optically stimulated luminescence (OSL) system developed for radiation protection dosimetry by Landauer, Inc., the InLight microStar reader, was tested for dosimetry procedures in radiotherapy. The system uses carbon-doped aluminum oxide, Al2O3:C, as a radiation detector material. Using this OSL system, a percent depth dose curve for 60Co gamma radiation was measured in solid water. Field size and SSD dependences of the detector response were also evaluated. The dose response relationship was investigated between 25 and 400 cGy. The decay of the response with time following irradiation and the energy dependence of the Al2O3:C OSL detectors were also measured. The results obtained using OSL dosimeters show good agreement with ionization chamber and diode measurements carried out under the same conditions. Reproducibility studies show that the response of the OSL system to repeated exposures is 2.5% (1sd), indicating a real possibility of applying the Landauer OSL commercial system for radiotherapy dosimetric procedures.

Patient-specific Monte Carlo dose calculations for (103)Pd breast brachytherapy.

This work retrospectively investigates patient-specific Monte Carlo (MC) dose calculations for (103)Pd permanent implant breast brachytherapy, exploring various necessary assumptions for deriving virtual patient models: post-implant CT image metallic artifact reduction (MAR), tissue assignment schemes (TAS), and elemental tissue compositions. Three MAR methods (thresholding, 3D median filter, virtual sinogram) are applied to CT images; resulting images are compared to each other and to uncorrected images. Virtual patient models are then derived by application of different TAS ranging from TG-186 basic recommendations (mixed adipose and gland tissue at uniform literature-derived density) to detailed schemes (segmented adipose and gland with CT-derived densities). For detailed schemes, alternate mass density segmentation thresholds between adipose and gland are considered. Several literature-derived elemental compositions for adipose, gland and skin are compared. MC models derived from uncorrected CT images can yield large errors in dose calculations especially when used with detailed TAS. Differences in MAR method result in large differences in local doses when variations in CT number cause differences in tissue assignment. Between different MAR models (same TAS), PTV [Formula: see text] and skin [Formula: see text] each vary by up to 6%. Basic TAS (mixed adipose/gland tissue) generally yield higher dose metrics than detailed segmented schemes: PTV [Formula: see text] and skin [Formula: see text] are higher by up to 13% and 9% respectively. Employing alternate adipose, gland and skin elemental compositions can cause variations in PTV [Formula: see text] of up to 11% and skin [Formula: see text] of up to 30%. Overall, AAPM TG-43 overestimates dose to the PTV ([Formula: see text] on average 10% and up to 27%) and underestimates dose to the skin ([Formula: see text] on average 29% and up to 48%) compared to the various MC models derived using the post-MAR CT images studied herein. The considerable differences between TG-43 and MC models underline the importance of patient-specific MC dose calculations for permanent implant breast brachytherapy. Further, the sensitivity of these MC dose calculations due to necessary assumptions illustrates the importance of developing a consensus modelling approach.

Clinical use of a commercial Monte Carlo treatment planning system for electron beams.

In 2002 we fully implemented clinically a commercial Monte Carlo based treatment planning system for electron beams. The software, developed by MDS Nordion (presently Nucletron), is based on Kawrakow's VMC++ algorithm. The Monte Carlo module is integrated with our Theraplan Plustrade mark treatment planning system. An extensive commissioning process preceded clinical implementation of this software. Using a single virtual 'machine' for each electron beam energy, we can now calculate very accurately the dose distributions and the number of MU for any arbitrary field shape and SSD. This new treatment planning capability has significantly impacted our clinical practice. Since we are more confident of the actual dose delivered to a patient, we now calculate accurate three-dimensional (3D) dose distributions for a greater variety of techniques and anatomical sites than we have in the past. We use the Monte Carlo module to calculate dose for head and neck, breast, chest wall and abdominal treatments with electron beams applied either solo or in conjunction with photons. In some cases patient treatment decisions have been changed, as compared to how such patients would have been treated in the past. In this paper, we present the planning procedure and some clinical examples.

Use of novel fibre-coupled radioluminescence and RADPOS dosimetry systems for total scatter factor measurements in small fields.

A dosimetry system based on Al2O3:C radioluminescence (RL), and RADPOS, a novel 4D dosimetry system using microMOSFETs, were used to measure total scatter factors, (S(c,p))(f(clin))(det), for the CyberKnife robotic radiosugery system. New Monte Carlo calculated correction factors are presented and applied for the RL detector response for the 5, 7.5 and 10 mm collimators in order to correct for the detector geometry and increased photoelectric cross section of Al2O3:C relative to water. For comparison, measurements were also carried out using a micro MOSFET, PTW60012 diode and GAFCHROMIC(®) film (EBT and EBT2). Uncorrected (S(c,p))(f(clin))(det) were obtained by taking the ratio of the detector response for each collimator to that for the 60 mm diameter reference field. Published Monte Carlo calculated correction factors were applied to the RADPOS, microMOSFET and diode detector measurements to yield corrected field factors, Ω(f(clin),f(msr))(Q(clin),Q(msr)), following the terminology of a recent formalism introduced for small and composite field relative dosimetry. With corrections, the RL measured Ω(f(clin),f(msr))(Q(clin),Q(msr)) were 0.656 ± 0.002, 0.815 ± 0.002 and 0.865 ± 0.003 for the 5, 7.5 and 10 mm collimators, respectively. This was in good agreement with RADPOS corrected field factors of 0.650 ± 0.010, 0.816 ± 0.024 and 0.867 ± 0.010 for the 5, 7.5 and 10 mm collimators, respectively. Both RL and RADPOS total scatter factors agreed within approximately two standard deviations of the GAFCHROMIC film values (average of EBT and EBT2) of 0.640 ± 0.006, 0.806 ± 0.007 and 0.859 ± 0.09. Corrected total scatter factors for all dosimetry systems agreed within one standard deviation for collimator sizes 10-60 mm. Our study suggests that the microMOSFET/RADPOS and optical fibre-coupled RL dosimetry system are well suited for total scatter factor measurements over the entire range of field sizes, provided that appropriate correction factors are applied for the collimator diameters smaller than 10 mm.

Commissioning and quality assurance of treatment planning computers.

The process of radiation therapy is complex and involves many steps. At each step, comprehensive quality assurance procedures are required to ensure the safe and accurate delivery of a prescribed radiation dose. This report deals with a comprehensive commissioning and ongoing quality assurance program specifically for treatment planning computers. Detailed guidelines are provided under the following topics: (a) computer program and system documentation and user training, (b) sources of uncertainties and suggested tolerances, (c) initial system checks, (d) repeated system checks, (e) quality assurance through manual procedures, and in vivo dosimetry, and (f) some additional considerations including administration and manpower requirements. In the context of commercial computerized treatment planning systems, uncertainty estimates and achievable criteria of acceptability are presented for: (a) external photon beams, (b) electron beams, (c) brachytherapy, and (d) treatment machine setting calculations. Although these criteria of acceptability appear large, they approach the limit achievable with most of today's treatment planning systems. However, developers of new or improved dose calculation algorithms should strive for the goal recommended by the International Commission of Radiation Units and Measurements of 2% in relative dose accuracy in low dose gradients or 2 mm spatial accuracy in regions with high dose gradients. For brachytherapy, the aim should be 3% accuracy in dose at distances of 0.5 cm or more at any point for any radiation source. Details are provided for initial commissioning tests and follow-up reproducibility tests. The final quality assurance for each patient is to perform an independent manual check of at least one point in the dose distributions, as well as the machine setting calculation. As a check of the overall treatment planning process, in vivo dosimetry should be performed on a select number of patients.

Measurements of the electron dose distribution near inhomogeneities using a plastic scintillation detector.

PURPOSE: Accurate measurement of the electron dose distribution near an inhomogeneity is difficult with traditional dosimeters which themselves perturb the electron field. We tested the performance of a new high resolution, water-equivalent plastic scintillation detector which has ideal properties for this application. METHODS AND MATERIALS: A plastic scintillation detector with a 1 mm diameter, 3 mm long cylindrical sensitive volume was used to measure the dose distributions behind standard benchmark inhomogeneities in water phantoms. The plastic scintillator material is more water equivalent than polystyrene in terms of its mass collision stopping power and mass scattering power. Measurements were performed for beams of electrons having initial energies of 6 and 18 MeV at depths from 0.2-4.2 cm behind the inhomogeneities. RESULTS: The detector reveals hot and cold spots behind heterogeneities at resolutions equivalent to typical film digitizer spot sizes. Plots of the dose distributions behind air, aluminum, lead, and formulations for cortical and inner bone-equivalent materials are presented. CONCLUSION: The plastic scintillation detector is suited for measuring the electron dose distribution near an inhomogeneity.

The survival of aerobic and anoxic human glioma and melanoma cells after irradiation at ultrahigh and clinical dose rates.

This in vitro study was undertaken to determine if ultrahigh dose rates could improve the radiation response of human tumors. Two cell lines, human glioma (U-87 MG), which is radioresistant, and human melanoma (HT-144), which is radiosensitive, were irradiated at ultrahigh and high dose rates under aerobic and anoxic conditions to determine if their oxygen enhancement ratios are modified by dose rate. In fact, the survival curves, and hence the oxygen enhancement ratios, were found to be independent of the dose rate. The oxygen enhancement ratio for glioma cells irradiated in plateau phase was 2.8 (+/- 0.3). The oxygen enhancement ratio was 2.7 (+/- 0.4) for melanoma cells in plateau phase and 2.8 (+/- 0.3) in exponential phase. These results indicate that there is no advantage in treating these tumors using ultrahigh dose rates instead of conventional dose rates.

Evaluation of a dual bias dual metal oxide-silicon semiconductor field effect transistor detector as radiation dosimeter.

A new type of direct reading semiconductor dosimeter has been investigated as a radiation detector for photon and electron therapy beams of various energies. The operation of this device is based on the measurement of the threshold voltage shift in a custom-built metal oxide-silicon semiconductor field effect transistor (MOSFET). This voltage is a linear function of absorbed dose. The extent of the linearity region is dependent on the voltage controlled operation during irradiation. Operating two MOSFETS at two different biases simultaneously during irradiation will result in sensitivity (V/Gy) reproducibility better than +/- 3% over a range in dose of 100 Gy and at a dose per fraction greater than 20 x 10(-2) Gy. The modes of operation give this device many advantages, such as continuous monitoring during irradiation, immediate reading, and permanent storage of total dose after irradiation. The availability and ease of use of these MOSFET detectors make them very promising in clinical dosimetry.

Measurement of Prepl Pwall factors in electron beams and in a 60Co beam for plane-parallel chambers.

This paper describes a method to measure the product of Prepl Pwall correction factors for ionization chambers and presents our measured values of Prepl Pwall for Markus plane-parallel chambers in electron beams. It is shown that the measured values of Prepl Pwall can be fitted to an equation, Prepl Pwall = c1 + c2 R50 + c3 (R50)2, for Markus chambers at the new reference depth for electron beams (6 MeV < or = nominal energy E < or = 20 MeV). We also present our measured values of Prepl Pwall for NACP and Markus chambers in a water phantom irradiated in a 60Co beam.

Measurements of electron beam peak scatter factors.

We have measured the peak scatter factor (PSF) for electron beams which in analogy to photon beams is defined as the ratio of the absorbed dose to water at the depth of dose maximum in a water phantom to the absorbed dose in free air at the same location for a given incident beam. In this study we have measured the PSFs as a function of the incident electron beam energy and field size. We have used both an RK ionization chamber and a silicon diode as radiation detectors. After applying the water/air stopping power ratio and fluence correction factors to the ionization chamber readings, the values of measured PSFs by using two different detectors are in a good agreement. The results show that the value of PSF increases with the increase of the field size and decreases with the increase of the beam energy. The range of the variation with field size between 2 x 2 and 10 x 10 cm2 is 1.24-1.42, 1.15-1.30, 1.15-1.17 for 6, 11, and 20 MeV beams, respectively.

Effects of gold and silver backings on the dose rate around an 125I seed.

Measurements of the effect of either gold or silver backing on the dose rate around an 125I seed were performed using a Therados RFA7 dosimetry system and a small diode detector which was 2.5 mm in diameter and 0.06 mm thick. It was found that the presence of the gold or silver backing modifies the diode response on the side of the 125I seed away from the backing. The effect depends on the backing material and the distance from the seed. There is a small increase close to the gold backing but a decrease further away. This decrease at distances greater than 10 mm from the seed is uniformly 10%, the same as found when the seed is backed by air. There is an increase of up to 25% observed with silver backing the seed and this increase remains significant more than 30 mm from the seed. When the response increases, the results are hard to interpret quantitatively because of variations in the diode response per unit dose with photon energy and extreme sensitivity to geometric changes. Nonetheless, except for the increase at close distances with the gold, the results are in agreement with EGS4 Monte Carlo photon transport simulations which are for a simplified geometry and account for x-ray fluorescence from the K-shell. Furthermore, the increase in the gold-backed case is qualitatively explained by Williamson's Monte Carlo calculations which take into account the L-shell fluorescent x-rays from gold.

GafChromic MD-55: investigated as a precision dosimeter.

GafChromic MD-55 is a fairly new, thin film dosimeter that develops a blue color (lambda max = 676 nm) when irradiated with ionizing radiation. The increase in absorbance is roughly proportional to the absorbed dose. In this study, GafChromic MD-55 was irradiated with 60Co gamma rays. A double irradiation method was used in which a dosimeter is given an unknown dose followed by a known, calibration dose. With this method, GafChromic MD-55 was used to measure doses in the vicinity of 6 Gy with an uncertainty of less than 1%. It was found that the measured optical density of GafChromic MD-55, as presently fabricated, is affected by the polarization of the analyzing light, an important consideration when using GafChromic MD-55 as a precision dosimeter. GafChromic MD-55 was found to consist of seven layers. The response to polarized light was measured for the whole dosimeter and for the three Mylar films which form part of the dosimeter.

The effect of dose rate dependence of p-type silicon detectors on linac relative dosimetry.

Cumulative radiation damage to silicon semiconductor diode detectors can induce dose rate dependent sensitivity, a concern in the pulsed beam of a linac. Two p-Si diode photon detectors were used in this study, diodes A and B. Both were preirradiated by the supplier to 5 kGy, with diode A receiving an estimated 8 kGy from measurements, and diode B, 25 kGy. At 6 MV, the PDD measured with diode B was lower (by 4.4% at a depth of 25 cm) than diode A. Using SSD to vary the dose per pulse from 0.02 to 0.64 mGy/pulse, diode A was dose rate independent (within 2%), while the sensitivity of diode B changed by 13%. Silicon diode detectors should be checked regularly against ionization chambers in the pulsed beam of a linac, especially older high-resistivity diodes that have accumulated dose from high-energy photon beams.

Clinical reference dosimetry: comparison between AAPM TG-21 and TG-51 protocols.

We compare the results of absorbed dose determined at reference conditions according to the AAPM TG-21 dose calibration protocol and the new AAPM TG-51 protocol. The AAPM TG-21 protocol for absorbed dose calibration is based on ionization chambers having exposure calibration factors for 60Co gamma rays, N(x). The new AAPM TG-51 dosimetry protocol for absorbed dose calibration is based on ionization chambers having 60Co absorbed dose-to-water calibration factor, N60Co(D,w). This study shows that the dose changes are within 1% for a cobalt beam, 0.5% for photon energies of 6 and 18 MV, and 2%-3% for electron beams with energies of 6 to 20 MeV. The chamber primary calibration factors, Nx and N60Co(D,w), are traceable to the Canadian primary standards laboratory (NRCC). We also present estimated dose changes between the two protocols when calibration factors are traceable to NIST in the United States.

A safe method of analysis for mechanical damage in spherical radioactive sources used in remote afterloading brachytherapy devices.

The examination of brachytherapy sources in a quality assurance program can result in significant exposure to staff. Even though exposure may be within prescribed limits attempts should be made to reduce such exposure to as low as reasonably achievable (ALARA concept). In this report, a system that can be used to examine small spherical radiation sources while minimizing exposure to staff is described.

Electron fluence correction factors for conversion of dose in plastic to dose in water.

In radiation dosimetry protocols, plastic is allowed as a phantom material for the determination of absorbed dose to water in electron beams. The electron fluence correction factor is needed in conversion of dose measured in plastic to dose in water. There are large discrepancies among recommended values as well as measured values of electron fluence correction factors when polystyrene is used as a phantom material. Using the Monte Carlo technique, we have calculated electron fluence correction factors for incident clinical beam energies between 5 and 50 MeV as a function of depth for clear polystyrene, white polystyrene and PMMA phantom materials and compared the results with those recommended in protocols as well as experimental values from published data. In the Monte Carlo calculations, clinical beams are simulated using the EGS4 user-code BEAM for a variety of medical accelerators. The study shows that our calculated fluence correction factor, phi pw, is a function of depth and incident beam energy Eo with little dependence on other aspects of beam quality. However the phi pw values at dmax are indirectly influenced by the beam quality since they vary with depth and dmax also varies with the beam quality. Calculated phi pw values at dmax are in a range of 1.005-1.045 for a clear polystyrene phantom, 1.005-1.038 for a white polystyrene phantom and 0.996-1.016 for a PMMA phantom. Our values of phi pw are about 1-2% higher than those determined according to the AAPM TG-25 protocol at dmax for clear or white polystyrene. Our calculated values of phi pw also explain some of the variations of measured data because of its depth dependence. A simple formula is derived which gives the electron fluence correction factor phi pw as a function of R50 at dmax or at the depth of 0.6R50-0.1 for any clinical electron beam with energy between 5 and 25 MeV for three plastics: clear polystyrene, white polystyrene and PMMA. The study also makes a careful distinction between phi pw and the corresponding IAEA Code of Practice quantity, hm.

Effects of changes in stopping-power ratios with field size on electron beam relative output factors.

Stopping-power ratios are a function of field size and vary with accelerators. To investigate how these variations affect relative output factor measurements made using ion chambers for electron beams, especially for small fields, (L/rho)air(water) is calculated using the Monte Carlo technique for different field sizes, beam energies, and accelerators and is compared to the data in TG-21 or TG-25, which are for mono-energetic broad beams. For very small field sizes defined by cutouts, if the change in (L/rho)air(water) with dmax is ignored (i.e., TG-25 is not carefully followed), there is an overestimate of relative output factors by up to 3%. Ignoring the field-size effect on stopping-power ratio adds an additional overestimate of up to one-half percent, and using mono-energetic stopping-power ratio data instead of realistic beam data gives another error, but in the opposite direction, of up to 0.7%. Due to the cancellation of these latter two errors, following TG-25 with (L/rho)air(water) data for broad mono-energetic beams will give the correct answer for the ROF measurement within 0.4% compared to using (L/rho)air(water) data for which the field-size effect is considered for realistic electron beams.

Evaluation of the first commercial Monte Carlo dose calculation engine for electron beam treatment planning.

The purpose of this study is to perform a clinical evaluation of the first commercial (MDS Nordion, now Nucletron) treatment planning system for electron beams incorporating Monte Carlo dose calculation module. This software implements Kawrakow's VMC++ voxel-based Monte Carlo calculation algorithm. The accuracy of the dose distribution calculations is evaluated by direct comparisons with extensive sets of measured data in homogeneous and heterogeneous phantoms at different source-to-surface distances (SSDs) and gantry angles. We also verify the accuracy of the Monte Carlo module for monitor unit calculations in comparison with independent hand calculations for homogeneous water phantom at two different SSDs. All electron beams in the range 6-20 MeV are from a Siemens KD-2 linear accelerator. We used 10,000 or 50,000 histories/cm2 in our Monte Carlo calculations, which led to about 2.5% and 1% relative standard error of the mean of the calculated dose. The dose calculation time depends on the number of histories, the number of voxels used to map the patient anatomy, the field size, and the beam energy. The typical run time of the Monte Carlo calculations (10,000 histories/cm2) is 1.02 min on a 2.2 GHz Pentium 4 Xeon computer for a 9 MeV beam, 10 x 10 cm2 field size, incident on the phantom 15 x 15 x 10 cm3 consisting of 31 CT slices and voxels size of 3 x 3 x 3 mm3 (total of 486,720 voxels). We find good agreement (discrepancies smaller than 5%) for most of the tested dose distributions. We also find excellent agreement (discrepancies of 2.5% or less) for the monitor unit calculations relative to the independent manual calculations. The accuracy of monitor unit calculations does not depend on the SSD used, which allows the use of one virtual machine for each beam energy for all arbitrary SSDs. In some cases the test results are found to be sensitive to the voxel size applied such that bigger systematic errors (>5%) occur when large voxel sizes interfere with the extensions of heterogeneities or dose gradients because of differences between the experimental and calculated geometries. Therefore, user control over voxelization is important for high accuracy electron dose calculations.

Monte Carlo investigation of electron beam output factors versus size of square cutout.

A major task in commissioning an electron accelerator is to measure relative output factors versus cutout size (i.e., cutout factors) for various electron beam energies and applicator sizes. We use the BEAM Monte Carlo code [Med Phys. 22, 503-524 (1995)] to stimulate clinical electron beams and to calculate the relative output factors for square cutouts. Calculations are performed for a Siemens MD2 linear accelerator with beam energies, 6, 9, 11, and 13 MeV. The calculated cutout factors for square cutouts in 10 X 10 cm2, 15 X 15 cm2, and 20 X 20 cm2 applicators at SSDs of 100 and 115 cm agree with the measurements made using a silicon diode within about 1% except for the smallest cutouts at SSD= 115 cm where they agree within 0.015. The details of each component of the dose, such as the dose from particles scattered off the jaws and the applicator, the dose from contaminant photons, the dose from direct electrons, etc., are also analyzed. The calculations show that inphantom side-scatter equilibrium is a major factor for the contribution from the direct component which usually dominates the output of a beam. It takes about 6 h of CPU time on a Pentium Pro 200 MHz computer to simulate an accelerator and additional 2 h to calculate the relative output factor for each cutout with a statistical uncertainty of 1%.