选取 RI-wB97M-V/def2-TZVP在ORCA中file orca_main/mainchk.cpp, line 5178的问题

Meso_Alt

本人在照猫画虎地尝试使用ORCA与Multiwfn来进行弱相互作用力（氢键）与光谱分析（目前并未加入水分子来形成氢键，仅添加了水溶剂模型），选用了 RI-wB97M-V/def2-TZVP计算级别时出现了如下问题：

   在使用Multiwfn生成ORCA输入文件后，将生成的1.inp输入进入ORCA。

    输入：orca 1.inp>1.out
    输出：[file orca_main/mainchk.cpp, line 5178]: Error (ORCA_MAIN): ... aborting the run




    Multiwfn操作：输入文件为5HT的.pdb分子模型文件；开启了SMD溶剂模型，选择了水作为溶剂；任务选择为优化+振动分析；计算级别选择RI-wB97M-V/def2-TZVP。

    使用ORCA 5.0.4 ,  Multiwfn 3.8(dev)

随帖附上.pdb模型、生成的.inp文件、输出的.out文件。

望诸位大佬指出我的操作哪里存在错误与疏忽，感谢。
本人初步猜测是该计算级别并不支持振动分析，因为删去1.inp中的振动分析任务代码“freq”后其可以单独运行优化任务，但是如果我只想做IR光谱，我应该如何选择方法、基组与任务？

附问：多个分子可以放进同一个模型文件中吗，譬如我想测算某些分子特定位置的水形成的氢键对IR光谱的影响诸如此类的？在相关论文中见到过类似操作，但不懂如何实现。

Meso_Alt

补充：RI-B3LYP-D3(BJ)/def2-TZVP(-f)和RI-B3LYP-D3(BJ)/def2-TZVP两种计算级别在同样的任务上会出现同样的问题。
我在怀疑TZVP基组不能进行振动分析？我水平过于有限并不能分析这是什么问题，恳请诸位大佬指教。

Meso_Alt

在.out文件里找到了报错，我先去查手册

Meso_Alt

在文章http://sobereva.com/490中找到了“RI-wB97M-V/def2-TZVP……由于这种色散校正形式，导致ORCA不支持其解析梯度，因此很不适合用于几何优化、振动分析。”原来是选错计算级别了。

sobereva

Meso_Alt 发表于 2024-7-9 01:28
在文章http://sobereva.com/490中找到了“RI-wB97M-V/def2-TZVP……由于这种色散校正形式，导致ORCA不支持 ...

这个老说法我早就改了，当前这个博文里没有这句话
ORCA 5开始已经有了wB97M-V的解析梯度，但没解析Hessian。由于几何优化之后一般都要做振动分析，所以还是不适合几何优化。

sobereva

Meso_Alt 发表于 2024-7-9 00:43
补充：RI-B3LYP-D3(BJ)/def2-TZVP(-f)和RI-B3LYP-D3(BJ)/def2-TZVP两种计算级别在同样的任务上会出现同样的 ...

只能说某基组是否适合做振动分析（准确度如何、是否划算），没有什么基组是不能做振动分析的

Meso_Alt

sobereva 发表于 2024-7-9 13:50
只能说某基组是否适合做振动分析（准确度如何、是否划算），没有什么基组是不能做振动分析的

理解了，感谢sob老师