24小时热门版块排行榜    

查看: 2902  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:18 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 01:04 by ZPa0EcMwuECS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 00:59 by winsaint
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见