Energy spectra, angular spread, fluence profiles and dose distributions of 6 and 18 MV photon beams: results of monte carlo simulations for a varian 2100EX accelerator

The purpose of this study is to provide detailed characteristics of incident photon beams for different field sizes and beam energies. This information is critical to the future development of accurate treatment planning systems. It also enhances our knowledge of radiotherapy photon beams. The EGS4 Monte Carlo code, BEAM, has been used to simulate 6 and 18 MV photon beams from a Varian Clinac-2100EX accelerator. A simulated realistic beam is stored in a phase space data file, which contains details of each particle's complete history including where it has been and where it has interacted. The phase space files are analysed to obtain energy spectra, angular distribution, fluence profile and mean energy profiles at the phantom surface for particles separated according to their charge and history. The accuracy of a simulated beam is validated by the excellent agreement between the Monte Carlo calculated and measured dose distributions. Measured depth-dose curves are obtained from depth-ionization curves by accounting for newly introduced chamber fluence corrections and the stopping-power ratios for realistic beams. The study presents calculated depth-dose components from different particles as well as calculated surface dose and contribution from different particles to surface dose across the field. It is shown that the increase of surface dose with the increase of the field size is mainly due to the increase of incident contaminant charged particles. At 6 MV, the incident charged particles contribute 7% to 21% of maximum dose at the surface when the field size increases from 10 x 10 to 40 x 40 cm2. At 18 MV, their contributions are up to 11% and 29% of maximum dose at the surface for 10 x 10 cm2 and 40 x 40 cm2 fields respectively. However, the fluence of these incident charged particles is less than 1% of incident photon fluence in all cases.     

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.     

Comparison of measured and Monte Carlo calculated dose distributions in inhomogeneous phantoms in clinical electron beams

Calculations of dose distributions in heterogeneous phantoms in clinical electron beams, carried out using the fast voxel Monte Carlo (MC) system XVMC and the conventional MC code EGSnrc, were compared with measurements. Irradiations were performed using the 9 MeV and 15 MeV beams from a Varian Clinac-18 accelerator with a 10 x 10 cm2 applicator and an SSD of 100 cm. Depth doses were measured with thermoluminescent dosimetry techniques (TLD 700) in phantoms consisting of slabs of Solid Water (SW) and bone and slabs of SW and lung tissue-equivalent materials. Lateral profiles in water were measured using an electron diode at different depths behind one and two immersed aluminium rods. The accelerator was modelled using the EGS4/BEAM system and optimized phase-space files were used as input to the EGSnrc and the XVMC calculations. Also, for the XVMC, an experiment-based beam model was used. All measurements were corrected by the EGSnrc-calculated stopping power ratios. Overall, there is excellent agreement between the corrected experimental and the two MC dose distributions. Small remaining discrepancies may be due to the non-equivalence between physical and simulated tissue-equivalent materials and to detector fluence perturbation effect correction factors that were calculated for the 9 MeV beam at selected depths in the heterogeneous phantoms.     

Monte Carlo evaluation of a photon pencil kernel algorithm applied to fast neutron therapy treatment planning

When dedicated software is lacking, treatment planning for fast neutron therapy is sometimes performed using dose calculation algorithms designed for photon beam therapy. In this work Monte Carlo derived neutron pencil kernels in water were parametrized using the photon dose algorithm implemented in the Nucletron TMS (treatment management system) treatment planning system. A rectangular fast-neutron fluence spectrum with energies 0-40 MeV (resembling a polyethylene filtered p(41)+Be spectrum) was used. Central axis depth doses and lateral dose distributions were calculated and compared with the corresponding dose distributions from Monte Carlo calculations for homogeneous water and heterogeneous slab phantoms. All absorbed doses were normalized to the reference dose at 10 cm depth for a field of radius 5.6 cm in a 30 x 40 x 20 cm3 water test phantom. Agreement to within 7% was found in both the lateral and the depth dose distributions. The deviations could be explained as due to differences in size between the test phantom and that used in deriving the pencil kernel (radius 200 cm, thickness 50 cm). In the heterogeneous phantom, the TMS, with a directly applied neutron pencil kernel, and Monte Carlo calculated absorbed doses agree approximately for muscle but show large deviations for media such as adipose or bone. For the latter media, agreement was substantially improved by correcting the absorbed doses calculated in TMS with the neutron kerma factor ratio and the stopping power ratio between tissue and water. The multipurpose Monte Carlo code FLUKA was used both in calculating the pencil kernel and in direct calculations of absorbed dose in the phantom.     

A virtual source model for Monte Carlo modeling of arbitrary intensity distributions

A photon virtual source model was developed for simulating arbitrary, external beam, intensity distributions using the Monte Carlo method. The source model consists of a photon fluence grid composed of a matrix of square elements, located 25-cm downstream from the linear accelerator target. Each particle originating from the fluence map is characterized by the seven phase space parameters, position (x, y, z), direction (u, v, w), and energy. The map was reconstructed from fluence and energy spectra acquired by modeling components of the linear accelerator treatment head using the Monte Carlo code MCNP4B. The effect of contaminant electrons is accounted for by the use of a sub-source derived from a phase-space simulation of a 25-MV linac treatment head using the code BEAM. The BEAM sub-source was incorporated into the MCNP4B phase-space model and is sampled using a field-size dependent sampling ratio. A Gaussian blurring kernel is convolved with the photon fluence map to account for the finite focal spot size and scattering effects from structures such as the flattening filter and MLC leaves. Depth dose and profile source calculations for 6-MV and 25-MV photon beams, for 5 x 5 cm2, 10 x 10 cm2, and 15 x 15 cm2 field sizes, are in good agreement with measurement and are well within acceptability criteria suggested by the AAPM Task Group Report No. 53. Irregular field calculations compared with film measurement and with a 3-D pencil beam algorithm show that the source model is capable of accurately simulating arbitrary MLC fields.     

Comparison of the Batho, ETAR and Monte Carlo dose calculation methods in CT based patient models

This paper shows the contribution that Monte Carlo methods make in regard to dose distribution calculations in CT based patient models and the role it plays as a gold standard to evaluate other dose calculation algorithms. The EGS4 based BEAM code was used to construct a generic 8 MV accelerator to obtain a series of x-ray field sources. These were used in the EGS4 based DOSXYZ code to generate beam data in a mathematical water phantom to set up a beam model in a commercial treatment planning system (TPS), CADPLAN V.2.7.9. Dose distributions were calculated with the Batho and ETAR inhomogeneity correction algorithms in head/sinus, lung, and prostate patient models for 2 x 2, 5 x 5, and 10 X 10 cm2 open x-ray beams. Corresponding dose distributions were calculated with DOSXYZ that were used as a benchmark. The dose comparisons are expressed in terms of 2D isodose distributions, percentage depth dose data, and dose difference volume histograms (DDVH's). Results indicated that the Batho and ETAR methods contained inaccuracies of 20%-70% in the maxillary sinus region in the head model. Large lung inhomogeneities irradiated with small fields gave rise to absorbed dose deviations of 10%-20%. It is shown for a 10 x 10 cm2 field that DOSXYZ models lateral scatter in lung that is not present in the Batho and ETAR methods. The ETAR and Batho methods are accurate within 3% in a prostate model. We showed how the performance of these inhomogeneity correction methods can be understood in realistic patient models using validated Monte Carlo codes such as BEAM and DOSXYZ.     

Benchmarking of the dose planning method (DPM) Monte Carlo code using electron beams from a racetrack microtron

A comprehensive set of measurements and calculations has been conducted to investigate the accuracy of the Dose Planning Method (DPM) Monte Carlo code for dose calculations from 10 and 50 MeV scanned electron beams produced from a racetrack microtron. Central axis depth dose measurements and a series of profile scans at various depths were acquired in a water phantom using a Scanditronix type RK ion chamber. Source spatial distributions for the Monte Carlo calculations were reconstructed from in-air ion chamber measurements carried out across the two-dimensional beam profile at 100 cm downstream from the source. The in-air spatial distributions were found to have full width at half maximum of 4.7 and 1.3 cm, at 100 cm from the source, for the 10 and 50 MeV beams, respectively. Energy spectra for the 10 and 50 MeV beams were determined by simulating the components of the microtron treatment head using the code MCNP4B. DPM calculations are on average within +/- 2% agreement with measurement for all depth dose and profile comparisons conducted in this study. The accuracy of the DPM code illustrated in this work suggests that DPM may be used as a valuable tool for electron beam dose calculations.     

Interface perturbation effects in high-energy electron beams

Near interfaces between two different media exposed to high-energy electron beams substantial dose and fluence perturbations due to backscatter can be observed. In this work, dose and fluence perturbations were studied for 4-19 MeV electron beams at backscatter interfaces of polystyrene, graphite, water, aluminium and lead. Measurements of relative dose using an NPL-designed thin-window plane-parallel ion chamber and a Markus ion chamber were performed to determine the effect of different interface materials and thicknesses. Results of Monte Carlo simulations with the EGSnrc code, including models of the ion chambers, were found to be in excellent agreement with the measurements. The well-documented increasing dose perturbation with increasing effective atomic number of the backscatter material and decreasing electron beam energy was confirmed. Simulations in a simplified slab geometry showed that, despite the decrease of average electron energy with depth in water, the dose perturbations decrease with increasing depth of the interface in water for all the materials in the study. This was ascribed to the change of the electron angular distribution with depth which has a different effect in water and in the presence of a high-Z interface. Electron fluence perturbations near a lead/water interface were found to cause small differences in unrestricted mass collision stopping power ratios, water to air. Effects of bremsstrahlung photons, characteristic photons and positrons from the backscattering material were found to be insignificant for electron interface dosimetry. When comparing simulations using EGSnrc and the older version of the same code, EGS4, underestimations of the dose perturbation effects of up to 7% were found when using the latter code to simulate 4 MeV electrons irradiating a lead/water interface. It is concluded that EGSnrc is a highly suitable tool for electron interface dosimetry studies.     

Clinical implementation of a Monte Carlo treatment planning system.

The purpose of this study was to implement the Monte Carlo method for clinical radiotherapy dose calculations. We used the EGS4/BEAM code to obtain the phase-space data for 6-20 MeV electron beams and 4, 6, and 15 MV photon beams for Varian Clinac 1800, 2100C, and 2300CD accelerators. A multiple-source model was used to reconstruct the phase-space data for both electron and photon beams, which retained the accuracy of the Monte Carlo beam data. The multiple-source model reduced the phase-space data storage requirement by a factor of 1000 and the accelerator simulation time by a factor of 10 or more. Agreement within 2% was achieved between the Monte Carlo calculations and measurements of the dose distributions in homogeneous and heterogeneous phantoms for various field sizes, source-surface distances, and beam modulations. The Monte Carlo calculated electron output factors were within 2% of the measured values for various treatment fields while the heterogeneity correction factors for various lung and bone phantoms were within 1% for photon beams and within 2% for electron beams. The EGS4/DOSXYZ Monte Carlo code was used for phantom and patient dose calculations. The results were compared to the dose distributions produced by a conventional treatment planning system and an intensity-modulated radiotherapy inverse-planning system. Significant differences (>5% in dose and >5 mm shift in isodose lines) were found between Monte Carlo calculations and the analytical calculations implemented in the commercial systems. Treatment sites showing the largest dose differences were for head and neck, lung, and breast cases.     

Secondary electron fluence perturbation by high-Z interfaces in clinical proton beams: a Monte Carlo study.

Fluence perturbation of secondary electrons from clinical proton beams (50-250 MeV) by thin high-Z planar interfaces was studied with Monte Carlo simulations. Starting from monoenergetic proton pencil beams, proton depth doses and proton fluence spectra were calculated, both in homogeneous water and near thin high-Z interfaces by using the proton transport Monte Carlo code PTRAN. This code was modified extensively to enable modelling of proton transport in non-homogeneous geometries. From the proton fluence spectra in water and in the interface materials, electron generation spectra were calculated analytically and were then used as input for an electron transport calculation with the Monte Carlo code EGS4/PRESTAII to obtain electron doses and electron fluence spectra. The interface materials used in the study were graphite, Al, Ti, Cu, Sn and Au. We found significant electron fluence perturbations on both sides of the planar interfaces, resulting in an electron dose increase upstream and a decrease downstream from the interfaces, with the magnitude of the effect depending strongly on the atomic number of the interface. For the most extreme case studied, 250 MeV protons and a gold interface, we obtained an electron dose increase of 41% upstream of the interface and a decrease of 15% downstream with both perturbations having a spatial extent of about 700 microm. The total dose perturbation due to this effect amounts to a 5% increase upstream and a 2% decrease downstream. A detailed analysis of dose and fluence perturbation is presented for a wide range of materials and proton energies.     

Model for calculating depth dose distributions for broad electron beams

Central axis electron beam depth dose distributions can be transformed by replacing dose by fluence and depth by a measure of angular dispersion. This transformation was applied to a set of broad beam central axis depth dose distributions calculated by Monte Carlo code for beams with initial energies ranging from 1 to 60 MeV in homogeneous media of water, aluminum, and copper. The resulting fluence distributions belong to a family of curves that can be parametrized by a single-valued function of initial beam energy and medium and can be used to calculate fluence distributions for accelerators. Fluence curves can be easily transformed to depth dose curves.     

Commissioning stereotactic radiosurgery beams using both experimental and theoretical methods

The purpose of this investigation is to study the feasibility of using an alternative method to commission stereotactic radiosurgery beams shaped by micro multi-leaf collimators by using Monte Carlo simulations to obtain beam characteristics of small photon beams, such as incident beam particle fluence and energy distributions, scatter ratios, depth-dose curves and dose profiles where measurements are impossible or difficult. Ionization chambers and diode detectors with different sensitive volumes were used in the measurements in a water phantom and the Monte Carlo codes BEAMnrc/DOSXYZnrc were used in the simulation. The Monte Carlo calculated data were benchmarked against measured data for photon beams with energies of 6 MV and 10 MV produced from a Varian Trilogy accelerator. The measured scatter ratios and cross-beam dose profiles for very small fields are shown to be not only dependent on the size of the sensitive volume of the detector used but also on the type of detectors. It is known that the response of some detectors changes at small field sizes. Excellent agreement was seen between scatter ratios measured with a small ion chamber and those calculated from Monte Carlo simulations. The values of scatter ratios, for field sizes from 6 x 6 mm2 to 98 x 98 mm2, range from 0.67 to 1.0 and from 0.59 to 1.0 for 6 and 10 MV, respectively. The Monte Carlo calculations predicted that the incident beam particle fluence is strongly affected by the X-Y-jaw openings, especially for small fields due to the finite size of the radiation source. Our measurement confirmed this prediction. This study demonstrates that Monte Carlo calculations not only provide accurate dose distributions for small fields where measurements are difficult but also provide additional beam characteristics that cannot be obtained from experimental methods. Detailed beam characteristics such as incident photon fluence distribution, energy spectra, including composition of primary and scattered photons, can be independently used in dose calculation models and to improve the accuracy of measurements with detectors with an energy-dependent response. Furthermore, when there are discrepancies between results measured with different detectors, the Monte Carlo calculated values can indicate the most correct result. The data set presented in this study can be used as a reference in commissioning stereotactic radiosurgery beams shaped by a BrainLAB m3 on a Varian 2100EX or 600C accelerator.     

Monte Carlo calculation of nine megavoltage photon beam spectra using the BEAM code

A recent paper analyzed the sensitivity to various simulation parameters of the Monte Carlo simulations of nine beams from three major manufacturers of commercial medical linear accelerators, ranging in energy from 4-25 MV. In this work the nine models are used: to calculate photon energy spectra and average energy distributions and compare them to those published by Mohan et al. [Med. Phys. 12, 592-597 (1985)]; to separate the spectra into primary and scatter components from the primary collimator, the flattening filter and the adjustable collimators; and to calculate the contaminant-electron fluence spectra and the electron contribution to the depth-dose curves. Notwithstanding the better precision of the calculated spectra, they are similar to those calculated by Mohan et al. The three photon spectra at 6 MV from the machines of three different manufacturers show differences in their shapes as well as in the efficiency of bremsstrahlung production in the corresponding target and filter combinations. The contribution of direct photons to the photon energy fluence in a 10 x 10 field varies between 92% and 97%, where the primary collimator contributes between 0.6% and 3.4% and the flattening filter contributes between 0.6% and 4.5% to the head-scatter energy fluence. The fluence of the contaminant electrons at 100 cm varies between 5 x 10(-9) and 2.4 x 10(-7) cm(-2) per incident electron on target, and the corresponding spectrum for each beam is relatively invariant inside a 10 x 10 cm2 field. On the surface the dose from electron contamination varies between 5.7% and 11% of maximum dose and, at the depth of maximum dose, between 0.16% and 2.5% of maximum dose. The photon component of the percentage depth-dose at 10 cm depth is compared with the general formula provided by AAPM's task group 51 and confirms the claimed accuracy of 2%.     

A systematic evaluation of air cavity dose perturbation in megavoltage x-ray beams

The EGS4 Monte Carlo radiation transport code was used to systematically study the dose perturbation near planar and cylindrical air cavities in a water medium irradiated by megavoltage x-ray beams. The variables of the problem included x-ray energy, cavity shape and dimension, and depth of the cavity in water. The Monte Carlo code was initially validated against published measurements and its results were found to agree within 2% with the published measurements. The study results indicate that the dose perturbation is strongly dependent on x-ray energy, field size, depth, and size of cavity in water. For example, the Monte Carlo calculations show dose reductions of 42% and 18% at 0.05 and 2 mm, respectively, beyond the air-water interface distal to the radiation source for a 3 cm thick air slab irradiated by a single 5x5 cm2 15 MV beam. The dose reductions are smaller for a parallel-opposed pair of 5x5 cm2 15 MV x-ray beams, being 21% and 11% for the same depths. The combined set of Monte Carlo calculations showed that the dose reduction near an air cavity is greater for: (a) Smaller x-ray field size, (b) higher x-ray energy, (c) larger air-cavity size, and (d) smaller depth in water where the air cavity is situated. A potential clinical application of these results to the treatment of prostate cancer is discussed.     

Bone and mucosal dosimetry in skin radiation therapy: a Monte Carlo study using kilovoltage photon and megavoltage electron beams

This study examines variations of bone and mucosal doses with variable soft tissue and bone thicknesses, mimicking the oral or nasal cavity in skin radiation therapy. Monte Carlo simulations (EGSnrc-based codes) using the clinical kilovoltage (kVp) photon and megavoltage (MeV) electron beams, and the pencil-beam algorithm (Pinnacle(3)?treatment planning system) using the MeV electron beams were performed in dose calculations. Phase-space files for the 105 and 220 kVp beams (Gulmay D3225 x-ray machine), and the 4 and 6?MeV electron beams (Varian 21 EX linear accelerator) with a field size of 5?cm diameter were generated using the BEAMnrc code, and verified using measurements. Inhomogeneous phantoms containing uniform water, bone and air layers were irradiated by the kVp photon and MeV electron beams. Relative depth, bone and mucosal doses were calculated for the uniform water and bone layers which were varied in thickness in the ranges of 0.5-2?cm and 0.2-1?cm. A uniform water layer of bolus with thickness equal to the depth of maximum dose (d(max)) of the electron beams (0.7?cm for 4 MeV and 1.5?cm for 6 MeV) was added on top of the phantom to ensure that the maximum dose was at the phantom surface. From our Monte Carlo results, the 4 and 6 MeV electron beams were found to produce insignificant bone and mucosal dose (<1%), when the uniform water layer at the phantom surface was thicker than 1.5?cm. When considering the 0.5?cm thin uniform water and bone layers, the 4 MeV electron beam deposited less bone and mucosal dose than the 6 MeV beam. Moreover, it was found that the 105 kVp beam produced more than twice the dose to bone than the 220 kVp beam when the uniform water thickness at the phantom surface was small (0.5?cm). However, the difference in bone dose enhancement between the 105 and 220 kVp beams became smaller when the thicknesses of the uniform water and bone layers in the phantom increased. Dose in the second bone layer interfacing with air was found to be higher for the 220 kVp beam than that of the 105 kVp beam, when the bone thickness was 1?cm. In this study, dose deviations of bone and mucosal layers of 18% and 17% were found between our results from Monte Carlo simulation and the pencil-beam algorithm, which overestimated the doses. Relative depth, bone and mucosal doses were studied by varying the beam nature, beam energy and thicknesses of the bone and uniform water using an inhomogeneous phantom to model the oral or nasal cavity. While the dose distribution in the pharynx region is unavailable due to the lack of a commercial treatment planning system commissioned for kVp beam planning in skin radiation therapy, our study provided an essential insight into the radiation staff to justify and estimate bone and mucosal dose.     

Electron beam modeling and commissioning for Monte Carlo treatment planning

A hybrid approach for commissioning electron beam Monte Carlo treatment planning systems has been studied. The approach is based on the assumption that accelerators of the same type have very similar electron beam characteristics and the major difference comes from the on-site tuning of the electron incident energy at the exit window. For one type of accelerator, a reference machine can be selected and simulated with the Monte Carlo method. A multiple source model can be built on the full Monte Carlo simulation of the reference beam. When commissioning electron beams from other accelerators of the same type, the energy spectra in the source model are tuned to match the measured dose distributions. A Varian Clinac 2100C accelerator was chosen as the reference machine and a four-source beam model was established based on the Monte Carlo simulations. This simplified beam model can be used to generate Monte Carlo dose distributions accurately (within 2%/2 mm compared to those calculated with full phase space data) for electron beams from the reference machine with various nominal energies, applicator sizes, and SSDs. Three electron beams were commissioned by adjusting the energy spectra in the source model. The dose distributions calculated with the adjusted source model were compared with the dose distributions calculated using the phase space data for these beams. The agreement is within 1% in most of cases and 2% in all situations. This preliminary study has shown the capability of the commissioning approach for handling large variation in the electron incident energy. The possibility of making the approach more versatile is also discussed.     

Pencil beam approach for correcting the energy dependence artifact in film dosimetry for IMRT verification

The higher sensitivity to low-energy scattered photons of radiographic film compared to water can lead to significant dosimetric error when the beam quality varies significantly within a field. Correcting for this artifact will provide greater accuracy for intensity modulated radiation therapy (IMRT) verification dosimetry. A procedure is developed for correction of the film energy-dependent response by creating a pencil beam kernel within our treatment planning system to model the film response specifically. Film kernels are obtained from EGSnrc Monte Carlo simulations of the dose distribution from a 1 mm diameter narrow beam in a model of the film placed at six depths from 1.5 to 40 cm in polystyrene and solid water phantoms. Kernels for different area phantoms (50 x 50 cm2 and 25 x 25 cm2 polystyrene and 30 x 30 cm2 solid water) are produced. The Monte Carlo calculated kernel is experimentally verified with film, ion chamber and thermoluminescent dosimetry (TLD) measurements in polystyrene irradiated by a narrow beam. The kernel is then used in convolution calculations to, predict the film response in open and IMRT fields. A 6 MV photon beam and Kodak XV2 film in a polystyrene phantom are selected to test the method as they are often used in practice and can result in large energy-dependent artifacts. The difference in dose distributions calculated with the film kernel and the water kernel is subtracted from film measurements to obtain a practically film artifact free IMRT dose distribution for the Kodak XV2 film. For the points with dose exceeding 5 cGy (11% of the peak dose) in a large modulated field and a film measurement inside a large polystyrene phantom at depth of 10 cm, the correction reduces the fraction of pixels for which the film dose deviates from dose to water by more than 5% of the mean film dose from 44% to 6%.     

A Monte Carlo investigation of fluence profiles collimated by an electron specific MLC during beam delivery for modulated electron radiation therapy

Modulated electron radiation therapy (MERT) is able to deliver conformal dose to shallow tumors while significantly reducing dose to distal structures and surrounding tissues. An electron specific multileaf collimator (eMLC) has been proposed and constructed as an effective means of delivering electron beams for MERT. The aim of this work is to apply the Monte Carlo method to investigate the fluence profiles collimated by the eMLC in order to achieve accurate beam delivery for MERT. In this work, the EGS4/BEAM code was used to simulate the eMLC collimated electron beams of 6-20 MeV generated from a Varian Clinac 2100C linear accelerator. An attempt was made to describe the fluence profiles with an analytic Sigmoid function. The function parameters were determined by the fittings of the Monte Carlo simulated fluence profiles. How the function parameters depend on the eMLC aperture size, the off-axis location, and the electron beam energy has been investigated. It has been found that the eMLC collimated fluence profiles are dependent on beam energy, while almost independent of leaf location or dimension of MLC aperture. There is little difference in the fluence profiles collimated by the leaf side and the leaf end for the straight-edged leaves. It is possible that these energy-dependent Sigmoid functions can serve as operators to account for the energy dependence of the eMLC collimated fluence profiles. These operators can be incorporated into the inverse planning algorithm to derive desired dose distributions using a set of electron beams of variable energy and field size suitable for delivery by the eMLC.     

Validation of calculations for electrons modulated with conventional photon multileaf collimators

Treating shallow tumors with a homogeneous dose while simultaneously minimizing the dose to distal critical organs remains a challenge in radiotherapy. One promising approach is modulated electron radiotherapy (MERT). Due to the scattering properties of electron beams, the commercially provided secondary and tertiary photon collimation systems are not conducive for electron beam delivery when standard source-to-surface distances are used. Also, commercial treatment planning systems may not accurately model electron-beam dose distributions when collimated without the standard applicators. However, by using the photon multileaf collimators (MLCs) to create segments to modulate electron beams, the quality of superficial tumor dose distributions may improve substantially. The purpose of this study is to develop and evaluate calculations for the narrow segments needed to modulate megavoltage electron beams using photon beam multileaf collimators. Modulated electron radiotherapy (MERT) will be performed with a conventional linear accelerator equipped with a 120 leaf MLC for 6-20 MeV electron beam energies. To provide a sharp penumbra, segments were delivered with short SSDs (70-85 cm). Segment widths (SW) ranging from 1 to 10 cm were configured for delivery and planning, using BEAMnrc Monte Carlo (MC) code, and the DOSXYZnrc MC dose calculations. Calculations were performed with voxel size of 0.2 x 0.2 x 0.1 cm3. Dosimetry validation was performed using radiographic film and micro- or parallel-plate chambers. Calculated and measured data were compared using technical computing software. Beam sharpness (penumbra) degraded with decreasing incident beam energy and field size (FS), and increasing SSD. A 70 cm SSD was found to be optimal. The PDD decreased significantly with decreasing FS. The comparisons demonstrated excellent agreement for calculations and measurements within 3%, 1 mm. This study shows that accurate calculations for MERT as delivered with existing photon MLC are feasible and allows the opportunity to take advantage of the dynamic leaf motion capabilities and control systems, to provide conformal dose distributions.