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

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.     

A computational implementation and comparison of several intensity modulated proton therapy treatment planning algorithms

The authors present a comparative study of intensity modulated proton therapy (IMPT) treatment planning employing algorithms of three-dimensional (3D) modulation, and 2.5-dimensional (2.5D) modulation, and intensity modulated distal edge tracking (DET) [A. Lomax, Phys. Med. Biol. 44, 185-205 (1999)] applied to the treatment of head-and-neck cancer radiotherapy. These three approaches were also compared with 6 MV photon intensity modulated radiation therapy (IMRT). All algorithms were implemented in the University of Florida Optimized Radiation Therapy system using a finite sized pencil beam dose model and a convex fluence map optimization model. The 3D IMPT and the DET algorithms showed considerable advantages over the photon IMRT in terms of dose conformity and sparing of organs at risk when the beam number was not constrained. The 2.5D algorithm did not show an advantage over the photon IMRT except in the dose reduction to the distant healthy tissues, which is inherent in proton beam delivery. The influences of proton beam number and pencil beam size on the IMPT plan quality were also studied. Out of 24 cases studied, three cases could be adequately planned with one beam and 12 cases could be adequately planned with two beams, but the dose uniformity was often marginally acceptable. Adding one or two more beams in each case dramatically improved the dose uniformity. The finite pencil beam size had more influence on the plan quality of the 2.5D and DET algorithms than that of the 3D IMPT. To obtain a satisfactory plan quality, a 0.5 cm pencil beam size was required for the 3D IMPT and a 0.3 cm size was required for the 2.5D and the DET algorithms. Delivery of the IMPT plans produced in this study would require a proton beam spot scanning technique that has yet to be developed clinically.     

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.     

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.     

Photon-beam subsource sensitivity to the initial electron-beam parameters

One limitation to the widespread implementation of Monte Carlo (MC) patient dose-calculation algorithms for radiotherapy is the lack of a general and accurate source model of the accelerator radiation source. Our aim in this work is to investigate the sensitivity of the photon-beam subsource distributions in a MC source model (with target, primary collimator, and flattening filter photon subsources and an electron subsource) for 6- and 18-MV photon beams when the energy and radial distributions of initial electrons striking a linac target change. For this purpose, phase-space data (PSD) was calculated for various mean electron energies striking the target, various normally distributed electron energy spread, and various normally distributed electron radial intensity distributions. All PSD was analyzed in terms of energy, fluence, and energy fluence distributions, which were compared between the different parameter sets. The energy spread was found to have a negligible influence on the subsource distributions. The mean energy and radial intensity significantly changed the target subsource distribution shapes and intensities. For the primary collimator and flattening filter subsources, the distribution shapes of the fluence and energy fluence changed little for different mean electron energies striking the target, however, their relative intensity compared with the target subsource change, which can be accounted for by a scaling factor. This study indicates that adjustments to MC source models can likely be limited to adjusting the target subsource in conjunction with scaling the relative intensity and energy spectrum of the primary collimator, flattening filter, and electron subsources when the energy and radial distributions of the initial electron-beam change.     

Verification of IMRT fields by film dosimetry

In intensity modulated radiation therapy (IMRT) the aim of an accurate conformal dose distribution is obtained through a complex process. This ranges from the calculation of the optimal distribution of fluence by the treatment planning system (TPS), to the dose delivery through a multilamellar collimator (MLC), with several segments per beam in the step and shoot approach. The above-mentioned consideration makes mandatory an accurate dosimetric verification of the IM beams. A high resolution and integrating dosimeter, like the radiographic film, permits one to simultaneously measure the dose in a matrix of points, providing a good means of obtaining dose distributions. The intrinsic limitation of film dosimetry is the sensitivity dependence on the field size and on the measurement depth. However, the introduction of a scattered radiation filter permits the use of a single calibration curve for all field sizes and measurement depths. In this paper the quality control procedure developed for dosimetric verification of IMRT technique is reported. In particular a system of film dosimetry for the verification of a 6 MV photon beam has been implemented, with the introduction of the scattered radiation filter in the clinical practice that permits one to achieve an absolute dose determination with a global uncertainty within 3.4% (1 s.d.). The film has been calibrated to be used both in perpendicular and parallel configurations. The work also includes the characterization of the Elekta MLC. Ionimetric independent detectors have been used to check single point doses. The film dosimetry procedure has been applied to compare the measured absolute dose distributions with the ones calculated by the TPS, both for test and clinical plans. The agreement, quantified by the gamma index that seldom reaches the 1.5 value, is satisfying considering that the comparison is performed between absolute doses.     

Energy spectra, sources, and shielding considerations for neutrons generated by a flattening filter-free Clinac

Neutron production is an unwanted result of high-energy radiation therapy and results in secondary exposure of patients and radiation therapists to radiation. Recent studies have shown that delivering therapy using a standard medical accelerator with the flattening filter removed may reduce neutron fluence by nearly 70% over the course of prostate intensity-modulated radiation therapy (IMRT). In the current study, the 197Au Bonner sphere technique was used to compare the neutron spectrum produced when the filter is present and when it is absent. In addition, the following was calculated: (1) the neutron-shielding parameters of source strength and ambient dose equivalent (H0) and (2) using the Monte Carlo technique, the sources of neutron production in the accelerator head. It was found that the neutron spectrum was nearly constant, regardless of the presence of the flattening filter; however, the total fluence and ambient dose equivalent over the course of prostate IMRT were more than 70% lower when the filter was removed. Similarly, shielding parameters were lower when the filter was removed. Finally, the primary collimator and jaws accounted for the majority of neutron production, both with and without the flattening filter; however, with the flattening filter removed, the upper jaw accounted for much more neutron production relative to when the filter was present. Ultimately, removal of the flattening filter may offer several clinical advantages, including a reduction in the dose from neutrons to the patient and to radiation personnel.     

Photon scatter in intensity modulating filters evaluated by first Compton scatter and Monte Carlo calculations and experiments in broad beams

High-atomic-number materials may be used as intensity modulating filters for inverse radiation treatment planning with photon beams. Such filters, when placed in a bremsstrahlung beam, attenuate the primary fluence, but also produce scattered photons that will reach the patient. To account for such effects in the optimization of photon beam intensities a semiempirical method based on narrow and broad beam transmission measurements was used to quantify the number of scattered photons produced in these filters. The method was verified by performing analytical calculations based on first scatter and a Monte Carlo simulation in 6 and 18 MV photon beams. The resultant experimental transmission ratios agree with calculations by these methods within 2 per cent under the experimental conditions investigated. The semiempirical method can thus be used as a basis for preliminary decision-making to select the proper material for intensity modulating filters and can provide a fast method to perform independent quality checks of the calculation accuracy of dose planning systems. Change in beam penetration is of less concern when treatments of target volumes at smaller depths are of interest. A 10 g cm(-2) thick filter made of low-melting-point alloy produces a change in percentage depth dose of less than 2 per cent for depths larger than 10 cm independent of field size. Similarly the scatter correction modifies the dose distribution by less than 5-10 per cent in most cases.     

Design of a fast multileaf collimator for radiobiological optimized IMRT with scanned beams of photons, electrons, and light ions

Intensity modulated radiation therapy is rapidly becoming the treatment of choice for most tumors with respect to minimizing damage to the normal tissues and maximizing tumor control. Today, intensity modulated beams are most commonly delivered using segmental multileaf collimation, although an increasing number of radiation therapy departments are employing dynamic multileaf collimation. The irradiation time using dynamic multileaf collimation depends strongly on the nature of the desired dose distribution, and it is difficult to reduce this time to less than the sum of the irradiation times for all individual peak heights using dynamic leaf collimation [Svensson et al., Phys. Med. Biol. 39, 37-61 (1994)]. Therefore, the intensity modulation will considerably increase the total treatment time. A more cost-effective procedure for rapid intensity modulation is using narrow scanned photon, electron, and light ion beams in combination with fast multileaf collimator penumbra trimming. With this approach, the irradiation time is largely independent of the complexity of the desired intensity distribution and, in the case of photon beams, may even be shorter than with uniform beams. The intensity modulation is achieved primarily by scanning of a narrow elementary photon pencil beam generated by directing a narrow well focused high energy electron beam onto a thin bremsstrahlung target. In the present study, the design of a fast low-weight multileaf collimator that is capable of further sharpening the penumbra at the edge of the elementary scanned beam has been simulated, in order to minimize the dose or radiation response of healthy tissues. In the case of photon beams, such a multileaf collimator can be placed relatively close to the bremsstrahlung target to minimize its size. It can also be flat and thin, i.e., only 15-25 mm thick in the direction of the beam with edges made of tungsten or preferably osmium to optimize the sharpening of the penumbra. The low height of the collimator will minimize edge scatter from glancing incidence. The major portions of the collimator leafs can then be made of steel or even aluminum, so that the total weight of the multileaf collimator will be as low as 10 kg, which may even allow high-speed collimation in real time in synchrony with organ movements. To demonstrate the efficiency of this collimator design in combination with pencil beam scanning, optimal radiobiological treatments of an advanced cervix cancer were simulated. Different geometrical collimator designs were tested for bremsstrahlung, electron, and light ion beams. With a 10 mm half-width elementary scanned photon beam and a steel collimator with tungsten edges, it was possible to make as effective treatments as obtained with intensity modulated beams of full resolution, i.e., here 5 mm resolution in the fluence map. In combination with narrow pencil beam scanning, such a collimator may provide ideal delivery of photons, electrons, or light ions for radiation therapy synchronized to breathing and other organ motions. These high-energy photon and light ion beams may allow three-dimensional in vivo verification of delivery and thereby clinical implementation of the BioArt approach using Biologically Optimized three-dimensional in vivo predictive Assay based adaptive Radiation Therapy [Brahme, Acta Oncol. 42, 123-126 (2003)].     

Delivery time comparison for intensity-modulated radiation therapy with/without flattening filter: a planning study

The treatment delivery time of intensity-modulated radiation therapy (IMRT) with a multileaf collimator (MLC) is generally longer than that of conventional radiotherapy. In theory, removing the flattening filter from the treatment head may reduce the beam-on time by enhancing the output dose rate, and then reduce the treatment delivery time. And in practice, there is a possibility of delivering the required fluence distribution by modulating the unflattened non-uniform fluence distribution. However, the reduction of beam-on time may be discounted by the increase of leaf-travel time and (or) verification-and-recording (V&R) time. Here we investigate the overall effect of flattening filter on the treatment delivery time of IMRT with MLCs implemented in the step and shoot method, as well as with compensators on six hybrid machines. We compared the treatment delivery time with/without flattening filter for ten nasopharynx cases and ten prostate cases by observing the variations of the ratio of the beam-on time, segment number, leaf-travel time and the treatment delivery time with dose rate, leaf speed and V&R time. The results show that, without the flattening filter, the beam-on time reduces for both static MLC and compensator-based techniques: the number of segments and the leaf-travel time increase slightly for the static MLC technique; the relative IMRT treatment delivery time decreases more with lower dose rate, higher leaf speed and shorter V&R overhead time. The absolute treatment delivery time reduction depends on the fraction dose. It is not clinically significant at a fraction dose of 2 Gy for the technique of removing the flattening filter, but becomes significant when the fraction dose is as high as that for radiosurgery.     

Photon beam characterization and modelling for Monte Carlo treatment planning

Photon beams of 4, 6 and 15 MV from Varian Clinac 2100C and 2300C/D accelerators were simulated using the EGS4/BEAM code system. The accelerators were modelled as a combination of component modules (CMs) consisting of a target, primary collimator, exit window, flattening filter, monitor chamber, secondary collimator, ring collimator, photon jaws and protection window. A full phase space file was scored directly above the upper photon jaws and analysed using beam data processing software, BEAMDP, to derive the beam characteristics, such as planar fluence, angular distribution, energy spectrum and the fractional contributions of each individual CM. A multiple-source model has been further developed to reconstruct the original phase space. Separate sources were created with accurate source intensity, energy, fluence and angular distributions for the target, primary collimator and flattening filter. Good agreement (within 2%) between the Monte Carlo calculations with the source model and those with the original phase space was achieved in the dose distributions for field sizes of 4 cm x 4 cm to 40 cm x 40 cm at source surface distances (SSDs) of 80-120 cm. The dose distributions in lung and bone heterogeneous phantoms have also been found to be in good agreement (within 2%) for 4, 6 and 15 MV photon beams for various field sizes between the Monte Carlo calculations with the source model and those with the original phase space.     

蒙特卡罗方法研究均整器对6MVX线能谱的影响

目的：用蒙特卡罗方法计算Varian Trilogy加速器无均整器条件下6MVX射线能谱。分析均整器对能谱的影响。方法：先用BEAMnrc分别模拟计算Varian Trilogy加速器6MV射线在具备均整器和无均整器条件下,方野边长分别为4cm、6cm、8cm、10cm、20cm和40cm时的相空间文件。以相空间文件为输入用BEAMDP分析光子能谱。结果：无均整器条件下光子注量增大,在光子能谱峰值附近最明显．射野边长为4cm时去掉均整器后光子注量增加的最多为6．284倍,随着射野增大增加倍数减小,射野边长为40cm时最小为2．398倍．无均整器条件下光子谱峰值能量降低,光子谱整体左移,平均能量明显减小。结论：去除均整器后,加速器的输出光子能谱发生较大变化。随之剂量特性发生改变,临床上可能产生一定的获益或未知情况,尚需要进一步的研究支持。     

Modelling lateral beam quality variations in pencil kernel based photon dose calculations

Standard treatment machines for external radiotherapy are designed to yield flat dose distributions at a representative treatment depth. The common method to reach this goal is to use a flattening filter to decrease the fluence in the centre of the beam. A side effect of this filtering is that the average energy of the beam is generally lower at a distance from the central axis, a phenomenon commonly referred to as off-axis softening. The off-axis softening results in a relative change in beam quality that is almost independent of machine brand and model. Central axis dose calculations using pencil beam kernels show no drastic loss in accuracy when the off-axis beam quality variations are neglected. However, for dose calculated at off-axis positions the effect should be considered, otherwise errors of several per cent can be introduced. This work proposes a method to explicitly include the effect of off-axis softening in pencil kernel based photon dose calculations for arbitrary positions in a radiation field. Variations of pencil kernel values are modelled through a generic relation between half value layer (HVL) thickness and off-axis position for standard treatment machines. The pencil kernel integration for dose calculation is performed through sampling of energy fluence and beam quality in sectors of concentric circles around the calculation point. The method is fully based on generic data and therefore does not require any specific measurements for characterization of the off-axis softening effect, provided that the machine performance is in agreement with the assumed HVL variations. The model is verified versus profile measurements at different depths and through a model self-consistency check, using the dose calculation model to estimate HVL values at off-axis positions. A comparison between calculated and measured profiles at different depths showed a maximum relative error of 4% without explicit modelling of off-axis softening. The maximum relative error was reduced to 1% when the off-axis softening was accounted for in the calculations.     

The photon sensitivity of a moderated activation neutron detector

A method is described for estimating the overresponse produced by neutrons generated in the hydrocarbon moderator of a neutron dosimetry system by photon reactions. It is shown that photon interactions in the moderator produce an apparent neutron fluence per unit absorbed dose of 1.40 X 10(6) n/cm2/Gy-x for a flattened 33-MV x-ray beam and 2.66 X 10(6) n/cm2/Gy-x for the same beam without the flattening filter in place. The photon response of the moderated activation neutron detector was found to decrease rapidly with end point energy of the x-ray beam in the range 33 to 18.7 MeV.     

Spectral reconstruction of high energy photon beams for kernel based dose calculations

A kernel-based dose computation method with finite-size pencil beams (FSPBs) requires knowledge of the photon spectrum. Published methods of indirect spectral measurements using transmission measurements through beam attenuators use mathematical fits with a large number of parameters and constraints. In this study, we examine a simple strategy for fitting transmission data that models important physical characteristics of photon beams produced in clinical linear accelerators. The shape of an unattenuated bremsstrahlung spectrum is known, varying linearly from a maximum at zero energy to a value of zero at a maximum energy. This unattenuated spectrum is altered primarily by absorption of low energy photons by the flattening filter, causing the true spectrum to roll off to zero at low photon energies. A fitting equation models this behavior and has these advantages over previous methods: (1) the equation describes the shape of a bremsstrahlung spectrum based on physical expectations; and (2) only three fit parameters are required with a single constraint. Results for 4 MV and 6 MV accelerators for central axis and off-axis beams show good agreement with the maximum, average and modal energies for known spectra. Previously published models, representations of beam fluence (energy fluence, dN/dE), experimental methods, and the fitting process are discussed.     

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.     

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.