Re:Error in mkemap(table overflow, in Japanese)

Posted on : November 18, 2004 (Thu) 17:51:45

by Kobayashi Kazuaki

>インプットファイルは次のようなものになります。
>(計算環境は OS:RedHatLinux 8, Compiler:Intel Fortran Compiler 7です)
>
>#----------------------- input data ----------------------------#
># go/ngo/dos/dsp/spc file name
> go data/fe_sp222k12
>
>#- primitive vector --------------------------------------#
># this aux option is only effective for newest versions of AkaiKKR
> aux
> 1.000000000000 0.000000000000 0.000000000000
> 0.000000000000 1.000000000000 0.000000000000
> 0.000000000000 0.000000000000 1.000000000000
> 3.011
>#L.C. of supercell(a.u.)
（中略）
>
>mkemap.fをとりあえず開いてみましたが、原因は特定できませんでした。。
>teble overflowの意味もよくわかりません。
>また、specx.fのパラメータも大きくしたり小さくしたりいろいろやってみましたが、まったくエラーが解消されません。とりあえず、現在のspecx.f中の指定パラメータは以下のとおりです。
> & (natmmx=16, ncmpmx=5, msizmx=150, mxlmx=3, nk1x=100, nk3x=1,
> & msex=35, ngmx=15, nrpmx=250, ngpmx=250, npmx=200, msr=4000)
>
>よければ、アドバイスをいただければと思います。よろしくお願いします。

　筆者は、当該プログラムを使用していないので的外れなもの
かもしれませんが、上記データにある、3.011 (a.u.)がスーパー
セルの格子定数なのでしょうか？。

　鉄(bcc)の非磁性での格子定数は、2.87 Åで、これはa.u.換算で
5 a.u.くらいになります。１６個のスーパーセルならその２倍の値
が必要で、10 a.u.以上となります。それと比べると、3.011 a.u.
（←これが本当にスーパーセルの格子定数なら）というのは非常に
小さ過ぎる値と思われます。もし、鉄１６個の系を3.011 a.u.の
格子定数で計算すれば、通常の条件でのバンド計算は、まず確実
に計算が途中で破綻すると考えられます。

　以上、何かのお役に立てれば幸いです。

(--@--は変なメイル対策です)