应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » [求助]vasp计算单个Ni原子能量不收敛问题
北京石油化工学院2026年研究生招生接收调剂公告
162/2返回列表上一页12
查看: 2198  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

ypf1449

银虫 (初入文坛)
K点取1X1X1是没有办法收敛的，我计算的另外一个半金属开始设置为1X1X1到2000步都没有收敛，K点增大以后很快就收敛了！！
赞一下 回复此楼
芳圆时态
11楼2008-07-02 10:05:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leongxf

银虫 (初入文坛)

计算结果

采取各位的建议，模拟了很多次，终于收敛了，也找到了不收敛的原因，分析一下，请各位多指教。
1. 刚开始的KEPOINTS，取1 1 1，是因为看到论坛上有贴的计算H原子能的例子，我test了这个列子，发现 1 1 1是收敛的，故Ni原子的也取一样的，后来不收敛，我把k点增加到10 10 10，还是不收敛，所以我分析可能是INCAR的参数设置问题。
根据楼上建议，分别尝试以下的INCAR（以下所有k点都为10 10 10 ）
1. 结果不收敛（使用PW91势）
ISTART = 0 ; ICHARG=2
PREC= n
ISMEAR = -5; SIGMA = 0.2;
ISPIN=1
IBRION=-1
NSW=0
2。收敛
ISTART = 0 ; ICHARG=2
PREC= n
ISMEAR = 0; SIGMA = 0.2;
ISPIN=2
VOSKOWN=1
IBRION=-1
NSW=0
3。收敛（使用PBE势）***********因为有人说用此势比较好，故尝试      
ISTART = 0 ; ICHARG=2
PREC= n
ISMEAR = 0; SIGMA = 0.2;
ISPIN=2
IBRION=-1
NSW=0
4。收敛（使用PW91势）
ISTART = 0 ; ICHARG=2
PREC= n
ISMEAR = 0; SIGMA = 0.2;
ISPIN=1
IBRION=-1
NSW=0

综上所述，以及我另外做的其他尝试，排除了poscar的影响，以及势函数的影响，发现SIGMA = 0.2是收敛的主要因素。当然我也在以上条件3，4 INCAR中其他条件未变的情况下，仅改变SIGMA = 0.01，发现不收敛了。所以我想
SIGMA = 0.2是收敛的主要因素

2，3，4计算中，4没有计算磁矩，2，3考虑自旋，但计算结果发现
E2=-0.31eV
E3=-0.54ev
2，3区别在于势函数不一样，所以采用不同的势函数计算了体材料Ni的单原子能分别为（我拟和了NI的晶格参数，所以结果都为取平衡晶格参数的值）
PW91 -5.416
PBE    -5.454
这样计算的NI的结合能是否为
PW91-E2=-5.106
PBE-E3=-4.914
试验上的数据参考《潜入原子方法理论及其在材料科学中的应用》为4.46ev，所以是不是我的计算还是有问题？还是我的结合能计算方法有问题？
赞一下(1人) 回复此楼
12楼2008-07-03 11:49:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mermaid1128

木虫 (著名写手)
★
wuchenwf(金币+1):thank you 1-27 16:05
k点太小 可以先用所有默认值计算
赞一下(10人) 回复此楼
13楼2008-07-03 14:25:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★
wuli8(金币+1,VIP+0):(*^__^*) 2-19 19:47
我以前计算Ti原子能量遇到不收敛情况，后来我是设置IBRION=2解决的。此外，Ni有磁性，采用自旋极化。
赞一下(10人) 回复此楼
14楼2009-02-19 14:35:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangyc

木虫 (正式写手)
应该和K点选取没有关系吧，增加K点数目只能改进描述原子间的相互作用，单原子计算中需要吗？
赞一下(2人) 回复此楼
15楼2009-02-19 16:16:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilyx8050

禁虫 (正式写手)
本帖内容被屏蔽
16楼2009-02-26 16:40:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leongxf 的主题更新
162/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298求调剂 +4 manman511 2026-04-05 4/200 2026-04-08 16:50 by tjzhao
[考研] 274求调剂求调剂 +10 Jachenbingoo 2026-04-06 13/650 2026-04-08 14:25 by zhq0425
[考研] 求调剂！生物与医药专硕 +9 逆转陆先生 2026-04-01 10/500 2026-04-08 10:37 by maolC
[考研] 338求调剂 +8 wxygxsaaaaa 2026-04-06 8/400 2026-04-08 06:58 by 无际的草原
[考研] 11408 325分 +3 jgtxuxgkx 2026-04-07 3/150 2026-04-07 23:10 by lbsjt
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-04-02 6/300 2026-04-07 18:45 by 求调剂zz
[考研] 工科 22408 267求推荐 +4 wanwan00 2026-04-05 5/250 2026-04-06 22:47 by chenzhimin
[考研] 0703求调剂383分 +9 W55j 2026-04-03 9/450 2026-04-06 06:50 by houyaoxu
[考研] 材料专硕(0856) 339分求调剂 +10 哈哈哈鹅哈哈哈 2026-04-04 10/500 2026-04-05 18:51 by 蓝云思雨
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-04 5/250 2026-04-05 14:06 by imissbao
[考研] 材料与化工306分找调剂 +12 沧海轻舟e 2026-04-03 13/650 2026-04-04 23:45 by lqwchd
[考研] 302求调剂一志愿华中师范大学 +8 小江小江江江 2026-04-02 8/400 2026-04-04 19:50 by 蓝云思雨
[考研] 一志愿沪9，求生物学调剂，326分 +6 刘墨墨 2026-04-04 6/300 2026-04-04 19:44 by 唐沐儿
[考研] 368求调剂 +5 今华习 2026-04-03 7/350 2026-04-04 18:47 by imissbao
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 338求调剂 +4 zzz，，r 2026-04-03 4/200 2026-04-03 16:39 by lijunpoly
[考研] 调剂 +7 祉岷. 2026-04-02 7/350 2026-04-03 09:11 by 花呗还欠600
[考研] 312求调剂 +6 小小墨123 2026-04-02 7/350 2026-04-03 07:32 by jsw79
[考研] 一志愿复旦材料，英一专硕，总分357调剂 +4 1050389037 2026-04-02 5/250 2026-04-02 21:40 by dongzh2009
[考研] 一志愿北京科技材料科学与工程288分，求调剂 +14 是辰啊 2026-04-02 14/700 2026-04-02 21:10 by dongzh2009
信息提示
关闭
请填处理意见
关闭