计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 请教Molpro计算振动频率的问题

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

查看数: 4158 | 评论数: 4 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

chenguangan

发布时间: 2021-12-30 11:14

正文摘要:

各位大佬，我在用molpro搜索三原子体系过渡态并计算振动频率时遇到了一个问题。程序成功搜索出了过渡态的结构，跟我通过计算预估的结构基本一致， &nb ...

回复 Reply

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就不会提示，

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

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

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

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

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

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

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

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

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