PHASE FAQ (Frequently Asked Questions)
1.インストール、実行
2.電子状態計算
- SCF 計算の途中で，電荷密度が負になるというエラーがでる
- cutoff wf やcutoff cd の値はどうやって決めればよいのですか？
- 入力ファイルでxctype を指定しますが，擬ポテンシャルデータ中にあるxctype と一致している必要がありますか？
- SCF の収束の条件を満たしていないように見えるのに，計算が終了するのはなぜですか。
- 計算結果が本当に収束したかどうかは，出力ファイルのどこを見れば確認できるのでしょうか。
1.インストール、実行
.rhost にマシン名とユーザーアカウント名とを書いていますか。
2.電子状態計算
(例：** WARN CHG.DEN < 0.0 AT 1 - 0.8220751D-1)
charge mesh の精度は充分でしょうか？電荷密度のカットオフを指定するcutoff cd の値を大きくして計算してみましょう。
(なお，SCF 計算の初期の段階にのみ，上記warning が現れる場合は，問題ありません。)
cutoff wf の値は用いる擬ポテンシャルに依存します。通常は，結晶の計算をし，格子定数や，体積弾性率の実験値を再現する様に値を決めます。通常のTM ポテンシャルの場合はcutoff cd ＝ 4 cutoff wf となります。ただし，擬ポテンシャルがPCC を使っている場合やウルトラソフトポテンシャルの場合は，一般にそれよりも大きくする必要があります。
合理的な計算をするために，原則的に必要です。入力でxctype を指定しなかった場合，擬ポテンシャルのファイル中に書かれているxctype の値が使われます。尚，複数の擬ポテンシャルが使われる場合，各擬ポテンシャルデータ中に書かれているxctype の種類が一致していないと，job は実行されません。
収束条件は，1 原子当たりのエネルギーに対して設定されています。これに対して，PHASE の標準出力結果であるoutput000 に出力されているのは全エネルギーなので，これを原子数で割ったものが，収束条件を満たしているかどうかを確認しなければなりません。
SCF 計算の結果が収束していれば，出力ファイルであるnfefn.data やnfdynm.data に何か書き込まれているはずです。原子に働く力の計算が収束したならば、continue.data の最後に，convergenceという記述が見つかるはずです。