Using Monte Carlo simulations to commission photon beam output factors—a feasibility study

Technical note: Evaluation of the dose enhancement effect for a novel transmission-type x-ray tube using the Monte Carlo method

PURPOSE

Transmission-target x-ray tubes generate more x-rays than reflection thick-target x-ray tubes. A transmission x-ray tube combined with radiosensitizers has a better radiation enhancement effect. This study investigated the feasibility of using a transmission x-ray tube with radiosensitizers in clinical radiotherapy and its effect on radiation dose enhancement.

METHODS

This study used MCNP6.2 to simulate the model of a transmission x-ray tube and Co-60 beam.   The radiation enhancement effect of radiosensitizers was examined with iodine-127 (I-127), radioiodinated iododeoxyuridine (IUdR), and gold nanoparticles (GNPs).

RESULTS

The study results showed that the dose enhancement factor (DEF) of the transmission x-ray tube with GNPs was 10.27, which was higher than that of I-127 (6.46) and IUdR (3.08). The DEF of the Co-60 beam with GNPs, I-127, and IUdR was 1.23, 1.19, and 1.2, respectively. The Auger electron flux of the transmission x-ray tube with GNPs was 1.19E+05 particles/cm2 .

CONCLUSIONS

This study found that a transmission x-ray tube with appropriate radiosensitizers could produce a high rate of Auger electrons to fulfill the radiation enhancement effect, and this procedure has the potential to become a radiotherapy modality.

Beam modeling and beam model commissioning for Monte Carlo dose calculation-based radiation therapy treatment planning: Report of AAPM Task Group 157

Dose calculation plays an important role in the accuracy of radiotherapy treatment planning and beam delivery. The Monte Carlo (MC) method is capable of achieving the highest accuracy in radiotherapy dose calculation and has been implemented in many commercial systems for radiotherapy treatment planning. The objective of this task group was to assist clinical physicists with the potentially complex task of acceptance testing and commissioning MC-based treatment planning systems (TPS) for photon and electron beam dose calculations. This report provides an overview on the general approach of clinical implementation and testing of MC-based TPS with a specific focus on models of clinical photon and electron beams. Different types of beam models are described including those that utilize MC simulation of the treatment head and those that rely on analytical methods and measurements. The trade-off between accuracy and efficiency in the various source-modeling approaches is discussed together with guidelines for acceptance testing of MC-based TPS from the clinical standpoint. Specific recommendations are given on methods and practical procedures to commission clinical beam models for MC-based TPS.

The comparison between 6 MV Primus LINAC simulation output using EGSnrc and commissioning data

Introduction

Monte Carlo calculation method is considered to be the most accurate method for dose calculation in radiotherapy. The purpose of this research is comparison between 6 MV Primus LINAC simulation output with commissioning data using EGSnrc and build a Monte Carlo geometry of 6 MV Primus LINAC as realistically as possible. The BEAMnrc and DOSXYZnrc (EGSnrc package) Monte Carlo model of the LINAC head was used as a benchmark.

Methods

In the first part, the BEAMnrc was used for the designing of the LINAC treatment head. In the second part, dose calculation and for the design of 3D dose file were produced by DOSXYZnrc. The simulated PDD and beam profile obtained were compared with that calculated using commissioning data. Good agreement was found between calculated PDD (1·1%) and beam profile using Monte Carlo simulation and commissioning data. After validation, TPR20,10, TMR and Spvalues were calculated in five different field.

Results

Good agreement was found between calculated values by using Monte Carlo simulation and commissioning data. Average differences for five field sizes in this approach is about 0·83% for Sp. for TPR20,10differences for field sizes 10×10 cm2is 0·29% and for TMR in five field sizes, the average value is ~1·6%.

Conclusion

In conclusion, the BEAMnrc and DOSXYZnrc codes package have very good accuracy in calculating dose distribution for 6 MV photon beam and it can be considered as a promising method for patient dose calculations and also the Monte Carlo model of primus linear accelerator built in this study can be used as method to calculate the dose distribution for cancer patients.

A multi-centre analytical study of small field output factor calculations in radiotherapy

An audit methodology was developed and applied for output factor (OF) calculations in radiotherapy. The auditees were asked to calculate OFs for field sizes from 10 × 10 cm² to 2 × 2 cm². Sixty five beams were audited; missing reference OFs were interpolated. The calculated OFs were in 73% of cases higher than the reference data. The smaller the field size, the higher the overestimations which were observed in the higher fraction of cases. Treatment planning systems generally overestimated OFs for small fields. The reference dataset helped radiotherapy centres to identify discrepancies which were higher than typical.

Implementation of a double Gaussian source model for the BEAMnrc Monte Carlo code and its influence on small fields dose distributions

The shape of the radiation source of a linac has a direct impact on the delivered dose distributions, especially in the case of small radiation fields. Traditionally, a single Gaussian source model is used to describe the electron beam hitting the target, although different studies have shown that the shape of the electron source can be better described by a mixed distribution consisting of two Gaussian components. Therefore, this study presents the implementation of a double Gaussian source model into the BEAMnrc Monte Carlo code. The impact of the double Gaussian source model for a 6 MV beam is assessed through the comparison of different dosimetric parameters calculated using a single Gaussian source, previously commissioned, the new double Gaussian source model and measurements, performed with a diode detector in a water phantom. It was found that the new source can be easily implemented into the BEAMnrc code and that it improves the agreement between measurements and simulations for small radiation fields. The impact of the change in source shape becomes less important as the field size increases and for increasing distance of the collimators to the source, as expected. In particular, for radiation fields delivered using stereotactic collimators located at a distance of 59 cm from the source, it was found that the effect of the double Gaussian source on the calculated dose distributions is negligible, even for radiation fields smaller than 5 mm in diameter. Accurate determination of the shape of the radiation source allows us to improve the Monte Carlo modeling of the linac, especially for treatment modalities such as IMRT, were the radiation beams used could be very narrow, becoming more sensitive to the shape of the source.

PACS number(s): 87.53.Bn, 87.55.K, 87.56.B‐, 87.56.jf

Monte Carlo validation of the TrueBeam 10XFFF phase-space files for applications in lung SABR

PURPOSE

To establish the clinical acceptability of universal Monte Carlo phase-space data for the 10XFFF (flattening filter free) photon beam on the Varian TrueBeam Linac, including previously unreported data for small fields, output factors, and inhomogeneous media. The study was particularly aimed at confirming the suitability for use in simulations of lung stereotactic ablative radiotherapy treatment plans.

METHODS

Monte Carlo calculated percent depth doses (PDDs), transverse profiles, and output factors for the TrueBeam 10 MV FFF beam using generic phase-space data that have been released by the Varian MC research team were compared with in-house measurements and published data from multiple institutions (ten Linacs from eight different institutions). BEAMnrc was used to create field size specific phase-spaces located underneath the jaws. Doses were calculated with DOSXYZnrc in a water phantom for fields ranging from 1 × 1 to 40 × 40 cm(2). Particular attention was paid to small fields (down to 1 × 1 cm(2)) and dose per pulse effects on dosimeter response for high dose rate 10XFFF beams. Ion chamber measurements were corrected for changes in ion collection efficiency (P(ion)) with increasing dose per pulse. MC and ECLIPSE ANISOTROPIC ANALYTICAL ALGORITHM (AAA) calculated PDDs were compared to Gafchromic film measurement in inhomogeneous media (water, bone, lung).

RESULTS

Measured data from all machines agreed with Monte Carlo simulations within 1.0% and 1.5% for PDDs and in-field transverse profiles, respectively, for field sizes >1 × 1 cm(2) in a homogeneous water phantom. Agreements in the 80%-20% penumbra widths were better than 2 mm for all the fields that were compared. For all the field sizes considered, the agreement between their measured and calculated output factors was within 1.1%. Monte Carlo results for dose to water at water/bone, bone/lung, and lung/water interfaces as well as within lung agree with film measurements to within 2.8% for 10 × 10 and 3 × 3 cm(2) field sizes. This represents a significant improvement over the performance of the ECLIPSE AAA.

CONCLUSIONS

The 10XFFF phase-space data offered by the Varian Monte Carlo research team have been validated for clinical use using measured, interinstitutional beam data in water and with film dosimetry in inhomogeneous media.

Development and reproducibility evaluation of a Monte Carlo-based standard LINAC model for quality assurance of multi-institutional clinical trials

Technical developments in radiotherapy (RT) have created a need for systematic quality assurance (QA) to ensure that clinical institutions deliver prescribed radiation doses consistent with the requirements of clinical protocols. For QA, an ideal dose verification system should be independent of the treatment-planning system (TPS). This paper describes the development and reproducibility evaluation of a Monte Carlo (MC)-based standard LINAC model as a preliminary requirement for independent verification of dose distributions. The BEAMnrc MC code is used for characterization of the 6-, 10- and 15-MV photon beams for a wide range of field sizes. The modeling of the LINAC head components is based on the specifications provided by the manufacturer. MC dose distributions are tuned to match Varian Golden Beam Data (GBD). For reproducibility evaluation, calculated beam data is compared with beam data measured at individual institutions. For all energies and field sizes, the MC and GBD agreed to within 1.0% for percentage depth doses (PDDs), 1.5% for beam profiles and 1.2% for total scatter factors (Scps.). Reproducibility evaluation showed that the maximum average local differences were 1.3% and 2.5% for PDDs and beam profiles, respectively. MC and institutions' mean Scps agreed to within 2.0%. An MC-based standard LINAC model developed to independently verify dose distributions for QA of multi-institutional clinical trials and routine clinical practice has proven to be highly accurate and reproducible and can thus help ensure that prescribed doses delivered are consistent with the requirements of clinical protocols.

Monte Carlo modeling of HD120 multileaf collimator on Varian TrueBeam linear accelerator for verification of 6X and 6X FFF VMAT SABR treatment plans

A Monte Carlo (MC) validation of the vendor-supplied Varian TrueBeam 6 MV flattened (6X) phase-space file and the first implementation of the Siebers-Keall MC MLC model as applied to the HD120 MLC (for 6X flat and 6X flattening filter-free (6X FFF) beams) are described. The MC model is validated in the context of VMAT patient-specific quality assurance. The Monte Carlo commissioning process involves: 1) validating the calculated open-field percentage depth doses (PDDs), profiles, and output factors (OF), 2) adapting the Siebers-Keall MLC model to match the new HD120-MLC geometry and material composition, 3) determining the absolute dose conversion factor for the MC calculation, and 4) validating this entire linac/MLC in the context of dose calculation verification for clinical VMAT plans. MC PDDs for the 6X beams agree with the measured data to within 2.0% for field sizes ranging from 2 × 2 to 40 × 40 cm2. Measured and MC profiles show agreement in the 50% field width and the 80%-20% penumbra region to within 1.3 mm for all square field sizes. MC OFs for the 2 to 40 cm2 square fields agree with measurement to within 1.6%. Verification of VMAT SABR lung, liver, and vertebra plans demonstrate that measured and MC ion chamber doses agree within 0.6% for the 6X beam and within 2.0% for the 6X FFF beam. A 3D gamma factor analysis demonstrates that for the 6X beam, > 99% of voxels meet the pass criteria (3%/3 mm). For the 6X FFF beam, > 94% of voxels meet this criteria. The TrueBeam accelerator delivering 6X and 6X FFF beams with the HD120 MLC can be modeled in Monte Carlo to provide an independent 3D dose calculation for clinical VMAT plans. This quality assurance tool has been used clinically to verify over 140 6X and 16 6X FFF TrueBeam treatment plans.

A Monte Carlo approach to validation of FFF VMAT treatment plans for the TrueBeam linac

PURPOSE

To commission and benchmark a vendor-supplied (Varian Medical Systems) Monte Carlo phase-space data for the 6 MV flattening filter free (FFF) energy mode on a TrueBeam linear accelerator for the purpose of quality assurance of clinical volumetric modulated arc therapy (VMAT) treatment plans. A method for rendering the phase-space data compatible with BEAMnrc/DOSXYZnrc simulation software package is presented.

METHODS

Monte Carlo (MC) simulations were performed to benchmark the TrueBeam 6 MV FFF phase space data that have been released by the Varian MC Research team. The simulations to benchmark the phase space data were done in three steps. First, the original phase space which was created on a cylindrical surface was converted into a format that was compatible with BEAMnrc. Second, BEAMnrc was used to create field size specific phase spaces located underneath the jaws. Third, doses were calculated with DOSXYZnrc in a water phantom for fields ranging from 1 × 1 to 40 × 40 cm(2). Calculated percent depth doses (PDD), transverse profiles, and output factors were compared with measurements for all the fields simulated. After completing the benchmarking study, three stereotactic body radiotherapy (SBRT) VMAT plans created with the Eclipse treatment planning system (TPS) were calculated with Monte Carlo. Ion chamber and film measurements were also performed on these plans. 3D gamma analysis was used to compare Monte Carlo calculation with TPS calculations and with film measurement.

RESULTS

For the benchmarking study, MC calculated and measured values agreed within 1% and 1.5% for PDDs and in-field transverse profiles, respectively, for field sizes >1 × 1 cm(2). Agreements in the 80%-20% penumbra widths were better than 2 mm for all the fields that were compared. With the exception of the 1 × 1 cm(2) field, the agreement between measured and calculated output factors was within 1%. It is of note that excellent agreement in output factors for all field sizes including highly asymmetric fields was achieved without accounting for backscatter into the beam monitor chamber. For the SBRT VMAT plans, the agreement between Monte Carlo and ion chamber point dose measurements was within 1%. Excellent agreement between Monte Carlo, treatment planning system and Gafchromic film dose distribution was observed with over 99% of the points in the high dose volume passing the 3%, 3 mm gamma test.

CONCLUSIONS

The authors have presented a method for making the Varian IAEA compliant 6 MV FFF phase space file of the TrueBeam linac compatible with BEAMnrc/DOSXYZnrc. After benchmarking the modified phase space against measurement, they have demonstrated its potential for use in MC based quality assurance of complex delivery techniques.

Clinical commissioning and use of the Novalis Tx linear accelerator for SRS and SBRT

The purpose of this study was to perform comprehensive measurements and testing of a Novalis Tx linear accelerator, and to develop technical guidelines for com-missioning from the time of acceptance testing to the first clinical treatment. The Novalis Tx (NTX) linear accelerator is equipped with, among other features, a high-definition MLC (HD120 MLC) with 2.5 mm central leaves, a 6D robotic couch, an optical guidance positioning system, as well as X-ray-based image guidance tools to provide high accuracy radiation delivery for stereotactic radiosurgery and stereotactic body radiation therapy procedures. We have performed extensive tests for each of the components, and analyzed the clinical data collected in our clinic. We present technical guidelines in this report focusing on methods for: (1) efficient and accurate beam data collection for commissioning treatment planning systems, including small field output measurements conducted using a wide range of detectors; (2) commissioning tests for the HD120 MLC; (3) data collection for the baseline characteristics of the on-board imager (OBI) and ExacTrac X-ray (ETX) image guidance systems in conjunction with the 6D robotic couch; and (4) end-to-end testing of the entire clinical process. Established from our clinical experience thus far, recommendations are provided for accurate and efficient use of the OBI and ETX localization systems for intra- and extracranial treatment sites. Four results are presented. (1) Basic beam data measurements: Our measurements confirmed the necessity of using small detectors for small fields. Total scatter factors varied significantly (30% to approximately 62%) for small field measurements among detectors. Unshielded stereotactic field diode (SFD) overestimated dose by ~ 2% for large field sizes. Ion chambers with active diameters of 6 mm suffered from significant volume averaging. The sharpest profile penumbra was observed for the SFD because of its small active diameter (0.6 mm). (2) MLC commissioning: Winston Lutz test, light/radiation field congruence, and Picket Fence tests were performed and were within criteria established by the relevant task group reports. The measured mean MLC transmission and dynamic leaf gap of 6 MV SRS beam were 1.17% and 0.36 mm, respectively. (3) Baseline characteristics of OBI and ETX: The isocenter localization errors in the left/right, posterior/anterior, and superior/inferior directions were, respectively, -0.2 ± 0.2 mm, -0.8 ± 0.2 mm, and -0.8 ± 0.4 mm for ETX, and 0.5 ± 0.7 mm, 0.6 ± 0.5 mm, and 0.0 ± 0.5 mm for OBI cone-beam computed tomography. The registration angular discrepancy was 0.1 ± 0.2°, and the maximum robotic couch error was 0.2°. (4) End-to-end tests: The measured isocenter dose differences from the planned values were 0.8% and 0.4%, measured respectively by an ion chamber and film. The gamma pass rate, measured by EBT2 film, was 95% (3% DD and 1 mm DTA). Through a systematic series of quantitative commissioning experiments and end-to-end tests and our initial clinical experience, described in this report, we demonstrate that the NTX is a robust system, with the image guidance and MLC requirements to treat a wide variety of sites - in particular for highly accurate delivery of SRS and SBRT-based treatments.

Monte Carlo simulation for Neptun 10 PC medical linear accelerator and calculations of output factor for electron beam

AIM

Exact knowledge of dosimetric parameters is an essential pre-requisite of an effective treatment in radiotherapy. In order to fulfill this consideration, different techniques have been used, one of which is Monte Carlo simulation.

MATERIALS AND METHODS

This study used the MCNP-4C to simulate electron beams from Neptun 10 PC medical linear accelerator. Output factors for 6, 8 and 10 MeV electrons applied to eleven different conventional fields were both measured and calculated.

RESULTS

The measurements were carried out by a Wellhofler-Scanditronix dose scanning system. Our findings revealed that output factors acquired by MCNP-4C simulation and the corresponding values obtained by direct measurements are in a very good agreement.

CONCLUSION

In general, very good consistency of simulated and measured results is a good proof that the goal of this work has been accomplished.

Spectral method for the correction of the Cerenkov light effect in plastic scintillation detectors: a comparison study of calibration procedures and validation in Cerenkov light-dominated situations

PURPOSE

The purposes of this work were: (1) To determine if a spectral method can accurately correct the Cerenkov light effect in plastic scintillation detectors (PSDs) for situations where the Cerenkov light is dominant over the scintillation light and (2) to develop a procedural guideline for accurately determining the calibration factors of PSDs.

METHODS

The authors demonstrate, by using the equations of the spectral method, that the condition for accurately correcting the effect of Cerenkov light is that the ratio of the two calibration factors must be equal to the ratio of the Cerenkov light measured within the two different spectral regions used for analysis. Based on this proof, the authors propose two new procedures to determine the calibration factors of PSDs, which were designed to respect this condition. A PSD that consists of a cylindrical polystyrene scintillating fiber (1.6 mm3) coupled to a plastic optical fiber was calibrated by using these new procedures and the two reference procedures described in the literature. To validate the extracted calibration factors, relative dose profiles and output factors for a 6 MV photon beam from a medical linac were measured with the PSD and an ionization chamber. Emphasis was placed on situations where the Cerenkov light is dominant over the scintillation light and on situations dissimilar to the calibration conditions.

RESULTS

The authors found that the accuracy of the spectral method depends on the procedure used to determine the calibration factors of the PSD and on the attenuation properties of the optical fiber used. The results from the relative dose profile measurements showed that the spectral method can correct the Cerenkov light effect with an accuracy level of 1%. The results obtained also indicate that PSDs measure output factors that are lower than those measured with ionization chambers for square field sizes larger than 25 x 25 cm2, in general agreement with previously published Monte Carlo results.

CONCLUSIONS

The authors conclude that the spectral method can be used to accurately correct the Cerenkov light effect in PSDs. The authors confirmed the importance of maximizing the difference of Cerenkov light production between calibration measurements. The authors also found that the attenuation of the optical fiber, which is assumed to be constant in the original formulation of the spectral method, may cause a variation of the calibration factors in some experimental setups.

Monte Carlo simulation of the effect of miniphantom on in-air output ratio

PURPOSE

The aim of the study was to quantify the effect of miniphantoms on in-air output ratio measurements, i.e., to determine correction factors for in-air output ratio.

METHODS

Monte Carlo (MC) simulations were performed to simulate in-air output ratio measurements by using miniphantoms made of various materials (PMMA, graphite, copper, brass, and lead) and with different longitudinal thicknesses or depths (2-30 g/cm2) in photon beams of 6 and 15 MV, respectively, and with collimator settings ranging from 3 x 3 to 40 x 40 cm2. EGSnrc and BEAMnrc (2007) software packages were used. Photon energy spectra corresponding to the collimator settings were obtained from BEAMnrc code simulations on a linear accelerator and were used to quantify the components of in-air output ratio correction factors, i.e., attenuation, mass energy absorption, and phantom scatter correction factors. In-air output ratio correction factors as functions of miniphantom material, miniphantom longitudinal thickness, and collimator setting were calculated and compared to a previous experimental study.

RESULTS

The in-air output ratio correction factors increase with collimator opening and miniphantom longitudinal thickness for all the materials and for both energies. At small longitudinal thicknesses, the in-air output ratio correction factors for PMMA and graphite are close to 1. The maximum magnitudes of the in-air output ratio correction factors occur at the largest collimator setting (40 x 40 cm2) and the largest miniphantom longitudinal thickness (30 g/cm2): 1.008 +/- 0.001 for 6 MV and 1.012 +/- 0.001 for 15 MV, respectively. The MC simulations of the in-air output ratio correction factor confirm the previous experimental study.

CONCLUSIONS

The study has verified that a correction factor for in-air output ratio can be obtained as a product of attenuation correction factor, mass energy absorption correction factor, and phantom scatter correction factor. The correction factors obtained in the present study can be used in studies involving in-air output ratio measurements using miniphantoms.

Application of a Monte Carlo-based method for total scatter factors of small beams to new solid state micro-detectors

The scope of this work was to apply a method for estimation of total scatter factors of the smallest beams of the Cyberknife radiosurgery system to newly available solid‐state detectors: the PTW 60008 diode, the SunNuclear EdgeDetector™ diode, and the Thomson and Nielsen TN502RDM micromosfet. The method is based on a consistency check between Monte Carlo simulation of the detectors and experimental results, and was described in a recent publication. Corrected total scatter factors were in excellent agreement with the findings of the former study. The results showed that the diodes tend to overestimate the total scatter factor of small beams, probably due to excessive scatter from the material surrounding the active layer. The correction factor for diodes and for the micromosfet, however, was found to be independent of the electron beam width. This is a desirable characteristic because it allows standard correction factors to be used for treatment units of the same type, without the need of case‐by‐case Monte Carlo simulation.

PACS numbers: 87.55.kh; 87.55.ne; 87.56.Fc.

Total scatter factors of small beams: a multidetector and Monte Carlo study

The scope of this study was to estimate total scatter factors (S(c,p)) of the three smallest collimators of the Cyberknife radiosurgery system (5-10 mm in diameter), combining experimental measurements and Monte Carlo simulation. Two microchambers, a diode, and a diamond detector were used to collect experimental data. The treatment head and the detectors were simulated by means of a Monte Carlo code in order to calculate correction factors for the detectors and to estimate total scatter factors by means of a consistency check between measurement and simulation. Results for the three collimators were: S(c,p) (5 mm) = 0.677 +/- 0.004, S(c,p) (7.5 mm) = 0.820 +/- 0.008, S(c,p) (10 mm) = 0.871 +/- 0.008, all relative to the 60 mm collimator at 80 cm source-to-detector distance. The method also allows the full width at half maximum of the electron beam to be estimated; estimations made with different collimators and different detectors were in excellent agreement and gave a value of 2.1 mm. Correction factors to be applied to the detectors for the measurement of S(c,p) were consistent with a prevalence of volume effect for the microchambers and the diamond and a prevalence of scattering from high-Z material for the diode detector. The proposed method is more sensitive to small variations of the electron beam diameter with respect to the conventional method used to commission Monte Carlo codes, i.e., by comparison with measured percentage depth doses (PDD) and beam profiles. This is especially important for small fields (less than 10 mm diameter), for which measurements of PDD and profiles are strongly affected by the type of detector used. Moreover, this method should allow S(c,p) of Cyberknife systems different from the unit under investigation to be estimated without the need for further Monte Carlo calculation, provided that one of the microchambers or the diode detector of the type used in this study are employed. The results for the diamond are applicable only to the specific detector that was investigated due to excessive variability in manufacturing.

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 calculations for absolute dosimetry to determine machine outputs for proton therapy fields

The prescribed dose in radiation therapy has to be converted into machine monitor units for patient treatment. This is done routinely for each spread-out Bragg peak (SOBP) field either by calibration measurements, by using analytical algorithms or by relying on empirical data. At the Northeast Proton Therapy Center, a monitor unit corresponds to a fixed amount of charge collected in a segmented transmission ionization chamber inside the treatment head. The goal of this work was to use a detailed Monte Carlo model of the treatment head to calculate the dose delivered to the patient as a function of ionization chamber reading, i.e. to yield absolute dose in patients in terms of machine monitor units. The results show excellent agreement with measurements. For 50 SOBP fields considered in this study, the mean absolute difference between the experimental and the calculated value is 1.5%, where approximately 50% of the fields agree within 1%. This is within the uncertainties of the data. The Monte Carlo method has advantages over analytical algorithms because it takes into account scattered and secondary radiation, does not rely on empirical parameters, and provides a tool to study the influence of parts of the treatment head on the ionization chamber reading. Compared to experimental methods the Monte Carlo method has the advantage of being able to verify the dose in the patient geometry.

Commissioning stereotactic radiosurgery beams using both experimental and theoretical methods

The purpose of this investigation is to study the feasibility of using an alternative method to commission stereotactic radiosurgery beams shaped by micro multi-leaf collimators by using Monte Carlo simulations to obtain beam characteristics of small photon beams, such as incident beam particle fluence and energy distributions, scatter ratios, depth-dose curves and dose profiles where measurements are impossible or difficult. Ionization chambers and diode detectors with different sensitive volumes were used in the measurements in a water phantom and the Monte Carlo codes BEAMnrc/DOSXYZnrc were used in the simulation. The Monte Carlo calculated data were benchmarked against measured data for photon beams with energies of 6 MV and 10 MV produced from a Varian Trilogy accelerator. The measured scatter ratios and cross-beam dose profiles for very small fields are shown to be not only dependent on the size of the sensitive volume of the detector used but also on the type of detectors. It is known that the response of some detectors changes at small field sizes. Excellent agreement was seen between scatter ratios measured with a small ion chamber and those calculated from Monte Carlo simulations. The values of scatter ratios, for field sizes from 6 x 6 mm2 to 98 x 98 mm2, range from 0.67 to 1.0 and from 0.59 to 1.0 for 6 and 10 MV, respectively. The Monte Carlo calculations predicted that the incident beam particle fluence is strongly affected by the X-Y-jaw openings, especially for small fields due to the finite size of the radiation source. Our measurement confirmed this prediction. This study demonstrates that Monte Carlo calculations not only provide accurate dose distributions for small fields where measurements are difficult but also provide additional beam characteristics that cannot be obtained from experimental methods. Detailed beam characteristics such as incident photon fluence distribution, energy spectra, including composition of primary and scattered photons, can be independently used in dose calculation models and to improve the accuracy of measurements with detectors with an energy-dependent response. Furthermore, when there are discrepancies between results measured with different detectors, the Monte Carlo calculated values can indicate the most correct result. The data set presented in this study can be used as a reference in commissioning stereotactic radiosurgery beams shaped by a BrainLAB m3 on a Varian 2100EX or 600C accelerator.

Absolute dose calculations for Monte Carlo simulations of radiotherapy beams

Monte Carlo (MC) simulations have traditionally been used for single field relative comparisons with experimental data or commercial treatment planning systems (TPS). However, clinical treatment plans commonly involve more than one field. Since the contribution of each field must be accurately quantified, multiple field MC simulations are only possible by employing absolute dosimetry. Therefore, we have developed a rigorous calibration method that allows the incorporation of monitor units (MU) in MC simulations. This absolute dosimetry formalism can be easily implemented by any BEAMnrc/DOSXYZnrc user, and applies to any configuration of open and blocked fields, including intensity-modulated radiation therapy (IMRT) plans. Our approach involves the relationship between the dose scored in the monitor ionization chamber of a radiotherapy linear accelerator (linac), the number of initial particles incident on the target, and the field size. We found that for a 10 x 10 cm2 field of a 6 MV photon beam, 1 MU corresponds, in our model, to 8.129 x 10(13) +/- 1.0% electrons incident on the target and a total dose of 20.87 cGy +/- 1.0% in the monitor chambers of the virtual linac. We present an extensive experimental verification of our MC results for open and intensity-modulated fields, including a dynamic 7-field IMRT plan simulated on the CT data sets of a cylindrical phantom and of a Rando anthropomorphic phantom, which were validated by measurements using ionization chambers and thermoluminescent dosimeters (TLD). Our simulation results are in excellent agreement with experiment, with percentage differences of less than 2%, in general, demonstrating the accuracy of our Monte Carlo absolute dose calculations.

Reference photon dosimetry data and reference phase space data for the 6MV photon beam from Varian Clinac 2100 series linear accelerators

The current study presents the reference photon dosimetry data (RPDD) and reference phase space data (RPSD) for the 6 MV photon beam from Varian 2100 series linear accelerators. The RPDD provide the basic photon dosimetry data, typically collected during the initial commissioning of a new linear accelerator, including output factors, depth dose data, and beam profile data in air and in water. The RPSD provide the full phase space information, such as position, direction, and energy for each particle generated inside the head of any particular linear accelerator in question. The dosimetric characteristics if the 6 MV photon beam from the majority of the aforementioned accelerators, which are unaltered from the manufacturer's original specifications, can be fully described with these two data sets within a clinically acceptable uncertainty (approximately +/-2 %). The current study also presents a detailed procedure to establish the RPDD and RPSD using measured data and Monte Carlo calculations. The RPDD were constructed by compiling our own measured data and the average data based on the analysis of more than 50 sets of measured data from the Radiological Physics Center (RPC) and 10 sets of clinical dosimetry data obtained from 10 different institutions participating in the RPC's quality assurance monitoring program. All the measured data from the RPC and the RPC-monitored institutions were found to be within a statistically tight range (i.e., 1sigma approximately 1% or less) for each dosimetric quantity. The manufacturer's standard data, except for in-air off-axis factors that are available only from the current study, were compared with the RPDD, showing that the manufacturer's standard data could also be used as the RPDD for the photon beam studied in this study. The RPSD were obtained from Monte Carlo calculations using the BEAMnrc/ DOSXYZnrc code system with 6.2 MeV (a spread of 3% full width at half maximum) and 1.0 mm full width at half maximum as the values of the energy and radial spread of a Gaussian electron pencil beam incident on the target, respectively. The RPSD were capable of generating Monte Carlo data that agreed with the RPDD within the acceptance criteria adopted in the current study (e.g., 1% or 1 mm for depth dose). A complete set of the RPDD and RPSD from the current study is available from the RPC website (http://rpc.mdanderson.org) or via mass storage media such as DVD or CD-ROM upon request.

An investigation of accelerator head scatter and output factor in air

Our purpose in this study was to investigate whether the Monte Carlo simulation can accurately predict output factors in air. Secondary goals were to study the head scatter components and investigate the collimator exchange effect. The Monte Carlo code, BEAMnrc, was used in the study. Photon beams of 6 and 18 MV were from a Varian Clinac 2100EX accelerator and the measurements were performed using an ionization chamber in a mini-phantom. The Monte Carlo calculated in air output factors was within 1% of measured values. The simulation provided information of the origin and the magnitude of the collimator exchange effect. It was shown that the collimator backscatter to the beam monitor chamber played a significant role in the beam output factors. However the magnitude of the scattered dose contributions from the collimator at the isocenter is negligible. The maximum scattered dose contribution from the collimators was about 0.15% and 0.4% of the total dose at the isocenter for a 6 and 18 MV beam, respectively. The scattered dose contributions from the flattening filter at the isocenter were about 0.9-3% and 0.2-6% of the total dose for field sizes of 4x4 cm2-40x40 cm2 for the 6 and 18 MV beam, respectively. The study suggests that measurements of head scatter factors be done at large depth well beyond the depth of electron contamination. The insight information may have some implications for developing generalized empirical models to calculate the head scatter.

A new approach in dose measurement and error analysis for narrow photon beams (beamlets) shaped by different multileaf collimators using a small detector

Dose measurement for narrow stereotactic beams and intensity modulation radiotherapy beamlets is difficult and error-prone due to the lack of lateral electron equilibrium. A small detector position error and finite sensitive volume as well as the nonfocus collimation could result in considerable (> 10%) measurement errors. A new method is introduced here to measure the dose and error components so that the accuracy and precision of the dose measurement can be improved. Based on superposition principle, we can create exactly the small field of interest by subtraction of a reference open field (O-field) and two strip fields (S-fields) from the sum of four quadrant fields (Q-fields). The position effect on the dose measurement is determined by the standard deviation of the four Q-field readings. The collimator leaf-edge effect (LEE) is quantified by the difference between the readings of the two S-fields using a detector that has very small sensitive volume. The detector-volume effect can be analytically estimated from the integrals of the dose distributions of the two S-fields over the detector volume. Using a pinpoint ion chamber (PTW N31006) and a stereotactic silicon diode detector (Scanditronix, DEB050), we have measured scatter factors (SF) and tissue-maximum ratios for 6-MV x-ray fields with sizes of 3 x 3 and 6 x 6 mm2 shaped by a BrainLAB micromultileaf collimator (microMLC) (M3), 4.4 x 4.4 and 8.8 x 8.8 mm2 shaped by a 3DLine double-focus MLC, and 5 x 5 and 10 x 10 mm2 by a Varian Millennium MLC. Our experimental results demonstrate that the large errors are often caused by a small setup error or measuring point displacement from the central ray of the beam. The LEE is almost independent of the depth but closely related to the field size and the type of MLC. The volume effect becomes significant when the detector diameter is comparable to the half size of the small fields. Application of the new method using different detectors had achieved less than 8.3% total experiment error for all of the small fields of interest except for the SF of the 3 x 3 mm2 field from the pinpoint ion chamber that has 15% volume effect. Importantly, the new method using a solid water phantom is clinically convenient and highly reproducible.