A new method of readout in radiochromic film dosimetry

Radiochromic film as a dosimetry medium offers several advantages in high-resolution radiography. A new technique of readout was developed to measure the optical density distributions of the film in purely directed light. This technique implements radiochromic film dosimetry near the film's absorption maximum by using a single-mode top-surface emitting laser diode (675.2 nm). The effective sensitivity of the film, compared with a helium-neon laser densitometer (632.8 nm), is increased approximately threefold. Good accuracy, high spatial resolution and simple assembly of the readout system is achieved. Beam profiles of the four final collimator helmets of a Leksell Gamma Knife (Elekta Inc., Sweden) were experimentally determined. Measured profiles and full-widths at half maximum are consistent with the computer generated data of the dose planning system (Kula 4.4, Elekta Inc., Sweden). The output factor of the 4 mm collimator (the smallest collimator with the steepest dose gradient), essential for the application of well defined doses, was checked. The measurements established an output factor of 826 +/- 9 that lies 9 +/- 1% lower than the adjusted one.     

Dosimetry studies with TLDs for stereotactic radiation techniques for intraocular tumours

Between March 1993 and January 1997, stereotactic radiation techniques were used to irradiate 66 intraocular tumour patients with the Gamma Knife (Leksell Gamma Knife, model B unit) at the University of Vienna, Austria. This study investigates the dosimetry for stereotactic irradiation of ocular structures. For the dosimetry program KULA 4.4, Gamma Knife stereotactic irradiation of the eye represents an extreme frontal skull position. In addition, irradiation of the eye may be performed in the usual supine position in exceptional cases only. With the patient in the prone position, the dose planning program has to calculate with a significantly large number of single-beam extrapolations. In our first experiment we measured the isocentre dose for eight different gamma-angle positions, both in prone and supine positions, using TLD measurements in an Alderson head phantom. We found a maximum deviation of +/- 1.6% using these individually calibrated TLDs. In the second experiment we examined the dose cross profiles for the two most frequently used treatment positions (supine position, gamma = 65 degrees, and prone position, gamma = 140 degrees). For this purpose we implanted a specially designed TLD array into the orbit of a human cadaver head. We found excellent agreement of the dose values measured for the isocentre as well as the posterior part of the eye with orbit with deviations of less than -2.7%. However, for the anterior part of the eye, deviations between computer-generated calculations and the TLD measurements were found to range up to -30%. These differences were noticed both for supine and prone positions. For the Gamma Knife stereotactic irradiation of ocular tumours or pathologies, precautions should be taken to avoid significant underdosage in the anterior part of the radiation field.     

TLD array for precise dose measurements in stereotactic radiation techniques

We developed a new TLD array for precise dose measurement and verification of the spatial dose distribution in small radiation targets. It consists of a hemicylindrical, tissue-equivalent rod made of polystyrene with 17 parallel moulds for an exact positioning of each TLD. The spatial resolution of the TLD array was evaluated using the Leskell spherical phantom. Dose planning was performed with KULA 4.4 under stereotactic conditions on axial CT images. In the Leksell gamma unit the TLD array was irradiated with a maximal dose of 10 Gy with an unplugged 14 mm collimator. The doses delivered to the TLDs were rechecked by diode detector and film dosimetry and compared to the computer-generated dose profile. We found excellent agreement of our measured values, even at the critical penumbra decline. For the 14 mm and 18 mm collimator and for the 11 mm collimator combination we compared the measured and calculated data at full width at half maximum. This TLD array may be useful for phantom or tissue model studies on the spatial dose distribution in confined radiation targets as used in stereotactic radiotherapy.     

On absorbed dose in narrow 60Co gamma-ray beams and dosimetry of the gamma knife

Using separate analytical functions describing primary dose, P0(dm,r), collimator scatter, Sc(r), and phantom scatter, TAR(d,r), an expression for absorbed dose in narrow 60Co gamma-ray beams is developed and each function is quantified: D(d,r) = P0(dm,r) Sc(r) TAR(d,r). The absorbed dose is calculated in beams as narrow as 0.2 cm in radius. Analytical and experimental results are compared using measured dose data for the Gamma Knife. Close agreement with experimental data is observed.     

The FE-lspd model for electron beam dosimetry

The FE-lspd model is a two-component electron beam model that distinguishes between electrons that can be described by small-angle transport theory and electrons that are too widely scattered for small-angle transport theory to be applicable. The two components are called the primary beam and the laterally scattered primary distribution (lspd). The primary beam component incorporates a simple version of the Fermi-Eyges model and dominates dose calculations at therapeutic depths. The lspd component corrects erros in the lateral spreading of the primary beam component, thereby improving the accuracy by which the FE-lspd model calculates dose distribution in blocked fields. Comparisons were made between dose profiles and central-axis depth dose distributions in small fields calculated by the FE-lspd, Fermi-Eyges and EGS4 Monte Carlo models for a 10 MeV beam in a homogeneous water phantom. The maximum difference between the dose calculated using the FE-lspd model and EGS4 Monte Carlo is about 6% at a field diameter of about 1 cm, and less than 2% for field sizes greater than 3 cm diameter. The maximum difference between the Fermi-Eyges and Monte Carlo calculations is about 18% at a field diameter of about 2.5 cm. A comparison was made with the central-axis depth dose distribution measured in water for a 3 cm diameter field in a 10 MeV clinical electron beam. The errors in the dose distribution were found to be less than 2% using the FE-lspd model but almost 18% using the Fermi-Eyges model. A comparison was also made with pencil beam profiles calculated using the second-order Fermi-Eyges transport model.     

Radiosurgical management of meningiomas

The management of residual, recurrent, or small skull base meningiomas is controversial. Stereotactic radiosurgery has emerged as an alternative treatment. We report our experience from September 1991 to August 1994 of treatment of 20 such patients [18 females -age 19-82 years, followed for 6-36 months (mean 15.5 months)] with the Leksell Gamma Knife. Nine patients were treated either with recurrent (2 patients-2 operations each) or residual tumor. Twelve patients had skull base, 3 optic nerve, 3 parasagittal, and 1 residual torcular tumor. Mean volume/diameter was 9.172 mm3/25 mm.     

Dose characteristics of in-house-built collimators for stereotactic radiotherapy with a linear accelerator

Dose characteristics of a stereotactic radiotherapy unit based on a standard Varian Clinac 4/100 4 MV linear accelerator, in-house-built Lipowitz collimators and the SMART stereotactic radiotherapy treatment planning software have been determined. Beam collimation is constituted from the standard collimators of the linear accelerator and a tertiary collimation consisting of a replaceable divergent Lipowitz collimator. Four collimators with isocentre diameters of 15, 25, 35 and 45 mm, respectively, were constructed. Beam characteristics were measured in air, acrylic or water with ionization chamber, photon diode, electron diode, diamond detector and film. Monte Carlo simulation was also applied. The radiation leakage under the collimators was less than 1% at 50 mm depth in water. Specific beam characteristics for each collimator were imported to SMART and dose planning with five non-coplanar converging 140 degrees arcs separated by 36 degrees angles was performed for treatment of a RANDO phantom. Dose verification was made with TLD and radiochromic film. The in-house-built collimators were found to be suitable for stereotactic radiotherapy and patient treatments with this system are in progress.     

Changes with age in the modulation of natural killer activity of murine leukocytes by gastrin-releasing peptide, neuropeptide Y and sulfated cholecystokinin octapeptide

From Jan. 1993 to Sept. 1995 23 patients suffering from brain metastases from renal cell carcinoma were treated with the Leksell Gamma Knife at the University of Vienna. At the time of diagnosis 13 patients had single and 10 patients presented with multiple metastatic lesions with a total of 44 metastases in MRI scans. Median tumour volume was 5500 cmm (range 100-24000 cmm). Predominant neurological symptoms and signs were different forms of hemiparesis, focal and generalized seizures, cognitive deficit, headache, dizziness, ataxia and CN XII paresis. Fourteen patients received Gamma Knife Radiosurgery (GKRS) with a median dose of 22 Gy (range 8-30 Gy) at the tumour margin. Nine patients underwent a combined treatment of a radiosurgical boost with a median dose of 18 Gy (range 10-22 Gy) at the tumour margin followed by Whole Brain Radiotherapy (total dose 30 Gy/2 weeks). In 20 patients tumour volume reduction up to 30% of the primary tumour volume was found after 4 weeks, evaluated on CT or MRI. A total remission was seen in 4 cases 3 months after GKRS. We achieved a local tumour control of 96%. Rapid neurological improvement after GKRS was seen in 17 patients. The median survival time was 11 months; the one-year actual survival in this unselected group was 48%. Five long term survivors were still alive, 18 patients had subsequently died, 15 of them of general tumour progression. GKRS induces a significant tumour remission accompanied by rapid neurological improvement and therefore provides the opportunity for extended high quality survival. Neither local tumour control was improved nor CNS relapse free survival was prolonged significantly by additional WBRT.     

Benchmarking the MCNP code for Monte Carlo modelling of an in vivo neutron activation analysis system

The Monte Carlo computer code MCNP (version 4A) has been used to develop a personal computer-based model of the Swansea in vivo neutron activation analysis (IVNAA) system. The model included specification of the neutron source (252Cf), collimators, reflectors and shielding. The MCNP model was 'benchmarked' against fast neutron and thermal neutron fluence data obtained experimentally from the IVNAA system. The Swansea system allows two irradiation geometries using 'short' and 'long' collimators, which provide alternative dose rates for IVNAA. The data presented here relate to the short collimator, although results of similar accuracy were obtained using the long collimator. The fast neutron fluence was measured in air at a series of depths inside the collimator. The measurements agreed with the MCNP simulation within the statistical uncertainty (5-10%) of the calculations. The thermal neutron fluence was measured and calculated inside the cuboidal water phantom. The depth of maximum thermal fluence was 3.2 cm (measured) and 3.0 cm (calculated). The width of the 50% thermal fluence level across the phantom at its mid-depth was found to be the same by both MCNP and experiment. This benchmarking exercise has given us a high degree of confidence in MCNP as a tool for the design of IVNAA systems.     

Monte Carlo simulation of a miniature, radiosurgery x-ray tube using the ITS 3.0 coupled electron-photon transport code

A miniature, interstitial x-ray generator has recently been developed and is currently undergoing clinical trials for the treatment of brain tumors. The maximum photon energy from this x-ray tube is 50 keV, although most of the initial testing has been carried out at 40 keV. Dose rates of up to 2 Gy/min in a water phantom at a distance of 10 mm from the tube tip are produced. In this paper we describe the modeling and simulation of x-ray production from this device using the ITS 3.0 Monte Carlo code. Verification of the simulation of x-ray production in the device was carried out by comparing predictions of spatial photon distribution, energy spectrum, and dose versus depth in water with experimentally obtained measurements. Agreement between the simulated results and experimental measurements was fairly good when comparing the angular distribution of photons emitted from the x-ray tube and very good when comparing dose rate versus depth in a water phantom. Discrepancies observed when comparing the calculated and measured estimates of characteristic line radiation were reduced by incorporation of a modification to the ITS code. Possible causes of the remaining discrepancy in bremsstrahlung intensity are discussed.     

Monte Carlo methods for the in vivo analysis of cisplatin using X-ray fluorescence

A Monte Carlo method has been used to model the measurement of cisplatin uptake with in vivo X-ray fluorescence. A user-code has been written for the EGS4 Monte Carlo system that incorporates linear polarisation and multiple element fluorescence extensions. The yield of fluorescent photons to the mainly Compton scattered background is computed for our detector arrangement. The detector consists of a mutually orthogonal arrangement of X-ray tube, aluminium polariser and high purity germanium scintillation detector. The influence of tube voltage on the minimum detectable concentration is modelled for 100 through 150 kVp X-radiation. The code is able to predict absorbed dose to the patient which will influence the optimal choice of tube voltage. The influence of alterations to collimator design and scatterer construction can also be examined. A minimum detectable concentration of 50 ppm is determined from measurements with a 115 kVp X-ray source and a 615 ppm cisplatin sample in a water phantom.     

Combined experimental and Monte Carlo verification of 137Cs brachytherapy plans for vaginal applicators

Dose rates in a phantom around a shielded and an unshielded vaginal applicator containing Selectron low-dose-rate 137Cs sources were determined by experiment and Monte Carlo simulation. Measurements were performed with thermoluminescent dosimeters in a white polystyrene phantom using an experimental protocol geared for precision. Calculations for the same set-up were done using a version of the EGS4 Monte Carlo code system modified for brachytherapy applications into which a new combinatorial geometry package developed by Bielajew was recently incorporated. Measured dose rates agree with Monte Carlo estimates to within 5% (1 SD) for the unshielded applicator, while highlighting some experimental uncertainties for the shielded applicator. Monte Carlo calculations were also done to determine a value for the effective transmission of the shield required for clinical treatment planning, and to estimate the dose rate in water at points in axial and sagittal planes transecting the shielded applicator. Comparison with dose rates generated by the planning system indicates that agreement is better than 5% (1 SD) at most positions. The precision thermoluminescent dosimetry protocol and modified Monte Carlo code are effective complementary tools for brachytherapy applicator dosimetry.     

Validation of Monte Carlo dose calculations near 125I sources in the presence of bounded heterogeneities

PURPOSE: Dose distributions around low energy (< 60 keV) brachytherapy sources, such as 125I, are known to be very sensitive to changes in tissue composition. Available 125I dosimetry data describe the effects of replacing the entire water medium by heterogeneous material. This work extends our knowledge of tissue heterogeneity effects to the domain of bounded tissue heterogeneities, simulating clinical situations. Our goals are three-fold: (a) to experimentally characterize the variation of dose rate as a function of location and dimensions of the heterogeneity, (b) to confirm the accuracy of Monte Carlo dose calculation methods in the presence of bounded tissue heterogeneities, and (c) to use the Monte Carlo method to characterize the dependence of heterogeneity correction factors (HCF) on the irradiation geometry. METHODS AND MATERIALS: Thermoluminescent dosimeters (TLD) were used to measure the deviations from the homogeneous dose distribution of an 125I seed due to cylindrical tissue heterogeneities. A solid water phantom was machined accurately to accommodate the long axis of the heterogeneous cylinder in the transverse plane of a 125I source. Profiles were obtained perpendicular to and along the cylinder axis, in the region downstream of the heterogeneity. Measurements were repeated at the corresponding points in homogeneous solid water. The measured heterogeneity correction factor (HCF) was defined as the ratio of the detector reading in the heterogeneous medium to that in the homogeneous medium at that point. The same ratio was simulated by a Monte Carlo photon transport (MCPT) code, using accurate modeling of the source, phantom, and detector geometry. In addition, Monte Carlo-based parametric studies were performed to identify the dependence of HCF on heterogeneity dimensions and distance from the source. RESULTS: Measured and calculated HCFs reveal excellent agreement (< or = 5% average) over a wide range of materials and geometries. HCFs downstream of 20 mm diameter by 10 mm thick hard bone cylinders vary from 0.12 to 0.30 with respect to distance, while for an inner bone cylinder of the same dimension, it varies from 0.72 to 0.83. For 6 mm diameter by 10 mm thick hard bone and inner bone cylinders, HCF varies 0.27-0.58 and 0.77-0.88, respectively. For lucite, fat, and air, the dependence of HCF on the 3D irradiation geometry was much less pronounced. CONCLUSION: Monte Carlo simulation is a powerful, convenient, and accurate tool for investigating the long neglected area of tissue composition heterogeneity corrections. Simple one dimensional dose calculation models that depend only on the heterogeneity thickness cannot accurately characterize 125I dose distributions in the presence of bone-like heterogeneities.     

Stereotactic radiosurgery for tentorial meningiomas

Radical microsurgical resection is the procedure of choice for tentorial meningiomas. Despite advances in microsurgery, tentorial meningiomas continue to challenge surgeons and patients. To evaluate the response of tentorial meningiomas, we evaluated 41 patients who had Gamma knife stereotactic radiosurgery during a 9 year period. Patient age varied from 32 to 79 years. Headache, trigeminal neuralgia, or facial paraesthesia were the most common presenting symptoms. Sensory deficits in the distribution of the trigeminal nerve were the most common finding. Eighteen patients (44%) had undergone between 1 and 5 (mean, 1.9) resections prior to radiosurgery; 23 had tumors diagnosed by neuroimaging. The average tumor diameter in this series was 20 mm. The maximum tumor dose varied from 24 to 40 Gy (mean, 30.5 Gy), and the tumor margin dose varied from 12 to 20 Gy (mean, 15.3 Gy). During the average follow-up interval of 3 years (range, 1-8 years), 19 patients had clinical improvement, 20 remained stable, and 2 patients deteriorated. Follow-up imaging showed a reduction in tumor size in 18 patients, no further tumor growth in 22, and an increase in tumor size in one (overall tumor control rate of 98%). Stereotactic radiosurgery using the Gamma Knife was a safe and effective primary or adjuvant treatment for patients with tentorial meningiomas.     

Experimental and Monte Carlo dosimetry of the Henschke applicator for high dose-rate 192Ir remote afterloading

We have performed extensive computational and experimental dosimetry of the Henschke applicator with respect to high dose-rate 192Ir brachytherapy using a GAMMAMED remote afterloader. Our goal was to generate clinically useful two- and three-dimensional look-up tables. Dose measurements of the Henschke applicator involved using TLD chips placed in a polystyrene phantom. Monte Carlo simulations were performed using the MCNP code. The computational models included the detailed geometry of 192Ir source, tandem tube, and shielded ovoid. The measured dose rates were corrected for the dependence of TLD sensitivity on the distance of measurement points from the source. Transit dose delivered during source extension to and retraction from a given dwell position was estimated by Monte Carlo simulations, and a correction was applied to the experimental values. For the applicator tandem, the ratio of dose rates obtained by MCNP to those measured by TLD chips ranges from 0.92 to 1.10 with an average of 0.98 and a standard deviation of 0.02. The measured and calculated dose rates at 1 cm on the transverse axis are 1.10 cGy U-1 h-1. For the shielded ovoid, the ratio ranges from 0.88 to 1.16 with an average of 1.00 and a standard deviation of 0.07. Causes of the discrepancy between the Monte Carlo and TLD results were identified. We found that the combined uncertainty of measured dose rates due to these causes is 5.6% for the applicator tandem and 8.4% for the shielded ovoid. Therefore, the results of the Monte Carlo simulation are considered to have been validated by the measurements within the uncertainty involved in the calculation and measurements.     

Signal profile measurements for evaluation of the volume-selection performance of ISIS

High-resolution signal profiles obtained with a test phantom were used in this study to evaluate the volume-selection performance of an implementation of ISIS (Image Selected In vivo Spectroscopy). The phantom simulated the brain with regard to volume and loading of coil. A remotely controlled, movable signal source inside the phantom was filled with orthophosphoric acid. Signal profiles of the volume of interest (VOI) were measured in three perpendicular directions. Special interest was focused on the transition zones, the position of the profiles, and the effects of off-resonance and T1 smearing. The transition zones were on average 5.6 mm wide and the full width at half maximum (FWHM) was 35 mm for a VOI of 40 x 40 x 40 mm3. The positions of the centre of the signal profiles were x = 3.2, y = -0.7 and z = 3.3 mm off-centre. The deviation of the volume position could be explained by off-resonance effects during imaging and spectroscopy. These data illustrate the importance of detailed knowledge of the volume-selection performance when attempting precision measurements using image-guided in vivo MRS.     

Design and dosimetric characteristics of a high dose rate remotely afterloaded endocavitary applicator system

PURPOSE: An applicator is described for endocavitary treatment of rectal cancers using a high dose rate (HDR) remote afterloading system with a single high-intensity 192Ir source as an alternative to the 50 kVp x-ray therapy contact unit most frequently used in this application. METHODS AND MATERIALS: The applicator consists of a tungsten-alloy collimator with a 45 degree beveled end, placed in a protoscope with an elliptical cross-section. The resultant 3 cm diameter circular treatment aperture, located in the beveled face of the proctoscope, is irradiated by circular array of dwell positions located about 6.5 mm from the applicator surface. This beveled end allows patients with posterior wall tumors to be treated in the dorsal lithotomy position. The dose-rate distributions about the applicator were determined using a combination of thermoluminescent dosimetry (TLD-100 detectors) and radiochromic film dose measurement techniques along with Monte Carlo dosimetry calculations. TLD-100 (3 x 3 x 0.9 mm3 chips) measurements were used to measure the distribution of dose over the proctoscope surface as well as the central axis dose-rate distribution. Relative radiochromic film measurements were used to measure off-axis ratios (flatness and penumbra width) within the treatment aperture. These data were combined with Monte Carlo simulation results to obtain the final dose distribution. RESULTS: The tungsten collimator successfully limits the dose to the tissue in contact with the proctoscope walls to less than 12% of the prescribed dose. These results indicate that the HDR applicator system has slightly more penetrating depth-dose characteristics than the most widely used contact therapy x-ray machine. Flatness characteristics of the two treatment delivery systems are comparable, although the HDR endocavitary applicator has a significantly wider penumbra. Finally, the HDR applicator has a lower surface dose rate (1.5-4 Gy/min of dwell time) compared to 9-10 Gy/min for the x-ray unit. CONCLUSIONS: An applicator system has been developed for endocavitary treatment of early stage rectal carcinoma that uses a single-stepping source HDR remote afterloading system as a radiation source. The advantages of the HDR-based system over x-ray therapy contact units currently used in this clinical application are (a) enhanced flexibility in applicator design and (b) widespread availability of single-stepping source HDR remote afterloading systems.     

Dose perturbation caused by high-density inhomogeneities in small beams in stereotactic radiosurgery

The influence of high-density tissue heterogeneities in small-diameter beams used in stereotactic radiosurgery has been investigated. Dose perturbation immediately behind aluminium sheets, used to simulate a high-density tissue inhomogeneity such as bone, was studied in a solid water phantom. Dose reduction factors (DRFs), which are the ratios of the dose in the presence of the inhomogeneity to dose in a uniform density solid water phantom, were measured with a diamond detector for three thicknesses of aluminium. DRFs exhibit dependence on both the inhomogeneity thickness and the beam diameter. The DRF decreases with inhomogeneity thickness. The DRF initially decreases with increase in the beam diameter from 12.5 to 25 mm. For fields greater than 25 mm, the DRFs are nearly constant. The commonly used algorithms such as the TAR ratio method underestimate the magnitude of the measured effect. A good agreement between these measurements and Monte Carlo calculations is obtained. The influence of the high-density inhomogeneity on the tissue maximum ratio (TMR) was also measured with the inhomogeneity at a fixed depth dmax from the entrance surface. The TMR is reduced for all detector-inhomogeneity distances investigated. The dose build-up phenomenon observed in the presence of low-density air inhomogeneity is absent in the presence of a high-density inhomogeneity. The beam width (defined by 50% dose points) immediately beyond the inhomogeneity is unaffected by the high-density inhomogeneity. However, the 90%-10% and 80%-20% dose penumbra widths and the dose outside the beam edge (beyond the 50% dose point) are reduced. This reduction in dose outside the beam edge is caused by the reduced range of the secondary radiation (photons and electrons) in the high-density medium.     

The application of correlated sampling to the computation of electron beam dose distributions in heterogeneous phantoms using the Monte Carlo method

Although the Monte Carlo method is capable of computing the dose distribution in heterogeneous phantoms directly, there are some advantages to computing a heterogeneity correction factor. If this approach is adopted there are savings in time using correlated sampling. This technique forces histories to have the same energy, position, direction and random number seed as incident on both the heterogeneous and homogeneous water phantom. This ensures that a history that has, by chance, travelled through only water in the heterogeneous phantom will have the same path as it would have through the homogeneous phantom, resulting in a reduced variance when a ratio of heterogeneous dose to homogeneous dose is formed. Metrics to describe the distributions of uncertainty, efficiency, and degree of correlation are defined. EGS4 Monte Carlo calculation of the dose distribution from a 20 MeV electron beam on water phantoms containing aluminum or air slab heterogeneities illustrate that this technique is the most efficient when the heterogeneity is deep within the phantom, but that improved efficiency can be realized even when the heterogeneity is at or near the surface. This is because some correlation between the two histories is retained despite passage through the heterogeneity.     

Monte Carlo calculations of electron beam output factors for a medical linear accelerator

The purpose of this study was to investigate the application of the Monte Carlo technique to the calculation and analysis of output factors for electron beams used in radiotherapy. The code EGS4/BEAM was used to obtain phase-space files for 6, 12 and 20 MeV clinical electron beams from a scattering-foil linac (Varian Clinac 2100C) for a clinically representative range of applicator and square or rectangular insert combinations. The source-to-surface distance used was 100 cm. The field sizes ranged from 1 x 1 cm2 to 20 x 20 cm2. These phase-space files were analysed to study the intrinsic beam characteristics and used as source input for relative dose and output factor computations in homogeneous water phantoms using the code EGS4/DOSXYZ. The calculated relative central-axis depth-dose and transverse dose profiles at various depths of clinical interest agreed with the corresponding measured dose profiles to within 2% of the maximum dose. Calculated output factors for the fields studied agreed with measured output factors to about 2%. This demonstrated that for the Varian Clinac 2100C linear accelerator, electron beam dose calculations in homogeneous water phantoms can be performed accurately at the 2% level using Monte Carlo simulations.