Clinical use of a commercial Monte Carlo treatment planning system for electron beams.

In 2002 we fully implemented clinically a commercial Monte Carlo based treatment planning system for electron beams. The software, developed by MDS Nordion (presently Nucletron), is based on Kawrakow's VMC++ algorithm. The Monte Carlo module is integrated with our Theraplan Plustrade mark treatment planning system. An extensive commissioning process preceded clinical implementation of this software. Using a single virtual 'machine' for each electron beam energy, we can now calculate very accurately the dose distributions and the number of MU for any arbitrary field shape and SSD. This new treatment planning capability has significantly impacted our clinical practice. Since we are more confident of the actual dose delivered to a patient, we now calculate accurate three-dimensional (3D) dose distributions for a greater variety of techniques and anatomical sites than we have in the past. We use the Monte Carlo module to calculate dose for head and neck, breast, chest wall and abdominal treatments with electron beams applied either solo or in conjunction with photons. In some cases patient treatment decisions have been changed, as compared to how such patients would have been treated in the past. In this paper, we present the planning procedure and some clinical examples.

Monte Carlo-aided dosimetry of the new Bebig IsoSeed 103Pd interstitial brachytherapy seed.

A new model 103Pd interstitial brachytherapy source, the IsoSeed 103Pd, was recently introduced by Bebig Isotopentechnik und Umweltdiagnostik GmbH for permanent implant applications. This study presents the first quantitative theoretical study of the seed's dosimetric quantities. Monte Carlo photon transport (MCPT) simulation techniques have been used to evaluate the dose-rate distributions around the model IsoSeed 103Pd source in liquid water and air phantoms. These results have been used to calculate and tabulate the anisotropy function, F(r, theta), radial dose function, g(r), and anisotropy factors, phi(r), and dose-rate constant as defined by AAPM Task Group 43 (TG-43) Report. Cartesian "away" and "along" tables, giving the dose rates per unit air-kerma strength in water in the range 0.1-3 cm distance around the seed have also been tabulated. The dose-rate constant, lambda, was evaluated by simulating the wide-angle, free-air chamber (WAFAC) calibration geometry recently implemented by NIST (National Institute of Standards and Technology) to realize the primary standard of air-kerma strength (SK,N99) for low-energy photon-emitting brachytherapy sources. The dose-rate constant has been found to be lambda=0.660+/-0.017 in units of dose-rate per unit air-kerma strength (cGy x h(-1) x U(-1)).

Monte Carlo-aided dosimetry of a new high dose-rate brachytherapy source.

In this article we introduce a new high-intensity 192Ir source design for use in a recently reengineered microSelectron-HDR remote afterloading device for high dose-rate (HDR) brachytherapy. The maximum rigid length and outer diameter of the new source are reduced to 4.95 and 0.90 mm, respectively, compared to 5.50 and 1.10 mm for the previous source design introduced in 1991. In addition, a smaller diameter and more flexible steel cable are used, allowing the source cable to negotiate smaller diameter catheters or more tortuously curved catheters. Using Monte Carlo photon transport simulation, the complete two-dimensional (2-D) dose-rate distribution is calculated over the 0.1-7 cm distance range and are presented both as conventional 2-D Cartesian lookup tables and in the formalism recommended by the American Association of Physicists in Medicine Task Group 43 (TG-43) Report. The dose distribution of this source is very similar to that of its predecessor, except near the source tip and in the shadow of the cable assembly, where differences of 5%-8% are apparent. The accuracy of various methods for extrapolating beyond the tubulated anisotropy functions to short distances is evaluated. It is demonstrated that linear extrapolation from the anisotropy functions defined by TG-43 accurately (+/- 2%) estimates dose rate at short and long distances lying outside the radial distance range of the original measured data from which the anisotropy and radial dose functions were derived. In contrast, the algorithm used on the vendor's planning system results in large calculation errors at distances less than 5 mm.

Analytical approach to heterogeneity correction factor calculation for brachytherapy.

In brachytherapy treatment planning, the effects of tissue and applicator heterogeneities are commonly neglected due to lack of accurate, general, and fast three-dimensional (3D) dose-computational algorithms. A novel approach, based on analytical calculation of scattered photon fluxes inside and around a disk-shaped heterogeneity, has been developed for use in the three-dimensional scatter-subtraction algorithm. Specifically, our model predicts the central-ray dose distribution for a collimated photon isotropic source or brachytherapy "minibeam" in the presence of a slab of heterogeneous material. The model accounts for the lateral dimensions, location, composition, density, and thickness of the heterogeneity using precalculated scatter-to-primary ratios (SPRs) for the corresponding homogeneous problem. The model is applicable to the entire brachytherapy energy range (25 to 662 keV) and to a broad range of materials having atomic numbers of 13 to 82, densities of 2.7 g.cm-3 (Al) to 21.45 g.cm-3 (Pt) and thicknesses up to 1 mean free path. For this range of heterogeneous materials, the heterogeneity correction factors (HCFs) vary from 0.09 to 0.75. The model underestimates HCF when multiple scattering prevails and overestimates HCF when absorption dominates. However, the analytic model agrees with Monte Carlo photon transport (MCPT) benchmark calculations within 1.8% to 10% for 125I, 169Yb, 192Ir, and 137Cs for a wide variety of materials, with the exception of Ag. For 125I shielded by Ag, where the mean discrepancy can exceed 25%, the error is due to K-edge characteristic x rays originating within the heterogeneity. The proposed approach provides reductions in CPU time required of 5 x 10(4)-10(5) and 100 in comparison with direct MCPT simulation and 1D numerical integration, respectively. The limitations of model applicability, as determined by the physical properties of heterogeneity material and accuracy required, are also discussed.

Evaluation of the first commercial Monte Carlo dose calculation engine for electron beam treatment planning.

The purpose of this study is to perform a clinical evaluation of the first commercial (MDS Nordion, now Nucletron) treatment planning system for electron beams incorporating Monte Carlo dose calculation module. This software implements Kawrakow's VMC++ voxel-based Monte Carlo calculation algorithm. The accuracy of the dose distribution calculations is evaluated by direct comparisons with extensive sets of measured data in homogeneous and heterogeneous phantoms at different source-to-surface distances (SSDs) and gantry angles. We also verify the accuracy of the Monte Carlo module for monitor unit calculations in comparison with independent hand calculations for homogeneous water phantom at two different SSDs. All electron beams in the range 6-20 MeV are from a Siemens KD-2 linear accelerator. We used 10,000 or 50,000 histories/cm2 in our Monte Carlo calculations, which led to about 2.5% and 1% relative standard error of the mean of the calculated dose. The dose calculation time depends on the number of histories, the number of voxels used to map the patient anatomy, the field size, and the beam energy. The typical run time of the Monte Carlo calculations (10,000 histories/cm2) is 1.02 min on a 2.2 GHz Pentium 4 Xeon computer for a 9 MeV beam, 10 x 10 cm2 field size, incident on the phantom 15 x 15 x 10 cm3 consisting of 31 CT slices and voxels size of 3 x 3 x 3 mm3 (total of 486,720 voxels). We find good agreement (discrepancies smaller than 5%) for most of the tested dose distributions. We also find excellent agreement (discrepancies of 2.5% or less) for the monitor unit calculations relative to the independent manual calculations. The accuracy of monitor unit calculations does not depend on the SSD used, which allows the use of one virtual machine for each beam energy for all arbitrary SSDs. In some cases the test results are found to be sensitive to the voxel size applied such that bigger systematic errors (>5%) occur when large voxel sizes interfere with the extensions of heterogeneities or dose gradients because of differences between the experimental and calculated geometries. Therefore, user control over voxelization is important for high accuracy electron dose calculations.

Dosimetric modeling of the microselectron high-dose rate 192Ir source by the multigroup discrete ordinates method.

The DANTSYS multigroup discrete ordinates computer code is applied to quantitatively estimate the absorbed dose rate distributions in the vicinity of a microSelectron 192Ir high-dose-rate (HDR) source in two-dimensional cylindrical R-Z geometry. The source is modeled in a cylindrical water phantom of diameter 20 cm and height 20 cm. The results are also used for evaluation of the Task Group 43 (TG-43) dosimetric quantities. The DANTSYS accuracy is estimated by direct comparisons with corresponding Monte Carlo results. Our 210-group photon cross section library developed previously, together with angular quadratures consisting of 36 (S16) to 210 (S40) directions and associated weights per octant, are used in the DANTSYS simulations. Strong ray effects are observed but are significantly mitigated through the use of DANTSYS's stochastic ray-tracing first collision source algorithm. The DANTSYS simulations closely approximate Monte Carlo estimates of both direct dose calculations and TG-43 dosimetric quantities. The discrepancies with S20 angular quadrature (55 directions and weights per octant) or higher are shown to be less than +/- 5% (about 2.5 standard deviations of Monte Carlo calculations) everywhere except for limited regions along the Z axis of rotational symmetry, where technical limitations in the DANTSYS first collision source implementation makes adequate suppression of ray effects difficult to achieve. The efficiency of DANTSYS simulations is compared with that of the EGS4 Monte Carlo code. It is demonstrated that even with the 210-group cross section library, DANTSYS achieves two-fold efficiency gains using the the S20 quadrature set. The potential of discrete ordinates method for further efficiency improvements is also discussed.