| 查看: 12313 | 回复: 187 | |||||||||||
| 【奖励】 本帖被评价155次，作者lihb734增加金币 121.4 个 | |||||||||||
| 本帖产生 1 个 QC强帖 ，点击这里进行查看 | |||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||||||
[资源]
【lihb个人文集】浅谈BSSE矫正之相互作用能的计算
|
|||||||||||
|
一、原理与理解 在量子化学中，由于采用的是有限基组近似，故在计算一个超分子比如二聚体AB的相互作用能时，会由于描述片段A的基组扩散到B上，反之亦然。这样，描述片段A,B,二聚体分子的基组就不再处于同一水平了——基函数的数目发生一定的变化。这样会导致计算二聚体能量时所用的基函数数目增多，人为降低了二聚体的能量。从而导致相互作用能的计算 DeltaE=Ea+Eb-Eab 偏大（取正值）。这种现象就成为基组重叠矫正误差（basis set ）也可以这么理解，由于之前提到的目前使用的基组是有限的，非完备的，根据变分原理：决定能量的自由度增加，能量降低。这样导致二聚体的能量与没有这种基组叠加误差的能量相比时偏低，导致DeltaE偏大。 我们也可以还一种思路来理解：我们知道，描述分子的基函数数目越多，得到的能量准确，这样我们认为Eab的计算是准确的。但Ea与Eb计算时，与Eab相比，所用的基组数目要少一些，故导致前二者能量偏高，这样算出来的DeltaE也偏大。 这是一种人为效应，所以需要矫正，于是1970年Boys等人(S. F. Boys and F. Bernardi, “Calculation of Small Molecular Interactions by Differences of Separate Total Energies - Some Procedures with Reduced Errors,” Mol. Phys., 19 (1970) 553.)提出了couterpoise矫正的思路，在单体上增加基函数，来降低这种人为的偏高，以达到矫正的目的。 二、高斯实例与计算 在高斯中计算DeltaE加入关键词Counterpoise就可以了，若有两个片段则指定：Counterpoise=2 即可。 输入文件例子： # UHF/6-21G(d) Counterpoise=2 Counterpoise on water dimer 1,2 1,2 0,1 O(Fragment=1) 0.00 0.00 0.00 O(Fragment=2) 0.00 0.00 2.98 H(Fragment=1) 0.49 0.76 -0.29 H(Fragment=1) 0.49 -0.76 -0.29 H(Fragment=2) -0.91 0.00 3.24 H(Fragment=2) -0.01 0.00 2.03 三、结果处理 成功计算以后，会得到七个相关的能量：搜索scf done会得到五个能量1-5，分别为：（数字是能量顺序的编号，便于下边的计算） 未校正的二聚体能量 1 校正后的片段1的能量 2 校正后的片段2的能量 3 未校正的片段1的能量 4 未校正的片段2的能量 5 两外两个是： 6 Counterpoise: corrected energy = -2510.716555327645 7 Counterpoise: BSSE energy = 0.004963814828 这两个分别为矫正之后二聚体的能量； BSSE重叠能，也就是片段A和片段B被高估的能量，或者二聚体被低估的能量。 这样我们计算矫正之后的相互作用能应该至少有三种思路： 1）2+3-1 此时把二聚体的能量看成是准确计算的，因为它的变分自由度最多。所以采用矫正的片段能量。 2）4+5-6 这时我们把片段的能量堪称是可以准确计算的，此时二聚体的能量计算偏低，故采用矫正的二聚体能量6. 3）4+5-1-7 其思路和1）差不多，这不过被高估的片段A ，B的能量最后才被扣除 因此，从上述可以看出，最偷懒的方法是找到前三个能量就可以简单进行计算。 四、注意 1）根据基组的完备度，显然基组越小，BSSE效应越明显，越有必要矫正；基组越大时，矫正的必要性越小，这也是Grimme主张的在进行DFT-D计算相互作用能时，采用TZ2P以上的基组可以考虑不矫正的原因； 2）弥散函数比极化函数更加扩散，这样弥散函数的BSSE效应要强于极化函数的； 3) 对于非SCF水平下（HF和DFT）的计算，请务必选择正确的能量，比如采用MP2计算时，此时scf done后给出的能量只是HF水平下的能量，真正的MP2水平下的能量是：EUMP2= 后边的数值； 4）计算过程中，通常先优化二聚体，或者直接从晶体结构中取二聚体的几何即可。不要试图去优化片段，这样得到能量不再是相互作用能，此时还包括了片段形变后的能量。 五、尾声 为了帮助大家熟悉BSSE和计算相互作用能。本人整理了一个excel表格，只要填入相应的数字就可以得到相关能量。并结合了一个例子输出列出这七个能量所在的行号。 六、免责声明 本文是根据我对量化基本理论的理解整理的，难免有措辞不当，理解不准确的地方。请您看完之后予以批评指正。谢谢！ [ Last edited by lihb734 on 2011-6-1 at 16:49 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 量子化学资料 | 我的量化学习 | 高斯/MOPAC | 知识能改变命运吗？ |
| 计算各种马克 | 经验 | 计算化学 | 个人收藏 |
| 出国留学 |
» 本帖已获得的红花（最新10朵）
yongma2008» 猜你喜欢
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有277人回复
-
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
-
英国利物浦大学屠昕和雍雪课题组招收CSC博士生
已经有27人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）
已经有0人回复
-
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人
已经有45人回复
-
英国贝尔法斯特女王大学招收 2026 年 CSC/BBSRC phd 机器学习 催化剂设计方向
已经有0人回复
-
法国里昂大学招生CSC 资助博士生
已经有0人回复
-
法国里昂大学接收CSC 资助访问学者
已经有0人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 谈谈谐振频率校正因子
已经有9人回复
- 零点能/结构优化关键词
已经有18人回复
- 是否可以不用freq算零点能
已经有14人回复
- 关于gaussian09中BSSE校正
已经有10人回复
- 如何计算氢键键能
已经有7人回复
- 这样算的高斯 准不准确
已经有16人回复
- 考虑BSSE的结构优化问题
已经有8人回复
- BSSE校正中 片段2的变化如何改写
已经有8人回复
- 我们用VASP计算材料的力学性能，是绝对零度的。有意义吗？正解者有红花相送哦
已经有8人回复
- 计算两个分子相互作用
已经有14人回复
- 高斯计算中离子电荷的处理
已经有15人回复
- CP校正的必要性？
已经有8人回复
- 【求助】关于BSSE校正：优化后进行BSSE校正码？
已经有6人回复
42楼2011-07-15 17:02:21
3楼2011-05-05 21:41:29
11楼2011-05-06 08:36:48
简单回复
8tao74楼
2011-05-05 21:41
回复
五星好评 
2011-05-05 21:43
回复
2011-05-05 21:43
回复
五星好评 
xuyy20107楼
2011-05-05 21:54
回复
五星好评
daiyulan8512楼
2011-05-06 08:45
回复
五星好评 顶一下，感谢分享！
duanya06113楼
2011-05-06 09:38
回复
五星好评
Quan.14楼
2011-05-06 10:50
回复
五星好评 顶一下，感谢分享！