Error in mkemap(table overflow)

Posted on : November 17, 2004 (Wed) 23:54:48

by takafumi.ihara

今年春のCMDでお世話になりました。井原と申します。その節はお世話になりました。
Fe（bcc）のスーパーセルの計算を試みたのですが、エラーが出てそこから進みません。

インプットファイルは次のようなものになります。
(計算環境は 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.)
# 1 Angstrom = 0.529177249 a.u.

#- edelt ewidth nrl/sra mjw/vbh/vwn mag/nmag/-mag/rvrs/kick ------#
0.0010 1.0000 sra vwn mag

#- init/1st/2nd quit/update 0/1/2/../t/l/m/h/u iteration pmix -#
2nd update 6 40 0.02000

#- number of type -----------------------------------------------#
1
# type components rtin field l_max Z concentration----------#
2a 1 0. 0.000 2 26 1.00
#Fe

# number of atoms------------------------------------------------#
16
#- position ---------------------------------------------- type -#
0.000000000000a 0.000000000000b 0.000000000000c 2a
0.500000000000a 0.500000000000b 0.500000000000c 2a
0.000000000000a 0.500000000000b 0.500000000000c 2a
0.500000000000a 0.000000000000b 0.500000000000c 2a
0.500000000000a 0.500000000000b 0.000000000000c 2a
0.000000000000a 0.000000000000b 0.500000000000c 2a
0.500000000000a 0.000000000000b 0.000000000000c 2a
0.000000000000a 0.500000000000b 0.000000000000c 2a
0.250000000000a 0.250000000000b 0.250000000000c 2a
0.250000000000a 0.750000000000b 0.750000000000c 2a
0.750000000000a 0.250000000000b 0.750000000000c 2a
0.750000000000a 0.750000000000b 0.250000000000c 2a
0.250000000000a 0.250000000000b 0.750000000000c 2a
0.250000000000a 0.750000000000b 0.250000000000c 2a
0.750000000000a 0.250000000000b 0.250000000000c 2a
0.750000000000a 0.750000000000b 0.750000000000c 2a

出たエラーは次のとおりです。
***err in mkemap...table overflows

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)

よければ、アドバイスをいただければと思います。よろしくお願いします。

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.の
格子定数で計算すれば、通常の条件でのバンド計算は、まず確実
に計算が途中で破綻すると考えられます。

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

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

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

Posted on : November 19, 2004 (Fri) 16:36:28

by takafumi ihara

早速のご指摘ありがとうございます．ご指摘通り，格子定数の設定が小さすぎました．私自身もいろいろなパッケージを使っており，今まで１６個スーパーセルの計算を何回かこなしてきたつもりなのですが，今回何を考えたのか5angstromを約０．５で割って約10a.u.とするべきところを，逆に掛けて今回のような格子定数を設定してしまいました．本当にお恥ずかしい限りです．
若輩者ですが，今後ともご教唆の程をよろしくお願いします．

_/_/_/ 井原隆文(Takafumi Ihara)
_/_/ _/ mailto:tihara0912@hotmail.com