24小时热门版块排行榜    

查看: 2913  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科又开始会评了,不知道这次命运如何 +6 雨打竹帘 2026-08-30 10/500 2026-08-31 09:29 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 面上意见出来了 +10 黄鸟于飞Chao 2026-08-29 19/950 2026-08-31 08:01 by blueearth171
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 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
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见