计算化学公社

标题: 求助有机TADF分子的SCF不收敛解决办法尝试 [打印本页]

作者
Author: lindaiyu 时间: 2026-3-14 09:58
标题: 求助有机TADF分子的SCF不收敛解决办法尝试
各位老师好，我在计算如下有机TADF分子的结构和频率时，最初选用#p opt freq b3lyp/6-311++g(d,p) scrf=(smd,solvent=toluene) ，SCF不收敛，看了帖子（http://sobereva.com/61）之后自己想了两种方法，分别计算了两个分子：
方法1.#p opt freq b3lyp/6-311++g(d,p) scrf=(smd,solvent=toluene) SCF(novaracc,noincfock) int=acc2e=12，按这个方法在超算平台结果已计算出来
方法2.#p opt freq b3lyp/6-311++g(d,p) scrf=(smd,solvent=toluene) scf=(conver=6,vshift=500)。还在计算ing...
tadf分子结构如下图，方法一对应分子元素为C、H、O、N，Cl，方法二对应元素为C、H、O、N、Cl、S。
(, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 4)

想问一下老师哪种方法要好一些，精度更高？

作者
Author: fatpig 时间: 2026-3-14 10:49
你这计算加弥散有啥意义，如果一定要加弥散，那给 Cl 单独加就行了，连 H 都加上的目的是啥...

加弥散通常就是SCF难收敛的一个重要原因...

作者
Author: lindaiyu 时间: 2026-3-14 10:53
本帖最后由 lindaiyu 于 2026-3-14 02:54 编辑

fatpig 发表于 2026-3-14 02:49
你这计算加弥散有啥意义，如果一定要加弥散，那给 Cl 单独加就行了，连 H 都加上的目的是啥...

加弥散通 ...

原来是这样啊，我不清楚哪里加弥散哪里不加，我再去看看社长那几篇基组泛函的文章，非常感谢！

作者
Author: zjxitcc 时间: 2026-3-14 15:52
本帖最后由 zjxitcc 于 2026-3-14 15:55 编辑

跟SCF不收敛没有关系，您搞错了问题，错症下药。
（1）结构非常不合理，有几个原子看起来距离过近。在GaussView里点击Clean一下，确保结构合理一些，然后用半经验量子化学方法GFN2-xTB或PM7先优化一遍，把优化好的结构拿来做DFT结构优化。
（2）基组非常不合理，急需阅读Sob老师这篇博文
谈谈量子化学中基组的选择http://sobereva.com/336
（3）B3LYP没有加色散校正，等于白算。急需阅读
谈谈“计算时是否需要加DFT-D3色散校正？” http://sobereva.com/413
DFT-D色散校正的使用 http://sobereva.com/210

作者
Author: sobereva 时间: 2026-3-15 06:14
opt freq不建议用SMD，下文说了
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）
对于甲苯这么低极性的溶剂，完全可以近似为真空。
opt freq任务也不建议用conver=6这么宽松，下文说了
解决SCF不收敛问题的方法
http://sobereva.com/61

此体系此任务加弥散完全多余，不仅无意义地徒增耗时还显著增加SCF收敛的难度，属于自讨苦吃。

但当前导致SCF不收敛的最关键原因还是初始结构忒烂，截图里都能看出很多地方自动判断的成键关系都明显不对，哪怕收敛了也没意义，纯粹是瞎算

作者
Author: lindaiyu 时间: 2026-3-16 10:29

sobereva 发表于 2026-3-14 22:14
opt freq不建议用SMD，下文说了
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

感谢sob老师，我把初始结构调好了，先用GFN2-xTB算一遍，然后用#p opt freq b3lyp/6-311g(d,p) em=GD3BJ scrf=(read,solvent=toluene) 算，非常感谢1

作者
Author: lindaiyu 时间: 2026-3-16 10:39

zjxitcc 发表于 2026-3-14 07:52
跟SCF不收敛没有关系，您搞错了问题，错症下药。
（1）结构非常不合理，有几个原子看起来距离过近。在Gaus ...

谢谢您，我把初始结构调好了，先用GFN2-xTB算，然后再用#p opt freq b3lyp/6-311g(d,p) em=GD3BJ scrf=(read,solvent=toluene) 算，非常感谢！

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