vasp自洽计算过程中不收敛的解决办法交流 - 第一性原理 - Vasp&MedeA

回帖置顶 ( 共有3个 )

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫

1ST强帖: 2
应助: 622 (博士)
金币: 6517.3
散金: 1221
红花: 111
帖子: 1345
在线: 517.7小时
虫号: 1730773
注册: 2012-04-01
性别: GG
专业: 理论和计算化学

★ ★
obaica: 金币+1, 3ks 2016-08-24 20:33:23
obaica: 回帖置顶 2016-08-24 22:58:37
蛋蛋小童鞋: 金币+1, 应助指数+1, 谢谢参与 2016-08-25 09:24:16

考虑两种情况:结构不合理，电子步换算法

发自小木虫Android客户端

5楼2016-08-24 20:07:17

obaica

禁言 (著名写手)

obaica: 回帖置顶 2016-08-24 22:57:37
obaica: 取消置顶 2016-08-24 22:58:44
obaica: 回帖置顶 2016-08-24 22:58:51

本帖内容被屏蔽

7楼2016-08-24 22:57:33

hs0218

金虫 (正式写手)

DDr.

1ST强帖: 1
应助: 290 (大学生)
金币: 874.6
红花: 39
帖子: 923
在线: 383.5小时
虫号: 1236662
注册: 2011-03-17
性别: GG
专业: 金属结构材料

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
obaica: 金币+5 2016-08-25 16:28:23
obaica: 回帖置顶 2016-08-25 16:29:11

不管是电子自洽还是原子优化自洽理论上不可能不收敛
原因很多你优化参数选择ISIF 以及AMIX等的选择
还有精度也有关系

路漫漫其修远兮

11楼2016-08-25 16:26:26

回帖支持 ( 显示支持度最高的前 50 名 )

magicmonk

至尊木虫 (著名写手)

应助: 172 (高中生)
金币: 11451.3
红花: 24
帖子: 1124
在线: 3108.8小时
虫号: 1191315
注册: 2011-01-17
专业: 理论和计算化学

★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
obaica: 金币+1, 3ks 2016-08-24 14:32:59
蛋蛋小童鞋: 金币+1, 谢谢参与 2016-08-25 09:25:06

自洽不收敛还真是少见，一般都是结构弛豫不收敛。如果是自洽的话，把波函数打开，再接着跑一次应该就可以了

我们最深的恐惧不是我们能力不够，我们最深的恐惧是能力超越了极限。

4楼2016-08-24 13:48:13

hs0218

金虫 (正式写手)

DDr.

1ST强帖: 1
应助: 290 (大学生)
金币: 874.6
红花: 39
帖子: 923
在线: 383.5小时
虫号: 1236662
注册: 2011-03-17
性别: GG
专业: 金属结构材料

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
obaica: 金币+5 2016-08-25 19:55:58

引用回帖:

12楼: Originally posted by obaica at 2016-08-25 16:29:00
你能够详细说一下：AMIX如何选择，我不太熟悉这个参数。还有ALGO这个该如何选择？谢谢你！...

ALGO=FAST 就可以了

AMIX的默认值是0.2 在某些情况下比较粗的优化时候需要降低到0.1甚至更小
至于ISIF优化，直接用3是很大的计算量，建议先用7找到平衡体积，然后用2基于平衡体积找原子平衡位置

vasp里面INCAR的很多参数都是需要测试的，尤其是针对于自己的材料体系，拿来直接用可以但是可能出现问题，所以多做测试！

还有，你应该多看看vasp manual，那里面对每个参数的解释非常详细。

路漫漫其修远兮

14楼2016-08-25 18:13:52

Feynmanux

至尊木虫 (文坛精英)

应助: 4 (幼儿园)
金币: 13985.3
红花: 2
帖子: 13311
在线: 222.4小时
虫号: 1157895
注册: 2010-11-28
性别: GG
专业: 凝聚态物性 II ：电子结构

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

帮顶，前几天恰好出现了这么个情况。

做优化到只要一个离子步就达到精度结束。isif =2 7 都是如此。

然而到了自洽，却怎么也不收敛，就在一个能量值附近起伏。

IALGO 48 38 都试过了，不收敛；改动 sigma，potim，不收敛。

最后没办法，将K点缩小了，居然收敛了。。。

18楼2016-09-04 13:21:31

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

专家经验: +224
1ST强帖: 2
应助: 47 (小学生)
金币: 3480.8
散金: 14
红花: 223
帖子: 2477
在线: 1253.7小时
虫号: 1369423
注册: 2011-08-16
性别: GG
专业: 金属材料的磨损与磨蚀
管辖: 第一性原理

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
obaica: 金币+5 2016-09-04 19:21:25

如果难收敛首要得考虑模型是否合理，其次看看计算参数选择是否合适、是否经过严格测试，最后实在不行再从算法下手，自洽中的电子步算法换几种之后都不行的话，就得考虑是不是结构本身的问题了。

检查一下之前离子弛豫过程中每一步的电子步是否都收敛了，如果离子弛豫过程中有电子步是没收敛的，之后自洽可能没法正常收敛。

不排除结构不合理的情况，如果结构不合理，离子弛豫前几步可能电子步不是正常收敛的，但是对于vasp，电子步不收敛计算并不会停止（类似于Gaussian的IOp(5/13=1)），这样的结果可能会有些问题。

我遇到过极少情况确实会遇到特定的泛函优化得到的晶格然后切面计算确实收敛不了的，不过这可能是很特殊的情况了。

不一定挂在论坛，计算问题问题欢迎留言。

23楼2016-09-04 14:45:36