oniomQM部分如何限制性优化以及消虚频

liurui

想在大富勒烯笼中放一个CH4固定一个CH来进行限制优化，其中富勒烯用的UFF力场，CH4用的wb97xd方法，发现报错了，尝试了几种帖子(http://sobereva.com/404)中的方法，仍然报如下图错误；此外，发现扫描和频率分开算不报错，但是存在虚频，输入命令为：#p opt=(recalc=4,calcall,verytight,maxcycle=10000) freq oniom(wb97xd/aug-cc-pvdz :uff) nosymm scf=(maxcycle=10000,NoVarAcc)  int=superfine,recalc=1,2,3,4,5都测试了一下，发现recalc=4的虚频要小一些，虚频为-32，想请问如何除掉这种小虚频？或者说如何实现oniom的限制性优化加频率计算呢？

sobereva

强烈不建议用ONIOM算这个，瞎折腾，ONIOM(QM:UFF)对相互作用的描述太烂了，算着玩玩也罢，文章肯定发不出去（用这个甚至还不如用GFN-xTB或PM6-D3H4）。以现在主流计算条件用DFT优化这种体系轻轻松松。B3LYP-D3(BJ)/6-311G**做这种体系优化足以得到不错结果，也花不了多少时间