优化过渡态没有收敛，但是虚频符合要求

mfdsrax2 · 发表于 Post on 2022-6-15 09:18:57

使用b3lyp-d3(bj)/6-31g** calcall优化过渡态，最终轻微震荡，也没有达到收敛限，但是由于是calcall，最后默认频率计算，我看了一下，有-41的虚频，并且符合反应方向，IRC正在算，感觉应该问题不大。那么这个结构能用吗？或者有哪些方法帮助满足收敛限？

gauss98 · 发表于 Post on 2022-6-15 10:04:04

根据经验，-41的虚频多半不是反应坐标，大概率是某个转动震动模式。

wzkchem5 · 发表于 Post on 2022-6-15 20:48:26

一个结构是过渡态的充要条件是梯度为0且虚频有且仅有一个。因为一般是前者容易满足而后者不容易满足，所以大家经常只强调虚频有几个的问题，但是你现在这个情况是后者满足但前者不满足，所以结构也是不能用的，最简单的解决方法是以最后一帧为初始结构，继续带着calcall优化，并且把maxstep设小一些。

sai77 · 发表于 Post on 2022-6-17 08:59:26

calcall跑了114步...这要是体系大点得跑多少时间

mfdsrax2 · 发表于 Post on 2022-6-19 14:40:22

wzkchem5 发表于 2022-6-15 20:48
一个结构是过渡态的充要条件是梯度为0且虚频有且仅有一个。因为一般是前者容易满足而后者不容易满足，所以 ...

经过再次计算，均方根受力已经收敛，最大受力非常接近收敛（当前值0.000453，收敛域0.000450）目前的关键词是#p opt=(calcall,ts,maxstep=2,gdiis,notrust,noeigen) b3lyp/6-31g(d,p) em=gd3bj
但是还是小幅度震荡，还应该怎么样继续调整，减少maxstep还是增加？或者尝试别的方法？

wzkchem5 · 发表于 Post on 2022-6-19 15:25:16

mfdsrax2 发表于 2022-6-19 07:40
经过再次计算，均方根受力已经收敛，最大受力非常接近收敛（当前值0.000453，收敛域0.000450）目前的关键 ...

可以加大格点试试

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

[Gaussian/gview] 优化过渡态没有收敛，但是虚频符合要求