小木虫 --- 700万学术达人喜爱的学术科研平台

首页 >> 第一原理 >>【求助】Wien2k计算孤立原子,急急急!!!!

【求助】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
。。。。。。
。。。。。。
。。。。。。(省略)

 返回小木虫查看更多

今日热帖

用户评论

人呢?:(:o

个人建议,最好建个20*20*20埃的超胞,把孤立原子放进去算即可

Originally posted by zdhlover at 2009-9-11 17:06:
个人建议,最好建个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原子测试一下,看看怎样能正确自洽收敛。谢谢。:),

不懂,帮不上忙,帮顶。但是我有个疑问:
为什么要算单个原子啊,不知道算出来的东西是否又意义?

Originally posted by zeoliters at 2009-9-12 16:00:
不懂,帮不上忙,帮顶。但是我有个疑问:
为什么要算单个原子啊,不知道算出来的东西是否又意义?
我主要是想算结合能, binding energy or cohesive energy.
不知您有何高见,请赐教。:)

猜你喜欢

学术必备
与600万学术达人在线互动!


扫描下载送金币