Scan后做频率分析没有虚频

gnli · 发表于 Post on 2022-4-15 22:55:40

本帖最后由 gnli 于 2022-4-15 23:00 编辑

大家好我想找C=O加氢反应的过渡态（金属中心的氢原子加到C上，N上的氢加到O上）
1.先扫描了C-H键的距离，得到了如下扫描曲线。
2.然后取能量最高点做频率计算，却没有虚频。
3.然后又在峰附近用很小的步长继续做了一个scan 但是峰值附近的结构还是没有虚频。做了Mn和Fe两个金属的计算，其他计算也是同样方法，频率计算都正常。
比较了这个结构和其他几个的结构，键长非常接近。所以很纳闷。
请问这种情况原因是什么，求大神指点。

zjxitcc · 发表于 Post on 2022-4-15 23:25:55

（1）先确认是否波函数是否稳定、过渡金属上的自旋布居是否符合期望，以免收敛到了错误的波函数解上；
（2）柔性扫描最高点通常是合理的过渡态初始结构，做freq有个明显虚频，但这是经验规则，不是物理定律，少数体系违背经验规则是可以的。

gnli · 发表于 Post on 2022-4-16 00:03:29

zjxitcc 发表于 2022-4-15 23:25
（1）先确认是否波函数是否稳定、过渡金属上的自旋布居是否符合期望，以免收敛到了错误的波函数解上；
（2 ...

做了scan峰值处结构的波函数稳定性测试：
#p stable 6-311+g(d,p) empiricaldispersion=gd3bj pbe1pbe
结果显示：“The wavefunction is stable under the perturbations considered.”

那请问出现这种非经验情况，应该怎么找到过渡态呢？

dnlx · 发表于 Post on 2022-4-16 11:09:31

这个协同反应的反应物和过渡态“没有局部带显著负电荷”（http://sobereva.com/119），是否需要去掉弥散函数试一下？

gnli · 发表于 Post on 2022-4-16 17:20:46

dnlx 发表于 2022-4-16 11:09
这个协同反应的反应物和过渡态“没有局部带显著负电荷”（http://sobereva.com/119），是否需要去掉弥散函 ...

不加弥散也不行，还是没有虚频

dnlx · 发表于 Post on 2022-4-16 19:13:14

gnli 发表于 2022-4-16 17:20
不加弥散也不行，还是没有虚频

如果是我的话，只能改用ORCA的NEB-TS去寻找过渡态。

gnli · 发表于 Post on 2022-4-17 02:41:08

本帖最后由 gnli 于 2022-4-17 02:43 编辑

ORCA的话还得现学个软件。。。用 opt=(ts,calcfc,noeigen)的话总是跑回反应物了
我现在把scan导出来的结构用最慢的opt=(ts,calcall,noeigen) 继续优化看看情况

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

[Gaussian/gview] Scan后做频率分析没有虚频