24小时热门版块排行榜    

查看: 2725  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 山东省面上项目限额评审 +3 石瑞0426 2026-03-19 3/150 2026-03-21 18:59 by RRRKKK
[考研] 求调剂 +5 十三加油 2026-03-21 5/250 2026-03-21 18:48 by 学员8dgXkO
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 生物学一志愿985,分数349求调剂 +3 zxts12 2026-03-21 3/150 2026-03-21 16:34 by 33来了真来了
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 北科281学硕材料求调剂 +5 tcxiaoxx 2026-03-20 5/250 2026-03-20 21:35 by laoshidan
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
信息提示
请填处理意见