Features of PHITS3.10
PHITS development team, Mar. 20191

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

NucleusLow ← Energy → High

1 TeV

JAMQMD+ GEM

3.0 GeV

Intra-nuclear cascade (INCL4.6) d+t
Evaporation (GEM)
3He

20 MeV

0.01 meVe- / e+Photon
1 TeV

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

NuclearData Library(JENDL-4.0)
+EGM

Muon

α

1 MeV
1 keV

Quantum
MolecularDynamics(JQMD)+ GEM
10 MeV/uIonization
ATIMA

VirtualPhotoNuclear
JAM/
JQMD
+
GEM
200 MeV

PhotoNuclear
EPDL97
orEGS5EGS5ATIMA+Original

1 keV

JAM/
JQMD
+GEM+JENDL+NRF

1 keV

Muonic **Trackatom + structure *Only in water
Capture 1 meV

Neutrino with energy below 20 MeV can cause interaction with 1H, 2H,and electron, though it is not written in this map

2

Major Upgraded Features in v3.08
PHITS 3.08 was distributed only through PHITS tutorial

Upgraded Points from v3.02


JAMQMD2 is developed for high-energy nucleus-nucleus interaction
Neutrino interaction model is developed & implemented[weight window] & [t-wwg] become applicable to xyz-mesh
Emission of Auger electrons, internal conversion electrons, & characteristic
X-rays can be considered in the RI source function
Several types of dose conversion coefficients are implemented
Source generation function in tetrahedral mesh is developed
Tetrahedral geometry created by CAD become readable in PHITS
Neutron and γ-ray emissions are considered at once in GEM
Detector response can be calculated even when [forced collisions] is used
Probability density of the number of interactions per history can becalculated using [t-interact] (former [t-star])
Several activation cross section data libraries used for [t-yield] and [tdchain] are implemented[t-heat] is merged into [t-deposit]
PHITS executable file in MPI version for Windows is implemented3

Development of JAMQMD version 2
Nucleus-nucleus interaction model for 3 GeV/n ~ 1 TeV/n
Useful for cosmic-ray and high-energy physics simulations

100GeV/n 12C

200200断面積 (mb)
Cross section

160160Exp (実測)
実験値(Cecchini et al. 1993)Exp.
jamqmd Ver1
JAMQMD
Ver.1 ＋GEM
ktuy
JAMQMD
Ver. 2 ＋GEM

8080μ, γ10km

120120Improved
精度改善

40400088ν91011129101112Atomic numberof residual nuclei
生成核原子番号Cross section of 14.6AGeV 28Si(27Al,x)1313nAirshower simulation

 Improve the accuracy of fragment yields
 Automatically used for nucleus-nucleus interaction above 3GeV/n*T.Ogawa et al., Phys. Rev. C., 98, 024611 (2018)4

Development of Neutrino Interaction Model*Applicable to neutrino interaction with 1H, 2H or e2H disintegration

νe →e+conversione (ν, ν’) e’
scattering
Neutrino can penetrate
Energy distribution of secondary
nuclear reactor
particles produced by 2H(νe,X)
 Consider only when ntrnore is set to 1 (D=0) in [parameter] Recommend to use this option in combination with [Forced Collisions] Interaction with heavier target & neutrino oscillation are under consideration*Some functions are available only after version 3.10

5xyz-mesh for [weight window] & [t-wwg]

Few particles canreach lateral edges
because cells aredivided only to thevertical direction

mesh = regmesh = xyz

Statistical uncertainties of [t-track] obtained using [t-wwg] with reg (left)
and xyz (right) meshes. Cells are divided only to the vertical direction.
(see lecture/advanced/weight in more detail)You do not have to divide a cell into many pieces
for defining different weight windows6

Major Upgraded Features in v3.10
PHITS 3.10 is distributed since April 2019

Upgraded Points from v3.08
 Procedure for coupling PHITS with thermal analysis software such as ANSYS
Fluent is established Sample files for using JENDL-4.0/HE are provided
 Function to read magnetic field maps written in xyz or r-z grid is developed [t-dchain] has become applicable to xyz-mesh Error and warning ID numbers are introduced for some messages
 Accuracy of electron track-structure mode is improved Detector resolution in any form given by user-defined function can be considered INC-ELF, an intra-nuclear cascade model developed in Kyushu Univ., is improved A new gshow option is introduced to visualize tetrahedral and voxel geometries Time for reading lattice structure is reduced by introducing compressed format
 Function to calculate water-equivalent dose is implemented in [t-deposit]
 Remaining batch number has become adjustable using batch.out
7Coupling with thermal analysis
 Create geometry in thermal analysis code such as ANSYS Fluent
 Meshing to tetrahedral mesh => Output in NASTRAN bulk data format
 Import => PHITS radiation transport => Output in OpenFOAM field data Import => Thermal analysis by ANSYS Fluent

PHITSFluent
NASTRANbulk dataformat

ProtonbeamFluent
OpenFOAM
Field dataformat

Temperature
Tetra-mesh geom.
Deposition energy
Read & write in general tetrahedral mesh geometry are possible
and thus seamless coupling to other physics becomes possible
8Sample files for using JENDL-4.0/HE
What is JENDL-4.0/HE?
 High-energy nuclear data library for neutrons and protons up to 200 MeV
 Precisely reproduce the DDX of nuclear reactions particularly for lighter targets
 Data for neutrons below 20 MeV are the same as JENDL-4.0
10Important notice!

29DDX (mb/sr/MeV)

Be(p,n) for Ep = 9.7 MeV at 15 deg

10100−20

 Data for only several elements areincluded in the PHITS package. Datafor other elements will be released via
JAEA website in future Event generator mode cannot be usedin combination with JENDL-4.0/HEExp. (Iwamoto et al.)JENDL−4.0/HE
INCL
48Emitted Neutron Energy (MeV)DDX of 9Be(p,n) reaction calculated by
PHITS using JENDL-4.0/HE or INCL Calculation of deposition energies isnot feasible using either [t-deposit] nor[t-heat]
Applications are limited to e.g. shielding
calculation and neutron source design
Please see phits/recommendation/jendlHES. Kunieda et al., JAEA-Conf 2016-004 (2016), https://wwwndc.jaea.go.jp/ftpnd/jendl/jendl40he.html

Function to Read Magnetic Field Map
 Both xyz and r-z grid fields can be read
 Both map- and list-type data can be read
 Applicable to symmetry and asymmetry fields to certain axes Particle motions are analyzed by Runge-Kutta methodwith magnetic field

Sourcewithout magnetic field

SourceElectron trajectories in metal target with (left) and without (right) magnetic field
Please see phits/utility/magmap in more detail
Magnetic field data are provided by Dr. Sakata @ Osaka Univ.

10