Analysis of various beamlet sizes for IMRT with 6 MV photons

Application of intensity modulated radiation therapy (IMRT) using multileaf collimation often requires the use of small beamlets to optimize the delivered radiation distribution. Small-beam dose distribution measurements were compared to dose distributions calculated using a commercial treatment planning system that models its data acquired using measurements from relatively large fields. We wanted to evaluate only the penumbra, percent depth-dose (PDD) and output model, so we avoided dose distribution features caused by rounded leaf ends and interleaf leakage by making measurements using the secondary collimators. We used a validated radiochromic film dosimetry system to measure high-resolution dose distributions of 6 MV photon beams. A commercial treatment planning system using the finite size pencil beam (FSPB) dose calculation algorithm was commissioned using measured central axis outputs from 4.0x4.0 to 40.0x40.0 cm2 beams and radiographic-film profile measurements of a 4.0x4.0 cm2 beam at twice the depth of maximum dose (dmax). Calculated dose distributions for square fields of 0.5x0.5 cm2, and 1.0x1.0 cm2, to 6.0x6.0 cm2, in 1.0x1.0 cm2, increments were compared against radiochromic film measurements taken with the film oriented parallel to the beam central axis in a water equivalent phantom. The PDD of the smaller field sizes exhibited behavior typical of small fields, namely a decrease in dmax with decreasing field size. The FSPB accurately modeled the depth-dose and central axis output for depths deeper than the nominal dmax of 1.5 cm plus 0.5 cm. The dose distribution in the build-up and penumbra regions was not accurately modeled for depths less than 2 cm, especially for the fields of 2.0x2.0 cm2 and smaller. Using the gamma function with 2 mm and 2% criteria, the dose model was shown to accurately predict the penumbra. While for single small beams the compared dose distributions passed the gamma function criteria, the clinical appropriateness of these criteria is not clear for a composite IMRT plan. Further investigation of the cumulative impact of the observed dose discrepancies is warranted. We speculate that the observed differences in the penumbra regions arise from some energy dependent artifact in the radiographic-film profiles used for commissioning. In the future, radiochromic film based commissioning might provide a more accurate data set for dose modeling.     

The accuracy of dose calculations by anisotropic analytical algorithms for stereotactic radiotherapy in nasopharyngeal carcinoma

Nasopharyngeal tumors are commonly treated with intensity-modulated radiotherapy techniques. For photon dose calculations, problems related to loss of lateral electronic equilibrium exist when small fields are used. The anisotropic analytical algorithm (AAA) implemented in Varian Eclipse was developed to replace the pencil beam convolution (PBC) algorithm for more accurate dose prediction in an inhomogeneous medium. The purpose of this study was to investigate the accuracy of the AAA for predicting interface doses for intensity-modulated stereotactic radiotherapy boost of nasopharyngeal tumors. The central axis depth dose data and dose profiles of phantoms with rectangular air cavities for small fields were measured using a 6 MV beam. In addition, the air-tissue interface doses from six different intensity-modulated stereotactic radiotherapy plans were measured in an anthropomorphic phantom. The nasopharyngeal region of the phantom was especially modified to simulate the air cavities of a typical patient. The measured data were compared to the data calculated by both the AAA and the PBC algorithm. When using single small fields in rectangular air cavity phantoms, both AAA and PBC overestimated the central axis dose at and beyond the first few millimeters of the air-water interface. Although the AAA performs better than the PBC algorithm, its calculated interface dose could still be more than three times that of the measured dose when a 2 × 2 cm(2) field was used. Testing of the algorithms using the anthropomorphic phantom showed that the maximum overestimation by the PBC algorithm was 20.7%, while that by the AAA was 8.3%. When multiple fields were used in a patient geometry, the dose prediction errors of the AAA would be substantially reduced compared with those from a single field. However, overestimation of more than 3% could still be found at some points at the air-tissue interface.     

Evaluation of Al2O3:C optically stimulated luminescence (OSL) dosimeters for passive dosimetry of high-energy photon and electron beams in radiotherapy

This article investigates the performance of Al2O3: C optically stimulated luminescence dosimeters (OSLDs) for application in radiotherapy. Central-axis depth dose curves and optically stimulated luminescence (OSL) responses were obtained in a water phantom for 6 and 18 MV photons, and for 6, 9, 12, 16, and 20 MeV electron beams from a Varian 21EX linear accelerator. Single OSL measurements could be repeated with a precision of 0.7% (one standard deviation) and the differences between absorbed doses measured with OSLDs and an ionization chamber were within +/- 1% for photon beams. Similar results were obtained for electron beams in the low-gradient region after correction for a 1.9% photon-to-electron bias. The distance-to-agreement values were of the order of 0.5-1.0 mm for electrons in high dose gradient regions. Additional investigations also demonstrated that the OSL response dependence on dose rate, field size, and irradiation temperature is less than 1% in the conditions of the present study. Regarding the beam energy/quality dependence, the relative response of the OSLD for 18 MV was (0.51 +/- 0.48)% of the response for the 6 MV photon beam. The OSLD response for the electron beams relative to the 6 MV photon beam. The OSLD response for the electron beams relative to the 6 MV photon beam was in average 1.9% higher, but this result requires further confirmation. The relative response did not seem to vary with electron energy at dmax within the experimental uncertainties (0.5% in average) and, therefore, a fixed correction factor of 1.9% eliminated the energy dependence in our experimental conditions.     

Analysis of the penumbra enlargement in lung versus the quality index of photon beams: a methodology to check the dose calculation algorithm

It is well known that considerable underdosage can occur at the edges of a tumor inside the lung because of the degradation of penumbra due to lack of lateral electronic equilibrium. Although present even at smaller energies, this phenomenon is more pronounced for higher energies. Apart from Monte Carlo calculation, most of the existing Treatment Planning Systems (TPSs) cannot deal at all, or with acceptable accuracy, with this effect. A methodology has been developed for assessing the dose calculation algorithms in the lung region where lateral electronic disequilibrium exists, based on the Quality Index (QI) of the incident beam. A phantom, consisting of layers of polystyrene and lung material, has been irradiated using photon beams of 4, 6, 15, and 20 MV. The cross-plane profiles of each beam for 5x5, 10x10, and 25x10 fields have been measured at the middle of the phantom with the use of films. The penumbra (20%-80%) and fringe (50%-90%) enlargement was measured and the ratio of the widths for the lung to that of polystyrene was defined as the Correction Factor (CF). Monte Carlo calculations in the two phantoms have also been performed for energies of 6, 15, and 20 MV. Five commercial TPS's algorithms were tested for their ability to predict the penumbra and fringe enlargement. A linear relationship has been found between the QI of the beams and the CF of the penumbra and fringe enlargement for all the examined fields. Monte Carlo calculations agree very well (less than 1% difference) with the film measurements. The CF values range between 1.1 for 4 MV (QI 0.620) and 2.28 for 20 MV (QI 0.794). Three of the tested TPS's algorithms could not predict any enlargement at all for all energies and all fields and two of them could predict the penumbra enlargement to some extent. The proposed methodology can help any user or developer to check the accuracy of its algorithm for lung cases, based on a simple phantom geometry and the QI of the incident beam. This check is very important especially when higher energies are used, as the inaccuracies in existing algorithms can lead to an incorrect choice of energy for lung treatment and consequently to a failure in tumor control.     

Dosimetric validation of the anisotropic analytical algorithm for photon dose calculation: fundamental characterization in water

In July 2005 a new algorithm was released by Varian Medical Systems for the Eclipse planning system and installed in our institute. It is the anisotropic analytical algorithm (AAA) for photon dose calculations, a convolution/superposition model for the first time implemented in a Varian planning system. It was therefore necessary to perform validation studies at different levels with a wide investigation approach. To validate the basic performances of the AAA, a detailed analysis of data computed by the AAA configuration algorithm was carried out and data were compared against measurements. To better appraise the performance of AAA and the capability of its configuration to tailor machine-specific characteristics, data obtained from the pencil beam convolution (PBC) algorithm implemented in Eclipse were also added in the comparison. Since the purpose of the paper is to address the basic performances of the AAA and of its configuration procedures, only data relative to measurements in water will be reported. Validation was carried out for three beams: 6 MV and 15 MV from a Clinac 2100C/D and 6 MV from a Clinac 6EX. Generally AAA calculations reproduced very well measured data, and small deviations were observed, on average, for all the quantities investigated for open and wedged fields. In particular, percentage depth-dose curves showed on average differences between calculation and measurement smaller than 1% or 1 mm, and computed profiles in the flattened region matched measurements with deviations smaller than 1% for all beams, field sizes, depths and wedges. Percentage differences in output factors were observed as small as 1% on average (with a range smaller than +/-2%) for all conditions. Additional tests were carried out for enhanced dynamic wedges with results comparable to previous results. The basic dosimetric validation of the AAA was therefore considered satisfactory.     

A Monte Carlo based three-dimensional dose reconstruction method derived from portal dose images

The verification of intensity-modulated radiation therapy (IMRT) is necessary for adequate quality control of the treatment. Pretreatment verification may trace the possible differences between the planned dose and the actual dose delivered to the patient. To estimate the impact of differences between planned and delivered photon beams, a three-dimensional (3-D) dose verification method has been developed that reconstructs the dose inside a phantom. The pretreatment procedure is based on portal dose images measured with an electronic portal imaging device (EPID) of the separate beams, without the phantom in the beam and a 3-D dose calculation engine based on the Monte Carlo calculation. Measured gray scale portal images are converted into portal dose images. From these images the lateral scattered dose in the EPID is subtracted and the image is converted into energy fluence. Subsequently, a phase-space distribution is sampled from the energy fluence and a 3-D dose calculation in a phantom is started based on a Monte Carlo dose engine. The reconstruction model is compared to film and ionization chamber measurements for various field sizes. The reconstruction algorithm is also tested for an IMRT plan using 10 MV photons delivered to a phantom and measured using films at several depths in the phantom. Depth dose curves for both 6 and 10 MV photons are reconstructed with a maximum error generally smaller than 1% at depths larger than the buildup region, and smaller than 2% for the off-axis profiles, excluding the penumbra region. The absolute dose values are reconstructed to within 1.5% for square field sizes ranging from 5 to 20 cm width. For the IMRT plan, the dose was reconstructed and compared to the dose distribution with film using the gamma evaluation, with a 3% and 3 mm criterion. 99% of the pixels inside the irradiated field had a gamma value smaller than one. The absolute dose at the isocenter agreed to within 1% with the dose measured with an ionization chamber. It can be concluded that our new dose reconstruction algorithm is able to reconstruct the 3-D dose distribution in phantoms with a high accuracy. This result is obtained by combining portal dose images measured prior to treatment with an accurate dose calculation engine.     

Surface dosimetry for oblique tangential photon beams: a Monte Carlo simulation study

The effect of beam obliquity on the surface relative dose profiles for the tangential photon beams was studied. The 6 and 15 MV photon beams with 4 x 4 and 10 x 10 cm2 field sizes produced by a Varian 21 EX linear accelerator were used. Phase-space models of the photon beams were created using Monte Carlo simulations based on the EGSnrc code, and were verified using film measurements. The relative dose profiles in the phantom skin, at 2 mm depth from the surface of the half-phantom geometry, or HPG, were calculated for increasing gantry angles from 270 to 280 deg clockwise. Relative dose profiles of a full phantom enclosing the whole tangential beam (full phantom geometry, or FPG) were also calculated using Monte Carlo simulation as a control for comparison. The results showed that, although the relative dose profiles in the phantom skin did not change significantly with an oblique beam using a FPG, the surface relative depth dose was increased for the HPG. In the HPG, with 6 MV photon beams and field size = 10 x 10 cm2, when the beam angle, starting from 270 deg, was increased from 1 to 3 deg, the relative depth doses in the phantom skin were increased from 68% to 79% at 10 cm depth. This increase in dose was slightly larger than the dose from 15 MV photon beams with the same field size and beam angles, where the relative depth doses in phantom skin were increased from 81% to 87% at 10 cm depth. A parameter called the percent depth dose (PDD) ratio, defined as the relative depth dose from the HPG to the relative depth dose from the FPG at a given depth along the phantom skin, was used to evaluate the effect of the phantom-air interface. It is found that the PDD ratio increased significantly when the beam angle was changed from zero to 1-3 degrees. Moreover, the PDD ratio, for a given field size, experienced a greater increase for 6 MV than for 15 MV. For the same photon beam energy, the PDD ratio increased more with a 4 x 4 cm2 field compared to 10 x 10 cm2. The results in this study will be useful for physicists and dosimetrists to predict the surface relative dose variations when using clinical tangential-like photon beams in radiation therapy.     

Mixing intensity modulated electron and photon beams: combining a steep dose fall-off at depth with sharp and depth-independent penumbras and flat beam profiles.

For application in radiotherapy, intensity modulated high-energy electron and photon beams were mixed to create dose distributions that feature: (a) a steep dose fall-off at larger depths, similar to pure electron beams, (b) flat beam profiles and sharp and depth-independent beam penumbras, as in photon beams, and (c) a selectable skin dose that is lower than for pure electron beams. To determine the required electron and photon beam fluence profiles, an inverse treatment planning algorithm was used. Mixed beams were realized at a MM50 racetrack microtron (Scanditronix Medical AB, Sweden), and evaluated by the dose distributions measured in a water phantom. The multileaf collimator of the MM50 was used in a static mode to shape overlapping electron beam segments, and the dynamic multileaf collimation mode was used to realize the intensity modulated photon beam profiles. Examples of mixed beams were generated at electron energies of up to 40 MeV. The intensity modulated electron beam component consists of two overlapping concentric fields with optimized field sizes, yielding broad, fairly depth-independent overall beam penumbras. The matched intensity modulated photon beam component has high fluence peaks at the field edges to sharpen this penumbra. The combination of the electron and the photon beams yields dose distributions with the characteristics (a)-(c) mentioned above.     

Portal dose image prediction for dosimetric treatment verification in radiotherapy. II. An algorithm for wedged beams

A method is presented for calculation of a two-dimensional function, T(wedge)(x,y), describing the transmission of a wedged photon beam through a patient. This in an extension of the method that we have published for open (nonwedged) fields [Med. Phys. 25, 830-840 (1998)]. Transmission functions for open fields are being used in our clinic for prediction of portal dose images (PDI, i.e., a dose distribution behind the patient in a plane normal to the beam axis), which are compared with PDIs measured with an electronic portal imaging device (EPID). The calculations are based on the planning CT scan of the patient and on the irradiation geometry as determined in the treatment planning process. Input data for the developed algorithm for wedged beams are derived from (the already available) measured input data set for transmission prediction in open beams, which is extended with only a limited set of measurements in the wedged beam. The method has been tested for a PDI plane at 160 cm from the focus, in agreement with the applied focus-to-detector distance of our fluoroscopic EPIDs. For low and high energy photon beams (6 and 23 MV) good agreement (approximately 1%) has been found between calculated and measured transmissions for a slab and a thorax phantom.     

Comparison of measured and Monte Carlo calculated dose distributions from the NRC linac

We have benchmarked photon beam simulations with the EGS4 user code BEAM [Rogers et al., Med. Phys. 22, 503-524 (1995)] by comparing calculated and measured relative ionization distributions in water from the 10 and 20 MV photon beams of the NRC linac. Unlike previous calculations, the incident electron energy is known independently to 1%, the entire extra-focal radiation is simulated, and electron contamination is accounted for. The full Monte Carlo simulation of the linac includes the electron exit window, target, flattening filter, monitor chambers, collimators, as well as the PMMA walls of the water phantom. Dose distributions are calculated using a modified version of the EGS4 user code DOSXYZ which additionally allows scoring of average energy and energy fluence in the phantom. Dose is converted to ionization by accounting for the (L/rho)water(air) variation in the phantom, calculated in an identical geometry for the realistic beams using a new EGS4 user code, SPRXYZ. The variation of (L/rho)water(air) with depth is a 1.25% correction at 10 MV and a 2% correction at 20 MV. At both energies, the calculated and the measured values of ionization on the central axis in the buildup region agree within 1% of maximum ionization relative to the ionization at 10 cm depth. The agreement is well within statistics elsewhere. The electron contamination contributes 0.35(+/- 0.02) to 1.37(+/- 0.03)% of the maximum dose in the buildup region at 10 MV and 0.26(+/- 0.03) to 3.14(+/- 0.07)% of the maximum dose at 20 MV. The penumbrae at 3 depths in each beam (in g/cm2), 1.99 (dmax, 10 MV only), 3.29 (dmax, 20 MV only), 9.79 and 19.79, agree with ionization chamber measurements to better than 1 mm. Possible causes for the discrepancy between calculations and measurements are analyzed and discussed in detail.     

Commissioning of a medical accelerator photon beam Monte Carlo simulation using wide-field profiles

A method for commissioning an EGSnrc Monte Carlo simulation of medical linac photon beams through wide-field lateral profiles at moderate depth in a water phantom is presented. Although depth-dose profiles are commonly used for nominal energy determination, our study shows that they are quite insensitive to energy changes below 0.3 MeV (0.6 MeV) for a 6 MV (15 MV) photon beam. Also, the depth-dose profile dependence on beam radius adds an additional uncertainty in their use for tuning nominal energy. Simulated 40 cm x 40 cm lateral profiles at 5 cm depth in a water phantom show greater sensitivity to both nominal energy and radius. Beam parameters could be determined by comparing only these curves with measured data.     

Monitor unit calculation on the beam axis of open and wedged asymmetric high-energy photon beams

An ESTRO booklet and a report of the Netherlands Commission on Radiation Dosimetry have been published recently describing empirical methods for monitor unit (MU) calculations in symmetrical high-energy photon beams. Both documents support the same basic ideas; firstly the separation of head scatter and volume scatter components and secondly the determination of head scatter quantities in a mini-phantom. Based on these ideas the methods previously described for MU calculations in symmetrical beams are extended to asymmetrical open and wedged beams in isocentric treatment conditions. All required dosimetric parameters (normalized head scatter factors, phantom scatter correction factors, wedge factors, off-axis ratios, quality index, and depth dose parameters) are determined as a function of beam axis position in order to study their off-axis dependence. Measurements are performed for 6 MV and 18 MV photon beams provided by two different dual-energy linear accelerators, a GE Saturne 42 and a Varian 2100 CD linac.     

Monitor unit calculations for wedged asymmetric photon beams

Algorithms for calculating monitor units (MUs) in wedged asymmetric high-energy photon beams as implemented in treatment planning systems have their limitations. Therefore an independent method for MU calculation is necessary. The aim of this study was to develop an empirical method to determine MUs for points at the centre of wedged fields, asymmetric in two directions. The method is based on the determination of an off-axis factor (OAF) that corrects for the difference in dose between wedged asymmetric and wedged symmetric beams with the same field size. Measurements were performed in a water phantom irradiated with 6 and 18 MV photon beams produced by Elekta accelerators, which are fitted with an internal motorized wedge that has a complex shape. The OAF perpendicular to the wedge direction changed significantly with depth for the 18 MV beam. Dose values measured for a set of 18 test cases were compared with those calculated with our method. The maximum difference found was 6.5% and in 15 cases this figure was smaller than 2.0%. The analytical method of Khan and the empirical method of Georg were also tested and showed errors up to 12.8%. It can be concluded that our simple formalism is able to calculate MUs in wedged asymmetric fields with an acceptable accuracy in most clinical situations.     

Response to high-energy photons of PTW31014 PinPoint ion chamber with a central aluminum electrode

Since its introduction the PinPoint (PTW-Freiburg) micro-ionization chamber has been proposed for relative dosimetry (output factors, depth dose curves, and beam profiles) as well as for determination of absolute dose of small high-energy photon beams. This paper investigates the dosimetric performance of a new design (type 31014) of the PinPoint ion chamber with a central aluminum electrode. The study included characterization of inherent and radiation-induced leakage, ion collection efficiency and polarity effect, relative response of the chamber, measurement of beam profiles, and depth dose curves. The 6 and 15 MV photon beams of a Varian 2100 C/D were considered. At the nominal operating voltage of 400 V the PinPoint type 31014 chamber was found to present a strong field size dependence of the polarity correction factor and an excess of the collected charge, which can lead to an underestimation of the collection efficiency if determined with the conventional "two-voltage" method. In comparison to the original PinPoint design (type 31006) the authors found for type 31014 chamber no overresponse to large-area fields if polarity correction is applied. If no correction is taken into consideration, the authors found the chamber's output to be inaccurate for large-area fields (0.5% accuracy limited up to the 12 x 12 and 20 x 20 cm2 field for the 6 and 15 MV beams, respectively), which is a direct consequence of the stem and polarity effects due to the chamber's very small sensitive volume (0.015 cc) and cable irradiation. Beam profiles and depth dose curves measured with type 31014 PinPoint chamber for small and medium size fields were compared to data measured with a 0.125 cc ion chamber and with high-resolution Kodak EDR2 films. Analysis of the penumbra (80%-20% distance) showed that the spatial resolution of type 31014 PinPoint ion chamber approaches (penumbra broadening < or = 0.6 mm) EDR2 film results.     

Properties of unflattened photon beams shaped by a multileaf collimator

Several studies have shown that removal of the flattening filter from the treatment head of a clinical accelerator increases the dose rate and changes the lateral profile in radiation therapy with photons. However, the multileaf collimator (MLC) used to shape the field was not taken into consideration in these studies. We therefore investigated the effect of the MLC on flattened and unflattened beams. To do this, we performed measurements on a Varian Clinac 21EX and MCNPX Monte Carlo simulations to analyze the physical properties of the photon beam. We compared lateral profiles, depth dose curves, MLC leakages, and total scatter factors for two energies (6 and 18 MV) of MLC-shaped fields and jaw-shaped fields. Our study showed that flattening filter-free beams shaped by a MLC differ from the jaw-shaped beams. Similar differences were also observed for flattened beams. Although both collimating methods produced identical depth dose curves, the penumbra size and the MLC leakage were reduced in the softer, unflattened beam and the total scatter factors showed a smaller field size dependence.     

Photon scatter in portal images: accuracy of a fluence based pencil beam superposition algorithm

The accuracy of a pencil beam algorithm to predict scattered photon fluence into portal imaging systems was studied. A data base of pencil beam kernels describing scattered photon fluence behind homogeneous water slabs (1-50 cm thick) at various air gap distances (0-100 cm) was generated using the EGS Monte Carlo code. Scatter kernels were partitioned according to particle history: singly-scattered, multiply-scattered, and bremsstrahlung and positron annihilation photons. Mean energy and mean angle with respect to the incident photon pencil beam were also scored. This data allows fluence, mean energy, and mean angular data for each history type to be predicted using the pencil beam algorithm. Pencil beam algorithm predictions for 6 and 24 MV incident photon beams were compared against full Monte Carlo simulations for several inhomogeneous phantoms, including approximations to a lateral neck, and a mediastinum treatment. The accuracy of predicted scattered photon fluence, mean energy, and mean angle was investigated as a function of air gap, field size, photon history, incident beam resolution, and phantom geometry. Maximum errors in mean energies were 0.65 and 0.25 MeV for the higher and lower energy spectra, respectively, and 15 degrees for mean angles. The ability of the pencil beam algorithm to predict scatter fluence decreases with decreasing air gap, with the largest error for each phantom occurring at the exit surface. The maximum predictive error was found to be 6.9% with respect to the total fluence on the central axis. By maintaining even a small air gap (approximately 10 cm), the error in predicted scatter fluence may be kept under 3% for the phantoms and beam energies studied here. It is concluded that this pencil beam algorithm is sufficiently accurate (using International Commission on Radiation Units and Measurements Report No. 24 guidelines for absorbed dose) over the majority of clinically relevant air gaps, for further investigation in a portal dose prediction algorithm.     

A comparative dosimetric study on tangential photon beams,intensity-modulated radiation therapy (IMRT) and modulated electron radiotherapy (MERT) for breast cancer treatment

Recently, energy- and intensity-modulated electron radiotherapy (MERT) has garnered a growing interest for the treatment of superficial targets. In this work. we carried out a comparative dosimetry study to evaluate MERT, photon beam intensity-modulated radiation therapy (IMRT) and conventional tangential photon beams for the treatment of breast cancer. A Monte Carlo based treatment planning system has been investigated, which consists of a set of software tools to perform accurate dose calculation, treatment optimization, leaf sequencing and plan analysis. We have compared breast treatment plans generated using this home-grown treatment optimization and dose calculation software forthese treatment techniques. The MERT plans were planned with up to two gantry angles and four nominal energies (6, 9, 12 and 16 MeV). The tangential photon treatment plans were planned with 6 MV wedged photon beams. The IMRT plans were planned using both multiple-gantry 6 MV photon beams or two 6 MV tangential beams. Our results show that tangential IMRT can reduce the dose to the lung, heart and contralateral breast compared to conventional tangential wedged beams (up to 50% reduction in high dose volume or 5 Gy in the maximum dose). MERT can reduce the maximum dose to the lung by up to 20 Gy and to the heart by up to 35 Gy compared to conventional tangential wedged beams. Multiple beam angle IMRT can significantly reduce the maximum dose to the lung and heart (up to 20 Gy) but it induces low and medium doses to a large volume of normal tissues including lung, heart and contralateral breast. It is concluded that MERT has superior capabilities to achieve dose conformity both laterally and in the depth direction, which will be well suited for treating superficial targets such as breast cancer.