请教有机分子用wB97M-V几何优化中遇到的问题

啃骨魔

各位专家好！
我想算一下四氰基乙烯的HOMO用于求亲核指数，由于目标分子用RI-wB97M-V def2-TZVP算的单点，想直接在这个级别优化加单点得了；结果失败：

WARNING: Gradients for DFT-NL dispersion correction are not available,
         therefore, geometry optimizations are not possible.
  ===>   Skipping actual calculation

请问这个“Gradients for DFT-NL dispersion correction are not available”是怎么回事啊？如何解决啊？

最后是B3LYPdef2-TZVP优化后再做的RI-wB97M-V def2-TZVP单点，EHOMO为-0.430521 a.u.

thanhtam

这个泛函没有解析梯度，不能做优化

sobereva

详谈Multiwfn产生ORCA量子化学程序的输入文件的功能
http://sobereva.com/490（http://bbs.keinsci.com/thread-13560-1-1.html）

这里明确说了

• RI-wB97M-V/def2-TZVP。关键词：wB97M-V def2-TZVP def2/J RIJCOSX strongSCF grid5 gridx5
根据Phys. Chem. Chem. Phys., 19, 32184 (2017)、Mol. Phys., 115, 2315 (2017)、J. Chem. Theory Comput., 15, 3610 (2019)等测试，ωB97M-V在除了双杂化泛函以外的泛函中，无论算主族还是含过渡金属的体系，无论算反应能、势垒还是弱相互作用，性能都是顶级的，在《简谈量子化学计算中DFT泛函的选择》我也提到了。这个泛函是meta-杂化泛函并且结合了VV10色散校正，也正是由于这种色散校正形式，导致ORCA不支持其解析梯度，因此很不适合用于几何优化、振动分析。众所周知，几何优化用的级别可以比单点要低，这点在《浅谈为什么优化和振动分析不需要用大基组》（http://sobereva.com/387）专门提了，因此推荐大家用相对便宜的B97-3c或RI-B3LYP-D3(BJ)/def2-TZVP(-f)来优化结构，最后用RI-wB97M-V/def2-TZVP来算单点。相对于BLYP、B3LYP这种较为粗放的泛函，ωB97M-V属于稍昂贵一些但精度明显更高的泛函，可以达到定量上较准确的结果，所以有理由用比起默认更好的DFT积分格点。而且根据我的经验此泛函对积分格点要求稍高，用太次的话有时候连SCF收敛都费劲。另外，由于这个泛函只拿来用于算单点，因此积分格点质量用得好一些也不会造成耗时过高。考虑到这些，用了grid5将DFT交换相关泛函的积分格点为比起默认更好的档次（与Gaussian的int=ultrafine接近），而gridx5代表将COSX部分的格点也设为比默认明显更好的档次。对于ωB97M-V这样的非双杂化泛函而言，对于一般问题用def2-TZVP就已经够用了，要求更高的话也可以考虑def2-QZVPP。如果你要算弱相互作用能的话也推荐用def2-QZVPP，若结合counterpoise校正还能更好一点。由于ORCA的默认的SCF收敛限相对于wB97M-V的精度水平来说过于宽松了些，因此同时用了strongSCF来让能量收敛到更高精度。

wuzhiyi

如果真的很小的话 可以尝试NUMGRAD

啃骨魔

谢谢各位老师指正！

rendong

这个泛函好像也不能用于TDDFT计算