Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 7-Day Trial for You or Your Team.

Learn More →

MCHP (Monte Carlo + Human Phantom): Platform to facilitate teaching nuclear radiation physics

MCHP (Monte Carlo + Human Phantom): Platform to facilitate teaching nuclear radiation physics Some concepts in nuclear radiation physics are abstract and intellectually demanding. In the present paper, an “MCHP platform” (MCHP was an acronym for Monte Carlo simula- OPENACCESS tions + Human Phantoms) was proposed to provide assistance to the students through visu- alization. The platform involved Monte Carlo simulations of interactions between ionizing Citation: Shahmohammadi Beni M, Watabe H, Krstic D, Nikezic D, Yu KN (2021) MCHP (Monte radiations and the Oak Ridge National Laboratory (ORNL) adult male human phantom. As Carlo + Human Phantom): Platform to facilitate an example to demonstrate the benefits of the proposed MCHP platform, the present paper teaching nuclear radiation physics. PLoS ONE investigated the variation of the absorbed photon dose per photon from a Cs source in 16(9): e0257638. https://doi.org/10.1371/journal. three selected organs, namely, brain, spine and thyroid of an adult male for concrete and pone.0257638 lead shields with varying thicknesses. The results were interesting but not readily compre- Editor: Mahendra Singh Dhaka, Mohanlal Sukhadia hensible without direct visualization. Graphical visualization snapshots as well as video clips University, INDIA of real time interactions between the photons and the human phantom were presented for Received: June 18, 2021 the involved cases, and the results were explained with the help of such snapshots and Accepted: September 4, 2021 video clips. It is envisaged that, if the platform is found useful and effective by the readers, Published: September 17, 2021 the readers can also propose examples to be gradually added onto this platform in future, Copyright:© 2021 Shahmohammadi Beni et al. with the ultimate goal of enhancing students’ understanding and learning the concepts in an This is an open access article distributed under the undergraduate nuclear radiation physics course or a related course. terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Introduction Data Availability Statement: All data files are available from Figshare: https://figshare.com/ Nuclear radiation physics is an interesting subject in physics but some concepts are abstract articles/media/MCHP_Monte_Carlo_Human_ and intellectually demanding. The main topics in nuclear radiation physics include the nature Phantom_Platform_to_facilitate_teaching_ of ionizing radiations, their interactions and detection. Notably, the stochastic nature of radio- nuclear_radiation_physics/13798268. activity as well as interactions between ionizing radiation with matter have made their direct Funding: The present work was supported by the visualization relatively tedious and challenging, and as a result many concepts involved in Special Grant for the Development of Virtual nuclear radiation physics are being taught nowadays using average properties and/or probabil- Teaching and Learning (VTL) no. 6430120 from ities. For example, Fig 1 shows the photon interaction cross-section data for some elements the University Grants Committee of Hong Kong 1 12 14 16 SAR (https://www.ugc.edu.hk/eng/ugc/about/ H, C, N and O that are major nuclei in a modeled human cell or organ. These were PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 1 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics press_speech_other/press/2021/pr20210121. obtained from the National Institute of Standards and Technology (NIST) XCOM library html) (awarded to K.N. Yu). The funder had no role homepage: https://physics.nist.gov/PhysRefData/Xcom/html/xcom1.html. The four different in study design, data collection and analysis, interactions shown in the figure include coherent, incoherent, photoelectric effect and pair decision to publish, or preparation of the production. In another example, in the discussion of the linear attenuation coefficient μ of manuscript. photons in a target material, which is defined as the product n × σ, where n is the number den- Competing interests: The authors have declared sity of atomic nuclei in the target material (i.e., number of atomic nuclei per unit volume), that no competing interests exist. while σ is the interaction cross section between the photons and the atomic nuclei. Such treat- ments of the topics might be the most convenient ones for teaching, but the general lack of visualization has inevitably made the concepts abstract, which might not have helped the stu- dents’ understanding and learning. In the present paper, we proposed an “MCHP platform” (MCHP is acronym for Monte Carlo simulations + Human Phantoms) for teaching such abstract concepts in an undergradu- ate nuclear radiation physics course (or a related course) through general provision of 1 12 14 16 Fig 1. Interaction cross-section data for H, C, N and O that are the major nuclei in a modeled human cell or organ. The four different interactions shown in the figure include coherent, incoherent, photoelectric effect and pair production. These were obtained from the National Institute of Standards and Technology (NIST) XCOM library homepage: https://physics.nist.gov/PhysRefData/Xcom/html/xcom1.html. https://doi.org/10.1371/journal.pone.0257638.g001 PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 2 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics Fig 2. The Oak Ridge National Laboratory (ORNL) adult male human phantom [1] used in the present work. https://doi.org/10.1371/journal.pone.0257638.g002 visualization. To further arouse the students’ interest, the platform involved the Oak Ridge National Laboratory (ORNL) adult male human phantom that we have previously developed [1] as shown in Fig 2, and the topics could be discussed in terms of, for examples, the radiation dose delivered to different organs in the human phantom in real time. In this way, the students could likely perceive that the discussions had a direct relevance to themselves. As mentioned above, interactions between ionizing radiation with matter are stochastic, so Monte Carlo sim- ulations were employed to achieve direct visualization. An example was also provided in this preliminary work to help teach the concepts of inter- actions between photons and matter, attenuation coefficient and radiation shielding against photons [2]. The materials and methods employed in working out this example would be described in section 2 below, while the results and discussion will be presented in section 3 below. It is envisaged that, if the platform is found useful and effective by the readers, the read- ers can also propose examples to be gradually added onto this platform in future, with the ulti- mate goal of enhancing students’ understanding and learning the concepts in an undergraduate nuclear radiation physics course or a related course. Methodology and numerical results We refer interested readers to the Refs. [1–4] regarding the modelling and irradiation of the ORNL adult male human phantom. In the present example, the Monte Carlo N-Particle (MCNP) code version 5 was used to simulate the transport and interaction of high energy pho- tons. Similar to our previous work [1], an isotropic point-like Cs radioactive source was used, which emitted photons with energy of 661.6 keV in all directions. The source to phantom distance was kept fixed at 100 cm from the face of phantom (please refer to Ref. [2] for the detailed setup). The two shielding materials considered in the present example were (1) con- crete and (2) lead, both with thicknesses of 10 and 20 cm [5]. Three main organs, namely, brain, spine and thyroid were used to demonstrate the impact of the ionizing radiation on the human body. The dose from the photons in each of these three selected organs were determined. Fig 3 shows the absorbed photon dose per photon in the brain, spine and thyroid of the adult male phantom for no shielding, 10 and 20 cm thick concrete and lead shielding. Interest- ingly, Fig 3 shows that the absorbed photon doses per photon in the organs would be reduced when a shielding material was inserted between the source and the phantom, and the absorbed photon doses per photon would be reduced more for a thicker shield. Furthermore, the lead shields were more effective than concrete shields. The energy of the employed Cs source is 661.6 keV and at this energy the mass attenuation coefficient for concrete and lead would be −2 −1 2 about 8.236×10 and 1.248×10 cm /g, respectively. In addition, the density of concrete and lead would be 2.35 and 11.35 g/cm , respectively, which leads to much larger linear attenuation PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 3 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics Fig 3. Absorbed photon dose per photon in the brain, spine and thyroid of the adult male phantom for no shielding, 10 and 20 cm thick concrete and lead shielding. The error bars represent the statistical uncertainties associated with Monte Carlo simulation. Some of the error bars are smaller than symbols. https://doi.org/10.1371/journal.pone.0257638.g003 coefficient (mass attenuation coefficient × density) for lead than concrete for shielding against gamma rays. This can be noticed from the obtained results. Although such results would usu- ally deem to be “expected”, the complete explanation might be more tedious. For example, it would not be trivial whether the reduction in the absorbed doses were due to the absorption of the photons due to photoelectric effect (and pair production as well for general consideration; pair production was not involved in the current example involving the photons from Cs) within the shielding material, or due to the absorption of photon energy due to Compton scat- tering within the shielding material. Furthermore, while photoelectric effect (and pair produc- tion for general consideration) will attenuate the number of photons reaching the human phantom (which could also be observed from the visualization in Fig 4 below) and will thus attenuate the absorbed dose in the human phantom, Compton scattering will reduce the energy of the photons reaching the human phantom (which could also be observed from the visualization in Fig 5 below) and whether that will attenuate the absorbed dose in the human phantom will further depend on the interaction cross section of the photons with the reduced energies [6, 7]. For the energy range of photons considered in the current example (where the original source was Cs), incoherent (Compton) scattering is a main photon-matter interaction mechanism [8–15] (see the cross-section data shown in Fig 2), and the dose will be delivered PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 4 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics Fig 4. Graphical representation of photon tracks for the cases with (a) no shielding, (b) 10 cm concrete, (c) 20 cm concrete, (d) 10 cm lead and (e) 20 cm lead shielding. Blue dots: positions of photon interactions; green lines: links between positions of multiple interaction (i.e., blue dots); green open circle: positions at which photons would not undergo further multiple scattering. Some symbols may have overlapped with one another. For a better understanding, the readers could cross check with the video clips, where the symbols are displayed in sequence. https://doi.org/10.1371/journal.pone.0257638.g004 PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 5 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 6 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics Fig 5. Graphical representation of the energy of interacting photons for the cases with (a) no shielding, (b) 10 cm concrete, (c) 20 cm concrete, (d) 10 cm lead and (e) 20 cm lead shielding. The lighter yellow dots represent photons with higher energy, while darker yellow dots represent photons with lower energy. https://doi.org/10.1371/journal.pone.0257638.g005 to the organs in the phantom mainly due to photoelectric effect and incoherent scattering. Unfortunately, the numerical absorbed photon doses alone could not help the students’ under- standing and learning on photon scattering and spreading, which could however be achieved with the help of visualization in Figs 4 and 5 below. To sum up, the results in Fig 3 were interesting, but it might not be straightforward to understand these results without direct visualization. In section 3 below, as a demonstration of the benefits from the direct visualization capability proposed in the current MCHP platform, graphical visualization snapshots as well as video clips of real time interactions between the photons and the human phantom were presented for the involved cases, and the results were explained with the help of such snapshots and video clips. Visualization and discussion To illustrate the added values of the currently proposed MCHP platform, two graphical visuali- zation snapshots for the five involved cases are shown in Figs 4 and 5, namely, (a) no shielding, (b) 10 cm concrete, (c) 20 cm concrete, (d) 10 cm lead and (e) 20 cm lead shielding (distance between the shield and the phantom was ~60 cm). The visualizations were prepared using Vised [16] that came in bundle with the MCNP package. In addition, video clips on interac- tions between the photons and the phantom were also prepared, corresponding to the cases with no shielding (Fig 4A), 10 cm thick concrete shield (Fig 4B) and 10 cm thick lead shield (Fig 4D). Interested readers can download these video clips from https://figshare.com/articles/ media/MCHP_Monte_Carlo_Human_Phantom_Platform_to_facilitate_teaching_nuclear_ radiation_physics/13798268. In Fig 4 and the video clips mentioned above, the blue dots repre- sented the positions of photon interactions, and the green lines connected the positions of multiple interaction (i.e., blue dots). Those positions at which the photons would not undergo further multiple scattering were covered with a green open circle. It is remarked here that some symbols in Fig 4 overlapped with one another because of the high densities of dots, lines and open circles. For a better understanding, the readers could cross check with the video clips, where the symbols were shown in sequence. Fig 4(A) is for the case with no shielding. Fig 4(B) and 4(C) represent the cases for 10 and 20 cm thick concrete shields, respectively, and the number of radiation interactions with human phantom has been significantly reduced. Fig 4(D) and 4(E) represent the cases for 10 and 20 cm thick lead shields, respectively, and almost no photons reached the human phantom. These visual- izations clearly demonstrated that absorption of photons due to photoelectric effect within the shielding material contributed to the reduction in the absorbed doses in the human phantom. On the other hand, Fig 5 shows the energy of interacting photons, using dots with different colors to represent photons with different energies. Fig 5(A) is for the case with no shielding, and the energies of the photons which interacted with the human phantom were uniform within the phantom and relatively high (as shown by lighter yellow color dots) when compared to Fig 5(B) and 5(C). Fig 5(B) and 5(C) represent the cases for 10 and 20 cm thick concrete shields, respectively, and the energies of a substantial portion of the photons which interacted with the human phantom became lower (as shown by darker yellow color dots) after having interacted with the shielding material. Fig 5(D) and 5(E) represent the cases for 10 and 20 cm thick lead shields, respectively, which show further reduction in the number and energy of the photons which interacted with the human phantom. These visualizations clearly demonstrated PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 7 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics that the absorption of photon energy due to Compton scattering within the shielding material also contributed to the reduction in the absorbed doses in the human phantom. The results in Figs 4 and 5 clearly revealed that the reduction in the absorbed doses per pho- ton provided by the shielding materials were due to both absorption of the photons due to photoelectric effect, and absorption of photon energy due to Compton scattering within the shielding material. Without these visualizations, it would not be trivial to make such conclu- sions. As described above, the reduction in the absorbed doses were traditionally “explained” by the mass attenuation coefficient or the linear attenuation coefficient (mass attenuation coefficient × density), but this explanation did not answer the question whether the reduction in the absorbed doses were due to the absorption of the photons or due to absorption of pho- ton energy within the shielding material. In particular, the radiation spreading of photons from the Cs source and their interactions with the phantom would enhance the understand- ing of students on the importance of incoherent scattering, which was a dominating photon- matter interaction mechanism over a large energy range for most commonly encountered nuclei present in the human body and other living organisms. The video clips provided a much clearer visualization of the radiation interactions with the phantom, the stochastic nature of these interactions and the spreading of photons. In addition, it is now directly observable that after incoherent scattering of photons with the surrounding air or the shielding material, the photons would no longer be mono-energetic. Without the visualization aids, the energy distribution could only be plotted for one domain at one time. In contrast, with the help of direct visualization, changes in the photon energy in multiple domains could be examined at the same time, such as surrounding air, shielding material and human phantom (e.g., see Fig 5B and 5C). From the above discussion, we genuinely felt that the MCHP platform proposed in the cur- rent paper could help generate “case studies” which could help undergraduate student learn the subject of nuclear radiation physics. The case studies could help explain concepts which were relatively difficult to understand using mathematical equations alone, particularly those which only described average properties, and could arouse the students’ interest since the stu- dents could easily perceive that the discussions had a direct relevance to themselves. If this MCHP platform is welcome by the community teaching nuclear radiation physics and if the platform can become popular in future, a website will be created to host the created case stud- ies, probably with the discussions as well, which can be used by members of the teaching com- munity. We envisage that the development of the MCHP platform can benefit both the teachers and students in a nuclear radiation physics course or a related course. Conclusions The present work introduced a novel teaching platform named MCHP (acronym for Monte Carlo + Human Phantoms). Through visualization, this platform helped students better understand interactions between nuclear radiation and the human body. In the present work, irradiation by a Cs point source through different shielding materials with different thick- nesses was used to demonstrate the use of the MCHP platform. Our previously developed ORNL human phantom model was employed. The MCHP platform would be very useful for pedagogical purposes. In future, we would incorporate other well-known examples in the field of nuclear radiation physics into the platform. Acknowledgments We acknowledge the support of the Neutron computer cluster from the Department of Phys- ics, City University of Hong Kong, for the computational work involved in this paper. PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 8 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics Author Contributions Conceptualization: Mehrdad Shahmohammadi Beni, Kwan Ngok Yu. Data curation: Mehrdad Shahmohammadi Beni. Formal analysis: Mehrdad Shahmohammadi Beni. Funding acquisition: Kwan Ngok Yu. Investigation: Mehrdad Shahmohammadi Beni. Methodology: Mehrdad Shahmohammadi Beni. Project administration: Kwan Ngok Yu. Resources: Hiroshi Watabe. Software: Dragana Krstic, Dragoslav Nikezic. Supervision: Hiroshi Watabe, Kwan Ngok Yu. Validation: Dragana Krstic, Dragoslav Nikezic, Kwan Ngok Yu. Visualization: Mehrdad Shahmohammadi Beni, Kwan Ngok Yu. Writing – original draft: Mehrdad Shahmohammadi Beni, Kwan Ngok Yu. Writing – review & editing: Hiroshi Watabe, Dragana Krstic, Dragoslav Nikezic, Kwan Ngok Yu. References 1. Krstic D, Nikezic D. Input files with ORNL—mathematical phantoms of the human body for MCNP-4B. Comput. Phys. Commun. 2007; 176: 33–37. 2. Akhlaghi P, Miri-Hakimabad H, Rafat-Motavalli L. Effects of shielding the radiosensitive superficial organs of ORNL pediatric phantoms on dose reduction in computed tomography. J. Med. Phys. 2014; 39: 238–246. https://doi.org/10.4103/0971-6203.144490 PMID: 25525312 3. Shahmohammadi Beni M, Ng CYP, Krstic D, Nikezic D, Yu KN. Conversion coefficients for determina- tion of dispersed photon dose during radiotherapy: NRUrad input code for MCNP. PLoS ONE 2017; 12: e0174836. https://doi.org/10.1371/journal.pone.0174836 PMID: 28362837 4. Shahmohammadi Beni M, Krstic D, Nikezic D, Yu KN. A comparative study on dispersed doses during photon and proton radiation therapy in pediatric applications. PLoS ONE 2021; 16: e0248300. https:// doi.org/10.1371/journal.pone.0248300 PMID: 33690664 5. Singer G, Laurence GC. Materials and Methods of X-ray Protection. Radiology 1946; 46: 69–76. https:// doi.org/10.1148/46.1.69 PMID: 21011652 6. Shahmohammadi Beni M, Krstic D, Nikezic D, Yu KN. Realistic dosimetry for studies on biological responses to X-rays andγ-rays. J. Radiat. Res. 2017; 58: 729–736. https://doi.org/10.1093/jrr/rrx019 PMID: 28444359 7. Shahmohammadi Beni M, Krstic D, Nikezic D, Yu KN. Modeling kV X-ray-induced coloration in radio- chromic films. Appl. Sci. 2018; 8: 106. 8. Shahmohammadi Beni M, Krstic D, Nikezic D, Yu KN. Monte Carlo studies on photon interactions in radiobiological experiments. PLoS ONE 2018; 13: e0193575. https://doi.org/10.1371/journal.pone. 0193575 PMID: 29561871 9. Shahmohammadi Beni M, Hau TC, Krstic D, Nikezic D, Yu KN. Monte Carlo studies on neutron interac- tions in radiobiological experiments. PLoS ONE 2017; 12: e0181281. https://doi.org/10.1371/journal. pone.0181281 PMID: 28704557 10. Nikezic D, Shahmohammadi Beni M, Krstic D, Yu KN. Characteristics of protons exiting from a polyeth- ylene converter irradiated by neutrons with energies between 1 keV and 10 MeV. PLoS ONE 2016; 11: e0157627. https://doi.org/10.1371/journal.pone.0157627 PMID: 27362656 11. Hubbell JH. Photon mass attenuation and energy-absorption coefficients. Int. J. Appl. Radiat. Isot. 1982; 33: 1269–1290. PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 9 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics 12. Xu Z, Ni K, Mao J, Luo T, Lin W. Monte Carlo Simulations Reveal New Design Principles for Efficient Nanoradiosensitizers Based on Nanoscale Metal–Organic Frameworks. Adv. Mater. 2021: 2104249. https://doi.org/10.1002/adma.202104249 PMID: 34432917 13. Frosio T, Bertreix P, Menaa N, Thomas S, Eberhardt H, Endres J. Photons fluence to local skin Dose coefficients and benchmark with three Monte-Carlo codes. Application to the computation of radioactive material transport limits. Appl. Radiat. Isot. 2021; 176: 109892. https://doi.org/10.1016/j.apradiso.2021. 109892 PMID: 34418729 14. Rostamani N, Khabaz R. Monte Carlo simulation estimates of absorbed dose in human organs due to the external exposure by decorative granite stones. Radiat. Phys. Chem. 2021; 189: 109702. 15. Shahmohammadi Beni M, Krstic D, Nikezic D, Yu KN. Medium-thickness-dependent proton dosimetry for radiobiological experiments. Sci. Rep. 2019; 9: 11577. https://doi.org/10.1038/s41598-019-48100-9 PMID: 31399622 16. The MCNP Vised website. Available from: http://www.mcnpvised.com/ PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 10 / 10 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png PLoS ONE Public Library of Science (PLoS) Journal

MCHP (Monte Carlo + Human Phantom): Platform to facilitate teaching nuclear radiation physics

Download PDF
10 pages

Loading next page...
 
/lp/public-library-of-science-plos-journal/mchp-monte-carlo-human-phantom-platform-to-facilitate-teaching-nuclear-YbgKIyC0sE

References (17)

Publisher
Public Library of Science (PLoS) Journal
Copyright
Copyright: © 2021 Shahmohammadi Beni et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: All data files are available from Figshare: https://figshare.com/articles/media/MCHP_Monte_Carlo_Human_Phantom_Platform_to_facilitate_teaching_nuclear_radiation_physics/13798268. Funding: The present work was supported by the Special Grant for the Development of Virtual Teaching and Learning (VTL) no. 6430120 from the University Grants Committee of Hong Kong SAR (https://www.ugc.edu.hk/eng/ugc/about/press_speech_other/press/2021/pr20210121.html) (awarded to K.N. Yu). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist.
eISSN
1932-6203
DOI
10.1371/journal.pone.0257638
Publisher site
See Article on Publisher Site

Abstract

Some concepts in nuclear radiation physics are abstract and intellectually demanding. In the present paper, an “MCHP platform” (MCHP was an acronym for Monte Carlo simula- OPENACCESS tions + Human Phantoms) was proposed to provide assistance to the students through visu- alization. The platform involved Monte Carlo simulations of interactions between ionizing Citation: Shahmohammadi Beni M, Watabe H, Krstic D, Nikezic D, Yu KN (2021) MCHP (Monte radiations and the Oak Ridge National Laboratory (ORNL) adult male human phantom. As Carlo + Human Phantom): Platform to facilitate an example to demonstrate the benefits of the proposed MCHP platform, the present paper teaching nuclear radiation physics. PLoS ONE investigated the variation of the absorbed photon dose per photon from a Cs source in 16(9): e0257638. https://doi.org/10.1371/journal. three selected organs, namely, brain, spine and thyroid of an adult male for concrete and pone.0257638 lead shields with varying thicknesses. The results were interesting but not readily compre- Editor: Mahendra Singh Dhaka, Mohanlal Sukhadia hensible without direct visualization. Graphical visualization snapshots as well as video clips University, INDIA of real time interactions between the photons and the human phantom were presented for Received: June 18, 2021 the involved cases, and the results were explained with the help of such snapshots and Accepted: September 4, 2021 video clips. It is envisaged that, if the platform is found useful and effective by the readers, Published: September 17, 2021 the readers can also propose examples to be gradually added onto this platform in future, Copyright:© 2021 Shahmohammadi Beni et al. with the ultimate goal of enhancing students’ understanding and learning the concepts in an This is an open access article distributed under the undergraduate nuclear radiation physics course or a related course. terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Introduction Data Availability Statement: All data files are available from Figshare: https://figshare.com/ Nuclear radiation physics is an interesting subject in physics but some concepts are abstract articles/media/MCHP_Monte_Carlo_Human_ and intellectually demanding. The main topics in nuclear radiation physics include the nature Phantom_Platform_to_facilitate_teaching_ of ionizing radiations, their interactions and detection. Notably, the stochastic nature of radio- nuclear_radiation_physics/13798268. activity as well as interactions between ionizing radiation with matter have made their direct Funding: The present work was supported by the visualization relatively tedious and challenging, and as a result many concepts involved in Special Grant for the Development of Virtual nuclear radiation physics are being taught nowadays using average properties and/or probabil- Teaching and Learning (VTL) no. 6430120 from ities. For example, Fig 1 shows the photon interaction cross-section data for some elements the University Grants Committee of Hong Kong 1 12 14 16 SAR (https://www.ugc.edu.hk/eng/ugc/about/ H, C, N and O that are major nuclei in a modeled human cell or organ. These were PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 1 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics press_speech_other/press/2021/pr20210121. obtained from the National Institute of Standards and Technology (NIST) XCOM library html) (awarded to K.N. Yu). The funder had no role homepage: https://physics.nist.gov/PhysRefData/Xcom/html/xcom1.html. The four different in study design, data collection and analysis, interactions shown in the figure include coherent, incoherent, photoelectric effect and pair decision to publish, or preparation of the production. In another example, in the discussion of the linear attenuation coefficient μ of manuscript. photons in a target material, which is defined as the product n × σ, where n is the number den- Competing interests: The authors have declared sity of atomic nuclei in the target material (i.e., number of atomic nuclei per unit volume), that no competing interests exist. while σ is the interaction cross section between the photons and the atomic nuclei. Such treat- ments of the topics might be the most convenient ones for teaching, but the general lack of visualization has inevitably made the concepts abstract, which might not have helped the stu- dents’ understanding and learning. In the present paper, we proposed an “MCHP platform” (MCHP is acronym for Monte Carlo simulations + Human Phantoms) for teaching such abstract concepts in an undergradu- ate nuclear radiation physics course (or a related course) through general provision of 1 12 14 16 Fig 1. Interaction cross-section data for H, C, N and O that are the major nuclei in a modeled human cell or organ. The four different interactions shown in the figure include coherent, incoherent, photoelectric effect and pair production. These were obtained from the National Institute of Standards and Technology (NIST) XCOM library homepage: https://physics.nist.gov/PhysRefData/Xcom/html/xcom1.html. https://doi.org/10.1371/journal.pone.0257638.g001 PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 2 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics Fig 2. The Oak Ridge National Laboratory (ORNL) adult male human phantom [1] used in the present work. https://doi.org/10.1371/journal.pone.0257638.g002 visualization. To further arouse the students’ interest, the platform involved the Oak Ridge National Laboratory (ORNL) adult male human phantom that we have previously developed [1] as shown in Fig 2, and the topics could be discussed in terms of, for examples, the radiation dose delivered to different organs in the human phantom in real time. In this way, the students could likely perceive that the discussions had a direct relevance to themselves. As mentioned above, interactions between ionizing radiation with matter are stochastic, so Monte Carlo sim- ulations were employed to achieve direct visualization. An example was also provided in this preliminary work to help teach the concepts of inter- actions between photons and matter, attenuation coefficient and radiation shielding against photons [2]. The materials and methods employed in working out this example would be described in section 2 below, while the results and discussion will be presented in section 3 below. It is envisaged that, if the platform is found useful and effective by the readers, the read- ers can also propose examples to be gradually added onto this platform in future, with the ulti- mate goal of enhancing students’ understanding and learning the concepts in an undergraduate nuclear radiation physics course or a related course. Methodology and numerical results We refer interested readers to the Refs. [1–4] regarding the modelling and irradiation of the ORNL adult male human phantom. In the present example, the Monte Carlo N-Particle (MCNP) code version 5 was used to simulate the transport and interaction of high energy pho- tons. Similar to our previous work [1], an isotropic point-like Cs radioactive source was used, which emitted photons with energy of 661.6 keV in all directions. The source to phantom distance was kept fixed at 100 cm from the face of phantom (please refer to Ref. [2] for the detailed setup). The two shielding materials considered in the present example were (1) con- crete and (2) lead, both with thicknesses of 10 and 20 cm [5]. Three main organs, namely, brain, spine and thyroid were used to demonstrate the impact of the ionizing radiation on the human body. The dose from the photons in each of these three selected organs were determined. Fig 3 shows the absorbed photon dose per photon in the brain, spine and thyroid of the adult male phantom for no shielding, 10 and 20 cm thick concrete and lead shielding. Interest- ingly, Fig 3 shows that the absorbed photon doses per photon in the organs would be reduced when a shielding material was inserted between the source and the phantom, and the absorbed photon doses per photon would be reduced more for a thicker shield. Furthermore, the lead shields were more effective than concrete shields. The energy of the employed Cs source is 661.6 keV and at this energy the mass attenuation coefficient for concrete and lead would be −2 −1 2 about 8.236×10 and 1.248×10 cm /g, respectively. In addition, the density of concrete and lead would be 2.35 and 11.35 g/cm , respectively, which leads to much larger linear attenuation PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 3 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics Fig 3. Absorbed photon dose per photon in the brain, spine and thyroid of the adult male phantom for no shielding, 10 and 20 cm thick concrete and lead shielding. The error bars represent the statistical uncertainties associated with Monte Carlo simulation. Some of the error bars are smaller than symbols. https://doi.org/10.1371/journal.pone.0257638.g003 coefficient (mass attenuation coefficient × density) for lead than concrete for shielding against gamma rays. This can be noticed from the obtained results. Although such results would usu- ally deem to be “expected”, the complete explanation might be more tedious. For example, it would not be trivial whether the reduction in the absorbed doses were due to the absorption of the photons due to photoelectric effect (and pair production as well for general consideration; pair production was not involved in the current example involving the photons from Cs) within the shielding material, or due to the absorption of photon energy due to Compton scat- tering within the shielding material. Furthermore, while photoelectric effect (and pair produc- tion for general consideration) will attenuate the number of photons reaching the human phantom (which could also be observed from the visualization in Fig 4 below) and will thus attenuate the absorbed dose in the human phantom, Compton scattering will reduce the energy of the photons reaching the human phantom (which could also be observed from the visualization in Fig 5 below) and whether that will attenuate the absorbed dose in the human phantom will further depend on the interaction cross section of the photons with the reduced energies [6, 7]. For the energy range of photons considered in the current example (where the original source was Cs), incoherent (Compton) scattering is a main photon-matter interaction mechanism [8–15] (see the cross-section data shown in Fig 2), and the dose will be delivered PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 4 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics Fig 4. Graphical representation of photon tracks for the cases with (a) no shielding, (b) 10 cm concrete, (c) 20 cm concrete, (d) 10 cm lead and (e) 20 cm lead shielding. Blue dots: positions of photon interactions; green lines: links between positions of multiple interaction (i.e., blue dots); green open circle: positions at which photons would not undergo further multiple scattering. Some symbols may have overlapped with one another. For a better understanding, the readers could cross check with the video clips, where the symbols are displayed in sequence. https://doi.org/10.1371/journal.pone.0257638.g004 PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 5 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 6 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics Fig 5. Graphical representation of the energy of interacting photons for the cases with (a) no shielding, (b) 10 cm concrete, (c) 20 cm concrete, (d) 10 cm lead and (e) 20 cm lead shielding. The lighter yellow dots represent photons with higher energy, while darker yellow dots represent photons with lower energy. https://doi.org/10.1371/journal.pone.0257638.g005 to the organs in the phantom mainly due to photoelectric effect and incoherent scattering. Unfortunately, the numerical absorbed photon doses alone could not help the students’ under- standing and learning on photon scattering and spreading, which could however be achieved with the help of visualization in Figs 4 and 5 below. To sum up, the results in Fig 3 were interesting, but it might not be straightforward to understand these results without direct visualization. In section 3 below, as a demonstration of the benefits from the direct visualization capability proposed in the current MCHP platform, graphical visualization snapshots as well as video clips of real time interactions between the photons and the human phantom were presented for the involved cases, and the results were explained with the help of such snapshots and video clips. Visualization and discussion To illustrate the added values of the currently proposed MCHP platform, two graphical visuali- zation snapshots for the five involved cases are shown in Figs 4 and 5, namely, (a) no shielding, (b) 10 cm concrete, (c) 20 cm concrete, (d) 10 cm lead and (e) 20 cm lead shielding (distance between the shield and the phantom was ~60 cm). The visualizations were prepared using Vised [16] that came in bundle with the MCNP package. In addition, video clips on interac- tions between the photons and the phantom were also prepared, corresponding to the cases with no shielding (Fig 4A), 10 cm thick concrete shield (Fig 4B) and 10 cm thick lead shield (Fig 4D). Interested readers can download these video clips from https://figshare.com/articles/ media/MCHP_Monte_Carlo_Human_Phantom_Platform_to_facilitate_teaching_nuclear_ radiation_physics/13798268. In Fig 4 and the video clips mentioned above, the blue dots repre- sented the positions of photon interactions, and the green lines connected the positions of multiple interaction (i.e., blue dots). Those positions at which the photons would not undergo further multiple scattering were covered with a green open circle. It is remarked here that some symbols in Fig 4 overlapped with one another because of the high densities of dots, lines and open circles. For a better understanding, the readers could cross check with the video clips, where the symbols were shown in sequence. Fig 4(A) is for the case with no shielding. Fig 4(B) and 4(C) represent the cases for 10 and 20 cm thick concrete shields, respectively, and the number of radiation interactions with human phantom has been significantly reduced. Fig 4(D) and 4(E) represent the cases for 10 and 20 cm thick lead shields, respectively, and almost no photons reached the human phantom. These visual- izations clearly demonstrated that absorption of photons due to photoelectric effect within the shielding material contributed to the reduction in the absorbed doses in the human phantom. On the other hand, Fig 5 shows the energy of interacting photons, using dots with different colors to represent photons with different energies. Fig 5(A) is for the case with no shielding, and the energies of the photons which interacted with the human phantom were uniform within the phantom and relatively high (as shown by lighter yellow color dots) when compared to Fig 5(B) and 5(C). Fig 5(B) and 5(C) represent the cases for 10 and 20 cm thick concrete shields, respectively, and the energies of a substantial portion of the photons which interacted with the human phantom became lower (as shown by darker yellow color dots) after having interacted with the shielding material. Fig 5(D) and 5(E) represent the cases for 10 and 20 cm thick lead shields, respectively, which show further reduction in the number and energy of the photons which interacted with the human phantom. These visualizations clearly demonstrated PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 7 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics that the absorption of photon energy due to Compton scattering within the shielding material also contributed to the reduction in the absorbed doses in the human phantom. The results in Figs 4 and 5 clearly revealed that the reduction in the absorbed doses per pho- ton provided by the shielding materials were due to both absorption of the photons due to photoelectric effect, and absorption of photon energy due to Compton scattering within the shielding material. Without these visualizations, it would not be trivial to make such conclu- sions. As described above, the reduction in the absorbed doses were traditionally “explained” by the mass attenuation coefficient or the linear attenuation coefficient (mass attenuation coefficient × density), but this explanation did not answer the question whether the reduction in the absorbed doses were due to the absorption of the photons or due to absorption of pho- ton energy within the shielding material. In particular, the radiation spreading of photons from the Cs source and their interactions with the phantom would enhance the understand- ing of students on the importance of incoherent scattering, which was a dominating photon- matter interaction mechanism over a large energy range for most commonly encountered nuclei present in the human body and other living organisms. The video clips provided a much clearer visualization of the radiation interactions with the phantom, the stochastic nature of these interactions and the spreading of photons. In addition, it is now directly observable that after incoherent scattering of photons with the surrounding air or the shielding material, the photons would no longer be mono-energetic. Without the visualization aids, the energy distribution could only be plotted for one domain at one time. In contrast, with the help of direct visualization, changes in the photon energy in multiple domains could be examined at the same time, such as surrounding air, shielding material and human phantom (e.g., see Fig 5B and 5C). From the above discussion, we genuinely felt that the MCHP platform proposed in the cur- rent paper could help generate “case studies” which could help undergraduate student learn the subject of nuclear radiation physics. The case studies could help explain concepts which were relatively difficult to understand using mathematical equations alone, particularly those which only described average properties, and could arouse the students’ interest since the stu- dents could easily perceive that the discussions had a direct relevance to themselves. If this MCHP platform is welcome by the community teaching nuclear radiation physics and if the platform can become popular in future, a website will be created to host the created case stud- ies, probably with the discussions as well, which can be used by members of the teaching com- munity. We envisage that the development of the MCHP platform can benefit both the teachers and students in a nuclear radiation physics course or a related course. Conclusions The present work introduced a novel teaching platform named MCHP (acronym for Monte Carlo + Human Phantoms). Through visualization, this platform helped students better understand interactions between nuclear radiation and the human body. In the present work, irradiation by a Cs point source through different shielding materials with different thick- nesses was used to demonstrate the use of the MCHP platform. Our previously developed ORNL human phantom model was employed. The MCHP platform would be very useful for pedagogical purposes. In future, we would incorporate other well-known examples in the field of nuclear radiation physics into the platform. Acknowledgments We acknowledge the support of the Neutron computer cluster from the Department of Phys- ics, City University of Hong Kong, for the computational work involved in this paper. PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 8 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics Author Contributions Conceptualization: Mehrdad Shahmohammadi Beni, Kwan Ngok Yu. Data curation: Mehrdad Shahmohammadi Beni. Formal analysis: Mehrdad Shahmohammadi Beni. Funding acquisition: Kwan Ngok Yu. Investigation: Mehrdad Shahmohammadi Beni. Methodology: Mehrdad Shahmohammadi Beni. Project administration: Kwan Ngok Yu. Resources: Hiroshi Watabe. Software: Dragana Krstic, Dragoslav Nikezic. Supervision: Hiroshi Watabe, Kwan Ngok Yu. Validation: Dragana Krstic, Dragoslav Nikezic, Kwan Ngok Yu. Visualization: Mehrdad Shahmohammadi Beni, Kwan Ngok Yu. Writing – original draft: Mehrdad Shahmohammadi Beni, Kwan Ngok Yu. Writing – review & editing: Hiroshi Watabe, Dragana Krstic, Dragoslav Nikezic, Kwan Ngok Yu. References 1. Krstic D, Nikezic D. Input files with ORNL—mathematical phantoms of the human body for MCNP-4B. Comput. Phys. Commun. 2007; 176: 33–37. 2. Akhlaghi P, Miri-Hakimabad H, Rafat-Motavalli L. Effects of shielding the radiosensitive superficial organs of ORNL pediatric phantoms on dose reduction in computed tomography. J. Med. Phys. 2014; 39: 238–246. https://doi.org/10.4103/0971-6203.144490 PMID: 25525312 3. Shahmohammadi Beni M, Ng CYP, Krstic D, Nikezic D, Yu KN. Conversion coefficients for determina- tion of dispersed photon dose during radiotherapy: NRUrad input code for MCNP. PLoS ONE 2017; 12: e0174836. https://doi.org/10.1371/journal.pone.0174836 PMID: 28362837 4. Shahmohammadi Beni M, Krstic D, Nikezic D, Yu KN. A comparative study on dispersed doses during photon and proton radiation therapy in pediatric applications. PLoS ONE 2021; 16: e0248300. https:// doi.org/10.1371/journal.pone.0248300 PMID: 33690664 5. Singer G, Laurence GC. Materials and Methods of X-ray Protection. Radiology 1946; 46: 69–76. https:// doi.org/10.1148/46.1.69 PMID: 21011652 6. Shahmohammadi Beni M, Krstic D, Nikezic D, Yu KN. Realistic dosimetry for studies on biological responses to X-rays andγ-rays. J. Radiat. Res. 2017; 58: 729–736. https://doi.org/10.1093/jrr/rrx019 PMID: 28444359 7. Shahmohammadi Beni M, Krstic D, Nikezic D, Yu KN. Modeling kV X-ray-induced coloration in radio- chromic films. Appl. Sci. 2018; 8: 106. 8. Shahmohammadi Beni M, Krstic D, Nikezic D, Yu KN. Monte Carlo studies on photon interactions in radiobiological experiments. PLoS ONE 2018; 13: e0193575. https://doi.org/10.1371/journal.pone. 0193575 PMID: 29561871 9. Shahmohammadi Beni M, Hau TC, Krstic D, Nikezic D, Yu KN. Monte Carlo studies on neutron interac- tions in radiobiological experiments. PLoS ONE 2017; 12: e0181281. https://doi.org/10.1371/journal. pone.0181281 PMID: 28704557 10. Nikezic D, Shahmohammadi Beni M, Krstic D, Yu KN. Characteristics of protons exiting from a polyeth- ylene converter irradiated by neutrons with energies between 1 keV and 10 MeV. PLoS ONE 2016; 11: e0157627. https://doi.org/10.1371/journal.pone.0157627 PMID: 27362656 11. Hubbell JH. Photon mass attenuation and energy-absorption coefficients. Int. J. Appl. Radiat. Isot. 1982; 33: 1269–1290. PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 9 / 10 PLOS ONE MCHP platform to facilitate teaching nuclear radiation physics 12. Xu Z, Ni K, Mao J, Luo T, Lin W. Monte Carlo Simulations Reveal New Design Principles for Efficient Nanoradiosensitizers Based on Nanoscale Metal–Organic Frameworks. Adv. Mater. 2021: 2104249. https://doi.org/10.1002/adma.202104249 PMID: 34432917 13. Frosio T, Bertreix P, Menaa N, Thomas S, Eberhardt H, Endres J. Photons fluence to local skin Dose coefficients and benchmark with three Monte-Carlo codes. Application to the computation of radioactive material transport limits. Appl. Radiat. Isot. 2021; 176: 109892. https://doi.org/10.1016/j.apradiso.2021. 109892 PMID: 34418729 14. Rostamani N, Khabaz R. Monte Carlo simulation estimates of absorbed dose in human organs due to the external exposure by decorative granite stones. Radiat. Phys. Chem. 2021; 189: 109702. 15. Shahmohammadi Beni M, Krstic D, Nikezic D, Yu KN. Medium-thickness-dependent proton dosimetry for radiobiological experiments. Sci. Rep. 2019; 9: 11577. https://doi.org/10.1038/s41598-019-48100-9 PMID: 31399622 16. The MCNP Vised website. Available from: http://www.mcnpvised.com/ PLOS ONE | https://doi.org/10.1371/journal.pone.0257638 September 17, 2021 10 / 10

Journal

PLoS ONE – Public Library of Science (PLoS) Journal

Published: Sep 17, 2021

There are no references for this article.