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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » [求助]vasp计算单个Ni原子能量不收敛问题
162/2返回列表上一页12
查看: 2656  |  回复: 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
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 14/700 2026-11-01 11:33 by Equinoxhua
[基金申请] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 11:09 by BCdvBAEpt3zm
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:56 by BCdvBAEpt3zm
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 10:53 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 10:27 by BCdvBAEpt3zm
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 X5bHxOULxFqb 2026-10-30 26/1300 2026-11-01 10:22 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 07:01 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 17/850 2026-11-01 06:48 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 3/150 2026-11-01 04:04 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 03:51 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
信息提示
关闭
请填处理意见
关闭