|
 
 |
| ORIGINAL ARTICLE |
|
| Year : 2012 | Volume
: 2
| Issue : 3 | Page : 144-148 |
|
Calculation of excess dose to the eye phantom due to a distanced shielding for electron therapy in head and neck cancers
Keyvan Jabbari1, Mahnaz Roayaei2, Hosein Saberi3
1 Department of Medical Physics and Engineering, School of Medicine; Medical Image and Signal Processing Research Center, Isfahan, Iran
2 Department of Radiation Oncology, Isfahan University of Medical Sciences, Isfahan, Iran
3 Department of Medical Physics and Engineering, School of Medicine, Isfahan, Iran
| Date of Submission | 06-Jun-2012 |
| Date of Acceptance | 19-Jul-2012 |
| Date of Web Publication | 20-Sep-2019 |
Correspondence Address:
Keyvan Jabbari
Department of Medical Physics and Engineering, School of Medicine, Isfahan University of Medical Sciences, Isfahan
Iran
 Source of Support: None, Conflict of Interest: None  | 2 |
DOI: 10.4103/2228-7477.112004
For superficial lesions, the electrons may be used for radiation therapy. The high energy photons and electrons are produced by a Linear accelerator (Linac). Many of electron fields need the shielding of normal or critical organs. The electron shields are usually lead slabs with few millimeter thicknesses which should be placed near the skin, less than 1 cm away from skin. In the inspection of patients setting in a clinic by a physicist, it was noted that in some cases the technician places the shields far away from skin in the way that the shadow of the field still matches the shielded area. This is due to a conceptual mistake in which one assumes that electrons travel in a straight line and matching the shadow of lead slab is enough for the shielding. This project is about Monte Carlo simulation of this case and dosimetry in which the excess dose to the tissue under the shield is calculated. In this study, BEAMnrc and DOSXYZnrc are used for simulation of the Linac and the electron shields. The water phantom as well as the Linac head (NEPTON Linac) is simulated in the electron mode. The simulation is performed in three various cases in which the lead shield is placed in distances of 1, 20, 40 cm from the surface of the phantom. In all cases, the edge of the shield is matched with the light field, so the shields get smaller as they move from the surface because of the divergence of the light field. The simulations were done in two energies, 6 and 13 MeV. The experiments also were done with EDR2 film dosimetry and the simulation results were validated using the experimental results. In all cases, the dose under the shield was normalized to the dose in the center. The dose of the normal organ under the shield was 2, 38, 43% with respect to the center for shield distances of 1, 20, 40 cm, respectively. So there is a considerable increasing of the dose due to the distanced shielding. In this work exact amount of the dose from this mistake (distanced shielding) is calculated and simulated. Keywords: Electron therapy, Monte Carlo, shielding
How to cite this article:
Jabbari K, Roayaei M, Saberi H. Calculation of excess dose to the eye phantom due to a distanced shielding for electron therapy in head and neck cancers. J Med Signals Sens 2012;2:144-8 |
How to cite this URL:
Jabbari K, Roayaei M, Saberi H. Calculation of excess dose to the eye phantom due to a distanced shielding for electron therapy in head and neck cancers. J Med Signals Sens [serial online] 2012 [cited 2023 Jun 1];2:144-8. Available from: https://www.jmssjournal.net/text.asp?2012/2/3/144/112004 |
| Introduction | |  |
Radiation therapy techniques are typically planned to conform to the volume of tissue needing irradiation. Therefore, shielding of the portion of the body outside the treatment field is essentially carried out by the machine collimators that define the field edges. After the surgery, radiation therapy prescription with incomplete margin results to improved local control and survival. [1],[2],[3],[4]
Shielding of the normal tissue is one of the important aspects of radiation therapy. Therapeutic procedures use higher energies in the megavoltage (MV) range. At these energies, the radiation is more penetrating. For high energy photons, to provide a 3% transmission, a thickness of approximately 6.5 cm of lead would be required for a photon beam with a typical half-value layer of 1.3 cm. [4]
Non-melanoma skin cancer is the most common type of the cancer despite of the low morbidity. Radiation therapy is one of the main options for treatment of this kind of cancer which has an equivalent outcome to surgery. For superficial lesions, the electrons may be used in for radiation therapy. [5],[6] Many of electron fields need the shielding of normal or critical organs. The electron shields are usually lead slabs with few millimeters thickness which should be place near the skin, less than 1 cm away from skin. The margin of the electron field should be between 1 and 2.5 cm for various sizes of tumors. For example, the usage of the small field sizes around the eye cause the inhomogeneity of the dose. In electron fields, three points are important: immobilization of the patient, perpendicular direction of the beam to the surface, and the shielding. [7],[8],[9],[10],[11],[12],[13] Therefore, protection of the normal tissue around the tumor is essential and the shield should not transmit more than 5% of the radiation. The shielding of critical organs such as eyes is more important. The internal shielding is inconvenient for the patient, and usually the external shield is used. On the other hand, the external shielding has few advantages such as the coverage of the lachrymal gland [14] and the ciliary protection.
In the inspection of patients setting in a clinic by authors (physicist and doctors), it was noticed that in some cases the technician places the shields far away from skin in the way that shadow of field still matches the shielded area [Figure 1]. This is due to a conceptual mistake in which one assumes that electrons travel in a straight line and matching the shadow of lead slab is enough for the shielding. This project is about Monte Carlo simulation of this case in which the excess dose to the tissue under the shield is calculated. The film dosimetry of this situation is also performed. | Figure 1: The lead shield used in electron mode, which is placed in various distances from the surface
Click here to view |
| Materials and Methods | | |
Unlike photons, electrons travel through a curve even in short distances. Therefore, the adjustment of the lead shields using the light field could cause an error. The light field is designed in a way that represents the straight line from the source to the surface. This design is appropriate for photon shielding since photons move through a straight line between each interaction. When an electron shield is placed far away from the surface, the electrons can reach the shielded organ due to a lateral scattering. It should be noted that in all cases the shield is adjusted to the light field so it looks like a right set up from the light field point of the view.
The Monte Carlo method is used to calculate the exact amount of the excess dose versus the distance of the shield from the surfaces. Monte Carlo technique is a statistical method based on random numbers. In this method, many particles with random manner are simulated to calculate the parameter of interest. After simulation of very large number of particles (photons or electrons depending on Linac mode) very accurate results can be produced. [15],[16],[17]
For simulation, the head of the Linac is simulated and the main components of the Linac head should be considered [Figure 2]. The primary electron beam first enters the head and passes the exit window. The primary electron beam is a narrow high energy electron pencil beam. This narrow beam cannot be used for treatment and it has to be spread over a wide area. For this reason there is a scattering foil after the collimator. The scattering foil is a lead slab with 0.5 mm thickness. | Figure 2: Simulation of the Linac's head using BEAMnrc. In electron mode, external applicators are attached to the head. The handmade shield should be placed after applicators, close to the surface
Click here to view |
In the path of the electron beam, we also have the ionization chambers which monitor the amount of the Linac output. The ionization chamber has three layers of aluminum, air and Kapton. At the end, we have secondary collimators; however, in electron mode, the applicator is attached to the head. The applicators extend to vicinity of the surface and generate sharp rectangular field sizes. There are many other components in the head such as mirror and light system which are not placed in the path of the electrons during the treatment and has no effect on the radiation dose.
In this study, BEAMnrc and DOSXYZnrc are used for simulation of the geometry. BEAMnrc is designed to simulate various linear accelerators. [18] The code is easy to learn the accuracy of this code has been evaluated in many different situations. [19],[20],[21],[22],[23],[24],[25],[26],[27],[28],[29],[30] The water phantom as well as the Linac head (NEPTON Linac) is simulated in the electron mode. The lead shields are designed in various distances from the surface. At the long distances from the surface, the shields are placed between the applicator slabs. This is exactly what happened in real clinical cases in which, for ease of work, the technicians placed the shields on the applicators instead of the surface of the patient. The shape of the shields in the simulation was similar to a real clinical case, which is a simple rectangular 5 × 5 cm 2 field with a 2 × 2 cm 2 shield at one corner of the field. The simulation is performed in 3 various cases in which the lead shield is placed in distances of 1, 20, 40 cm from the surface of the phantom. In all cases, the edge of the shield is matched with the light field, so the shields get smaller as they move from the surface because of the divergence of the light field. The simulations were done in two energies, 6 and 13 MeV.
The output file of the BEAMnrc includes the entire characteristics of the exited electrons which is available after the run. This file is called Phase-space file and it can be used as an input file for DOSXYZnrc for dose calculation in the phantom. [31],[32] All these Monte Carlo codes are available for free in NRC website (National Research Council of Canada). The DOSEXYZ code using exit Phase-space file, transports all of the exited particles from the Linac head in the phantom. The output file of the code includes the absorbed dose in each voxel of the phantom. This output file is imported to in-house MATLAB software and the absorbed dose and isodose line can be derived for final evaluation.
The experiments were done with EDR2 film dosimetry and the simulation results were validated using the experimental results. EDR2 films are one of best means for two-dimensional high energy dosimetry [33],[34] The calibration curve of this film is linear in a wide dose range, between 5 and 600 cGy. [35] For experiment the films were placed in the surface of a PMMA phantom under the Linac. The size of the phantom was 30 cm in each dimension. The applicators were also attached to the gantry since the Linac is used in the electron mode. The electron shields were made from lead sheets with 3 mm thickness. For the shield placed 1 cm from the surface the shape of the field was rectangle with size of 5×5 cm 2 which had a 2×2 cm 2 shield at one corner. For larger distances of the shield from the surface which is related to the wrong set-up, the general shape of the shield is the same but the sizes are reduced because of the divergence of the light field. The dose delivery is calculated in the way that 100 cGy is delivered to the center of the field. Before the experiments the calibration of the film was performed considering the recommendations of Zeu et al. and Childress et al. and the calibration curve was derived for various dose values. [36],[37],[38]
| Results and Discussion | | |
For simulation of the BEAMnrc, the run were done with minimum 40 million histories. The number of the history determines the number of the primary electrons that starts from the beginning of the Linac head. As mentioned before, the output of the BEAMnrc software is called phase-space file which is in a special format and includes all information of the particles in the related plane. In one step for evaluation of the results, the phase-space files are converted to the intensity map of the particles using in-house MATLAB software. MATLAB with its powerful image processing toolbox is the appropriate software for evaluation of the dosimetry results. [39],[40],[41],[42],[43] The MATLAB software using BEAMdp extracts the spatial distribution of the particles in the phase-space file which is placed at the end of the Linac head and applicators [Figure 3]. It should be noted that small rectangular shield in the right corner of the field resemble the eye shield in the real clinical case. | Figure 3: Intensity map generated from BEAMnrc code for the case in which the shield is placed 1, 20, and 40 cm from the surface. (5 × 5 cm2, 6 MeV). For large distances, the increasing of the dose outside the field is clearly visible
Click here to view |
The results of film dosimetry are illustrated in [Figure 4] in which d represents the distance of the shield from the surface. The Linac is in electron mode and the energy of the electrons is 6 MeV. As illustrated, for d = 1 cm, there is clear sharp penumbra and the area under the shield is protected. For the case of d = 20 cm the dose spreads out into the shielded region. This effect gets worse with increasing the distance of the shield from the surface. | Figure 4: Film dosimetry of the electron field in which the shield is placed 1 cm and 20 cm from the surface (6 MeV and 13 MeV electrons). For the case in which the shield is placed 20 cm from the surface, the light field still matches the shielded area
Click here to view |
The lateral scattering of the electron decreases for high energies since the electrons tend to move in forward direction. For this reason the amount of the excess dose decreases versus energy increasing, however we can still see the excess dose in the shielded region. This is illustrated in [Figure 4] for 13 MeV electrons. In the left figure, the shield is placed 1 cm away from the surface, in which we have again a radiation field with sharp edges. In the right, the distance of the shield is 20 cm from the surface and there is still considerable dose to the area under the shield.
For film development and processing, it is recommended that irradiated films should not developed sooner than 1 h from the time of irradiation. This point is required to get stable results for dosimetry. Therefore, all the films were kept in the cover and they were processed and developed 1 day after the experiment. It is also recommended that the development situation would be similar to the situation at the time of the film calibration. Therefore, the same equipment in one radiology center is used for processing and development of the films for film calibration as well as dosimetry experiments.
In all cases the dose under the shield were normalized to the dose in the center. The dose of the normal organ under the shield was 2%, 38%, 43% with respect to the center for shield distances of 1, 20, 40 cm, respectively. So there is a considerable increase of the dose due to distanced shielding.
[INLINE:1][INLINE:2]
The electron shielding has a different concept compare to photon shields. The electrons, because of the multiple scattering and continues lost of the energy, move through a curve. Even in the air, the path of the electron cannot be assumed as straight line. This is a known fact that the electron should be placed close to the surface. In this project the amount of the excess dose is calculated due to the distanced shielding. This shielding is due to a conceptual mistake of technician who applies the electron shield similar to photon shield far away from the skin. In this work exact amount of the dose from this mistake is calculated and simulated. For photon mode, the shields are usually placed in a tray which is more than 20 cm away from the surface however that electron shields should not be placed more than 1 cm from the skin of the patient.
| References | | |
| 1. | Perez CA. Management of incompletely excised carcinoma of the skin. Int J Radiat Oncol Biol Phys 1991;20:903-4.  |
| 2. | Wong JR, Wang CC. Radiation therapy in the management of cutaneous malignancies. Clin Dermatol 2001;19:348-53. |
| 3. | Petsuksiri J, Frank SJ, Garden AS, Ang KK, Morrison WH, Chao KS, et al. Outcomes after radiotherapy for squamous cell carcinoma of the eyelid. Cancer 2008;112:111-8. |
| 4. | Kwan W, Wilson D, Moravan V. Radiotherapy for locally advanced basal cell and squamous cell carcinomas of the skin. Int J Radiat Oncol Biol Phys 2004;60:406-11. |
| 5. | Tsao MN, Tsang RW, Liu FF, Panzarella T, Rotstein L. Radiotherapy management for squamous cell carcinoma of the nasal skin: The Princess Margaret Hospital experience. Int J Radiat Oncol Biol Phys 2002;52:973-9. |
| 6. | Amdur RJ, Kalbaugh KJ, Ewald LM, Parsons JT, Mendenhall WM, Bova FJ, et al. Radiation therapy for skin cancer near the eye: Kilovoltage x-rays versus electrons. Int J Radiat Oncol Biol Phys 1992;23:769-79. |
| 7. | Kwok SK, Ho PC, Leung SF, Gandhi S, Lee VW, Lam DS, et al. An analysis of the incidence and risk factors of developing severe keratopathy in eyes after megavoltage external beam irradiation. Ophthalmology 1998;105:2051-5. |
| 8. | Meli JA, Vitali PE. A simple eye-lens shielding technique compatible with independent jaws of a linear accelerator. Med Dosim 1992;17:221-3. |
| 9. | Mukundan S Jr, Wang PI, Frush DP, Yoshizumi T, Marcus J, Kloeblen E, et al. MOSFET dosimetry for radiation dose assessment of bismuth shielding of the eye in children. AJR Am J Roentgenol 2007;188:1648-50. |
| 10. | Butson MJ, Cheung T, Yu PK, Price S, Bailey M. Measurement of radiotherapy superficial X-ray dose under eye shields with radiochromic film. Phys Med 2008;24:29-33. |
| 11. | Dirican B, Oysul K, Beyzadeoglu M, Surenkok S, Pak Y. Experimental lens-sparing optimization in therapeutic orbital irradiation with electron beams. Neoplasma 2004;51:390-4. |
| 12. | Reft C, Rash C, Dabrowski J, Roeske JC, Hallahan D. Eye shielding for patients treated with total body irradiation. Med Dosim 1996;21:73-8. |
| 13. | Rustgi SN. Dose distribution under external eye shields for high energy electrons. Int J Radiat Oncol Biol Phys 1986;12:141-4. |
| 14. | Baker CR, Luhana F, Thomas SJ. Absorbed dose behind eye shields during kilovoltage photon radiotherapy. Br J Radiol 2002;75:685-8. |
| 15. | Chetty IJ, Curran B, Cygler JE, DeMarco JJ, Ezzell G, Faddegon BA, et al. Report of the AAPM Task Group No. 105: Issues associated with clinical implementation of Monte Carlo-based photon and electron external beam treatment planning. Med Phys 2007;34:4818-53. |
| 16. | Verhaegen F, Seuntjens J. Monte Carlo modelling of external radiotherapy photon beams. Phys Med Biol 2003;48:R107-64. |
| 17. | Kawrakow I. The effect of Monte Carlo statistical uncertainties on the evaluation of dose distributions in radiation treatment planning. Phys Med Biol 2004;49:1549-56. |
| 18. | Downes P, Jarvis R, Radu E, Kawrakow I, Spezi E. Monte Carlo simulation and patient dosimetry for a kilovoltage cone-beam CT unit. Med Phys 2009;36:4156-67. |
| 19. | Rogers DW, Faddegon BA, Ding GX, Ma CM, We J, Mackie TR. BEAM: A Monte Carlo code to simulate radiotherapy treatment units. Med Phys 1995;22:503-24. |
| 20. | Taylor RE, Rogers DW. EGSnrc Monte Carlo calculated dosimetry parameters for 192Ir and 169Yb brachytherapy sources. Med Phys 2008;35:4933-44. |
| 21. | Taylor RE, Rogers DW. An EGSnrc Monte Carlo-calculated database of TG-43 parameters. Med Phys 2008;35:4228-41. |
| 22. | Faddegon BA, Kawrakow I, Kubyshin Y, Perl J, Sempau J, Urban L. The accuracy of EGSnrc, Geant4 and PENELOPE Monte Carlo systems for the simulation of electron scatter in external beam radiotherapy. Phys Med Biol 2009;54:6151-63. |
| 23. | Hasenbalg F, Fix MK, Born EJ, Mini R, Kawrakow I. VMC++ versus BEAMnrc: A comparison of simulated linear accelerator heads for photon beams. Med Phys 2008;35:1521-31. |
| 24. | Jabbari K, Sarfehnia A, Podgorsak EB, Seuntjens JP. Monte Carlo feasibility study of orthogonal bremsstrahlung beams for improved radiation therapy imaging. Phys Med Biol 2007;52:1171-84. |
| 25. | Rogers DW, Kawrakow I. Monte Carlo calculated correction factors for primary standards of air kerma. Med Phys 2003;30:521-32. |
| 26. | Rogers DW. Fifty years of Monte Carlo simulations for medical physics. Phys Med Biol 2006;51:R287-301. |
| 27. | Ma CM, Faddegon BA, Rogers DW, Mackie TR. Accurate characterization of Monte Carlo calculated electron beams for radiotherapy. Med Phys 1997;24:401-16. |
| 28. | Ma CM, Jiang SB. Monte Carlo modeling of electron beams from medical accelerators. Phys Med Biol 1999;44:R157-89. |
| 29. | Ma CM, Mok E, Kapur A, Pawlicki T, Findley D, Brain S, et al. Clinical implementation of a Monte Carlo treatment planning system. Med Phys 1999;26:2133-43. |
| 30. | Sheikh-Bagheri D, Rogers DW, Ross CK, Seuntjens JP. Comparison of measured and Monte Carlo calculated dose distributions from the NRC linac. Med Phys 2000;27:2256-66. |
| 31. | Mainegra-Hing E, Reynaert N, Kawrakow I. Novel approach for the Monte Carlo calculation of free-air chamber correction factors. Med Phys 2008;35:3650-60. |
| 32. | Sarfehnia A, Jabbari K, Seuntjens J, Podgorsak EB. Experimental verification of beam quality in high-contrast imaging with orthogonal bremsstrahlung photon beams. Med Phys 2007;34:2896-906. |
| 33. | Ju SG, Han Y, Kum O, Cheong KH, Shin EH, Shin JS, et al. Comparison of film dosimetry techniques used for quality assurance of intensity modulated radiation therapy. Med Phys 2010;37:2925-33. |
| 34. | Sathiyan S, Ravikumar M, Varatharaj C, Supe SS, Keshava SL. IMRT Implementation and patient specific dose verification with film and ion chamber array detectors. Gulf J Oncolog 2010;1:20-7. |
| 35. | Novotny J Jr, Bhatnagar JP, Quader MA, Bednarz G, Lunsford LD, Huq MS. Measurement of relative output factors for the 8 and 4 mm collimators of Leksell Gamma Knife Perfexion by film dosimetry. Med Phys 2009;36:1768-74. |
| 36. | Zhu XR, Jursinic PA, Grimm DF, Lopez F, Rownd JJ, Gillin MT. Evaluation of Kodak EDR2 film for dose verification of intensity modulated radiation therapy delivered by a static multileaf collimator. Med Phys 2002;29:1687-92. |
| 37. | Childress NL, Salehpour M, Dong L, Bloch C, White RA, Rosen II. Dosimetric accuracy of Kodak EDR2 film for IMRT verifications. Med Phys 2005;32:539-48. |
| 38. | Childress NL, Rosen II. Effect of processing time delay on the dose response of Kodak EDR2 film. Med Phys 2004;31:2284-8. |
| 39. | Jabbari K, Keall P, Seuntjens J. Considerations and limitations of fast Monte Carlo electron transport in radiation therapy based on precalculated data. Med Phys 2009;36:530-40. |
| 40. | Abella V, Miro R, Juste B, Verdu G. Comparison between voxelized, volumized and analytical phantoms applied to radiotherapy simulation with Monte Carlo. Conf Proc IEEE Eng Med Biol Soc 2009;2009:1258-61. |
| 41. | Chow JC, Leung MK, Islam MK, Norrlinger BD, Jaffray DA. Evaluation of the effect of patient dose from cone beam computed tomography on prostate IMRT using Monte Carlo simulation. Med Phys 2008;35:52-60. |
| 42. | Mukumoto N, Tsujii K, Saito S, Yasunaga M, Takegawa H, Yamamoto T, et al. A preliminary study of in-house Monte Carlo simulations: An integrated Monte Carlo verification system. Int J Radiat Oncol Biol Phys 2009;75:571-9. |
| 43. | Georg D, Nyholm T, Olofsson J, Kjaer-Kristoffersen F, Schnekenburger B, Winkler P, et al. Clinical evaluation of monitor unit software and the application of action levels. Radiother Oncol 2007;85:306-15. |
| Authors | | |
 Keyvan Jabbari received B.S. degree in physics from Isfahan University of Technology, Isfahan, IRAN in 1999, M.S. degree in medical radiation engineering from Amir Kabir University, Tehran in 2002, M.S. degree in medical physics from University of Manitoba , Winnipeg, Canada, in 2004. Then he got Ph.D. of Medical Physics from McGill University, Montreal, Canada in 2008.
Since then he is assistant professor in Department of Medical Physics and Medical Engineering , School of Medicine, Isfahan University of Medical Sciences. He is also clinical medical physicist in Isfahan Milad Hospital and Seyed-o-Shohada cancer center. His research interest is fast Monte Carlo in treatment planning, dose calculation and plan optimization.
 Mahnaz Roayaei received specialist M.D. degree from Isfahan University of Medical Sciences, Isfahan, IRAN in 2004. She is radiation therapist in Seyedo Shohada hospital. She is currently assistant professor in School of Medicine, Isfahan University of Medical Sciences. Her research interest is Treatment planning for radiation therapy.
 Hosein Saberi received B.Sc degree in physics from Qom University, Qom, Iran (2005), M.Sc degree in medical physics from Tehran University of Medical Science, Tehran, Iran (2008). Currently he is doing PhD in Isfahan University of Medical Sciences, Isfahan, Iran. His research interest is Monte Carlo in treatment planning and Medical Image processing.
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
|
Search Pubmed for