Error (supercell3*3*3 Calculation)

Posted on : November 24, 2004 (Wed) 21:58:08

by takafumi ihara

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

インプットファイルは次のようなものになります。
(計算環境は OS:RedHatLinux 8, Compiler:Intel Fortran Compiler 7です)

#----------------------- input data ----------------------------#
# go/ngo/dos/dsp/spc file name
go data/05

#- 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
16.129
#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 5 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------------------------------------------------#
54
#- position ---------------------------------------------- type -#
0.00000000a 0.00000000b 0.00000000c 2a
0.33333333a 0.00000000b 0.00000000c 2a
＜以下，原子54個(typeは全て2a)分の座標の続き；略＞

出たエラーは次のとおりです。
セグメントエラー spec3 < in/05 > out/05 (coreを出力しました)

k点の数を減らしてみたり，specx.fのパラメータも大きくしたり小さくしたりいろいろやってみましたが、まったくエラーが解消されません。とりあえず、現在のspec*.f中の指定パラメータは以下のとおりです。
& (natmmx=54, ncmpmx=1, msizmx=150, mxlmx=3, nk1x=40, nk3x=21,
& msex=201, ngmx=15, nrpmx=250, ngpmx=250, npmx=100, msr=400)

身勝手な質問だとは承知していますが，よろしければアドバイスをいただければと思います。よろしくお願いします。