计算化学公社

标题: 杂多酸阴离子结构优化SCF不收敛 [打印本页]

作者
Author: zch22076027 时间: 2026-6-25 14:39
标题: 杂多酸阴离子结构优化SCF不收敛
本帖最后由 zch22076027 于 2026-6-25 14:41 编辑

各位老师好，本人在利用Gaussian16进行SiW11O39的优化时，遇到了SCF不收敛问题，研读了sob老师的解决SCF不收敛问题的方法 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)，以及令Gaussian 16中SCF未收敛到默认收敛限也能继续做后续计算的方法 - 思想家公社的门口：量子化学·分子模拟·二次元等几篇博文尝试了其中几种方法后，依旧难以收敛，现在不知道怎么调整了，请各位老师帮忙检查一下！

作者
Author: zjxitcc 时间: 2026-6-25 16:30
绝对不要使用IOp(7/127=-99)，您已中毒，错误理解博文。不要用conver=6。

您针对的是一个真实体系，体系里往往有抗衡离子（对当前体系，也就是阳离子），应当把对应的阳离子补上去，不要使用-8电荷进行计算。如果是在溶液中，还应当恰当设置隐式溶剂模型，不要偷懒当成气相做计算。

作者
Author: zch22076027 时间: 2026-6-25 17:09
本帖最后由 zch22076027 于 2026-6-25 17:14 编辑

zjxitcc 发表于 2026-6-25 16:30
绝对不要使用IOp(7/127=-99)，您已中毒，错误理解博文。不要用conver=6。

您针对的是一个真实体系，体系 ...

好的老师。感谢您的回复。这个体系的阴离子是一个四己基氢氧化铵（C24H53NO），带1个正电荷，如果要全部平衡的话需要8个阳离子，将导致体系太大。并且合作者认为不能缩短碳链，因为碳链的长度对实验性能很重要。因此现在硬着头皮在计算。这不是在溶液中的，所以没有加隐式溶剂模型。

作者
Author: zjxitcc 时间: 2026-6-25 17:17
本帖最后由 zjxitcc 于 2026-6-25 17:35 编辑

zch22076027 发表于 2026-6-25 17:09
好的老师。感谢您的回复。这个体系的阴离子是一个四己基氢氧化铵（C24H53NO），带1个正电荷，如果要全部 ...

我仍然建议把阳离子加上，然后使用半经验方法GFN2-xTB做结构优化和频率分析（xTB程序中有很多半经验方法，GFN2-xTB是GFN系列里目前最准的）。

如果后续还需其他计算，可以选择合适的泛函和基组，用量化程序ORCA做DFT单点计算，并开启RI加速。ORCA在RI加速的情况下，算几百个原子、上千个原子体系的DFT单点犹如砍瓜切菜（注意，请不要回复我“没有ORCA”、“没有ORCA版权”、“只有Gaussian版权”之类的话，因为ORCA是学术免费的量化软件，在审稿人眼中看来这是你已经拥有的东西）。

作者
Author: zch22076027 时间: 2026-6-25 17:40

zjxitcc 发表于 2026-6-25 17:17
我仍然建议把阳离子加上，然后使用半经验方法GFN2-xTB做结构优化和频率分析（xTB程序中有很多半经验方法 ...

好的，再次感谢您的帮助！我还想请教一下，不收敛的原因是因为负电荷太多了，还是初始的结构不太合理。

作者
Author: zjxitcc 时间: 2026-6-25 17:45

zch22076027 发表于 2026-6-25 17:40
好的，再次感谢您的帮助！我还想请教一下，不收敛的原因是因为负电荷太多了，还是初始的结构不太合理。

负电荷太多；在有负电荷的情况下，基组却没有考虑弥散函数。

作者
Author: zch22076027 时间: 2026-6-25 22:44

zjxitcc 发表于 2026-6-25 17:45
负电荷太多；在有负电荷的情况下，基组却没有考虑弥散函数。

我之前也试过给O加弥散，但是也没有效果。总之谢谢您。我先按照您说的试试吧。

作者
Author: zjxitcc 时间: 2026-6-25 23:26

zch22076027 发表于 2026-6-25 22:44
我之前也试过给O加弥散，但是也没有效果。总之谢谢您。我先按照您说的试试吧。

只加弥散函数是不够的，需要多条措施并举，例如保持xqc，增大SCF迭代圈数，以及上面所说的绝对不要使用IOp(7/127=-99)等等。另外，最好只给O加弥散函数，而不是给所有原子都加弥散函数。

作者
Author: zch22076027 时间: 2026-6-25 23:39

zjxitcc 发表于 2026-6-25 23:26
只加弥散函数是不够的，需要多条措施并举，例如保持xqc，增大SCF迭代圈数，以及上面所说的绝对不要使用IO ...

好的谢谢老师。

作者
Author: wzkchem5 时间: 2026-6-26 09:26

zch22076027 发表于 2026-6-25 17:09
好的老师。感谢您的回复。这个体系的阴离子是一个四己基氢氧化铵（C24H53NO），带1个正电荷，如果要全部 ...

如果碳链的长度对实验性能很重要，碳链长度不对都不行，那完全不加阳离子岂不是更不行

作者
Author: zch22076027 时间: 2026-7-3 10:21

wzkchem5 发表于 2026-6-26 09:26
如果碳链的长度对实验性能很重要，碳链长度不对都不行，那完全不加阳离子岂不是更不行

感谢老师回复。当时没考虑清楚，想的是阴阳离子分开计算，体系小点，但是没考虑难收敛的问题。现在想起来确实有问题。尝试众多方法无果后，换软件计算了

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