求助高斯计算离子液体与R116之间的结合能，计算很慢还会出现不收敛的情况

Flies68

各位老师们，求助一下关于高斯计算离子液体与R116之间结合能的问题。如图，是我的构型与gif文件（离子液体与R116都分别进行了优化），但计算很慢还会出现不收敛的情况。是我构型出了问题还是该用其他方法来计算。
当然我也试过用溶剂化的方法，根据http://sobereva.com/431中给出的方法很容易便可算出溶解自由能。但其给出的介电常数都是平均值，计算后也不知道R116气体与离子液体的哪一部分存在结合力，计算后得出R116气体的溶解自由能也无从考证。
各位大佬们有什么高见吗？欢迎留言

睡觉不打鼾

请仔细阅读：
谈谈量子化学中基组的选择 - 量子化学 (Quantum Chemistry) - 计算化学公社
“单点能计算、几何优化（包括优化到过渡态）、振动分析是最常见的任务类型。几何优化和振动分析对基组的敏感度远低于单点能的计算，而这些任务的耗时又远高于单点计算，所以几何优化绝对用不着大基组，诸如用6-31G**这样偏小的基组一般就可以，至多用到6-311G**或def-TZVP这样中等尺寸基组也就够了，用诸如def2-TZVP这样档次的基组已经属于明显浪费；等优化完，算单点能的时候再把基组提升一、两个档，是明智的做法。这属于量子化学研究者的最基本常识，完全不必担心有哪个内行人会对此有质疑，详见《浅谈为什么优化和振动分析不需要用大基组》（http://sobereva.com/387）。”

以及：
谈谈弥散函数和“月份”基组 - 思想家公社的门口：量子化学·分子模拟·二次元
“所以，需要对体系加弥散函数时，强烈不建议用诸如6-311++G**，而建议用6-311+G**。对于Dunning系列基组要加弥散的情况，也建议用后文说的月份基组而不是aug-cc-pVnZ。即便研究氢键键能也不需要给氢加弥散函数，测试表明这对结果的影响微乎其微。但对于计算(超)极化率、氢带明显负电（如LiH）等情况，给氢加弥散函数确实有重要意义，此时需要给氢加弥散。”

还有：
解决SCF不收敛问题的方法 - 思想家公社的门口：量子化学·分子模拟·二次元

另外你这离子液体对吗？

mizu-bai

建模都有问题啦，先想清楚自己要研究的是什么体系

sobereva

那叫gjf不叫gif

说清楚是SCF不收敛还是几何优化不收敛，仔细看下文
解决SCF不收敛问题的方法
http://sobereva.com/61
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

没事甭给氢加弥散，下文明确说了
谈谈弥散函数和“月份”基组
http://sobereva.com/119

用区区四核，还是windows版，做耗时的opt freq任务，还不知道监控计算过程，没法不慢

另外，算这种问题还应当做构型搜索，否则结果绝对招内行的质疑，完全不是随便自己摆个初始构型优化就行的。看
Molclus主页
http://www.keinsci.com/research/molclus.html
Molclus FAQ
http://sobereva.com/730（http://bbs.keinsci.com/thread-50349-1-1.html）
使用molclus程序做团簇构型搜索和分子构象搜索
http://bbs.keinsci.com/thread-577-1-1.html
genmer：生成团簇初始构型结合molclus做团簇结构搜索的超便捷工具
http://bbs.keinsci.com/thread-2369-1-1.html

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “高斯计算离子液体与R116之间的结合能” 改了，以后务必注意！

Flies68

mizu-bai 发表于 2026-4-21 08:00
建模都有问题啦，先想清楚自己要研究的是什么体系

老师您好，我是新手小白，对建模不太熟悉，烦请您告诉我一下这个建模的具体问题是什么，有没有什么推荐的学习规范建模的渠道？