计算化学公社

标题: 关于能量震荡不收敛的问题 [打印本页]

作者
Author: lastzealot 时间: 2019-9-8 12:14
标题: 关于能量震荡不收敛的问题
最近在用qe做一些第一性原理计算的工作，经常遇到能量不收敛的情况，增加步数也没效果。
网上查了一下，smear效果不错，但是有人不大建议使用，原因是会存在分数占据，使能量不准确。
所以我想请问下卢老师：在面对scf不收敛的情况，一般都有哪些办法？怎样的选择顺序呢？
谢谢

作者
Author: lastzealot 时间: 2019-9-8 12:39
以前用高斯的时候，看过卢老师关于解决不收敛的帖子，很实用。所以，来请教下在第一性原理计算方面这个问题解决方法。

作者
Author: biogon 时间: 2019-9-8 21:33
增大步数肯定不行的，scf不收敛的博文里面不都说了，虽然程序不一样但是本质一样的，解决方法手册一般都有的

作者
Author: 卡开发发 时间: 2019-9-8 22:10
本帖最后由卡开发发于 2019-9-10 21:56 编辑

lastzealot 发表于 2019-9-8 12:39
以前用高斯的时候，看过卢老师关于解决不收敛的帖子，很实用。所以，来请教下在第一性原理计算方面这个问题 ...

其实可以用QE对照下，有些能实现的可以做一下类比，如：

1、格点方面问题引起的不收敛总可以调整动能截断和密度截断。
2、mixing_ndim，即DIIS大小可以进行调整，设置成1这样DIIS应该就不起作用了。
3、mixing_beta，即damping的权重系数。
4、展宽问题，smear是可以帮助收敛，但固体计算中smear主要的作用不是为了加速SCF，而是为了减少k点使用；另外smear确实会影响电子结构和能量，所以需要严格控制电子熵的大小。体系就有金属性的时候你不用smear显然也不合理，至少力算出来不合理，甚至SCF也难收敛。
5、startingwfc，可以更改波函数初猜的模式，不过我不清楚不同动能截断或k点是否能被读取进来或者有其他的限制，我没这么用过。

另外，mixing_mode没有对应的，我猜测可能是和Broyden混合的参数有关系，对于电子分布极为均匀或极为不均匀的情况也许还可以更换这个选项。

作者
Author: Marisa 时间: 2019-9-9 15:36
体系是否有磁性，磁性体系收敛较为困难

作者
Author: lastzealot 时间: 2019-9-10 21:19

Marisa 发表于 2019-9-9 15:36
体系是否有磁性，磁性体系收敛较为困难

有磁性，MnO2，Mn有一个单电子

作者
Author: lastzealot 时间: 2019-9-10 21:51

卡开发发发表于 2019-9-8 22:10
其实可以用QE对照下，有些能实现的可以做一下类比，如：

1、格点方面问题引起的不收敛总可以调整动能 ...

谢谢大神，非常给力，祝您节日快乐

作者
Author: Marisa 时间: 2019-9-11 15:27

lastzealot 发表于 2019-9-10 21:19
有磁性，MnO2，Mn有一个单电子

https://mp.weixin.qq.com/s?__biz ... ne=23&srcid=#rd
可以看看这个帖子，类推下到QE需要调节那些参数

欢迎光临计算化学公社 (http://ccc.keinsci.com/)