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.
Nippon 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.
Satoh, Daiki; Furuta, Takuya; Takahashi, Fumiaki; Lee, C.*; Bolch, W. E.*
Journal of Nuclear Science and Technology, 54(9), p.1018 - 1027, 2017/09
The personal dose equivalent was calculated for the public (newborns; 1-, 5-, 10-, and 15-year-old children; and adults) in an environment contaminated with radioactive cesium (Cs and Cs) distributed in a soil at specific depths of 0.0, 0.5, 2.5, 5.0, 10.0, and 50.0 g/cm. Monte Carlo calculations were performed using pediatric and adult computational phantoms incorporated into a particle and heavy ion transport code system (PHITS). Compared with the effective dose and ambient dose equivalent at a height of 100 cm above the ground, the personal dose equivalent was found to provide an acceptable assessment for the effective dose and did not exceed the ambient dose equivalent in the environmental radiation field, while the personal dose equivalent values increased for younger subjects. The weighted-integral method to obtain the personal dose equivalent for a volumetric source was applied to the analysis of exponential radioactive cesium distributions in the soil observed in Fukushima, and the calculation results successfully reproduced the measured data.
Furuta, Takuya; Sato, Tatsuhiko; Han, M. C.*; Yeom, Y. S.*; Kim, C. H.*; Brown, J. L.*; Bolch, W. E.*
Physics in Medicine & Biology, 62(12), p.4798 - 4810, 2017/06
A new function to treat tetrahedral-mesh geometry, a type of polygon-mesh geometry, was implemented in the Particle and Heavy Ion Transport code Systems (PHITS). Tetrahedral-mesh is suitable to describe complex geometry including curving shapes. In addition, construction of three-dimensional geometry using CAD software becomes possible with file format conversion. We have introduced a function to create decomposition maps of tetrahedral-mesh objects at the initial process so that the computational time for transport process can be reduced. Owing to this function, transport calculation in tetrahedral-mesh geometry can be as fast as that for the geometry in voxel-mesh with the same number of meshes. Due to adaptability of tetrahedrons in size and shape, dosimetrically equivalent objects can be represented by tetrahedrons with much fewer number of meshes compared with the voxels. For dosimetric calculation using computational human phantom, significant acceleration of the computational speed, about 4 times, was confirmed by adopting the tetrahedral mesh instead of the voxel.
Sato, Tatsuhiko; Iwamoto, Yosuke; Hashimoto, Shintaro; Ogawa, Tatsuhiko; Furuta, Takuya; Abe, Shinichiro; Kai, Takeshi; Matsuda, Norihiro; Iwase, Hiroshi*; Niita, Koji*
Hoshasen, 43(2), p.55 - 58, 2017/05
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 registered researchers and technicians have used this system for various applications such as accelerator design, radiation shielding and protection, medical physics, and space- and geo-sciences. This paper summarizes the physics models and functions recently implemented in PHITS, between versions 2.52 and 2.88.
Iwamoto, Yosuke; Sato, Tatsuhiko; Hashimoto, Shintaro; Ogawa, Tatsuhiko; Furuta, Takuya; Abe, Shinichiro; Kai, Takeshi; Matsuda, Norihiro; Hosoyamada, Ryuji*; Niita, Koji*
Journal of Nuclear Science and Technology, 54(5), p.617 - 635, 2017/05
We performed a benchmark study for 58 cases using the recent version 2.88 of the Particle and Heavy Ion Transport code System (PHITS) in the following fields: particle production cross-sections for nuclear reactions, neutron transport calculations, and electro-magnetic cascade. This paper reports details for 22 cases. In cases of nuclear reactions with energies above 100 MeV and electro-magnetic cascade, overall agreements were found to be satisfactory. On the other hand, PHITS did not reproduce the experimental data for an incident proton energy below 100 MeV, because the intranuclear cascade model INCL4.6 in PHITS is not suitable for the low-energy region. For proton incident reactions over 100 MeV, PHITS did not reproduce fission product yields due to the problem of high-energy fission process in the evaporation model GEM. To overcome these inaccuracies, we are planning to incorporate a high-energy version of the evaluated nuclear data library JENDL-4.0/HE, and so on.
Igaku Butsuri, 37(3), p.177 - 180, 2017/00
Gel dosimeters are gelated aqueous solutions of radiosensitive substances and useful tools to visualize three-dimensional dose distribution of radiations. A study which uses gel dosimeters to quantify the dose distribution in biological samples is explained. In this study, the measured distribution was compared with the simulation result of PHITS to check the accuracy of Monte Carlo simulation in particle therapy situations. The comparison shows that the measured complex distal edge structure of the dose distribution reflecting the particle transport in inhomogeneous biological sample was reproduced by simulation with a precision of 2 mm difference. The result also shows that gel dosimeters are powerful tools to visualize three-dimensional dose distribution and have potential to be used for various purposes such as quality assurance of treatment beams as well as accuracy check of simulations.
Iwamoto, Yosuke; Sato, Tatsuhiko; Niita, Koji*; Hashimoto, Shintaro; Ogawa, Tatsuhiko; Furuta, Takuya; Abe, Shinichiro; Kai, Takeshi; Matsuda, Norihiro; Iwase, Hiroshi*; et al.
JAEA-Conf 2016-004, p.63 - 69, 2016/09
A general purpose Monte Carlo Particle and Heavy Ion Transport code System, PHITS, is being developed through the collaboration of several institutes. PHITS can deal with the transport of nearly all particles, including neutrons, protons, heavy ions, photons, and electrons, over wide energy ranges using various nuclear reaction models and data libraries. PHITS users apply the code to various research and development fields such as nuclear technology, accelerator design, medical physics, and cosmic-ray research. This presentation briefly summarizes the physics models implemented in PHITS, and introduces some new models such as muon-induced nuclear reaction model and a de-excitation model EBITEM. We will also present the radiation damage cross sections for materials, PKA spectra and kerma factors calculated by PHITS under the IAEA-CRP activity titled "Primary radiation damage cross section."
Satoh, Daiki; Furuta, Takuya; Takahashi, Fumiaki; Endo, Akira; Lee, C.*; Bolch, W. E.*
Journal of Nuclear Science and Technology, 53(1), p.69 - 81, 2016/01
To estimate effective doses for the public exposed to external radiation from radioactive cesium (Cs and Cs) deposited on the ground by the Fukushima nuclear accident, we calculate the conversion coefficients for converting activity concentration to effective dose rate by using the Particle and Heavy Ion Transport code System (PHITS). The data were produced from different age groups within the public (newborns; 1-, 5-, 10-, and 15-year-old children; and adults) for the situations in which radioactive cesium is distributed uniformly in the soil over a planar area and at specific depths of 0.0, 0.5, 2.5, 5.0, 10.0, and 50.0 g/cm. On the basis of the results, we also derive the conversion coefficients for exponentially distributed volumetric sources. In addition, we obtain the conversion coefficients that give the effective dose accumulated over the first and second months, the first year, and over a lifetime (50 years) because of the contamination remaining on the ground. These calculations indicate that the conversion coefficients to obtain the effective dose rate are higher for the younger ages compared with adults but do not exceed the ambient dose equivalent rate. Furthermore, we find that the difference between the calculated effective dose rates according to the International Commission on Radiological Protection (ICRP) 1990 and 2007 Recommendations is small (7% maximum) for a ground contamination of radioactive cesium.
Furuta, Takuya; Hashimoto, Shintaro; Sato, Tatsuhiko
Igaku Butsuri, 36(1), p.50 - 54, 2016/00
An application of Particle and Heavy Ion Transport Code System; PHITS, for medical physics is to simulate treatment planning of radiation therapy. Treatment planning simulation is conducted by constructing patient geometry from patient CT data, calculating radiation transport of external beam, and deducing dose distribution inside patient body. However, it is not easy to extract information such as patient location and CT value distribution from patient CT data or to construct complex accelerator geometry in PHITS format. Therefore, we developed two user assistance programs, DICOM2PHITS and PSFC4PHITS. DICOM2PHITS is a program to construct the voxel PHITS simulation geometry from patient CT DICOM image data. PSFC4PHITS is a program to convert the IAEA phase-space file data to PHITS format to be used as simulation source of PHITS. We explain these two programs by showing some applications in this article.
Furuta, Takuya; Takahashi, Fumiaki
Radiation Protection Dosimetry, 167(4), p.392 - 398, 2015/12
A new approach was proposed to compute radiation dose to objects of interest in infinitely expanded radiation fields with Monte Carlo transport codes. The particles, which were emitted from the location far from the interested object, could be effectively taken into account by setting reflection boundaries at the borders of the computational area in the newly developed method. Here, the positions and momenta of the particles at each reflection were recorded and used as the sources to simulate particles from the exterior region of the computational area. The validity of new method was checked by radiation transport calculations of objects on the infinitely expanded contaminated ground surface. The results in our new approach agreed within statistical errors of the simulation with those in the conventional approach considering a very large computational area. In addition, an efficiency in computational time was compared between the new and conventional approaches.
Sato, Tatsuhiko; Niita, Koji*; Matsuda, Norihiro; Hashimoto, Shintaro; Iwamoto, Yosuke; Furuta, Takuya; Noda, Shusaku; Ogawa, Tatsuhiko; Iwase, Hiroshi*; Nakashima, Hiroshi; et al.
Annals of Nuclear Energy, 82, p.110 - 115, 2015/08
The general purpose Monte Carlo Particle and Heavy Ion Transport code System, PHITS, is being developed through a collaboration of several institutes in Japan and Europe. The Japan Atomic Energy Agency is responsible for managing the entire project. PHITS can deal with the transport of nearly all particles, including neutrons, protons, heavy ions, photons, and electrons, over wide energy ranges using various nuclear reaction models and data libraries. This paper briefly summarizes the physics models implemented in PHITS, and introduces some important functions useful for particular purposes, such as an event generator mode and beam transport functions.