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

A novel analytical method for computing dose from kilovoltage beams used in Image-Guided radiation therapy

PURPOSE

A modified convolution/superposition algorithm is proposed to compute dose from the kilovoltage beams used in IGRT. The algorithm uses material-specific energy deposition kernels instead of water-energy deposition kernels.

METHODS

Monte Carlo simulation was used to model the Elekta XVI unit and determine dose deposition characteristics of its kilovoltage beams. The dosimetric results were compared with ion chamber measurements. The dose from the kilovoltage beams was then computed using convolution/superposition along with material-specific energy deposition kernels and compared with Monte Carlo and measurements. The material-specific energy deposition kernels were previously generated using Monte Carlo.

RESULTS

The obtained gamma indices (using 2%/2mm criteria for 95% of points) were lower than 1 in almost all instances which indicates good agreement between simulated and measured depth doses and profiles. The comparisons of the algorithm with measurements in a homogeneous solid water slab phantom, and that with Monte Carlo in a head and neck CT dataset produced acceptable results. The calculated point doses were within 4.2% of measurements in the homogeneous phantom. Gamma analysis of the calculated vs. Monte Carlo simulations in the head and neck phantom resulted in 94% of points passing with a 2%/2mm criteria.

CONCLUSIONS

The proposed method offers sufficient accuracy in kilovoltage beams dose calculations and has the potential to supplement the conventional megavoltage convolution/superposition algorithms for dose calculations in low energy range.

Increasing efficiency of BEAMnrc-simulated Co-60 beams using directional source biasing

PURPOSE

This study describes the implementation of a directional source biasing (DSB) scheme for efficiently simulating Cobalt-60 treatment heads using the BEAMnrc Monte Carlo code. Previous simulation of Co-60 beams with BEAMnrc was impractical because of the time required to track photons not directed into the treatment field and to simulate secondary charged particles.

METHODS

In DSB, efficiency is increased by splitting each photon emitted by the Co-60 source a user-defined number of times. Only those split primary photons directed into a user-defined splitting field (encompassing the treatment field) are sampled, yielding many low-weight photons directed into the field. Efficiency can be further increased by taking advantage of radial symmetry at the top of the treatment head to reduce the number of split primary photons tracked in this portion. There is also an option to generate contaminant electrons in DSB.

RESULTS

The DSB scheme in BEAMnrc increases the photon fluence calculation efficiency in a 10 × 10 cm(2) Co-60 beam by a factor of 1800 with a concurrent increase in contaminant electron fluence calculation efficiency by a factor of 1200. Implementation of DSB in beampp, a C++ code for accelerator simulations based on EGSnrc and the C++ class library, egspp, increases photon fluence efficiency by a factor of 2800 and contaminant electron fluence efficiency by a factor of 1600. Optimum splitting numbers are in the range of 20,000-40,000. For dose calculations in a water phantom (0.5 × 0.5 × 0.5 cm(3) voxels) this translates into a factor of ∼400 increase in dose calculation efficiency (all doses > 0.5 × Dmax). An example calculation of the ratio of dose to water to dose to chamber (the basis of the beam quality correction factor) to within 0.2% in a realistic chamber using a full simulation of a Co-60 treatment head as a source indicates the practicality of Co-60 simulations with DSB.

CONCLUSIONS

The efficiency improvement resulting from DSB makes Monte Carlo commissioning of Co-60 beams and calculation of beam quality correction factors feasible.

Investigation of practical approaches to evaluating cumulative dose for cone beam computed tomography (CBCT) from standard CT dosimetry measurements: a Monte Carlo study

A function called Gx(L) was introduced by the International Commission on Radiation Units and Measurements (ICRU) Report-87 to facilitate measurement of cumulative dose for CT scans within long phantoms as recommended by the American Association of Physicists in Medicine (AAPM) TG-111. The Gx(L) function is equal to the ratio of the cumulative dose at the middle of a CT scan to the volume weighted CTDI (CTDIvol), and was investigated for conventional multi-slice CT scanners operating with a moving table. As the stationary table mode, which is the basis for cone beam CT (CBCT) scans, differs from that used for conventional CT scans, the aim of this study was to investigate the extension of the Gx(L) function to CBCT scans. An On-Board Imager (OBI) system integrated with a TrueBeam linac was simulated with Monte Carlo EGSnrc/BEAMnrc, and the absorbed dose was calculated within PMMA, polyethylene (PE), and water head and body phantoms using EGSnrc/DOSXYZnrc, where the body PE body phantom emulated the ICRU/AAPM phantom. Beams of width 40-500 mm and beam qualities at tube potentials of 80-140 kV were studied. Application of a modified function of beam width (W) termed Gx(W), for which the cumulative dose for CBCT scans f (0) is normalized to the weighted CTDI (CTDIw) for a reference beam of width 40 mm, was investigated as a possible option. However, differences were found in Gx(W) with tube potential, especially for body phantoms, and these were considered to be due to differences in geometry between wide beams used for CBCT scans and those for conventional CT. Therefore, a modified function Gx(W)100 has been proposed, taking the form of values of f (0) at each position in a long phantom, normalized with respect to dose indices f 100(150)x measured with a 100 mm pencil ionization chamber within standard 150 mm PMMA phantoms, using the same scanning parameters, beam widths and positions within the phantom. f 100(150)x averages the dose resulting from a CBCT scan over the 100 mm length. Like the Gx(L) function, the Gx(W)100 function showed only a weak dependency on tube potential at most positions for the phantoms studied. The results were fitted to polynomial equations from which f (0) within the longer PMMA, PE, or water phantoms can be evaluated from measurements of f 100(150)x. Comparisons with other studies, suggest that these functions may be suitable for application to any CT or CBCT scan acquired with stationary table mode.

Characterization of an orthovoltage biological irradiator used for radiobiological research

Orthovoltage irradiators are routinely used to irradiate specimens and small animals in biological research. There are several reports on the characteristics of these units for small field irradiations. However, there is limited knowledge about use of these units for large fields, which are essential for emerging large-field irregular shape irradiations, namely total marrow irradiation used as a conditioning regimen for hematological malignancies. This work describes characterization of a self-contained Orthovoltage biological irradiator for large fields using measurements and Monte Carlo simulations that could be used to compute the dose for in vivo or in vitro studies for large-field irradiation using this or a similar unit. Percentage depth dose, profiles, scatter factors, and half-value layers were measured and analyzed. A Monte Carlo model of the unit was created and used to generate depth dose and profiles, as well as scatter factors. An ion chamber array was also used for profile measurements of flatness and symmetry. The output was determined according to AAPM Task Group 61 guidelines. The depth dose measurements compare well with published data for similar beams. The Monte Carlo-generated depth dose and profiles match our measured doses to within 2%. Scatter factor measurements indicate gradual variation of these factors with field size. Dose rate measured by placing the ion chamber atop the unit's steel plate or solid water indicate enhanced readings of 5 to 28% compared with those measured in air. The stability of output over a 5-year period is within 2% of the 5-year average.

A Monte Carlo investigation of cumulative dose measurements for cone beam computed tomography (CBCT) dosimetry

Many studies have shown that the computed tomography dose index (CTDI100) which is considered as a main dose descriptor for CT dosimetry fails to provide a realistic reflection of the dose involved in cone beam computed tomography (CBCT) scans. Several practical approaches have been proposed to overcome drawbacks of the CTDI100. One of these is the cumulative dose concept. The purpose of this study was to investigate four different approaches based on the cumulative dose concept: the cumulative dose (1) f(0,150) and (2) f(0,∞) with a small ionization chamber 20 mm long, and the cumulative dose (3) f100(150) and (4) f100(∞) with a standard 100 mm pencil ionization chamber. The study also aimed to investigate the influence of using the 20 and 100 mm chambers and the standard and the infinitely long phantoms on cumulative dose measurements. Monte Carlo EGSnrc/BEAMnrc and EGSnrc/DOSXYZnrc codes were used to simulate a kV imaging system integrated with a TrueBeam linear accelerator and to calculate doses within cylindrical head and body PMMA phantoms with diameters of 16 cm and 32 cm, respectively, and lengths of 150, 600, 900 mm. f(0,150) and f100(150) approaches were studied within the standard PMMA phantoms (150 mm), while the other approaches f(0,∞) and f100(∞) were within infinitely long head (600 mm) and body (900 mm) phantoms. CTDI∞ values were used as a standard to compare the dose values for the approaches studied at the centre and periphery of the phantoms and for the weighted values. Four scanning protocols and beams of width 20-300 mm were used. It has been shown that the f(0,∞) approach gave the highest dose values which were comparable to CTDI∞ values for wide beams. The differences between the weighted dose values obtained with the 20 and 100 mm chambers were significant for the beam widths <120 mm, but these differences declined with increasing beam widths to be within 4%. The weighted dose values calculated within the infinitely long phantoms with both the chambers for the beam widths ≤140 were within 3% of those within the standard phantoms, but the differences rose to be within 15% at wider beams. By comparing the approaches studied in this investigation with other methodologies taking into account the efficiency of the approach as a dose descriptor and the simplicity of the implementation in the clinical environment, the f(0,150) method may be the best for CBCT dosimetry combined with the use of correction factors.

Monte Carlo model of the scanning beam digital x-ray (SBDX) source

The scanning-beam digital x-ray (SBDX) system has been developed for fluoroscopic imaging using an inverse x-ray imaging geometry. The SBDX system consists of a large-area x-ray source with a multihole collimator and a small detector. The goal of this study was to build a Monte Carlo (MC) model of the SBDX source as a useful tool for optimization of the SBDX imaging system in terms of its hardware components and imaging parameters. The MC model of the source was built in the EGSnrc/BEAMnrc code and validated using the DOSXYZnrc code and Gafchromic film measurements for 80, 100, and 120 kV x-ray source voltages. The MC simulated depth dose curves agreed with measurements to within 5%, and beam profiles at three selected depths generally agreed within 5%. Exposure rates and half-value layers for three voltages were also calculated from the MC simulations. Patient skin-dose per unit detector-dose was quantified as a function of patient size for all three x-ray source voltages. The skin-dose to detector-dose ratio ranged from 5-10 for a 20 cm thick patient to 1 × 10(3)-1 × 10(5) for a 50 cm patient for the 120 and 80 kV beams, respectively. Simulations of imaging dose for a prostate patient using common imaging parameters revealed that skin-dose per frame was as low as 0.2 mGy.

A PENELOPE-based system for the automated Monte Carlo simulation of clinacs and voxelized geometries-application to far-from-axis fields

PURPOSE

Two new codes, PENEASY and PENEASYLINAC, which automate the Monte Carlo simulation of Varian Clinacs of the 600, 1800, 2100, and 2300 series, together with their electron applicators and multileaf collimators, are introduced. The challenging case of a relatively small and far-from-axis field has been studied with these tools.

METHODS

PENEASY is a modular, general-purpose main program for the PENELOPE Monte Carlo system that includes various source models, tallies and variance-reduction techniques (VRT). The code includes a new geometry model that allows the superposition of voxels and objects limited by quadric surfaces. A variant of the VRT known as particle splitting, called fan splitting, is also introduced. PENEASYLINAC, in turn, automatically generates detailed geometry and configuration files to simulate linacs with PENEASY. These tools are applied to the generation of phase-space files, and of the corresponding absorbed dose distributions in water, for two 6 MV photon beams from a Varian Clinac 2100 C∕D: a 40 × 40 cm(2) centered field; and a 3 × 5 cm(2) field centered at (4.5, -11.5) cm from the beam central axis. This latter configuration implies the largest possible over-traveling values of two of the jaws. Simulation results for the depth dose and lateral profiles at various depths are compared, by using the gamma index, with experimental values obtained with a PTW 31002 ionization chamber. The contribution of several VRTs to the computing speed of the more demanding off-axis case is analyzed.

RESULTS

For the 40 × 40 cm(2) field, the percentages γ(1) and γ(1.2) of voxels with gamma indices (using 0.2 cm and 2% criteria) larger than unity and larger than 1.2 are 0.2% and 0%, respectively. For the 3 × 5 cm(2) field, γ(1) = 0%. These figures indicate an excellent agreement between simulation and experiment. The dose distribution for the off-axis case with voxels of 2.5 × 2.5 × 2.5 mm(3) and an average standard statistical uncertainty of 2% (1σ) is computed in 3.1 h on a single core of a 2.8 GHz Intel Core 2 Duo processor. This result is obtained with the optimal combination of the tested VRTs. In particular, fan splitting for the off-axis case accelerates execution by a factor of 240 with respect to standard particle splitting.

CONCLUSIONS

PENEASY and PENEASYLINAC can simulate the considered Varian Clinacs both in an accurate and efficient manner. Fan splitting is crucial to achieve simulation results for the off-axis field in an affordable amount of CPU time. Work to include Elekta linacs and to develop a graphical interface that will facilitate user input is underway.

Verification of a commercial implementation of the Macro-Monte-Carlo electron dose calculation algorithm using the virtual accelerator approach

In this work, the accuracy of the implementation of the Macro Monte Carlo electron dose calculation algorithm into the radiation therapy treatment planning system Eclipse is evaluated. This implementation - called eMC - uses a particle source based on the Rotterdam Initial Phase-Space model. A three-dimensional comparison of eMC calculated dose to dose distributions resulting from full treatment head simulations with the Monte Carlo code package EGSnrc is performed using the 'virtual accelerator' approach. Calculated dose distributions are compared for a homogeneous tissue equivalent phantom and a water phantom with air and bone inhomogeneities. The performance of the eMC algorithm in both phantoms can be considered acceptable within the 2%/2 mm Gamma index criterion. A systematic underestimation of dose by the eMC algorithm within the air inhomogeneity is found.

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 simulation and patient dosimetry for a kilovoltage cone-beam CT unit

PURPOSE

The purpose of this work is to characterize the x-ray volume imager (XVI), the cone-beam computed tomography (CBCT) unit mounted on the Elekta Synergy linac, with F1 bowtie filter and to calculate the three-dimensional dose delivered to patients using volumetric acquisition.

METHODS

The XVI is modeled in detail using a new Monte Carlo (MC) code, BEAMPP, under development at the National Research Council Canada. In this investigation, a new component module is developed to accurately model the unit's bowtie filter used in conjunction with the available beam collimators at the clinical energy of 120 kV. The modeling is compared against percentage depth dose (PDD) and profile measurements. Kilovoltage radiation beams' phase space files are also analyzed. The authors also describe a method for the absolute dose calibration of the MC model of the CBCT unit when used in a clinical volumetric acquisition mode. Finally, they calculate three-dimensional patient dose from CBCT image acquisition in three clinical cases of interest: Pelvis, lung, and head and neck.

RESULTS

The agreement between measurement and MC is shown to be very good: Within +/- 2% for the PDD and within +/- 3.5% inside the radiation field for all the collimators with the F1 bowtie filter. A full account of the absolute calibration method is given and dose calculation is validated against ion chamber measurements in different locations of a plastic phantom. Calculations and experiments agree within +/- 2% or better in both at the center and the periphery of the phantom, with worst agreement of 4.5% at the surface of the phantom and for one specific combination of collimator and filter. Patient dose from CBCT scan reveals that dose to tissue is between 2 and 2.5 cGy for a pelvis or a lung full acquisition. For H&N dose to tissue is 5 cGy, with the unit presets used in this work. Dose to bony structures can be two to three times higher than dose to tissue.

CONCLUSIONS

The XVI CBCT unit has been fully modeled including the F1 bowtie filter. Absolute dose distribution from the unit has been successfully validated. Full MC patient dose calculation has shown that the three-dimensional dose distribution from CBCT is complex. Patient dose from CBCT exposure cannot be completely accounted for by using a numerical factor as an estimate of the dose at the center of the body. Furthermore, additional dose to bone should be taken into account when adopting any IGRT strategy and weighed vs the unquestionable benefits of the technique in order to optimize treatment. Full three-dimensional dose calculation is recommended if patient dose from CBCT is to be integrated in any adaptive planning strategy.

An enhanced HOWFARLESS option for DOSXYZnrc simulations of slab geometries

The Monte Carlo code DOSXYZnrc is a valuable instrument for calculating absorbed dose within a three-dimensional Cartesian geometry. DOSXYZnrc includes several variance reduction techniques used to increase the efficiency of the Monte Carlo calculation. One such technique is HOWFARLESS which is used to increase the efficiency of beam commissioning calculations in homogeneous phantoms. The authors present an enhanced version of HOWFARLESS which extends the application to include phantoms inhomogeneous in one dimension. When the enhanced HOWFARLESS was used, efficiency increases as high as 14 times were observed without any loss in dose accuracy. The efficiency gains of an enhanced HOWFARLESS simulation was found to be dependent on both slab geometry and slab density. As the number of two-dimensional voxel layers per slab increases, so does the efficiency gain. Also, as the mass density of a slab is decreased, the efficiency gains increase.