A monte carlo approach to electron contamination sources in the Saturne-25 and -41

The various components of the accelerator treatment head act as sources of contaminating electrons. The presence of contamination electrons increases the surface dose, which deteriorates the skin-sparing effect. The present study examines the sources of this 'contamination', the influence on the surface dose and the shape of the build-up curve. The Monte Carlo simulation of two linear accelerators, Saturne-25 and -41, allowed us to study the influence of electron contamination in various therapeutic energies and in different geometries. The Saturne-25 and -41 cover a wide range of therapeutic energies with nominal energies 12/23 MV and 6/15 MV, respectively. The analysis of the results shows that at a source-to-surface distance of 100 cm and a wide opening of the collimators, the main sources of contaminating electrons are the flattening filter and the air below it. The contribution of the secondary contamination electrons on the surface dose is 16% for 6 MV and 12 MV, 6% for 15 MV and 17% for 23 MV. The energy spectra of electrons coming from the flattening filter and the air below it are completely different. The air produces electrons of low energies. The mean energies of these spectra vary from 1 MeV to 2 MeV depending on the nominal energy of the photon beam. The secondary electrons generated by the flattening filter produce a wide energy spectrum with mean energies of the same order of the bremsstrahlung spectrum. The flattening filter absorbs the secondary electrons generated in the target, the primary collimator and the air inside the head.     

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%.     

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.     

Monte Carlo study of backscatter in a flattening filter free clinical accelerator

In conventional linear accelerators, the flattening filter provides a uniform lateral dose profile. In intensity modulated radiation therapy applications, however, the flatness of the photon field and hence the presence of a flattening filter, is not necessary. Removing the filter may provide some advantages, such as faster treatments and smaller out-of-field doses to the patients. In clinical accelerators the backscattered radiation dose from the collimators must be taken into account when the dose to the target volume in the patient is being determined. In the case of a conventional machine, this backscatter is known to great precision. In a flattening filter free accelerator, however, the amount of backscatter may be different. In this study we determined the backscatter contribution to the monitor chamber signal in a flattening filter free clinical accelerator (Varian Clinac 21EX) with Monte Carlo simulations. We found that with the exception of very small fields in the 18-MV photon mode, the contribution of backscattered radiation to the monitor signal did not differ from that of conventional machines with a flattening filter. Hence, a flattening filter free clinical accelerator would not necessitate a different backscatter correction.     

Monte Carlo study of photon fields from a flattening filter-free clinical accelerator

In conventional clinical linear accelerators, the flattening filter scatters and absorbs a large fraction of primary photons. Increasing the beam-on time, which also increases the out-of-field exposure to patients, compensates for the reduction in photon fluence. In recent years, intensity modulated radiation therapy has been introduced, yielding better dose distributions than conventional three-dimensional conformal therapy. The drawback of this method is the further increase in beam-on time. An accelerator with the flattening filter removed, which would increase photon fluence greatly, could deliver considerably higher dose rates. The objective of the present study is to investigate the dosimetric properties of 6 and 18 MV photon beams from an accelerator without a flattening filter. The dosimetric data were generated using the Monte Carlo programs BEAMnrc and DOSXYZnrc. The accelerator model was based on the Varian Clinac 2100 design. We compared depth doses, dose rates, lateral profiles, doses outside collimation, total and collimator scatter factors for an accelerator with and without a flatteneing filter. The study showed that removing the filter increased the dose rate on the central axis by a factor of 2.31 (6 MV) and 5.45 (18 MV) at a given target current. Because the flattening filter is a major source of head scatter photons, its removal from the beam line could reduce the out-of-field dose.     

Monte Carlo study of a Cyberknife stereotactic radiosurgery system

This study investigated small-field dosimetry for a Cyberknife stereotactic radiosurgery system using Monte Carlo simulations. The EGSnrc/BEAMnrc Monte Carlo code was used to simulate the Cyberknife treatment head, and the DOSXYZnrc code was implemented to calculate central axis depth-dose curves, off-axis dose profiles, and relative output factors for various circular collimator sizes of 5 to 60 mm. Water-to-air stopping power ratios necessary for clinical reference dosimetry of the Cyberknife system were also evaluated by Monte Carlo simulations. Additionally, a beam quality conversion factor, kQ, for the Cyberknife system was evaluated for cylindrical ion chambers with different wall material. The accuracy of the simulated beam was validated by agreement within 2% between the Monte Carlo calculated and measured central axis depth-dose curves and off-axis dose profiles. The calculated output factors were compared with those measured by a diode detector and an ion chamber in water. The diode output factors agreed within 1% with the calculated values down to a 10 mm collimator. The output factors with the ion chamber decreased rapidly for collimators below 20 mm. These results were confirmed by the comparison to those from Monte Carlo methods with voxel sizes and materials corresponding to both detectors. It was demonstrated that the discrepancy in the 5 and 7.5 mm collimators for the diode detector is due to the water non-equivalence of the silicon material, and the dose fall-off for the ion chamber is due to its large active volume against collimators below 20 mm. The calculated stopping power ratios of the 60 mm collimator from the Cyberknife system (without a flattening filter) agreed within 0.2% with those of a 10 X 10 cm2 field from a conventional linear accelerator with a heavy flattening filter and the incident electron energy, 6 MeV. The difference in the stopping power ratios between 5 and 60 mm collimators was within 0.5% at a 10 cm depth in water. Furthermore, kQ values for the Cyberknife system were in agreement within 0.3% with those of the conventional 6 MV-linear accelerator for the cylindrical ion chambers with different wall material.     

Radiosurgery with unflattened 6-MV photon beams

One of the major drawbacks to doing stereotactic radiosurgery with a linear accelerator is the long time required to deliver the target dose. Single fractions of 25 Gy delivered at the isocenter and at depth in the skull may require beam times in excess of 15 min for a typical linear accelerator with a maximum dose rate of 250 cGy/min in tissue. In an effort to decrease the treatment time for this technique, the flattening filter has been removed from an AECL Therac-6 linear accelerator and the characteristics of the resulting beam have been measured. Flatness is acceptable for the field sizes used with this technique and the dose rate is increased by a factor of 2.75.     

A semianalytical method for the design of a linac x-ray beam flattening filter

The purpose of this study was to design an improved flattening filter for a Therac 20 medical linear accelerator. Profiles of the 18-MV x-ray beam produced by this accelerator measured along the diagonal of a 40 X 40 cm field at a depth of 5 cm were measured, and it was found that there were regions near the corners of the field where the dose was 109% of the central axis dose. An iterative algorithm for designing flattening filters was developed which required, as input, precise measurements of the following data: the unflattened primary beam profile, the fraction of the beam due to contamination radiation arising from interactions of primary photons with the flattening filter and the collimator assemblies, and the attenuation of the primary photons in water and lead as a function of angle from the central axis of the beam. A new flattening filter was designed and profiles of the beam were measured at a number of depths. These measurements showed that the beam was flattened to within +/- 1% out to 24 cm along the diagonal of a 40 X 40 cm field at a depth of 5 cm.     

A Monte Carlo multiple source model applied to radiosurgery narrow photon beams

Monte Carlo (MC) methods are nowadays often used in the field of radiotherapy. Through successive steps, radiation fields are simulated, producing source Phase Space Data (PSD) that enable a dose calculation with good accuracy. Narrow photon beams used in radiosurgery can also be simulated by MC codes. However, the poor efficiency in simulating these narrow photon beams produces PSD whose quality prevents calculating dose with the required accuracy. To overcome this difficulty, a multiple source model was developed that enhances the quality of the reconstructed PSD, reducing also the time and storage capacities. This multiple source model was based on the full MC simulation, performed with the MC code MCNP4C, of the Siemens Mevatron KD2 (6 MV mode) linear accelerator head and additional collimators. The full simulation allowed the characterization of the particles coming from the accelerator head and from the additional collimators that shape the narrow photon beams used in radiosurgery treatments. Eight relevant photon virtual sources were identified from the full characterization analysis. Spatial and energy distributions were stored in histograms for the virtual sources representing the accelerator head components and the additional collimators. The photon directions were calculated for virtual sources representing the accelerator head components whereas, for the virtual sources representing the additional collimators, they were recorded into histograms. All these histograms were included in the MC code, DPM code and using a sampling procedure that reconstructed the PSDs, dose distributions were calculated in a water phantom divided in 20000 voxels of 1 x 1 x 5 mm3. The model accurately calculates dose distributions in the water phantom for all the additional collimators; for depth dose curves, associated errors at 2sigma were lower than 2.5% until a depth of 202.5 mm for all the additional collimators and for profiles at various depths, deviations between measured and calculated values were less than 2.5% or 1 mm.     

A flattening filter free photon treatment concept evaluation with Monte Carlo

In principle, the concept of flat initial radiation-dose distribution across the beam is unnecessary for intensity modulated radiation therapy. Dynamic leaf positioning during irradiation could appropriately adjust the fluence distribution of an unflattened beam that is peaked in the center and deliver the desired uniform or nonuniform dose distribution. Removing the flattening filter could lead to reduced treatment time through higher dose rates and reduced scatter, because there would be substantially less material in the beam; and possibly other dosimetric and clinical advantages. This work aims to evaluate the properties of a flattening filter free clinical accelerator and to investigate its possible advantages in clinical intensity modulated radiation therapy applications by simulating a Varian 2100-based treatment delivery system with Monte Carlo techniques. Several depth-dose curves and lateral dose distribution profiles have been created for various field sizes, with and without the flattening filter. Data computed with this model were used to evaluate the overall quality of such a system in terms of changes in dose rate, photon and electron fluence, and reduction in out-of-field stray dose from the scattered components and were compared to the corresponding data for a standard treatment head with a flattening filter. The results of the simulations of the flattening filter free system show that a substantial increase in dose rate can be achieved, which would reduce the beam on time and decrease the out-of-field dose for patients due to reduced head-leakage dose. Also close to the treatment field edge, a significant improvement in out-of-field dose could be observed for small fields, which can be attributed to the change in the photon spectra, when the flattening filter is removed from the beamline.     

A virtual photon energy fluence model for Monte Carlo dose calculation

The presented virtual energy fluence (VEF) model of the patient-independent part of the medical linear accelerator heads, consists of two Gaussian-shaped photon sources and one uniform electron source. The planar photon sources are located close to the bremsstrahlung target (primary source) and to the flattening filter (secondary source), respectively. The electron contamination source is located in the plane defining the lower end of the filter. The standard deviations or widths and the relative weights of each source are free parameters. Five other parameters correct for fluence variations, i.e., the horn or central depression effect. If these parameters and the field widths in the X and Y directions are given, the corresponding energy fluence distribution can be calculated analytically and compared to measured dose distributions in air. This provides a method of fitting the free parameters using the measurements for various square and rectangular fields and a fixed number of monitor units. The next step in generating the whole set of base data is to calculate monoenergetic central axis depth dose distributions in water which are used to derive the energy spectrum by deconvolving the measured depth dose curves. This spectrum is also corrected to take the off-axis softening into account. The VEF model is implemented together with geometry modules for the patient specific part of the treatment head (jaws, multileaf collimator) into the XVMC dose calculation engine. The implementation into other Monte Carlo codes is possible based on the information in this paper. Experiments are performed to verify the model by comparing measured and calculated dose distributions and output factors in water. It is demonstrated that open photon beams of linear accelerators from two different vendors are accurately simulated using the VEF model. The commissioning procedure of the VEF model is clinically feasible because it is based on standard measurements in air and water. It is also useful for IMRT applications because a full Monte Carlo simulation of the treatment head would be too time-consuming for many small fields.     

Energy and angular distributions of photons from medical linear accelerators

For accurate three-dimensional treatment planning, new models of dose calculations are being developed which require the knowledge of the energy spectra and angular distributions of the photons incident on the surface of the patient. Knowledge of the spectra is also useful in other applications, including the design of filters and beam modifying devices and determination of factors to convert ionization chamber measurements to dose. We have used Monte Carlo code (EGS) to compute photon spectra for a number of different linear accelerators. Both the target and the flattening filter have been accurately modeled. We find the mean photon energy to have a value lower than the generally perceived value of one-third the maximum energy. As expected, the spectra become softer as the distance from the central axis increases. Verification of the spectra is performed by computing dose distributions and half-value layers in water using the calculated spectra and comparing the results with measured data. We also examined the angular distributions of photons incident on the surface of the phantom. In currently used models of dose computations, it is assumed that the angular distribution of photons with respect to fan lines emanating from the source is negligible. Although the angular spread of photons with respect to the incident direction has been found to be small, its contribution to the diffuseness of the beam boundaries is significant.     

Reference dosimetry condition and beam quality correction factor for CyberKnife beam

This article is intended to improve the certainty of the absorbed dose determination for reference dosimetry in CyberKnife beams. The CyberKnife beams do not satisfy some conditions of the standard reference dosimetry protocols because of its unique treatment head structure and beam collimating system. Under the present state of affairs, the reference dosimetry has not been performed under uniform conditions and the beam quality correction factor kQ for an ordinary 6 MV linear accelerator has been temporally substituted for the kQ of the CyberKnife in many sites. Therefore, the reference conditions and kQ as a function of the beam quality index in a new way are required. The dose flatness and the error of dosimeter reading caused by radiation fields and detector size were analyzed to determine the reference conditions. Owing to the absence of beam flattening filter, the dose flatness of the CyberKnife beam was inferior to that of an ordinary 6 MV linear accelerator. And if the absorbed dose is measured with an ionization chamber which has cavity length of 2.4, 1.0 and 0.7 cm in reference dosimetry, the dose at the beam axis for a field of 6.0 cm collimator was underestimated 1.5%, 0.4%, and 0.2% on a calculation. Therefore, the maximum field shaped with a 6.0 cm collimator and ionization chamber which has a cavity length of 1.0 cm or shorter were recommended as the conditions of reference dosimetry. Furthermore, to determine the kQ for the CyberKnife, the realistic energy spectrum of photons and electrons in water was simulated with the BEAMnrc. The absence of beam flattening filter also caused softer photon energy spectrum than that of an ordinary 6 MV linear accelerator. Consequently, the kQ for ionization chambers of a suitable size were determined and tabulated as a function of measurable beam quality indexes in the CyberKnife beam.     

A Monte Carlo study of a flattening filter-free linear accelerator verified with measurements

A Monte Carlo model of an Elekta Precise linear accelerator has been built and verified by measured data for a 6 and 10 MV photon beam running with and without a flattening filter in the beam line. In this study the flattening filter was replaced with a 6 mm thick copper plate, provided by the linac vendor, in order to stabilize the beam. Several studies have shown that removal of the filter improves some properties of the photon beam, which could be beneficial for radiotherapy treatments. The investigated characteristics of this new beam included output, spectra, mean energy, half value layer and the origin of scattered photons. The results showed an increased dose output per initial electron at the central axis of 1.76 and 2.66 for the 6 and 10 MV beams, respectively. The number of scattered photons from the accelerator head was reduced by (31.7 ± 0.03)% (1 SD) for the 6 MV beam and (47.6 ± 0.02)% for the 10 MV beam. The photon energy spectrum of the unflattened beam was softer compared to a conventional beam and did not vary significantly with the off-axis distance, even for the largest field size (0-20 cm off-axis).     

Dosimetric properties of photon beams from a flattening filter free clinical accelerator

Basic dosimetric properties of 6 MV and 18 MV photon beams from a Varian Clinac 21EX accelerator operating without the flattening filter have been measured. These include dose rate data, depth dose dependencies and lateral profiles in a water phantom, total scatter factors and transmission factors of a multileaf collimator. The data are reviewed and compared with measurements for the flattened beams. The unflattened beams have the following: a higher dose rate by factors of 2.3 (6 MV) and 5.5 (18 MV) on the central axis; lower out-of-field dose due to reduced head scatter and softer spectra; less variation of the total scatter factor with field size; and less variation of the shape of lateral dose profiles with depth. The findings suggest that with a flattening filter free accelerator better radiation treatments can be developed, with shorter delivery times and lower doses to normal tissues and organs.     

Radiation characteristics of helical tomotherapy

Helical tomotherapy is a dedicated intensity modulated radiation therapy (IMRT) system with on-board imaging capability (MVCT) and therefore differs from conventional treatment units. Different design goals resulted in some distinctive radiation field characteristics. The most significant differences in the design are the lack of flattening filter, increased shielding of the collimators, treatment and imaging operation modes and narrow fan beam delivery. Radiation characteristics of the helical tomotherapy system, sensitivity studies of various incident electron beam parameters and radiation safety analyses are presented here. It was determined that the photon beam energy spectrum of helical tomotherapy is similar to that of more conventional radiation treatment units. The two operational modes of the system result in different nominal energies of the incident electron beam with approximately 6 MeV and 3.5 MeV in the treatment and imaging modes, respectively. The off-axis mean energy dependence is much lower than in conventional radiotherapy units with less than 5% variation across the field, which is the consequence of the absent flattening filter. For the same reason the transverse profile exhibits the characteristic conical shape resulting in a 2-fold increase of the beam intensity in the center. The radiation leakage outside the field was found to be negligible at less than 0.05% because of the increased shielding of the collimators. At this level the in-field scattering is a dominant source of the radiation outside the field and thus a narrow field treatment does not result in the increased leakage. The sensitivity studies showed increased sensitivity on the incident electron position because of the narrow fan beam delivery and high sensitivity on the incident electron energy, as common to other treatment systems. All in all, it was determined that helical tomotherapy is a system with some unique radiation characteristics, which have been to a large extent optimized for intensity modulated delivery.     

Fast Fourier transform convolution calculations of x-ray isodose distributions in homogeneous media

Convolution concepts were implemented using the discrete fast Fourier transform (FFT) to model the three-dimensional dose distribution due to x-rays produced by a medical linear accelerator. Convolution kernels were employed that had been calculated by Mackie using the EGS4 Monte Carlo code. The EGS4 code was also used to estimate initially the spectrum by simulating the production, filtering, and flattening of the beam in the collimator of the linear accelerator. The continuous bremsstrahlung spectrum was modeled using five discrete energies. The more subtle field-size effects of collimator scattering on the spectrum were obtained by calculating corrections to the spectral components using a least-squares search technique. Dose distributions were obtained using FFT convolutions of the kernels for each energy with the spectrally weighted fluence distributions for that energy. The dose distributions were compared with isodose distributions measured in a water phantom. The agreement was generally found to be better than 1% on the central axis. The calculation time for a single three-dimensional beam was approximately 20 min using a VAX/750 without an array processor. Methods were explored to reduce the calculation time using similar hardware, and estimates were made of how to reduce the calculation time using a more sophisticated computer system.     

Sensitivity of megavoltage photon beam Monte Carlo simulations to electron beam and other parameters

The BEAM code is used to simulate nine photon beams from three major manufacturers of medical linear accelerators (Varian, Elekta, and Siemens), to derive and evaluate estimates for the parameters of the electron beam incident on the target, and to study the effects of some mechanical parameters like target width, primary collimator opening, flattening filter material and density. The mean energy and the FWHM of the incident electron beam intensity distributions (assumed Gaussian and cylindrically symmetric) are derived by matching calculated percentage depth-dose curves past the depth of maximum dose (within 1% of maximum dose) and off-axis factors (within 2sigma at 1% statistics or less) with measured data from the AAPM RTC TG-46 compilation. The off-axis factors are found to be very sensitive to the mean energy of the electron beam, the FWHM of its intensity distribution, its angle of incidence, the dimensions of the upper opening of the primary collimator, the material of the flattening filter and its density. The off-axis factors are relatively insensitive to the FWHM of the electron beam energy distribution, its divergence and the lateral dimensions of the target. The depth-dose curves are sensitive to the electron beam energy, and to its energy distribution, but they show no sensitivity to the FWHM of the electron beam intensity distribution. The electron beam incident energy can be estimated within 0.2 MeV when matching either the measured off-axis factors or the central-axis depth-dose curves when the calculated uncertainties are about 0.7% at the 1 sigma level. The derived FWHM (+/-0.1 mm) of the electron beam intensity distributions all fall within 1 mm of the manufacturer specifications except in one case where the difference is 1.2 mm.     

Generation and modelling of megavoltage photon beams for contrast-enhanced radiation therapy

Contrast-enhanced radiation therapy (CERT) is a treatment approach involving the irradiation of tumours containing high atomic number (Z) contrast media, using low-quality x-ray beams. This work describes the experimental generation of x-ray beams using a linear accelerator with low-Z target materials (beryllium and aluminium), in order to produce photon energy spectra appropriate for CERT. Measurements were made to compare the experimental beams to conventional linear accelerator photon beams in terms of per cent depth dose. Monte Carlo simulation was used to model the generation of each beam, and models were validated against experimental measurement. Validated models were used to demonstrate changes in photon spectra as well as to quantify the variation of tumour dose enhancement with iodinated contrast medium concentration in a simulated tumour volume. Finally, the ratio of the linear attenuation coefficient for iodinated contrast medium relative to water was determined experimentally as a function of iodine concentration. Beams created with low-Z targets show significant changes in energy spectra compared to conventional beams. For the 4 MeV/Be beam, for example, 33% of photons have energies below 60 keV. Measurements and calculation show that both the linear attenuation coefficient ratio and dose enhancement factor (DEF) increase most rapidly at concentrations below 46 mg I ml(-1). There is a significant dependence of DEF on electron energy and a lesser dependence on target material. The 4 MeV/Be beam is the most promising in terms of magnitude of DEF - for example, DEF values of 1.16 and 1.29 are obtained for concentrations of 20 mg I ml(-1) and 50 mg I ml(-1), respectively. DEF will increase or decrease, respectively, for shallower or deeper tumours at a rate of approximately 1.1% cm(-1). In summary, we show that significant dose enhancement is possible by altering the linear accelerator target and filtration, but the magnitude is highly dependent on contrast medium concentration.