计算化学公社

标题: 请问Gaussian计算主客体分子单点能：用何基组泛函才能将分子间作用力考虑进去 [打印本页]

作者
Author: nmhy0616 时间: 4 hour ago
标题: 请问Gaussian计算主客体分子单点能：用何基组泛函才能将分子间作用力考虑进去
本帖最后由 nmhy0616 于 2026-5-21 17:05 编辑

各位大佬好，新手求助，近期在计算一个主客体体系的单点能。主体为四阳离子紫精分子的环蕃，客体为含萘环的长链烷醚，主体需套在客体上形成一个轮烷结构。

从理论分析，四阳离子紫精环蕃应套在客体的萘环位置处，此时体系的自由能（Free Energy）应最低。这是因为紫精为缺电子结构（电子受体），萘环为富电子结构（电子供体），两者间的静电相互作用使之稳定。实验单晶结构也证实该构象最稳定。

然而在计算中，我首先使用以下设置：

# b3lyp/6-311+g(d,p) scrf=(smd,solvent=acetonitrile)

计算结果很不理想：主体位于客体萘环位置时，自由能反而是最高而非最低。与老师讨论后，认为该计算可能未合理考虑分子间静电作用。随后尝试了：

# M06-2X/6-311+G(d,p) scrf=(smd,solvent=acetonitrile)

得到的结果部分符合预期的趋势，但主体位于客体萘环位置时仍非是自由能最低点。

请问各位大佬，针对这种主客体之间的静电/供体-受体相互作用，是否有更合适的泛函和溶剂模型，能够更充分地描述其相互作用力，从而正确计算出FES最低点与实验一致？

第一次发帖，如果有不恰当的地方还请各位海涵，感谢指教！

作者
Author: zjxitcc 时间: 4 hour ago
（1）“能隙差”是什么？HOMO-LUMO gap大小与主客体结合能大小没有必然关系。

（2）B3LYP如果没有加色散校正，去做这种计算，结果大概率是完全错误的，这是做计算基本常识。我的意思是一开始就不应该提交计算，而不是算完之后再觉得有问题。如果您有给B3LYP加色散校正，却没有在帖子里写出来，也是不对的做法。建议详细阅读
谈谈“计算时是否需要加DFT-D3色散校正？”http://sobereva.com/413

（3）“整体仍然不令人满意”是什么意思？怎样是满意的？优化出来的主客体复合物结构与实验晶体结构差多少Angstrom以内吗？应当予以解释。

作者
Author: pal 时间: 2 hour ago
除了色散校正，标题写的是单点能，所以结构优化是如何进行的并没有说明

作者
Author: nmhy0616 时间: 1 hour ago

zjxitcc 发表于 2026-5-21 14:19
（1）“能隙差”是什么？HOMO-LUMO gap大小与主客体结合能大小没有必然关系。

（2）B3LYP如果没有加 ...

感谢老师的回答

我的描述确实存在严重谬误。我们实际关注的并非HOMO-LUMO能隙，而是以主客体自由能面（Free Energy Surface, FES, 单位 kcal/mol）为纵轴、以主体套在客体的不同相对位置为横轴的能量曲线。我之前错误地认为体系能隙差与自由能之间存在某种换算关系，这是不正确的认知。HOMO-LUMO能隙大小与主客体结合能大小并无必然联系，感谢老师指出这一问题。

第二个问题，B3LYP确实没有加色散校正，感谢老师指出的问题，我们对理论计算模拟方面涉猎确实较少，一些基础的计算理论常识了解不多，感谢老师分享的《谈谈“计算时是否需要加DFT-D3色散校正？”》我会认真学习这篇文章，避免今后再出现类似的低级错误。

第三个问题，整体仍不令人满意的原因是，单晶实验结果显示，主体套在客体萘环上构成最稳定构象，因此该结构所对应的自由能面上应处于全局最低点。然而在M06-2X泛函下的计算结果中，主体套在萘环位置虽已不再是最高的能量点，但依然不是全局最低点。我所期望的判定标准是：该结构对应的能量点应恰好落在FES的全局最小处。

作者
Author: nmhy0616 时间: 1 hour ago

pal 发表于 2026-5-21 15:42
除了色散校正，标题写的是单点能，所以结构优化是如何进行的并没有说明

感谢老师的回答

关于结构优化的操作补充说明如下：
我们首先对主体和客体分别进行了结构优化，采用的计算方法均为：
“#p opt freq b3lyp/6-311g(d,p) scrf=(solvent=acetonitrile)”
之后将优化好的主体与客体直接组合，进行单点能计算。

我们也曾将主体客体结合后进行过结构优化（但是并未通过分子动力学模拟进行分子对接，实在是条件有限），优化的过程中出现问题，就是在优化的时候主体在客体上的位置会随着结构优化而改变，暂时没有找到解决办法，后来组里老师说可以直接进行单点能的计算，不需要在专门优化，就没再深究这个问题

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