计算化学公社

标题: 高斯计算离子液体与R116之间的结合能 [打印本页]

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

作者
Author: 睡觉不打鼾 时间: halfhour ago
请仔细阅读：
谈谈量子化学中基组的选择 - 量子化学 (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不收敛问题的方法 - 思想家公社的门口：量子化学·分子模拟·二次元

另外你这离子液体对吗？

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