Monte Carlo dosimetry of a novel Yttrium-90 disc source for episcleral brachytherapy

Abstract

Purpose: To calculate the dose distribution using Monte Carlo simulations for a novel high-dose-rate Yttrium-90 (Y-90) disc source recently developed for episcleral brachytherapy and provide a lookup table for treatment planning. Methods: Monte Carlo simulations were performed to calculate the in-water dose distribution of the Y-90 disc source using the “GATE”, a software based on the “Geant4” Monte Carlo simulation toolkit developed by the international OpenGATE collaboration. The geometry of this novel beta source, its capsule, and the surrounding water medium were accurately modeled in the simulation input files. The standard Y-90 element beta spectrum from ICRU 72 was used, and the physics processes for beta and photon interactions with matters were all included. The dose distribution of this Y-90 disc source was measured in a separate study using Gafchromic EBT-3 films and the results were reported elsewhere. To match the setup of the experiment, a Gafchromic EBT-3 film was also included in the simulation geometry. The simulated dose profiles were exported from the 3D dose distribution results and compared with the measured dose profiles. Transverse dose profiles at different distances from the seed surface were also obtained to study the lateral coverage of the source. Results: The measured percent depth dose (PDD) curves along the central axis perpendicular to the surface of the Y-90 disc were constructed from the experimental and simulated data, and normalized to the reference point at 1 mm from the source capsule. Both PDD curves agreed well up to 4 mm from the source surface (maximum difference ± 10%) but deviated from each other beyond 4 mm. The deviation might be caused by the increased measurement uncertainty in the low-dose region. The dose rate at the reference point calculated from the Monte Carlo simulation was 1.09 cGy/mCi-s and agreed very well with the measured dose rate of 1.05 cGy/mCi-s. If the 80% isodose line is selected as the lateral coverage, the lateral dose coverage is maximal (∼4.5 mm) at the plane next to the source surface, and gradually decreases with the increasing distance, approaching 3.5 mm when the plane is 5 mm from the 6-mm diameter source surface. Conclusion: Monte Carlo simulations were successfully performed to confirm the measured PDD curve of the novel Y-90 disc source. This simulation work laid a solid foundation for characterizing the full dosimetry parameters of this source for episcleral brachytherapy applications.