计算化学公社

标题: ORCA几何优化输出文件中SCF迭代疑问 [打印本页]

作者
Author: xinyanli 时间: 2022-7-14 14:59
标题: ORCA几何优化输出文件中SCF迭代疑问
我用Multiwfn产生了ORCA的几何优化输入文件，利用ORCA 5.0进行几何优化也正常结束了。但是关于其中输出文件的SCF迭代时，第一步迭代出现如下图所示的警告。这是什么意思呢？而且第一步迭代花了整整十个小时，而后面的每步迭代基本只需要半小时左右，是初始结构不合理吗？这对最终的优化结果有影响吗？

作者
Author: hebrewsnabla 时间: 2022-7-14 15:59
后面的每一步都会读取前一步收敛的轨道作为初猜，所以比第一步快是正常的。

第一个Warning 只是表示开启TRAH，这是默认策略，没什么不正常的，只是说明这个体系SCF比较难收敛。但是考虑到收敛耗时如此之长，并且多次negative hessian warning，也许你需要检查一下波函数稳定性。结构不合理是SCF难收敛的一种可能的原因，但是并不一定是这个原因，光凭你给出的信息并不能判断结构合不合理。

作者
Author: xinyanli 时间: 2022-7-14 19:28

hebrewsnabla 发表于 2022-7-14 15:59
后面的每一步都会读取前一步收敛的轨道作为初猜，所以比第一步快是正常的。

第一个Warning 只是表示开启 ...

好的，谢谢您。那我想问下波函数稳定性该怎么检查呢？

作者
Author: sobereva 时间: 2022-7-14 20:09

xinyanli 发表于 2022-7-14 19:28
好的，谢谢您。那我想问下波函数稳定性该怎么检查呢？

带着以下关键词，以在收敛后做波函数稳定性测试确保稳妥
%scf
STABPerform true
end
详细看手册

作者
Author: xinyanli 时间: 2022-7-15 10:25

sobereva 发表于 2022-7-14 20:09
带着以下关键词，以在收敛后做波函数稳定性测试确保稳妥
%scf
STABPerform true

那请问老师我是应该测试初始结构的波函数稳定性吗？还是最后优化好的结构呢？

作者
Author: zjxitcc 时间: 2022-7-15 10:47
本帖最后由 zjxitcc 于 2022-7-15 10:51 编辑

xinyanli 发表于 2022-7-15 10:25
那请问老师我是应该测试初始结构的波函数稳定性吗？还是最后优化好的结构呢？

如果你已经做过了常规计算，那么此时最快的、最直接的做法是检验优化好的结构的波函数稳定性，如果是稳定的那就罢了；如果你还没开始算，则应一开始就检验波函数稳定性，确保读取了稳定的波函数进行结构优化；优化完成后，仍要检查波函数稳定性。
另外，如果是含d或f单电子的过渡金属分子，可能有多个SCF解，甚至多个SCF稳定解。即使你检验了、发现波函数是稳定的，也不能保证是能量最低的、合理的解。

作者
Author: xinyanli 时间: 2022-7-15 11:09

zjxitcc 发表于 2022-7-15 10:47
如果你已经做过了常规计算，那么此时最快的、最直接的做法是检验优化好的结构的波函数稳定性，如果是稳定 ...

好的，明白了。谢谢老师。

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