|Year : 2010 | Volume
| Issue : 1 | Page : 23-32
Monte Carlo simulation of a multi-leaf collimator design for telecobalt machine using BEAMnrc code
Komanduri M Ayyangar1, M Dinesh Kumar2, Pradush Narayan1, Fenedit Jesuraj1, MR Raju1
1 International Cancer Center, Mahatma Gandhi Memorial Medical Trust Hospital, Bhimavaram, AP-534 204, India
2 Yashoda Hospital, Secunderabad, India
|Date of Submission||26-May-2009|
|Date of Decision||14-Jul-2009|
|Date of Acceptance||28-Aug-2009|
|Date of Web Publication||7-Jan-2010|
Komanduri M Ayyangar
Mahatma Gandhi Memorial Medical Trust Hospital, Peda Amiram, Bhimavaram, AP-534 204
Source of Support: None, Conflict of Interest: None
| Abstract|| |
This investigation aims to design a practical multi-leaf collimator (MLC) system for the cobalt teletherapy machine and check its radiation properties using the Monte Carlo (MC) method. The cobalt machine was modeled using the BEAMnrc Omega-Beam MC system, which could be freely downloaded from the website of the National Research Council (NRC), Canada. Comparison with standard depth dose data tables and the theoretically modeled beam showed good agreement within 2%. An MLC design with low melting point alloy (LMPA) was tested for leakage properties of leaves. The LMPA leaves with a width of 7 mm and height of 6 cm, with tongue and groove of size 2 mm wide by 4 cm height, produced only 4% extra leakage compared to 10 cm height tungsten leaves. With finite 60 Co source size, the interleaf leakage was insignificant. This analysis helped to design a prototype MLC as an accessory mount on a cobalt machine. The complete details of the simulation process and analysis of results are discussed.
Keywords: BEAMnrc, Cobalt-60, multi-leaf collimator, Monte Carlo
|How to cite this article:|
Ayyangar KM, Kumar M D, Narayan P, Jesuraj F, Raju M R. Monte Carlo simulation of a multi-leaf collimator design for telecobalt machine using BEAMnrc code. J Med Phys 2010;35:23-32
|How to cite this URL:|
Ayyangar KM, Kumar M D, Narayan P, Jesuraj F, Raju M R. Monte Carlo simulation of a multi-leaf collimator design for telecobalt machine using BEAMnrc code. J Med Phys [serial online] 2010 [cited 2019 Jul 21];35:23-32. Available from: http://www.jmp.org.in/text.asp?2010/35/1/23/58780
| Introduction|| |
Clinical radiotherapy using 60 Co machines are a cost-effective alternative to low energy linear accelerators.  Addition of MLC in these machines will make treatment more versatile for implementing 3D CRT plans. There are nearly 300 cobalt machines in India and many more in the world. Most of these machines are similar in design to Theratron-80 models (original design by Atomic Energy Canada Limited (AECL)).
Monte Carlo (MC) simulation can be used to test beams from radiotherapy machines of various designs in terms of their adequacy of shielding and beam transport. For a diverging type collimator Han et al,  simulated a cobalt beam from Theratron-780C model teletherapy machine using MC computations with Stanford Electron Gamma Shower (EGS) code. They used an original fluence of 2 million photons generated with 2 cm diameter teletherapy source, with depleted uranium collimator. The MC method allows for realistic and accurate results.  Instead of actually doing the experiment with expensive equipment, a virtual machine can be simulated and virtual experiments can be performed using MC methods provided the problem is defined accurately and accurate data are used.
MC is vital for development in radiation oncology physics. MC calculations have been used to generate basic data for external beam treatment planning. Finite size pencil beam point spread functions derived from MC data were used in convolution algorithms and analytical anisotropic algorithms. In brachytherapy, basic TG43 data were generated by MC and verified by experiment. In this work, we applied MC methodology to validate the efficacy of a prototype design of MLC for a telecobalt machine.
The cobalt-60 teletherapy machine has been modeled in the past by various investigators. ,,,,,, Complete cobalt source was simulated by Mora et al,  using BEAM code as early as in 1999. Earlier, in 1988, Rogers  used EGS to calculate the cobalt source in a limited study. The source size used by Rogers was 2 cm diameter. Subsequently, smaller cobalt sources with higher specific activity were available. In 2004, Al-Basheer  has conducted a detailed study of the properties of Theratron telecobalt machine using MCNP code. However, he used a source size of 1.5 cm diameter and 3 cm height. More recently there were other studies ,, of simulating Theratron sources using MC modeling.
In this paper, we present the building and testing of cobalt-60 teletherapy machine for virtual simulation purposes, calculating depth dose and cross profiles, and generating the leakage data for a typical MLC design.
| Materials and Methods|| |
The current simulation was performed for a Phoenix model cobalt machine manufactured by Kirloskar Technologies, Harihar, India. This is an isocentric machine with source to axis (SAD) distance of 80 cm and source diaphragm distance (SDD) of 45 cm. In addition, a block tray holder with SDD of 55 cm is available. This machine is of similar design to a Theratron -80 60 Co machine (M/s Theratronics, Canada).
EGS and OMEGA BEAM Software
The MC code used in this paper is known as the OMEGA BEAM (Ottawa-Madison Electron Gamma Algorithm) which was originally developed collaboratively by NRC Canada and University of Wisconsin through a grant awarded by NIH in 1994.  The BEAM code was specifically developed for radiation therapy beams. It simulates - targets, flattening filters, scattering foils, mirrors, jaws, applicators, etc. The original BEAM code has undergone several revisions and is now called BEAMnrc. The version used in this paper was revised as of February 2007. The BEAM code requires the EGS software as prerequisite, is now called EGSnrcMP, where MP signifies multi-platform. The EGSnrcMP is an extended and improved version of the earlier EGS4 package and incorporates many improvements in the implementation of the condensed history technique for the simulation of charged particle transport and latest low energy cross sections. (See EGSnrcMP link below).
Free versions of the software are available from http://www.irs.inms.nrc.ca/EGSnrc/EGSnrc.html (EGSnrcMP) and http://www.irs.inms.nrc.ca/BEAM/beamhome.html none (BEAMnrc).
Documentation was available in PDF format and was a part of the installation package. Installation packages contained the following modules:
60 Co Source Configuration and Beam Transport
- EGSnrcMP - the software that tracks the electron-gamma transport
- BEAMnrc - builds the treatment machines
- BEAMdp - analyses phase space files created by BEAMnrc
- DOSXYZnrc - a transport code to process phase space files and compute dose distributions in CT matrix or phantoms.
The source size used in the current simulation was 1.5 cm diameter and 2 cm height. The active source was surrounded by steel cladding and a lead shield except for an opening for the beam to exit. Figure 1a shows a diagram of the source housing used for simulation in this work.
All the simulations were carried out on a Dell PC with 2 GHz speed and 512 MB RAM. The current simulation was a three-step process. In the first step, simulation was carried out only for the source, source housing and primary collimator. Results of the simulation were stored as a phase space file at the end of the primary collimator. The primary collimator opening was enough to accommodate the largest field size at the isocenter. The phase space file was scored over a diameter of 10 cm and a total number of 1.7 billion histories were collected in about 160 hours of computational time. This phase space file typically was of size 1.0 Gb. [Figure 1]b shows the cobalt source housing and primary collimator modeled in step I of the MC simulation. A detailed description of the input file used for this step is given in appendix.
In step II, simulation was performed using the phase space file generated in step I through the inner and outer jaws and to the top of water phantom placed at 80 cm from source plane. Step II had to be repeated for each new field size. Typically, for a 10 cm x 10 cm field simulation, five billion histories are processed and computational time is around 40 hours. This computation resulted in generation of a phase space file on the top of the water phantom. This file typically was of size 1.2 Gb. In addition, during this process, BEAMnrc calculated the central axis depth dose data in the phantom. [Figure 2] shows the complete simulation geometry used in step II with true distances and dimensions of the inner collimator on the Phoenix cobalt machine.
The transport code automatically calculated the percent standard deviation of the mean dose at each scored zone. Selection of the number of histories was based on reducing this value to less than 1%. In step III the phase space file from the previous step was given as input to DOSXYZnrc program which continues the transport of phase space particles into the phantom defined as 3D voxels of size 2.5 x 2.5 x 5.0 mm 3 . The 5 mm grid size was in the direction of depth. Typically this process computed 4.7 billion histories in 116 hours. From the 3D dose matrix generated by DOSXYZnrc, depth dose and beam profile data were derived. In addition, cross sectional data was extracted to generate isodose distribution in planes perpendicular to the beam.
Our current design of the MLC for our cobalt machine consisted of 20 leaf pairs of 15 cm leaf length to cover an area of 20 cm x 20 cm field size with a 5 cm leaf over travel from the central axis. The leaf banks would be mounted on 1 cm thick Lucite plate, which was at a distance of 55 cm from the source, attached to the tray holder with two aluminum brackets. The total weight of the collimator with leaves made of low melting point alloy (LMPA), also known as Cerrobend, would be about 30 kg. The design of the proposed MLC is general enough such that when the collimator is built, it can be fitted on to similar machines.
We first constructed a wooden model with the help of local carpenters. The model allowed engineers to grasp the functionality and foresee engineering issues. [Figure 3] shows a picture of this wooden model. In the current study, prototype leaves were made with both LMPA and lead metal.
In the current design, the leaf height was arbitrarily chosen as 6 cm based on 2% leakage for lead. The leaf thickness (width) was 7 mm projecting a size of 1 cm at the isocenter. The leaf had a rounded edge with a 9.25 cm radius of curvature. The tongue and grooves were 2 mm wide and 4 cm height. [Figure 4] shows the cross section and other dimensions of the leaf. The leaf had an additional groove in the bottom to help make the leaf travel on a rail. In addition, the tongue, the groove and the aluminum bracket helps to hold it in position.
This paper addresses the radiation characteristics of this design using MC method. For simulation, the MLC leaf was at a distance of 48 cm from the source. To compare the current MLC design with a near perfect MLC at the same position, we have simulated a tungsten leaf MLC with a leaf height of 10 cm and a source distance of 45 cm. This would have a transmission of the primary radiation beam of only 0.002 %. This tungsten collimator did not have any tongues and grooves in the MC modeling but had rounded leaves. Since it was at the same position of the tray holder it would have the fluence modifications and penumbral changes and even leakage through the closed leaves. To study the leakage characteristics, in this simulation, MLC leaves were adjusted to a 4 cm x 4 cm opening while the jaws were fixed at 10 cm x 10 cm size. While the simulation of an irregularly shaped field would be more realistic, proper analysis cannot be made without the availability of accurate measured data and hence the choice of a regular small field inside the 10 cm x 10 cm field. In addition, arbitrary shaped irregular fields cannot have predictable penumbra because of complex scatter contributions.
| Results|| |
[Figure 5] shows a comparison of 10 cm x 10 cm beam spectrum at isocenter with data from Mora et al.  It can be seen that there is close agreement between the two spectra. Deviations could be due to differences between simulation dimensions of the source and housing. [Figure 6] compares 10 cm x 10 cm field depth dose calculation results of direct BEAMnrc simulation with 5 billion histories compared with standard depth dose tables for cobalt beams from BJR25.  The averaged agreement was within 1%. [Figure 7] shows a comparison of 5 cm x 5 cm and 15 cm x 15 cm field depth dose calculations using 1 billion histories with BJR25  data. The agreement for 15 cm x 15 cm data was within 1%. However, the deviation with 5 cm x 5 cm data was nearly 2%.
[Figure 8] shows isodose distributions for a 10 cm x10 cm field at 4.75 cm and 14.75 cm depths. In comparison with shallow depth, there was a blooming effect of the isodose lines at larger depth as expected. [Figure 9] shows beam profiles for 10 cm x 10 cm field at depths of d max , 5 cm, 10 cm, 15 cm and 20 cm depths. Since the voxel size was relatively large, there was noise in the data due to lack of resolution and lack of adequate statistics. Hence, smoothed profile data has been presented.
[Figure 10] shows typical leaf position for a 4 cm x 4 cm opening for the tungsten leaves of 10 cm height and when the inner and outer jaws were set to 10 cm x 10 cm field.
[Figure 11] shows the comparison of isodose distributions for the tungsten leaf MLC and the LMPA leaf MLC. It can be seen that the LMPA distribution has a wider 5% isodose surface compared to the tungsten plot. Roughly, the 5% isodose has spread around the 4 x 4 field size by nearly 1.5 cm all around. There was not much difference in the isodoses in the central 4 x 4 area. More accurate observations can be done by plotting the profiles as described below.
[Figure 12] shows the radiation characteristics, plotted as profiles, of the proposed LMPA leaves compared with tungsten leaves. From the x-axis profile it was observed additional radiation leakage due to the LMPA collimator was about 4% outside the defined field. As shown in the profiles of the y-axis, the leakage of LMPA compared with tungsten adds an additional 10% leakage in the area of leaf end transmission.
| Discussion|| |
As shown in [Figure 6],[Figure 7], the percent depth dose data obtained from the current simulation closely matched the data from BJR25. The BJR25  data for depth dose are considered a good average standard for these machines. Differences if any from the BJR25 data can be attributed to differences in individual machines. A comparison of differences in scatter factors was recently presented by Senthilkumar and Ramakrishnan  and previously by Sharma et al.  According to these studies, although scatter factors differ by up to 2% due to electron contamination, the depth dose data are supposed to be very similar if not identical. The deviation of 2% from BJR25  data we have observed with 5 cm x 5 cm data needs further investigation especially for small field sizes. This could be important especially for irregular MLC fields.
Most Linear accelerators are now fitted with MLC. The advantages of MLC in limiting the radiation field to the tumor shape are well known. A recent article by Adams and Warrington  compared several conformal and IMRT plans between cobalt and linear accelerator. They have demonstrated that high quality radiotherapy treatments for cobalt units would be feasible if the cobalt unit were to be fitted with proper beam blocking and compensating systems.
In addition, the use of MLC in principle, allows the cobalt unit to be used for intensity modulated radiation therapy (IMRT). There was only one commercial MLC known as WIMRT (WIMRT DMLC by TOPSLANE International, Clearwater, Florida, USA. Distributed by Acceletronics, Mumbai, India) from Germany, designed for linacs, that was advertised as suitable for cobalt machines. This was designed for radiosurgery and IMRT applications with 3 mm leaf width and can cover only 12 cm of width.
There are no known publications on any MLC for cobalt machines except on an indigenous development from India.  The paper demonstrated that cobalt conformal plans with MLC can match similar plans from linac. This MLC produced a limited projected field size of only 13 cm x 13 cm at the isocenter.
Although there are many ways to design MLC, the current design aims to retrofit the collimator to any cobalt machine quickly without requiring alteration of the machine. It entails the MLC mounted on the tray holder. This puts a limitation on the design and hence compromises the properties of the MLC. Also the collimator would be bulky compared to a collimator built into the machine at the time of machine fabrication rather than retrofit.
Most MLC designs use tungsten leaves. However, tungsten is expensive and hard to machine. Since most radio surgery cones are made using LMPA, this material was chosen for the prototype design of the leaves. Our experience with LMPA has been that very smooth surface and tongue and groove construction was possible by carefully pouring the blocks and further machining on a milling machine. The material can be sprayed with a thin plastic coat for additional durability. When lead is used for the 6 cm thick MLC leaves, only 2% of the primary radiation of cobalt-60 with an average energy of 1.25 MeV would be transmitted; whereas with LMPA, a 3.5% transmission would result. Although lead offers this advantage, the metal is relatively soft. A steel foil could be used as a covering for lead leaves. Although steel foil covering will slightly increase the transmission penumbra, its incorporation may not adversely affect the overall radiation characteristics of the beam. However, making a tight steel covering can be an engineering task.
It was encouraging to find in this study that there was no visible increase in penumbra of the beam when the beam profiles shapes of tungsten and LMPA were compared as can be seen from [Figure 11],[Figure 12]. However, LMPA showed a leakage of nearly 4% assessed from the blocked area.
When two opposite rounded leaves join together, the junction would have less material to attenuate the primary beam. This created a large amount of leakage but was very narrow in extent. This leaf end transmission is usually only about 2% for linear accelerators due to the very small source size. For the current simulation with the cobalt teletherapy machine, this effect is significant. This was as expected and can be eliminated by using the primary collimator to block the leaf gap form opposite leaves. Since cobalt machines do not have asymmetric jaws, the leaf-end transmission can only be partially reduced.
In conclusion, it was possible to generate accurate data for treatment planning purposes using the MC approach. The MC simulation has proved helpful to evaluate MLC design for cobalt-60 teletherapy machine. Without the use of such simulation it was not possible to assess the radiation properties of the MLC design especially near beam edges where critical structures could lie.
While this study focused on finding only leakage properties, further studies are needed to determine photon fluence modifications made by the MLC for various field configurations. After completion of the proposed MLC fabrication, another MC simulation would be desirable, using the exact fabricated parameters, to compare the simulation with measurements.
| Acknowledgments|| |
The authors wish to thank Dr. Rabi Raja Singh, Associate Professor of Medical Physics, Christian Medical College, Vellore and Dr. A. R. Reddy, Retired Director of Defence Research Labortatory, Jodhpur, for their valuable comments and corrections.
| References|| |
|1.||Ravichandran R. Has the time come for doing away with Cobalt-60 teletherapy for cancer treatments. J Med Phys 2009;34:63-5. |
|2.||Han K, Ballon D, Chui C, Mohan R. Monte Carlo simulation of cobalt-60 beam. Med Phys 1987;14:414-8. |
|3.||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. |
|4.||Mora GM, Maio A, Rogers DW. Monte Carlo simulation of a typical 60 Co therapy source. Med Phys 1999;26:2494-502. |
|5.||Rogers DW, Ewart GM, Bielajew AF, Van Dyk G. Calculation of electron contamination in a 60 Co therapy beam. Proceedings of the IAEA International Symposium on Dosimetry in Radiotherapy, IAEA, Vienna, 1988, Vol. 1, pp.303-312. |
|6.||Ahmad Al-Basheer. A New Approach For Reducing Scattered Photons And Electron Contamination In Cobalt-60 Therapy Beam. M.S. Thesis, University Of Florida, 2004. Available from: http://etd.fcla.edu/UF/UFE0004410/albasheer_a.pdf. [cited in 2004]. |
|7.||de Luelmo SC. Characterization of the 60Co therapy unit Siemens Gammatron 1 using BEAMnrc Monte Carlo simulations. Thesis for Master of Science in Medical Radiation Physics, Stockholm University, Sweden: 2006. |
|8.||Mirσ R, Soler J, Gallardo S, Campayo JM, Dνez S, Verdϊ G. MCNP simulation of a Theratron 780 radiotherapy unit. Radiat Prot Dosimetry 2005;116:65-8. |
|9.||Sichani BT, Sohrabpour M. Monte Carlo dose calculations for radiotherapy machines: Theratron 780-C teletherapy case study. Phys Med Biol 2004;49:807-18. |
|10.||Mesbahi A, Allahverdi M, Gheraati H. Monte Carlo dose calculations in conventional thorax fields for 60Co photons. Radiat Med 2005;23:341-50. |
|11.||Central axis depth dose data for use in radiotherapy. Br J Radiol 1996; Supplement 25. |
|12.||Senthilkumar S, Ramakrishnan V. Design of mini phantom and measurement of cobalt-60 beam data parameters. J Med Physics 2008;33:100-7. |
|13.||Sharma SD, Tripathi UB, Kumar P, Govinda Rajan KN, Bhatt BC. Measured head scatter correction factors for tele-cobalt units. J Med Physics 2000;25:6-12. |
|14.||Adams EJ, Warrington AP. A comparison between cobalt and linear accelerator-based treatment plans for conformal and intensity-modulated radiotherapy. Br J Radiol 2008;81:304-10. |
|15.||Singh IR, Ravindran BP, Ayyangar KM. Design and development of motorized MLC for telecobalt unit. Technol Cancer Res Treat 2006;5:597-605. |
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7], [Figure 8], [Figure 9], [Figure 10], [Figure 11], [Figure 12]
|This article has been cited by|
||Monte Carlo simulation based study of a proposed multileaf collimator for a telecobalt machine
| || Sahani, G., Dutt Sharma, S., Dash Sharma, P.K., Sharma, D.N., Hussain, S.A. |
| ||Medical Physics. 2013; 40(2): art-021705 |
||Monte Carlo simulation based study of a proposed multileaf collimator for a telecobalt machine
| ||G. Sahani,Sunil Dutt Sharma,P. K. Dash Sharma,D. N. Sharma,S. A. Hussain |
| ||Medical Physics. 2013; 40(2): 021705 |
|[Pubmed] | [DOI]|