Monte Carlo simulations of the differential beam hardening effect of a flattening filter on a therapeutic x-ray beam

An optimized filter design approach for enhancing imaging quality in industrial linear accelerator

BACKGROUND

The polychromatic X-rays generated by a linear accelerator (Linac) often result in noticeable hardening artifacts in images, posing a significant challenge to accurate defect identification. To address this issue, a simple yet effective approach is to introduce filters at the radiation source outlet. However, current methods are often empirical, lacking scientifically sound metrics.

OBJECTIVE

This study introduces an innovative filter design method that optimizes filter performance by balancing the impact of ray intensity and energy on image quality.

MATERIALS AND METHODS

Firstly, different spectra under various materials and thicknesses of filters were obtained using GEometry ANd Tracking (Geant4) simulation. Subsequently, these spectra and their corresponding incident photon counts were used as input sources to generate different reconstructed images. By comprehensively comparing the intensity differences and noise in images of defective and non-defective regions, along with considering hardening indicators, the optimal filter was determined.

RESULTS

The optimized filter was applied to a Linac-based X-ray computed tomography (CT) detection system designed for identifying defects in graphite materials within high-temperature gas-cooled reactor (HTR), with defect dimensions of 2 mm. After adding the filter, the hardening effect reduced by 22%, and the Defect Contrast Index (DCI) reached 3.226.

CONCLUSION

The filter designed based on the parameters of Average Difference (AD) and Defect Contrast Index (DCI) can effectively improve the quality of defect images.

Characterization of the radiation beam of a tomotherapy equipment with MCNP

This work aims to model and characterize the radiation beam of one Accuray tomotherapy equipment using the Monte Carlo Code MCNP5 (Monte Carlo N-Particle). This tomotherapy equipment is used for delivering high doses of radiation in tumor regions to kill cancer cells and shrink the tumor during radiation therapy of cancer patients, however, the radiation can damage surrounding areas and nearby organs at risk (OAR) if the radiation field is not well delimited. In particular, intensity-modulated radiotherapy treatments (IMRT) with tomotherapy equipment offer great benefits to patients allowing treatment of tumor regions without affecting surrounding areas and OAR. Nowadays, it is well known that a correct simulation of transport of radiation in tomotherapy equipment facilitates considerably the estimation of ideal doses in the tumor, surrounding regions, and OAR. For that reason, in this work, we simulated the geometry of the 6 MV ACCURAY Tomotherapy equipment of the CECAN using the MCNP5. The model includes a TomoLINAC consisting of an electron source that emits Gaussian distribution particles with an average energy of 5.7 MeV and width of 0.3 MeV. The emitted particles impact the tungsten target and pass through primary collimators and jaws that define the irradiation field in the isocenter. To validate the geometry and radiation transport in the TomoLINAC the curves of depth dose percentage (PDD) estimated by simulation and the curves measured experimentally were tuned. In the same way, the simulated transverse and longitudinal profiles were compared with the experimental results. In addition, a comparison between the qualities of the radiation beam characterized with MCNP and measured experimentally in CECAN showed a deviation of 1%. For the simulations, cylindrical detectors located inside a water phantom were considered and it was employed the tally *F8. A good agreement was observed between the PDD's curves obtained from the simulation and those measured experimentally for a field of 5 × 10 cm2 in the isocenter and SSD (distance from the source to the surface) of 85 cm. Also, the comparison between the simulated and experimental transverse profiles obtained at 1.5 cm, 10 cm and 15 cm depth with a radiation field of 5 × 40 cm2 showed very good agreement. The longitudinal profiles were estimated with the same depths as the transverse ones, but for each of them, the openings of the jaws were 5.0 cm, 2.5 cm and 1.0 cm in the longitudinal direction, which corresponds to the direction in which the patient's table moves. The comparison between the simulated and experimental longitudinal profiles showed good concordance too. Once the radiation beam of the ACCURAY tomotherapy equipment had been characterized, experimental dose measurements were made using a Cheese phantom and two A1SL ionization chambers. These results obtained experimentally were compared with those estimated with MCNP for a field of 5 × 40 cm2 at the isocenter and SAD of 85 cm and, it was concluded that both results were similar considering the regions of uncertainty. Finally, we must highlight that the modeling and characterization of the radiation beam of CECAN's ACCURAY tomotherapy equipment can be a key tool for dose estimations in different cancer treatment plans and future research.

Study of Dosimetric Properties of Flattening Filter Free Photon Beam Passing through Cadmium Free Compensator Alloy

BACKGROUND

This study aims to investigate radiation beam geometry of Cyberknife beam and change in dosimetric characteristics of six megavoltage (6MV) flattening filter free (FFF) beam after passing through high density cadmium free compensator alloy.

MATERIAL AND METHODS

In this experimental study, changes in FFF beam dosimetric characteristics after passing through compensator alloy was measured. Transmitted intensity of FFF beam was measured in air by an ion chamber at a source to detector distance (SDD) of 800mm. Extended SDD measurement also has been performed at a distance of 1270mm to analyze scattering due to compensator. Linear attenuation coefficient (µ_eff) was measured for cadmium free compensator alloy using simple exponential attenuation model. Percentage depth doses (PDDs) have been measured by a radiation field analyzer with compensator material to observe the beam hardening and change in surface doses and depth doses.

RESULTS

Linear attenuation coefficient of compensator alloy was measured 0.042 (Standard Deviation ±0.00099) mm^-1 and it was found that there is no change with increase in collimator size. Even after increasing distance source from detector, µ_eff has no change. PDDs were found to increase with thickness of compensator. PDD from a 60mm collimator size increased by 5% and 6% at a depth of 100mm and 200mm, respectively in water. PDD also increased with collimator size less significantly. Surface dose was found to decrease with increase in compensator thickness.

CONCLUSION

Cyberknife beam has been found to be narrow beam geometry. FFF beam contains lesser scattered photons. Presence of high density compensator filters out the soft x-ray photon causes significant dosimetric changes.

Monte Carlo Study of Unflattened Photon Beams Shaped by Multileaf Collimator

INTRODUCTION

This study investigates basic dosimetric properties of unflattened 6 MV photon beam shaped by multileaf collimator and compares them with those of flattened beams.

MATERIALS AND METHODS

Monte Carlo simulation model using BEAM code was developed for a 6MV photon beam based on Varian Clinic 600 unique performance linac operated with and without a flattening filter in beam line. Dosimetric features including lateral profiles, central axis depth dose, photon and electron spectra were calculated for flattened and unflattened cases, separately.

RESULTS

An increase in absolute depth dose with a factor of more than 2.4 was observed for unflattened beam which was dependent on depth. PDDs values were found to be lower for unflattened beam for all field sizes. Significant decrease in calculated mlc leakage was observed when the flattening filter was removed from the beam line. The total scatter factor, S_CP was found to show less variation with field sizes for unflattened beam indicating a decrease in head scatter. The beam profiles for unflattened case are found to have lower relative dose value in comparison with flattened beam near the field edge, and it falls off faster with distance.

CONCLUSION

Our study showed that increase in the dose rate and lower peripheral dose could be considered as realistic advantages for unflattened 6MV photon beams.

Energy spectrum and dose enhancement due to the depth of the Lipiodol position using flattened and unflattened beams

Aim

Lipiodol was used for stereotactic body radiotherapy combining trans arterial chemoembolization. Lipiodol used for tumour seeking in trans arterial chemoembolization remains in stereotactic body radiation therapy. In our previous study, we reported the dose enhancement effect in Lipiodol with 10× flattening-filter-free (FFF). The objective of our study was to evaluate the dose enhancement and energy spectrum of photons and electrons due to the Lipiodol depth with flattened (FF) and FFF beams.

Methods

FF and FFF for 6 MV beams from TrueBeam were used in this study. The Lipiodol (3 × 3 × 3 cm³) was located at depths of 1, 3, 5, 10, 20, and 30 cm in water. The dose enhancement factor (DEF) and the energy fluence were obtained by Monte Carlo calculations of the particle and heavy ion transport code system (PHITS).

Results

The DEFs at the centre of Lipiodol with the FF beam were 6.8, 7.3, 7.6, 7.2, 6.1, and 5.7% and those with the FFF beam were 20.6, 22.0, 21.9, 20.0, 12.3, and 12.1% at depths of 1, 3, 5, 10, 20, and 30 cm, respectively, where Lipiodol was located in water. Moreover, spectrum results showed that more low-energy photons and electrons were present at shallow depth where Lipiodol was located in water. The variation in the low-energy spectrum due to the depth of the Lipiodol position was more explicit with the FFF beam than that with the FF beam.

Conclusions

The current study revealed variations in the DEF and energy spectrum due to the depth of the Lipiodol position with the FF and FFF beams. Although the FF beam could reduce the effect of energy dependence due to the depth of the Lipiodol position, the dose enhancement was overall small. To cause a large dose enhancement, the FFF beam with the distance of the patient surface to Lipiodol within 10 cm should be used.

Study of unflattened photon beams shaped by multileaf collimator using BEAMnrc code

Purpose

In our study basic dosimetric properties of a flattening filter free 6 MV photon beam shaped by multileaf collimators (MLC) is examined using the Monte Carlo (MC) method.

Methods and Materials

BEAMnrc code was used to make a MC simulation model for 6 MV photon beam based on Varian Clinic 600 unique performance linac, operated with and without a flattening filter in beam line. Dosimetric features including central axis depth dose, beam profiles, photon and electron spectra were calculated and compared for flattened and unflattened cases.

Results

Dosimetric field size and penumbra were found to be smaller for unflattened beam, and the decrease in field size was less for MLC shaped in comparison with jaw-shaped unflattened beam. Increase in dose rate of >2·4 times was observed for unflattened beam indicating a shorter beam delivery time for treatment. MLC leakage was found to decrease significantly when the flattening filter was removed from the beam line. The total scatter factor showed slower deviation with field sizes for unflattened beam indicating a reduced head scatter.

Conclusions

Our study demonstrated that improved accelerator characteristics can be achieved by removing flattening filter from beam line.

Measurement and comparison of head scatter factor for 7 MV unflattened (FFF) and 6 MV flattened photon beam using indigenously designed columnar mini phantom

AIM

To measure and compare the head scatter factor for 7 MV unflattened and 6 MV flattened photon beam using a home-made designed mini phantom.

BACKGROUND

The head scatter factor (Sc) is one of the important parameters for MU calculation. There are multiple factors that influence the Sc values, like accelerator head, flattening filter, primary and secondary collimators.

MATERIALS AND METHODS

A columnar mini phantom was designed as recommended by AAPM Task Group 74 with high and low atomic number material for measurement of head scatter factors at 10 cm and d max dose water equivalent thickness.

RESULTS

The Sc values measured with high-Z are higher than the low-Z mini phantoms observed for both 6MV-FB and 7MV-UFB photon energies. Sc values of 7MV-UFB photon beams were smaller than those of the 6MV-FB photon beams (0.6-2.2% (Primus), 0.2-1.4% (Artiste) and 0.6-3.7% (Clinac iX (2300CD))) for field sizes ranging from 10 cm × 10 cm to 40 cm × 40 cm. The SSD had no influence on head scatter for both flattened and unflattened beams. The presence of wedge filters influences the Sc values. The collimator exchange effects showed that the opening of the upper jaw increases Sc irrespective of FF and FFF.

CONCLUSIONS

There were significant differences in Sc values measured for 6MV-FB and unflattened 7MV-UFB photon beams over the range of field sizes from 10 cm × 10 cm to 40 cm × 04 cm. Different results were obtained for measurements performed with low-Z and high-Z mini phantoms.

Comparison of Head Scatter Factor for 6MV and 10MV flattened (FB) and Unflattened (FFF) Photon Beam using indigenously Designed Columnar Mini Phantom

To measure and compare the head scatter factor for flattened (FB) and unflattened (FFF) of 6MV and 10MV photon beam using indigenously designed mini phantom. A columnar mini phantom was designed as recommended by AAPM Task Group 74 with low and high atomic number materials at 10 cm (mini phantom) and at approximately twice the depth of maximum dose water equivalent thickness (brass build-up cap). Scatter in the accelerator (Sc) values of 6MV-FFF photon beams are lesser than that of the 6MV-FB photon beams (0.66-2.8%; Clinac iX, 2300CD) and (0.47-1.74%; True beam) for field sizes ranging from 10 × 10 cm² to 40 × 40 cm². Sc values of 10MV-FFF photon beams are lesser (0.61-2.19%; True beam) than that of the 10MV-FB photons beams for field sizes ranging from 10 × 10 cm² to 40 × 40 cm². The SSD had no influence on head scatter for both flattened and unflattened beams and irrespective of head design of the different linear accelerators. The presence of field shaping device influences the Sc values. The collimator exchange effect reveals that the opening of the upper jaw increases Sc irrespective of FB or FFF photon beams and different linear accelerators, and it is less significant in FFF beams. Sc values of 6MV-FB square field were in good agreement with that of AAPM, TG-74 published data for Varian (Clinac iX, 2300CD) accelerator. Our results confirm that the removal of flattening filter decreases in the head scatter factor compared to flattened beam. This could reduce the out-of-field dose in advanced treatment delivery techniques.

Photon energy-modulated radiotherapy: Monte Carlo simulation and treatment planning study

PURPOSE

To demonstrate the feasibility of photon energy-modulated radiotherapy during beam-on time.

METHODS

A cylindrical device made of aluminum was conceptually proposed as an energy modulator. The frame of the device was connected with 20 tubes through which mercury could be injected or drained to adjust the thickness of mercury along the beam axis. In Monte Carlo (MC) simulations, a flattening filter of 6 or 10 MV linac was replaced with the device. The thickness of mercury inside the device varied from 0 to 40 mm at the field sizes of 5 × 5 cm(2) (FS5), 10 × 10 cm(2) (FS10), and 20 × 20 cm(2) (FS20). At least 5 billion histories were followed for each simulation to create phase space files at 100 cm source to surface distance (SSD). In-water beam data were acquired by additional MC simulations using the above phase space files. A treatment planning system (TPS) was commissioned to generate a virtual machine using the MC-generated beam data. Intensity modulated radiation therapy (IMRT) plans for six clinical cases were generated using conventional 6 MV, 6 MV flattening filter free, and energy-modulated photon beams of the virtual machine.

RESULTS

As increasing the thickness of mercury, Percentage depth doses (PDD) of modulated 6 and 10 MV after the depth of dose maximum were continuously increased. The amount of PDD increase at the depth of 10 and 20 cm for modulated 6 MV was 4.8% and 5.2% at FS5, 3.9% and 5.0% at FS10 and 3.2%-4.9% at FS20 as increasing the thickness of mercury from 0 to 20 mm. The same for modulated 10 MV was 4.5% and 5.0% at FS5, 3.8% and 4.7% at FS10 and 4.1% and 4.8% at FS20 as increasing the thickness of mercury from 0 to 25 mm. The outputs of modulated 6 MV with 20 mm mercury and of modulated 10 MV with 25 mm mercury were reduced into 30%, and 56% of conventional linac, respectively. The energy-modulated IMRT plans had less integral doses than 6 MV IMRT or 6 MV flattening filter free plans for tumors located in the periphery while maintaining the similar quality of target coverage, homogeneity, and conformity.

CONCLUSIONS

The MC study for the designed energy modulator demonstrated the feasibility of energy-modulated photon beams available during beam-on time. The planning study showed an advantage of energy-and intensity modulated radiotherapy in terms of integral dose without sacrificing any quality of IMRT plan.

Dosimetric properties of a beam quality-matched 6 MV unflattened photon beam

The purpose of this study was to report the characteristics of an equivalent quality unflattened (eqUF) photon beam in clinical implementation and to provide a generalized method to describe unflattened (UF) photon beam profiles. An unflattened photon beam with a beam quality equivalent to the corresponding flat 6 MV photon beam (WF) was obtained by removing the flattening filter from a Siemens ONCOR Avant‐Garde linear accelerator and adjusting the photon energy. A method independent from the WF beam profile was presented to describe UF beam profiles and other selected beam characteristics were examined. The short‐term beam stability was examined by dynamic beam profiles, recorded every 0.072 s in static and gated delivery, and the long‐term stability was evidenced by the five‐year clinical quality assurance records. The dose rate was raised fivefold using the eqUF beam. The depth of maximum dose (dmax) shifted 3 mm deeper, but the percent depth dose beyond dmax was very similar to that of the WF beam. The surface dose and out‐of‐field dose were lower, but the penumbra was slightly wider. The variation in head scatter and phantom scatter with changes in field size was smaller; the variation in the profile shape with change in depth was also smaller. The eqUF beam is stable 0.072 s after the beam is turned on, and the five‐year beam stability was comparable to that of the WF beam. A fivefold dose rate increase was observed in the eqUF beam with similar beam characteristics to other reported UF beam data except for a deeper dmax and a slightly wider penumbra. The initial and long‐term stability of the eqUF beam profile is on parity with the WF beam. The UF beam profile can be described in the generalized method independently without relying on the WF beam profile.

PACS number: 87.55.de

Assessment of computerized treatment planning system accuracy in calculating wedge factors of physical wedged fields for 6 MV photon beams

Wedge filters are commonly used in external beam radiotherapy to achieve a uniform dose distribution within the target volume. The main objective of this study was to investigate the accuracy of the beam modifier algorithm of Theraplan plus (TPP version 3.8) treatment planning system and to confirm that either the beam hardening, beam softening and attenuation coefficients along with wedge geometry and measured wedge factor at single depth and multiple fields sizes can be the replacement of wedged profile and wedged cross-sectional data or not. In this regard the effect of beam hardening and beam softening was studied with physical wedges for 6 MV photons. The Normalized Wedge Factors (NWFs) were measured experimentally as well as calculated with the Theraplan plus, as a function of depth and field size in a water phantom for 15°, 30°, 45°, and 60° wedge filters. The beam hardening and softening was determined experimentally by deriving the required coefficients for all wedge angles. The TPP version 3.8 requires wedge transmission factor at single depth and multiple field sizes. Without incorporating the hardening and softening coefficients the percent difference between measured and calculated NFWs was as high as 7%. After the introduction of these parameters into the algorithm, the agreement between measured and TPP (V 3.8) calculated NWFs were improved to within 2 percent for various depths. Similar improvement was observed in TPP version 3.8 while calculating NWFs for various field sizes when the required coefficients were adjusted. In conclusion, the dose calculation algorithm of TPP version 3.8 showed good accuracy for a 6 MV photon beam provided beam hardening and softening parameters are taken into account. From the results, it is also concluded that, the beam hardening, beam softening and attenuation coefficients along with wedge geometry and measured wedge factor at single depth and multiple fields sizes can be the replacement of wedged profile and wedged cross-sectional data in the TPS. The study also indicated that by ignoring the beam softening and beam hardening will result in an inaccurate dose to the target volume of the patient.

Characterization of the theorectical radiation dose enhancement from nanoparticles

Recently, nanoparticles have been considered as a method of providing radiation dose enhancement in tumors. In order to quantify this affect, a dose enhancement factor (DEF) is defined that represents the ratio of the dose deposited in tumor with nanoparticles, divided by the dose deposited in the tumor without nanoparticles. Materials with atomic numbers (Z) ranging from 25 to 90 are considered in this analysis. In addition, the energy spectrum for a number of external beam x-ray sources and common radionuclides are evaluated. For a nanoparticle concentration of 5 mg/ml, the DEF is < 1.05 for Co-60, Ir-192, Au-198, Cs-137, 6, 18, and 25 MV x-rays for all materials considered. However, relatively large increases in the DEF are observed for 50, 80, 100, and 140 KVp x-rays as well as Pd-103 and I-125. The DEF increases for all sources as Z varies from 25-35. From Z = 40-60, the DEF plateaus or slightly decreases. For higher Z materials (Z>70), the DEF increases and is a maximum for the highest Z materials. High atomic number nanoparticles coupled with low energy external beam x-rays or brachytherapy sources offer the potential of significantly enhancing the delivered dose.

Multibeam tomotherapy: a new treatment unit devised for multileaf collimation, intensity-modulated radiation therapy

A fully integrated system for treatment planning, application, and verification for automated multileaf collimator (MLC) based, intensity-modulated, image-guided, and adaptive radiation therapy (IMRT, IGRT and ART, respectively) is proposed. Patient comfort, which was the major development goal, will be achieved through a new unit design and short treatment times. Our device for photon beam therapy will consist of a new dual energy linac with five fixed treatment heads positioned evenly along one plane but one electron beam generator only. A minimum of moving parts increases technical reliability and reduces motion times to a minimum. Motion is allowed solely for the MLCs, the robotic patient table, and the small angle gantry rotation of +/- 36 degrees. Besides sophisticated electron beam guidance, this compact setup can be built using existing modules. The flattening-filter-free treatment heads are characterized by reduced beam-on time and contain apertures restricted in one dimension to the area of maximum primary fluence output. In the case of longer targets, this leads to a topographic intensity modulation, thanks to the combination of "step and shoot" MLC delivery and discrete patient couch motion. Owing to the limited number of beam directions, this multislice cone beam serial tomotherapy is referred to as "multibeam tomotherapy." Every patient slice is irradiated by one treatment head at any given moment but for one subfield only. The electron beam is then guided to the next head ready for delivery, while the other heads are preparing their leaves for the next segment. The "Multifocal MLC-positioning" algorithm was programmed to enable treatment planning and optimize treatment time. We developed an overlap strategy for the longitudinally adjacent fields of every beam direction, in doing so minimizing the field match problem and the effects of possible table step errors. Clinical case studies show for the same or better planning target volume coverage, better organ-at-risk sparing, and comparable mean integral dose to the normal tissue a reduction in treatment time by more than 50% to only a few minutes in comparison to high-quality 3-D conformal and IMRT treatments. As a result, it will be possible to incorporate features for better patient positioning and image guidance, while sustaining reasonable overall treatment times at the same time. The virtual multibeam tomotherapy design study TOM'5-CT contains a dedicated electron beam CT (TOM'AGE) and an objective optical topometric patient positioning system (TOPOS). Thanks to the wide gantry bore of 120 cm and slim gantry depths of 70 cm, patients can be treated very comfortably, in all cases tumor-isocentrically, as well as with noncoplanar beam arrangements as in stereotactic radiosurgery with a couch rotation of up to +/- 54 degrees. The TOM'5 treatment unit on which this theoretical concept is based has a stand-alone depth of 40 cm and an outer diameter of 245 cm; the focus-isocenter distance of the heads is 100 cm with a field size of 40 cm x 7 cm and 0.5 cm leaves, which operate perpendicular to the axis of table motion.

Dosimetric characteristics of unflattened 6MV photon beams of a clinical linear accelerator: A Monte Carlo study

The purpose of this study is to evaluate the dosimetric properties of a flattening filter free 6 MV photon beam. The 6 MV photon beam of a Varian Clinac 21 EX linac was modeled using the MCNP4C Monte Carlo (MC) code. Dosimetric features including central axis absorbed doses, beam profiles and photon energy spectra were calculated for flattened and unflattened 6 MV photon beams. A substantial increase in the dose rate was seen for the unflattened beam, which was decreased with field size and depth. The penumbra width was decreased less than 0.2 mm (about 5%) and a 25% decrease in out-of-field dose was observed for the unflattened beam. The photon energy spectra were softer for the unflattened beam and the mean energies of spectra were higher for smaller field size. Our study showed that increase in the dose rate and lower out-of-field dose could be considered as practical advantages for unflattened 6 MV beams.

Dosimetric properties of a flattening filter-free 6-MV photon beam: a Monte Carlo study

PURPOSE

The dosimetric features of an unflattened 6-MV photon beam of an Elekta SL-25 linac was calculated by the Monte Carlo (MC) method.

MATERIAL AND METHODS

The head of the Elekta SL-25 linac was simulated using the MCNP4C MC code. The accuracy of the model was evaluated using measured dosimetric features, including depth dose values and dose profiles in a water phantom. The flattening filter was then removed, and beam dosimetric properties were calculated by the MC method and compared with those of the flattened photon beam.

RESULTS

Our results showed a significant (twofold) increase in the dose rate for all field sizes. Also, the photon beam spectra for an unflattened beam were softer, which led to a steeper reduction in depth doses. The decrease in the out-of-field dose and increase in the contamination electrons and a buildup region dose were the other consequences of removing the flattening filter.

CONCLUSION

Our study revealed that, for recent radiotherapy techniques, the use of multileaf collimators for beam shaping removing the flattening filter could offer some advantages, including an increased dose rate and decreased out-of-field dose.

Optimization of photon beam flatness for radiation therapy

In this work, we investigate the relation between lateral fluence/dose distributions and photon beam uniformity, possibly identifying ways to improve these characteristics. The calculations included treatment head scatter properties associated with three common types of linear accelerators in order to study their impact on the results. For 6 and 18 MV photon beams the lateral fluence distributions were optimized with respect to the resulting calculated flatness, as defined by the International Electrotechnical Commission (IEC), at 10 cm depth in six different field sizes. The limits proposed by IEC for maximum dose ratios ('horns') at the depth of dose maximum have also been accounted for in the optimization procedure. The conclusion was that typical head scatter variations among different types of linear accelerators have a very limited effect on the optimized results, which implies that the existing differences in measured off-axis dose distributions are related to non-equivalent optimization objectives. Finally, a comparison between the theoretically optimized lateral dose distributions and corresponding dose measurements for the three investigated accelerator types was performed. Although the measured data generally fall within the IEC requirements the optimized distributions show better results overall for the evaluated uniformity parameters, indicating that there is room for improved flatness performance in clinical photon beams.

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.

Dose uncertainties in photon pencil kernel calculations at off-axis positions

The purpose of this study was to investigate the specific problems associated with photon dose calculations in points located at a distance from the central beam axis. These problems are related to laterally inhomogeneous energy fluence distributions and spectral variations causing a lateral shift in the beam quality, commonly referred to as off-axis softening (OAS). We have examined how the dose calculation accuracy is affected when enabling and disabling explicit modeling of these two effects. The calculations were performed using a pencil kernel dose calculation algorithm that facilitates modeling of OAS through laterally varying kernel properties. Together with a multi-source model that provides the lateral energy fluence distribution this generates the total dose output, i.e., the dose per monitor unit, at an arbitrary point of interest. The dose calculation accuracy was evaluated through comparisons with 264 measured output factors acquired at 5, 10, and 20 cm depth in four different megavoltage photon beams. The measurements were performed up to 18 cm from the central beam axis, inside square fields of varying size and position. The results show that calculations including explicit modeling of OAS were considerably more accurate, up to 4%, than those ignoring the lateral beam quality shift. The deviations caused by simplified head scatter modeling were smaller, but near the field edges additional errors close to 1% occurred. When enabling full physics modeling in the dose calculations the deviations display a mean value of -0.1%, a standard deviation of 0.7%, and a maximum deviation of -2.2%. Finally, the results were analyzed in order to quantify and model the inherent uncertainties that are present when leaving the central beam axis. The off-axis uncertainty component showed to increase with both off-axis distance and depth, reaching 1% (1 standard deviation) at 20 cm depth.

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.

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.

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.

ORANGE: a Monte Carlo dose engine for radiotherapy

This study presents data for the verification of ORANGE, a fast MCNP-based dose engine for radiotherapy treatment planning. In order to verify the new algorithm, it has been benchmarked against DOSXYZ and against measurements. For the benchmarking, first calculations have been done using the ICCR-XIII benchmark. Next, calculations have been done with DOSXYZ and ORANGE in five different phantoms (one homogeneous, two with bone equivalent inserts and two with lung equivalent inserts). The calculations have been done with two mono-energetic photon beams (2 MeV and 6 MeV) and two mono-energetic electron beams (10 MeV and 20 MeV). Comparison of the calculated data (from DOSXYZ and ORANGE) against measurements was possible for a realistic 10 MV photon beam and a realistic 15 MeV electron beam in a homogeneous phantom only. For the comparison of the calculated dose distributions and dose distributions against measurements, the concept of the confidence limit (CL) has been used. This concept reduces the difference between two data sets to a single number, which gives the deviation for 90% of the dose distributions. Using this concept, it was found that ORANGE was always within the statistical bandwidth with DOSXYZ and the measurements. The ICCR-XIII benchmark showed that ORANGE is seven times faster than DOSXYZ, a result comparable with other accelerated Monte Carlo dose systems when no variance reduction is used. As shown for XVMC, using variance reduction techniques has the potential for further acceleration. Using modern computer hardware, this brings the total calculation time for a dose distribution with 1.5% (statistical) accuracy within the clinical range (less then 10 min). This means that ORANGE can be a candidate for a dose engine in radiotherapy treatment planning.

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.

An MCNP-based model of a medical linear accelerator x-ray photon beam

The major components in the x-ray photon beam path of the treatment head of the VARIAN Clinac 2300 EX medical linear accelerator were modeled and simulated using the Monte Carlo N-Particle radiation transport computer code (MCNP). Simulated components include x-ray target, primary conical collimator, x-ray beam flattening filter and secondary collimators. X-ray photon energy spectra and angular distributions were calculated using the model. The x-ray beam emerging from the secondary collimators were scored by considering the total x-ray spectra from the target as the source of x-rays at the target position. The depth dose distribution and dose profiles at different depths and field sizes have been calculated at a nominal operating potential of 6 MV and found to be within acceptable limits. It is concluded that accurate specification of the component dimensions, composition and nominal accelerating potential gives a good assessment of the x-ray energy spectra.

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.

Monte Carlo modelling of external radiotherapy photon beams

An essential requirement for successful radiation therapy is that the discrepancies between dose distributions calculated at the treatment planning stage and those delivered to the patient are minimized. An important component in the treatment planning process is the accurate calculation of dose distributions. The most accurate way to do this is by Monte Carlo calculation of particle transport, first in the geometry of the external or internal source followed by tracking the transport and energy deposition in the tissues of interest. Additionally, Monte Carlo simulations allow one to investigate the influence of source components on beams of a particular type and their contaminant particles. Since the mid 1990s, there has been an enormous increase in Monte Carlo studies dealing specifically with the subject of the present review, i.e., external photon beam Monte Carlo calculations, aided by the advent of new codes and fast computers. The foundations for this work were laid from the late 1970s until the early 1990s. In this paper we will review the progress made in this field over the last 25 years. The review will be focused mainly on Monte Carlo modelling of linear accelerator treatment heads but sections will also be devoted to kilovoltage x-ray units and 60Co teletherapy sources.

[Improvement of the accuracy of the Laplace transform method for the determination of radiotherapy spectra of clinical linear accelerators]

The present study focused on the reconstruction of the bremsstrahlung spectrum of a clinical linear accelerator from the measured transmission curve, with the aim of improving the accuracy of this method. The essence of the method was the analytic inverse Laplace transform of a parameter function fitted to the measured transmission curve. We tested known fitting functions, however they resulted in considerable fitting inaccuracy, leading to inaccuracies of the bremsstrahlung spectrum. In order to minimise the fitting errors, we employed a linear combination of n equations with 2n-1 parameters. The fitting errors are now considerably smaller. The measurement of the transmission function requires that the energy-dependent detector response is taken into account. We analysed the underlying physical context and developed a function that corrects for the energy-dependent detector response. The factors of this function were experimentally determined or calculated from tabulated values.

A Monte Carlo study on internal wedges using BEAM

To do calculations for wedged photon beams with the NRC Monte Carlo simulation package BEAM, a new Component Module for wedges called WEDGE has been designed and built. After an initial series of benchmarks using monoenergetic photon beams as well as realistic 6 MV and 10 MV beams, it was found, that the new CM did work fine for the large wedge (maximum field size 30 x 40 cm2) of the Elekta SL-linac. The next step was to calculate dose distributions and output factors for a range of wedged fields with field size from 3 x 3 cm2 to 30 x 30 cm2. Results from these simulations have been compared to measurements. Calculated values for the reference wedge transmission factor and the relative wedge transmission factors were within 1.5% from the measured data. Dose distributions showed an identical behavior; both depth-dose curves as well as cross profiles were within 1.5% from measured data, usually even better. Despite the increased mean energy, there was no indication that, as a result, the phantom scatter output factors will change for a 10 MV photon beam. It was found that by adding a wedge the contributions for the different sources of head scatter changed considerably as compared to the open fields, apart from the additional scatter from the wedge. Another consequence of inserting a wedge was an increase in the mean energy of both primary and scattered radiation with 0.3 MV and 0.7 MV, respectively, for all wedged fields with respect to the corresponding open fields. Despite the statistical uncertainty in the calculated data, which is in the same order of magnitude as the effect to be determined, it was possible to derive reliable data for the beam hardening from the calculated dose distributions. Only for the smallest field (field size 3 x 3 cm2) a large difference between the measured and calculated beam hardening factor was found due to the relative large voxel size of 1 x 1 x 1 cm3 compared to the field size. For a description of the influence of a wedge on a photon beam, the results of this study strongly support the use of a reference wedge transmission factor (determined under reference conditions) in combination with a relative wedge transmission factor. The product of these variables should replace the collimator scatter output factor used in open fields. The influence on the dose distribution should be incorporated by using the (field size dependent) beam hardening. The ultimate solution will be to make this beam hardening depending on the actual position in the radiation field, as the photon energy varies over the field (holds also for open fields).

Monte Carlo modelling of electron beams from medical accelerators

Monte Carlo simulation of radiation transport is considered to be one of the most accurate methods of radiation therapy dose calculation. With the rapid development of computer technology, Monte Carlo based treatment planning for radiation therapy is becoming practical. A basic requirement for Monte Carlo treatment planning is a detailed knowledge of the radiation beams from medical accelerators. A practical approach to obtain the above is to perform Monte Carlo simulation of radiation transport in the medical accelerator. Additionally, Monte Carlo modelling of the treatment machine head can also improve our understanding of clinical beam characteristics, help accelerator design and improve the accuracy of clinical dosimetry by providing more realistic beam data. This paper summarizes work over the past two decades on Monte Carlo simulation of clinical electron beams from medical accelerators.

Dose calculations for external photon beams in radiotherapy

Dose calculation methods for photon beams are reviewed in the context of radiation therapy treatment planning. Following introductory summaries on photon beam characteristics and clinical requirements on dose calculations, calculation methods are described in order of increasing explicitness of particle transport. The simplest are dose ratio factorizations limited to point dose estimates useful for checking other more general, but also more complex, approaches. Some methods incorporate detailed modelling of scatter dose through differentiation of measured data combined with various integration techniques. State-of-the-art methods based on point or pencil kernels, which are derived through Monte Carlo simulations, to characterize secondary particle transport are presented in some detail. Explicit particle transport methods, such as Monte Carlo, are briefly summarized. The extensive literature on beam characterization and handling of treatment head scatter is reviewed in the context of providing phase space data for kernel based and/or direct Monte Carlo dose calculations. Finally, a brief overview of inverse methods for optimization and dose reconstruction is provided.

Calculating photon beam characteristics with Monte Carlo techniques

This study describes the results from a simulation of a 10 MV photon beam from a medical linac using the BEAM code. To check the quality of the generated photon beam, the characteristics of this beam (depth dose curve, cross profiles, and output factors) have been calculated and compared to measured data. By splitting up the radiation head in two parts, the target section and the collimator section calculation times were long, but acceptable when aiming at phase space files containing some 5 million particles. Given the number of particles evaluated, the accuracy of all data was around 2%. Analysis of the phase space files for different field size supports results from previous studies about contaminant particles and sources for scattered radiation for photon beams from medical linacs. The total scatter output factor Scp as well as the collimator scatter output factor Sc have been calculated within 2% of measurements. Also, the ratio between dose at a reference point for the full scatter situation and the no-scatter situation has been calculated correctly. All depth dose curves and cross profiles have also been calculated correctly, although with only moderate statistics. Improvements are possible by increasing the number of particles in the simulations (up to 50 million for the largest field size) at least 4-8 times, although calculation times will increase with the same factor. Nevertheless, the method proved itself as reliable. Still, the accuracy should be improved to 1% or better. This is necessary as we plan to use Monte Carlo simulations to benchmark three-dimensional radiotherapy planning systems. By increasing the number of particles in the phase space files and subsequently increasing the number of particles in each simulation, this 1% accuracy will be achieved. The easy way to increase the number of particles in a simulation by increasing the number of times phase space files, which were already recycled ten times, are reused from ten times (this study) to 40 times or more will not work, as it introduces artifacts, especially in the cross profiles.

Analytical considerations of beam hardening in medical accelerator photon spectra

Beam hardening is a well-known phenomenon for therapeutic accelerator beams passing through matter in narrow beam geometry. This study assesses quantitatively the magnitude of beam hardening of therapeutic beams in water. A formal concept of beam hardening is proposed which is based on the decrease of the mean attenuation coefficient with depth. On the basis of this concept calculations of beam hardening effects are easily performed by means of a commercial spreadsheet program. Published accelerator spectra and the tabulated values of attenuation coefficients serve as input for these calculations. It is shown that the mean attenuation coefficient starts at depth zero with an almost linear decrease and then slowly levels off to a limit value. A similar behavior is found for the beam hardening coefficient. A physically reasonable, semianalytical model is given which fits the data better than previously published functions. The energy dependence of the initial attenuation coefficient is evaluated and shown. It fits well to published experimental data. The initial beam hardening coefficient, however, shows no energy dependence. Its mean value (eta0) approximately 0.006 cm(-1)) is also in close agreement to the measured data.

Validation of Monte Carlo generated phase-space descriptions of medical linear accelerators

The accuracy of Monte Carlo codes in dose calculation systems relies on the correctness of the input data. Monte Carlo calculations are performed to generate phase-space descriptions of the Varian 2100C accelerator at 6 and 18 MeV. Before these data can be reliably used as the input for dose calculations in patients, they must be properly validated. This validation consists of three different stages: validation of the coding of the geometry, validation of the user code for the Monte Carlo code, and validation of calculated results. Geometric validation is performed by isolating and testing treatment head components independently. The user code is checked by testing for energy conservation and the variance reduction schemes incorporated into the user code are checked by comparison of results calculated with and without their employment. Validation of the phase-space description is performed by calculation of depth dose curves and lateral profiles for dose deposition in phantom, with difference plots used to illustrate any discrepancies. Calculated and experimental in-phantom output is also determined. After complete validation, the calculated data can then be reliably used as the input for dose calculations.

An MCNP-based model of a linear accelerator x-ray beam

The Monte Carlo N-Particle radiation transport computer code (MCNP) has been employed on a personal computer to develop a simple model simulating the major components within the beam path of a linear accelerator radiation head, namely the electron target, primary conical collimator, beam flattening filter, wedge filter and the secondary collimators. The model was initially used to calculate the energy spectra and angular distributions of the x-ray beam for the Philips SL 75/5 linear accelerator, in a plane immediately beneath the flattening filter. These data were subsequently used as a 'source' of x-rays at the target position, to assess the emergent beam from the secondary collimators. The depth dose distributions and dose profiles at constant depth for various field sizes have been calculated for a nominal operating potential of 4 MV and found to be within acceptable limits. It is concluded that the technique may be used to calculate the energy spectra of any linear accelerator upon specification of the component dimensions, materials and nominal accelerating potential. It is anticipated that this work will serve as the basis of a quality control tool for linear accelerators and treatment planning systems.

The flatness of Siemens linear accelerator x-ray fields

The primary definer for Siemens MXE and MDX linear accelerators projects a circular opening with a radius of 25 cm at 100 cm from the target. Our measurements of photon beam profiles, however, indicate that the photon fluence drops to 95% of the central axis value at a radius of 18 cm. The flattening filter for these machines projects a flattened field size that is much smaller than the primary definer would allow. The clinical implications of this mismatch for large rectangular fields and for fields defined by asymmetric jaws are discussed and solutions are considered. A large field flattener was designed for our Siemens MXE 6 MV beam using Monte Carlo simulation of the treatment head and water phantom. The accuracy required of source and geometry details for dose distributions calculation is presented. The key parameters are the mean energy and focal spot size of the electron beam incident on the exit window, the material composition, and thickness profile of the exit window, target, flattener, and primary collimator, and the position of the primary collimator relative to the target. Profiles were more sensitive than central axis depth doses to simulation details. The beam energy and primary collimator position were selected to achieve good agreement between measured and calculated dose distributions. The flattener we designed with Monte Carlo was machined from brass and mounted on our MXE treatment unit. Measurements demonstrate that the large field flattener extends the useful radius of the field out to 22 cm, right into the penumbra cast by the primary collimator.