24小时热门版块排行榜    

查看: 2869  |  回复: 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一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 04:07 by 6vVgjDL4CnGu
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 KXLV3nuBVBY7 2026-08-07 3/150 2026-08-08 03:27 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见