ORIGINAL ARTICLE Year : 2015  Volume : 40  Issue : 1  Page : 1317 Monte Carlobased dose calculation for ^{32} P patch source for superficial brachytherapy applications Sridhar Sahoo^{1}, Selvam T Palani^{1}, SK Saxena^{1}, D A R Babu^{2}, A Dash^{2}, ^{1} Radiological Physics and Advisory Division, Bhabha Atomic Research Centre, Mumbai, Maharashtra, India ^{2} Isotope Production and Applications Division, Bhabha Atomic Research Centre, Mumbai, Maharashtra, India Correspondence Address: Skin cancer treatment involving ^{32} P source is an easy, less expensive method of treatment limited to small and superficial lesions of approximately 1 mm deep. Bhabha Atomic Research Centre (BARC) has indigenously developed ^{32} P nafionbased patch source (1 cm Χ 1 cm) for treating skin cancer. For this source, the values of dose per unit activity at different depths including dose profiles in water are calculated using the EGSnrcbased Monte Carlo code system. For an initial activity of 1 Bq distributed in 1 cm ^{2} surface area of the source, the calculated central axis depth dose values are 3.62 Χ 10 ^{10} GyBq ^{1} and 8.41 Χ 10 ^{11} GyBq ^{1} at 0.0125 and 1 mm depths in water, respectively. Hence, the treatment time calculated for delivering therapeutic dose of 30 Gy at 1 mm depth along the central axis of the source involving 37 MBq activity is about 2.7 hrs.
Introduction Basal cell carcinoma is one of the most common skin cancers, occurs mostly in middle aged people, and is more probable for the fair complexion people. [1] The treatment modalities for skin cancers are surgical excision, radiotherapy and chemotherapy. Each treatment modality has its own advantages and disadvantages. Removing the affected area by surgical excision is usually preferred in many cases, but the recurrence rates after treatment are high. Radiotherapy treatment using external beam therapy is too expensive and it also delivers unnecessary dose to underlying normal tissues. Chemotherapy has its own side effects. Mould or superficial brachytherapy is a promising alternative treatment method for such skin cancers, where highenergy beta emitting radionuclides such as 32 P, 90 Sr/ 90 Y, 188 Re are used to overcome the disadvantages of radiotherapy and surgery. In superficial brachytherapy, prescribed dose can be delivered to the affected area without excessive damage to the neighboring normal tissues. This technique is simple, less trauma to patients, and less expensive as compared to external beam therapy. Lee et al., introduced the treatment of skin cancer and Bowen's disease using beta emitting 165 Hoimpregnated patch sources. [2] Successful tumor destruction was observed both in animal and human studies. Mukherjee et al., in their studies evaluated 90 Y skin patches and 188 Re radioactive bandages for therapy of superficial tumors in mice. [3],[4] Treatment of skin cancer using 188 Relabeled paper patches has been reported by Jeong et al. [5] Pandey et al., reported the use of 32 P cellulosebased adsorbent paper skin patches to control the tumor regression in C57BL6 mice bearing melanoma. [6] Park et al., [7] studied the use of 32 P ophthalmic applicator after pterygium and glaucoma surgeries. They demonstrated that dose distributions obtained using the 32 P source is beneficial for reducing the incidence rate of radiationinduced cataract and it can deliver therapeutic doses to the surface of the conjunctiva while sparing the lens better than the 90 Sr/ 90 Y applicators. Xu et al., [8] investigated the therapeutic effects of the chromic phosphate particlebased 32 P source in a rabbit VX 2 lung tumor animal model and found that the tumor volume significantly decreased after implantation of source particle. Salguerio et al., designed 32 P brachytherapy patch source (1 mm in height × 5 mm in dia.) for skin diseases using phosphoric acid and chromic phosphate in combination with natural rubber or silicone and evaluated its therapeutic efficacy. [9],[10] They reported arrest of tumor growth and complete regression of tumor in some cases with 40 Gy of singledose scheme in animal studies. They estimated the dose rate at selected depths (0.0001, 0.01, 4 and 7.5 mm) using the Monte Carlobased MCNP5 code. [10],[11] The activity per unit area considered in their calculations was 10.6 MBqcm 2 . The surface area of the source was 0.196 cm 2 . Hence, the total activity of the source considered in their work was 2.081 MBq. We repeated their study using the DOSRZnrc usercode. [12] The dose rate values showed a good agreement for 0.0001 and 0.01 mm depths. For 4 and 7.5 mm depths, the published values were higher by a factor of about 22 and 3.6 × 10 4 , respectively. We concluded that this large discrepancy in the dose rate values at 4 and 7.5 mm depths published by Salguerio et al., [9] was due to possible systematic error in their Monte Carlo calculations. [13] 32 P is a suitable radioisotope for such therapeutic application due to many advantages over other beta emitting radioisotopes. It is a pure beta emitter with maximum energy of 1.71 MeV. Its halflife is 14.2 days. Hence it is less hazardous material from transportation, storage, and waste disposal point of view. The maximum range of 32 P beta particle in soft tissue is 8 mm (the average range is 3 mm). [14] Due to its short range, there will be negligible radiation dose to the underlying healthy normal tissues and bone. Isotope Production and Applications Division, BARC has indigenously developed nafionzirconium phosphate filmbased 32 P patch source for superficial brachytherapy applications. [15] A nafion117 membrane of thickness 100 μm is treated with ZrOCl 2 solution, and subsequently dipped in orthophosphoric acid. These radioactive 32 P patches are cut in to 1 cm × 1 cm sizes and then subsequently laminated with thermoplastic polyurethane sheets of thickness 40 μm. The above preparation method is robust, inexpensive and reproducible and complies with the safety standard stipulated by Atomic Energy Regulatory Board, India. [16] The detailed preparation of nafion117 patches is explained by Saxena et al. [15] The present study is aimed at calculating central axis depth dose and dose profiles in water phantom for the indigenously developed 32 Pnafionbased patch source. For this purpose, the EGSnrcbased Monte Carlo code system is used. [17] Based on the calculated dose rate data, the treatment time to deliver a therapeutic dose of 30 Gy at reference depth is also calculated, as per the IAEAtecdoc1274. [18] Materials and Methods Monte Carlo calculations DOSXYZnrc usercode [19] of the EGSnrccode system [17] is used to calculate central axis depth doses and dose profiles in the unit density water medium for simulating the 1 cm × 1 cm 32 Pnafionpatch source. The 1 cm × 1 cm 32 Pnafionpatch source is positioned on 2 × 2 ×2 cm 3 water phantom. The thickness of source is 100 μm. The geometry and coordinate system used in the Monte Carlo calculations is shown in [Figure 1]. The elemental composition and density of phosphorousloaded zirconiumnafion117 composite membrane used in the Monte Carlo calculation is given in [Table 1]. [15]{Figure 1}{Table 1} The 32 P beta spectrum [see [Figure 2]] needed for the Monte Carlo calculation is based on ICRU Report No. 56. [20] In the Monte Carlo calculations, it is considered that the source particles are uniformly distributed in the nafion patch of dimensions 1 cm × 1 cm × 100 μm. The water phantom was divided in to voxels of dimension of 0.25 × 0.25 × 0.25 mm 3 for generating dose profiles. Dose distributions in water are scored in these voxels. Separate simulation is carried out to score central axis depth dose by using bigger voxel dimensions (2 × 2× 0.25 mm 3 ).{Figure 2} In the Monte Carlo calculations, all secondary particles such as knock on electrons and secondary bremsstrahlung photons produced by primary source electrons are completely followed in the simulation. The PEGS4 dataset needed for EGSnrc calculations is based on the XCOM compilations. [21] The PEGS4 dataset is generated by setting AE = ECUT = 0.521 MeV and AP = PCUT = 0.01 MeV, where the parameters AE and AP are the low energy thresholds for the production of knockon electrons and secondary bremsstrahlung photons, respectively. All Monte Carlo simulations utilized the PRESTAII algorithm. The electron stepsize parameter is set at ESTEP = 0.25. To increase the speed of the calculations, electron range rejection technique was used by setting ESAVE = 2 MeV. Auxiliary simulation was also carried without range rejection. The calculations suggest that using the range rejection with ESAVE = 2 MeV, improves the efficiency of the calculations by about 40%. In order to know the effect of boundary crossing algorithm on computational time, both PRESTAI and EXACT boundary crossing algorithms were used in the calculations. The study shows that the using the PRESTAI option results in improving the efficiency of the calculations by a factor 2 when compared to using the EXACT boundarycrossing algorithm. This observation is consistent with the findings by Walters and Kawrakow in their EGSnrcbased Monte Carlo study involving radiotherapy electron beams. [22] Results and Discussion The variation of the dose values per unit activity (GyBq 1 ) as a function of depth (mm) in water for the 32 Pnafionpatch source is shown in [Figure 3]. The dose decreases rapidly with increasing depth in water. Central axis dose at 4 mm depth in water is only 0.08% of the central axis surface dose. Such a rapid decrease in dose will result in better sparing of the normal tissues.{Figure 3} [Table 2] compares the values of central axis depth dose per unit activity (GyBq 1 ) of 32 Pnafionpatch source with the corresponding values of 32 Psiliconpatch [9] for different depths in water. Higher dose rate values are observed in the case of 32 Psiliconpatch source, because the radioactivity is distributed in lesser surface area (0.196 cm 2 ) as compared to 32 Pnafionpatch source, where surface area is 1 cm 2 .{Table 2} For treatment time calculation, 1 mm depth from the surface along the central axis of the source is considered as reference depth. [18] The value of dose in water calculated at 1 mm from the source surface is 8.41 × 10 11 GyBq 1 . Hence, the time required to deliver a therapeutic dose of 30 Gy for a 37 MBq of radioactivity distributed in 1 cm 2 of 32 Pnafionpatch source is about 2.7 hours. [Figure 4] presents the dose rate profiles along the xaxis of 32 Pnafionpatch source for three different depths z = 0.5, 1 and 2 mm from the source surface. [Figure 5] presents normalized dose values along the xaxis at depth 1 mm. The central axis dose value at 1 mm depth is used for normalization. Dose rate value at 3.5 mm away from the central axis is about 91% of the central axis value. Whereas dose rate at 5 mm away from the central axis is only 50% of the central axis value. [Figure 6] [Figure 7] [Figure 8] show isodose profiles of the 32 Pnafionpatch source at depths of 0.5, 1, and 2 mm. About 3.253.5 mm distance around the central axis is covered by about 90% isodose line for depths of 0.5, 1, and 2 mm. Hence, the 32 Pnafionpatch source is effective for treatment of approximately 6.57.0 mm diameter lesions.{Figure 4}{Figure 5}{Figure 6}{Figure 7}{Figure 8} Conclusions Dose distributions for the indigenously developed 1 cm × 1 cm 32 Pnafion skin patch source are calculated using the Monte Carlobased EGSnrc code system. The calculated treatment time for delivering therapeutic dose of 30 Gy at 1 mm depth along the central axis of the source involving 37 MBq activity is about 2.7 hrs. This source is effective for treatment of approximately 6.57.0 mm diameter lesions. Acknowledgments The authors would like to thank Dr. D. N. Sharma, Formerly Director, Health, Safety and Environment Group, Bhabha Atomic Research Centre for his encouragement and support for this work. References


