计算化学公社

标题: 请教Molpro计算振动频率的问题 [打印本页]

作者
Author: chenguangan 时间: 2021-12-30 11:14
标题: 请教Molpro计算振动频率的问题
各位大佬，我在用molpro搜索三原子体系过渡态并计算振动频率时遇到了一个问题。

         程序成功搜索出了过渡态的结构，跟我通过计算预估的结构基本一致，

         但是在紧接着计算振动频率的时候，计算出来是两个虚频一个实频，并且out文件里警告梯度不等于零。

         我就没弄明白，都已经搜索出来这个静态点了，为什么计算振频的时候梯度会不等于零呢？

         有没有懂的大佬，请教教我！谢谢！

作者
Author: RAL 时间: 2021-12-30 11:27
原因在下文中提到了
《Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法》
http://bbs.keinsci.com/thread-633-1-1.html

作者
Author: fatpig 时间: 2021-12-30 11:34
估计你并没有找到真正的过渡态。

虽然你写了nosym，但又把两个P-H键设成一样，这其实就把过渡态的结构限制在了C2v，只是波函数使用C1。

所以可以再加一个变量，令两个键长不相等，再试试。

作者
Author: chenguangan 时间: 2021-12-30 15:03

fatpig 发表于 2021-12-30 11:34
估计你并没有找到真正的过渡态。

虽然你写了nosym，但又把两个P-H键设成一样，这其实就把过渡态的结构限 ...

大佬你好，

因为我这个分子是PH2，这个过渡态肯定是个C2V结构的。

我担心搜出来是个不对称的咋办，我按你说的先搜索试试吧。

我刚刚把H和P的基组都改成aug-cc-pV(Q+d)Z算了一下，频率就对了。。。

我看sobereva大神的帖子说H不用加弥散函数，而且H肯定用不着加d吧。。。

我整个势能面拟合用了13000+势能点，全是H没加弥散函数的，拟合出来的平衡结构和振频也跟其他人差不多，不知道为啥搜个过渡态那么多问题

作者
Author: chenguangan 时间: 2021-12-30 15:23

RAL 发表于 2021-12-30 11:27
原因在下文中提到了
《Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法》
http://bbs.keins ...

多谢大佬！

好像是这样的，我刚看了一下，结构优化梯度小于10-3的时候就收敛了，计算频率的时候大于10-4就会提示梯度不为零，小于10-4就不会提示，

其实梯度已经很小了，结果还能对上

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