Physical Models > Map of Models


     Physical Models included in PHITS

Note: The followings are description for PHITS2.64!

Above figure summarizes the physics models recommended for use in PHITS for simulating nuclear and atomic collisions. The intra-nuclear cascade models JAM [1] and INCL4.6 [2] are employed for simulating the dynamic stage of hadron-induced nuclear reactions in the high and intermediate energy regions, respectively. Other intra-nuclear cascade models, such as INC-ELF [3] and a modified version of BERTINI [4], can be also selected in this region. On the other hand, the quantum molecular dynamics model JQMD [5] is generally employed for nucleus-induced reactions, because consideration of the nuclear forces between every combination of nucleons is very important in the simulation. Simulations for deuteron-, triton-, 3He-, and -induced reactions at intermediate energies are the exceptional cases where INCL4.6 is recommended to be employed, because INCL4.6 can deal with light-ion induced reactions and is much faster than JQMD. The evaporation and fission model GEM [6] is adopted for simulating the static stage for both hadron- and nucleus-induced reactions. The energy losses of charged particles, except for electrons, are calculated using the SPAR [7] or ATIMA [8] codes with the continuous slowing down approximation.

Nuclear and atomic data libraries are generally used for simulating low-energy neutron-induced nuclear reactions and photo- and electro-atomic interactions, respectively. These data libraries are written in the ACE format, which is the same format as that adopted in MCNP [9], but different from that of original evaluated-data files, such as JENDL-4.0 [10]. Thus, it is necessary to develop PHITS-readable data libraries by compiling the original evaluated-data files using nuclear data processing systems such as NJOY99 [11]. The neutron and photon data libraries included in the PHITS package were compiled from JENDL-4.0 and its sub-library based on EPDL97 [12], respectively. The electron library was developed based on EEDL [13] and ITS3.0 [14], in addition to EPDL97. Photo-nuclear reactions can be treated only for the lower energy regions in which the giant resonance is the dominant reaction mechanism.

The highest energies of particles recommended to be simulated using PHITS are 200 GeV for hadrons, 100 GeV/n for nucleus, 100 MeV for electrons and positrons, and 100 GeV for photons. Note that PHITS can deal with the transport of particles above these energies, but the accuracy of simulations for such high-energy particles has not been validated. On the other hand, the lowest particle energy that can be simulated using PHITS is 1 keV, except for neutrons, because the continuous slowing down approximation cannot be applied to the transport simulation of charged particles below this energy. Event-by-event analysis of each ionization is necessary for such low-energy particle transport simulation, but this process is not incorporated in PHITS because of its long computational time. Thus, the minimum spatial resolution reliably achieved by PHITS simulations is approximately 10-4 g/cm2, e.g. 1 micron in liquid water. For neutrons, PHITS can handle transport down to 10-5 eV using the neutron data library.

Reference
  1. Y. Nara et al. "Relativistic nuclear collisions at 10A GeV energies from p+Be to Au+Au with the hadronic cascade model", Phys. Rev. C, 61, 024901 (2000).
  2. A. Boudard et al. "New potentialities of the Liege intranuclear cascade model for reactions induced by nucleons and light charged particles", Phys. Rev. C, 87, 014606 (2013).
  3. Y. Sawada et al. "Intranuclear cascade with emission of light fragment code implemented in the transport code system PHITS", Nucl. Instrum. Meth. B, 291, 38-44 (2012).
  4. H. Takada et al. "An Upgraded Version of the Nucleon Meson Transport Code: NMTC/JAERI97", JAERI-Data/Code 98-005, Japan Atomic Energy Research Institute (1998).
  5. K. Niita et al. "Analysis of the (N,Xn') Reactions by Quantum Molecular-Dynamics Plus Statistical Decay Model", Phys. Rev. C, 52, 2620-2635 (1995).
  6. S. Furihata, "Statistical analysis of light fragment production from medium energy proton-induced reactions", Nucl. Instrum. Meth. B, 171, 251-258 (2000).
  7. T.W. Armstrong et al. "SPAR, a FORTRAN program for computing stopping powers and ranges for muons, charged pions, protons and heavy ions", ORNL-4869, Oak Ridge National Laboratory (1973).
  8. H. Geissel, C. Scheidenberger, P. Malzacher, J. Kunzendorf, H. Weick, ATIMA, GSI
  9. X-5 Monte Carlo Team, MCNP - A General N-Particle Transport Code, Verion 5, Volume I: Overview and Theory, LA-UR-03-1987, Los Alamos National Laboratory (2003).
  10. K. Shibata et al."JENDL-4.0: A New Library for Nuclear Science and Engineering", J. Nucl. Sci. Technol., 48, 1-30 (2011).
  11. R.E. MacFarlane et al., The NJOY nuclear data processing system version 91, LA-12740-M, Los Alamos National Laboratory (1994).
  12. D.E. Cullen et al., EPDL97: The Evaluated Photon Data Library, '97 Version, UCRL-50400, 6, Rev. 5, Lawrence Livermore National Laboratory (1997).
  13. S.T. Perkins et al., Tables and Graphs of Electron-Interaction Cross Sections from 10 eV to 100 GeV Derived from the LLNL Evaluated Electron Data Library (EEDL), Z = 1 - 100, UCRL-50400, 31, Lawrence Livermore National Laboratory (1991).
  14. J.A. Halbleib et al., ITS Version 3.0: The integrated TIGER series of coupled electron/photon Monte Carlo transport codes, SAND91-1634, Sandia National Laboratories (1992).