Head scatter off-axis for megavoltage x rays

Detailed Monte Carlo analysis of the secondary photons coming out of the therapeutic X-ray beam of linear accelerator

External photon beam radiotherapy is often used in tumor treatment. The photons are generated from the target which had stricken by the primary electron beam (incident particles). The photon beam contains the primary photons coming directly from the target and secondary photons coming from the photon interactions with head component materials (scattered photons). Altogether is thereafter used in radiotherapy treatment. This Monte Carlo study aims to investigate and evaluate the secondary radiations (photons) in terms of fluence, energy fluence, spectral distribution, mean energy and angular spread distribution.

The secondary photons, which contributed in radiotherapy treatment, are examined and evaluated in number (fluence) and energy. At the phantom surface, the secondary photons originated in the whole linac head are mainly coming from the primary collimator. In 0.45% of secondary photons coming from the whole linac head, the primary collimator contributes by 86% and they are more energetic. However, the flattening filter and the secondary collimator contribute together by less than 14% and their photons are less energetic and then can deteriorate the beam dosimetry quality. To improve the radiotherapy treatment quality, the number of photons of low energy should be as low as possible in the clinical beam. Our work can be a basic investigation to use in the improvement of linac head configuration and specially the beam modifiers.

Part I: Out-of-field dose mapping for 6X and 6X-flattening-filter-free beams on the TrueBeam for extended distances

PURPOSE

With increasing cancer treatment success rates, many patients go on to live long, productive lives following recovery. Therefore, minimizing potential side effects due to dose outside the treated field is becoming a significant consideration in radiation therapy. With many potential treatment configurations available, it is important to quantify how out-of-field dose varies with common variables such as distance from isocenter, couch angle, jaw size, and flattening-filter setting. The accurate quantification of out-of-field dose at extended distances could also benefit researchers and detector developers. While data exist for out-of-field dose from older linear accelerator (Linac) models, the phenomenon has not been described for the latest generation of machines, such as the Varian TrueBeam. The purpose of this study was to comprehensively quantify out-of-field dose for the Varian TrueBeam Linac low energy photons in a wide range of positions and treatment geometries.

METHOD AND MATERIALS

Out-of-field doses were measured using two phantom setups: (a) A large volume ion chamber with a buildup sleeve to quantify head leakage and collimator scatter background dose; and (b) A farmer ion chamber in solid water to incorporate phantom scatter in addition to collimator scatter, and head leakage background dose. In both cases, the ion chamber was positioned with its length along the slowly varying transverse direction (perpendicular to the radial from isocenter). Doses were measured for four symmetric jaw settings (2 × 2 cm² , 4 × 4 cm² , 10 × 10 cm² , and 20 × 20 cm² ) for a range of distances from the isocenter (0-100 cm). The angular dependence of the out-of-field dose was measured using four different angles: 0°, 45°, 90°, and 135° with respect to the in-plane direction. All measurements were performed for both 6X and 6X-flattening-filter-free (FFF) beams.

RESULTS

The lowest out-of-field doses were observed at 60 cm away from isocenter in both in-plane and cross-plane directions for fields smaller than 10 × 10 cm² . Out-of-field dose decreased with decreasing jaw size (a factor of 4.7 for 6X-FFF and a factor of 3.1 for 6X going from 20 × 20 cm² to 2 × 2 cm² at 60 cm from isocenter in the in-plane direction). The 6X-FFF beam produced out-of-field doses as low as 64% of the 6X beam.

CONCLUSION

This study presents a comprehensive description of 6X and 6X-FFF out-of-field doses on a Varian TrueBeam Linac including measurements at a range of positions, angles, and jaw settings and with and without phantom scatter.

Monitor unit calculations for external photon and electron beams: Report of the AAPM Therapy Physics Committee Task Group No. 71

A protocol is presented for the calculation of monitor units (MU) for photon and electron beams, delivered with and without beam modifiers, for constant source-surface distance (SSD) and source-axis distance (SAD) setups. This protocol was written by Task Group 71 of the Therapy Physics Committee of the American Association of Physicists in Medicine (AAPM) and has been formally approved by the AAPM for clinical use. The protocol defines the nomenclature for the dosimetric quantities used in these calculations, along with instructions for their determination and measurement. Calculations are made using the dose per MU under normalization conditions, D'0, that is determined for each user's photon and electron beams. For electron beams, the depth of normalization is taken to be the depth of maximum dose along the central axis for the same field incident on a water phantom at the same SSD, where D'0 = 1 cGy/MU. For photon beams, this task group recommends that a normalization depth of 10 cm be selected, where an energy-dependent D'0 ≤ 1 cGy/MU is required. This recommendation differs from the more common approach of a normalization depth of dm, with D'0 = 1 cGy/MU, although both systems are acceptable within the current protocol. For photon beams, the formalism includes the use of blocked fields, physical or dynamic wedges, and (static) multileaf collimation. No formalism is provided for intensity modulated radiation therapy calculations, although some general considerations and a review of current calculation techniques are included. For electron beams, the formalism provides for calculations at the standard and extended SSDs using either an effective SSD or an air-gap correction factor. Example tables and problems are included to illustrate the basic concepts within the presented formalism.

Multisource modeling of flattening filter free (FFF) beam and the optimization of model parameters

PURPOSE

With the introduction of flattening filter free (FFF) linear accelerators to radiation oncology, new analytical source models for a FFF beam applicable to current treatment planning systems is needed. In this work, a multisource model for the FFF beam and the optimization of involved model parameters were designed.

METHODS

The model is based on a previous three source model proposed by Yang et al. ["A three-source model for the calculation of head scatter factors," Med. Phys. 29, 2024-2033 (2002)]. An off axis ratio (OAR) of photon fluence was introduced to the primary source term to generate cone shaped profiles. The parameters of the source model were determined from measured head scatter factors using a line search optimization technique. The OAR of the photon fluence was determined from a measured dose profile of a 40 x 40 cm2 field size with the same optimization technique, but a new method to acquire gradient terms for OARs was developed to enhance the speed of the optimization process. The improved model was validated with measured dose profiles from 3 x 3 to 40 x 40 cm2 field sizes at 6 and 10 MV from a TrueBeam STx linear accelerator. Furthermore, planar dose distributions for clinically used radiation fields were also calculated and compared to measurements using a 2D array detector using the gamma index method.

RESULTS

All dose values for the calculated profiles agreed with the measured dose profiles within 0.5% at 6 and 10 MV beams, except for some low dose regions for larger field sizes. A slight overestimation was seen in the lower penumbra region near the field edge for the large field sizes by 1%-4%. The planar dose calculations showed comparable passing rates (> 98%) when the criterion of the gamma index method was selected to be 3%/3 mm.

CONCLUSIONS

The developed source model showed good agreements between measured and calculated dose distributions. The model is easily applicable to any other linear accelerator using FFF beams as the required data include only the measured PDD, dose profiles, and output factors for various field sizes, which are easily acquired during conventional beam commissioning process.

Validation of the Eclipse AAA algorithm at extended SSD

The accuracy of dose calculations at extended SSD is of significant importance in the dosimetric planning of total body irradiation (TBI). In a first step toward the implementation of electronic, multi-leaf collimator compensation for dose inhomogeneities and surface contour in TBI, we have evaluated the ability of the Eclipse AAA to accurately predict dose distributions in water at extended SSD. For this purpose, we use the Eclipse AAA algorithm, commissioned with machine-specific beam data for a 6 MV photon beam, at standard SSD (100 cm). The model was then used for dose distribution calculations at extended SSD (179.5 cm). Two sets of measurements were acquired for a 6 MV beam (from a Varian linear accelerator) in a water tank at extended SSD: i) open beam for 5 x 5, 10 x 10, 20 x 20 and 40 x 40 cm2 field sizes (defined at 179.5 cm SSD), and ii) identical field sizes but with a 1.3 cm thick acrylic spoiler placed 10 cm above the water surface. Dose profiles were acquired at 5 cm, 10 cm and 20 cm depths. Dose distributions for the two setups were calculated using the AAA algorithm in Eclipse. Confidence limits for comparisons between measured and calculated absolute depth dose curves and normalized dose profiles were determined as suggested by Venselaar et al. The confidence limits were within 2% and 2 mm for both setups. Extended SSD calculations were also performed using Eclipse AAA, commissioned with Varian Golden beam data at standard SSD. No significant difference between the custom commissioned and Golden Eclipse AAA was observed. In conclusion, Eclipse AAA commissioned at standard SSD can be used to accurately predict dose distributions in water at extended SSD for 6 MV open beams.

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.

Determination of depth and field size dependence of multileaf collimator transmission in intensity-modulated radiation therapy beams

Intensity‐modulated radiation therapy (IMRT) plans for the treatment of large and complex volumes may contain a relatively large contribution from multileaf collimator (MLC) transmission. In such cases, comprehensive characterization of direct and scatter MLC transmission is important. We designed a set of tests (open beam, closed static MLC, and dynamic MLC gap) to determine dosimetric MLC properties as a function of field size and depth at the central axis.

We developed a generalized model of MLC transmission to account for direct MLC transmission, MLC scatter, beam hardening, and leaf‐end transmission (dosimetric gap). The model is consistent with the beam model used in IMRT optimization. We tested the model for extreme asymmetric fields relevant for large targets and for split IMRT fields. We applied our MLC scatter estimation formula to clinically relevant cases and showed that MLC scatter is contributing an undesired background dose. This contribution is relatively large, especially in low‐dose regions. (For instance, a uniform extra dose may dramatically increase normal‐lung toxicity in thorax treatment.) For complex IMRT of large‐volume targets, we found direct MLC transmission dose to be as high as 30%, and MLC scatter, up to 10% within the target volume for the selected cases. We identified that the dose discrepancies between the IMRT planning system [Eclipse (Varian Medical Systems, Palo Alto, CA)] and ionization chamber measurements (inside and outside of the field) are attributable to an inadequate model of MLC transmission in the planning system (constant‐value model).

In the present study, we measured MLC transmission properties for Varian 6EX (6 MV) and 21EXs (6 and 10 MV) linear accelerators; however, the experimental method and theoretical model are more general.

PACS number: 87.53.‐j

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.

A novel approach to accurate portal dosimetry using CCD-camera based EPIDs

A new method for portal dosimetry using CCD camera-based electronic portal imaging devices (CEPIDs) is demonstrated. Unlike previous approaches, it is not based on a priori assumptions concerning CEPID cross-talk characteristics. In this method, the nonsymmetrical and position-dependent cross-talk is determined by directly imaging a set of cross-talk kernels generated by small fields ("pencil beams") exploiting the high signal-to-noise ratio of a cooled CCD camera. Signal calibration is achieved by imaging two reference fields. Next, portal dose images (PDIs) can be derived from electronic portal dose images (EPIs), in a fast forward-calculating iterative deconvolution. To test the accuracy of these EPI-based PDIs, a comparison is made to PDIs obtained by scanning diode measurements. The method proved accurate to within 0.2+/-0.7% (1 SD), for on-axis symmetrical and asymmetrical fields with different field widths and homogeneous phantom thicknesses, off-axis Alderson thorax fields and a strongly modulated IMRT field. Hence, the proposed method allows for fast, accurate portal dosimetry. In addition, it is demonstrated that the CEPID cross-talk signal is not only induced by optical photon reflection and scatter within the CEPID structure, but also by high-energy back-scattered radiation from CEPID elements (mirror and housing) towards the fluorescent screen.

Influence of ion chamber response on in-air profile measurements in megavoltage photon beams

This article presents an investigation of the influence of the ion chamber response, including buildup caps, on the measurement of in-air off-axis ratio (OAR) profiles in megavoltage photon beams using Monte Carlo simulations with the EGSnrc system. Two new techniques for the calculation of OAR profiles are presented. Results of the Monte Carlo simulations are compared to measurements performed in 6, 10 and 25 MV photon beams produced by an Elekta Precise linac and shown to agree within the experimental and simulation uncertainties. Comparisons with calculated in-air kerma profiles demonstrate that using a plastic mini phantom gives more accurate air-kerma measurements than using high-Z material buildup caps and that the variation of chamber response with distance from the central axis must be taken into account.

Output ratio in air for MLC shaped irregular fields

For accurate monitor unit calculation, it is important to calculate the output ratio in air, Sc, for an irregular field shaped by MLC. We have developed an algorithm to calculate Sc based on an empirical model [Med. Phys. 28, 925-937 (2001)] by projecting each leaf position to the isocenter plane. Thus it does not require the exact knowledge of the head geometry. Comparisons were made for three different types of MLC: those with MLC replacing the inner collimator jaws; those with MLC replacing the outer collimator jaws; and those with MLC as a tertiary attachment. When the MLC leaf positions are substantially different from the secondary collimators (or the rectangular field encompassing the irregular field), one observes an up to 5% difference in the value of head-scatter correction factor, HCF, defined as the ratio of output ratio in air between the MLC shaped irregular field and that of the rectangular field encompassing the irregular field. No collimator exchange effect was observed for rectangular fields shaped by MLC (e.g., 5x30 and 30x5 cm2 diagonal) when the secondary collimators are fixed, unlike that for the rectangular fields shaped by the inner and outer collimator jaws, where it can be 1-2%. For the same MLC shaped irregular field, the value of Sc increases from the Elekta, to the Siemens, to the Varian accelerators, with an up to 4% difference. The calculation agrees with measurement to within 1.2% for points both on and off the central-axis. The fitting parameters used in the algorithm are derived from measurements for square field sizes on the central-axis.