请问Gaussian计算主客体分子单点能：用何基组泛函才能将分子间作用力考虑进去

nmhy0616

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

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

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

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

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

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

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

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

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

zjxitcc

（1）“能隙差”是什么？HOMO-LUMO gap大小 与 主客体结合能大小 没有必然关系。

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

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

pal

除了色散校正，标题写的是单点能，所以结构优化是如何进行的并没有说明

nmhy0616

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的全局最小处。

nmhy0616

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

感谢老师的回答

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

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

sobereva

nmhy0616 发表于 2026-5-21 17:01
感谢老师的回答

我的描述确实存在严重谬误。我们实际关注的并非HOMO-LUMO能隙，而是以主客体自由能面 ...

不要盲目拿理论计算孤立体系的情况和晶体衍射实验得到的结构对比，不具备严格可比性，仔细看
实验测定分子结构的方法以及将实验结构用于量子化学计算需要注意的问题
http://sobereva.com/569

B3LYP-D3(BJ)结合6-311+G(d,p)或6-311G(d,p)做优化，完全能合理描述分子间静电作用和色散作用

pal

nmhy0616 发表于 2026-5-21 17:15
感谢老师的回答

关于结构优化的操作补充说明如下：

组合后的结构在优化后会出现明显的变化说明当前结构不是局部最优，能量显然不会是最低的，去做单点也不会把他变成最低的。应该先加上色散校正后看看优化后的情况

nmhy0616

sobereva 发表于 2026-5-21 23:39
不要盲目拿理论计算孤立体系的情况和晶体衍射实验得到的结构对比，不具备严格可比性，仔细看
实验测定分 ...

好的老师，感谢老师回答

nmhy0616

pal 发表于 2026-5-22 09:23
组合后的结构在优化后会出现明显的变化说明当前结构不是局部最优，能量显然不会是最低的，去做单点也不会 ...

好的老师，我再去尝试一下，感谢老师的回答