Kuga, Naoya*; Shiiba, Takuro*; Sato, Tatsuhiko; Hashimoto, Shintaro; Kuroiwa, Yasuyoshi*
Journal of Nuclear Science and Technology, 10 Pages, 2023/00
This study aims to verify the accuracy of PHITS in terms of photon and electron transport and provide essential data for its application in clinical dosimetry in high-energy photon beam therapy. Percentage depth dose (PDD), beam profiles, and output factor (OPF) in a water phantom with various field sizes created by a Clinac 21EX linear accelerator were measured using an ionization chamber. Experimental setups were precisely reproduced by PHITS version 3.24, and the percentage differences (%Diff) between the measured and calculated data were evaluated. The average %Diff of PDDs obtained from PHITS and measurement were within 10% and 2% in the build-up and fall-off regions, respectively. For beam profiles, the average %Diff in the plateau region was within 3%; the differences between the calculated and measured distances from the central axis to 50% dose level were within 2 mm. These differences were lower than their tolerance levels. The consistency between the PHITS and EGSnrc was better; their %Diff was within 1% in most cases. The concurrence between the PHITS and measurement shown in this study demonstrates the potential clinical application of PHITS in high-energy photon beam therapy, given its similar dose calculation accuracy compared with EGSnrc.
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.
Walter, H.*; Colonna, M.*; Cozma, D.*; Danielewicz, P.*; Ko, C. M.*; Kumar, R.*; Ono, Akira*; Tsang, M. Y. B*; Xu, J.*; Zhang, Y.-X.*; et al.
Progress in Particle and Nuclear Physics, 125, p.103962_1 - 103962_90, 2022/07
Transport models are the main method to obtain physics information on the nuclear equation of state and in-medium properties of particles from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions to reach consistent conclusions from the same type of physical model. To this end, calculations under controlled conditions of physical input and set-up were performed by the various participating codes. These included both calculations of nuclear matter in a periodic box, which test individual ingredients of a transport code, and calculations of complete collisions of heavy ions. Over the years, five studies were performed within this project. They show, on one hand, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known under the names of Boltzmann-Uehling-Uhlenbeck (BUU) and Quantum Molecular Dynamics (QMD) type codes. On the other hand, there still exist substantial differences when these codes are applied to real heavy-ion collisions. The results of transport simulations of heavy-ion collisions will have more significance if codes demonstrate that they can verify benchmark calculations such as the ones studied in these evaluations.
Iwamoto, Yosuke; Hashimoto, Shintaro; Sato, Tatsuhiko; Matsuda, Norihiro; Kunieda, Satoshi; elik, Y.*; Furutachi, Naoya*; Niita, Koji*
Journal of Nuclear Science and Technology, 59(5), p.665 - 675, 2022/05
A benchmark study of PHITS3.24 has been conducted using neutron-shielding experiments listed in the Shielding Integral Benchmark Archive and Database. Five neutron sources were selected, which are generated from (1) 43- and 68-MeV proton-induced reaction on a thin lithium target, (2) 52-MeV proton-induced reaction on a thick graphite target, (3) 590-MeV proton-induced reaction on a thick lead target, (4) 500-MeV proton-induced reaction on a thick tungsten target, and (5) 800-MeV proton-induced reaction on a thick tantalum target. For all cases, overall agreements in the results are satisfactory when using the JENDL-4.0/HE to simulate neutron- and proton-induced reactions up to 200 MeV. However, discrepancies using PHITS default settings are observed in the results. For an accurate neutron-shielding design for accelerator facilities, using JENDL-4.0/HE in the particle and heavy-ion transport code system calculation is favorable.
Nagai, Yasuki*; Kawabata, Masako*; Hashimoto, Shintaro; Tsukada, Kazuaki; Hashimoto, Kazuyuki*; Motoishi, Shoji*; Saeki, Hideya*; Motomura, Arata*; Minato, Futoshi; Ito, Masatoshi*
Journal of the Physical Society of Japan, 91(4), p.044201_1 - 044201_10, 2022/04
Recently, Lu is considered as one of the most important medical RIs for treating neuroendocrine tumors. A plan to produce Lu with high purity by using enriched Yb samples with irradiation of deuteron beams in accelerators has been discussed. However, since the other Yb isotopes contained in the Yb sample interacts with deuterons, Lu isotopes other than Lu are produced as impurities. Since the purity of Lu is important for medical use, a method to evaluate the impurity of Lu has been required. In this study, we proposed a new method to estimate production yields of each Lu isotopes in Yb samples with arbitrary isotopic compositions by using excitation functions of Yb()Lu reactions and the particle transport calculation code PHITS. The method plays an important role in discussing the isotopic composition of enriched samples to produce high-purity Lu using accelerators.
Kawabata, Masako*; Motoishi, Shoji*; Ota, Akio*; Motomura, Arata*; Saeki, Hideya*; Tsukada, Kazuaki; Hashimoto, Shintaro; Iwamoto, Nobuyuki; Nagai, Yasuki*; Hashimoto, Kazuyuki*
Journal of Radioanalytical and Nuclear Chemistry, 330(3), p.913 - 922, 2021/12
Both Cu and Cu are promising radionuclides in nuclear medicine. Production yields of these radionuclides were quantified by irradiating 55.4 g of natural zinc with accelerator neutrons. Clinically suitable Cu and Cu yields were estimated by experimental based numerical simulations using 100 g of enriched Zn and Zn, respectively, and elevated neutron fluxes from 40 MeV, 2 mA deuterons. A combined thermal- and resin-separation method was developed to isolate Cu and Cu from zinc, resulting in 73% separation efficiency and 97% zinc recovery. Such methods can provide large scale production of Cu and Cu for clinical applications.
Sato, Tatsuhiko; Hashimoto, Shintaro; Inaniwa, Taku*; Takada, Kenta*; Kumada, Hiroaki*
International Journal of Radiation Biology, 97(10), p.1450 - 1460, 2021/10
The stochastic microdosimetric kinetic (SMK) model is one of the most sophisticated and precise models used in the estimation of the relative biological effectiveness of carbon-ion radiotherapy (CRT) and boron neutron capture therapy (BNCT). Through the introduction of Taylor expansion (TE) or fast Fourier transform (FFT), we developed two simplified SMK models and implemented them into the Particle and Heavy Ion Transport code System (PHITS). This study enables the instantaneous calculation of the equieffective dose for CRT and BNCT, considering their cellular-scale dose heterogeneities. Treatment-planning systems that use the improved PHITS as a dose-calculation engine are under development.
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.
Hashimoto, Shintaro; Sato, Tatsuhiko
EPJ Web of Conferences, 239, p.03015_1 - 03015_4, 2020/09
Particle transport simulation codes based on the Monte Carlo technique have been successfully applied to shielding calculations in accelerator facilities. Estimation of not only statistical uncertainties, which depend on the number of trials, but also systemic uncertainties, which are caused by uncertainty of total cross section models, is required to confirm the reliability of the simulation results. We evaluated unclear quantities of internal parameters included in the total cross section model by the KALMAN code, which is based on the least squares technique, comparing with experimental data of the total cross section. The uncertainties in the total cross sections obtained by the new model are comparable to the experimental errors. In the present study, the systematic uncertainty included in the simulation results can be estimated by performing the transport calculations with variation of the internal parameters within their unclear quantities.
Iwamoto, Yosuke; Meigo, Shinichiro; Hashimoto, Shintaro
Journal of Nuclear Materials, 538, p.152261_1 - 152261_9, 2020/09
The displacements-per-atom (dpa) is widely used as an exposure unit to predict the operating lifetime of materials in radiation environments. Because the athermal-recombination-corrected dpa (arc-dpa) model is a more realistic model than the standard Norgertt-Robinson-Torrens (NRT) model, new evaluation of radiation damage will be performed using the arc-dpa model as a standard. In this work, the recent arc-dpa model of various materials are incorporated in PHITS, and the rescaling factors (NRT-dpa/arc-dpa) over a wide energy range are reported. For neutron incidences with the energy spectrum determined in selected nuclear facilities and proton incidences with energies of 600 MeV-50 GeV, the rescaling factor for each material is independent of these irradiation conditions with almost the same value for each material. Our findings will be beneficial for rescaling the NRT-dpa model used for radiation damage applications over a wide energy region.
Tsukada, Kazuaki; Nagai, Yasuki*; Hashimoto, Shintaro; Minato, Futoshi; Kawabata, Masako*; Hatsukawa, Yuichi*; Hashimoto, Kazuyuki*; Watanabe, Satoshi*; Saeki, Hideya*; Motoishi, Shoji*
Journal of the Physical Society of Japan, 89(3), p.034201_1 - 034201_7, 2020/03
We found anomalously large yields of Ga, Ga, Zn and Cu by neutron irradiation on a ZnO sample in a polyethylene shield. Neutron beams are generated from the Be() reaction for 50 MeV deuterons. The yields obtained were more than 20 times larger than those in the unshielded sample. On the other hand, the yields of Ga, Ga, Zn and Cu from a metallic Zn sample and the yields of Cu, Ni and Zn from the ZnO and Zn samples were almost insensitive to the shield conditions. This finding would provide us a unique capability of accelerator neutrons to simultaneously produce a large amount of several radioisotopes, including proton induced reaction products, by using a single sample. The experimental data were compared with the yields estimated by using the Particle and Heavy Ion Transport code System and the result was discussed.
Nakayama, Shinsuke; Iwamoto, Osamu; Iwamoto, Nobuyuki; Hashimoto, Shintaro
Kaku Deta Nyusu (Internet), (123), p.53 - 59, 2019/06
The 2019 International Conference on Nuclear Data for Science and Technology (ND2019) was held at the China National Convention Center on May 19-24, 2019. The series of the ND conferences are the largest conferences in nuclear data research field that are held every three years. In this paper, as a part of the conference reports of ND2019, the authors gave summaries of the presentations on nuclear data evaluation and theory conducted at the conference.
Hashimoto, Shintaro; Sato, Tatsuhiko
Journal of Nuclear Science and Technology, 56(4), p.345 - 354, 2019/04
Particle transport simulations based on the Monte Carlo method have been applied to shielding calculations. Estimation of not only statistical uncertainty related to the number of trials but also systematic one induced by unclear physical quantities is required to confirm the reliability of calculated results. In this study, we applied a method based on analysis of variance to shielding calculations. We proposed random- and three-condition methods. The first one determines randomly the value of the unclear quantity, while the second one uses only three values: the default value, upper and lower limits. The systematic uncertainty can be estimated adequately by the random-condition method, though it needs the large computational cost. The three-condition method can provide almost the same estimate as the random-condition method when the effect of the variation is monotonic. We found criterion to confirm convergence of the systematic uncertainty as the number of trials increases.
Ogawa, Tatsuhiko; Sato, Tatsuhiko; Hashimoto, Shintaro; Niita, Koji*
Physical Review C, 98(2), p.024611_1 - 024611_15, 2018/08
Particle production by nucleus-nucleus reactions in the energy range from GeV to TeV is substantially important for safety evaluation in heavy ion accelerators and evaluation of space radiation dose. A lot of models and theories have been studied. In the models developed in the past, interaction between nucleons were dependent on the reference frame; therefore the moving incident nucleus and the target nucleus at rest transferred to the common frame were disintegrated. Previously, intentional bias was introduced to the calculation algorithms to supplement stability but residual nucleus mass and secondary particle production was underestimated. In this study, a reaction model JAMQMD was developed, in which intra-nucleon interaction was described in a frame-independent way. This model can reproduce the stability of nuclei regardless of the reference frame and the yield of residual nuclei as well as secondary particles including deuterons. JQMD Ver.2 developed 3 years ago can simulate nucleus-nucleus reactions up to 3 GeV/nucleon; therefore the development of JAMQMD is the doorway to simulate nucleus-nucleus reactions regardless of the incident energy. JAMQMD is an useful model for not only radiation protection studies but also analysis of fundamental physics studies.
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.
Ogawa, Tatsuhiko; Hashimoto, Shintaro; Sato, Tatsuhiko; Niita, Koji*; Kamae, Tsuneyoshi*
Proceedings of 3rd International Conference on Particle Physics and Astrophysics (ICPPA 2017) (Internet), p.391 - 398, 2018/04
In the reactions of cosmic-ray heavy ions producing high-energy gamma-rays, heavy ions are directed to the ground therefore prompt gamma-rays from projectile fragments are boosted by Doppler effect and observed on the ground. Recently, a lot of experiments pay attention to such high energy gamma-rays. In order to simulate such high-energy gammas, event-by-event simulation of fragmentation reaction which determines the excitation energy and angular momentum of the produced fragment, and de-excitation simulation based on the nuclear structure data. Such models are available in the general-purpose radiation transport simulation code PHITS. JAERI Quantum Molecular Dynamics model was recently updated to accurately simulate charge and mass distribution of fragments. EBITEM, which was released recently, can simulate gamma deexcitation after evaporation based on the excitation energy and angular momentum. Thus latest PHITS can accurately simulate production of -rays attributed to cosmic-ray heavy ions. This study shows a new approach to reproduce -rays by cosmic-ray heavy ions.
Minato, Futoshi; Tsukada, Kazuaki; Sato, Nozomi*; Watanabe, Satoshi*; Saeki, Hideya*; Kawabata, Masako*; Hashimoto, Shintaro; Nagai, Yasuki*
Journal of the Physical Society of Japan, 86(11), p.114803_1 - 114803_6, 2017/11
We have measured the yield of Mo, the mother nuclide of Tc used in nuclear diagnostic procedure. Mo was produced by Mo(,)Mo using neutrons with thermal energy up to about 40 MeV, provided by C(,). The Mo yield agrees with an estimated yield with the use of the latest data of C(,) and the evaluated cross section given in the JENDL. Based on this, a new calculation was carried out to produce Mo to seek for a good economical condition. Various conditions such as the MoO sample mass, the distance between the carbon target and the sample, the radius of the deuteron beam, and the neutron irradiation time were considered. The calculated Mo yield indicates that about 30% of the Mo demand in Japan can be fulfilled with a single accelerator. The elusion of Tc from the Mo twice per day would meet about 50% of the Mo demand.
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.