Differential dose contributions on total dose distribution of 125I brachytherapy source

This work provides an improvement of the approach using Monte Carlo simulation for the Amersham Model 6711 ¹²⁵I brachytherapy seed source, which is well known by many theoretical and experimental studies. The source which has simple geometry was researched with respect to criteria of AAPM Tg-43 Report. The approach offered by this study involves determination of differential dose contributions that come from virtual partitions of a massive radioactive element of the studied source to a total dose at analytical calculation point. Some brachytherapy seeds contain multi-radioactive elements so the dose at any point is a total of separate doses from each element. It is momentous to know well the angular and radial dose distributions around the source that is located in cancerous tissue for clinical treatments. Interior geometry of a source is effective on dose characteristics of a distribution. Dose information of inner geometrical structure of a brachytherapy source cannot be acquired by experimental methods because of limits of physical material and geometry in the healthy tissue, so Monte Carlo simulation is a required approach of the study. EGSnrc Monte Carlo simulation software was used. In the design of a simulation, the radioactive source was divided into 10 rings, partitioned but not separate from each other. All differential sources were simulated for dose calculation, and the shape of dose distribution was determined comparatively distribution of a single-complete source. In this work anisotropy function was examined also mathematically.     

A Monte Carlo evaluation for effects of probable dimensional uncertainties of low dose rate brachytherapy seeds on dose

The aim of this study is to determine effects of size deviations of brachytherapy seeds on two dimensional dose distributions around the seed. Although many uncertainties are well known, the uncertainties which stem from geometric features of radiation sources are weakly considered and predicted. Neither TG-43 report which is not completely in common consensus, nor individual scientific MC and experimental studies include sufficient data for geometric uncertainties. Sizes of seed and its components can vary in a manufacturing deviation. This causes geometrical uncertainties, too. In this study, three seeds which have different geometrical properties were modeled using EGSnrc-Code Packages. Seeds were designed with all their details using the geometry package. 5% deviations of seed sizes were assumed. Modified seeds were derived from original seed by changing sizes by 5%. Normalizations of doses which were calculated from three kinds of brachytherapy seed and their derivations were found to be about 3%–20%. It was shown that manufacturing differences of brachytherapy seed cause considerable changes in dose distribution.     

Comparison of planned and measured rectal dose in-vivo during high dose rate Cobalt-60 brachytherapy of cervical cancer

Cobalt-60 (Co-60) is a relatively new source for the application of high-dose rate (HDR) brachytherapy. Radiation dose to the rectum is often a limiting factor in achieving the full prescribed dose to the target during brachytherapy of cervical cancer. The aim of this study was to measure radiation doses to the rectum in-vivo during HDR Co-60 brachytherapy. A total of eleven HDR brachytherapy treatments of cervical cancer were recruited in this study. A series of diodes incorporated in a rectal probe was inserted into the patient's rectum during each brachytherapy procedure. Real-time measured rectal doses were compared to calculated doses by the treatment planning system (TPS). The differences between calculated and measured dose ranged from 8.5% to 41.2%. This corresponds to absolute dose differences ranging from 0.3 Gy to 1.5 Gy. A linear relationship was observed between calculated and measured doses with linear regression R² value of 0.88, indicating close association between the measured and calculated doses. In general, absorbed doses for the rectum as calculated by TPS were observed to be higher than the doses measured using the diode probe. In-vivo dosimetry is an important quality assurance method for HDR brachytherapy of cervical cancer. It provides information that can contribute to the reduction of errors and discrepancies in dose delivery. Our study has shown that in-vivo dosimetry is feasible and can be performed to estimate the dose to the rectum during HDR brachytherapy using Co-60.     

Dosimetric characterization of a novel 90Y source for use in the conformal superficial brachytherapy device

PurposeTo characterize the dose distribution in water of a novel beta-emitting brachytherapy source for use in a Conformal Superficial Brachytherapy (CSBT) device.Methods and materialsYttrium-90 (⁹⁰Y) sources were designed for use with a uniquely designed CSBT device. Depth dose and planar dose measurements were performed for bare sources and sources housed within a 3D printed source holder. Monte Carlo simulated dose rate distributions were compared to film-based measurements. Gamma analysis was performed to compare simulated and measured dose rates from seven ⁹⁰Y sources placed simultaneously using the CSBT device.ResultsThe film-based maximum measured surface dose rate for a bare source in contact with the surface was 3.35 × 10^–7 cGy s⁻¹ Bq⁻¹. When placed in the source holder, the maximum measured dose rate was 1.41 × 10^–7 cGy s⁻¹ Bq⁻¹. The Monte Carlo simulated depth dose rates were within 10% or 0.02 cm of the measured dose rates for each depth of measurement. The maximum film surface dose rate measured using a seven-source configuration within the CSBT device was 1.78 × 10⁻⁷ cGy s⁻¹ Bq⁻¹. Measured and simulated dose rate distribution of the seven-source configuration were compared by gamma analysis and yielded a passing rate of 94.08%. The gamma criteria were 3% for dose-difference and 0.07056 cm for distance-to-agreement. The estimated measured dose rate uncertainty was 5.34%.Conclusions⁹⁰Y is a unique source that can be optimally designed for a customized CSBT device. The rapid dose falloff provided a high dose gradient, ideal for treatment of superficial lesions. The dose rate uncertainty of the ⁹⁰Y-based CSBT device was within acceptable brachytherapy standards and warrants further investigation.     

Impact of lung density on isolated lung tumor dose in VMAT using inline MR-Linac

PurposeThis study investigated the impact of lung density on the isolated lung tumor dose for volumetric modulated arc therapy (VMAT) in an inline magnetic resonance linear accelerator (MR-Linac) using the Monte Carlo (MC) simulation.MethodsCT images of the thorax phantoms with lung tumors of 1, 2, and 3 cm diameters were converted into voxel-base phantoms with lung densities of 0.1, 0.2, and 0.3 g/cm³, respectively. The dose distributions were calculated for partial-arc VMAT. The dose distributions were compared using dose differences, dose volume histograms, and dose volume indices.ResultsIn all cases, the inline magnetic field significantly enhanced the lung tumor dose compared to that at 0 T. For the 1 cm lung tumor, the inline magnetic field of 1 T increased the minimum dose of 95% of the Planning target volume (PTV D₉₅) by 14.0% in 0.1 g/cm³ lung density as compared to that in 0.3 g/cm³ at 0 T. In contrast, at 0 and 0.5 T, the PTV D₉₅ in 0.3 g/cm³ lung density was larger than that in lung density of 0.1 g/cm³. For the 2 cm lung tumor, a similar tendency to 1 cm was observed, whereas the dose impact of lung density was smaller than that for 1 cm. For the 3 cm lung tumor, the lung tumor dose was independent of lung density at 0.5 T and 1.0 T.ConclusionThe inline MR-Linac with the magnetic field over 1 T can enhance the PTV D₉₅ for VMAT regardless of the lung density.     

Validation of 4D Monte Carlo dose calculations using a programmable deformable lung phantom

PurposeTo validate the accuracy of 4D Monte Carlo (4DMC) simulations to calculate dose deliveries to a deforming anatomy in the presence of realistic respiratory motion traces. A previously developed deformable lung phantom comprising an elastic tumor was modified to enable programming of arbitrary motion profiles. 4D simulations of the dose delivered to the phantom were compared with the measurements.MethodsThe deformable lung phantom moving with irregular breathing patterns was irradiated using static and VMAT beam deliveries. Using the RADPOS 4D dosimetry system, point doses were measured inside and outside the tumor. Dose profiles were acquired using films along the motion path of the tumor (S-I). In addition to dose measurements, RADPOS was used to record the motion of the tumor during dose deliveries. Dose measurements were then compared against 4DMC simulations with EGSnrc/4DdefDOSXYZnrc using the recorded tumor motion.ResultsThe agreements between dose profiles from measurements and simulations were determined to be within 2%/2 mm. Point dose agreements were within 2σ of experimental and/or positional/dose reading uncertainties. 4DMC simulations were shown to accurately predict the sensitivity of delivered dose to the starting phase of breathing motions. We have demonstrated that our 4DMC method, combined with RADPOS, can accurately simulate realistic dose deliveries to a deforming anatomy moving with realistic breathing traces. This 4DMC tool has the potential to be used as a quality assurance tool to verify treatments involving respiratory motion. Adaptive treatment delivery is another area that may benefit from the potential of this 4DMC tool.     

Mobile shielding evaluation on the fetal dose during a breast radiotherapy using Monte Carlo simulation

PurposeEvaluation of the out-of-field dose is an important aspect in radiotherapy. Due to the fetus radiosensitivity, this evaluation becomes even more conclusive when the patient is pregnant. In this work, a linear accelerator Varian Clinac 2100c operating at 6 MV, a pregnant anthropomorphic phantom (Maria), and different shields added above the abdominal region of the phantom were used for the analysis based on MCNPX. Methods: The simulations were performed for the medial and lateral projections, using either an open field collimation (10×16 cm²) or a multileaf collimator. The added shields (M1 and M2) were designed based on models proposed by Stovall et al. [1], intending to reduce the deposited dose on the fetus and related structures. Results: The presence of the shields showed to be effective in reducing the doses on the fetus, amniotic sac, and placenta, for example. A reduction of about 43% was found in the dose on the fetus when M2 was added, using the open field collimation, in comparison with the situation with no shield, being the lateral projection the main responsible for the dose. The use of MLC significatively reduced the doses in different structures, including on the fetus and amniotic sac, for example, in comparison to the open field situation. A slight increment on the dose in organs such as the eyes, thyroid and brain was found in both collimation systems, due to the presence of the shields. The contribution of the leakage radiation from the tube head of the linear accelerator was found to be in the order of µGy, being reduced by the presence of the M2 shield. Conclusion: Using the shields showed to be an essential feature in order to reduce the dose not only on the fetus, but also in important structures responsible to its development.     

Dosimetric characterization of the 60Co BEBIG Co0.A86 high dose rate brachytherapy source using PENELOPE

⁶⁰Co sources are being used as an alternative to ¹⁹²Ir sources in high dose rate brachytherapy treatments. In a recent document from AAPM and ESTRO, a consensus dataset for the ⁶⁰Co BEBIG (model Co0.A86) high dose rate source was prepared by using results taken from different publications due to discrepancies observed among them. The aim of the present work is to provide a new calculation of the dosimetric characteristics of that ⁶⁰Co source according to the recommendations of the AAPM and ESTRO report. Radial dose function, anisotropy function, air-kerma strength, dose rate constant and absorbed dose rate in water have been calculated and compared to the results of previous works. Simulations using the two different geometries considered by other authors have been carried out and the effect of the cable density and length has been studied.     

Fetal dose measurements and shielding efficiency assessment in a custom setup of 192Ir brachytherapy for a pregnant woman with breast cancer

PurposeTo assess the radiation dose to the fetus of a pregnant patient undergoing high-dose-rate (HDR) ¹⁹²Ir interstitial breast brachytherapy, and to design a new patient setup and lead shielding technique that minimizes the fetal dose.MethodsRadiochromic films were placed between the slices of an anthropomorphic phantom modeling the patient. The pregnant woman was seated in a chair with the breast over a table and inside a leaded box. Dose variation as a function of distance from the implant volume as well as dose homogeneity within a representative slice of the fetal position was evaluated without and with shielding.ResultsWith shielding, the peripheral dose after a complete treatment ranged from 50 cGy at 5 cm from the caudal edge of the breast to <0.1 cGy at 30 cm. The shielding reduces absorbed dose by a factor of two near the breast and more than an order of magnitude beyond 20 cm. The dose is heterogeneous within a given axial plane, with variations from the central region within 50%. Interstitial HDR ¹⁹²Ir brachytherapy with breast shielding can be more advantageous than external-beam radiotherapy (EBRT) from a radiation protection point of view, as long as the distance to the uterine fundus is higher than about 10 cm. Furthermore, the weight of the shielding here proposed is notably lower than that needed in EBRT.ConclusionsShielded breast brachytherapy may benefit pregnant patients needing localized radiotherapy, especially during the early gestational ages when the fetus is more sensitive to ionizing radiation.     

Monte Carlo investigation on the effect of air gap under bolus in post-mastectomy radiotherapy

ObjectivesTo investigate the dosimetric effect of air gaps under bolus on skin dose for left-sided post-mastectomy radiotherapy with loco regional involvement.MethodsEight patients were planned retrospectively with volume modulated arc therapy (VMAT) and conventional static Field-in-Field (FinF) methods. Three different setups were applied for the 5-mm bolus over the chest wall having 0, 5 or 10 mm air gap under the bolus. The dose calculation was performed using Monte Carlo (MC) simulation. In addition, Analytic Anisotropic Algorithm (AAA) was used to demonstrate the differences observed in clinical setting.ResultsThe investigated air gaps under the bolus had minimal effect on surface dose for FinF plans (relative difference ≤ 2.6%), whereas for VMAT plans the surface dose decreased 13.6% when compared to the case with no air gap. In both FinF and VMAT, the largest differences between AAA and MC were seen at the surface where AAA underestimated the dose by 1.5 Gy (p < 0.05) on average; while the dose in the target volume excluding the surface was relatively similar being on average 0.3 Gy (p > 0.05) larger with AAA than with MC calculations.ConclusionsThe surface dose was significantly lower with VMAT technique than with FinF technique. Possible air gaps under the bolus reduced the surface dose significantly further for VMAT but not for FinF treatments, which may have clinical impact on recurrence rate. AAA was shown to underestimate the surface dose when compared to MC calculation.     

Simulation and measurement of microbeam dose distribution in lung tissue

Microbeam radiation therapy (MRT), a so far preclinical method in radiation oncology, modulates treatment doses on a micrometre scale. MRT uses treatment fields with a few ten micrometre wide high dose regions (peaks) separated by a few hundred micrometre wide low dose regions (valleys) and was shown to spare tissue much more effectively than conventional radiation therapy at similar tumour control rates. While preclinical research focused primarily on tumours of the central nervous system, recently also lung tumours have been suggested as a potential target for MRT.This study investigates the effect of the lung microstructure, comprising air cavities of a few hundred micrometre diameter, on the microbeam dose distribution in lung. In Monte Carlo simulations different models of heterogeneous lung tissue are compared with pure water and homogeneous air–water mixtures. Experimentally, microbeam dose distributions in porous foam material with cavity sizes similar to the size of lung alveoli were measured with film dosimetry at the European Synchrotron Radiation Facility (ESRF) in Grenoble, France.Simulations and experiments show that the microstructure of the lung has a huge impact on the local doses in the microbeam fields. Locally, material inhomogeneities may change the dose by a factor of 1.7, and also average peak and valley doses substantially differ from those in homogeneous material.Our results imply that accurate dose prediction for MRT in lung requires adequate models of the lung microstructure. Even if only average peak and valley doses are of interest, the assumption of a simple homogeneous air–water mixture is not sufficient. Since anatomic information on a micrometre scale are unavailable for clinical treatment planning, alternative methods and models have to be developed.     

Tolerance of the vaginal vault to high-dose rate brachytherapy and concomitant chemo-pelvic irradiation: Long-term perspective

Aim/background

We sought to determine the tolerance level and complication rates of the vaginal vault to combined high-dose-rate intra-cavitary brachytherapy with concomitant chemo-radiotherapy.

Patients and methods

A retrospective review of medical records of all the patients who received definitive chemo-radiotherapy for cervical cancer between 1998 and 2002 was undertaken. The records were reviewed for doses and for radiation-associated early and late sequelae of the vagina, rectum and bladder. Cumulative biological effective dose was calculated for two reference vaginal surface points.

Results

Fifty patients were included. Average age at diagnosis was 54 years. Median follow-up was 59 months. There were no recorded instances of acute grade IV toxicity. Maximal high-dose-rate vaginal surface dose (upper central point) was 103 Gy, and maximal brachytherapy lateral surface dose was 70 Gy. Maximal cumulative biological effective dose for the lateral surface reference point was 465.5 Gy₃, and the maximal cumulative biological effective dose for the superior reference point was 878.6 Gy₃. There were no cases of vaginal necrosis or fistulas, and no cases of grade IV late vaginal, rectal or bladder toxicity. No correlation was found between the maximal vaginal surface dose and vaginal, rectal or bladder toxicity.

Conclusions

The maximal surface HDR brachytherapy dose of 103 Gy and the maximal cBED of 878.6 Gy₃ were not associated with fistula or necrosis or other grade 3–4 vaginal complications. Concomitant chemo-radiotherapy, including pelvic radiotherapy and high-dose-rate intracavitary brachytherapy, is relatively safe for cervical cancer patients.     

Monte Carlo dose calculation of GZP6 60Co stepping source based on a matrix shift technique

Background

As a routine method for stepping source simulation, a Monte Carlo program is run according to the number of steps and then the summation of dose from each run is taken to obtain total dose distribution. This method is time consuming.

Aim

As an alternative method, a matrix shift based technique was applied to simulate a stepping source for brachytherapy.

Materials and methods

The stepping source of GZP6 brachytherapy unit was simulated. In a matrix shift method, it is assumed that a radiation source is stationary and instead the data matrix is shifted based on the number of steps. In this study, by running MCNPX program for one point and calculation of the dose matrix using the matrix shift method, the isodose curves for the esophageal cancer tumor lengths of 4 and 6 cm were obtained and compared with the isodose curves obtained by running MCNPX programs in each step position separately (15 and 23 steps for esophageal cancer tumor lengths of 4 and 6 cm, respectively).

Results

The difference between the two dose matrixes for the stepping and matrix shift methods based on the average dose differences are 3.85 × 10⁻⁴ Gy and 5.19 × 10⁻⁴ Gy for treatment length of 4 cm and 6 cm, respectively. Dose differences are insignificant and these two methods are equally valid.

Conclusions

The matrix shift method presented in this study can be used for calculation of dose distribution for a brachytherapy stepping source as a quicker tool compared to other routine Monte Carlo based methods.     

A Monte Carlo study on dose distribution validation of GZP6 60Co stepping source

Aim

Stepping source in brachytherapy systems is used to treat a target lesion longer than the effective treatment length of the source. Cancerous lesions in the cervix, esophagus and rectum are examples of such a target lesion.

Background

In this study, the stepping source of a GZP6 afterloading intracavitary brachytherapy unit was simulated using Monte Carlo (MC) simulation and the results were used for the validation of the GZP6 treatment planning system (TPS).

Materials and methods

The stepping source was simulated using MCNPX Monte Carlo code. Dose distributions in the longitudinal plane were obtained by using a matrix shift method for esophageal tumor lengths of 8 and 10 cm. A mesh tally has been employed for the absorbed dose calculation in a cylindrical water phantom. A total of 5 × 10⁸ photon histories were scored and the MC statistical error obtained was at the range of 0.008–3.5%, an average of 0.2%.

Results

The acquired MC and TPS isodose curves were compared and it was shown that the dose distributions in the longitudinal plane were relatively coincidental. In the transverse direction, a maximum dose difference of 7% and 5% was observed for tumor lengths of 8 and 10 cm, respectively.

Conclusion

Considering that the certified source activity is given with ±10% uncertainty, the obtained difference is reasonable. It can be concluded that the accuracy of the dose distributions produced by GZP6 TPS for the stepping source is acceptable for its clinical applications.     

Real-time estimation of surface dose based on incident air kerma in diagnostic radiology

PurposeTo estimate the surface dose in diagnostic radiology in real time based on the relationship between the incident air kerma and the surface dose.MethodsThe air kerma for 20 X-ray beams with tube voltages of 50–140 kV and a half-value layer (HVL) of 2.27–9.65 mm Al was measured using an ionization chamber. The beam quality was classified based on the quality indexes (QIs) of 0.4, 0.5, and 0.6, which are defined as the ratio of the effective energy to the maximum energy corresponding to the tube potential. The surface dose for 20 X-ray beams was evaluated based on the measured air kerma, backscatter factor, and ratio of the mass–energy absorption coefficients of water to air, which were calculated using the Monte Carlo method. Finally, the relationship between the air kerma and the surface dose was investigated for X-ray beams with the specific QI values.ResultsThe surface dose at a water phantom was represented by a linear approximation of R² > 0.98, with the air kerma, regardless of the X-ray beam quality. The surface dose estimated based on a linear approximation with the air kerma indicated an agreement within 8% with that evaluated by the chamber measurements at HVL > 3.4 mm Al.ConclusionIt is possible to estimate the surface dose in real time using the linear relationship between the incident air kerma and the surface dose regardless of the X-ray beam quality by accepting ±10% uncertainty in the surface dose estimation.