2. Version change log

Version 3.36 (2026-03-31)

  • Parallel computation capability has been implemented in DCHAIN, enabling independent parallel execution for each region or mesh. Both OpenMP and MPI parallelization are supported, while hybrid (OpenMP+MPI) parallelization is not supported. Parallel execution can be controlled by specifying $OMP or $MPI at the beginning of the input file.

Version 3.35 (2025-03-28)

  • Several bugs in the PHITS RI-source [source] output mode have been fixed. In addition, a new function to output [source] with proj = photon as well as all, electron, positron, and alpha was implemented. Note that the meaning of PHITSOUT has been changed in accordance with this improvement.

  • The default value of IGSDEF and IGSORG is now set to 0.

  • Some activation cross section libraries have been updated to include the data for \((n,n')\) reactions.

  • The activation cross-section of JENDL-5 has been implemented. INXSLIB=100,101,102,120 changes.

Version 3.25 (2021-07-12)

  • New input parameters: IFISYD, INFYLIB, INFYERG, ISFYLIB, HNFYLIB, HSFYLIB, ERGTHFA, ERGFAFU, RADLEV, ILCHAIN

  • DCHAIN can now handle neutron-induced and spontaneous fission yields (previously disabled since 2015) with GEF model-calculated yields and (optionally) supplemented with evaluated data yields where available (with newly added libraries ENDF/B-VIII.0, JENDL/FPY-2011, JEFF-3.3, and hybrid combinations of these).

  • New option to output PHITS RI-source [Source] sections, allowing for \(\gamma\), e\(^-\), e\(^+\), and \(\alpha\) decay emissions, at exact energies, to be simulated (see new PHITSOUT options).

  • The maximum number of decay chains per nuclide (ICHAIN) and maximum length of chains (ILCHAIN) can now be controlled independently.

  • Manual Sections 4.6.3 and 4.7 relating to fission added.

Version 3.24 (2021-03-22)

  • Missing files for quantities (neutron flux/nuclide yields) set to be ignored by JMODE no longer result in the calculation aborting.

  • Output \(^*\).pht files with photon [Source] sections are now formatted as multi-sources rather than individual sources when using xyz mesh geometries; see Section 4.13.

  • The ordering of hybrid activation cross section libraries was changed to place higher priority on the dedicated activation libraries; see the footnote in Section 4.6.2.

  • Manual Sections 4.1, 4.2, 4.9, and 4.13 added.

Version 3.23 (2021-01-18)

  • The new main DCHAIN reference paper [RMA+20] has been published.

  • New input parameters: IMTCARD, IMTCNUM, IMTCMETA, THMATND, IDOSUNIT, IANGBPWR, ANGELOUT_NUCLIDES

  • New \(^*\).mat output file added containing PHITS-formatted materials cards.

  • Photon effective dose rate units (dose, area, and time parts) can now be changed.

  • Data for individual nuclides and the beam schedule can be added to ANGEL plots.

  • Regions in ANGELOUT_REGION can now be listed by number IREGS or name HRGCMM.

  • Example calculation added in Section 4.12 of this manual.

  • Maximum number of materials in an xyz mesh problem increased from 10 to 100.

  • Fixed bug in \(^*\).yld output where target nuclide inventories were double-counted.

Version 3.21 (2020-09-14)

  • New input parameters: IXSRALL, IWRCHNUC, IDOSECF, IREDUFMT, IRDONCE, IPLOTMODE, IPLOTAXIS, FOAMOUT, FOAMVALS

  • Added ENDF/B-VII.1, ENDF/B-VIII.0, JENDL-4.0, JEFF-3.3, FENDL/A-3.0, EAF-2010, BROND-3.1, CENDL-3.1, TENDL-2017, and hybrid activation cross section libraries. (files added in version 3.18, INXSLIB values updated in 3.21)

  • Added ICRP 116 effective dose and ICRP 74 ambient dose equivalent coefficients.

  • Added support for tetrahedral mesh geometries (mesh = tet).

  • Added 2D visualization options for xyz mesh geometries.

  • Added propagation of neutron flux statistical uncertainties from \(^*\).dtrk file, completing implementation of statistical uncertainty propagation from PHITS to DCHAIN.

  • Introduced reduced format for [T-Track] and [T-Yield] output from [T-Dchain].

  • Greatly improved efficiency, especially for large meshes (old behavior still available)

  • Added ability to track only specific nuclides in \(^*\).dcs file.

  • Source bounding box in \(^*\).pht is now automatically defined for xyz meshes.

  • New \(^*\).csv and \(^*\).foam output files added for tetrahedral meshes.

  • New \(^*\).pxy.ang, \(^*\).pxz.ang, and \(^*\).pyz.ang output files added for xyz grid meshes.

  • Maximum allowable number of time steps increased from 100 to 1000.

  • The number of allowable regions is now unlimited (limited by system memory).

Version 3.14 (2019-08-18)

  • New input parameters: INXSLIB, IDCYLIB, HNXSLIB, HDCYLIB, IWRTCHN, CHRLVTH, IWRCHDT, IWRCHSS, IERTDCHO, HENMTCF

  • Added new ENDF/B-VIII.0, JENDL/AD-2017, and hybrid decay library options

  • Added JENDL/AD-2017 activation cross section library

  • New \(^*\).dcs output file added containing decay chain schemes and reaction rates.

  • Added ability to propagate nuclide yield statistical uncertainties from \(^*\).dyld file