计算化学公社

标题: 求助：使用SCAN泛函，SCF计算无法收敛的问题 [打印本页]

作者
Author: 浅墨 时间: 2023-12-20 09:55
标题: 求助：使用SCAN泛函，SCF计算无法收敛的问题
本帖最后由浅墨于 2023-12-20 09:55 编辑

我想测试一下使用SCAN泛函进行计算，在刚开始跑cp2k的时候，scf收敛很正常，计算了几十步之后，scf计算步数满了也无法到达收敛值，增加scf步数也没用，附件是我的inp，log和ener文件，可以看到在72步之后收敛就出现问题了，想请教下是什么原因，谢谢。
1. 体系：含65个水分子的体相系统
2.inp文件：使用Multiwfn生成
3.赝势文件：使用SCAN泛函对应的赝势，是从hutter的github库中下载的：https://github.com/juerghutter/GTH/blob/master/SCAN/POTENTIAL

作者
Author: 喝杯热水 时间: 2023-12-20 10:49
https://pubs.acs.org/doi/full/10.1021/acs.jpclett.0c02405
这是r2SCAN的论文。里面写了SCAN泛函存在数值收敛困难的问题。

作者
Author: 喝杯热水 时间: 2023-12-20 10:51
主要表现就是SCAN在一些点，会出现剧烈波动的问题。所以SCAN本身数值收敛性就不太好。

作者
Author: 喝杯热水 时间: 2023-12-20 10:53
一定要用metaGGA级别SCAN家族的泛函的话，我建议使用r2SCAN代替，里面对数值收敛性问题解决很好。

作者
Author: 浅墨 时间: 2023-12-20 11:31
好的我用r2SCAN泛函做一下测试，感谢分享

作者
Author: goldNAN 时间: 2023-12-20 18:41
&XC_GRID
XC_DERIV NN10_SMOOTH
XC_SMOOTH_RHO NN10
&END XC_GRID

作者
Author: sobereva 时间: 2023-12-21 03:51
跑这种体系MD用SCAN或r2SCAN没有什么好处，比PBE还贵不少
用PBE-D3(BJ)足矣

Multiwfn里都能直接选r2SCAN，根本没有任何理由用数值稳定性差的SCAN
(, 下载次数 Times of downloads: 13)

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