Features of PHITS3.27
PHITS development team, Mar. 2022

1

Number of PHITS users since 2019

Top 10 countries
Country#usersJapan

1709

Indonesia

179

United States

162

Argentina

89

Myanmar

87

South Korea

84

France73

Philippines

60*@03/14/2022, based on the location of users’ affiliation

China

57

3,082 new users from 59 countries in three years!

Germany/Italy

44

Number of PHITS users in each country since 2019*https://phits.jaea.go.jp/usermap/PHITS_map_userbase.html

2Map of Models Recommended to Use in PHITS
NeutronProton, Pion(other hadrons)

NucleusLow ← Energy → High

1 TeV

e- / e+Photon
1 TeV

1 TeV/uIntra-nuclear cascade (JAM)+ Evaporation (GEM)

JAMQMD+ GEM

3.0 GeV

Intra-nuclear cascade (INCL4.6) d+t
Evaporation (GEM)
200 MeV
20 MeV

Muon

3He

JENDL-4.0/He

Quantum
Molecular
Dynamics(JQMD)+ GEM

VirtualPhotoNuclear
JAM/
JQMD
+GEM
200 MeV

PhotoNuclear
EPDL97
or
EGS5
EGS5α

JAM/
JQMD
+GEM
+
JENDL
+
NRF

10 MeV/u
ATIMA
Nuclear
1 MeV+
Data Library
Original
Stopping power (ATIMA)(JENDL-4.0) 1 keV or track structure (KURBUC)
1 keV 1 keV
Muonic **Track+atom + structure *Only in water
EGM
0.01 meVCapture

1 meV

JENDL-4.0/HE is recommended to use for shielding design of accelerators**Nuclear reaction models are still recommended to use for the other cases3

Major Upgraded Features in ver. 3.27
Upgraded Points from ver. 3.24
 The generalized track structure mode has been upgraded and renamed
as ITSART Functions to read nuclear data libraires for deuteron, alpha particle, andphoton in ACE format have been implemented*1)
 A function to automatically determine the maximum library energy ofeach nucleus from address file (xsdir) has been implemented*2)
 A conversion program from the cross sections contained in the EXFOR
database to [frag data] has been developed*2)
 GUI of RT-PHITS has been developed The limitation of the number of particle types (part) specified in a tally
was removed*2)
 A event generator mode using high-energy nuclear data library has
been developed*2)
 Memory used in [t-yield] and [t-dchain] has been substantially reduced*2)*1) Under support of ASTOM R&D and NAIS
*2) Under support of RIST4

Major Upgraded Features in ver. 3.27
Upgraded Points from ver. 3.24
 The RI source function has been improved to generate all kinds ofradiation emitted from decay when proj = all is set*1)
 The anatally function has become applicable to most tallies*1)
 Treatment of thermal scattering libraries has been improved, withsupport for continuous energy thermal scattering libraries, mixed elastic
thermal scattering libraries and libraries containing SANS models*2)
 The Low-Earth Orbit (LEO) mode has been added to the cosmic-ray
source mode

*1) Under support of CCSE of JAEA
*2) Under the support of Dr. JI Marquez Damian at European Spallation Source

5Generalized track structure mode ITSART Track structure mode for arbitrary ions and materials is implemented
 Every collision of protons/Ions is calculated not only in water
Trajectories of 1 MeV protons and secondary electrons in Si

Track structure mode ITSART

Default : ATIMA+[Delta-ray]

How to use?
 In [Track structure], if mID of the cell is -1, ITSART is used.
 In [Track structure], even if mID of the cell is 1, ITSART is used outside theenergy/particle range of ETS and KURBUC
 Chemical forms can be considered by specifying, for example,chem = H2O 1.0 N2 2.0
in [Material] section.
*T. Ogawa, Y.Hirata, Y.Matsuya, T.Kai, Scientific Reports, 24401 (2021)

6

Functions to read d, α & photo-nuclear libraries
 Data should be in ACE format (the same as neutron and proton libraries) Deuteron data library JENDL/DEU-2020* is contained in XS folder α & photon-nuclear data libraries should be downloaded from elsewhereNeutron energy spectra on 9Be(left) and C(right) thick target[frag data] format is replaced by ACE format since version 3.27
*S. Nakayama et al., J. Nucl. Sci. Technol. 58(7), 805-821 (2021)

7

Automatic determination of [data max][parameters] setting Proton library: dmax(1)>0
 HE neutron library: dmax(2) > 20
 deuteron library: dmax(15) > 0
 α-particle library: dmax(18) > 0
 Photonuclear library:dpnmax > 0 & ipnint = 1

Comprehensive benchmark ofPHITS+JENDL-4.0/HE for shielding

Yes
Automatically search the library withextension specified by “lib” parameter inxsdir for each nucleus

Yes

NoUse thelibrary Conventional library (neutron)
 Nuclear reaction model (others)

Neutron fluence from p(52MeV) on C
Y. Iwamoto et al., J. Nucl. Sci. Technol.,
DOI:10.1080/00223131.2021.1993372

 High energy nuclear data library can be easily used in PHITS JENDL-4.0/HE is adopted in the recommendation setting for shielding etc.

8Exfor data converter: Exfor2fragWhat is Exfor2frag?
 Convert the experimental cross section data contained in the EXFOR
database* to [frag data] file Exactly reproduce the experimental data in the PHITS simulation* https://www-nds.iaea.org/exfor/Convert

EXFOR file (left) for 16O(p,x)1H is converted into [frag data] file (right)
 In [frag data] section, set a generated fragdata.dat file with opt=5
See phits/utility/exfor2frag in more detail

9Extension from DICOM2PHITS to RT-PHITSRadioTherapy package based on PHITSDICOM data（RT-Plan, PET-Image, CT-Image, RT-Structure, RT-Dose）
Radiotherapy

Radionuclide Therapyor
RT-PHITS（Plan2PHITS, PET2PHITS, CT2PHITS, PHITS2DICOM）

General
DICOM software

Under development

orIntermediate files (Beam geom., Source info, Patient geom., Tally setting)
Detailed analysis

InputCreater4PHITS• 3D view• DVH analysis

PHITS input

Phase Space File

PHITSphits/utility/RT-phits

Dosedistribution

10

Development of RT-PHITS GUI (β ver.）RadioTherapy package based on PHITS User support function for RT-PHITS by a Python tkinter GUI
 It allows distribution analysis by interactive interface

Display activity dist. In PET images

Display dose dist. on CT images

Interactive setting

Display beam
device geometry（under development）

DVH analysis

See phits/utility/RT-phits/RTphitsGUI.pptx for more details

11Upcoming Futures
We are planning to …

 Improve the track structure mode Extension of the precise mode applicable to other elements/compounds
 Improvement of ITSART to be capable of handling excitation
 Development of coupling modules connecting to material & life sciences Improve affinity to nuclear data library Inclusion of JENDL-5.0 and photonuclear data libraries

 Develop user support functions
 Special editor for making PHITS input file Improvement of RT-PHITS

 Improve accuracy and nuclear reaction model Improvement of JQMD ver. 2.0 to be faster and more accurate Improvement in the evaluation methods for both statistical &systematic uncertainty
 Comprehensive V&V12