3.3. Card 3: Calculation options

Table 3.3.1 IDIVS PHITS: same name

Value	Description
(D=50)	Number of time subdivisions of each irradiation time step.

Table 3.3.2 ICHAIN PHITS: same name

Value	Description
(D=100 [1])	Maximum number of linear chains considered per nuclide.

Table 3.3.3 ILCHAIN PHITS: same name

Value	Description
(D=ICHAIN, DP=100)	Maximum length of each linear decay chain.

Table 3.3.4 ITDECS PHITS: same name

Value	Description
(D=0, DP=1; requires JMODE=0,2)	Option for considering depletion/burnup of target nuclides by the continuous primary beam irradiation [2].
= 0	Target nuclides are not depleted from beam interactions.
= 1	Beam-induced target nuclide burnup is taken into account.

Table 3.3.5 ITDECN PHITS: same name

Value	Description
(D=0, DP=1; requires JMODE=1,2)	Option for considering depletion/burnup of target nuclides by transmutations from neutrons below 20 MeV [3].
= 0	Target nuclides are not depleted from neutron reactions.
= 1	Neutron-induced target nuclide burnup is taken into account.

Table 3.3.6 ISOMTR PHITS: same name

Value	Description
(D=0, DP=2,0 [4])	Option for treatment of isomers in the nuclide yield file.
= 0	Isomers are considered distinct nuclides from their ground state.
= 1	Isomers are treated as if they were the ground state.
= 2	All isomers and the ground state are treated as being produced with equal probability.

Table 3.3.7 IXSRALL PHITS: same name

Value	Description
(D=0, DP=1)	Option for accelerating calculation time by reading and storing all neutron reaction cross section data into memory at once rather than reading and then discarding each cross section once per region.
= 0	\(\sigma\) data is read once per region (slower for multiple regions but consumes less memory; “old” methodology).
= 1	\(\sigma\) data is read once and saved for calculations in all regions (faster for multiple regions but consumes more memory).

Table 3.3.8 IFISYD PHITS: same name

Value	Description
(D=1, DP=0,2 [5])	Option for \(<\)20 MeV neutron-induced fission yields.
= 0	No effect (these fissions not considered).
= 1	Fission yields from induced fission reactions are taken into account. (requires JMODE=1,2)
= 2	Fission yields from induced and spontaneous fission reactions are taken into account. If JMODE= -1 or 0, only spontaneous fission yields are taken into account.

Note: It is strongly recommended to increase ICHAIN from its default value of 100 to 500\(\sim\)1000 instead (demonstrated in Sec. 4.7) if your problem contains fissioning nuclides as fission greatly increases the number of possible decay and reaction pathways/chains producing each nuclide.

[1]

[T-Dchain] increases default ICHAIN to 500 if IFISYD\(>\)0 and any A\(>\)220 target nuclides are found.

[2]

This option toggles burnup calculations using the negative yields present in the *.dyld file generated by PHITS to deplete target nuclides and using Pearlstein systematics [Pea87] with the ebeam and amp parameters (assuming a proton beam) to reduce the beam flux as the target is depleted.

[3]

When a neutron reaction occurs, one nuclide is transmuted into another. This option allows for the user to either ignore or consider destruction/removal of target nuclides in these reactions.

[4]

If the PHITS parameter igamma is set to 3, as recommended, the PHITS default for ISOMTR is 0. Note igamma=2 is the parameter’s default value in PHITS if not specified, setting ISOMTR= 2 by default.

[5]

[T-Dchain] sets IFISYD to 0 by default; if any target nuclides have A\(>\)220, it is set to 2 instead.