请教Molpro计算振动频率的问题

chenguangan · 发表于 Post on 2021-12-30 11:14:19

各位大佬，我在用molpro搜索三原子体系过渡态并计算振动频率时遇到了一个问题。

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

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

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

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

RAL · 发表于 Post on 2021-12-30 11:27:46

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

fatpig · 发表于 Post on 2021-12-30 11:34:58

估计你并没有找到真正的过渡态。

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

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

chenguangan · 发表于 Post on 2021-12-30 15:03:55

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没加弥散函数的，拟合出来的平衡结构和振频也跟其他人差不多，不知道为啥搜个过渡态那么多问题

chenguangan · 发表于 Post on 2021-12-30 15:23:16

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

多谢大佬！

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Molpro] 请教Molpro计算振动频率的问题

浏览过的版块