Rowan M. Thomson

Report of the Task Group 186 on model-based dose calculation methods in brachytherapy beyond the TG-43 formalism: Current status and recommendations for clinical implementation

The charge of Task Group 186 (TG-186) is to provide guidance for early adopters of model-based dose calculation algorithms (MBDCAs) for brachytherapy (BT) dose calculations to ensure practice uniformity. Contrary to external beam radiotherapy, heterogeneity correction algorithms have only recently been made available to the BT community. Yet, BT dose calculation accuracy is highly dependent on scatter conditions and photoelectric effect cross-sections relative to water. In specific situations, differences between the current water-based BT dose calculation formalism (TG-43) and MBDCAs can lead to differences in calculated doses exceeding a factor of 10. MBDCAs raise three major issues that are not addressed by current guidance documents: (1) MBDCA calculated doses are sensitive to the dose specification medium, resulting in energy-dependent differences between dose calculated to water in a homogeneous water geometry (TG-43), dose calculated to the local medium in the heterogeneous medium, and the intermediate scenario of dose calculated to a small volume of water in the heterogeneous medium. (2) MBDCA doses are sensitive to voxel-by-voxel interaction cross sections. Neither conventional single-energy CT nor ICRU∕ICRP tissue composition compilations provide useful guidance for the task of assigning interaction cross sections to each voxel. (3) Since each patient-source-applicator combination is unique, having reference data for each possible combination to benchmark MBDCAs is an impractical strategy. Hence, a new commissioning process is required. TG-186 addresses in detail the above issues through the literature review and provides explicit recommendations based on the current state of knowledge. TG-43-based dose prescription and dose calculation remain in effect, with MBDCA dose reporting performed in parallel when available. In using MBDCAs, it is recommended that the radiation transport should be performed in the heterogeneous medium and, at minimum, the dose to the local medium be reported along with the TG-43 calculated doses. Assignments of voxel-by-voxel cross sections represent a particular challenge. Electron density information is readily extracted from CT imaging, but cannot be used to distinguish between different materials having the same density. Therefore, a recommendation is made to use a number of standardized materials to maintain uniformity across institutions. Sensitivity analysis shows that this recommendation offers increased accuracy over TG-43. MBDCA commissioning will share commonalities with current TG-43-based systems, but in addition there will be algorithm-specific tasks. Two levels of commissioning are recommended: reproducing TG-43 dose parameters and testing the advanced capabilities of MBDCAs. For validation of heterogeneity and scatter conditions, MBDCAs should mimic the 3D dose distributions from reference virtual geometries. Potential changes in BT dose prescriptions and MBDCA limitations are discussed. When data required for full MBDCA implementation are insufficient, interim recommendations are made and potential areas of research are identified. Application of TG-186 guidance should retain practice uniformity in transitioning from the TG-43 to the MBDCA approach.

Comparison of dose calculation methods for brachytherapy of intraocular tumors

PURPOSE To investigate dosimetric differences among several clinical treatment planning systems (TPS) and Monte Carlo (MC) codes for brachytherapy of intraocular tumors using 125I or 103Pd plaques, and to evaluate the impact on the prescription dose of the adoption of MC codes and certain versions of a TPS (Plaque Simulator with optional modules). METHODS Three clinical brachytherapy TPS capable of intraocular brachytherapy treatment planning and two MC codes were compared. The TPS investigated were Pinnacle v8.0dp1, BrachyVision v8.1, and Plaque Simulator v5.3.9, all of which use the AAPM TG-43 formalism in water. The Plaque Simulator software can also handle some correction factors from MC simulations. The MC codes used are MCNP5 v1.40 and BrachyDose/EGSnrc. Using these TPS and MC codes, three types of calculations were performed: homogeneous medium with point sources (for the TPS only, using the 1D TG-43 dose calculation formalism); homogeneous medium with line sources (TPS with 2D TG-43 dose calculation formalism and MC codes); and plaque heterogeneity-corrected line sources (Plaque Simulator with modified 2D TG-43 dose calculation formalism and MC codes). Comparisons were made of doses calculated at points-of-interest on the plaque central-axis and at off-axis points of clinical interest within a standardized model of the right eye. RESULTS For the homogeneous water medium case, agreement was within approximately 2% for the point- and line-source models when comparing between TPS and between TPS and MC codes, respectively. For the heterogeneous medium case, dose differences (as calculated using the MC codes and Plaque Simulator) differ by up to 37% on the central-axis in comparison to the homogeneous water calculations. A prescription dose of 85 Gy at 5 mm depth based on calculations in a homogeneous medium delivers 76 Gy and 67 Gy for specific 125I and 103Pd sources, respectively, when accounting for COMS-plaque heterogeneities. For off-axis points-of-interest, dose differences approached factors of 7 and 12 at some positions for 125I and 103Pd, respectively. There was good agreement (approximately 3%) among MC codes and Plaque Simulator results when appropriate parameters calculated using MC codes were input into Plaque Simulator. Plaque Simulator and MC users are perhaps at risk of overdosing patients up to 20% if heterogeneity corrections are used and the prescribed dose is not modified appropriately. CONCLUSIONS Agreement within 2% was observed among conventional brachytherapy TPS and MC codes for intraocular brachytherapy dose calculations in a homogeneous water environment. In general, the magnitude of dose errors incurred by ignoring the effect of the plaque backing and Silastic insert (i.e., by using the TG-43 approach) increased with distance from the plaques central-axis. Considering the presence of material heterogeneities in a typical eye plaque, the best method in this study for dose calculations is a verified MC simulation.

Monte Carlo dosimetry for and eye plaque brachytherapy

A Monte Carlo study of dosimetry for eye plaque brachytherapy is performed. BrachyDose, an EGSnrc user code which makes use of Yegins multi-geometry package, is used to fully model I125 (model 6711) and Pd103 (model 200) brachytherapy seeds and the standardized plaques of the Collaborative Ocular Melanoma Study (COMS). Three-dimensional dose distributions in the eye region are obtained. In general, dose to water is scored; however, the implications of replacing water with eye tissues are explored. The effect of the gold alloy (Modulay) backing is investigated and the dose is found to be sensitive to the elemental composition of the backing. The presence of the silicone polymer (Silastic) seed carrier results in substantial dose decreases relative to water, particularly for Pd103. For a 20mm plaque with a Modulay backing and Silastic insert, fully loaded with 24 seeds, the dose decrease relative to water is of the order of 14% for I125 and 20% for Pd103 at a distance of 1cm from the inner sclera along the plaques central axis. For the configurations of seeds used in COMS plaques, interseed attenuation is a small effect within the eye region. The introduction of an air interface results in a dose reduction in its vicinity which depends on the plaques position within the eye and the radionuclide. Introducing bone in the eyes vicinity also causes dose reductions. The dose distributions in the eye for the two different radionuclides are compared and, for the same prescription dose, Pd103 generally offers a lower dose to critical normal structures. BrachyDose is sufficiently fast to allow full Monte Carlo dose calculations for routine clinical treatment planning.

Monte Carlo dosimetry for 125I and 103Pd eye plaque brachytherapy.

A Monte Carlo study of dosimetry for eye plaque brachytherapy is performed. BrachyDose, an EGSnrc user code which makes use of Yegins multi-geometry package, is used to fully model I125 (model 6711) and Pd103 (model 200) brachytherapy seeds and the standardized plaques of the Collaborative Ocular Melanoma Study (COMS). Three-dimensional dose distributions in the eye region are obtained. In general, dose to water is scored; however, the implications of replacing water with eye tissues are explored. The effect of the gold alloy (Modulay) backing is investigated and the dose is found to be sensitive to the elemental composition of the backing. The presence of the silicone polymer (Silastic) seed carrier results in substantial dose decreases relative to water, particularly for Pd103. For a 20mm plaque with a Modulay backing and Silastic insert, fully loaded with 24 seeds, the dose decrease relative to water is of the order of 14% for I125 and 20% for Pd103 at a distance of 1cm from the inner sclera along the plaques central axis. For the configurations of seeds used in COMS plaques, interseed attenuation is a small effect within the eye region. The introduction of an air interface results in a dose reduction in its vicinity which depends on the plaques position within the eye and the radionuclide. Introducing bone in the eyes vicinity also causes dose reductions. The dose distributions in the eye for the two different radionuclides are compared and, for the same prescription dose, Pd103 generally offers a lower dose to critical normal structures. BrachyDose is sufficiently fast to allow full Monte Carlo dose calculations for routine clinical treatment planning.

A generic high‐dose rate 192Ir brachytherapy source for evaluation of model‐based dose calculations beyond the TG‐43 formalism

PURPOSE In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) (192)Ir source and a virtual water phantom were designed, which can be imported into a TPS. METHODS A hypothetical, generic HDR (192)Ir source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic (192)Ir source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra(®) Brachy with advanced collapsed-cone engine (ACE) and BrachyVision ACUROS™ ]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and PENELOPE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201)(3) voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR (192)Ir source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods. RESULTS TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ace algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 ± 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ace at clinically relevant distances. CONCLUSIONS A hypothetical, generic HDR (192)Ir source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.

Monte Carlo dosimetry for 125I and 103Pd eye plaque brachytherapy with various seed models

A Monte Carlo study of dosimetry for eye plaque brachytherapy is performed. BrachyDose, an EGSnrc user code which makes use of Yegins multi-geometry package, is used to fully model 125I (model 6711) and 103Pd (model 200) brachytherapy seeds and the standardized plaques of the Collaborative Ocular Melanoma Study (COMS). Three-dimensional dose distributions in the eye region are obtained. In general, dose to water is scored; however, the implications of replacing water with eye tissues are explored. The effect of the gold alloy (Modulay) backing is investigated and the dose is found to be sensitive to the elemental composition of the backing. The presence of the silicone polymer (Silastic) seed carrier results in substantial dose decreases relative to water, particularly for 103Pd. For a 20 mm plaque with a Modulay backing and Silastic insert, fully loaded with 24 seeds, the dose decrease relative to water is of the order of 14% for 125I and 20% for 103Pd at a distance of 1 cm from the inner sclera along the plaques central axis. For the configurations of seeds used in COMS plaques, interseed attenuation is a small effect within the eye region. The introduction of an air interface results in a dose reduction in its vicinity which depends on the plaques position within the eye and the radionuclide. Introducing bone in the eyes vicinity also causes dose reductions. The dose distributions in the eye for the two different radionuclides are compared and, for the same prescription dose, 103Pd generally offers a lower dose to critical normal structures. BrachyDose is sufficiently fast to allow full Monte Carlo dose calculations for routine clinical treatment planning.

Sci—Sat AM(2): Brachy — 05: Fast Monte Carlo Dose Calculations for Brachytherapy with BrachyDose

A fast dose calculation algorithm called BrachyDose has been developed for brachytherapy applications. BrachyDose is based on the EGSnrc code system for simulating radiation transport. Complex geometries are modelled through the superposition of basic geometric entities (spheres, cuboids, cylinders, and cones) using Yegins multi-geometry package; the phantom geometry may be defined using a CT dataset. A database of brachytherapy sources has been developed and benchmarked, as has a database of eye plaque applicators. BrachyDose scores collision kerma, which is equivalent to absorbed dose for most situations of interest, using a tracklength estimator. The phase space of particles emitted from brachytherapy sources may be generated with BrachyDose and used in subsequent simulations to avoid the repeated simulation of particle transport within sources. A particle recycling feature has been implemented for multisource configurations in which the first source acts as a particle generator; particles emitted from this source are reinitiated at each source location. Dose calculations for prostate permanent implants achieving 2% average uncertainty in the prostate region take less than 30 seconds in (2 mm)3 voxels on a single 3.0 GHz Woodcrest core; calculation times for eye plaque therapy are on the order of three minutes in (0.5 mm)3 voxels. These calculation times are sufficiently fast for routine clinical treatment planning. A graphical user interface (GUI) for BrachyDose has been developed. Working towards clinical implementation, efforts are underway to integrate data in the DICOM-RT format with BrachyDose.

Model-based dose calculations for 125I lung brachytherapy

PURPOSE Model-baseddose calculations (MBDCs) are performed using patient computed tomography (CT) data for patients treated with intraoperative (125)I lung brachytherapy at the Mayo Clinic Rochester. Various metallic artifact correction and tissue assignment schemes are considered and their effects on dose distributions are studied. Dose distributions are compared to those calculated under TG-43 assumptions. METHODS Dose distributions for six patients are calculated using phantoms derived from patient CT data and the EGSnrc user-code BrachyDose. (125)I (GE Healthcare/Oncura model 6711) seeds are fully modeled. Four metallic artifact correction schemes are applied to the CT data phantoms: (1) no correction, (2) a filtered back-projection on a modified virtual sinogram, (3) the reassignment of CT numbers above a threshold in the vicinity of the seeds, and (4) a combination of (2) and (3). Tissue assignment is based on voxel CT number and mass density is assigned using a CT number to mass density calibration. Three tissue assignment schemes with varying levels of detail (20, 11, and 5 tissues) are applied to metallic artifact corrected phantoms. Simulations are also performed under TG-43 assumptions, i.e., seeds in homogeneous water with no interseed attenuation. RESULTS Significant dose differences (up to 40% for D(90)) are observed between uncorrected and metallic artifact corrected phantoms. For phantoms created with metallic artifact correction schemes (3) and (4), dose volume metrics are generally in good agreement (less than 2% differences for all patients) although there are significant local dose differences. The application of the three tissue assignment schemes results in differences of up to 8% for D(90); these differences vary between patients. Significant dose differences are seen between fully modeled and TG-43 calculations with TG-43 underestimating the dose (up to 36% in D(90)) for larger volumes containing higher proportions of healthy lung tissue. CONCLUSIONS Metallic artifact correction is necessary for accurate application of MBDCs for lung brachytherapy; simpler threshold replacement methods may be sufficient for early adopters concerned with clinical dose metrics. Rigorous determination of voxel tissue parameters and tissue assignment is required for accurate dose calculations as different tissue assignment schemes can result in significantly different dose distributions. Significant differences are seen between MBDCs and TG-43 dose distributions with TG-43 underestimating dose in volumes containing healthy lung tissue.PURPOSE Model-baseddose calculations (MBDCs) are performed using patient computed tomography (CT) data for patients treated with intraoperative125 I lung brachytherapy at the Mayo Clinic Rochester. Various metallic artifact correction and tissue assignment schemes are considered and their effects on dose distributions are studied. Dose distributions are compared to those calculated under TG-43 assumptions. METHODS Dose distributions for six patients are calculated using phantoms derived from patient CT data and the EGSnrc user-code BrachyDose.125 I (GE Healthcare/Oncura model 6711) seeds are fully modeled. Four metallic artifact correction schemes are applied to the CT data phantoms: (1) no correction, (2) a filtered back-projection on a modified virtual sinogram, (3) the reassignment of CT numbers above a threshold in the vicinity of the seeds, and (4) a combination of (2) and (3). Tissue assignment is based on voxel CT number and mass density is assigned using a CT number to mass density calibration. Three tissue assignment schemes with varying levels of detail (20, 11, and 5 tissues) are applied to metallic artifact corrected phantoms. Simulations are also performed under TG-43 assumptions, i.e., seeds in homogeneous water with no interseed attenuation. RESULTS Significant dose differences (up to 40% for D90 ) are observed between uncorrected and metallic artifact corrected phantoms. For phantoms created with metallic artifact correction schemes (3) and (4), dose volume metrics are generally in good agreement (less than 2% differences for all patients) although there are significant local dose differences. The application of the three tissue assignment schemes results in differences of up to 8% for D90 ; these differences vary between patients. Significant dose differences are seen between fully modeled and TG-43 calculations with TG-43 underestimating the dose (up to 36% in D90 ) for larger volumes containing higher proportions of healthy lung tissue. CONCLUSIONS Metallic artifact correction is necessary for accurate application of MBDCs for lung brachytherapy; simpler threshold replacement methods may be sufficient for early adopters concerned with clinical dose metrics. Rigorous determination of voxel tissue parameters and tissue assignment is required for accurate dose calculations as different tissue assignment schemes can result in significantly different dose distributions. Significant differences are seen between MBDCs and TG-43 dose distributions with TG-43 underestimating dose in volumes containing healthy lung tissue.

Modified COMS plaques for 125I and 103Pd iris melanoma brachytherapy.

PURPOSE Novel plaques are used to treat iris melanoma at the Mayo Clinic Rochester. The plaques are a modification of the Collaborative Ocular Melanoma Study (COMS) 22 mm plaque design with a gold alloy backing, outer lip, and silicone polymer insert. An inner lip surrounds a 10 mm diameter cutout region at the plaque center. Plaques span 360°, 270°, and 180° arcs. This article describes dosimetry for these plaques and others used in the treatment of anterior eye melanomas. METHODS AND MATERIALS The EGSnrc user-code BrachyDose is used to perform Monte Carlo simulations. Plaques and seeds are fully modeled. Three-dimensional dose distributions for different plaque models, TG-43 calculations, and (125)I (model 6711) and (103)Pd (model 200) seeds are compared via depth-dose curves, tabulation of doses at points of interest, and isodose contours. RESULTS Doses at points of interest differ by up to 70% from TG-43 calculations. The inner lip reduces corneal doses. Matching plaque arc length to tumor extent reduces doses to eye regions outside the treatment area. Maintaining the same prescription dose, (103)Pd offers lower doses to critical structures than (125)I, with the exception of the sclera adjacent to the plaque. CONCLUSION The Mayo Clinic plaques offer several advantages for anterior eye tumor treatments. Doses to regions outside the treatment area are significantly reduced. Doses differ considerably from TG-43 predictions, illustrating the importance of complete Monte Carlo simulations. Calculations take a few minutes on a single CPU, making BrachyDose sufficiently fast for routine clinical treatment planning.

Model-based dose calculations for COMS eye plaque brachytherapy using an anatomically realistic eye phantom.

PURPOSE To investigate the effects of the composition and geometry of ocular media and tissues surrounding the eye on dose distributions for COMS eye plaque brachytherapy with(125)I, (103)Pd, or (131)Cs seeds, and to investigate doses to ocular structures. METHODS An anatomically and compositionally realistic voxelized eye model with a medial tumor is developed based on a literature review. Mass energy absorption and attenuation coefficients for ocular media are calculated. Radiation transport and dose deposition are simulated using the EGSnrc Monte Carlo user-code BrachyDose for a fully loaded COMS eye plaque within a water phantom and our full eye model for the three radionuclides. A TG-43 simulation with the same seed configuration in a water phantom neglecting the plaque and interseed effects is also performed. The impact on dose distributions of varying tumor position, as well as tumor and surrounding tissue media is investigated. Each simulation and radionuclide is compared using isodose contours, dose volume histograms for the lens and tumor, maximum, minimum, and average doses to structures of interest, and doses to voxels of interest within the eye. RESULTS Mass energy absorption and attenuation coefficients of the ocular media differ from those of water by as much as 12% within the 20-30 keV photon energy range. For all radionuclides studied, average doses to the tumor and lens regions in the full eye model differ from those for the plaque in water by 8%-10% and 13%-14%, respectively; the average doses to the tumor and lens regions differ between the full eye model and the TG-43 simulation by 2%-17% and 29%-34%, respectively. Replacing the surrounding tissues in the eye model with water increases the maximum and average doses to the lens by 2% and 3%, respectively. Substituting the tumor medium in the eye model for water, soft tissue, or an alternate melanoma composition affects tumor dose compared to the default eye model simulation by up to 16%. In the full eye model simulations, the average dose to the lens is larger by 7%-9% than the dose to the center of the lens, and the maximum dose to the optic nerve is 17%-22% higher than the dose to the optic disk for all radionuclides. In general, when normalized to the same prescription dose at the tumor apex, doses delivered to all structures of interest in the full eye model are lowest for(103)Pd and highest for (131)Cs, except for the tumor where the average dose is highest for (103)Pd and lowest for (131)Cs. CONCLUSIONS The eye is not radiologically water-equivalent, as doses from simulations of the plaque in the full eye model differ considerably from doses for the plaque in a water phantom and from simulated TG-43 calculated doses. This demonstrates the importance of model-based dose calculations for eye plaque brachytherapy, for which accurate elemental compositions of ocular media are necessary.