Sublet, J.-Ch.*; Bondarenko, I. P.*; Bonny, G.*; Conlin, J. L.*; Gilbert, M. R.*; Greenwood, L. R.*; Griffin, P. J.*; Helgesson, P.*; Iwamoto, Yosuke; Khryachkov, V. A.*; et al.
European Physical Journal Plus (Internet), 134(7), p.350_1 - 350_50, 2019/07
Nuclear reaction with nuclear data is the origin of defects produced by cascade damage in irradiated materials. Therefore, it is important to consider nuclear reaction correctly for calculations of the damage energy of Primary Knock on Atom (PKA) and the number of Displacement Per Atom (DPA). Here, radiation damage metrics considering nuclear reaction enables us to simulate transport of each defect and clustering defects in the irradiated material. This paper reviews the theory of nuclear reaction and damage energy and describes the latest methodologies about uncertainty propagation and quantification in nuclear data and damage calculations based on molecular dynamics.
Panikkath, P.*; Otsuka, Naohiko*; Iwamoto, Yosuke; Mohanakrishnan, P.*
European Physical Journal A, 55(6), p.91_1 - 91_9, 2019/06
The thermal neutron capture cross section of Ga was evaluated to be around 4.7 b in many nuclear data libraries. On the other hand, it was evaluated to be 3.71 b in JENDL-4.0. Since these cross section data were obtained by the activation measurement with a monitor foil to determine the absolute value of cross section, data depend on the monitor foil adopted in the measurement. In this work, we measured the thermal neutron capture cross section of Ga with a few monitor foils. In the experiment, neutrons decelerated in concrete shield were irradiated with a sample with Ga foil and Au and Mn monitor foils. A covariance analysis was performed to obtain the off-diagonal weighted mean of the results determined with the two monitor reactions. As a result, the thermal neutron capture cross section of Ga was 4.050.27 b and is close to the value of JENDL-4.0.
Tsai, P.-E.; Iwamoto, Yosuke; Hagiwara, Masayuki*; Sato, Tatsuhiko; Ogawa, Tatsuhiko; Satoh, Daiki; Abe, Shinichiro; Ito, Masatoshi*; Watabe, Hiroshi*
Proceedings of 2017 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC 2017) (Internet), 3 Pages, 2018/11
The energy spectra of primary knock-on atoms (PKAs) are essential for radiation damage assessment in design of accelerator facilities. However up to date the experimental data are still limited, due to the poor mass resolution and the high measurement threshold energies in the conventional setup of nuclear physics experiments using solid state detectors, which are typically above a few MeV/nucleon. In this study, a novel detection system consisting of two time detectors and one dE-E energy detector is proposed and being constructed to measure the PKA spectra. The system and detector design was based on Monte Carlo simulations by using the PHITS code. The PHITS simulations show that the system is able to distinguish the PKA isotopes above 0.2-0.3 MeV/nucleon for A=2030 amu; the PKA mass identification thresholds decrease to 0.1 MeV/nucleon for PKAs lighter than 20 amu. The detection system will be tested in the summer of 2017, and the test results will be presented at the conference.
Iwamoto, Yosuke; Yoshida, Makoto*; Yoshiie, Toshimasa*; Satoh, Daiki; Yashima, Hiroshi*; Matsuda, Hiroki; Meigo, Shinichiro; Shima, Tatsushi*
Journal of Nuclear Materials, 508, p.195 - 202, 2018/09
To validate the displacement damage model in radiation transport codes used for the estimation of radiation damages at accelerator facilities, we measured electrical resistance increase of aluminum and copper induced by radiation defects under the cryogenic 200 MeV proton irradiation. The irradiation device had the structure to cool two irradiation samples at same time using thermal conductance. The aluminum and copper wire with 250 m diameter was sandwiched between two AlN plates with excellent thermal conductivity and electrical insulation. As a result, temperature of irradiation samples was kept at below 5 K under proton irradiation with beam intensity below 3 nA. The experimental displacement cross section agreed with calculated results with defect production efficiency.
Iwamoto, Yosuke; Matsuda, Hiroki; Meigo, Shinichiro; Satoh, Daiki; Nakamoto, Tatsushi*; Yoshida, Makoto*; Ishi, Yoshihiro*; Kuriyama, Yasutoshi*; Uesugi, Tomonori*; Yashima, Hiroshi*; et al.
Proceedings of 61st ICFA Advanced Beam Dynamics Workshop on High-Intensity and High-Brightness Hadron Beams (HB 2018) (Internet), p.116 - 121, 2018/07
The radiation damage model in the radiation transport code PHITS has been developed to calculate the basic data of the radiation damage including the energy of the target Primary Knock on Atom (PKA). For the high-energy proton incident reactions, a target PKA created by the secondary particles was more dominant than a target PKA created by the projectile. To validate the radiation damage model in metals irradiated by 100 MeV protons, we developed a proton irradiation device with a Gifford-McMahon cryocooler to cryogenically cool wire samples. By using this device, the defect-induced electrical resistivity changes related to the DPA cross section of copper and aluminum were measured under irradiation with 125 and 200 MeV protons at cryogenic temperature. A comparison of the experimental data with the calculated results indicates that the DPA cross section with defect production efficiencies provide better quantitative descriptions.
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.
Meigo, Shinichiro; Matsuda, Hiroki; Iwamoto, Yosuke; Iwamoto, Hiroki; Hasegawa, Shoichi; Maekawa, Fujio; Yoshida, Makoto*; Ishida, Taku*; Makimura, Shunsuke*; Nakamoto, Tatsushi*
Proceedings of 9th International Particle Accelerator Conference (IPAC '18) (Internet), p.499 - 501, 2018/06
no abstracts in English
Shigyo, Nobuhiro*; Iwase, Hiroshi*; Iwamoto, Yosuke; Sato, Tatsuhiko
Nippon Genshiryoku Gakkai-Shi, 60(5), p.294 - 298, 2018/05
For calculations of neutron productions and radiation damages in nuclear design of Accelerator Driven System (ADS) and exposure doses in radiation cancer therapy, radiation transport simulations with the high-energy nuclear data library evaluated with experimental data play an important role. Experimental data are needed for validation of nuclear reaction models in calculation codes. In this paper, we explain examples of nuclear data measurements in high energy region and the progress of the Particle and Heavy-Ion Transport code System (PHITS) developed in Japan.
Nuclear Instruments and Methods in Physics Research B, 419, p.32 - 37, 2018/03
In this study, the Monte Carlo displacement damage calculation method in the Particle and Heavy-Ion Transport code System (PHITS) was improved to calculate displacements per atom (DPA) due to irradiation by electrons and rays. For the damage due to electrons and rays, PHITS simulates electromagnetic cascades using the EGS5 algorithm and calculates DPA values using the recoil energies. A comparison of DPA values calculated by PHITS and the Monte Carlo assisted Classical Method (MCCM) reveals that they were in good agreement for -ray irradiations of silicon and iron at energies that were less than 10 MeV. Above 10 MeV, PHITS can calculate DPA values not only for electrons but also for charged particles produced by photonuclear reactions. For irradiation of 90-cm-thick carbon by protons with energies of more than 30 GeV, the ratio of the secondary electron DPA values to the total DPA values is more than 10% and increases with an increase in incident energy.
Araki, Shohei*; Watanabe, Yukinobu*; Kitajima, Mizuki*; Sadamatsu, Hiroki*; Nakano, Keita*; Kin, Tadahiro*; Iwamoto, Yosuke; Satoh, Daiki; Hagiwara, Masayuki*; Yashima, Hiroshi*; et al.
EPJ Web of Conferences (Internet), 146, p.11027_1 - 11027_4, 2017/09
In recently years, deuteron-induced reaction is considered to produce the neutron source for application fields such as radiation damage fusion materials and boron neutron capture therapy. However, as the experimental data are not sufficient at incident energies above 60 MeV, the theoretical models are not validated. Therefore, we measured the double differential cross sections (DDXs) for Li, Be and C at 100 MeV at the Research Center for Nuclear Physics in Osaka University. The DDXs were measured at 6 angles (025 and neutron energy was determined by a time of flight method. Three different-size NE213 liquid organic scintillators located at a distance of 7 m, 24 m and 74 m respectively were adopted as neutron detectors. In the measured DDXs, a broad peak due to deuteron breakup process was observed at approximately half of the deuteron incident energy. The DDXs calculated by PHITS did not reproduce the experimental ones due to lack of theoretical model.
Iwamoto, Yosuke; Sato, Tatsuhiko; Satoh, Daiki; Hagiwara, Masayuki*; Yashima, Hiroshi*; Masuda, Akihiko*; Matsumoto, Tetsuro*; Iwase, Hiroshi*; Shima, Tatsushi*; Nakamura, Takashi*
EPJ Web of Conferences (Internet), 153, p.08019_1 - 08019_3, 2017/09
To develop 100-400 MeV quasi-monoenergetic neutron field, we measured neutron and unexpected -ray energy spectra of the Li(p,n) reaction with 80-389 MeV protons in the 100-m time-of-flight (TOF) tunnel at the Research Center for Nuclear Physics (RCNP). Neutron energy spectra with energies above 3 MeV were measured by the TOF method and energy spectra with energies above 0.1 MeV were measured by the automatic unfolding function of the radiation dose monitor DARWIN. For neutron spectra, the contribution of peak intensity to the total intensity integrated with energies above 3 MeV varied between 0.38 and 0.48. For -ray spectra, high-energetic -rays at around 70 MeV originated from the decay of were observed over 200 MeV. For the 246-MeV proton incident reaction, the contribution of -ray dose to neutron dose is negligible because the ratio of -ray to neutron is 0.014.
Theis, C.*; Carbonez, P.*; Feldbaumer, E.*; Forkel-Wirth, D.*; Jaegerhofer, L.*; Pangallo, M.*; Perrin, D.*; Urscheler, C.*; Roesler, S.*; Vincke, H.*; et al.
EPJ Web of Conferences (Internet), 153, p.08018_1 - 08018_5, 2017/09
At CERN, gas-filled ionization chambers PTW-34031 (PMI) are commonly used in radiation fields including neutrons, protons and -rays. A response function for each particle is calculated by the radiation transport code FLUKA. To validate a response function to high energy neutrons, benchmark experiments with quasi mono-energetic neutrons have been carried out at RCNP, Osaka University. For neutron irradiation with energies below 200 MeV, very good agreement was found comparing the FLUKA simulations and the measurements. In addition it was found that at proton energies of 250 and 392 MeV, results calculated with neutron sources underestimate the experimental data due to a non-negligible gamma component originating from the target Li(p,n)Be reaction.
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.
Matsumoto, Tetsuro*; Masuda, Akihiko*; Nishiyama, Jun*; Iwase, Hiroshi*; Iwamoto, Yosuke; Satoh, Daiki; Hagiwara, Masayuki*; Yashima, Hiroshi*; Yashima, Hiroshi*; Shima, Tatsushi*; et al.
EPJ Web of Conferences (Internet), 153, p.08016_1 - 08016_3, 2017/09
Neutron energy spectra behind concrete and iron shields were measured for quasi-monoenergetic neutrons above 200 MeV using a Bonner sphere spectrometer (BSS). Quasi-monoenergetic neutrons were produced by the Li(p,xn) reaction with 246-MeV and 389-MeV protons. The response function of BSS was also measured at neutron energies from 100 MeV to 387 MeV. In data analysis, the measured response function was used and the multiple neutron scattering effect between the BSS and the shielding material was considered. The neutron energy spectra behind the concrete and iron shields were obtained by the unfolding method using the MAXED code. Ambient dose equivalents were obtained as a function of a shield thickness successfully. For the case of the 244 MeV neutron incidence, the multiple neutron scattering effect on the effective dose is large under 50 cm thickness of the concrete shield.
Mares, V.*; Trinkl, S.*; Iwamoto, Yosuke; Masuda, Akihiko*; Matsumoto, Tetsuro*; Hagiwara, Masayuki*; Satoh, Daiki; Yashima, Hiroshi*; Shima, Tatsushi*; Nakamura, Takashi*
EPJ Web of Conferences (Internet), 153, p.08020_1 - 08020_3, 2017/09
To validate response of an extended range Bonner Sphere Spectrometer (ERBSS) with He proportional counter, neutron energy spectra were measured using an ERBSS in the quasi-mono-energetic neutron field at the Research Center for Nuclear Physics (RCNP). Using 100 MeV and 296 MeV proton beams, neutron fields with nominal peak energies of 96 MeV and 293 MeV were generated via Li(p,n)Be reactions. The energy spectra were measured at a distance of 35 m from the target. To deduce the corresponding neutron spectra from thermal to the nominal maximum energy, the ERBSS data were unfolded using the MSANDB unfolding code. At high energies, the neutron spectra were also measured by means of the TOF method using NE213 organic liquid scintillators. The agreement between ERBSS and TOF neutron spectra above 5 MeV is very good. Comparison in terms of ambient dose equivalent, H(10) between ERBSS and TOF values for both proton energies shows very good agreement.
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.
Iwamoto, Yosuke; Ogawa, Tatsuhiko
Nuclear Instruments and Methods in Physics Research B, 396, p.26 - 33, 2017/04
It is important to validate the calculations of PKA spectra that are used to estimate radiation damage in materials. Here, the PKA spectra of fission-relevant materials were calculated using the Particle and Heavy Ion Transport code System (PHITS) and also using the data processing code NJOY with the nuclear data libraries. The heating number, which is the integral of the PKA spectrum, was also calculated using PHITS and compared with data extracted from the data libralies. From analyzing the PKA spectra, we found that the energy and angular recoil distributions were incorrect for Ge, As, Y, and Ag in the ENDF/B-VII.1 library. From analyzing the heating number, we found that the data extracted from the ACE file of TENDL2015 for all elements were problematic in the neutron capture region. However, PHITS+TENDL2015 can calculate PKA spectra and heating numbers correctly.
Masuda, Akihiko*; Matsumoto, Tetsuro*; Iwamoto, Yosuke; Hagiwara, Masayuki*; Satoh, Daiki; Sato, Tatsuhiko; Iwase, Hiroshi*; Yashima, Hiroshi*; Nakane, Yoshihiro; Nishiyama, Jun*; et al.
Nuclear Instruments and Methods in Physics Research A, 849, p.94 - 101, 2017/03
Quasi-monoenergetic high-energy neutron fields induced by Li(p,n) reactions are used for the response evaluation of neutron-sensitive devices. The quasi-monoenergetic high-energy field consists of high-energy monoenergetic peak neutrons and unwanted continuum neutrons down to the low-energy region. A two-angle differential method has been developed to compensate for the effect of the continuum neutrons in the response measurements. In this study, the two-angle differential method was demonstrated for Bonner sphere detectors, which are typical examples of moderator-based neutron-sensitive detectors, to investigate the method's applicability and its dependence on detector characteristics. Through this study, the adequacy of the two-angle differential method was experimentally verified, and practical suggestions were made pertaining to this method.