specxのパラメータ設定

Posted on : October 07, 2004 (Thu) 20:47:24

by Norihito Suzuki

赤井先生、小谷先生。

たいへんご無沙汰しています。元阿久津研にいた鈴木悟仁です。第2、3回のCMDでもお世話になりました。
最近、酸化ガリウムの計算をしてみようと、チャレンジしていますがエラーが出て止まってしまいます。

XtalEditで次のインプットファイルを作成し、ewidthを2.0に変更しspecxに読み込ませました。
(計算環境は FreeBSD 5.2上にIntel Fortran Compiler 8.0.046で実行ファイルを作成しました。)

go O_Ga.inp
aux
0.500000000000 0.124284546198 0.000000000000
0.500000000000 -0.124284546198 0.000000000000
-0.112318989958 0.000000000000 0.460751013572
23.111348840320

0.0010 2.000 sra vwn nmag
2nd update t 40 0.02000
2
O 1 0. 0.000 2 8 1.00
Ga 1 0. 0.000 3 31 1.00

10
0.090400000000a 0.090400000000b -0.205200000000c Ga
-0.090400000000a -0.090400000000b 0.205200000000c Ga
0.341400000000a 0.341400000000b -0.314300000000c Ga
-0.341400000000a -0.341400000000b 0.314300000000c Ga
0.167400000000a 0.167400000000b 0.101100000000c O
-0.167400000000a -0.167400000000b -0.101100000000c O
0.495700000000a 0.495700000000b 0.253300000000c O
-0.495700000000a -0.495700000000b -0.253300000000c O
0.827900000000a 0.827900000000b 0.436500000000c O
-0.827900000000a -0.827900000000b -0.436500000000c O

するとまず次のエラーが出ました。
***msg in specx...lmx may be too large, check it.

そのまま計算を続けてくれるかと思いきや

:
:
***** self-consistent iteration starts *****
O_6Ga_4
***err in zroots...m too large


と表示した後、止まりました。
まず、最初のspecx.fでのエラーメッセージは138行目を見ると
"lmx" ではなく "mxl" のtypoでしょうか?

zroots.fでのエラーは良くわかりません。
specx.fのパラメータもいろいろいじってみましたが良くわかりませんでした。


natmmx=36, ncmpmx=5, msizmx=132, mxlmx=3, nk1x=24, nk3x=1,
msex=35, ngmx=15, nrpmx=250, ngpmx=250, npmx=200, msr=400

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

Re:parameter setting of specx (in English)

Posted on : October 13, 2004 (Wed) 19:42:34

by Norihito Suzuki

I have cleared some problems, but I still have another problems.

> ***msg in specx...lmx may be too large, check it.

This message was disappeared by setting mxlmx=4. But mxlmx=4 means the program consider
s, p, d, and f-orbital. Despite I consider only Ga and O, I wonder why the case mxlmx=3
caused above message.