过渡态优化，用calcall结合小步长仍不收敛，求助

h840473807 · 发表于 Post on 2017-12-25 16:05:59

用B3LYP/def2-SVP等级结合TS方法优化一个过渡态，calcfc下能量振荡，无法收敛,接着用calcall,仍然振荡无法收敛，然后结合小步长，maxstep=5,直到maxstep=1(均带有notrust)，均无法收敛，出现振荡，取出calcall,maxstep=1优化中振荡的两个结构做freq分析，发现有均有两个虚频，主虚频和我的反应坐标对应，还有一个小虚频无法消除。请问我接下来可以用什么方法来消除虚频，优化收敛？

h840473807 · 发表于 Post on 2017-12-25 16:09:44

附件中bts开头为优化文件，freq开头为振动分析文件。freqbts6e和freqbts6e7分别为图中优化曲线（对应bts6e）的第6，第7个结构

sobereva · 发表于 Post on 2017-12-25 17:09:45

尝试gdiis（其实你这体系应该加上D3校正，配体间的相互作用还是比较明显的）

ikea1984 · 发表于 Post on 2017-12-25 17:37:07

有反应物/产物坐标没，感觉QST2比较容易收敛，CalcAll+小步长仍然很难

h840473807 · 发表于 Post on 2017-12-25 20:44:03

@sobereva 多谢老师，用了gdiis优化，8步就收敛了，请问加入D3校正，会增加很多计算量吗

h840473807 · 发表于 Post on 2017-12-25 20:44:43

@ikea1984 多谢解答

sobereva · 发表于 Post on 2017-12-25 22:19:24

h840473807 发表于 2017-12-25 20:44
@sobereva 多谢老师，用了gdiis优化，8步就收敛了，请问加入D3校正，会增加很多计算量吗

对计算量没任何影响

yanmu · 发表于 Post on 2018-1-14 21:32:48

我想问一下你是用的# opt=(calcall,notrust,ts，Maxstep=1)找到的过渡态对应两个虚频，后来直接在此基础上加了gdiis消除虚频了吗？我最近遇到了类似的问题，想请教一下，谢谢你

sobereva · 发表于 Post on 2018-1-15 00:05:28

yanmu 发表于 2018-1-14 21:32
我想问一下你是用的# opt=(calcall,notrust,ts，Maxstep=1)找到的过渡态对应两个虚频，后来直接在此基础上 ...

GDIIS没有直接的“消虚频”的功效，只不过有助于势能面平缓的情况下的几何优化收敛
应当怎么合理处理看
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278

h840473807 · 发表于 Post on 2018-1-15 12:29:00

yanmu 发表于 2018-1-14 21:32
我想问一下你是用的# opt=(calcall,notrust,ts，Maxstep=1)找到的过渡态对应两个虚频，后来直接在此基础上 ...

我用opt=(calcall,notrust,ts，Maxstep=1)不收敛，改用opt=(calcfc,notrust,ts，gdiis)收敛但仍有虚频率，最后用opt=(calcfc,notrust,ts，gdiis，tight) int=ultrafine 终于消除了虚频

yanmu · 发表于 Post on 2018-1-15 15:30:14

h840473807 发表于 2018-1-15 12:29
我用opt=(calcall,notrust,ts，Maxstep=1)不收敛，改用opt=(calcfc,notrust,ts，gdiis)收敛但仍有虚频率 ...

好的，谢谢你，我试一下看看可不可以

avatar_man · 发表于 Post on 2022-12-20 09:22:11

我的也是加上gdiis之后，虚频就消失了，请问其他的结构用不用全采用 gdiis重新优化？已经算了七八个中间态了，全改gdiis工作量有点大啊！

LittlePupil · 发表于 Post on 2022-12-20 21:52:15

avatar_man 发表于 2022-12-20 09:22
我的也是加上gdiis之后，虚频就消失了，请问其他的结构用不用全采用 gdiis重新优化？已经算了七八个中间 ...

不需要，只有当前理论方法对应的势能面和之前不同时才需要重算，而gdiis是一种优化算法，只改变体系在势能面上的运动方式，不会改变势能面。

红烧肉盖饭 · 发表于 Post on 2023-2-23 19:02:29

您好，我想问一下如何得到优化曲线，我把未收敛的log文件下载导入gaussview，在结论里看不到优化曲线

h840473807 · 发表于 Post on 2023-3-6 16:29:53

红烧肉盖饭发表于 2023-2-23 19:02
您好，我想问一下如何得到优化曲线，我把未收敛的log文件下载导入gaussview，在结论里看不到优化曲线

先打开GV用open打开.log文件，同时勾选窗口左下角的Read Intermediate Geometries即可

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

[Gaussian/gview] 过渡态优化，用calcall结合小步长仍不收敛，求助