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. Both condense history and track-structure methods were implemented to handle the electron and positron transport, although the latter is reliable only for simulations in liquid water. In addition, several usersupportive 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. Owing to the continuous improvement and promotion of the code, the number of registered users has exceeded 3,000, and it is being used in diverse areas of study, including accelerator design, radiation shielding and protection, medical physics, and cosmic-ray research. 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.

An upgraded version of the Particle and Heavy Ion Transport code System, PHITS2.52, was developed and released to the public. The new version has been greatly improved from the previously released version, PHITS2.24, in terms of not only the code itself but also the contents of its package, such as the attached data libraries. In the new version, a higher accuracy of simulation was achieved by implementing several latest nuclear reaction models. The reliability of the simulation was improved by modifying both the algorithms for the electron-, positron-, and photon-transport simulations and the procedure for calculating the statistical uncertainties of the tally results. Estimation of the time evolution of radioactivity became feasible by incorporating the activation calculation program DCHAIN-SP into the new package. The efficiency of the simulation was also improved as a result of the implementation of shared-memory parallelization and the optimization of several time-consuming algorithms. Furthermore, a number of new user-support tools and functions that help users to intuitively and effectively perform PHITS simulations were developed and incorporated. Due to these improvements, PHITS is now a more powerful tool for particle transport simulation applicable to various research and development fields, such as nuclear technology, accelerator design, medical physics, and cosmic-ray research.

The nuclear liquid-gas phase transition of the system in ideal thermal equilibrium is studied with antisymmetrized molecular dynamics. The time evolution of a many-nucleon system confined in a container is solved for a long time to get a microcanonical ensemble of a given energy and volume. The temperature and the pressure are extracted from this ensemble and the caloric curves are constructed. The present work is the first time that a microscopic dynamical model which describes nuclear multifragmentation reactions well is directly applied to get the nuclear caloric curve. The obtained constant pressure caloric curves clearly show the characteristic feature of the liquid-gas phase transition, namely negative heat capacity (backbending), which is expected for the phase transition in finite systems.

We performed a benchmark study for 58 cases (22 cases reported in this paper and 36 cases reported in online as supplementary materials of this paper) using the recent version (version 2.88) of the Particle and Heavy-Ion Transport code System (PHITS) in the following fields: (1) particle production cross-sections for nuclear reactions from 20 MeV to 1 GeV, (2) thick-target neutron yields and neutron shielding, (3) depth-dose distribution in water using 12 C beam, and (4) electron and photon transportation over a wide-energy range from keV to GeV. Overall agreements were found to be sufficiently satisfactory; however, several discrepancies are observed, particularly in particle productions with energies below 100 MeV, neutron production for 7 Li(p,n) 7 Be, and photonuclear reactions. To overcome these inaccuracies and to further improve the code, it will be necessary to incorporate a high-energy version of the evaluated nuclear data library JENDL-4.0/HE and the photonuclear data file JENDL-PD in the PHITS package.