[ Material ] セクション ================================================== このセクションでは体系を構成する物質を定義します。 ここで定義した物質は、[cell]セクションをはじめとする幾つかのセクションで参照されます。 書式は以下の通りです。物質番号の後に、物質を構成する元素(核種)とその組成比を交互に書きます。 .. code-block:: text [ Material ] MAT[物質番号] 元素 組成比 元素 組成比 ... ... 物質番号nは1から99999まで選択でき、MAT[n]あるいはMnの2種類の表現で指定することができます。 ただし、MAT[n]において、MATと[の間には空白文字を入れることができないのでご注意ください。 コメント文字として\#, \%, !, \$の4文字が使用できます。 バージョン2.88までは行頭の c （直後に半角スペースがある場合）もコメント文字として使用することができていましたが、それより後のバージョンではデフォルトでは使用できなくなりました。 もし、コメント文字として使用したい場合は、[parameters]セクションにおいてicommat=1としてください。 また、ファイルのインクルード文、定数定義などは、セクションの途中でも用いることが出来ます。 元素(核種)の定義 -------------------------------------------------- 上記書式の元素の箇所には、定義したい物質の元素記号を書いてください。 核種を指定する場合は、質量数を加えた表式を書いてください。 208Pb, Pb-208, 82208のどの形式でも可能です。 例えば水素は、1H, H-1, 1001と指定します。 元素記号のみ、あるいは質量数をゼロにした場合(例えばPb, 82000)は、その元素の同位体を天然存在比で展開します。 ただし、JENDL-4.0に含まれない核種に関しては展開されませんのでご注意ください。 PHITSパッケージに含まれているJENDL-4.0以外の核データライブラリーを使用する場合は、断面積ディレクトリファイルfile(7) (D=xsdir)にそのライブラリ付属のxsdirの情報を追加してください。 もし、20MeV以下の中性子データライブラリーについて、JENDL-4.0以外のものを使用する場合は、物質番号毎にそのライブラリーのidを :numref:`tbl-mat-para` のNLIBにより指定してください。 ライブラリーのidは、ライブラリー番号(2桁の数字)とデータのクラス(1文字のアルファベット)を組み合わせて定義されており、例えばJENDL-4.0の中性子ライブラリーであれば50cとなります。 中性子以外にも、PLIB, ELIB, HLIBを指定することで、それぞれ光子、電子、陽子のライブラリーを指定して利用することができます。 また、物質番号毎ではなく、核種毎に使用するデータライブラリーを指定する場合は、その核種の指定を208Pb.50c, Pb-208.50c, 82208.50cのような形式で行ってください。 ただし、基本的にPHITSでは、中性子などの入射粒子に対して各々1種類のライブラリーを指定するようにし、その指定はNLIB, PLIB, ELIB, HLIBにより行ってください。 もし、2種類目のライブラリーを指定したい場合は、中性子に関してのみ可能ですので、NLIBによる指定を行った上で、ライブラリーを変更したい核種を208Pb.50cのような形式で定義してください。 また、各物質の後に拡張子でライブラリを指定できるのは中性子のみです。 たとえば、:math:`\mathrm{^6Li}` と :math:`\mathrm{^7Li}` について、陽子のJENDL/4.0-HEライブラリを利用する場合、HLIBを用いて、 .. code-block:: text 6Li 0.001 7Li 0.999 hlib=51h としてください。 なお、idによるライブラリーの指定がない場合は、file(7)中のリストを上から検索し、その核種で最初に該当したライブラリーを使用します。 どの核データが使用されたのかを知りたい場合は、[parameters]においてkmout=1とすることで、サマリーの出力ファイルfile(6) (D=phits.out)にその情報が書き出されます。 準安定核種を指定する場合は、原子番号と質量数を並べた5桁の整数形式を使用してください。 ただし、質量数には50を加えます。 例えば、Am-244m を指定する場合は 95294 としてください。 対応する核種の核データライブラリが存在する場合は、それを参照します。 物質組成比の定義 -------------------------------------------------- 上記書式の組成比の箇所には、その物質番号で定義した物質に含まれる元素(核種)の組成比を書いてください。 正の数で書いた場合は原子数比、負の数で書いた場合は質量比の意味となります。 例えば、水( :math:`\mathrm{H_2O}` )を定義したい場合は、 .. code-block:: text MAT[1] H 2 O 1 あるいは、 .. code-block:: text MAT[1] H -2/18 O -16/18 と指定します。 水の分子量がおよそ18であることから、水に含まれる水素と酸素の質量比がそれぞれ約2/18と16/18となるためです。 粒子輸送計算時に使用される物質の密度は、[cell]セクションで定義します。 ただし、[t-deposit], [t-deposit2], [t-let], [t-sed]セクションにおいてletmatを指定する場合は、本セクションの組成比の箇所で物質密度を定義する必要があります。 その際、正の数の場合は粒子密度(単位は[ :math:`10^{24}\,\mathrm{atoms/cm^3}` ])、負の数の場合は質量密度(単位は[ :math:`\mathrm{g/cm^3}` ])として定義します。 物質パラメータ -------------------------------------------------- 核データライブラリを使う領域に関して、物質毎のパラメータをkeyword=valueの形で指定します。 挿入場所は、その物質番号のサブセクション中なら自由です。 指定できるパラメータとその意味は以下の通りです。 .. rst-class:: no-caption-number .. list-table:: **gas** :header-rows: 0 * - 値 - 説明 * - | (D=0) - | 電子の阻止能の密度効果。 * - | =0 - | 液体もしくは固体の計算に適する。 * - | =1 - | 気体の計算に適する。 .. rst-class:: no-caption-number .. list-table:: **estep** :header-rows: 0 * - 値 - 説明 * - | - | 電子輸送のサブステップの数。 * - | =n - | サブステップの数をnにする。nがbuild-inのデフォルト値より小さい時は無視される。 .. _tbl-mat-para: .. rst-class:: no-caption-number .. list-table:: **nlib** :header-rows: 0 * - 値 - 説明 * - | - | デフォルトの中性子ライブラリー番号。 * - | =id - | idをデフォルトの中性子ライブラリーとする。 .. rst-class:: no-caption-number .. list-table:: **plib** :header-rows: 0 * - 値 - 説明 * - | - | デフォルトの光子ライブラリー番号。 * - | =id - | idをデフォルトの光子ライブラリーとする。 .. rst-class:: no-caption-number .. list-table:: **elib** :header-rows: 0 * - 値 - 説明 * - | - | デフォルトの電子ライブラリー番号。 * - | =id - | idをデフォルトの電子ライブラリーとする。 .. rst-class:: no-caption-number .. list-table:: **hlib** :header-rows: 0 * - 値 - 説明 * - | - | デフォルトの陽子ライブラリー番号。 * - | =id - | idをデフォルトの陽子ライブラリーとする。 .. rst-class:: no-caption-number .. list-table:: **cond** :header-rows: 0 * - 値 - 説明 * - | (D=0) - | 伝導体設定。 * - | <0 - | 絶縁体。 * - | =0 - | 1個でも絶縁体があれば絶縁体、それ以外は伝導体。 * - | >0 - | 1個でも伝導体があれば伝導体。 阻止能の指定 -------------------------------------------------- PHITSでは、通常、物質の阻止能はndedxパラメータで指定するモデル(ATIMA、SPARなど)で計算します。 しかし、それらのモデルでは精度が十分ではない場合や、ICRU90で定義した阻止能を使いたい場合は、阻止能のデータベースを直接読み込ませることができます。 その指定方法は .. code-block:: text dedxfile = ファイル名 で指定できます。 この阻止能データベースファイルは、file(29)で指定したフォルダ(初期設定はfile(1)/data/dedx)に格納しておく必要があります。 また、様々な物質に対する陽子及び :math:`\alpha` 粒子に対してPSTAR及びASTAR [#]_ で計算した阻止能データベースを上記フォルダに格納していますので、必要に応じて指定してご使用ください。 なお、本機能はATIMAで阻止能を計算する放射線に対してのみ有効となります。 また、エネルギー分散(nedisp)や角度分散(nspred)は、従来と同様に機能します。 .. [#] https://www.nist.gov/pml/stopping-power-range-tables-electrons-protons-and-helium-ions 阻止能データベースファイルのフォーマット ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ dedxfileで指定する阻止能データベースファイルでは、\$と\#がコメントとして利用できます。 また、PHITS入力ファイルとほぼ同じく大文字小文字の区別はありません。 指定できるパラメータはunitとkfのみで、それぞれ阻止能を与えるエネルギーの単位と放射線の種類を表します。 .. code-block:: text unit = 1: MeV (全運動エネルギー、デフォルト) = 2: MeV/u(atomic mass unitあたりの運動エネルギー) = 3: MeV/n(核子あたりの運動エネルギー) unitは、最初のkfコードを指定する前に指定してください。 kfは対象となる粒子のkfコードを指定し、その後に阻止能データベースを .. code-block:: text エネルギー 阻止能 順にエネルギーの昇順で定義してください。 ここで、エネルギーはunitで指定した単位、阻止能は :math:`\mathrm{MeV\,g^{-1}\,cm^{2}}` の単位で与えてください。 阻止能の単位は変更できませんのでご注意ください。 指定したエネルギー範囲外の放射線やkfコードを指定しなかった放射線に対しては、ATIMAで計算した値がそのまま利用されます。 複数の放射線に対する阻止能を定義する場合は、1つのデータベースファイルにまとめてください。 複数用意したデータベースファイルを1つの物質に対して使用することはできません。 :math:`S(\alpha,\beta)` の指定 -------------------------------------------------- 低エネルギー中性子の輸送の際、熱中性子散乱データ :math:`S(\alpha,\beta)` のライブラリーが必要になる場合があります。 特に、熱エネルギー付近で挙動が大きく変わりますのでご注意ください。 :math:`S(\alpha,\beta)` ライブラリーは、 .. code-block:: text MTn 物質ID で指定できます。 ここで、nはその物質の物質番号です。 物質IDは、xsdirに書かれているID番号(lwtr.20tなど)です。 例えば、常温の水の場合は、 .. code-block:: text M1 H 2.0 O 1.0 MT1 lwtr.20t のように指定します。 各データの物質や温度情報は/XS/tsl/tsl-tableに書かれていますので、そちらをご参照下さい。 飛跡構造解析計算用の化学形の指定 -------------------------------------------------- 飛跡構造解析計算を行う場合、物質の化学形が重要です。 電子の飛跡構造解析や、KURBUCを用いた陽子や炭素の飛跡構造解析は、水の断面積を密度によって規格化することによって物質の断面積を計算しますので、化学形の指定は不要です。 しかし、ITSART (Ion Track Structure calculation model for Arbitrary Radiation and Targets)を用いた陽子・イオンの飛跡構造解析の場合、構成元素をもとに断面積やイオン化ポテンシャルを計算することに加え、以下の物質はさらに化合物であることを考慮できます。 :math:`\mathrm{H_2O, CO_2, NH_2, NH_3, SF_6, TeF_6, CH_4, CH_3, }` :math:`\mathrm{C_2H_2, C_2H_4, C_2H_6, C_6H_6,(CH_3)_2NH}` 化合物としての化学形を考慮する場合、 .. code-block:: text M1 H 2.0 O 1.0 chem = H2O のように表記します。 化合物や純物質の混合物の場合、物質量(モル)の比で .. code-block:: text M1 H 1.62 O 0.01 N 1.6 C 0.4 Ar 0.1 chem = H2O 0.01 N2 0.8 CH4 0.4 Ar 0.1 のように表記します。 Chem内の物質量比は自動で規格化されるので、絶対値は重要ではありません。 例題 -------------------------------------------------- 幾つかの例題を下に示します。 [material]セクションの例題(1) .. code-block:: text [ Material ] MAT[ 1 ] 1H 1.0000000E-04 208Pb 1.7238000E-02 204Pb 4.6801000E-04 206Pb 7.9430000E-03 207Pb 7.2838000E-03 MAT[ 2 ] 1H 1.0000000E-09 14N 4.6801000E-05 16O 7.9430000E-06 元素(核種)と組成比の順番は、元素-組成比がデフォルトで、これを逆にしたいときは、 [material]セクションの例題(2) .. code-block:: text [ Material ] den nuc MAT[ 1 ] 1.0000000E-04 1H 1.7238000E-02 208Pb 4.6801000E-04 204Pb 7.9430000E-03 206Pb 7.2838000E-03 207Pb MAT[ 2 ] 1.0000000E-09 1H 4.6801000E-05 14N 7.9430000E-06 16O のように、den nuc を入れて定義します。 [material]セクションの例題(3) .. code-block:: text [ Material ] m1 80196.49c 5.9595d-5 80198.49c 3.9611d-3 80199.49c 6.7025d-3 80200.49c 9.1776d-3 80201.49c 5.2364d-3 80202.49c 1.1863d-2 80204.49c 2.2795d-3 $ ...Be... m3 4009.37c 1.2362E-1 mt3 be.01 $ ...h2o (25C)... m4 1001.37c 6.6658d-2 8016.37c 3.3329d-2 mt4 lwtr.01 $ ...b4c (natural boron; 25%-density)... m5 6012.37c 6.8118d-3 5011.37c 2.1825d-2 $ ...liquid-h2 (20K)... m6 1001.49c 3.1371d-2 1011.49c 1.0457d-2 mt6 orthoh.00 parah.00