Monte Carlo calculation of the dose distributions of two 106Ru eye applicators

Monte Carlo simulation of tilted contact plaque brachytherapy placement for juxtapapillary retinoblastoma

Background

The 106-Ruthenium contact plaque applicator is utilized for the treatment of intraocular tumor within a thickness of less than 6 mm. If anything obstructs the placement of the plaque applicator, the treatment is generally difficult because the applicator has to be temporarily located just on the opposite side of the retinal tumor. Furthermore, the plaque applicator edge of approximately 1 mm does not contain ¹⁰⁶Ru, estimating the delivered radiation dose for eccentric tumor is challenging because the lateral dose profile is inadequately provided by the manufacture’s certification. This study aims to simulate tumor coverage of the tilted applicator placement for treating an infant with juxtapapillary retinoblastoma and to achieve the effective treatment.

Case presentation

We present an infant with retinoblastoma whose tumor involved macular and was invading just temporal side of the optic disc. Additionally, posterior staphyloma was induced by a series of previous treatments, making it more difficult to treat the standard plaque placement. Thus, the applicator type of CCA was intentionally tilted to the eyeball and the distance between the posterior edge of the applicator and the eyeball had to be then equal to or more than 2 mm based on the dose distribution of the applicator calculated using Monte Carlo simulation to minimize damage to surrounding tissues while covering the tumor. It was then comparable to the certification and previous reports. Based on the acquired dose distribution, the optimal placement of the applicator was derived from varying the distance between the applicator’s edge and the eyeball, and the distance was then determined to be 2 mm. In this case, the minimum dose rate in the tumor was 25.5 mGy/min, and the time required to deliver the prescribed dose was 26.2 h. Therefore, the tilted ¹⁰⁶Ru plaque applicator placement could deliver the required dose for the treatment. The physical examination revealed no active tumor as a result of the treatment.

Conclusions

Optimizing the placement of the ¹⁰⁶Ru plaque applicator, it was possible to guarantee that the prescribed dose will be delivered to the tumor even if the standard placement is not possible for the juxtapapillary tumor.

Dose characteristics of Au-198 eye brachytherapy applicator: A Monte Carlo study

PURPOSE

The use of ocular plaques is a promising treatment option for eye melanoma brachytherapy. Although several studies have been done on various ocular plaques, little is known about the dose characterization of ¹⁹⁸Au plaque.

MATERIALS AND METHOD

The full mathematical model of the eye phantom, tumor, ¹⁰⁶Ru/¹⁰⁶Rh CCA, and ¹⁹⁸Au plaque were simulated using the Monte Carlo MCNPX code. The dose distribution was measured in the plaque's central axis direction, and a dose profile was also measured at a distance of 2.5 mm from the plaque surface.

RESULTS

The findings showed that ¹⁹⁸Au plaque has superior dosimetric characteristics than CCA plaque for tumors with a thickness of greater than 3.5 mm, while CCA plaque is better for tumors with a thickness of less than 3.5 mm. The dose to the sclera and choroid is higher in the case of CCA plaque, while the dose to the organs at risk (lens and optic nerve) is greater in the case of ¹⁹⁸Au applicator. In the case of ¹⁹⁸Au plaque, however, the dose to sensitive organs was within their permissible dose range.

CONCLUSION

In the treatment of medium and large tumors, ¹⁹⁸Au plaque is more successful than CCA plaque. It can produce a much more homogeneous lateral dose profile in the target. In the treatment of dome-shaped tumors, ¹⁹⁸Au plaque may be more successful than CCA plaque. As a result, the tumor's shape influences the plaque type selection.

Monte Carlo Computation of Dose-Volume Histograms in Structures at Risk of an Eye Irradiated with Heterogeneous Ruthenium-106 Plaques

Background/Aims

The aim of this work is to compare Monte Carlo simulated absorbed dose distributions obtained from ¹⁰⁶Ru eye plaques, whose heterogeneous emitter distribution is known, with the common homogeneous approximation. The effect of these heterogeneities on segmented structures at risk is analyzed using an anthropomorphic phantom.

Methods

The generic CCA and CCB, with a homogeneous emitter map, and the specific CCA1364 and CCB1256 ¹⁰⁶Ru eye plaques are modeled with the Monte Carlo code PENELOPE. To compare the effect of the heterogeneities in the segmented volumes, cumulative dose-volume histograms are calculated for different rotations of the aforementioned plaques.

Results

For the cornea, the CCA with the equatorial placement yields the lowest absorbed dose rate while for the CCA1364 in the same placement the absorbed dose rate is 33% higher. The CCB1256 with the hot spot oriented towards the cornea yields the maximum dose rate per unit of activity while it is 44% lower for the CCB.

Conclusions

Dose calculations based on a homogeneous distribution of the emitter substance yield the lowest absorbed dose in the analyzed structures for all plaque placements. Treatment planning based on such calculations may result in an overdose of the structures at risk.

Pre and post operative radiation protection in Ru-106 brachytherapy ophthalmic plaque surgery and related material shielding properties

Pre and post-operative exposure levels of medical staff and people from public in intra-operative Ru-106 ophthalmic brachytherapy are reported, together with attenuation properties of selected shielding materials. In particular radiation exposure of workers during plaque transportation and during medical assistance of implanted plaque patient was measured. Taking into account dose rates and considering standard assistance procedure of hospitalized patients, the exposure of medical staff and people of the public were evaluated for a given workload. In order to provide tools to optimize radiation protection, considering social and economic aspects due to possible hospital discharge or hospital stay, the attenuation properties of common shielding materials (lead, concrete, red brick, PMMA and gypsum) were measured, considering both narrow and broad beam setups. The eye was simulated using a water equivalent phantom and plaque was fixed on it. All measurements were performed with calibrated survey meters. Results were compared with numerical simulation of bremsstrahlung X-ray radiation spectra emitted from patient eye. Exposure levels measured at 1 m distance in front of the implanted eye are 0.05 µSv/h/MBq, at 10 cm from patient head, 0.44 µSv/h/MBq (plaque side), 0.4 µSv/h/MBq (front), 0.25 µSv/h/MBq (lateral, opposed to plaque), 0.2 µSv/h/MBq (back). Average exposure levels, under conservative assumptions, for medical staff is 17 µSv/patient and less than 23 µSv/patient for careers and comforters. TVLs in lead and concrete are about 1.6 cm and 11.5 cm respectively.

Monte Carlo Simulation of the Treatment of Uveal Melanoma Using Measured Heterogeneous 106Ru Plaques

Background/Aims

Ruthenium plaques are used for the treatment of ocular tumors. The aim of this work is the comparison between simulated absorbed dose distributions tallied in an anthropomorphic phantom, obtained from ideal homogeneous plaques, and real eye plaques in which the actual heterogeneous distribution of ¹⁰⁶Ru was measured. The placement of the plaques with respect to the tumor location was taken into consideration to optimize the effectiveness of the treatment.

Methods

The generic CCA and CCB, and the specific CCA1364 and CCB1256 ¹⁰⁶Ru eye plaques were modeled with the Monte Carlo code PENELOPE. To compare the suitability of each treatment for an anterior, equatorial and posterior tumor location, cumulative dose-volume histograms for the tumors and structures at risk were calculated.

Results

Eccentric placements of the plaques, taking into account the inhomogeneities of the emitter map, can substantially reduce the dose delivered to structures at risk while maintaining the prescribed dose at the tumor apex.

Conclusions

The emitter map distribution of the plaque and the computerized tomography of the patient used in a Monte Carlo simulation allow an accurate determination of the plaque position with respect to the tumor with the potential to reduce the dose to sensitive structures.

Monte Carlo Estimation of Absorbed Dose Distributions Obtained from Heterogeneous 106Ru Eye Plaques

BACKGROUND

The distribution of the emitter substance in ¹⁰⁶Ru eye plaques is usually assumed to be homogeneous for treatment planning purposes. However, this distribution is never homogeneous, and it widely differs from plaque to plaque due to manufacturing factors.

METHODS

By Monte Carlo simulation of radiation transport, we study the absorbed dose distribution obtained from the specific CCA1364 and CCB1256 ¹⁰⁶Ru plaques, whose actual emitter distributions were measured. The idealized, homogeneous CCA and CCB plaques are also simulated.

RESULTS

The largest discrepancy in depth dose distribution observed between the heterogeneous and the homogeneous plaques was 7.9 and 23.7% for the CCA and CCB plaques, respectively. In terms of isodose lines, the line referring to 100% of the reference dose penetrates 0.2 and 1.8 mm deeper in the case of heterogeneous CCA and CCB plaques, respectively, with respect to the homogeneous counterpart.

CONCLUSIONS

The observed differences in absorbed dose distributions obtained from heterogeneous and homogeneous plaques are clinically irrelevant if the plaques are used with a lateral safety margin of at least 2 mm. However, these differences may be relevant if the plaques are used in eccentric positioning.

Gold nanoparticle-based brachytherapy enhancement in choroidal melanoma using a full Monte Carlo model of the human eye

The effects of gold nanoparticles (GNPs) in 125I brachytherapy dose enhancement on choroidal melanoma are examined using the Monte Carlo simulation technique. Usually, Monte Carlo ophthalmic brachytherapy dosimetry is performed in a water phantom. However, here, the compositions of human eye have been considered instead of water. Both human eye and water phantoms have been simulated with MCNP5 code. These simulations were performed for a fully loaded 16 mm COMS eye plaque containing 13 125I seeds. The dose delivered to the tumor and normal tissues have been calculated in both phantoms with and without GNPs. Normally, the radiation therapy of cancer patients is designed to deliver a required dose to the tumor while sparing the surrounding normal tissues. However, as the normal and cancerous cells absorbed dose in an almost identical fashion, the normal tissue absorbed radiation dose during the treatment time. The use of GNPs in combination with radiotherapy in the treatment of tumor decreases the absorbed dose by normal tissues. The results indicate that the dose to the tumor in an eyeball implanted with COMS plaque increases with increasing GNPs concentration inside the target. Therefore, the required irradiation time for the tumors in the eye is decreased by adding the GNPs prior to treatment. As a result, the dose to normal tissues decreases when the irradiation time is reduced. Furthermore, a comparison between the simulated data in an eye phantom made of water and eye phantom made of human eye composition, in the presence of GNPs, shows the significance of utilizing the composition of eye in ophthalmic brachytherapy dosimetry Also, defining the eye composition instead of water leads to more accurate calculations of GNPs radiation effects in ophthalmic brachytherapy dosimetry.

Monte Carlo Simulation of the Treatment of Eye Tumors with (106)Ru Plaques: A Study on Maximum Tumor Height and Eccentric Placement

BACKGROUND/AIMS

Ruthenium plaques are used for the treatment of ocular tumors. There is, however, a controversy regarding the maximum treatable tumor height. Some advocate eccentric plaque placement, without a posterior safety margin, to avoid collateral damage to the fovea and optic disc, but this has raised concerns about marginal tumor recurrence. There is a need for quantitative information on the spatial absorbed dose distribution in the tumor and adjacent tissues. We have overcome this obstacle using an approach based on Monte Carlo simulation of radiation transport.

METHODS

CCA and CCB (106)Ru plaques were modeled and their geometry embedded in a computerized tomography scan of the eye of a patient. Different tumor sizes and locations were simulated with the general-purpose Monte Carlo code PENELOPE.

RESULTS

Cumulative dose-volume histograms were obtained for the tumors and the tissues at risk considered. Plots of isodose lines for both plaques were obtained in a computerized tomography study.

CONCLUSIONS

Ruthenium eye plaques are an adequate treatment option for tumors up to around 5 mm in height. According to our results, assuming a correct placement of the plaque, a tumor of 6.5 mm apical height is about the maximum size that can be treated safely with the large CCB plaque.

A new human eye model for ophthalmic brachytherapy dosimetry

The present work proposes a new mathematical eye model for ophthalmic brachytherapy dosimetry. This new model includes detailed description of internal structures that were not treated in previous works, allowing dose determination in different regions of the eye for a more adequate clinical analysis. Dose calculations were determined with the MCNP-4C Monte Carlo particle transport code running n parallel environment using PVM. The Amersham CKA4 ophthalmic applicator has been chosen and the depth dose distribution has been determined and compared to those provide by the manufacturer. The results have shown excellent agreement. Besides, absorbed dose values due to both 125I seeds and 60Co plaques were obtained for each one of the different structures which compose the eye model and can give relevant information in eventual clinical analyses.

Dosimetry characterization of 32P intravascular brachytherapy source wires using Monte Carlo codesPENELOPEandGEANT4

Monte Carlo calculations using the codes PENELOPE and GEANT4 have been performed to characterize the dosimetric parameters of the new 20 mm long catheter-based 32P beta source manufactured by the Guidant Corporation. The dose distribution along the transverse axis and the two-dimensional dose rate table have been calculated. Also, the dose rate at the reference point, the radial dose function, and the anisotropy function were evaluated according to the adapted TG-60 formalism for cylindrical sources. PENELOPE and GEANT4 codes were first verified against previous results corresponding to the old 27 mm Guidant 32P beta source. The dose rate at the reference point for the unsheathed 27 mm source in water was calculated to be 0.215 +/- 0.001 cGy s(-1) mCi(-1), for PENELOPE, and 0.2312 +/- 0.0008 cGy s(-1) mCi(-1), for GEANT4. For the unsheathed 20 mm source, these values were 0.2908 +/- 0.0009 cGy s(-1) mCi(-1) and 0.311 0.001 cGy s(-1) mCi(-1), respectively. Also, a comparison with the limited data available on this new source is shown. We found non-negligible differences between the results obtained with PENELOPE and GEANT4.

Dosimetry for radiocolloid therapy of cystic craniopharyngiomas

The dosimetry for radiocolloid therapy of cystic craniopharyngiomas is investigated. Analytical calculations based on the Loevinger and the Berger formulas for electrons and photons, respectively, are compared with Monte Carlo simulations. The role of the material of which the colloid introduced inside the craniopharyngioma is made of as well as that forming the cyst wall is analyzed. It is found that the analytical approaches provide a very good description of the simulated data in the conditions where they can be applied (i.e., in the case of a uniform and infinite homogeneous medium). However, the consideration of the different materials and interfaces produces a strong reduction of the dose delivered to the cyst wall in relation to that predicted by the Loevinger and the Berger formulas.

Characterization of a high-dose-rate 90Sr-90Y source for intravascular brachytherapy by using the Monte Carlo code PENELOPE

Radiation treatment with catheter-based beta-emitter sources is currently under clinical trial to prevent restenosis. In the present paper, we address the characterization of the high-dose-rate 90Sr-90Y seeds of the Beta-Cath system supplied by Novoste Corporation, one of the commercially available sources for intravascular brachytherapy. The Monte Carlo code PENELOPE has been used to simulate the transport of electrons emitted by the encapsulated 90Sr-90Y seeds. The calculated radial dose function and anisotropy function for a single seed in water are compared with simulation results from other authors. Regarding g(r), the present result lies between the ITS3 and EGS4 curves, being somewhat closer to ITS3, while in the case of F(r, theta) some differences appear for certain angular intervals and radial distances. In order to put the observed differences into perspective, we have calculated radial doses for point isotropic sources in water. Our results for 0.5 and 1 MeV electrons are in good agreement with simulations using EGSnrc, and an excellent agreement is obtained with ITS for point 90Sr-90Y emitters. Dose distributions in water are calculated for source 'trains' consisting of 1, 2, 3, 4, 5, 9 and 12 seeds. The dose at the source midplane is enhanced if the number of seeds is up to 4, and saturates for trains with 5 or more seeds. We also compare the dose distribution obtained by simply adding the contributions of individual seeds with the simulation of the complete source train. It is found that both calculation procedures yield essentially the same result for distances greater than about 2 mm. Finally, the contribution of bremsstrahlung photons to the dose is briefly analysed.