ORCA的波函数稳定性分析离纳入隐式溶剂化还有多远？

啊不错的飞过海

依ORCA 5.0.3手册说法，ORCA波函数稳定性分析模块尚未实现溶剂化模型支持[1]；
但同手册里，依赖电子hessian的TRAH SCF[2][3]与过程类似的TDDFT[4]均有至少是C-PCM级别的溶剂化模型支持，TDDFT还支持SMD。
ORCA是否已可以做带隐式溶剂化的波函数稳定性分析，还是仍待开发相关功能？

类似的情况也出现在RI支持上，TRAH支持RIJK密度拟合[2][3]（虽说不如RIJCOSX快），但波函数稳定性分析不支持。

[1]ORCA 5.0.3手册 9.10，p527
[2]ORCA 5.0.3手册 8.1.5，p111
[3]ORCA 5.0.3手册 9.7.7，p513-514
[4]ORCA 5.0.3手册 9.26.6，p798

wzkchem5

我们内部的程序版本已经支持对CPCM、SMD做稳定性分析了，所以ORCA 6.0会支持。估计还有几个月就发布了。
RIJK还不支持稳定性分析，主要是因为用RIJK的人越来越少了，COSX越来越快，导致RIJK比RIJCOSX快的情况越来越少，所以没有投入太多精力在RIJK上面

zjxitcc

就目前来说，可行的workaround有例如：
（1）ORCA算完之后导出轨道到Gaussian，使用Gaussian在CPCM下使用同样的泛函和基组进行波函数稳定性分析，虽然Gaussian的SCF会稍微迭代几圈，但这不会改变波函数本质，往往也不影响波函数稳定性。对应的小程序是ORCA的orca_2mkl，MOKIT的mkl2fch或molden2fch。这里所说的迭代几圈是由于积分格点不一致、RI近似等因素造成的，并非传轨道的小程序有什么近似。
（2）Gaussian在同样的计算机别下进行计算和波函数稳定性分析，传轨道进ORCA做计算，不再检验波函数稳定性。对应的小程序是MOKIT的fch2mkl和ORCA的orca_2mkl。