计算化学公社

标题: 优化过渡态没有收敛，但是虚频符合要求 [打印本页]

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

作者
Author: gauss98 时间: 2022-6-15 10:04
根据经验，-41的虚频多半不是反应坐标，大概率是某个转动震动模式。

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

作者
Author: sai77 时间: 2022-6-17 08:59
calcall跑了114步...这要是体系大点得跑多少时间

作者
Author: mfdsrax2 时间: 2022-6-19 14:40

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还是增加？或者尝试别的方法？
(, 下载次数 Times of downloads: 12)
(, 下载次数 Times of downloads: 17)

作者
Author: wzkchem5 时间: 2022-6-19 15:25

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

可以加大格点试试

欢迎光临计算化学公社 (http://bbs.keinsci.com/)