计算化学公社

 找回密码 Forget password
 注册 Register
Views: 164|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 求助:Gaussian计算相互作用使用CP法进行BSSE矫正是否还必须加弥散函数?

[复制链接 Copy URL]

5

帖子

1

威望

288

eV
积分
313

Level 3 能力者

各位老师好,最近我在进行双三氟甲烷磺酰亚胺锂(锂离子和盐阴离子)和碳酸乙烯碳酸乙烯亚乙酯(VEC)等一系统不同碳酸酯的相互作用能的计算。对此我有一些疑问想向各位老师请教:

我所研究的体系并不需要考虑分子内的相互作用,目前采用的是def2-TZVP。计算阴离子和碳酸酯的复合体系因为阴离子存在应该是要求加弥散函数的,目前考虑的是用带D后缀的def2系列(def2-TZVPD),但是加弥散后基本上不太算的动(CPU为9950x,目前也很难进行设备升级或上超算完成)


学习sob老师发布的大体系弱相互作用计算的解决之道http://sobereva.com/214”时候发现其中提到:“之所以不含弥散函数的小基组在计算弱相互作用上结果差,其实关键因素就在于BSSE。普遍提倡在计算弱相互作用时加上弥散函数,其实从实效上讲就是在于减轻BSSE。换句话说,对于分子间相互作用的计算,弥散函数起到的作用和CP实际上是相仿佛的。.....所以说,如果用CP来解决掉BSSE问题,那么对于分子间相互作用计算就不太需要弥散函数了(虽然再加上弥散的话结果还能再改善一点点)。


这是否说明即使不加弥散函数,单独靠Counterpoise矫正,在def2-TZVP这种较大的3-zeta的基组(谈谈BSSE校正与Gaussian对它的处理http://sobereva.com/46中提到:“对于3-zeta档次基组,def2-TZVP和6-311+G(2d,p)算弱相互作用能能及格,但考虑counterpoise后才能足够理想”),也能得到定性正确的可接受的结果?还是说在阴离子存在的体系无论什么情况下都一定要加弥散函数?如果是后者的话,请教各位老师对我目前的情况是否有比较好的解决方案,能解决opt和频率计算过程中耗时难以接受的问题

127

帖子

0

威望

2765

eV
积分
2892

Level 5 (御坂)

2#
发表于 Post on 5 day ago | 只看该作者 Only view this author
你先用 def2-svpd 或者 ma-svp 优化和算频率,后续用def2-tzvpd 或者 ma-tzvp 算单点就行了...

5

帖子

1

威望

288

eV
积分
313

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 5 day ago | 只看该作者 Only view this author
fatpig 发表于 2026-2-11 10:03
你先用 def2-svpd 或者 ma-svp 优化和算频率,后续用def2-tzvpd 或者 ma-tzvp 算单点就行了...

感谢fatpig老师的回复,这个情况我之前考虑过,主要担心的点如下:相互作用的结果应该跟几何优化的结果关联很大吧,如果用2-zeta基组做几何优化,是否会导致就算后续CP单点换成3-zeta加弥散的基组,其结果也会有较大的偏差?还是说做了BSSE矫正后,就算全程(包括CP单点)是2-zeta基组,也能得到相对合理的数值呢?
上述这个问题之前在学习sob老师博文的时候就没有太明确,还请各位老师解惑

409

帖子

3

威望

2608

eV
积分
3077

Level 5 (御坂)

4#
发表于 Post on 5 day ago | 只看该作者 Only view this author
本帖最后由 Stardust0831 于 2026-2-11 15:22 编辑
Monsoon 发表于 2026-2-11 15:05
感谢fatpig老师的回复,这个情况我之前考虑过,主要担心的点如下:相互作用的结果应该跟几何优化的结果关 ...

结构对基组的依赖性比较低,2zeta下算的结构大部分时候都够用了。
带D后缀的def2系列(def2-TZVPD)是给极化率、偶极矩这种高度需要弥散函数的场合用的。
def2-TZVP本身就有一定半弥散特征了,在orca中结合近乎免费的gCP矫正来考虑bsse是最优选择。
如果有负电荷特别集中的区域,可以用ma-def2-TZVP结合counterpoise,也足够完美。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
SharkYYX2025 + 2 谢谢

查看全部评分 View all ratings

5

帖子

1

威望

288

eV
积分
313

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 5 day ago | 只看该作者 Only view this author
Stardust0831 发表于 2026-2-11 15:17
结构对基组的依赖性比较低,2zeta下算的结构大部分时候都够用了。
带D后缀的def2系列(def2-TZVPD)是给 ...

感谢老师解惑!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-16 20:22 , Processed in 0.155782 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list