求助：使用SCAN泛函，SCF计算无法收敛的问题

浅墨 · 发表于 Post on 2023-12-20 09:55:12

本帖最后由浅墨于 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

喝杯热水 · 发表于 Post on 2023-12-20 10:49:12

https://pubs.acs.org/doi/full/10.1021/acs.jpclett.0c02405
这是r2SCAN的论文。里面写了SCAN泛函存在数值收敛困难的问题。

喝杯热水 · 发表于 Post on 2023-12-20 10:51:28

主要表现就是SCAN在一些点，会出现剧烈波动的问题。所以SCAN本身数值收敛性就不太好。

喝杯热水 · 发表于 Post on 2023-12-20 10:53:39

一定要用metaGGA级别SCAN家族的泛函的话，我建议使用r2SCAN代替，里面对数值收敛性问题解决很好。

浅墨 · 发表于 Post on 2023-12-20 11:31:57

好的我用r2SCAN泛函做一下测试，感谢分享

goldNAN · 发表于 Post on 2023-12-20 18:41:14

&XC_GRID
XC_DERIV NN10_SMOOTH
XC_SMOOTH_RHO NN10
&END XC_GRID

sobereva · 发表于 Post on 2023-12-21 03:51:22

跑这种体系MD用SCAN或r2SCAN没有什么好处，比PBE还贵不少
用PBE-D3(BJ)足矣

Multiwfn里都能直接选r2SCAN，根本没有任何理由用数值稳定性差的SCAN

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助：使用SCAN泛函，SCF计算无法收敛的问题