计算化学公社

标题: 单点能计算时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，我去学习一下怎么检查波函数稳定性。另外还想请 ...

可以。

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