计算化学公社

标题: 请问在使用Gaussian09做(ZnO- CH3(CH2)7OH)4的几何优化时一直不收敛该怎么解决 [打印本页]

作者
Author: yaoyaoyao哥 时间: 2025-6-10 10:38
标题: 请问在使用Gaussian09做(ZnO- CH3(CH2)7OH)4的几何优化时一直不收敛该怎么解决
首先使用了opt=(tight,cartesian)进行优化，报错l9999；
然后在上一步的结构基础上使用opt=GDIIS进行优化，后又报错l502；
之后看了sob老师的文章SCF不收敛时如何解决，又尝试使用以下信息在上一步的结构基础上进行计算#p M052X/GenECP SCF=VShift=500 opt freq，还是报错l502，以下文件是此步骤的输出文件。
之后我还尝试使用B3LYP进行计算，还是报错了l9999，请问各位老师还有什么方法可以解决。

作者
Author: sobereva 时间: 2025-6-10 12:12
绝对不该用M05-2X算这个，认真看
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）
改用TPSSh-D3(BJ)或PBE0-D3(BJ)

认真看
解决SCF不收敛问题的方法
http://sobereva.com/61
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164
要认真、细致、完整看全文，不要断章取义

“报错了l9999”这种是极为幼稚的提问方式。不懂怎么恰当描述问题就上传输入文件和压缩后的输出文件
但凡用了GenECP，要么上传输入文件，要么准确描述清楚用的基组和赝势，否则显然根本没有提供关键信息

作者
Author: yaoyaoyao哥 时间: 2025-6-10 14:08

sobereva 发表于 2025-6-10 12:12
绝对不该用M05-2X算这个，认真看
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs ...

好的谢谢老师，我再认真学习一下，使用的基组和赝势为6-31+G(d,p)和LANL2DZ

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-6-10 16:49
本帖最后由 Uus/pMeC6H4-/キ于 2025-6-10 17:06 编辑

我觉得这个Zn20和Zn24仅三配位的结构不太对，或许像下面这样让四个锌和四个羟基形成团簇的结构会合理些。图为CYLview制作。
(, 下载次数 Times of downloads: 16)
用PBE0-D3(BJ)泛函结合def2-TZVP与def2-SVP混合全电子基组，在S4点群对称性下做几何优化的输入输出文件如下所示。其中def2-TZVP用于核心的锌、氧等1-16号原子，def2-SVP用于周围四个正辛基的碳与氢。此处柔性分子因最终最大受力足够小而判定为几何优化收敛，未做进一步振动分析。写错的title行来自前一步优化，忽略就好
(, 下载次数 Times of downloads: 2)
编辑：上面的压缩包是tar -cvJf压缩的，但是后缀从xz打错成gz了（尴尬）解压的时候应用tar -xvJf来处理。

作者
Author: yaoyaoyao哥 时间: 2025-6-10 17:00

Uus/pMeC6H4-/キ发表于 2025-6-10 16:49
我觉得这个Zn20和Zn24仅三配位的结构不太对，或许像下面这样让四个锌和四个羟基形成团簇的结构会合理些。图 ...

好的，谢谢老师

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