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.     

Partially wedged radiation beams

To increase dose homogeneity within certain radiotherapy targets, we defined a partially wedged radiation beam as a beam with wedge modification in one part of the field only. Partially wedged beams may be beneficial in cases with curved surfaces inside parts of the beam only, where they may compensate for missing tissue and/or for variations in depth to the target region. Possible sites suitable for partially wedged beams include urinary bladder and tangential breast irradiation. Customized partially wedged beams were delivered applying dynamic collimation techniques. Two different linear detector arrays, a semiconductor diode array and an ionization chamber array, were used independently in the same standard water tank to verify that the partially wedged beams were delivered according to the definition. Dose calculations of partial wedge fields were implemented in our treatment planning system and compared with the measured dose distributions. We re-planned a representative treatment plan for both advanced urinary bladder cancer and tangential breast irradiation using partially wedged beams. For both patients the target dose homogeneity was improved, and the doses to surrounding critical normal tissues were reduced.     

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.     

Characteristics of photon beams from Philips SL25 linear accelerators

The Philips SL25 accelerator is a multimodality machine offering asymmetric collimator jaws and a new type of beam bending and transport system. It produces photon beams, nominally at 6 and 25 MV, and a scattered electron beam with nine selectable energies between 4 and 22 MeV. Dosimetric characteristics for the 6- and 25-MV photon beams are presented with respect to field flatness, surface and depth dose characteristics, isodose distribution, field size factors for both open and wedged fields, and narrow beam transmission data in different materials.     

Testing of the analytical anisotropic algorithm for photon dose calculation

The analytical anisotropic algorithm (AAA) was implemented in the Eclipse (Varian Medical Systems) treatment planning system to replace the single pencil beam (SPB) algorithm for the calculation of dose distributions for photon beams. AAA was developed to improve the dose calculation accuracy, especially in heterogeneous media. The total dose deposition is calculated as the superposition of the dose deposited by two photon sources (primary and secondary) and by an electron contamination source. The photon dose is calculated as a three-dimensional convolution of Monte-Carlo precalculated scatter kernels, scaled according to the electron density matrix. For the configuration of AAA, an optimization algorithm determines the parameters characterizing the multiple source model by optimizing the agreement between the calculated and measured depth dose curves and profiles for the basic beam data. We have combined the acceptance tests obtained in three different departments for 6, 15, and 18 MV photon beams. The accuracy of AAA was tested for different field sizes (symmetric and asymmetric) for open fields, wedged fields, and static and dynamic multileaf collimation fields. Depth dose behavior at different source-to-phantom distances was investigated. Measurements were performed on homogeneous, water equivalent phantoms, on simple phantoms containing cork inhomogeneities, and on the thorax of an anthropomorphic phantom. Comparisons were made among measurements, AAA, and SPB calculations. The optimization procedure for the configuration of the algorithm was successful in reproducing the basic beam data with an overall accuracy of 3%, 1 mm in the build-up region, and 1%, 1 mm elsewhere. Testing of the algorithm in more clinical setups showed comparable results for depth dose curves, profiles, and monitor units of symmetric open and wedged beams below dmax. The electron contamination model was found to be suboptimal to model the dose around dmax, especially for physical wedges at smaller source to phantom distances. For the asymmetric field verification, absolute dose difference of up to 4% were observed for the most extreme asymmetries. Compared to the SPB, the penumbra modeling is considerably improved (1%, 1 mm). At the interface between solid water and cork, profiles show a better agreement with AAA. Depth dose curves in the cork are substantially better with AAA than with SPB. Improvements are more pronounced for 18 MV than for 6 MV. Point dose measurements in the thoracic phantom are mostly within 5%. In general, we can conclude that, compared to SPB, AAA improves the accuracy of dose calculations. Particular progress was made with respect to the penumbra and low dose regions. In heterogeneous materials, improvements are substantial and more pronounced for high (18 MV) than for low (6 MV) energies.     

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.     

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.     

Backscatter towards the monitor ion chamber in high-energy photon and electron beams: charge integration versus Monte Carlo simulation

In some linear accelerators, the charge collected by the monitor ion chamber is partly caused by backscattered particles from accelerator components downstream from the chamber. This influences the output of the accelerator and also has to be taken into account when output factors are derived from Monte Carlo simulations. In this work, the contribution of backscattered particles to the monitor ion chamber response of a Varian 2100C linac was determined for photon beams (6, 10 MV) and for electron beams (6, 12, 20 MeV). The experimental procedure consisted of charge integration from the target in a photon beam or from the monitor ion chamber in electron beams. The Monte Carlo code EGS4/BEAM was used to study the contribution of backscattered particles to the dose deposited in the monitor ion chamber. Both measurements and simulations showed a linear increase in backscatter fraction with decreasing field size for photon and electron beams. For 6 MV and 10 MV photon beams, a 2-3% increase in backscatter was obtained for a 0.5 x 0.5 cm2 field compared to a 40 x 40 cm2 field. The results for the 6 MV beam were slightly higher than for the 10 MV beam. For electron beams (6, 12, 20 MeV), an increase of similar magnitude was obtained from measurements and simulations for 6 MeV electrons. For higher energy electron beams a smaller increase in backscatter fraction was found. The problem is of less importance for electron beams since large variations of field size for a single electron energy usually do not occur.     

Optimized treatment planning for prostate cancer comparing IMPT,VHEET and 15 MV IMXT

The merits of intensity-modulated very-high energy electron therapy (VHEET) and intensity-modulated proton therapy (IMPT) in relation to intensity-modulated x-ray therapy (IMXT) with respect to the treatment of the prostate have been quantified. Optimized dose distributions were designed for 5-11 beams of 250 MeV VHEET and 15 MV IMXT as well as 1-9 beam ports of IMPT. In the case of the comparison between 250 MeV VHEET and 15 MV IMXT, it was found that the quality of target coverage achievable with VHEET was comparable to or sometimes better than that provided by IMXT. However, VHEET provided an improvement over IMXT in the dose sparing of the sensitive structures and normal tissues. Compared to IMXT, VHEET decreased the mean rectal dose and bladder dose by up to 10% of the prescribed target dose, while reducing by up to 12% of the prescribed target dose the integral dose to normal tissues. In quantifying the merits of IMPT relative to IMXT, it was found that using intensity-modulated proton beams for inverse planning instead of intensity-modulated photon beams improved target dose homogeneity by up to 1.3% of the prescribed target dose, while reducing the mean rectal dose, bladder dose, and normal tissue integral dose by up to 27%, 30% and 28% of the prescribed target dose respectively. The comparison of optimized planning for IMPT and VHEET showed that the quality of target coverage achievable with IMPT is comparable to or better (by up to 1.3% of the prescribed target dose) than that provided by VHEET. Compared to VHEET, IMPT delivered a mean rectal dose and a bladder dose that was lower by up to 17% and 23% of prescribed target dose respectively, and also reduced the integral dose to normal tissues by up to 17% of the prescribed target dose. These results indicate that of the three modalities the greatest dose escalation will be possible with IMPT, then VHEET, and then IMXT. It follows that IMPT will result in the highest probability of complication-free tumour control, while IMXT will provide the lowest probability.     

Optimal electron and combined electron and photon therapy in the phase space of complication-free cure

The possibility of using intensity-modulated high-energy electrons beams alone or in combination with photon beams to treat tumours located at depths from 5 cm to 25 cm has been investigated. A radiobiologically based optimization algorithm using the probability of complication-free tumour control has been used to calculate the optimal dose distributions. Two different target volumes have been used; one advanced cervical cancer with locally involved lymph nodes and one astrocytoma in the upper brain hemisphere. Treatments with only electron beams and also combinations between electron and photon beams have been investigated. The dependence of the expected treatment outcome on the beam energy and directions was investigated, and to some extent on the number of beam portals. It is shown that the beam direction intervals resulting in a high expected treatment outcome increase with increasing electron energy and also with some electron-photon combinations. For an eccentrically placed, not too deeply situated tumour surrounded by sensitive normal tissue it is shown that the expected treatment outcome can be improved by using electron beams in combination with photon beams compared with using two photon beams, and using two electron beams results in almost as high an expected treatment outcome. The possibility of improving the dose conformity from electron beams by adding photon fields parallel or orthogonal to the electron beams is demonstrated.     

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.     

The characterization of unflattened photon beams from a 6 MV linear accelerator

Commissioning data have been measured for an Elekta Precise linear accelerator running at 6 MV without a flattening filter with the aim of studying the effects of flattening filter removal on machine operation and beam characterization. Modern radiotherapy practice now routinely relies on the use of fluence modifying techniques such as IMRT, i.e. the active production of non-flat beams. For these techniques the flattening filter should not be necessary. It is also possible that the increased intensity around the central axis associated with unflattened beams may be useful for conventional treatment planning by acting as a field-in-field or integrated boost technique. For this reason open and wedged field data are presented. Whilst problems exist in running the machine filter free clinically, this paper shows that in many ways the beam is actually more stable, exhibiting almost half the variation in field symmetry for changes in steering and bending currents. Dosimetric benefits are reported here which include a reduction in head scatter by approx. 70%, decreased penumbra (0.5 mm), lower dose outside of the field edge (11%) and a doubling in dose rate (2.3 times for open and 1.9 times for wedged fields). Measurements also show that reduced scatter also reduces leakage radiation by approx. 60%, significantly lowering whole body doses. The greatest benefit of filter-free use is perceived to be for IMRT where increased dose rate combined with reduced head scatter and leakage radiation should lead to improved dose calculation, giving simpler, faster and more accurate dose delivery with reduced dose to normal tissues.     

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.     

Photonuclear dose calculations for high-energy photon beams from Siemens and Varian linacs

The dose from photon-induced nuclear particles (neutrons, protons, and alpha particles) generated by high-energy photon beams from medical linacs is investigated. Monte Carlo calculations using the MCNPX code are performed for three different photon beams from two different machines: Siemens 18 MV, Varian 15 MV, and Varian 18 MV. The linac head components are simulated in detail. The dose distributions from photons, neutrons, protons, and alpha particles are calculated in a tissue-equivalent phantom. Neutrons are generated in both the linac head and the phantom. This study includes (a) field size effects, (b) off-axis dose profiles, (c) neutron contribution from the linac head, (d) dose contribution from capture gamma rays, (e) phantom heterogeneity effects, and (f) effects of primary electron energy shift. Results are presented in terms of absolute dose distributions and also in terms of DER (dose equivalent ratio). The DER is the maximum dose from the particle (neutron, proton, or alpha) divided by the maximum photon dose, multiplied by the particle quality factor and the modulation scaling factor. The total DER including neutrons, protons, and alphas is about 0.66 cSv/Gy for the Siemens 18 MV beam (10 cm x 10 cm). The neutron DER decreases with decreasing field size while the proton (or alpha) DER does not vary significantly except for the 1 cm x 1 cm field. Both Varian beams (15 and 18 MV) produce more neutrons, protons, and alphas particles than the Siemens 18 MV beam. This is mainly due to their higher primary electron energies: 15 and 18.3 MeV, respectively, vs 14 MeV for the Siemens 18 MV beam. For all beams, neutrons contribute more than 75% of the total DER, except for the 1 cm x 1 cm field (approximately 50%). The total DER is 1.52 and 2.86 cSv/Gy for the 15 and 18 MV Varian beams (10 cm x 10 cm), respectively. Media with relatively high-Z elements like bone may increase the dose from heavy charged particles by a factor 4. The total DER is sensitive to primary electron energy shift. A Siemens 18 MV beam with 15 MeV (instead of 14 MeV) primary electrons would increase by 40% the neutron DER and by 210% the proton + alpha DER. Comparisons with measurements (neutron yields from different materials and neutron dose equivalent) are also presented. Using the NCRP risk assessment method, we found that the dose equivalent from leakage neutrons (at 50-cm off-axis distance) represent 1.1, 1.1, and 2.0% likelihood of fatal secondary cancer for a 70 Gy treatment delivered by the Siemens 18 MV, Varian 15 MV, and Varian 18 MV beams, respectively.     

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.     

A 4 MV flattening filter-free beam: commissioning and application to conformal therapy and volumetric modulated arc therapy

Recent studies have indicated that radiotherapy treatments undertaken on a flattening filter-free (FFF) linear accelerator have a number of advantages over treatments undertaken on a conventional linear accelerator. In addition, 4 MV photon beams may give improved isodose coverage for some treatment volumes at air/tissue interfaces, compared to when utilizing the clinical standard of 6 MV photons. In order to investigate these benefits, FFF beams were established on an Elekta Beam Modulator linear accelerator for 4 MV photons. Commissioning beam data were obtained for open and wedged fields. The measured data were then imported into a treatment planning system and a beam model was commissioned. The beam model was optimized to improve dose calculations at shallow, clinically relevant depths. Following verification, the beam model was utilized in a treatment planning study, including volumetric modulated arc therapy, for a selection of lung, breast/chest wall and larynx patients. Increased dose rates of around 800 MU min(-1) were recorded for open fields (relative to 320 MU min(-1) for filtered open fields) and reduced head scatter was inferred from output factor measurements. Good agreement between planned and delivered dose was observed in verification of treatment plans. The planning study indicated that with a FFF beam, equivalent (and in some cases improved) isodose profiles could be achieved for small lung and larynx treatment volumes relative to 4 MV filtered treatments. Furthermore, FFF treatments with wedges could be replicated using open fields together with an 'effective wedge' technique and isocentre shift. Clinical feasibility of a FFF beam was therefore demonstrated, with beam modelling, treatment planning and verification being successfully accomplished.     

Characteristics of the high-energy photon beam of a 25-MeV accelerator

The CGR Saturne 25 is an isocentrically mounted standing wave medical linear accelerator that produces dual-energy photon beams and a scanned electron beam with six selectable energies between 4 and 25 MeV. The highest energy photon beam is nominally referred to as 23 MV. For this beam the mean energy of the accelerated electron beam on the 1.3 radiation length (4 mm) tungsten x-ray target is found to be approximately 21 MeV, with the energy acceptance stated to be +/- 5%. The electron beam traverses a 270 degrees bending magnet upstream of the x-ray production target. The resulting bremsstrahlung beam passes through a combination steel and lead flattening filter, 4-cm maximum thickness. Dosimetric data for the 23-MV beam are presented with respect to rectangular field output factor, depth of maximum dose as a function of field size, surface and buildup dose, central axis percent depth dose, tissue-phantom ratios, beam profile, applicability of inverse square, and block transmission. Some data are also presented on the effect of different flattening filter designs on apparent beam energy.