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.
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.
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.
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.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.