Isotope News, (787), p.20 - 23, 2023/06
Carbon ion radiotherapy has an advantage over conventional radiotherapy such that its superior dose concentration on the tumor helps to reduce unwanted dose to surrounding normal tissues. Nevertheless, a little dose to normal tissues, which is a potential risk of secondary cancer, is still unavoidable. In the current dose assessment, however, only assessment around target volume is performed for the tumor control and prevention of acute radiation injury of fatal organs. We therefore developed a system called RT-PHITS for CIRT to reproduce the carbon ion radiotherapy including the production and transport of secondary particles based on treatment planning data using PHITS. Using this system, whole-body dose assessment of patients in the past carbon ion radiotherapy can be performed. By comparing the dose assessment to the epidemiologic records of the patients, the relation between dose exposure of non-target organs and incidence of side effects such as secondary cancer will be elucidated.
Sabri, A. H. A.*; Tajudin, S. M.*; Aziz, M. Z. A.*; Furuta, Takuya
Radiological Physics and Technology, 16(1), p.109 - 117, 2023/03
The spatial distributions of photon dose rates in a brachytherapy room with an Iridium-192 were simulated by using the particle and heavy ion transport code system (PHITS). A geometry of the brachytherapy room with concrete walls and a sliding lead door was reproduced by tracing the existing room in Advanced Medical and Dental Institute at the Universiti Sains Malaysia in Penang. The simulation results were confirmed by comparing to the measured results using a thermoluminescent dosimeter. The simulation study suggested that an additional layer of 3-mm thick lead at the side wall of the entrance will efficiently reduce the dose outside the entrance due to the photons leaked from the edge of the entrance. Simulation with replacing the source with Cobalt-60 was also conducted and revealed the dose level outside the room was too high compared to regulatory value in the current room configuration.
Sato, Tatsuhiko; Iwamoto, Yosuke; Hashimoto, Shintaro; Ogawa, Tatsuhiko; Furuta, Takuya; Abe, Shinichiro; Kai, Takeshi; Matsuya, Yusuke; Matsuda, Norihiro; Hirata, Yuho; et al.
Journal of Nuclear Science and Technology, 9 Pages, 2023/00
The Particle and Heavy Ion Transport code System (PHITS) is a general-purpose Monte Carlo radiation transport code that can simulate the behavior of most particle species with energies up to 1 TeV (per nucleon for ions). Its new version, PHITS3.31, was recently developed and released to the public. In the new version, the compatibility with high-energy nuclear data libraries and the algorithm of the track-structure modes have been improved. In this paper, we summarize the upgraded features of PHITS3.31 with respect to the physics models, utility functions, and application software introduced since the release of PHITS3.02 in 2017.
Furuta, Takuya; Koba, Yusuke*; Hashimoto, Shintaro; Chang, W.*; Yonai, Shunsuke*; Matsumoto, Shinnosuke*; Ishikawa, Akihisa*; Sato, Tatsuhiko
Physics in Medicine & Biology, 67(14), p.145002_1 - 145002_15, 2022/07
Carbon ion radiotherapy has an advantage over conventional radiotherapy such that its superior dose concentration on the tumor helps to reduce unwanted dose to surrounding normal tissues. Nevertheless, a little dose to normal tissues, which is a potential risk of secondary cancer, is still unavoidable. The Monte Carlo simulation is a good candidate for the tool to assess secondary cancer risk, including the contributions of secondary particles produced by nuclear reactions. We therefore developed a new dose reconstruction system implementing PHITS as the engine. In this system, the PHITS input is automatically created from the DICOM data sets recorded in the treatment planning. The developed system was validated by comparing to experimental dose distribution in water and treatment plan on an anthropomorphic phantom. This system will be used for retrospective studies using the patient data in National Institute for Quantum and Science and Technology.
Katsuyama, Jinya; Yamaguchi, Yoshihito; Nemoto, Yoshiyuki; Furuta, Takuya; Kaji, Yoshiyuki
Proceedings of ASME 2022 Pressure Vessels and Piping Conference (PVP 2022) (Internet), 9 Pages, 2022/07
Igaku Butsuri, 41(4), P. 194, 2021/12
Number of medical uses of Particle and Heavy Ion Transport code System (PHITS) has been increased due to the recent high demands of medical use of radiations. The summary of such research works was described in the review article on medical application of Particle and Heavy Ion Transport code System PHITS published in Radiological Physics and Technology in 2021. There was a request from the editorial board of Japan Society of Medical Physics (JSMP) for writing an introductory article of this article in their internal journal. The research works on medical applications described in the review article, useful functions for medical application in PHITS, and newly opened user forum of PHITS have been introduced.
Tominaga, Masahide*; Nagayasu, Yukari*; Sasaki, Motoharu*; Furuta, Takuya; Hayashi, Hiroaki*; Oita, Masataka*; Nishiyama, Yuichi*; Haga, Akihiro*
Radiological Physics and Technology, 14(4), p.381 - 389, 2021/12
Due to recent advance of diagnostic radiology, the increase of diagnostic radiation exposure to patient becomes problem. Diagnostic Reference Levels has been released to optimized the radiation exposure to patients in Japan recently. The evaluation of entrance surface dose (ESD) is recommended to assess the dose level for general X-ray examination. The ESD can be easily evaluated by multiplying the backscatter factor of the patient body on the free-in-air air kerma. The air kerma free-in-air value used to estimate ESD may contain X-rays scattered from obstacles located at the time of measurement, which may induce non-minor error in assessments. We therefore studied the influence of scattered X-rays on air kerma measurement under various environments (distances, field sizes, and materials). It was found that the dependence on the X-ray energy and field size was different for different materials. The X-ray contamination can be ignored for all the materials when the distance to the scatterer exceeds 35 cm.
Chang, W.*; Koba, Yusuke*; Furuta, Takuya; Yonai, Shunsuke*; Hashimoto, Shintaro; Matsumoto, Shinnosuke*; Sato, Tatsuhiko
Journal of Radiation Research (Internet), 62(5), p.846 - 855, 2021/09
With the aim of developing a revaluation tool of treatment plan in carbon-ion radiotherapy using Monte Carlo (MC) simulation, we propose two methods; one is dedicated to identify realistic-tissue materials from a CT image with satisfying the well-calibrated relationship between CT numbers and stopping power ratio (SPR) provided by TPS, and the other is to estimate dose to water considering the particle- and energy-dependent SPR between realistic tissue materials and water. We validated these proposed methods by computing depth dose distribution in homogeneous and heterogeneous phantoms composed of human tissue materials and water irradiated by a 400 MeV/u carbon beam with 8 cm SOBP using a MC simulation code PHITS and comparing with results of conventional treatment planning system (TPS). Our result suggested that use of water as a surrogate of real tissue materials, which is adopted in conventional TPS, is inadequate for dose estimation from secondary particles because their production rates cannot be scaled by SPR of the primary particle in water. We therefore concluded that the proposed methods can play important roles in the reevaluation of the treatment plans in carbon-ion radiotherapy.
Furuta, Takuya; Sato, Tatsuhiko
Radiological Physics and Technology, 14(3), p.215 - 225, 2021/09
Number of the PHITS users has steadily increased since 2010 from when it is officially counted. Among them, increase of new users in medical physics is outstanding. Many research works in medical physics using PHITS have been published and the applications are widely spread in different fields such as applications to different types of radiotherapy, shielding calculations of medical facilities, application to radiation biology, and research and development of medical tools. In this article, we will introduce useful functions for medical application in PHITS by referring to examples of various medical applications.
Sato, Tatsuhiko; Furuta, Takuya; Liu, Y.*; Naka, Sadahiro*; Nagamori, Shushi*; Kanai, Yoshikatsu*; Watabe, Tadashi*
EJNMMI Physics (Internet), 8, p.4_1 - 4_16, 2021/01
An individual dosimetry system including the function for calculating EQDX was developed based on PHITS coupled with the microdosimetric kinetic model. It enables us to predict the therapeutic and side effects of TAT based on the clinical data largely available from conventional external radiotherapy.
Satoh, Daiki; Nakayama, Hiromasa; Furuta, Takuya; Yoshihiro, Tamotsu*; Sakamoto, Kensaku
PLOS ONE (Internet), 16(1), p.e0245932_1 - e0245932_26, 2021/01
In this study, we developed a simulation code named SIBYL, which estimates external gamma-ray doses at ground level from radionuclides distributed nonuniformly in atmosphere and on ground. SIBYL can combine with the local-scale atmospheric dispersion model LOHDIM-LES, and calculate the dose distributions according to the map of the activity concentrations simulated by LOHDIM-LES. To apply the SIBYL code to emergency responses of nuclear accidents, the time-consuming three-dimensional radiation transport simulations were performed in advance using the general-purpose Monte Carlo code PHITS, and then the results were compiled to the database for the SIBYL's dose calculations. Moreover, SIBYL can consider the dose attenuation by obstacles and the changes of terrain elevations. To examine the accuracy of SIBYL, typical five cases including Kr emission from a ventilation shaft and Cs dispersion inside urban area were investigated. The results of SIBYL agreed within 10% with those of PHITS at the most of target locations. Furthermore, the calculation speed was approximately 100 times faster than that of PHITS.
Toigawa, Tomohiro; Tsubata, Yasuhiro; Kai, Takeshi; Furuta, Takuya; Kumagai, Yuta; Matsumura, Tatsuro
Solvent Extraction and Ion Exchange, 39(1), p.74 - 89, 2021/00
Absorbed-dose estimation is essential for evaluation of the radiation feasibility of minor-actinide-separation processes. We propose a dose-evaluation method based on radiation permeability, with comparisons of heterogeneous structures seen in the solvent-extraction process, such as emulsions forming in the mixture of the organic and aqueous phases. A demonstration of radiation-energy-transfer simulation is performed with a focus on the minor-actinide-recovery process from high-level liquid waste with the aid of the Monte Carlo radiation-transport code PHITS. The simulation results indicate that the dose absorbed by the extraction solvent from alpha ray depends upon the emulsion structure, and that from beta and gamma ray depends upon the mixer-settler-apparatus size. Non-negligible contributions of well-permeable gamma rays were indicated in terms of the plant operation of the minor-actinide-separation process.
Ratliff, H.; Matsuda, Norihiro; Abe, Shinichiro; Miura, Takamitsu*; Furuta, Takuya; Iwamoto, Yosuke; Sato, Tatsuhiko
Nuclear Instruments and Methods in Physics Research B, 484, p.29 - 41, 2020/12
Carter, L. M.*; Crawford, T. M.*; Sato, Tatsuhiko; Furuta, Takuya; Choi, C.*; Kim, C. H.*; Brown, J. L.*; Bolch, W. E.*; Zanzonico, P. B.*; Lewis, J. S.*
Journal of Nuclear Medicine, 60(12), p.1802 - 1811, 2019/12
Voxel human phantoms have been used for internal dose assessment. More anatomically accurate representation become possible for skins or layer tissues owing to recent developments of advanced polygonal mesh-type phantoms and thus internal dose assessment using those advanced phantoms are desired. However, the Monte Carlo transport calculation by implementing those phantoms require an advanced knowledge for the Monte Carlo transport codes and it is only limited to experts. We therefore developed a tool, PARaDIM, which enables users to conduct internal dose calculation with PHITS easily by themselves. With this tool, a user can select tetrahedral-mesh phantoms, set radionuclides in organs, and execute radiation transport calculation with PHITS. Several test cases of internal dosimetry calculations were presented and usefulness of this tool was demonstrated.
Furuta, Takuya; El Basha, D.*; Iyer, S. S. R.*; Correa Alfonso, C. M.*; Bolch, W. E.*
Journal of Radiological Protection, 39(3), p.825 - 837, 2019/09
Despite large variation of human eye, only one computational eye model has been adopted in almost all the radiation transport simulation studies. We thus adopted a new scalable and deformable eye model and studied the radiation exposure by electrons, photons, and neutrons in the standard radiation fields such as AP, PA, RLAT, ROT, by using Monte Carlo radiation transport code PHITS. We computed the radiation exposure for 5 eye models (standard, large, small, myopic, hyperopic) and analyzed influence of absorbed dose in ocular structures on eye size and shape. Dose distribution of electrons is conformal and therefore the absorbed doses in ocular structures depend on the depth location of each ocular structure. We thus found a significant variation of the absorbed doses for each ocular structure for electron exposure due to change of the depth location affected by eye size and shape. On the other hand only small variation was observed for photons and neutrons exposures because of less conformal dose distribution of those particles.
Yeom, Y. S.*; Han, M. C.*; Choi, C.*; Han, H.*; Shin, B.*; Furuta, Takuya; Kim, C. H.*
Health Physics, 116(5), p.664 - 676, 2019/05
Recently, Task Group 103 of the ICRP developed the mesh-type reference computational phantoms (MCRPs), which are planned for use in future ICRP dose coefficient calculation. Performance of major Monte Carlo particle transport codes (Geant4, MCNP6, and PHITS) were tested with MCRP. External and internal exposure of various particles and energies were calculated and the computational times and required memories were compared. Additionally calculation for voxel-mesh phantom was also conducted so that the influence of different mesh-representation in each code was studied. Memory usage of MRCP was as large as 10 GB with Geant4 and MCNP6 while it is much less with PHITS (1.2 GB). In addition, the computational time required for MRCP tends to increase compared to voxel-mesh phantoms with Geant4 and MCNP6 while it is equal or tends to decrease with PHITS.
ANS RPSD 2018; 20th Topical Meeting of the Radiation Protection and Shielding Division of ANS (CD-ROM), 5 Pages, 2018/08
Recently we introduced a function in PHITS to treat tetrahedral-mesh geometry. Tetrahedral-mesh geometry is a structure composed of combination of tetrahedrons and able to construct complex objects. Tetrahedral-mesh objects can be obtained by converting polygon data using mesh generation software such as TetGen. We also implemented a function in PHITS to export tally results into the format of the three-dimensional visualization software ParaView. TetGen is able to convert the polygon data into ParaView format. Together with these tools, three-dimensional analysis can be realized for PHITS simulation using a polygon objects.
Sato, Tatsuhiko; Iwamoto, Yosuke; Hashimoto, Shintaro; Ogawa, Tatsuhiko; Furuta, Takuya; Abe, Shinichiro; Kai, Takeshi; Tsai, P.-E.; Matsuda, Norihiro; Iwase, Hiroshi*; et al.
Journal of Nuclear Science and Technology, 55(6), p.684 - 690, 2018/06
We have upgraded many features of the Particle and Heavy Ion Transport code System (PHITS) and released the new version as PHITS3.02. The accuracy and the applicable energy ranges of the code were greatly improved and extended, respectively, owing to the revisions to the nuclear reaction models and the incorporation of new atomic interaction models. In addition, several user-supportive functions were developed, such as new tallies to efficiently obtain statistically better results, radioisotope source-generation function, and software tools useful for applying PHITS to medical physics. In this paper, we summarize the basic features of PHITS3.02, especially those of the physics models and the functions implemented after the release of PHITS2.52 in 2013.
Hashimoto, Shintaro; Sato, Tatsuhiko; Iwamoto, Yosuke; Ogawa, Tatsuhiko; Furuta, Takuya; Abe, Shinichiro; Niita, Koji*
Kaku Deta Nyusu (Internet), (120), p.26 - 34, 2018/06
Particle and heavy-ion transport code system PHITS has been used for calculations of radiation shielding in accelerator facilities. PHITS describes physical phenomena induced by radiation as combination of transport and collision processes. The collision process including nuclear reactions is simulated by the three-step calculation: a generation of a reaction, pre-equilibrium, and compound processes. In the simulation, many physics models are used. This report explains roles of the models in PHITS and shows their developments we recently performed.
El Basha, D.*; Furuta, Takuya; Iyer, S. S. R.*; Bolch, W. E.*
Physics in Medicine & Biology, 63(10), p.105017_1 - 105017_13, 2018/05
With recent changes in the recommended annual limit on eye lens exposures to ionizing radiation by International Commission on Radiological Protection, there is considerable interest in predictive computational dosimetry models of the human eye and its various ocular structures. Several computational eye models to date have been constructed for this purpose but they are typically constructed of nominal size and of a roughly spherical shape associated with the emmetropic eye. We therefore constructed a geometric eye model that is both scalable (allowing for changes in eye size) and deformable (allowing for changes in eye shape), and that is suitable for use in radiation transport studies of ocular exposures and radiation treatments of eye disease. As an example, electron and photon anterior-posterior radiation transport with the constructed eye model was conducted and analyzed resultant energy-dependent dose profiles. Due to anterior-posterior irradiation, the energy dose response was shifted to higher energy for a larger-size eye or an axially deformed eye in prolate shape because the structures were located in deeper depth compared to the normal eye.