计算化学公社

标题: 求助，水解反应过渡态结构优化报错 [打印本页]

作者
Author: ReinForceKazuha 时间: 2025-6-5 18:05
标题: 求助，水解反应过渡态结构优化报错
本帖最后由 ReinForceKazuha 于 2025-6-5 18:13 编辑

各位老师大家好，我是Gaussian的入门者，现在想复现该文献的水解反应机理，之前按照sob老师的指导设置了输入文件，但是运行了三个小时后输出文件报错，提示结构优化部属不够导致不收敛（能量没有震荡），除此之外，我在output文件中发现“ Linear search not attempted -- option 19 set.”语句多次出现，计算耗时比较长。我把gjf输入文件作为附件上传了（output文件提示无法上传这么大的附件，截图了部分），想麻烦各位老师批评指正我在计算过程中的不当的地方，我该怎么做才能得到正确的过渡态呢，谢谢！
[attach]114124[/attach]

(, 下载次数 Times of downloads: 14) (, 下载次数 Times of downloads: 14) (, 下载次数 Times of downloads: 15)

作者
Author: chands 时间: 2025-6-5 18:59
本帖最后由 chands 于 2025-6-5 19:02 编辑

# opt=(calcfc,ts,noeigentest) freq 6-311++g(d,p) scrf=(smd,solvent=water)
geom=connectivity rm062x

复制代码

你把6-311++g(d,p)换成6-311+g(d,p)试试。社长说了，通常不用给H上弥散（氢负离子除外），可以试试IEFPCM（也就是Gaussian默认的溶剂模型），

作者
Author: ReinForceKazuha 时间: 2025-6-5 21:14

chands 发表于 2025-6-5 18:59
你把6-311++g(d,p)换成6-311+g(d,p)试试。社长说了，通常不用给H上弥散（氢负离子除外），可以试试IEFP ...

好的非常感谢老师！我再试试！

作者
Author: ReinForceKazuha 时间: 2025-6-5 21:19

chands 发表于 2025-6-5 18:59
你把6-311++g(d,p)换成6-311+g(d,p)试试。社长说了，通常不用给H上弥散（氢负离子除外），可以试试IEFP ...

老师您好，我按照您的方式进行尝试，7分钟后就报错了，提示这个
Optimization stopped.
-- Number of steps exceeded, NStep= 40
-- Flag reset to prevent archiving.
请问是需要添加关键词opt=maxcycle=100进行最大步数增加，还是我的初始构型就存在大问题呢？我是按照文献来摆的初始构型
谢谢老师！

作者
Author: chands 时间: 2025-6-5 21:30

ReinForceKazuha 发表于 2025-6-5 21:19
老师您好，我按照您的方式进行尝试，7分钟后就报错了，提示这个
Optimization stopped.
-- Number ...

试试社长的方法http://sobereva.com/164

作者
Author: sobereva 时间: 2025-6-6 00:47
如置顶的新社员必读贴和论坛首页的公告栏所示，试图上传超过500KB的文本型文件上传前必须先压缩再上传，以节约论坛空间、节约下载耗时。

关键词极其糟糕，M06-2X结合SMD做opt freq是经典的自找难收敛、出虚频的组合，我已经在论坛里说过很多次

初始构型也太烂，O8和H6哪能离得那么近，整体结构和文献里的明显有极大的差异

作者
Author: ReinForceKazuha 时间: 2025-6-6 17:09
本帖最后由 ReinForceKazuha 于 2025-6-6 17:12 编辑

sobereva 发表于 2025-6-6 00:47
如置顶的新社员必读贴和论坛首页的公告栏所示，试图上传超过500KB的文本型文件上传前必须先压缩再上传，以 ...

非常感谢sob老师的批评指正，那我按照文献所给出的泛函和基组进行计算看能否得到结果，同时对初始构型再次进行调整。
因为文献展示的是优化后的过渡态，opt=ts是需要自己摆初始构型，所以我在这里犯了一点错误，感谢sob老师，我再试试！

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