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.
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 and 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 and 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 (Internet), 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 and 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.
Furuta, Takuya; Maeyama, Takuya*; Ishikawa, Kenichi*; Fukunishi, Nobuhisa*; Fukasaku, Kazuaki*; Takagi, Shu*; Noda, Shigeho*; Himeno, Ryutaro*; Hayashi, Shinichiro*
Physics in Medicine and Biology, 60(16), p.6531 - 6546, 2015/08
Low reproducibility of dose distribution in inhomogeneous regions such as soft matter near bones is known with the simple dose analysis currently adopted in treatment planning of particle cancer therapy. Therefore a treatment planning system based on Monte Carlo simulation having better accuracy is highly desired. In order to assess the simulation accuracy of a Monte Carlo simulation code in situations closely related to medical application, we performed a comparison of dose distribution in a biological sample obtained by experiment and that by simulation. In particular, we irradiate a carbon beam on a biological sample composed of fresh chicken meat and bones, with a PAGAT gel dosimeter placed behind it, and compare the complex dose distribution in the gel dosimeter created by the beam passing through the inhomogeneous sample. Monte Carlo simulation using PHITS code was conducted by reconstructing the biological sample from its computed tomography images. The simulation accurately reproduced the experimental distal edge structure of the dose distribution with an accuracy under about 2 mm.
Furuta, Takuya; Takahashi, Fumiaki
Journal of Nuclear Science and Technology, 52(6), p.897 - 904, 2015/06
Dependence of the radiation dose reduction on sizes and materials of buildings was studied by numerical analyses using the Monte Carlo simulation code PHITS. Three building models were developed: the wooden house, open space concrete building, and thin-wall building, to study the effect of building size and construction material on dose reduction inside these structures. Here the floor area size of the building models were varied to clarify the influence of building configuration on dose reduction. The results demonstrated that the dose rates inside the buildings linearly decreased with increasing floor area on a logarithmic scale for all types of buildings. The calculation of dose distribution inside a building clarified that the distance from outer walls were a determinant factor for the dose rate at each position in the building. The obtained tendency was ensured by comparing to data for the dose reduction of typical buildings in Japan.