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.
Han, M. C.*; Yeom, Y. S.*; Lee, H. S.*; Shin, B.*; Kim, C. H.*; Furuta, Takuya
Physics in Medicine & Biology, 63(9), p.09NT02_1 - 09NT02_9, 2018/05
The multi-threading computation performances of the Geant4, MCNP6, and PHITS codes were evaluated using three tetrahedral-mesh phantoms with different complexity. Photon and neutron transport simulations were conducted and the initialization time, calculation time, and memory usage were measured as a function of the number of threads N used in the simulation. The initialization time significantly increases with the complexity of the phantom, but not much with the number of the threads. For the calculation time, Geant4 showed good parallelization efficiency with multi-thread computation (30 times speed-up factor for N = 40) adopting the private tallies while saturation of the speed-up factor were observed in MCNP6 and PHITS (10 and a few times for N = 40) due to the time delay for the sharing tallies. On the other hand, Geant4 requires larger memory specification and the memory usage rapidly increases with the number of threads compared to MCNP6 or PHITS. It is notable that when compared to the other codes, the memory usage of PHITS is much smaller, regardless of both the complexity of the phantom and the number of the threads.
Nihon Hoken Butsuri Gakkai Homu Peji (Internet), 2 Pages, 2017/10
The current standard of the international radiation protection is determined through life span study of atomic bomb survivors of Hiroshima and Nagasaki, and individual radiation dose estimates for each survivor calculated by Dosimetry System 2002 (DS02). Dose estimates in DS02 are based on input data such as location and shielding condition of the survivor. Therefore accuracy of dose estimates and accordingly that of radiation risk evaluation are largely affected by the accuracy of the input data. In this review, we present an article DS02R1 discussed improvement of the input data of individual atomic bomb survivors for DS02 and its consequence keeping the core system unchanged.
Sato, Tatsuhiko; Niita, Koji*; Iwamoto, Yosuke; Hashimoto, Shintaro; Ogawa, Tatsuhiko; Furuta, Takuya; Abe, Shinichiro; Kai, Takeshi; Matsuda, Norihiro; Okumura, Keisuke; et al.
EPJ Web of Conferences, 153, p.06008_1 - 06008_6, 2017/09
Particle and Heavy Ion Transport code System, PHITS, has been developed under the collaboration of several institutes in Japan and Europe. It can deal with the transport of nearly all particles up to 1 TeV (per nucleon for ion) using various nuclear reaction models and data libraries. More than 2,500 researchers and technicians have used the code for a variety of applications such as accelerator design, radiation shielding and protection, medical physics, and space and geosciences. This paper briefly summarizes physics models and functions newly implemented in PHITS between versions 2.52 and 2.82.