24小时热门版块排行榜    

查看: 2678  |  回复: 6

hlgzhl

金虫 (小有名气)

[求助] vasp计算单个铁原子能量

大家好,小弟现在用vasp计算单个Fe原子能量出现问题,incar文件是按照说明书设置的,并且设置ISPIN=2,POSCAR选取15*15*15的胞,把原子Fe原子放进去,K点取1。但每次提交任务后,根本不进行循环计算,也不提示任何错误信息。已经对incar文件调试了,均无效。
特此想大家询问一下,是否有遇到过这种情况,是如何解决的呢?
希望得到好心人的帮助!

发自小木虫Android客户端
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

西风月123

木虫 (著名写手)

我前两天也算了,单个Fe原子放在一个10A的立方盒子中,我用的赝势就是Fe不加任何尾缀的,K点用的也只1*1*1的,但我当时发现的结果就是开自旋ISPIN=2时,电子步好像要跑多于60多步才能收敛。
2楼2016-11-27 20:53:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hlgzhl

金虫 (小有名气)

引用回帖:
2楼: Originally posted by 西风月123 at 2016-11-27 20:53:49
我前两天也算了,单个Fe原子放在一个10A的立方盒子中,我用的赝势就是Fe不加任何尾缀的,K点用的也只1*1*1的,但我当时发现的结果就是开自旋ISPIN=2时,电子步好像要跑多于60多步才能收敛。

谢谢,可能Fe原子能量不容易收敛吧。可我的压根就没进行运算,直接停了。能看一下你的incar文件吗,多谢啊。

发自小木虫Android客户端
3楼2016-11-27 21:06:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

西风月123

木虫 (著名写手)

引用回帖:
3楼: Originally posted by hlgzhl at 2016-11-27 21:06:41
谢谢,可能Fe原子能量不容易收敛吧。可我的压根就没进行运算,直接停了。能看一下你的incar文件吗,多谢啊。
...

我明天给你贴上去吧,包括最后的能量, 我们可以一起交流
4楼2016-11-27 21:11:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hlgzhl

金虫 (小有名气)

引用回帖:
4楼: Originally posted by 西风月123 at 2016-11-27 21:11:41
我明天给你贴上去吧,包括最后的能量, 我们可以一起交流...

好的,多麻烦你了。

发自小木虫Android客户端
5楼2016-11-27 21:44:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

西风月123

木虫 (著名写手)

引用回帖:
5楼: Originally posted by hlgzhl at 2016-11-27 21:44:55
好的,多麻烦你了。
...

INCAR:
SYSTEM=ATOM

NCORE=4

PREC=H
LREAL=.FALSE.
ALGO=FAST
#ENCUT=520

ISTART=0
ICHARG=2

#NSW=120
NELM=200
LMAXMIX=4
EDIFF=1E-6
#EDIFFG=-0.01
ISMEAR=0
SIGMA=0.05
LORBIT=11
NEDOS=2000
#POTIM=0.1
#ISIF=3
#IBRION=2

ISYM=0
ISPIN=2


#LDAU= .TRUE.
#LDAUTYPE= 2
#LDAUL= -1    2
#LDAUU=  0    3.9
#LDAUJ=  0    0  
#LDAUPRINT= 2

#LWAVE=F
不用优化,最后energy  without entropy=       -3.43252748
你自己试一下,有问题可以再交流。
算出结果来,麻烦把总能贴出来看一下,谢谢!
6楼2016-11-28 16:45:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hlgzhl

金虫 (小有名气)

引用回帖:
6楼: Originally posted by 西风月123 at 2016-11-28 16:45:16
INCAR:
SYSTEM=ATOM

NCORE=4

PREC=H
LREAL=.FALSE.
ALGO=FAST
#ENCUT=520

ISTART=0
ICHARG=2

#NSW=120
NELM=200
LMAXMIX=4
EDIFF=1E-6
#EDIFFG=-0.01
ISMEAR=0
SIGMA=0.05
LORBIT=11
NED ...

很抱歉,最近比较忙。没来的及回复。我计算的结果是-3.395494eV。谢谢您的热心帮助!

发自小木虫Android客户端
7楼2016-12-14 15:25:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hlgzhl 的主题更新
信息提示
请填处理意见