对碳材料吸附有机分子的体系进行结构优化（Gaussian16）的基组选择求助

我是路飞

    各位老师们好，我研究的体系是碳材料（简化模型为七元苯环，会分别引入不同官能团羧基、羟基等）分别吸附两种染料分子（亚甲基蓝阳离子、胭脂红阴离子）。我的最终目标是计算吸附能、结合能，能隙，以及绘制静电势图，以分析其吸附机理。目前我的计算条件是自己的电脑（win系统 14核 内存16GB）。现在想对吸附材料单体、染料分子、吸附模型分别进行结构优化，但对具体过程的基组选择不确定，同时对是否要对整个结构优化过程采用同样的计算方法存在疑问，具体情况如下。
    因为我的体系中的涉及弱相互作用，学习sob老师谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的 - 量子化学 (Quantum Chemistry) - 计算化学公社，决定采用B3LYP-D3(BJ)。
    同时我所研究体系中有阴离子，学习谈谈弥散函数和“月份”基组 - 思想家公社的门口：量子化学·分子模拟·二次元后，决定采用基组6-31+(d,p)。
    最终我的想法是，首先在对吸附材料以及染料分子这样的单体进行结构优化时都采用基组：DFT/B3LYP/6-31+(d,p)，再对吸附模型进行结构优化时采用基组：DFT/B3LYP-D3(BJ)/6-31+(d,p)。但是因为计算吸附能时同时涉及吸附质、基底、整个体系三者的能量，我想问倘如整个结构优化过程没有采用同样的基组是否会对实现我的最终目标产生影响，即是否需要整个结构优化过程采用DFT/B3LYP-D3(BJ)/6-31+(d,p)？
    请各位前辈百忙中解答我以上困惑，同时希望能指出我过程中存在的错误，十分感谢。
同时附上吸附质的结构图：

-HMDS-

建议再看一下下面的博文
http://sobereva.com/272

从我个人的经验出发，你算共轭体系的弱相互作用，B3LYP-D3显然不是最合适的，相对好用的是wB97系列的基组，如果用Gaussian的话就是wB97XD。

H的极化函数应该是不必要的。

比较合适的计算级别是wB97XD/6-31+G(d)

结构优化全程是不需要一直用大基组的，用小基组做初步优化再提高精度是常规操作。保证最后得到结构的计算基组方法保持一致就可以。

我是路飞

-HMDS- 发表于 2026-8-4 17:13
建议再看一下下面的博文
http://sobereva.com/272

好的，感谢老师指教

我是路飞

-HMDS- 发表于 2026-8-4 17:13
建议再看一下下面的博文
http://sobereva.com/272

老师，我还想问一下，Gaussian里wB97XD应该不能色散校正，我现在的目的是要得到我这个吸附体系的几何优化构型，并且计算吸附能，绘制静电势分布图。倘若仅仅在Gaussian使用wB97XD而不加色散校正，最终可否实现我的目的。倘若不可以，老师建议采用哪些方式以实现目的。

-HMDS-

我是路飞 发表于 2026-8-4 13:54
老师，我还想问一下，Gaussian里wB97XD应该不能色散校正，我现在的目的是要得到我这个吸附体系的几何优化 ...

你猜猜wB97XD的D是指的什么。不过我平时不用Gaussian，不太知道他的D是哪个版本，貌似是D2。D3和D4在ORCA才支持，不过wB97本身作为一个长程分离泛函，底子就已经比B3LYP好不少了。

sobereva

我是路飞 发表于 2026-8-4 20:54
老师，我还想问一下，Gaussian里wB97XD应该不能色散校正，我现在的目的是要得到我这个吸附体系的几何优化 ...

wB97XD自带DFT-D2校正，根本不需要而且也不能再加色散校正
wB97XD用于当前体系优化毫无问题，包括优化整体和各个部分
注意基组名必须严格按规范写，6-31+G(d,p)绝对不能写成6-31+(d,p)

我是路飞

-HMDS- 发表于 2026-8-4 21:06
你猜猜wB97XD的D是指的什么。不过我平时不用Gaussian，不太知道他的D是哪个版本，貌似是D2。D3和D4在ORCA ...

理解了，感谢老师

我是路飞

sobereva 发表于 2026-8-5 00:56
wB97XD自带DFT-D2校正，根本不需要而且也不能再加色散校正
wB97XD用于当前体系优化毫无问题，包括优化整 ...

感谢老师指教，书写规范问题日后一定注意