Features of PHITS3.00
PHITS development team, Aug. 20171

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

NucleusLow ← Energy → High

1 TeV

Quantum
Molecular
3.0 GeV
Dynamics
Intra-nuclear cascade (INCL4.6) d (JQMD)++tEvaporation (GEM)
Evaporation
3He(GEM)
20 MeVα0.1 meV

e- / e+Photon
1 TeV

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

NuclearData Library(JENDL-4.0)

Muon

1 MeV
1 keV

VirtualPhotoNuclear
JAM/JQMD+GEMPhotoNuclear
EPDL97
orEGS5
EGS5200 MeV

10 MeV/uIonization
ATIMA

Event generator mode
all secondary particles are specified

1 keV*Trackstructure

1 meV

JAM/QMD+GEM+JENDL+NRF

1 keV*Only in water Transport photons and low-energy neutrons in default setting
 Track-structure mode was implemented (applicable only to e-/e+ in water)2

Major Upgraded Features in v2.88
Upgrade Points from v2.82
 Implement Weight Window Generator function
 Develop output option for a 3D-viewer ParaView Set ATIMA for default stopping power calculation model Improve “sum tally” function
 Improve muon and pion nuclear interaction models
 Implement radioisotope source function
 Develop JAMQMD
 Implement neutron decay process
 Revise some bugs in the EGS5 algorithm

3Weight Window Generator
What can do with Weight Window Generator [t-wwg]?
 Automatically determine the effective settings of [weight window] Help users to easily perform deep penetration calculation
1st trial
No [Weight Window]

2nd trial
Using [Weight Window]

3rd trial
Revised [Weight Window][t-track]
Relative Error
Deep penetration calculation using [t-wwg] & [Weight Window](history numbers are the same for all trial)
4Output Function for ParaView
What’s ParaView (http://www.paraview.org)?
 Open-source, multiplatform data analysis and visualization application
 Capable of drawing tally outputs in 3D picture & animation
How to activate the function?
Simply add “vtkout = 1” in a tally with “mesh = xyz”

ICRP voxel phantom
visualized by ParaView

Sample picture of ParaView

See lecture¥advance¥paraview in more detail5

RI Source Function
What can do with RI source function?
 generate photon sources with energy spectra of radioisotope (RI)
decay by simply specifying the activity and name of the RIs
 consider activities from daughter nuclides by setting a decay time use nuclear decay database DECDC*
Mono & continuous energy source
RI source(e-type = 22)(e-type = 28)[ Source ]totfact = 1.0…
e-type = 28
ni = 1
152Eu 100.
norm = 0
dtime = -10.0

dtime:
decay time

[ Source ]$ totfact =160.51…
e-type = 28
$ 1.218E-01 1.218E-01 100*2.867E-01$ 1.257E-01 1.257E-01 100*1.611E-04… repeated$ 1.698E+00 1.698E+00 100*5.907E-05$ 1.769E+00 1.769E+00 100*9.585E-05$ne = 132

Input formats for representing 152Eu
with activity of 100 Bq

Photon flux from 152Eu*A. Endo et al., JAERI 1347 (2005); equivalent to ICRP Publication 107

6Recent Upgrades: PHITS3.00
Not Registered in OECD/NEA yet

Upgrade Points from v2.88
 Development of JAMQMD2.0
 Development of track-structure mode Consideration of electro-magnetic field in EGS5 mode Extension of RI source function applicable to α and β decays Development of xyz-mesh source generation function
 Implementation of automatic termination function by CPU timeand standard deviation
 Development of volume calculation tally [t-volume] Implementation of sum-up function in [t-deposit]
 Extension of reaction channel in [counter], [t-star], & [t-product] Implementation of epsout = 2 option to draw error bar in ANGEL Introduction of weight window bias function [t-wwbg] & [ww bias]
 Revision of input format
7

Track-structure Mode
What is track-structure mode All ionization and excitation events are explicitly considered→ In contrast to continuous slowing down approximation (CSDA) Particle trajectory in microscopic scale (~ nm order) can be calculated→ 1μm resolution is minimum when CSDA is used
 Computational time becomes extremely long

Trajectories of 4 keV electrons
simulated by EGS5 mode (upper)
and track-structure mode (lower)(see recommendation/trackstructure)

Length unit can be changed by ANGEL parameter “cmum”

Currently applicable only to e- & e+ in water, but plan to extend to other cases

8Without electromagnetic field

10MV/cm

Positrons

10MV/cm

Electro-magnetic field in EGS mode

Electrons

水中に100MeV電子を入射した際の電子・陽電子フラックス
 Definition of the electro-magnetic fields is the same as previously→ Define in [magnetic field] or [electromagnetic field],
Set imagnf / ielctf = 1 in [parameters] (see lecture/advanced/option)

 Effective only when EGS mode is activated Be careful for generating too many secondary particles when highelectric field is applied

9Revision of Input Format
Major revision points
 “c” becomes unusable as comment marks in [material]→ If you set “icommat=1” in [parameters], “c” is regarded as comment
 “file(1)” is introduced to specify the PHITS installation folder name. If youset this parameter, you do not have to specify other input file names→ Each “file” parameter overwrite the setting when they are specified A new parameter “nucdata” is introduced to automatically adjust emin(2)and dmax(2) parameters suitable for JENDL-4.0→ Basically you do not have to specify emin(2) & dmax(2) “negs” becomes the control parameter for electromagnetic transport→ negs = -1: Transport only photon based on EPDL97 (default)
negs = 0: Ignore electron, positron, and photon transportnegs = 1: Transport electron, positron, and photon based on EGS5
 “s-type” in [source] needs not to be changed by energy type→ automatically judge mono-energetic or energy-distributed sources*Basically, you do not have to modify your input file suitable for new version*Exceptions: Define elements in [material] whose nuclear data are not evaluated → nucdata = 0
Use “c” as a comment remark in [material] → icommat = 1
10