3.5. カード4: 出力条件設定

このカードの従来の入力パラメータ IPRTB1、IPRTB2、RPRTB2、IPRTB3 は廃止され、現在は機能しない（旧カード8と、それに依存する ANCZA および HANCZA も同様）。

表 3.5.1 IYILD PHITS: same name

値	説明
(D=0, DP=2)	最終出力時刻における領域別の核種生成量ファイル（*.yld）の出力を指定する [1]。
= 0	出力しない。
= 1	生成核種の最終インベントリを出力する。
= 2	カード9a/9bに記載した標的核種の初期インベントリを含め、全核種の最終インベントリを出力する。

表 3.5.2 IGGRP PHITS: same name

値	説明
(D=3, IMODE\(\geq\)2が必要)	出力する\(\gamma\)線スペクトルのエネルギー群構造を選択する。
= 0	無効（IMODE\(<\)2の場合の既定値）。
= 1	HILO86の22群構造 (20.0\(\sim\)0.01 [MeV])
= 2	BERMUDAの41+3群構造 (20.0\(\sim\)0.0 [MeV])
= 3	VITAMIN-Jの42群構造 (50.0\(\sim\)0.001 [MeV])
= 4	THIDA-2の54群構造 (3.0\(\sim\)0.03 [MeV])
= 5	THIDA-2の78群構造 (7.6\(\sim\)0.03 [MeV])
= 6	FNSの40群構造 (14.0\(\sim\)0.01 [MeV])
= 99	FISPACT3の23群構造 (20.0\(\sim\)0.01 [MeV])

表 3.5.3 IBETAP PHITS: same name

値	説明
(D=1)	陽電子（\(\beta^+\)）の放出に伴う消滅光子を\(\gamma\)線スペクトルへ含めるかどうかを指定する。
= 0	消滅ガンマ線を考慮しない。
= 1	放出陽電子1個当たり2個の消滅光子を考慮する。

表 3.5.4 ACMIN PHITS: same name

値	説明
(D=0.0)	各出力ファイルへ記載する核種の放射能またはインベントリの下限値を指定する（DP=10\(^{-20}\)）。
= 0.0	下限値を全放射能の10\(^{-10}\)倍へ自動設定する。
\(>\) 0.0	放射能の下限値 [Bq]を指定する。
\(<\) 0.0	値の絶対値を核種数密度の下限値 [atoms/cm\(^3\)]として使用する。

表 3.5.5 ISTABL PHITS: same name

値	説明
(D=0)	出力ファイルへ安定核種を記載するかどうかを指定する。
= 0	出力しない。
= 1	出力する。

表 3.5.6 IGSDEF PHITS: same name

値	説明
(D=0)	\(\gamma\)線スペクトルをMCNPの線源定義（SDEF）カードとして*.gsdへ出力するかどうかを指定する。IMODE\(\geq\)2かつIGGRP\(>\)0の場合に有効である。
= 0	出力しない。
= 1	出力する。

表 3.5.7 IGSORG PHITS: same name

値	説明
(D=0)	\(\gamma\)線スペクトルの起源核種を*.gsoへ出力するかどうかを指定する。IMODE\(\geq\)2かつIGGRP\(>\)0の場合に有効である。
= 0	出力しない。
= 1	出力する。

表 3.5.8 IWRTCHN PHITS: same name

値	説明
(D=0, DP=1)	崩壊系列情報ファイル（*.dcs）を出力するかどうかを指定する。関連する各崩壊系列と、それが照射・冷却・出力時刻から決まる各計算時間ステップの核種インベントリ変化へ与える寄与を記載する。詳細は 4.10節 を参照すること。
= 0	出力しない。
= 1	出力する。

表 3.5.9 CHRLVTH PHITS: same name

値	説明
(D=-1.0, IWRTCHN=1が必要)	各計算時間ステップにおいて、崩壊系列を*.dcsへ出力するためのインベントリまたは放射能変化量の下限条件を指定する [2]。
\(<\) 0.0	値の絶対値を核種インベントリ変化量の下限値 [atoms/cm\(^3\)]として使用する。
= 0.0	下限値を設定せず、DCHAINが構築した全系列を出力する。
\(>\) 0 & \(\leq\) 10	核種インベントリの相対変化量を下限値として使用する。
\(>\) 10.0	核種の絶対放射能変化量の下限値 [Bq]として使用する。

表 3.5.10 IWRCHDT PHITS: same name

値	説明
(D=0, IWRTCHN=1が必要)	各崩壊系列を構成する過程が終端核種のインベントリ変化へ与える寄与を*.dcsへ出力するかどうかを指定する。
= 0	出力しない。
= 1	出力する。

表 3.5.11 IWRCHSS PHITS: same name

値	説明
(D=0, IWRTCHN=1が必要)	各照射期間における最初のIDIVS-1サブステップについても、崩壊系列データを*.dcsへ出力するかどうかを指定する [3]。
= 0	出力しない。
= 1	出力する。

表 3.5.12 IWRCHNUC PHITS: not settable

値	説明
(D=0, IWRTCHN=1が必要)	*.dcsへ出力する崩壊系列の終端核種を限定する。
= 0	すべての終端核種を通常どおり出力対象とする。
= \(n\)	*.dcsへ出力する終端核種を\(n\)個に限定する。この行の直後の\(n\)行に、カード9b に示すTGNZASまたはHTGNZAS形式で核種名を1行ずつ記載する [4]。

表 3.5.13 IDOSECF PHITS: same name

値	説明
(D=1)	各出力ファイルに記載する光子線量を計算する際に使用する、光子フルエンスから実効線量\(E\)または周辺線量当量\(H^*(10)\)への換算係数を選択する [5]。
= 0	旧既定値。出典が不明で、値に問題があると考えられる。
= 1	AP（前方－後方照射）のフルエンス－実効線量換算係数。ICRP 116 [ICR10]。
= 2	PA（後方－前方照射）のフルエンス－実効線量換算係数。ICRP 116 [ICR10]。
= 3	LLAT（左側方照射）のフルエンス－実効線量換算係数。ICRP 116 [ICR10]。
= 4	RLAT（右側方照射）のフルエンス－実効線量換算係数。ICRP 116 [ICR10]。
= 5	ROT（回転照射）のフルエンス－実効線量換算係数。ICRP 116 [ICR10]。
= 6	ISO（等方照射）のフルエンス－実効線量換算係数。ICRP 116 [ICR10]。
= 7	フルエンス－周辺線量当量\(H^*(10)\)換算係数。ICRP 74 [IInternationalCoRProtectionInternationalCoRUnits96] [6]。

表 3.5.14 IDOSUNIT PHITS: same name

値	桁	説明
(D=333, \(\mu\)Sv\(\cdot\)m\(^2\)/h)		各出力ファイルに記載する光子実効線量\(E\)または周辺線量当量\(H^*(10)\)の率 [7] の単位を選択する。3桁の数値XYZで指定し、各桁はそれぞれ線量、面積、時間の単位を表す。たとえば121はpSv\(\cdot\)cm\(^2\)/sを表す。
= 1	X［線量単位］	pSv
= 2	X［線量単位］	nSv
= 3	X［線量単位］	\(\mu\)Sv
= 4	X［線量単位］	mSv
= 5	X［線量単位］	Sv
= 1	Y［面積単位］	mm\(^2\)
= 2	Y［面積単位］	cm\(^2\)
= 3	Y［面積単位］	m\(^2\)
= 1	Z［時間単位］	1/s
= 2	Z［時間単位］	1/min
= 3	Z［時間単位］	1/h
= 4	Z［時間単位］	1/day

表 3.5.15 IMTCARD PHITS: same name

値	説明
(D=0, DP=1)	核種インベントリをPHITS/MCNP形式の[Material]セクションとして*.matへ出力するかどうかを指定する。
= 0	出力しない。
= 1	全出力時刻を1ファイルへ出力する。
= 2	ITOUTの各出力時刻を別ファイルへ出力する。PHITSのinfl:機能で挿入ファイルとして使用する場合に便利である。

表 3.5.16 IMTCNUM PHITS: same name

値	説明
(D=0, IMTCARD=1が必要)	各領域の物質番号 [8] を全出力時刻で一定にするか、出力時刻 ITOUT ごとに増加させるかを選択する。
= 0	全時刻で物質番号を一定にする。
= 1	出力時刻 ITOUT ごとに物質番号を調整し、増加させる。

表 3.5.17 IMTCMETA PHITS: same name

値	説明
(D=0, DP=2, IMTCARD=1が必要)	*.matの[Material]セクションへ出力する準安定核種の扱いを指定する。
= 0	準安定核種を出力しない。
= 1	MCNP形式の\(ZZZAAA\)値 [9] を使用する。\(AAA^\prime = AAA + 300 + 100\cdot M\)とし、\(M\)は準安定・核異性状態を表す1または2である [10]。
= 2	すべての準安定核種のインベントリを基底状態のインベントリへ加算する。

表 3.5.18 THMATND PHITS: same name

値	説明
(D=-1.0, DP=1.0E-6, IMTCARD=1が必要)	核種を*.matの[Material]セクションへ含めるために必要な数密度の下限値を指定する。
\(\leq\) 0	値の絶対値を数密度の下限値 [atoms/cm\(^3\)]として使用する。
\(>\) 0 & \(\leq\) 1	全原子数密度に対する割合を下限値として使用する。

[1]

単位体積当たりの全原子数に対する割合が\(10^{-12}\)以上の核種だけを出力する。IYILD=1では初期標的核種を全原子数に含めず、IYILD=2では含める。

[2]

各核種のインベントリは、その核種へ至る生成、反応、崩壊経路の寄与を合計して求める。これらは線形の「系列」を構成する。CHRLVTHは、各系列が核種インベントリの変化へどの程度寄与すれば出力対象となるかを制御する。

[3]

各照射期間はIDIVS（既定値50）個の計算サブステップへ分割される。IWRCHSSは、最後のサブステップだけを表示するか、すべてのサブステップを表示するかを指定する。詳細は 4.8節 を参照すること。

[4]

核種名には質量数が必要である。たとえばFe-55と26055.0は使用できるが、Feと26000.0は使用できない。

[5]

DCHAINの線量率出力の意味と導出は、4.9節 で説明する。

[6]

10 MeV以下の光子にはICRP 74 [IInternationalCoRProtectionInternationalCoRUnits96]、10 MeVを超える光子にはPelliccioni [Pel00]の係数を使用する。切替点には不連続がある。

[7]

DCHAINの線量率出力の意味と導出は、4.9節 で説明する。

[8]

PHITSの物質番号は5桁に制限される。IMTCNUM=0では、通常、regメッシュの物質番号をセル番号に、その他のメッシュでは領域インデックスに合わせる。IMTCNUM=1では、出力時刻のインデックスを物質番号の末尾へ追加する。領域数、出力時刻数、セル番号の大きさによって5桁に収まらない場合は、可能な範囲で割当方法を調整する。

[9]

この形式を正しく扱えるのはMCNPだけである。PHITSでは非常に重い核種として解釈される。現在、PHITSの[Material]セクションは準安定核種を直接指定できない。

[10]

MCNPでは1、2、3、4を使用できるが、DCHAINが現在扱うのは第1・第2核異性状態だけである。