Efficient photon beam dose calculations using DOSXYZnrc with BEAMnrc

Geometrical splitting technique to improve the computational efficiency in Monte Carlo calculations for proton therapy

PURPOSE

To present the implementation and validation of a geometrical based variance reduction technique for the calculation of phase space data for proton therapy dose calculation.

METHODS

The treatment heads at the Francis H Burr Proton Therapy Center were modeled with a new Monte Carlo tool (TOPAS based on Geant4). For variance reduction purposes, two particle-splitting planes were implemented. First, the particles were split upstream of the second scatterer or at the second ionization chamber. Then, particles reaching another plane immediately upstream of the field specific aperture were split again. In each case, particles were split by a factor of 8. At the second ionization chamber and at the latter plane, the cylindrical symmetry of the proton beam was exploited to position the split particles at randomly spaced locations rotated around the beam axis. Phase space data in IAEA format were recorded at the treatment head exit and the computational efficiency was calculated. Depth-dose curves and beam profiles were analyzed. Dose distributions were compared for a voxelized water phantom for different treatment fields for both the reference and optimized simulations. In addition, dose in two patients was simulated with and without particle splitting to compare the efficiency and accuracy of the technique.

RESULTS

A normalized computational efficiency gain of a factor of 10-20.3 was reached for phase space calculations for the different treatment head options simulated. Depth-dose curves and beam profiles were in reasonable agreement with the simulation done without splitting: within 1% for depth-dose with an average difference of (0.2 ± 0.4)%, 1 standard deviation, and a 0.3% statistical uncertainty of the simulations in the high dose region; 1.6% for planar fluence with an average difference of (0.4 ± 0.5)% and a statistical uncertainty of 0.3% in the high fluence region. The percentage differences between dose distributions in water for simulations done with and without particle splitting were within the accepted clinical tolerance of 2%, with a 0.4% statistical uncertainty. For the two patient geometries considered, head and prostate, the efficiency gain was 20.9 and 14.7, respectively, with the percentages of voxels with gamma indices lower than unity 98.9% and 99.7%, respectively, using 2% and 2 mm criteria.

CONCLUSIONS

The authors have implemented an efficient variance reduction technique with significant speed improvements for proton Monte Carlo simulations. The method can be transferred to other codes and other treatment heads.

An investigation into the use of MMCTP to tune accelerator source parameters and testing its clinical application

This paper presents an alternative method to tune Monte Carlo electron beam parameters to match measured data using a minimal set of variables in order to reduce the model setup time prior to clinical implementation of the model. Monte Carlo calculations provide the possibility of a powerful treatment planning verification technique. The nonstandardized and nonautomated process of tuning the required accelerator model is one of the reasons for delays in the clinical implementation of Monte Carlo techniques. This work aims to establish and verify an alternative tuning method that can be carried out in a minimal amount of time, allowing it to be easily implemented in a clinical setting by personnel with minimal experience with Monte Carlo methods. This tuned model can then be incorporated into the MMCTP system to allow the system to be used as a second dose calculation check for IMRT plans. The technique proposed was used to establish the primary electron beam parameters for accelerator models for the Varian Clinac 2100 6 MV photon beam using the BEAMnrc Monte Carlo system. The method is intended to provide a clear, direct, and efficient process for tuning an accelerator model using readily available clinical quality assurance data. The tuning provides a refined model, which agrees with measured dose profile curves within 1.5% outside the penumbra or 3 mm in the penumbra, for square fields with sides of 3 cm up to 30 cm. These models can then be employed as the basis for Monte Carlo recalculations of dose distributions, using the MMCTP system, for clinical treatment plans, providing an invaluable assessment tool. This was tested on six IMRT plans and compared to the measurements performed for the pretreatment QA process. These Monte Carlo values for the average dose to the chamber volume agreed with measurements to within 0.6%.

PACS number: 87.55.km

Is wax equivalent to tissue in electron conformal therapy planning? A Monte Carlo study of material approximation introduced dose difference

With CT-based Monte Carlo (MC) dose calculations, material composition is often assigned based on the standard Hounsfield unit ranges. This is known as the density threshold method. In bolus electron conformal therapy (BolusECT), the bolus material, machineable wax, would be assigned as soft tissue and the electron density is assumed equivalent to soft tissue based on its Hounsfield unit. This study investigates the dose errors introduced by this material assignment. BEAMnrc was used to simulate electron beams from a Trilogy accelerator. SPRRZnrc was used to calculate stopping power ratios (SPR) of tissue to wax, SPR (tissue) (wax), and tissue to water, SPR(tissue) (water), for 6, 9, 12, 15, and 18 MeV electron beams, of which 12 and 15MeV beams are the most commonly used energies in BolusECT. DOSXYZnrc was applied in dose distribution calculations in a tissue phantom with either flat wax slabs of various thicknesses or a wedge-shaped bolus on top. Dose distribution for two clinical cases, a chest wall and a head and neck, were compared with the bolus material treated as wax or tissue. The SPR(tissue) (wax) values for 12 and 15MeV beams are between 0.935 and 0.945, while the SPR(tissue) (water) values are between 0.990 and 0.991. For a 12 MeV beam, the dose in tissue immediately under the bolus is overestimated by 2.5% for a 3 cm bolus thickness if the wax bolus is treated as tissue. For 15 MeV beams, the error is 1.4%. However, in both clinical cases the differences in the PTV DVH is negligible. Due to stopping power differences, dose differences of up to 2.5% are observed in MC simulations if the bolus material is misassigned as tissue in BolusECT dose calculations. However, for boluses thinner than 2 cm that are more likely encountered in practice, the error is within clinical tolerance.

Monte Carlo modeling of a Novalis Tx Varian 6 MV with HD-120 multileaf collimator

A Monte Carlo model of the Novalis Tx linear accelerator equipped with high-definition multileaf collimator (HD-120 HD-MLC) was commissioned using ionization chamber measurements in water. All measurements in water were performed using a liquid filled ionization chamber. Film measurements were made using EDR2 film in solid water. Open rectangular fields defined by the jaws or the HD-MLC were used for comparison against measurements. Furthermore, inter- and intraleaf leakage calculated by the Monte Carlo model was compared against film measurements. The statistical uncertainty of the Monte Carlo calculations was less than 1% for all simulations. Results for all regular field sizes show an excellent agreement with commissioning data (percent depth-dose curves and profiles), well within 1% of difference in the relative dose and 1 mm distance to agreement. The computed leakage through HD-MLCs shows good agreement with film measurements. The Monte Carlo model developed in this study accurately represents the new Novalis Tx Varian linac with HD-MLC and can be used for reliable patient dose calculations.

Internal scatter, the unavoidable major component of the peripheral dose in photon-beam radiotherapy

In clinical photon beams, the dose outside the geometrical field limits is produced by photons originating from (i) head leakage, (ii) scattering at the beam collimators and the flattening filter (head scatter) and (iii) scattering from the directly irradiated region of the patient or phantom (internal scatter). While the first two components can be modified, e.g. by reinforcement of shielding components or by re-modeling the filter system, internal scatter remains an unavoidable contributor to the peripheral dose. Its relative magnitude compared to the other components, its numerical variation with beam energy, field size and off-axis distance as well as its spectral distribution are evaluated in this study. We applied a detailed Monte Carlo (MC) model of our 6/15 MV Siemens Primus linear accelerator beam head, provided with ideal head leakage shielding conditions (multi-leaf collimator without gaps) to assess the head scatter contribution. Experimental values obtained under real shielding conditions were used to evaluate the head leakage contribution. It was found that the MC-computed internal scatter doses agree with the results of our previous measurements, that internal scatter is the major contributor to the peripheral dose in the near periphery while head leakage prevails in the far periphery, and that the lateral decline of the internal scatter dose can be represented by the sum of two exponentials, with an asymptotic tenth value of 18 to 19 cm. Internal scatter peripheral doses from various elementary beams are additive, so that their sum increases approximately in proportion with field size. The ratio between normalized internal scatter doses at 6 and 15 MV is approximately 2:1. The energy fluence spectra of the internal scatter component at all points of interest outside the field have peaks near 500 keV. The fact that the energy-shifted internal scatter constitutes the major contributor to the dose in the near periphery has a general bearing for dosimetry, i.e. for energy-dependent detector responses and dose conversion factors, for the relative biological effectiveness and for second primary malignancy risk estimates in the peripheral region.

The development and experimental evaluation of a simple analytical model for the TPR in the build-up region of megavoltage photon beams

PURPOSE

The purpose of this work was to develop a simple analytical model for the tissue phantom ratio (TPR) in the build-up region of megavoltage photon beams and to experimentally evaluate the model under a variety of clinically relevant field configurations.

METHODS

Considering electron contamination and primary photons as the main components of the absorbed dose in the build-up region, an analytic expression for the TPR was derived. The electron contamination component was addressed with a biexponential function; the primary photon component was treated as nonlocal energy transport, i.e., assuming the energy deposited by secondary electrons can be described by a biexponential mode similar to that of contaminating electrons. The model contains five independent constants, which were fitted experimentally. The accuracy of the model was evaluated by comparing its results with in-phantom measurements taken on square, rectangular, irregular, and wedged fields, for 6 and 15 MV photon beams on a GE-Saturne 41 accelerator. More specifically, the accuracy of the model was quantified using the gamma index with 2% dose and 2 mm spatial tolerances as described by Low et al. [Med. Phys. 25, 656-661 (1998)].

RESULTS

For square cerrobende blocked fields, the maximum recorded gamma indices were 0.42 and 0.54 for 6 and 15 MV beams, respectively. For "I" shaped fields, the corresponding maxima were 0.64 and 0.52, respectively, while for "cross" shaped fields they were 0.42 and 0.76. For rectangular 10 × 30 cm fields, the corresponding maxima were 0.32 and 0.42, and for 7 × 20 cm fields, they were 0.70 and 0.35, respectively. For square 10 × 10 cm and 15 × 15 cm fields with an acrylic tray, the maxima were 0.57 and 0.45 for 6 MV and 0.32 and 0.77 for 15 MV beams, respectively. For a 10 × 10 cm 60° wedged field, the maxima were 0.53 and 0.33 for 6 and 15 MV beams, respectively. In all examined cases of irregular, rectangular, square (with and without tray), and wedged fields, the gamma index was less than unity. Thus, the model correctly predicted TPR in all cases, using the defined criteria.

CONCLUSIONS

A simple analytical model for the TPR in the build-up region was developed and evaluated experimentally. The model's predicted TPR values were compared with physical measurements for irregular, square (with and without tray), rectangular, and wedged fields, for 6 and 15 MV photon beams. In every case examined, the results of the model agreed with the experimental measurements based on specific quantitative agreement criteria. The model appears useful for predicting the TPR in the build-up region of megavoltage beams for different types of fields, in different configurations.

A Monte Carlo pencil beam scanning model for proton treatment plan simulation using GATE/GEANT4

This work proposes a generic method for modeling scanned ion beam delivery systems, without simulation of the treatment nozzle and based exclusively on beam data library (BDL) measurements required for treatment planning systems (TPS). To this aim, new tools dedicated to treatment plan simulation were implemented in the Gate Monte Carlo platform. The method was applied to a dedicated nozzle from IBA for proton pencil beam scanning delivery. Optical and energy parameters of the system were modeled using a set of proton depth-dose profiles and spot sizes measured at 27 therapeutic energies. For further validation of the beam model, specific 2D and 3D plans were produced and then measured with appropriate dosimetric tools. Dose contributions from secondary particles produced by nuclear interactions were also investigated using field size factor experiments. Pristine Bragg peaks were reproduced with 0.7 mm range and 0.2 mm spot size accuracy. A 32 cm range spread-out Bragg peak with 10 cm modulation was reproduced with 0.8 mm range accuracy and a maximum point-to-point dose difference of less than 2%. A 2D test pattern consisting of a combination of homogeneous and high-gradient dose regions passed a 2%/2 mm gamma index comparison for 97% of the points. In conclusion, the generic modeling method proposed for scanned ion beam delivery systems was applicable to an IBA proton therapy system. The key advantage of the method is that it only requires BDL measurements of the system. The validation tests performed so far demonstrated that the beam model achieves clinical performance, paving the way for further studies toward TPS benchmarking. The method involves new sources that are available in the new Gate release V6.1 and could be further applied to other particle therapy systems delivering protons or other types of ions like carbon.

Experimental measurements and Monte Carlo simulations of dose perturbation around a nonradioactive brachytherapy seed due to 6- and 18-MV photons

PURPOSE

Radioactive seeds used in permanent prostate brachytherapy are composed of high-Z metals and may exceed 100 in a patient. If supplemental external beam treatment is administered afterward, the seeds may cause substantial dose perturbation, which is being investigated in this article.

METHODS AND MATERIALS

Film measurements using 6-MV beam were primarily carried out using Kodak XV2 film layered above and below a nonradioactive iodine-125 ((125)I) seed. Monte Carlo simulations were carried out using DOSXYZnrc. Other experimental comparisons looked at changing beam energy, depth, and field size, including two opposing fields' pair. Effect of multiple seeds spatially spaced 0.5cm vertically was also studied.

RESULTS

For a single (125)I seed, on XV film, there is a localized dose enhancement of 6.3% upstream and -10.9% downstream. With two opposing fields, a cold spot around the seed of ∼3% was noticed. Increasing beam energy and field size decreased the magnitude of this effect, whereas the effect was found to increase with the increasing Z of material. DOSXYZnrc and EBT-2 film verified maximum dose enhancement of +15% upstream and -20% downstream of the (125)I seed surface.

CONCLUSIONS

In general, the dose perturbation because of the seeds was spatially limited to ∼2mm upstream and ∼5mm downstream to the incident beam. Similar to other heterogeneities, the seeds perturbation depends on incident beam energy, field size, and its Z. With multiple seeds spatially apart and multiple radiation fields routinely used in external beam radiotherapy, the cumulative effect may not result in clinically significant dose perturbation.

Validation and automation of the DYNJAWS component module of the BEAMnrc Monte Carlo code

The purpose of this work is to validate and automate the use of DYNJAWS; a new component module (CM) in the BEAMnrc Monte Carlo (MC) user code. The DYNJAWS CM simulates dynamic wedges and can be used in three modes; dynamic, step-and-shoot and static. The step-and-shoot and dynamic modes require an additional input file defining the positions of the jaw that constitutes the dynamic wedge, at regular intervals during its motion. A method for automating the generation of the input file is presented which will allow for the more efficient use of the DYNJAWS CM. Wedged profiles have been measured and simulated for 6 and 10 MV photons at three field sizes (5 cm × 5 cm, 10 cm × 10 cm and 20 cm × 20 cm), four wedge angles (15°, 30°, 45° and 60°), at d (max) and at 10 cm depth. Results of this study show agreement between the measured and the MC profiles to within 3% of absolute dose or 3 mm distance to agreement for all wedge angles at both energies and depths. The gamma analysis suggests that dynamic mode is more accurate than the step-and-shoot mode. The DYNJAWS CM is an important addition to the BEAMnrc code and will enable the MC verification of patient treatments involving dynamic wedges.

A Monte Carlo investigation of contaminant electrons due to a novel in vivo transmission detector

A novel transmission detector (IBA Dosimetry, Germany) developed as an IMRT quality assurance tool, intended for in vivo patient dose measurements, is studied here. The goal of this investigation is to use Monte Carlo techniques to characterize treatment beam parameters in the presence of the detector and to compare to those of a plastic block tray (a frequently used clinical device). Particular attention is paid to the impact of the detector on electron contamination model parameters of two commercial dose calculation algorithms. The linac head together with the COMPASS transmission detector (TRD) was modeled using BEAMnrc code. To understand the effect of the TRD on treatment beams, the contaminant electron fluence, energy spectra, and angular distributions at different SSDs were analyzed for open and non-open (i.e. TRD and block tray) fields. Contaminant electrons in the BEAMnrc simulations were separated according to where they were created. Calculation of surface dose and the evaluation of contributions from contaminant electrons were performed using the DOSXYZnrc user code. The effect of the TRD on contaminant electrons model parameters in Eclipse AAA and Pinnacle(3) dose calculation algorithms was investigated. Comparisons of the fluence of contaminant electrons produced in the non-open fields versus open field show that electrons created in the non-open fields increase at shorter SSD, but most of the electrons at shorter SSD are of low energy with large angular spread. These electrons are out-scattered or absorbed in air and contribute less to surface dose at larger SSD. Calculated surface doses with the block tray are higher than those with the TRD. Contribution of contaminant electrons to dose in the buildup region increases with increasing field size. The additional contribution of electrons to surface dose increases with field size for TRD and block tray. The introduction of the TRD results in a 12% and 15% increase in the Gaussian widths used in the contaminant electron source model of the Eclipse AAA dose algorithm. The off-axis coefficient in the Pinnacle(3) dose calculation algorithm decreases in the presence of TRD compared to without the device. The electron model parameters were modified to reflect the increase in electron contamination with the TRD, a necessary step for accurate beam modeling when using the device.

Relation between lineal energy distribution and relative biological effectiveness for photon beams according to the microdosimetric kinetic model

Our cell survival data showed the obvious dependence of RBE on photon energy: The RBE value for 200 kV X-rays was approximately 10% greater than those for mega-voltage photon beams. In radiation therapy using mega-voltage photon beams, the photon energy distribution outside the field is different with that in the radiation field because of a large number of low energy scattering photons. Hence, the RBE values outside the field become greater. To evaluate the increase in RBE, the method of deriving the RBE using the Microdosimetric Kinetic model (MK model) was proposed in this study. The MK model has two kinds of the parameters, tissue-specific parameters and the dose-mean lineal energy derived from the lineal energy distributions measured with a Tissue-Equivalent Proportional Counter (TEPC). The lineal energy distributions with the same geometries of the cell irradiations for 200 kV X-rays, (60)Co γ-rays, and 6 MV X-rays were obtained with the TEPC and Monte Carlo code GEANT4. The measured lineal energy distribution for 200 kV X-rays was quite different from those for mega-voltage photon beams. The dose-mean lineal energy of 200 kV X-rays showed the greatest value, 4.51 keV/µm, comparing with 2.34 and 2.36 keV/µm for (60)Co γ-rays and 6 MV X-rays, respectively. By using the results of the TEPC and cell irradiations, the tissue-specific parameters in the MK model were determined. As a result, the RBE of the photon beams (y(D): 2~5 keV/µm) in arbitrary conditions can be derived by the measurements only or the calculations only of the dose-mean lineal energy.

Mathematical modeling of the radiation dose received from photons passing over and through shielding walls in a PET/CT suite

Given that the financial cost of shielding PET/CT suites can be substantial, it has become increasingly important to be able to accurately assess the thickness of shielding required for barriers and whether it is necessary to extend such shielding all the way to the ceiling. The overall shielding requirement for a PET/CT installation must take into account both 511 keV gamma ray emissions from PET scans and lower energy x-ray scatter from CT scans. This paper deals with the overall impact of emissions from both modalities. Radiation exposure from both scatter over shielding barriers as well as transmission through these barriers is taken into account. A series of simulations of the dose received by a person positioned behind a shielding barrier in a typical PET/CT scanning suite were carried out using both Monte Carlo and analytical models. The transmission through lead barriers was found to be very dependent on the geometry of the radiation source and the resulting energy spectrum of the emitted radiation. The transmission from a patient source was found to be around half of that from a small vial and also half of that reported previously using parallel beams of mono-energetic radiation. For PET emissions, the dose from scatter over the barrier at waist height is relatively small but may have to be taken into account if the design dose limit is low. Shielding from floor to ceiling is probably not warranted in most instances for PET gamma emissions; in PET/CT installations, however, a thinner layer of shielding may need to extend to the ceiling of the imaging room to limit x-ray scatter over the wall from the CT unit.

Four-dimensional dosimetry validation and study in lung radiotherapy using deformable image registration and Monte Carlo techniques

Thoracic cancer treatment presents dosimetric difficulties due to respiratory motion and lung inhomogeneity. Monte Carlo and deformable image registration techniques have been proposed to be used in four-dimensional (4D) dose calculations to overcome the difficulties. This study validates the 4D Monte Carlo dosimetry with measurement, compares 4D dosimetry of different tumor sizes and tumor motion ranges, and demonstrates differences of dose-volume histograms (DVH) with the number of respiratory phases that are included in 4D dosimetry. BEAMnrc was used in dose calculations while an optical flow algorithm was used in deformable image registration and dose mapping. Calculated and measured doses of a moving phantom agreed within 3% at the center of the moving gross tumor volumes (GTV). 4D CT image sets of lung cancer cases were used in the analysis of 4D dosimetry. For a small tumor (12.5 cm³) with motion range of 1.5 cm, reduced tumor volume coverage was observed in the 4D dose with a beam margin of 1 cm. For large tumors and tumors with small motion range (around 1 cm), the 4D dosimetry did not differ appreciably from the static plans. The dose-volume histogram (DVH) analysis shows that the inclusion of only extreme respiratory phases in 4D dosimetry is a reasonable approximation of all-phase inclusion for lung cancer cases similar to the ones studied, which reduces the calculation in 4D dosimetry.

Fast, accurate photon beam accelerator modeling using BEAMnrc: a systematic investigation of efficiency enhancing methods and cross-section data

In this work, an investigation of efficiency enhancing methods and cross-section data in the BEAMnrc Monte Carlo (MC) code system is presented. Additionally, BEAMnrc was compared with VMC++, another special-purpose MC code system that has recently been enhanced for the simulation of the entire treatment head. BEAMnrc and VMC++ were used to simulate a 6 MV photon beam from a Siemens Primus linear accelerator (linac) and phase space (PHSP) files were generated at 100 cm source-to-surface distance for the 10 x 10 and 40 x 40 cm2 field sizes. The BEAMnrc parameters/techniques under investigation were grouped by (i) photon and bremsstrahlung cross sections, (ii) approximate efficiency improving techniques (AEITs), (iii) variance reduction techniques (VRTs), and (iv) a VRT (bremsstrahlung photon splitting) in combination with an AEIT (charged particle range rejection). The BEAMnrc PHSP file obtained without the efficiency enhancing techniques under study or, when not possible, with their default values (e.g., EXACT algorithm for the boundary crossing algorithm) and with the default cross-section data (PEGS4 and Bethe-Heitler) was used as the "base line" for accuracy verification of the PHSP files generated from the different groups described previously. Subsequently, a selection of the PHSP files was used as input for DOSXYZnrc-based water phantom dose calculations, which were verified against measurements. The performance of the different VRTs and AEITs available in BEAMnrc and of VMC++ was specified by the relative efficiency, i.e., by the efficiency of the MC simulation relative to that of the BEAMnrc base-line calculation. The highest relative efficiencies were approximately 935 (approximately 111 min on a single 2.6 GHz processor) and approximately 200 (approximately 45 min on a single processor) for the 10 x 10 field size with 50 million histories and 40 x 40 cm2 field size with 100 million histories, respectively, using the VRT directional bremsstrahlung splitting (DBS) with no electron splitting. When DBS was used with electron splitting and combined with augmented charged particle range rejection, a technique recently introduced in BEAMnrc, relative efficiencies were approximately 420 (approximately 253 min on a single processor) and approximately 175 (approximately 58 min on a single processor) for the 10 x 10 and 40 x 40 cm2 field sizes, respectively. Calculations of the Siemens Primus treatment head with VMC++ produced relative efficiencies of approximately 1400 (approximately 6 min on a single processor) and approximately 60 (approximately 4 min on a single processor) for the 10 x 10 and 40 x 40 cm2 field sizes, respectively. BEAMnrc PHSP calculations with DBS alone or DBS in combination with charged particle range rejection were more efficient than the other efficiency enhancing techniques used. Using VMC++, accurate simulations of the entire linac treatment head were performed within minutes on a single processor. Noteworthy differences (+/- 1%-3%) in the mean energy, planar fluence, and angular and spectral distributions were observed with the NIST bremsstrahlung cross sections compared with those of Bethe-Heitler (BEAMnrc default bremsstrahlung cross section). However, MC calculated dose distributions in water phantoms (using combinations of VRTs/AEITs and cross-section data) agreed within 2% of measurements. Furthermore, MC calculated dose distributions in a simulated water/air/water phantom, using NIST cross sections, were within 2% agreement with the BEAMnrc Bethe-Heitler default case.

Kilovoltage beam Monte Carlo dose calculations in submillimeter voxels for small animal radiotherapy

PURPOSE

Small animal conformal radiotherapy (RT) is essential for preclinical cancer research studies and therefore various microRT systems have been recently designed. The aim of this paper is to efficiently calculate the dose delivered using our microRT system based on a microCT scanner with the Monte Carlo (MC) method and to compare the MC calculations to film measurements.

METHODS

Doses from 2-30 mm diameter 120 kVp photon beams deposited in a solid water phantom with 0.2 x 0.2 x 0.2 mm3 voxels are calculated using the latest versions of the EGSnrc codes BEAMNRC and DOSXYZNRC. Two dose calculation approaches are studied: a two-step approach using phase-space files and direct dose calculation with BEAMNRC simulation sources. Due to the small beam size and submillimeter voxel size resulting in long calculation times, variance reduction techniques are studied. The optimum bremsstrahlung splitting number (NBRSPL in BEAMNRC) and the optimum DOSXYZNRC photon splitting (Nsplit) number are examined for both calculation approaches and various beam sizes. The dose calculation efficiencies and the required number of histories to achieve 1% statistical uncertainty--with no particle recycling--are evaluated for 2-30 mm beams. As a final step, film dose measurements are compared to MC calculated dose distributions.

RESULTS

The optimum NBRSPL is approximately 1 x 10(6) for both dose calculation approaches. For the dose calculations with phase-space files, Nsplit varies only slightly for 2-30 mm beams and is established to be 300. Nsplit for the DOSXYZNRC calculation with the BEAMNRC source ranges from 300 for the 30 mm beam to 4000 for the 2 mm beam. The calculation time significantly increases for small beam sizes when the BEAMNRC simulation source is used compared to the simulations with phase-space files. For the 2 and 30 mm beams, the dose calculations with phase-space files are more efficient than the dose calculations with BEAMNRC sources by factors of 54 and 1.6, respectively. The dose calculation efficiencies converge for beams with diameters larger than 30 mm.

CONCLUSIONS

A very good agreement of MC calculated dose distributions to film measurements is found. The mean difference of percentage depth dose curves between calculated and measured data for 2, 5, 10, and 20 mm beams is 1.8%.

Monte Carlo-based adaptive EPID dose kernel accounting for different field size responses of imagers

The aim of this study is to present an efficient method to generate imager-specific Monte Carlo (MC)-based dose kernels for amorphous silicon-based electronic portal image device dose prediction and determine the effective backscattering thicknesses for such imagers. EPID field size-dependent responses were measured for five matched Varian accelerators from three institutions with 6 MV beams at the source to detector distance (SDD) of 105 cm. For two imagers, measurements were made with and without the imager mounted on the robotic supporting arm. Monoenergetic energy deposition kernels with 0-2.5 cm of water backscattering thicknesses were simultaneously computed by MC to a high precision. For each imager, the backscattering thickness required to match measured field size responses was determined. The monoenergetic kernel method was validated by comparing measured and predicted field size responses at 150 cm SDD, 10 x 10 cm2 multileaf collimator (MLC) sliding window fields created with 5, 10, 20, and 50 mm gaps, and a head-and-neck (H&N) intensity modulated radiation therapy (IMRT) patient field. Field size responses for the five different imagers deviated by up to 1.3%. When imagers were removed from the robotic arms, response deviations were reduced to 0.2%. All imager field size responses were captured by using between 1.0 and 1.6 cm backscatter. The predicted field size responses by the imager-specific kernels matched measurements for all involved imagers with the maximal deviation of 0.34%. The maximal deviation between the predicted and measured field size responses at 150 cm SDD is 0.39%. The maximal deviation between the predicted and measured MLC sliding window fields is 0.39%. For the patient field, gamma analysis yielded that 99.0% of the pixels have gamma < 1 by the 2%, 2 mm criteria with a 3% dose threshold. Tunable imager-specific kernels can be generated rapidly and accurately in a single MC simulation. The resultant kernels are imager position independent and are able to predict fields with varied incident energy spectra and a H&N IMRT patient field. The proposed adaptive EPID dose kernel method provides the necessary infrastructure to build reliable and accurate portal dosimetry systems.

A method to improve accuracy and precision of water surface identification for photon depth dose measurements

The objective of this study is to present a method to reduce the setup error inherent in clinical depth dose measurements and, in doing so, to improve entrance dosimetry measurement reliability. Ionization chamber (IC) depth dose measurements are acquired with the depth scan extended into the air above the water surface. An inflection region is obtained in each resulting percent depth ionization (PDI) curve that can be matched against other measurements or to an inflection region obtained from an analogous Monte Carlo (MC) simulation. Measurements are made with various field sizes for the 6 and 18 MV photon beams, with and without a Pb foil in the beam, to determine the sensitivity of the dose inflection region to the beam conditions. The offset between reference and test data set inflection regions is quantified using two separate methods. When comparing sets of measured data, maxima in the second derivative of ionization are compared. When comparing measured data to MC simulation, the offset that minimizes the sum of squared differences between the reference and test curves in the ionization inflection region is found. These methods can be used to quantify the offset between an initial setup (test) position and the true surface (reference) position. The ionization inflection location is found to be insensitive to changes in field size, electron contamination, and beam energy. Data from a single reference condition should be sufficient to identify the surface location. The method of determining IC offsets is general and should be applicable to any IC and other radiation sources. The measurement method could reduce the time and effort required in the initial IC setup at a water surface as setup errors can be corrected offline. Given a reliable set of reference data to compare with, this method could increase the ability of quality assurance (QA) measurements to detect discrepancies in beam output as opposed to discrepancies in IC localization. Application of the measurement method standardizes the procedure for localizing cylindrical ICs at a water surface and thereby improves the reliability of measurements taken with these devices at all depths.

Monte Carlo simulation of the dose to nuclear medicine staff wearing protective garments

The literature contains both endorsements of, and advice against, the use of protective apparel in nuclear medicine procedures. The main issues usually centre around: Whether the shielding which can be provided by a protective garment light enough to wear (0 to 0.6 mm lead equivalent at the gamma energies commonly encountered in nuclear medicine) is enough to warrant its use; and (more recently); Whether the dose enhancement behind the protective garment from electron scatter in lead is sufficient to be of concern. In this work, the Monte Carlo code EGSnrc was used to investigate the effectiveness of lead of thicknesses of 0 to 0.6 mm, in shielding staff from photons of energies of 140 and 511 keV. Furthermore, dose escalation behind the lead was investigated. Reasonable dose reductions are obtained at 140 keV with protective garments of 0.5 mm lead equivalence. This perhaps warrants their use, in certain circumstances. At 511 keV, the reduction in dose is less than 10%, and their use is probably not justified (given the weight that has to be carried) from an ALARA point of view. It should be noted here that protective garments designed for X-ray shielding will generally not have the same lead equivalence at the gamma energies used in nuclear medicine. It should also be noted that protective garments which do not contain lead do not always attenuate as much as their stated lead equivalence claims. Dose escalation does occur, but the depth of penetration of the scattered electrons beyond the exit side of the lead shielding is such that it is highly unlikely that a significant dose would be delivered to viable tissue in wearers of protective garments.

Electron beam therapy at extended source-to-surface distance: a Monte Carlo investigation

Electron‐beam therapy is used to treat superficial tumors at a standard 100 cm source‐to‐surface distance (SSD). However, certain clinical situations require the use of an extended SSD. In the present study, Monte Carlo methods were used to investigate clinical electron beams, at standard and non‐standard SSDs, from a Siemens Oncor Avant Garde (Siemens Healthcare, Erlangen, Germany) linear accelerator (LINAC). The LINAC treatment head was modeled in BEAMnrc for electron fields 5 cm in diameter and 10×10 cm, 15×15 cm, and 20×20 cm; for 6 MeV, 9 MeV, and 12 MeV; and for 100 cm, 110 cm, and 120 cm SSD. The DOSXYZnrc code was used to calculate extended SSD factors and dose contributions from various parts of the treatment head.

The main effects of extended SSD on water phantom dose distributions were verified by Monte Carlo methods. Monte Carlo–calculated and measured extended SSD factors showed an average difference of ±1.8%. For the field 5 cm in diameter, the relative output at extended SSD declined more rapidly than it did for the larger fields. An investigation of output contributions showed this decline was mainly a result of a rapid loss of scatter dose reaching the dmax point from the lower scrapers of the electron applicator. The field 5 cm in diameter showed a reduction in dose contributions; the larger fields generally showed an increased contribution from the scrapers with increase in SSD. Angular distributions of applicator‐scattered electrons have shown a large number of acute‐angle electron tracks contributing to the output for larger field sizes, explaining the shallow output reduction.

PACS numbers: 87.53.Wz, 87.53.Vb, 87.53.Hv

Treatment planning for a small animal using Monte Carlo simulation

The development of a small animal model for radiotherapy research requires a complete setup of customized imaging equipment, irradiators, and planning software that matches the sizes of the subjects. The purpose of this study is to develop and demonstrate the use of a flexible in-house research environment for treatment planning on small animals. The software package, called DOSCTP, provides a user-friendly platform for DICOM computed tomography-based Monte Carlo dose calculation using the EGSnrcMP-based DOSXYZnrc code. Validation of the treatment planning was performed by comparing the dose distributions for simple photon beam geometries calculated through the Pinnacle3 treatment planning system and measurements. A treatment plan for a mouse based on a CT image set by a 360-deg photon arc is demonstrated. It is shown that it is possible to create 3D conformal treatment plans for small animals with consideration of inhomogeneities using small photon beam field sizes in the diameter range of 0.5-5 cm, with conformal dose covering the target volume while sparing the surrounding critical tissue. It is also found that Monte Carlo simulation is suitable to carry out treatment planning dose calculation for small animal anatomy with voxel size about one order of magnitude smaller than that of the human.

VMC++ versus BEAMnrc: a comparison of simulated linear accelerator heads for photon beams

BEAMnrc, a code for simulating medical linear accelerators based on EGSnrc, has been bench-marked and used extensively in the scientific literature and is therefore often considered to be the gold standard for Monte Carlo simulations for radiotherapy applications. However, its long computation times make it too slow for the clinical routine and often even for research purposes without a large investment in computing resources. VMC++ is a much faster code thanks to the intensive use of variance reduction techniques and a much faster implementation of the condensed history technique for charged particle transport. A research version of this code is also capable of simulating the full head of linear accelerators operated in photon mode (excluding multileaf collimators, hard and dynamic wedges). In this work, a validation of the full head simulation at 6 and 18 MV is performed, simulating with VMC++ and BEAMnrc the addition of one head component at a time and comparing the resulting phase space files. For the comparison, photon and electron fluence, photon energy fluence, mean energy, and photon spectra are considered. The largest absolute differences are found in the energy fluences. For all the simulations of the different head components, a very good agreement (differences in energy fluences between VMC++ and BEAMnrc <1%) is obtained. Only a particular case at 6 MV shows a somewhat larger energy fluence difference of 1.4%. Dosimetrically, these phase space differences imply an agreement between both codes at the <1% level, making VMC++ head module suitable for full head simulations with considerable gain in efficiency and without loss of accuracy.

A "HOWFARLESS" option to increase efficiency of homogeneous phantom calculations with DOSXYZnrc

This paper describes a "HOWFARLESS" transport option, which has been added to DOSXYZnrc to increase the efficiency of beam commissioning calculations in homogeneous phantoms. The algorithm speeds up charged particle transport by only considering the distance to the extreme outer boundaries of the phantom, thus eliminating the need to stop at voxel boundaries. Dose is deposited by approximating the total curved charged particle steps by two straight-line steps joined at a hinge point. Good agreement with normal simulations is achieved at all beam energies and for all practical maximum step lengths with a 1:1 mixture of approximations based on the initial position/ direction of the particle and on its final position/direction. Use of the "HOWFARLESS" option in phantom calculations for 6 and 18 MV photon beams (10 x 10 cm2 and 40 x 40 cm2 fields) from BEAMnrc-simulated accelerators increases the efficiency at the optimum photon splitting number by a factor of 2.9-5.4 when the exact EGSnrc boundary crossing algorithm (BCA) is used and by 51%-89% when the faster PRESTA-I BCA is employed. The efficiency gain due to the "HOWFARLESS" transport option increases with increasing beam energy and decreases with increasing field/dose voxel size. Efficiency improvement is greater when the efficiency of the particle source itself is not a factor, and in such cases the "HOWFARLESS" option improves the DOSXYZnrc efficiency by up to a factor of 13.1 (exact BCA) or 3.5 (PRESTA-I BCA) for photon beams, and up to a factor of 17.2 (exact BCA) or 5.2 (PRESTA-I BCA) for electron beams.

Dose calculation validation of VMC++ for photon beams

The radiation therapy specific Voxel Monte Carlo (VMC+ +) dose calculation algorithm achieves a dramatic improvement in MC dose calculation efficiency for radiation therapy treatment planning dose evaluation compared with other MC algorithms. This work aims to validate VMC+ + for radiation therapy photon beam planning. VMC++ was validated with respect to the well-benchmarked EGS-based DOSXYZnrc by comparing depth dose and lateral profiles for field sizes ranging from 1 X 1 to 40 x 40 cm(2) for 6 and 18 MV beams in a homogeneous water phantom and in a simulated bone-lung-bone phantom. Patient treatment plan dose distributions were compared for five prostate plans and five head-and-neck (H/N) plans, all using intensity-modulated radiotherapy beams. For all tests, the same incident particles were used in both codes to isolate differences due to modeling of the radiation source. Voxel-by-voxel observed differences were analyzed to distinguish between systematic and purely statistical differences. Dose-volume-histogram-derived dose indices were compared for the patient plans. For the homogeneous water phantom and the bone-lung-bone phantom, the depth dose curve predicted by VMC+ + agreed with that predicted by DOSXYZnrc within expected statistical uncertainty in all voxels except the surface voxel of the water phantom, where VMC+ + predicted a lower dose. When the electron cutoff parameter was decreased for both codes, the surface voxel agreed within expected statistical uncertainty. For prostate plans, the most severe difference between the codes resulted in 55% of the voxels showing a systematic difference of 0.32% of maximum dose. For H/N plans, the largest difference observed resulted in 2% of the voxels showing a systematic difference of 0.98% of maximum dose. For the prostate plans, the most severe difference in the planning target volume D95 was 0.4%, the rectum D35 was 0.2%, the rectum DI7 was 0.2%, the bladder D50 was 0.3% and the bladder D25 was 0.3%. For the H/N plans, the most severe difference in the gross tumor volume D98 was 0.4%, the clinical target volume D90 was 0.2%, the nodes D90 was 0.2%, the parotids D95 was 0.8%, and the cord D2 was 0.8%. All of these differences are clinically insignificant. VMC++ showed an average efficiency gain over DOSXYZnrc of at least an order of magnitude without introducing significant systematic bias. VMC + + can be used for photon beam MC patient dose computations without a clinically significant loss in accuracy.

Monte Carlo study of Si diode response in electron beams

Silicon semiconductor diodes measure almost the same depth-dose distributions in both photon and electron beams as those measured by ion chambers. A recent study in ion chamber dosimetry has suggested that the wall correction factor for a parallel-plate ion chamber in electron beams changes with depth by as much as 6%. To investigate diode detector response with respect to depth, a silicon diode model is constructed and the water/silicon dose ratio at various depths in electron beams is calculated using EGSnrc. The results indicate that, for this particular diode model, the diode response per unit water dose (or water/diode dose ratio) in both 6 and 18 MeV electron beams is flat within 2% versus depth, from near the phantom surface to the depth of R50 (with calculation uncertainty <0.3%). This suggests that there must be some other correction factors for ion chambers that counter-balance the large wall correction factor at depth in electron beams. In addition, the beam quality and field-size dependence of the diode model are also calculated. The results show that the water/diode dose ratio remains constant within 2% over the electron energy range from 6 to 18 MeV. The water/diode dose ratio does not depend on field size as long as the incident electron beam is broad and the electron energy is high. However, for a very small beam size (1 X 1 cm(2)) and low electron energy (6 MeV), the water/diode dose ratio may decrease by more than 2% compared to that of a broad beam.