计算化学公社

标题: 过渡态优化，用calcall结合小步长仍不收敛，求助 [打印本页]

作者
Author: h840473807 时间: 2017-12-25 16:05
标题: 过渡态优化，用calcall结合小步长仍不收敛，求助
用B3LYP/def2-SVP等级结合TS方法优化一个过渡态，calcfc下能量振荡，无法收敛,接着用calcall,仍然振荡无法收敛，然后结合小步长，maxstep=5,直到maxstep=1(均带有notrust)，均无法收敛，出现振荡，取出calcall,maxstep=1优化中振荡的两个结构做freq分析，发现有均有两个虚频，主虚频和我的反应坐标对应，还有一个小虚频无法消除。请问我接下来可以用什么方法来消除虚频，优化收敛？ (, 下载次数 Times of downloads: 53) (, 下载次数 Times of downloads: 43)

作者
Author: h840473807 时间: 2017-12-25 16:09
附件中bts开头为优化文件，freq开头为振动分析文件。freqbts6e和freqbts6e7分别为图中优化曲线（对应bts6e）的第6，第7个结构

作者
Author: sobereva 时间: 2017-12-25 17:09
尝试gdiis（其实你这体系应该加上D3校正，配体间的相互作用还是比较明显的）

作者
Author: ikea1984 时间: 2017-12-25 17:37
有反应物/产物坐标没，感觉QST2比较容易收敛，CalcAll+小步长仍然很难

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

作者
Author: h840473807 时间: 2017-12-25 20:44
@ikea1984 多谢解答

作者
Author: sobereva 时间: 2017-12-25 22:19

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

对计算量没任何影响

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

作者
Author: sobereva 时间: 2018-1-15 00:05

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

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

作者
Author: h840473807 时间: 2018-1-15 12:29

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 终于消除了虚频

作者
Author: yanmu 时间: 2018-1-15 15:30

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

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

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

作者
Author: LittlePupil 时间: 2022-12-20 21:52

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

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

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

作者
Author: h840473807 时间: 2023-3-6 16:29

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

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

作者
Author: lurensan 时间: 2023-10-12 09:34

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

你好，请问用gdiis时是不能用maxstep缩短步长吗？

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