计算化学公社

标题: 希望了解Gaussian当中stable任务和opt+freq任务最后输出的单点能有差异的原因 [打印本页]

作者
Author: philartist 时间: 2026-2-3 14:54
标题: 希望了解Gaussian当中stable任务和opt+freq任务最后输出的单点能有差异的原因
本帖最后由 philartist 于 2026-2-3 15:06 编辑

B45.gjf这个文件包含多个步骤，分别是opt+freq和stable检验波函数稳定性，任务结束后发现没有虚频，但是波函数存在RHF-->UHF的不稳定性，所以通过B45A1.gjf任务试图解决这个问题，奇怪的地方出现在这个任务的输出文件当中。

B45A1.gjf这个文件包含三个步骤：1.读取B45.chk中的结构和波函数，并使用stable=opt稳定波函数，新的chk文件保存为B45A1.chk。2.读取B45A1.chk当中的结构和波函数，做几何优化和振动分析，新的波函数覆盖B45A1.chk。3.读取覆盖后的B45A1.chk当中的结构和波函数，并用stable关键词检验最终的波函数是否稳定。

然而奇怪的事情是，第1步将波函数稳定下来之后单点能为'HF=-2397.8858973'，第2步为opt+freq，在Gaussian当中分为两个子步骤，单点能均为'HF=-2397.8870953'，与第1步不同，第3步是检验波函数稳定性，最后的单点能又变回'HF=-2397.8858973'.

这种情况出现的不太多，其他结构大部分时候stable和opt+freq任务给出的单点能还是一样的，少部分时候会遇到类似的情况。

尽管这两个结果差异很小，大约相差0.03 eV，但我还是想要了解一下通过stable和stable=opt这两个关键词计算后的结果给出的单点能为啥和通过opt和freq这两个关键词计算后的结果给出的单点能不一样？原则上哪个更靠谱一点？

B45A1的输入输出文件在附件当中。

作者
Author: wal 时间: 2026-2-3 15:32
opt的link没写%chk，事实上第三个stable任务读取的还是第一个任务stable过后的波函数

作者
Author: philartist 时间: 2026-2-3 17:39

wal 发表于 2026-2-3 15:32
opt的link没写%chk，事实上第三个stable任务读取的还是第一个任务stable过后的波函数

OMG，确实漏了一个%chk，大意了，多谢多谢！

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