PHITSホームページ戻る場合は → https://phits.jaea.go.jp/indexj.html
1. FAQ¶
1.1. パラメータ設定関連¶
- Q1.1: Ver.2.88以前で動いていたインプットファイルがそれ以降のversionで動きません。
A1.1: Ver.2.89以降、ミスを誘発しやすい仕様をいくつか変更しています。基本的には、従来のインプットをそのまま使えるように設計していますが、場合によってはインプットの修正が必要になる場合があります。以下、その例です。
Ver.2.89以降では天然炭素の記号Cとの混同を防ぐため[material]で行頭のc_(_は半角スペース)をコメント文字として利用できません(それ以外のセクションでは引き続き利用可能です)。[material]でエラーが出るようになった場合は、コメント文字として使っているc_を$や#に変更するか、[parameters]セクションにicommat=1と追記してください。
光子や電子の輸送のみを行う場合、従来は中性子の核データがなくても問題ありませんでしたが、Ver.2.93以降は自動的に中性子の核データを利用する設定になっているため"There is no cross section table(s) in xsdir"というエラーが出てしまいます。この場合は、[parameters]にnucdata=0と追記して核データを使わない設定に変更してください。
Ver.2.96以降では、[parameters]セクションが2つ以上ある場合、エラーを出して実行を止めるようにしました。2つ以上定義している場合は、1つに統合してください。
Ver.2.96以降、icntlに応じて不要なタリーを自動的にoffにするようにしました。自動的にoffにされたタリーセクション内にsetやinflコマンドがある場合、それらも無効になりますのでご注意ください。なお、Ver.2.96以降、offされたセクション内にsetやinflコマンドがある場合はwarningが出力されます。
- Q1.2: 電子・陽電子が輸送されません。
A1.2: 電子・陽電子の計算には時間が掛かるためデフォルトでは輸送しない設定になっています。電子・陽電子を輸送したい場合は、negs=1,2としてEGS5モードを利用するようにしてください。
- Q1.3: 核反応モデルに関して、デフォルトの設定は最も精度がよい設定ですか?
A1.3: 基本的にはそうなるように設定していますが、重イオン核反応や高エネルギー核反応により生成される残留核の収率を精度よく計算したい場合は、irqmd=1やismm=1のオプションを利用することを勧めています。ただし、これらのオプションを有効にすると、計算時間が長くなってしまう場合があります。
- Q1.4: イベントジェネレータモードは、どのような場合に使えばよいのでしょうか?
A1.4: イベントジェネレータモード(以下、e-mode)を使った方がよい計算は、検出器の応答関数計算や、半導体ソフトエラー発生率の計算など、イベント毎の情報が必要となる計算です。具体的には、[t-deposit], [t-let], [t-yield]などを使う計算では、e-mode を使った方がよい場合が多いです。逆に、使わない方がよい計算は、中性子束やガンマ線束だけを求める遮蔽計算などです。[t-track], [t-cross]などしか利用しない場合は、一般的にe-mode を使う必要はありません。詳しくはマニュアル5.2.27「Event Generator Mode」をご参照ください。なお、version2.76からイベントジェネレータモードを使う場合の奨励値はe-mode=2となっていますので、ご注意ください。
- Q1.5: 統計誤差導出方法(istdevの設定)はどういう条件で切り替えれば良いですか?
A1.5: バッチ分散モード(istdev=-1 or 1)しか選択できないメモリ共有型並列計算の場合を除き、基本的に全ヒストリー数(maxcas*maxbch)で統計誤差を評価するヒストリー分散モード(istdev=-2 or 2)の方を推奨します。しかしながら、ヒストリー分散モードでは、ヒストリーの計算が終了する毎に誤差情報を記録するため、タリーのメッシュ数が大きい場合(xyzメッシュでnxやnyを数1000ずつ設定するなど)には、計算時間が異常に長くなる可能性があります。そのような条件で計算する際は、バッチ分散モードに切り替えて、maxbchを十分に大きく(10以上)設定するようにしてください。
- Q1.6: JENDL以外の核データをPHITSで使うことは可能でしょうか?
A1.6: ACEフォーマット(MCNPが読み込める形式)であれば可能です。ただし、20 MeV以下の中性子ライブラリとそれ以外で手順が変わります。
- 20 MeV以下の中性子ライブラリの場合
新しく入手した核データに付録されているxsdirファイルをテキストエディタで開き、核データのアドレス情報に関する部分(例えば1001.80c 0.999167 xdata/endf71x/H/1001.710nc 0 1 4 17969 0 0 2.5301E-08)を全てコピーする。
PHITSのdataフォルダ内にあるxsdir.jndをテキストエディタで開き、その180行目付近にあるdirectoryの下に、コピーした内容を貼り付ける。上に書かれている方が優先されるため、最初のJENDL-4.0のアドレス情報(1001.50c 0.999167 neu/H_001.j40n_300K 0 1 1 30720 0 0 2.585E-08)より上にペーストしてください。
PHITSのXSフォルダ内に、上記アドレス情報が一致するように新しく入手した核データファイルをコピーする。例えば、上記のアドレス情報であれば、XSフォルダ内にxdata/endf71x/Hフォルダを作って1001.710ncファイルをその中にコピーする。
- その他のライブラリ(陽子・重陽子・光核反応・高エネルギー中性子ライブラリなど)の場合
新しく入手した核データに付録されているxsdirファイルをテキストエディタで開き、核データのアドレス情報に関する部分(例えば1001.80c 0.999167 xdata/endf71x/H/1001.710nc 0 1 4 17969 0 0 2.5301E-08)を全てPHITSのdataフォルダ内にあるxsdir.jndの末尾にコピー&ペーストする。
PHITSのXSフォルダ内に、上記アドレス情報が一致するように新しく入手した核データファイルをコピーする。例えば、上記のアドレス情報であれば、XSフォルダ内にxdata/endf71x/Hフォルダを作って1001.710ncファイルをその中にコピーする。
[parameters]セクションでlib(i)により、使うライブラリを明示的に指定する。例えば、上記のアドレス情報であれば、lib(2)=80cと指定する。また、使用するライブラリの最大エネルギーをdmax(i)で設定する。なお、指定された核種に対するライブラリーがない場合、警告メッセージが表示された上で、核反応モデルが自動的に利用される。
なお、最新版ENDF(ENDF/B-VIII.0)のACEフォーマットファイルは下記よりダウンロード可能です。 https://nucleardata.lanl.gov/ace/lib80x
この中には、アドレスファイルとしてxsdir及びxsdir_2.0が含まれますが、xsdirの方をお使い下さい。また、xsdir.jndにコピーした後、2行に跨がるデータがある場合は、「+」を「ptable」に一括置換してください。
- Q1.7: 光核反応により生成する粒子(中性子・放射性核種)の統計を効果的にためる方法はありますか?
A1.7: pnimulパラメータを使うことにより、他の反応(コンプトン散乱など)と比べて光核反応を起こしやすくすることができます。例えば、pnimul=2.0とすると、光核反応を引き起こす確率を2倍にする代わりに、光核反応で生成した粒子のウェイトを1/2にします。ウェイトを変化させるので、イベント毎の情報が必要な場合は使えません。また、pnimulパラメータを大きくしすぎると(100以上)、光子・電子の挙動がずれてきてしまいますので、本パラメータを使う場合は、光子・電子の挙動(フラックス・線量など)が大きく変動しないことを確認しながら使ってください。
- Q1.8: 再開始計算を行おうとすると、タリーを変更していないにも関わらず"Error: inconsistent tally parameters"というエラーが出て止まります。対処方法はありますか?
A1.8: この問題は、桁落ちの関係で新旧タリーのパラメータが異なるとPHITSが認識してしまうことにより発生します。タリーで使うパラメータを複雑な数式などで定義した場合によく発生します。この問題が生じた場合は[parameters]セクションでireschk = 1と設定し、タリーの一致性をチェックしないモードで再開始計算を実行して下さい。
1.2. エラー、コンパイル関連¶
- Q2.1: コンパイルがうまくできません。
A2.1: PHITSのコンパイルに関しては、基本的にはマニュアル「10. コンパイル」をご参照ください。PHITSは、Intel Fortran 11.1以降及びgfortran 4.7以降でコンパイル可能です。OneAPIのIntel Fortranコンパイラなら"phits/document/OneAPI-Fortran-manual.pdf" の通り、無償で入手できます。なお、PHITS 事務局が奨励するコンパイラは、計算速度の観点から Intel fortranです(Q2.6 参照)。
- Q2.2: 実行中にAccess violationもしくはSegmentation faultというエラーが起きてしまいます。
A2.2: PHITSのメモリ関連バグの可能性が高いです。問題が発生したインプットファイルの1行目に$DBG=1と追記してPHITSをデバッグモードで実行してください。何かエラーメッセージが出たら、そのメッセージを実行したPC環境の情報とともにPHITS事務局までご連絡ください。自分でコンパイルしたPHITS実行ファイルをお使いの方は、デバッグモード(makefile のDEBUG? = falseをtrueに変更して再コンパイル)で実行して、そのエラーメッセージをお教えください。
- Q2.3: 外部ファイルを取り込むためのコマンド「infl:」が使えません。
A2.3: infl:コマンドを使うためには、インプットファイルの1行目に「file=インプットファイル名」と書いておく必要があります。もしくは、通常のPHITS 用インプットファイルの他に、別のファイルを準備し、そのファイルの先頭行に「file=通常PHITS のインプットファイル名」と書いて実行することも可能です。
- Q2.4: Windowsでは実行できるインプットファイルがLINUX やUNIXでは実行できません。
A2.4: いろいろな原因が考えられますが、WindowsとLINUX系(UNIX含む)では、改行コードが違いファイルの転送に失敗している可能性があります。通常、WindowsからLINUX系にファイルを転送する際、FTPなどのプロトコルを使うと思いますが、その転送の際、アスキーモードとバイナリーモードというのがあり、アスキーファイル(PHITS のインプットファイルなど)は、アスキーモードで転送しないと、ファイルが破損してしまう可能性があります。お使いの FTPソフトの設定をご確認ください。なお、Windowsで作られたバイナリファイル(PHITS の実行ファイルなど)をLINUX系に転送しても、通常、正しく動作しません。
- Q2.5: Cygwin上でPHITS は動作しますか?
A2.5: はい。PHITS のmakefileにあるCygwinの項目をご参照ください。
- Q2.6: Intel fortranとgfortran、どちらがPHITSをコンパイルするのに適していますか?
A2.6: 基本的には、Intel fortranを使うことをお勧めいたします。理由は、計算速度が速いためです。gfortranでは、最適化オプションを付けると正しく動作しない場合があるため、奨励の最適化オプションは-O0(すなわち最適化しない)です。したがって、平均で3-5倍くらいの計算速度の差があると考えています。Intel fortranは2020年頃以降、oneAPIツールキットとして無償で提供されています。
- Q2.7: メモリ分散型並列(MPI)と共有型並列(OpenMP)はどのように使い分けるのがよいでしょうか?
A2.7: MPIプロトコルがインストールされたコンピュータであれば、基本的にMPIの方が計算速度は速くなります。ただし、通常のWindowsやMacにはMPIプロトコルがインストールされていませんので、簡単に並列計算を実行したいのであれば、OpenMPを使うことをお勧めしています。また、ボクセルファントムや高エネルギー核データライブラリなどの巨大なメモリを使用する計算では、MPIだとメモリ容量をオーバーしてしまう場合があります。その場合は、MPIとOpenMPを組み合わせたハイブリッド並列を実施する必要があります。その際、メモリの許す限りできるだけMPI並列数を増やした方が計算速度は速くなります。ただし、ジョブコントロール用にもOpenMP並列数分のコアを使用しますので、MPI並列数の設定にはご注意ください。例えば、128コアの計算機を使ってOpenMP並列数8で計算する場合、指定できるMPI並列数は、128/8-1 = 15、となります。
- Q2.8: DCHAINを実行した際に、「指定されたファイルが見つかりません」といったエラーが表示されます。
A2.8: DCHAINのインプットファイルやinflコマンドを含むPHITSのインプットファイルが全角文字を含むフォルダ名の中に置かれている場合、適切に動作しません。その場合は、半角文字のみのフォルダ名に変更してください。
- Q2.9: Lost particleが頻発します。どう対処したらよいでしょうか?
A2.9: いろいろな理由が考えられますが、最も多いのは、ジオメトリに2重定義や未定義領域がある場合です。[t-gshow]やicntl=8オプションを使ってジオメトリを描画し、黒や紫の部分がないか確認してください。また、ジオメトリが正しく定義されていても、以下のケースに該当する場合はLost particleが頻発する可能性がありますので、入力ファイルの修正が必要となります。
- 線源発生位置もしくは線源の進行方向とセルの境界面が完全に一致している場合。
この場合、線源の発生位置に少し幅を持たせてください。例えば、x=0に線源とセル境界面がある場合、線源の発生位置をx0=x1=0とするのではなく、x0=0, x1=0.001のようにしてください。
- Lattice構造を[transform]で回転させ、かつLattice構造を入れる外枠セルとLattice構造そのものの大きさが完全に一致している場合。
この場合、Lattice構造を入れる外枠セルを少しだけ小さくしてください。詳細は"phits/lecture/advanced/voxel"をご参照ください。
- [weight window]のxyzメッシュ境界面とセルの境界面が完全に一致している場合。
この場合、[weight window]のxyzメッシュを少しずらして定義してください。また、[t-wwg]を使って[weight window]を作成する場合は、[t-wwg]で定義するxyzメッシュを少しずらしてください。特に、パラメータを使ってジオメトリとxyzメッシュ範囲を連動させている場合はご注意ください。なお、version 3.34より、[weight window]及び[t-wwg]のxyzメッシュ最大値・最小値が整数で定義された場合は、その整数から最小値は-1*deltxyz、最大値は+2*deltxyzずらすようにしています。
なお、上記いずれにも該当しない場合でも、桁落ちの関係でどうしてもLost particleが出てしまう場合があります。その場合、おおよその目安として、数1000ヒストリに1回程度の頻度であれば、無視して問題ありません。それ以上の頻度であれば、PHITS事務局までご連絡ください。
- Q2.10: LinuxでPHITSのOpenMP並列計算ができません。ライブラリlibiomp5.soが必要というメッセージが出ます。
A2.10: OpenMPによる並列計算のためにはlibiomp5.soのライブラリのインストールが必要になる場合があります。この場合、各ディストリビューションの手順に従いlibomp-devをインストールしてください。例えばUbuntuの場合は、
sudo apt-get install libomp-dev
でインストールできます。これにより、/usr/lib/x86_64-linux-gnu/libomp.so.5のライブラリがインストールされます。次に、このディレクトリに移動
cd /usr/lib/x86_64-linux-gnu/
し、"libiomp5.so にシンボリックリンクを作ります。
sudo ln -s libomp.so.5 libiomp5.so
これにより、"libiomp5.so のライブラリが認識されるようになり、PHITSによるメモリ共有型並列計算が実行できるようになります。
- Q2.11: LinuxでPHITSを実行しようとすると、GLIBC_2.34 not foundの様なエラーがでます。
A2.11: PHITSパッケージに付属しているLinux用の実行ファイルは、新しいバージョンのGLIBCライブラリを使用しています。このため、古いLinuxのOS(例: Ubuntu 20.04LTS)では、このようなエラーが出て実行できません。この場合、OSのアップデートを検討していただくか、ご自身の環境で、PHITSのソースを再コンパイルしていただく必要があります。再コンパイルにはFortranコンパイラが必要となりますので、Intel FortranやGfortranなどのFortranコンパイラを用意してください。PHITSのソースコードはphits/srcフォルダに含まれています。makefileの中のENVFLAGSの値をご自身の環境に合わせて変更し、makeコマンドを実行することで、簡単にPHITSの再コンパイルをすることができますマニュアル「10. コンパイル」参照)。
- Q2.12: 計算の異常終了や無限ループなど、明らかにPHITSのバグと思われる症状が発生します。どうしたら良いでしょうか?
A2.12: 問題が発生したインプットファイルの1行目に$DBG=1と追記してPHITSをデバッグモードで実行してください。何かエラーメッセージが出たら、そのメッセージ、実行したPC環境の情報、可能ならば実行に必要な入力ファイル一式をPHITS事務局までお送りください。その際、無限ループの場合は以下の手順に従って1ヒストリーで問題を再現する初期乱数を特定してからファイルをお送りいただけると大変ありがたいです。
PHITS が無限ループに陥ったときに生成された batch.outを開く。
"next initial random seed"と書かれた次の行(例えば以下)をコピーし、入力ファイルの[parameters]セクションに貼り付け。
bitrseed = 1111111100000000100001111000001011001100011000101101000011001000
[parameters]セクション内でidmprijk=1を設定。
編集した入力ファイルを用いてPHITSを再実行。
rijkdmp.inpファイル(MPI並列時はrijkdmp***.inp、***はPE番号)に書き込まれたbitrseedをコピーし、[parameters]セクションに貼り付け。
再度PHITSを実行し、1ヒストリー目で無限ループに陥るかどうか確認。
1.3. タリー関連¶
- Q3.1: [t-heat]と[t-deposit]は何が違うのでしょうか?
A3.1: 共に放射線による付与エネルギー(発熱)を出力するタリーです。ただし、[t-heat]でしか実行できなかったカーマ近似による計算は、version 3.05以降、[t-deposit]でも可能となったため、現在[t-heat]の利用は奨励していません。
- Q3.2: [t-track]で重イオンのtrack length を計算すると、結果が不自然なのですが。
A3.2: 重イオンのエネルギーは、入射エネルギーや切断エネルギーに対しては核子あたりのエネルギー[MeV/n]で指定しますが、タリーの中でエネルギーを指定する場合は、核子あたりではなく、全エネルギー[MeV]で指定します。ただし、[parameters]においてiMeVperu=1とすることで、タリーにおけるエネルギーの単位をMeV/nに変えることができます。
- Q3.3: [t-let]、[t-sed]の結果が不自然なのですが。
A3.3: letmat でLET を計算する物質を指定する際、その物質がPHITS の体系内で使われている場合、密度はその体系内で使用している密度となります。例えば、水蒸気内の発熱量に対するLET 分布を、液体の水(1 \({\rm g/cm}^3\))に対するLETの関数として計算したい場合、[ Material ]で2つの水を定義し、1つをPHITSの体系内で使用し、もう1つをletmat で使用する必要があります。また、letmat で使用する物質の密度は、原子数密度の絶対値(1H 6.893407e-2 16O 3.446704e-2)を指定してください。詳しくは、奨励設定ファイルのParticleTherapyをご参照下さい。また、[t-let]の場合、LETメッシュ幅を細かく設定しすぎると結果がガタガタしてしまいますので、あまり細かいメッシュ幅を設定しないようにしてください(1桁10メッシュ程度まで)。
- Q3.4: タリー結果の誤差はどのように計算されているのでしょうか?
A3.4: バージョン2.50より、正しく計算できるように修正しました。詳しくはマニュアル「5.2.2 ヒストリー数、配列の大きさ」をご参照ください。
- Q3.5: 並列計算でdumpモードは使えるでしょうか?
A3.5: 基本的には使えます。ただし、メモリ共有型(OpenMP)並列の場合は、1つのファイルに全てのノードからの情報を書き込むため、書き込む順番がヒストリー番号毎に並ばなくなってしまいます。一方、メモリ分散型(MPI)並列の場合は、dumpファイル名にコア番号付けたファイルを出力しますので、そのような問題は起きません。
- Q3.6: mesh = regとmesh = xyzを用いて同じ直方体内のフラックスや発熱量を計算した場合に、それぞれの結果が異なるのですが。
A3.6: タリー結果を単位体積あたりに規格化して出力する場合、xyzやr-zメッシュではタリー領域の体積を自動的に計算して出力しますが、regメッシュの場合は、[volume]セクションで体積を指定しない限り体積を1 \({\rm cm}^3\) として出力します。例えば、mesh = regで[volume]セクションがない場合、[t-track]のunit=1[1/cm2/source]とunit=4[cm/source]は同じ値を与えます。
- Q3.7: [t-deposit]でoutput=depositとした場合に入射エネルギーよりも高いイベントが発生するのはなぜでしょうか?
A3.7: 発熱する核反応(Q値が正)が発生すると、入射粒子よりも2次粒子の合計エネルギーが高くなるため、入射エネルギーよりも付与エネルギーが大きくなる場合があります。また、xyzメッシュもしくはr-zメッシュかつエネルギー分散を考慮した場合(nedisp \(\ne\) 0)、エネルギー分散の考慮方法の仕様により、イベント毎の情報が正しく計算できなくなってしまいます。そのような場合は、mesh = regとして、[cell]セクションで領域を細かく分割して計算してください。
- Q3.8: axis = xzで表示する際、横軸と縦軸を入れ替えることは可能でしょうか?
A3.8: 可能です。まず、該当タリーのx-typeとz-typeに関連するパラメータを入れ替えます。次に下記3行をそのタリーに追加します。
trcl = 0 0 0 0 0 1 0 1 0 1 0 0 x-txt = x [cm] y-txt = z [cm]
これで図が90度回転します。
- Q3.9: タリーを増やすと計算時間が長くなるのですが対策方法はありますか?
A3.9: メッシュ数の多い場合などメモリ使用量の多いタリーを使うと、ヒストリー毎のメモリ初期化の計算時間がボトルネックになる場合があります。その場合は、istdev=1(バッチ分散)やitalsh=1(タリーメモリの共有化、OpenMPのみ)を試してみてください。また、version 2.94から3.24の間は、[t-dchain]及び[t-yield]を使うと計算時間が最大で5倍くらいになってしまいますので、これらのタリーを使った計算で計算時間にお困りの方は、PHITS事務局までご連絡ください。
- Q3.10: Sum overの値に対する統計誤差が知りたいのですが、簡単に出力できるでしょうか?
A3.10: バージョン3.35より、計算できるようになりました。
- Q3.11: stdcutを使って計算を打ちきる際、全てのメッシュではなく特定のメッシュに対する統計誤差のみチェックする方法はあるでしょうか?
A3.11: 残念ながらありません。一部のデータのみチェックしたい場合は、メッシュをチェックしたい領域に限定した新たなタリーを設定して、そちらにstdcutを定義してください。
- Q3.12: PHITSのタリーとMCNPのタリーとの対応を教えてください。
A3.12: 以下の通りです。
F1: Surface Current → [t-cross] with output=current
F2: Surface Flux → [t-cross] with output=flux
F4: Track Length cell flux → [t-track]
F4 with FM card: Reaction number → [t-interact],[t-yield],[t-product]
F5: Flux at a point and ring → [t-point]
F6: Energy Deposition → [t-deposit] with output=dose
F8: Pulse Height → [t-deposit] with output=deposit
- Q3.13: [t-cross]でmesh=regで計算した結果と同じ領域をmesh=xyzやr-zで定義して計算した結果が異なるのですが。
A3.13: [t-cross]でmesh=regの場合は、タリー面の表面積(area)を定義しないと絶対値が \({\rm /cm}^2\) に規格化されません。また、mesh=r-zの時にメッシュの側面が円柱状のあるセルの側面と一致する場合や、mesh=xyzの時に座標変換によって存在するセルの境界面とタリー面が一致する場合は、桁落ち誤差により数え落としが発生してしまうので、そのような場合はメッシュサイズを微妙に変化させて境界面と一致しないようにしてください。
1.4. 線源設定関連¶
- Q4.1: 等方線源ソース(s-type=9, dir = -all)のときの規格化はどのようにすればよいですか?
A4.1: PHITS で作る等方線源は、指定した球の中でフラックスが一定値: \(1/\pi/r_1^2\) ( \(r_1\) は球の内径)となるように作られています。したがって、PHITS の結果(/source)を単位フルエンス当たりに変換するには、結果に \(\pi r_1^2\) を乗じてください。また、等方線源の作成にはウェイトの概念を利用していますので、検出器の応答関数などイベント毎の分布計算には利用できませんので、ご注意ください。そのような計算が必要な方は、dir = isoとしてください。
- Q4.2: PHITS で認識する線源領域が、実際の領域とずれるのですが。
A4.2: 面線源や点線源を定義する際、その面や点が、Cellで使うsurfaceと完璧に一致している場合、PHITS の中で認識する線源領域が実際の領域とずれる場合があります。このような状態を避けるため、面線源や点線源は、cellで使う面とは少しずらして作成するようにしてください。また、特定の面に完璧に平行なビームを発生させた場合も同様のエラーが起きる可能性がありますのでご注意下さい。
- Q4.3: ADVANTGが作成する線源エネルギースペクトルバイアス関数(MCNPのSBカード)をPHITSで利用することは可能ですか?
A4.3: 可能です。下記の例のように、再現したいスペクトルをe-type=4で定義した上、p-typeで各エネルギービンから実際に発生させる粒子数比を定義してください。
e-type=4 ne = 3 0.0 2.0 1.0 1.5 2.0 1.0 3 p-type=1 1.0 2.0 3.0
1.5. PHIG-3D関連¶
- Q5.1: WindowsでPHIG-3Dを起動しようとすると"VCOMP140.DLLが見つからない"というエラーが出て止まってしまいます。
A5.1: Microsoft Visual C++ 2015 再頒布可能パッケージを次のウェブページからダウンロードしてインストールしてみてください。 https://www.microsoft.com/en-US/download/details.aspx?id=53840 これでもエラーが解消されない場合は、事務局までご連絡ください。
- Q5.2: MacでPHIG-3Dを起動しようとすると"開発元を検証できないため開けません"というエラーが出てしまいます。
A5.2: コントロールキーを押しながら右クリックしてショートカットメニューの「開く」を選んで起動してください。また、「システム環境設定」->「セキュリティとプライバシー」から許可を与えることも可能です。
- Q5.3: 描画ボタンを押すと"Not enough free memory"というエラーが出ます。
A5.3: PHITS本体を含む他のアプリなどがメモリを大量に使用している可能性がありますので、もし他のアプリが動いているならそれらを停止してから再度試してみてください。それでもうまく行かない場合は、左のタブにある「Setting」をクリックして「解像度(サンプリングレート)」の「点数/セル」の数値を小さくしてください。
1.6. その他¶
- Q6.1: 規制庁への加速器施設や廃棄物保管庫等の許認可申請において、PHITSによる遮蔽計算の結果を使用することができるでしょうか?
A6.1: RI法の申請においては、PHITSが認められています。 J-PARCをはじめとして、放医研の新治療棟、 北大の陽子線治療施設、神奈川がんセンターの重粒子線治療施設等で、実績が多数あります。一方、炉規法に関しては、申請において認められるかどうかに未確認です。許可取得例をご存知でしたら、事務局までご一報ください。これまでにユーザーの皆様からいただいた情報は、次のウェブページにまとめています。 https://phits.jaea.go.jp/image/shinsei-PHITS.pdf
- Q6.2: PHITSで得られた解析結果を商用目的で利用することは可能でしょうか?
A6.2: 基本的には可能です。例えば、PHITSの解析結果を元に新しい検出器を開発して販売することや、外部機関から加速器の遮蔽設計を請け負って対価を得ることも問題ありません。ただし、PHITSを含む商品(例えばPHITSの実行ファイルが含まれる治療計画システムなど)を開発して販売する場合は、原子力機構と商用利用契約を結んでいただく必要があります。そのような計画がある場合は、早めに事務局までお知らせ下さい。
- Q6.3: 放射線に関連する物理現象でPHITSで扱えないものは何でしょうか?
A6.3: PHITSは個々の放射線の挙動をモンテカルロ法を用いて解析するため、複数の放射線による影響(電流や放射線自身が作る電場の影響など)を考慮することができません。また、物理過程の後に引き起こされる化学過程や材質の変化は模擬することができません。詳しくは、マニュアル「Appendix B PHITSで扱えない物理現象リスト」をご参照ください。
- Q6.4: インストーラを削除してしまったため、新しいPCにPHITSをインストールできません。どうしたらよいでしょうか?
A6.4: ホームページの「連絡・質問」ページから「新バージョン請求」を選択して「使用しているPHITSのバージョン」欄に「新規インストール希望」と書いてPHITS事務局にご連絡ください。ユーザー登録の有無を確認した上で、登録済みであれば最新版を無料でお送りいたします。
- Q6.5: 大規模な並列計算を実行したいのですが、PHITSがインストールされた共用のスパコンはありますか?
A6.5: PHITSがインストールされた共用スパコンは、公益財団法人 計算科学振興財団(通称FOCUS)(https://www.j-focus.or.jp/)より有償で利用可能です。なお、このスパコンにインストールされたPHITSを使う場合でも、PHITSのユーザー登録は各自で必要となります。