计算化学公社

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

[Gaussian/gview] 为什么将一个结构拆成三个部分计算后,Total Basis Functions之和与整体不一致

[复制链接 Copy URL]

69

帖子

0

威望

471

eV
积分
540

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 typ_gauss 于 2025-11-20 19:23 编辑

请教大家一个问题:
整体结构141个原子,含有C、H、O、N、P、Br5种元素,整体采用B3LYP-D3BJ/6-311+G**(SDD for Br)/SMD,DCM的方法计算,其Total Basis Functions为2075。
将其拆解为3个部分后,每个部分分别含有86、41、14个原子,其中14个原子的部分含Br,因此只对这部分采用B3LYP-D3BJ/6-311+G**(SDD for Br)/SMD,DCM的方法,其他两个部分采用B3LYP-D3BJ/6-311+G**/SMD,DCM,结果三个部分的Total Basis Functions加和为1531,与整体的计算完全不一致。
这是为什么呢,是因为混合基组的原因吗?整体我是用的b3lyp/genecp,而拆解后由于两个组分不含有溴,则直接使用b3lyp/6-311+G**,是这个原因吗?谢谢!
用B3LYP-D3BJ/6-31+G*(SDD for Br),则整体和部分拆解的能对上。


6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

2#
发表于 Post on 2025-11-20 03:02:27 | 只看该作者 Only view this author
分清楚笛卡尔和球谐型基函数,看:
谈谈5d、6d型d壳层基函数与它们在Gaussian中的标识
http://sobereva.com/51

Gaussian里用gen或genecp的时候总是用球谐型基函数,直接指定某些基组时对d壳层默认用笛卡尔型基函数
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

69

帖子

0

威望

471

eV
积分
540

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-11-20 14:38:18 | 只看该作者 Only view this author
sobereva 发表于 2025-11-20 03:02
分清楚笛卡尔和球谐型基函数,看:
谈谈5d、6d型d壳层基函数与它们在Gaussian中的标识
http://sobereva.c ...

sob老师,感谢回复,我的关键词是加了5d的,并且后面我换了个没有Br元素的结构,直接用的b3lyp/6-311+g**,没有用genecp了,它依旧不一致,感觉是6-311+g**的原因。
比如整体部分用的# sp b3lyp/6-311+g** empiricaldispersion=gd3bj 5d scrf=(smd,solvent=dichloromethane) geom=check guess=read
拆开的部分也是用的# sp b3lyp/6-311+g** empiricaldispersion=gd3bj 5d scrf=(smd,solvent=dichloromethane) geom=check guess=read
然后计算出来发现basis functions整体和部分加和仍然不一致,但看了您的文章后,我又统计了一下primitive gaussians,这个确实是一致的,如果换成6-31g*则basis functions和primitive gaussians则都一致

1208

帖子

6

威望

2432

eV
积分
3760

Level 5 (御坂)

傻傻的木瓜

4#
发表于 Post on 2025-11-20 16:34:11 | 只看该作者 Only view this author
对照手册https://gaussian.com/basissets/?tabid=1与讨论帖http://bbs.keinsci.com/thread-19095-1-1.html来看,在5d/6d和直接指定/按混合基组指定上有坑的是2-zeta的6-31G(及其带极化和弥散的版本),而我没记错的话3-zeta的6-311G(及其带极化和弥散的版本)默认就是5d的。

这里更关键的是基组有弥散函数,可能在某个阈值下对不同体系去除线性相关基函数的数量不同,这就得观察输出文件里原本的基函数数目(以NBasis标记)和去除线性依赖后剩下的基函数数目(以NBsUse标记)了。
√546=23.36664289109

69

帖子

0

威望

471

eV
积分
540

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2025-11-20 19:17:07 | 只看该作者 Only view this author
本帖最后由 typ_gauss 于 2025-11-20 19:22 编辑

使用 b3lyp/6-311+g**的结果
片段1:
NBasis=  1308 RedAO= T EigKep=  1.07D-06  NBF=  1308
NBsUse=  1289 1.00D-06 EigRej=  7.21D-07 NBFU=  1289
1308 basis functions,  2074 primitive gaussians,

片段2:
NBasis=   566 RedAO= T EigKep=  2.73D-06  NBF=   566
NBsUse=   566 1.00D-06 EigRej= -1.00D+00 NBFU=   566
566 basis functions,   888 primitive gaussians,

片段3:
NBasis=   604 RedAO= T EigKep=  1.03D-06  NBF=   604
NBsUse=   596 1.00D-06 EigRej=  8.41D-07 NBFU=   596
604 basis functions,   962 primitive gaussians,


整体:
NBasis=  2478 RedAO= T EigKep=  1.02D-06  NBF=  2478
NBsUse=  2433 1.00D-06 EigRej=  9.80D-07 NBFU=  2433
2478 basis functions,  3924 primitive gaussians,

不好意思,给大家添麻烦了,我之前统计错了,使用b3lyp/6-311+g**的basis functions是一致的,抱歉。之前应该是genecp造成的,我有些地方加了5d,有些没加,有些混乱,所以写基组还是保持一致比较好,感谢大家!您提到的NBsUse确实值得思考,片段的加和(1289+566+596=2451)和整体(2433)不一致.

本版积分规则 Credits rule

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

GMT+8, 2026-1-24 08:47 , Processed in 0.226231 second(s), 21 queries , Gzip On.

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