计算化学公社

标题: 单点能计算时conver=6和conver=8的结果能否比较 [打印本页]

作者
Author: czhang8 时间: 2025-4-8 06:09
标题: 单点能计算时conver=6和conver=8的结果能否比较
各位大佬好！小弟最近为了DIAS分析交了一部分单点能任务，有一部分出现了l502报错。按照sob老师博文（http://sobereva.com/61）里的方法：“G09/16对单点计算默认是SCF=tight，相当于SCF=conver=8，然而此时对密度矩阵的收敛要求有点太严了，往往不容易达到，而且达到的时候能量通常都早已经收敛到极高精度了。对于计算单点能的目的可放心降到conver=6，相当于把收敛标准放宽100倍，通常收敛时能量变化已经非常小了。”加了scf（novaracc,conver=6)后确实收敛了。于是小弟有了另一个顾虑：其他部分的单点能计算用默认设置自动收敛了，也就是 conver=8，其结果能量值能否与用conver=6的结果放在一起进行比较？

如需要统一标准，以下两个方法是否均可行呢？
1）对其他部分同一加上scf（novaracc,conver=6)指令进行计算
2）对默认设定下难以收敛的任务，以此运算结果的波函数为初猜，进一步进行conver=8的运算

第一次发帖，如果有不恰当的地方还请各位海涵，感谢指教！

作者
Author: sobereva 时间: 2025-4-8 06:40
对于普通泛函算单点能，可以比较

作者
Author: zjxitcc 时间: 2025-4-8 07:07
正常收敛限的任务去读取scf(conver=6)任务的轨道，收敛后再比较。如果是开壳层体系，注意检查波函数稳定性。

作者
Author: czhang8 时间: 2025-4-8 07:11
感谢各位！确实是开壳层的计算，用的b3lyp-d3 和def2-tzvp，我去学习一下怎么检查波函数稳定性。另外还想请问一下，加了XQC/QC才达到收敛的任务，和默认算法的任务，能量之间可以比较吗？

作者
Author: zjxitcc 时间: 2025-4-8 12:56

czhang8 发表于 2025-4-8 07:11
感谢各位！确实是开壳层的计算，用的b3lyp-d3 和def2-tzvp，我去学习一下怎么检查波函数稳定性。另外还想请 ...

可以。

作者
Author: ZiyuSong 时间: 2026-4-30 09:32

sobereva 发表于 2025-4-8 06:40
对于普通泛函算单点能，可以比较

老师我在根据mcpb.py对一个从酶促反应里挖出来的含过渡金属的开壳团簇体系（200+原子），含有较高负电（-5）进行一系列gaussian计算，小模型可以收敛（opt,fre计算），但用来计算mk-resp电荷的大模型试了很多方法scf都无法收敛，如果我的目的只是单纯为了得到这个体系里金属的参数跑md，可不可以在计算大模型的mk-resp电荷时使用conver=5？观察了scf log其实只要十来步就能达到这个能级但后面就开始震荡了

作者
Author: sobereva 时间: 2026-4-30 11:18

ZiyuSong 发表于 2026-4-30 09:32
老师我在根据mcpb.py对一个从酶促反应里挖出来的含过渡金属的开壳团簇体系（200+原子），含有较高负电（- ...

至少=6
该怎么解决怎么解决
解决SCF不收敛问题的方法
http://sobereva.com/61

作者
Author: ZiyuSong 时间: 2026-4-30 23:14

sobereva 发表于 2026-4-30 11:18
至少=6
该怎么解决怎么解决
解决SCF不收敛问题的方法

好的，谢谢sob老师

作者
Author: ZiyuSong 时间: 2026-4-30 23:15

sobereva 发表于 2026-4-30 11:18
至少=6
该怎么解决怎么解决
解决SCF不收敛问题的方法

好的，谢谢sob老师

作者
Author: ZiyuSong 时间: 2026-5-13 04:58

sobereva 发表于 2026-4-30 11:18
至少=6
该怎么解决怎么解决
解决SCF不收敛问题的方法

老师我发现这个体系优化时scf在第一步优化前的收敛过程中rmsdp可以收敛得比较顺利但maxdp在缓慢下降的趋势中出现了多次震荡（在e-3-e-4）而且在默认圈数耗尽前也没法收敛，我现在想的解决办法一是增大默认圈数，二是减少macconv提早进入xqc，不知道这样合理吗？或者先设置成conver=6进行优化再用默认值优化一遍？大概哪种解决方式更合理？

作者
Author: sobereva 时间: 2026-5-13 16:22

ZiyuSong 发表于 2026-5-13 04:58
老师我发现这个体系优化时scf在第一步优化前的收敛过程中rmsdp可以收敛得比较顺利但maxdp在缓慢下降的趋 ...

增大默认圈数有用的前提是能看到有收敛趋势

没事甭用QC，博文里明确说了缺点和风险，博文里那么多更值得优先考虑的做法干嘛不用

作者
Author: ZiyuSong 时间: 2026-5-13 22:47

sobereva 发表于 2026-5-13 16:22
增大默认圈数有用的前提是能看到有收敛趋势

没事甭用QC，博文里明确说了缺点和风险，博文里那么多更值 ...

是这样的老师，我确实仔细看了那个博文里的所有方法，适合我这个体系的我也加上去了。现在的的泛函本身没有弥散(6311g**)，所以在scf里就只加了fermi没有使用novaracc,noincfock,int=acc2e=12因为软件是g16).然后因为体系里有过渡金属ni且观察到gap较小,所以提升了vshift=500，隐式溶剂使用了ether是对应一篇文献里的参数。这个体系是用mcpb.py截出来的酶促反应中心含有个卟啉大环，我冻结了除alpha carbon之外的几个外周的o，我附上了我用的routine，希望老师能抽空看一眼对于我这个体系还有什么更好的解决办法，因为目前尝试的所有没有qc的方法确实没法收敛

作者
Author: sobereva 时间: 2026-5-14 01:09

ZiyuSong 发表于 2026-5-13 22:47
是这样的老师，我确实仔细看了那个博文里的所有方法，适合我这个体系的我也加上去了。现在的的泛函本身没 ...

弥散函数是对于基组而言的不是对于泛函而言的

QC容易收敛到不稳定波函数，如果非要用这个，最后记得做波函数稳定性测试

别乱用Guess=Mix，根本不是这种场合用的

也可以尝试PBE0，由于gap更大因此SCF可能比TPSSh更容易收敛

作者
Author: ZiyuSong 时间: 2026-5-15 04:00

sobereva 发表于 2026-5-14 01:09
弥散函数是对于基组而言的不是对于泛函而言的

QC容易收敛到不稳定波函数，如果非要用这个，最后记得做 ...

好的，谢谢老师

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