5.16. [ Data Max ] section

dmax(i), which is the maximum energy of library use for i-th particle, can be defined for each nucleus in materials in this section. Maximum 6 [Data Max] sections are allowed to be used in a input file.

The format is as follows.

[Data Max]
  part = neutron proton
  mat   nucleus  dmax
  all   Fe       20
  5     all      50
  3     56Fe     150

Particle is defined in the first line as part=. From version 3.27, the upper limits of the nuclear data libraries for each nuclide are automatically determined from the address file (xsdir) coupled with dmax(i) and lib(i), where i is the particle type. Therefore, it is not necessary to specified [data max] section by users if they want to use all high-energy nuclear data libraries such as JENDL-4.0/HE defined in their xsdir file.

Three columns, mat, nucleus, and dmax, can be used. If you want to change the order of (mat) (nucleus), set as nucleus mat. You can use the skip operator non. In the mat column, material numbers can be specified and all eans all materials. In the nucleus column, you can specify nucleus as 56Fe and 26056 type. You can use Fe or 26000, which specifies all isotopes of Fe. You can also use all. By the number (MeV) in the dmax column, you define the maximum energy of library use for the nucleus.

If the same nucleus is defined in the [data max] sections, the latest definition has priority in an input file.

The values of dmax(1) and dmax(2) defined by [parameters] section should be the maximum value in [data max] section.

When dmax is set to higher than the upper limit of the library energy, PHITS simply uses the cross section for the highest energy up to dmax. For example, when dmax is set to 200 in the case of using the default library (JENDL-4.0), the cross section at 20 MeV is continuously used up to 200 MeV. When you use two (or more) libraries have different upper limits from each other, you should set the limit for each nucleus in this section.

If kmout=1 is specified in [parameters] section, the values of dmax for each nucleus in the materials are shown in the output file, file(6) (D=phits.out).