过渡态搜索结果有两个虚频，如何消掉其中一个？

yaoyuan0711 · 发表于 Post on 2022-2-12 16:06:34

各位老师好，我用Gaussian09搜索如下反应（N上的H转移到O上，C与N成键）的过渡态时，结果出现了两个虚频：第一个-1593是对应于反应路径的，第二个-962对应于甲基上3个氢原子的旋转。根据sob老师的帖子http://sobereva.com/278，我先加上了int=superfine还是两个虚频，然后取最后结构朝着-962的虚频分别旋转了0.3，0.4，结果跑到反应物和产物去了。这样的情况该怎么做？
另外，这个结构如果不加隐式溶剂模型的话，很容易找到只有一个虚频的过渡态。
搜索过渡态的关键词：#p opt=(calcfc,ts,noeigen) freq b3lyp/6-311g(d,p) scrf=(solvent=water,pcm) em=gd3bj int=superfine

计算结果log文件.zip (683.21 KB, 下载次数 Times of downloads: 9)

yaoyuan0711 · 发表于 Post on 2022-2-14 11:25:02

喵星大佬发表于 2022-2-14 11:17
不用特意纠结那个玩意，用了，可能有改进，不用，影响也不大，定性还是定量基本上都轮不到那点热力学校正 ...

好的，谢谢您的指导

喵星大佬 · 发表于 Post on 2022-2-14 11:17:32

本帖最后由喵星大佬于 2022-2-14 11:19 编辑

yaoyuan0711 发表于 2022-2-14 11:12
谢谢！但还有一个问题，b3lyp/def2tzvp能找到频率校正因子，def2svp好像没有。。。

不用特意纠结那个玩意，用了，可能有改进，不用，影响也不大，定性还是定量基本上都轮不到那点热力学校正量上的微小影响，其他能产生大的多的误差的地方注意好才是真的

yaoyuan0711 · 发表于 Post on 2022-2-14 11:12:07

fatpig 发表于 2022-2-14 10:32
理论上就是这样。

我们组以前都是用pople+小赝势，自从换了def2-svp/def2-tzvp，嘿，真香...

谢谢！但还有一个问题，b3lyp/def2tzvp能找到频率校正因子，def2svp好像没有。。。

喵星大佬 · 发表于 Post on 2022-2-14 10:50:09

fatpig 发表于 2022-2-14 10:32
理论上就是这样。

我们组以前都是用pople+小赝势，自从换了def2-svp/def2-tzvp，嘿，真香...

瓜丝用户，一般目的(反应能势垒之类的，包括激发态)，DFT，基本上6-31G*(*)，def-TZVP，def2-TZVP三个级别就差不多了，对于含阴离子的最多也是最多到may/jul-cc-pVTZ，基本用不到其他的基组

fatpig · 发表于 Post on 2022-2-14 10:32:21

yaoyuan0711 发表于 2022-2-14 08:49
一换基组就成了，不过我得把所有计算都重新算了。。。
是pople基组没有def2系列的基组构建的好吗？so ...

理论上就是这样。

我们组以前都是用pople+小赝势，自从换了def2-svp/def2-tzvp，嘿，真香...

yaoyuan0711 · 发表于 Post on 2022-2-14 10:09:23

喵星大佬发表于 2022-2-14 09:50
我的意思是，如果体系所有都跟你现在这个东西这么大，用def2-tzvp去算去振动分析，也不是不能接受，基本 ...

好的，谢谢您的建议

喵星大佬 · 发表于 Post on 2022-2-14 09:50:07

本帖最后由喵星大佬于 2022-2-14 09:51 编辑

yaoyuan0711 发表于 2022-2-14 08:57
您意思是说对于我当前的体系，有一个额外的-962的虚频也能接受？因为原子数不多，想计算的更精确点，用了 ...

我的意思是，如果体系所有都跟你现在这个东西这么大，用def2-tzvp去算去振动分析，也不是不能接受，基本上这个大小的基组做到7，80个原子的振动分析/找过渡态之类的就会很吃力了，可以用def-tzvp，在瓜丝里写作tzvp，这个基组跟6-311G**基本差不多，稍微大一丢丢

yaoyuan0711 · 发表于 Post on 2022-2-14 09:47:47

sai77 发表于 2022-2-14 09:35
内存设定是%mem=xxxGB

好的，谢谢您了，以前没注意这个问题

sai77 · 发表于 Post on 2022-2-14 09:35:02

yaoyuan0711 发表于 2022-2-14 08:54
感谢您的建议！我在输入文件一直写了%nprocshared=24，但是还是出现了您说的这种情况，这是什么原因？

内存设定是%mem=xxxGB

yaoyuan0711 · 发表于 Post on 2022-2-14 08:57:50

喵星大佬发表于 2022-2-14 08:53
这个基组相当大，尤其是振动分析会非常费劲，不过就这个结构来看还好体系不算大也不是完全不能接受就是了

您意思是说对于我当前的体系，有一个额外的-962的虚频也能接受？因为原子数不多，想计算的更精确点，用了较大的基组，后续用PWPB95/def2qzvpp计算单点、

yaoyuan0711 · 发表于 Post on 2022-2-14 08:54:40

ahxb 发表于 2022-2-12 20:59
插个和你的问题无关的建议，你的输入文件里似乎没有写%nprocshared要内存，程序只用了32MB内存，但你又用了 ...

感谢您的建议！我在输入文件一直写了%nprocshared=24，但是还是出现了您说的这种情况，这是什么原因？

喵星大佬 · 发表于 Post on 2022-2-14 08:53:39

这个基组相当大，尤其是振动分析会非常费劲，不过就这个结构来看还好体系不算大也不是完全不能接受就是了

yaoyuan0711 · 发表于 Post on 2022-2-14 08:49:46

本帖最后由 yaoyuan0711 于 2022-2-14 09:00 编辑

fatpig 发表于 2022-2-12 16:34
基组换成 def2-tzvp 试试...

一换基组就成了，不过我得把所有计算都重新算了。。。
是pople基组没有def2系列的基组构建的好吗？sob老师在http://sobereva.com/336提到过这一点

yaoyuan0711 · 发表于 Post on 2022-2-14 08:48:37

北大-陶豫发表于 2022-2-12 16:25
甲基转60度试试

谢谢！我按照您说的分别旋转+60度和-60度，-60度搜索出来虚频减小到190了，我应该继续旋转？

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

[新手求助] 过渡态搜索结果有两个虚频，如何消掉其中一个？

评分 Rate

评分 Rate

评分 Rate

评分 Rate