计算化学公社

标题: Mopac如何消虚频 [打印本页]

作者
Author: Shannon 时间: 2015-1-26 07:01
标题: Mopac如何消虚频
请教一个关于用mopac2012优化产生虚频的问题。
用PM7优化了一个35个原子的体系，但是用force关键词算频率的时候总是会产生虚频。
优化时候用的是这个方法：PM7 eps=78.54 recalc=0 FLEPO GNORM=0.0 DMAX=0.05 CHARGE=0
计算频率用的是：PM7 eps=78.54 Force LET CHARGE=0 试过很多关键词，还是不知道该如何消掉虚频……剩余两个虚频，一个 -200，一个 -100。
这个虚频是因为mopac没有解析二阶导数的原因吗？怎么都消不掉。
想用Mopac算出来的带有数目不对的虚频的震动频率，修改删除掉小虚频，然后用过渡态理论TST 计算反应速率，不知道结果可信吗？

作者
Author: sobereva 时间: 2015-1-26 08:20
优化时写上PRECISE关键词增加计算精度

作者
Author: Shannon 时间: 2015-1-26 08:30

sobereva 发表于 2015-1-26 08:20
优化时写上PRECISE关键词增加计算精度

precise用过了，也还是存留有虚频。我怀疑是因为体系内有几个受阻转动，可能是它们导致的顽固虚频。不知道能否在输出结果里面把negative eigenvalue的值直接改成正值，然后和 ORCA还有高斯结果拼到一起，塞到kisthelp里面算反应速度（感觉这样有点篡改数据，通不过审稿人的样子）。

作者
Author: sobereva 时间: 2015-1-26 09:22
可以把输入输出文件发给MOPAC作者J.J.P.Stewart问他怎么处理，他回信很快。
才35个原子，没必要用半经验

作者
Author: Shannon 时间: 2015-1-26 22:06

sobereva 发表于 2015-1-26 09:22
可以把输入输出文件发给MOPAC作者J.J.P.Stewart问他怎么处理，他回信很快。
才35个原子，没必要用半经验

谢谢sob，我多收集点数据，然后问问Stewart该怎么办

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