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.     

Application of Monte Carlo calculations for validation of a treatment planning system in high dose rate brachytherapy

AimThe accuracy of treatment planning systems is of vital importance in treatment outcomes in brachytherapy. In the current study the accuracy of dose calculations of a high dose rate (HDR) brachytherapy treatment planning system (TPS) was validated using the Monte Carlo method.Materials and methodsThree ⁶⁰Co sources of the GZP6 afterloading brachytherapy system were modelled using MCNP4C Monte Carlo (MC) code. The dose distribution around all the sources was calculated by MC and a dedicated treatment planning system. The results of both methods were compared.ResultsThere was good agreement (<2%) between TPS and MC calculated dose distributions except at a point near the sources (<1 cm) and beyond the tip of the sources.ConclusionsOur study confirmed the accuracy of TPS calculated dose distributions for clinical use in HDR brachytherapy.     

Calculation of dose distribution near an innovative concentric balloon catheter for endovascular brachytherapy

Purpose: Using a radioactive solution-filled catheter for intravascular irradiation has the potential problem of chemical and radiological toxicity in the case of a balloon rupture. In order to reduce this risk, an innovative concentric balloon catheter was developed.Methods and Materials: The concentric balloon was made by inner and outer balloons filled with saline and radioactive solution, respectively. The optimal inner radius was determined by comparing the dose rate reduction vs. the volume reduction for various inner and outer radii for 188Re, 32P, and 90Y solutions.Results: For a balloon with an outer radius of 1.5 mm, there was no advantage of a concentric balloon. For balloons with outer radii of 3.0 and 5 mm, the optimal inner radius was 1.5 and 3 mm, respectively.Conclusions: With the newly designed concentric balloon, the risk of toxicity can be reduced while keeping the dose rate high enough so that the treatment times within tolerable limits are still maintained.     

Radiation dose around a PET scanner installation: Comparison of Monte Carlo simulations,analytical calculations and experimental results

PurposeMonte Carlo study of radiation transmission around areas surrounding a PET room.MethodsAn extended population of patients administered with ¹⁸F-FDG for PET-CT investigations was studied, collecting air kerma rate and gamma ray spectra measurements at a reference distance. An MC model of the diagnostic room was developed, including the scanner and walls with variable material and thickness. MC simulations were carried out with the widely used code GEANT4.ResultsThe model was validated by comparing simulated radiation dose values and gamma ray spectra produced by a volumetric source with experimental measurements; ambient doses in the surrounding areas were assessed for different combinations of wall materials and shielding and compared with analytical calculations, based on the AAPM Report 108.In the range 1.5–3.0 times of the product between the linear attenuation coefficient and thickness of an absorber (μ x), it was observed that the effectiveness of different combinations of shielding is roughly equivalent. An extensive tabulation of results is given in the text.ConclusionsThe validation tests performed showed a satisfactory agreement between the simulated and expected results. The simulated dose rates incident on, and transmitted by the walls in our model of PET scanner room, are generally in good agreement with analytical estimates performed using the AAPM Publication No. 108 method. This provides an independent confirmation of AAPM's approach. Even in this specific field of application, GEANT4 proved to be a relevant and accurate tool for dosimetry estimates, shielding evaluation and for general radiation protection use.     

A Monte Carlo study on dose distribution evaluation of Flexisource 192Ir brachytherapy source

Aim

The aim of this study is to evaluate the dose distribution of the Flexisource ¹⁹²Ir source.

Background

Dosimetric evaluation of brachytherapy sources is recommended by task group number 43 (TG. 43) of American Association of Physicists in Medicine (AAPM).

Materials and methods

MCNPX code was used to simulate Flexisource ¹⁹²Ir source. Dose rate constant and radial dose function were obtained for water and soft tissue phantoms and compared with previous data on this source. Furthermore, dose rate along the transverse axis was obtained by simulation of the Flexisource and a point source and the obtained data were compared with those from Flexiplan treatment planning system (TPS).

Results

The values of dose rate constant obtained for water and soft tissue phantoms were equal to 1.108 and 1.106, respectively. The values of the radial dose function are listed in the form of tabulated data. The values of dose rate (cGy/s) obtained are shown in the form of tabulated data and figures. The maximum difference between TPS and Monte Carlo (MC) dose rate values was 11% in a water phantom at 6.0 cm from the source.

Conclusion

Based on dosimetric parameter comparisons with values previously published, the accuracy of our simulation of Flexisource ¹⁹²Ir was verified. The results of dose rate constant and radial dose function in water and soft tissue phantoms were the same for Flexisource and point sources. For Flexisource ¹⁹²Ir source, the results of TPS calculations in a water phantom were in agreement with the simulations within the calculation uncertainties. Furthermore, the results from the TPS calculation for Flexisource and MC calculation for a point source were practically equal within the calculation uncertainties.     

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.     

A Monte Carlo approach to calculate dose distribution around the lineal brachytherapy sources.

Monte Carlo means the statistical methods used to solve stochastic or deterministic physical or mathematical problems. The use of this technique for solving deterministic problems is indicated whenever there is no analytical solution for the problem or when this solution is very difficult or when the use of numerical methods requires an excessive amount of CPU-time. This paper describes the application of Monte Carlo methods using a computer program, called ISODOSE, which calculates isodose curves around linear radioactive sources to be used in brachytherapy treatment-planning. Brachytherapy is a special form of cancer treatment in which the radioactive source is placed near or inside the tumor, in order to produce cancer tissue necrosis. The program is written in C language, and the results will be compared to similar data obtained from a commercially available program at Hospital do Cancer de Recife (Brazil), radiation therapy institute. The use of Monte Carlo techniques in the ISODOSE program allows for plotting isodose curves around linear sources, and it is especially more precise near the borders of the source (singularities), taking less time than it would be possible by using deterministic methods. The ISODOSE program can be used for brachytherapy planning in small clinics in Recife and adjacent areas.     

Comparison between beta radiation dose distribution due to LDR and HDR ocular brachytherapy applicators using GATE Monte Carlo platform

Eye applicators with 90Sr/90Y and 106Ru/106Rh beta-ray sources are generally used in brachytherapy for the treatment of eye diseases as uveal melanoma. Whenever, radiation is used in treatment, dosimetry is essential. However, knowledge of the exact dose distribution is a critical decision-making to the outcome of the treatment. The Monte Carlo technique provides a powerful tool for calculation of the dose and dose distributions which helps to predict and determine the doses from different shapes of various types of eye applicators more accurately. The aim of this work consisted in using the Monte Carlo GATE platform to calculate the 3D dose distribution on a mathematical model of the human eye according to international recommendations. Mathematical models were developed for four ophthalmic applicators, two HDR 90Sr applicators SIA.20 and SIA.6, and two LDR 106Ru applicators, a concave CCB model and a flat CCB model. In present work, considering a heterogeneous eye phantom and the chosen tumor, obtained results with the use of GATE for mean doses distributions in a phantom and according to international recommendations show a discrepancy with respect to those specified by the manufacturers. The QC of dosimetric parameters shows that contrarily to the other applicators, the SIA.20 applicator is consistent with recommendations. The GATE platform show that the SIA.20 applicator present better results, namely the dose delivered to critical structures were lower compared to those obtained for the other applicators, and the SIA.6 applicator, simulated with MCNPX generates higher lens doses than those generated by GATE.     

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.     

Measurement-based model of a wide-bore CT scanner for Monte Carlo dosimetric calculations with GMCTdospp software

The aim of this work was to create a model of a wide-bore Siemens Somatom Sensation Open CT scanner for use with GMCTdospp, which is an EGSnrc-based software tool dedicated for Monte Carlo calculations of dose in CT examinations.The method was based on matching spectrum and filtration to half value layer and dose profile, and thus was similar to the method of Turner et al. (Med. Phys. 36, pp. 2154–2164). Input data on unfiltered beam spectra were taken from two sources: the TASMIP model and IPEM Report 78. Two sources of HVL data were also used, namely measurements and documentation. Dose profile along the fan-beam was measured with Gafchromic RTQA-1010 (QA+) film. Two-component model of filtration was assumed: bow-tie filter made of aluminum with 0.5 mm thickness on central axis, and flat filter made of one of four materials: aluminum, graphite, lead, or titanium.Good agreement between calculations and measurements was obtained for models based on the measured values of HVL. Doses calculated with GMCTdospp differed from the doses measured with pencil ion chamber placed in PMMA phantom by less than 5%, and root mean square difference for four tube potentials and three positions in the phantom did not exceed 2.5%. The differences for models based on HVL values from documentation exceeded 10%. Models based on TASMIP spectra and IPEM78 spectra performed equally well.     

Dual-energy contrast-enhanced digital mammography: Glandular dose estimation using a Monte Carlo code and voxel phantom

PurposeTo estimate the mean glandular dose of contrast enhanced digital mammography, using the EGSnrc Monte Carlo code and female adult voxel phantom.MethodsAutomatic exposure control of full field digital mammography system was used for the selection of the X-ray spectrum and the exposure settings for dual energy imaging. Measurements of the air-kerma and of the half value layers were performed and a Monte Carlo simulation of the digital mammography system was used to compute the mean glandular dose, for breast phantoms of various thicknesses, glandularities and for different X-ray spectra (low and high energy).ResultsFor breast phantoms of 2.0–8.0 cm thick and 0.1–100% glandular fraction, CC view acquisition, from AEC settings, can result in a mean glandular dose of 0.450 ± 0.022 mGy −2.575 ± 0.033 mGy for low energy images and 0.061 ± 0.021 mGy – 0.232 ± 0.033 mGy for high energy images. In MLO view acquisition mean glandular dose values ranged between 0.488 ± 0.007 mGy – 2.080 ± 0.021 mGy for low energy images and 0.065 ± 0.012 mGy – 0.215 ± 0.010 mGy for high energy images.ConclusionThe low kV part of contrast enhanced digital mammography is the main contributor to total mean glandular breast dose. The results of this study can be used to provide an estimated mean glandular dose for individual cases.     

FRoG dose computation meets Monte Carlo accuracy for proton therapy dose calculation in lung

PurposeTo benchmark and evaluate the clinical viability of novel analytical GPU-accelerated and CPU-based Monte Carlo (MC) dose-engines for spot-scanning intensity-modulated-proton-therapy (IMPT) towards the improvement of lung cancer treatment.MethodsNine patient cases were collected from the CNAO clinical experience and The Cancer Imaging Archive-4D-Lung-Database for in-silico study. All plans were optimized with 2 orthogonal beams in RayStation (RS) v.8. Forward calculations were performed with FRoG, an independent dose calculation system using a fast robust approach to the pencil beam algorithm (PBA), RS-MC (CPU for v.8) and general-purpose MC (gp-MC). Dosimetric benchmarks were acquired via irradiation of a lung-like phantom and ionization chambers for both a single-field-uniform-dose (SFUD) and IMPT plans. Dose-volume-histograms, dose-difference and γ-analyses were conducted.ResultsWith respect to reference gp-MC, the average dose to the GTV was 1.8% and 2.3% larger for FRoG and the RS-MC treatment planning system (TPS). FRoG and RS-MC showed a local γ-passing rate of ~96% and ~93%. Phantom measurements confirmed FRoG’s high accuracy with a deviation < 0.1%.ConclusionsDose calculation performance using the GPU-accelerated analytical PBA, MC-TPS and gp-MC code were well within clinical tolerances. FRoG predictions were in good agreement with both the full gp-MC and experimental data for proton beams optimized for thoracic dose calculations. GPU-accelerated dose-engines like FRoG may alleviate current issues related to deficiencies in current commercial analytical proton beam models. The novel approach to the PBA implemented in FRoG is suitable for either clinical TPS or as an auxiliary dose-engine to support clinical activity for lung patients.     

Monte Carlo dose verification of VMAT treatment plans using Elekta Agility 160-leaf MLC

In this study, we verified volumetric modulated arc therapy (VMAT) plans in an Elekta Synergy system with an integrated Agility 160-leaf multileaf collimator (MLC) by comparing them with Monte Carlo (MC)-calculated dose distributions using the AAPM TG-119 structure sets. The head configuration of the linear accelerator with the integrated MLC was simulated with the EGSnrc/BEAMnrc code. Firstly, the dosimetric properties of the MLC were evaluated with the MC technique and film measurements. Next, VMAT plans were created with the Pinnacle³ treatment planning system (TPS) for four regions in the AAPM TG-119 structures. They were then verified by comparing them with MC-calculated dose distributions using dose volume histograms (DVHs) and three-dimensional (3D) gamma analysis. The MC simulations for the Agility MLC dosimetric properties were in acceptable agreement with measurements. TPS-VMAT plans using TG-119 structure sets agreed with MC dose distributions within 2% in the comparison of D₉₅ in planning target volumes (PTVs) evaluated from DVHs. In contrast, higher dose regions such as D₂₀, D₁₀, and D₅ in PTVs for TPS tended to be smaller than MC values. This tendency was particularly noticeable for mock head and neck with complicated structures. In 3D gamma analysis, the passing rates with 3%/3mm criteria in PTVs were ≥99%, except for mock head and neck (89.5%). All passing rates for organs at risk (OARs) were in acceptable agreement of >96%. It is useful to verify dose distributions of PTVs and OARs in TPS-VMAT plans by using MC dose calculations and 3D gamma analysis.     

Radial secondary electron dose profiles and biological effects in light-ion beams based on analytical and Monte Carlo calculations using distorted wave cross sections

To speed up dose calculation, an analytical pencil-beam method has been developed to calculate the mean radial dose distributions due to secondary electrons that are set in motion by light ions in water. For comparison, radial dose profiles calculated using a Monte Carlo technique have also been determined. An accurate comparison of the resulting radial dose profiles of the Bragg peak for (1)H(+), (4)He(2+) and (6)Li(3+) ions has been performed. The double differential cross sections for secondary electron production were calculated using the continuous distorted wave-eikonal initial state method (CDW-EIS). For the secondary electrons that are generated, the radial dose distribution for the analytical case is based on the generalized Gaussian pencil-beam method and the central axis depth-dose distributions are calculated using the Monte Carlo code PENELOPE. In the Monte Carlo case, the PENELOPE code was used to calculate the whole radial dose profile based on CDW data. The present pencil-beam and Monte Carlo calculations agree well at all radii. A radial dose profile that is shallower at small radii and steeper at large radii than the conventional 1/r(2) is clearly seen with both the Monte Carlo and pencil-beam methods. As expected, since the projectile velocities are the same, the dose profiles of Bragg-peak ions of 0.5 MeV (1)H(+), 2 MeV (4)He(2+) and 3 MeV (6)Li(3+) are almost the same, with about 30% more delta electrons in the sub keV range from (4)He(2+)and (6)Li(3+) compared to (1)H(+). A similar behavior is also seen for 1 MeV (1)H(+), 4 MeV (4)He(2+) and 6 MeV (6)Li(3+), all classically expected to have the same secondary electron cross sections. The results are promising and indicate a fast and accurate way of calculating the mean radial dose profile.     

The FLUKA Monte Carlo code coupled with an OER model for biologically weighted dose calculations in proton therapy of hypoxic tumors

IntroductionThe increased radioresistance of hypoxic cells compared to well-oxygenated cells is quantified by the oxygen enhancement ratio (OER). In this study we created a FLUKA Monte Carlo based tool for inclusion of both OER and relative biological effectiveness (RBE) in biologically weighted dose (ROWD) calculations in proton therapy and applied this to explore the impact of hypoxia.MethodsThe RBE-weighted dose was adapted for hypoxia by making RBE model parameters dependent on the OER, in addition to the linear energy transfer (LET). The OER depends on the partial oxygen pressure (pO₂) and LET. To demonstrate model performance, calculations were done with spread-out Bragg peaks (SOBP) in water phantoms with pO₂ ranging from strongly hypoxic to normoxic (0.01–30 mmHg) and with a head and neck cancer proton plan optimized with an RBE of 1.1 and pO₂ estimated voxel-by-voxel using [¹⁸F]-EF5 PET. An RBE of 1.1 and the Rørvik RBE model were used for the ROWD calculations.ResultsThe SOBP in water had decreasing ROWD with decreasing pO₂. In the plans accounting for oxygenation, the median target doses were approximately a factor 1.1 lower than the corresponding plans which did not consider the OER. Hypoxia adapted target ROWDs were considerably more heterogeneous than the RBE_1.1-weighted doses.ConclusionWe realized a Monte Carlo based tool for calculating the ROWD. Read-in of patient pO₂ and estimation of ROWD with flexibility in choice of RBE model was achieved, giving a tool that may be useful in future clinical applications of hypoxia-guided particle therapy.     

Affordable compensator for IMRT delivery designed for low-and-middle income countries,a Monte Carlo study

We propose a novel cost-effective compensator that can be used to facilitate access to IMRT in low-and-middle income countries. The compensator has the advantages of simplicity, less downtime, increased reliability and less impact on treatment quality from patient motion during treatment. Moreover, the system can be used with either a cobalt-60 unit or linear accelerator. In this Monte Carlo study, the dosimetric properties of the new compensator design have been evaluated. Results were obtained for different field sizes of cobalt teletherapy machine, and the dose was scored at 0.5 cm depth in a water phantom. The effects of compensator thickness, filling material type and shape, and field size were identified. Furthermore, the percentage depth dose and beam profiles for various field sizes and at different depths were obtained. Beam profiles show no significant signature of the beads relative to a solid compensator; in addition, they exhibit a better flatness while preserving symmetry for all field sizes. A reusable bead-based compensator appears to be feasible, and provides dose distribution similar to a solid compensator with low cost and no hazards. Our results avail the ongoing efforts to expand the reach to IMRT in low- and middle-income countries.