【求助】Wien2k计算孤立原子，急急急！！！！
请问各位虫友：
怎样用wien2k去计算孤立原子的性质，如孤立原子的能量等。
我的struct文件：Fe-atom.struct
我用这样的struct文件（附件为结构文件，大家帮忙算算看），我算了不收敛，请大家帮我看看，这样算孤立原子对吗？（我知道vasp和siesta这样做可以）
Fe-atom
P LATTICE,NONEQUIV.ATOMS: 1
MODE OF CALC=RELA unit=ang
28.345903 28.345903 28.345903 90.000000 90.000000 90.000000
ATOM 1: X=0.00000000 Y=0.00000000 Z=0.00000000
MULT= 1 ISPLIT= 2
Fe NPT= 781 R0=0.00005000 RMT= 14.0000 Z: 26.0
LOCAL ROT MATRIX: 1.0000000 0.0000000 0.0000000
0.0000000 1.0000000 0.0000000
0.0000000 0.0000000 1.0000000
48 NUMBER OF SYMMETRY OPERATIONS
-1 0 0 0.00000000
0-1 0 0.00000000
0 0-1 0.00000000
。。。。。。
。。。。。。
。。。。。。（省略）
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
个人建议，最好建个20*20*20埃的超胞，把孤立原子放进去算即可
我前面已经说了，我现在已经按照15*15*15 Ang的超包在算了，但是最后自洽循环三次次次就退出了。提示如下：
cycle 3 (Wed Sep 9 16:44:41 CST 2009) (38/97 to go)
> lapw0 (16:44:41) 73.504u 5.448s 1:19.03 99.8% 0+0k 0+0io 0pf+0w
> lapw1 -up (16:46:00) 0.196u 0.040s 0:00.23 100.0% 0+0k 0+0io 0pf+0w
> lapw1 -dn (16:46:01) 0.196u 0.040s 0:00.23 100.0% 0+0k 0+0io 0pf+0w
> lapw2 -up (16:46:01) WARNING: EF not accurate, new emin,emax,NE-min,NE-max -8.398185287026881E-002
-8.398185096983876E-002 13.0000000000260 15.0000000000300
4.480u 0.120s 0:04.60 100.0% 0+0k 0+0io 0pf+0w
> stop error
我的问题是：请您帮我把附件中的Fe原子测试一下，看看怎样能正确自洽收敛。谢谢
，
不懂，帮不上忙，帮顶。但是我有个疑问：
为什么要算单个原子啊，不知道算出来的东西是否又意义？
我主要是想算结合能, binding energy or cohesive energy.
不知您有何高见，请赐教。