计算化学公社

标题: Gaussian几何优化达到四个YES后，为什么又多优化一步 [打印本页]

作者
Author: Xulx2022 时间: 2025-1-9 13:56
标题: Gaussian几何优化达到四个YES后，为什么又多优化一步
各位老师好，我在对有机分子和Au6配合物进行几何优化，程序正常停止后，输出文件一共优化了12步。从文件里我看到第11步已经出现了四个yes，但是程序没有停止，会再优化1步，第12步是两个yes，两个no。
输入卡关键词如下：
#p opt=(calcall,verytight,gdiis,maxstep=1,notrust,maxcyc=200) freq=ram
尝试替换换tight，verytight，nosymm，maxstep=1到5，int=superfine，sym=loose，结果都是这样的。

是什么原因导致的呢？是否能直接用第11帧的结构进行后续计算？
结构截图在下面、输出文件在附件里。

作者
Author: 倪宇晴 时间: 2025-1-9 14:05
多的一步是因为关键词里的 Freq。都opt=calcall了，如果只是做优化的话，freq不写也没关系。

作者
Author: Xulx2022 时间: 2025-1-9 14:37

倪宇晴发表于 2025-1-9 14:05
多的一步是因为关键词里的 Freq。都opt=calcall了，如果只是做优化的话，freq不写也没关系。

计算频率出现了不收敛，我应该怎么做？

作者
Author: 吴双 时间: 2025-1-9 16:27

Xulx2022 发表于 2025-1-9 14:37
计算频率出现了不收敛，我应该怎么做？

可以看看sob老师的这个帖子，http://sobereva.com/61

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-1-9 16:38
本帖最后由 Uus/pMeC6H4-/キ于 2025-7-2 18:37 编辑

单说问题本身，可以看下文。
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.keinsci.com/thread-633-1-1.html）

但是好奇一下，这个体系里平面三角形的Au6是哪来的，如果是切的金111表面那相对于分子尺寸过小不能反映晶面环境，如果就是要六个金原子的团簇那又感觉不像是构型搜索后众多极小点中能量最低的。编辑：现在回想起来，或许是要做表面增强拉曼光谱？这已经在http://bbs.keinsci.com/thread-19724-1-1.html有类似讨论了。

作者
Author: 倪宇晴 时间: 2025-1-10 20:47

Xulx2022 发表于 2025-1-9 14:37
计算频率出现了不收敛，我应该怎么做？

opt=verytight 把收敛限变严格了，我觉得可以算收敛了。可以试一下删掉verytight 和int=superfine。
另外，按照Gaussian 手册的介绍opt=verytight是搭配int=ultrafine用的。

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