| 查看: 5669 | 回复: 30 | |||||||||||||||
yjcmwgk禁虫 (文坛精英)
密度泛函·小卒
|
[交流]
【小卒个人文集】关于A+B生成AB的反应能和相互作用能的区别与计算方式的详细解释已有10人参与
|
|
这几天我的一个师妹说“不懂△E(interaction)和△E(reaction)的区别”，后来在QQ群上，又有人问我这个问题。 我寻思了一下，这个问题还真是“初学者不会，但是高手不屑一顾”的问题。小卒我作为“中手”，来系统解答这个问题吧。请高手斧正，请初学者批评。 其实这个问题是很简单的。讲述如下： （1）先对各个结构进行定义： 单独优化A，得到的结构是A1，其能量为E(A1)；单独优化B，得到的结构是B1，其能量为E(B1)。 优化复合物AB，得到的结构是X，对X进行BSSE校正后，得到其能量为E(X)。 （2）然后在结构X中，抽取原先属于A的那些原子），计算单点能，得到E(A2);抽取原先属于B的那些原子，计算单点能，得到E(B2)。 （3）计算： △E(reaction)=E(X)-E(A1)-E(B1)； △E(interaction)=E(X)-E(A2)-E(B2)； △E(A的变形)=E(A2)-E(A1)； △E(B的变形)=E(B2)-E(B1)； △E(reaction)=△E(interaction)+△E(A的变形)+△E(B的变形) （4）解释： 我们用图1和图2来解释这个问题。 图1表示了反应前和反应后的能量变化过程。2号线代表A分子在反应过程中，由于结构的变化，造成的分子能量的升高，也就是△E(A的变形)；3号线代表B分子在反应过程中，由于结构的变化，造成的分子能量的升高，也就是△E(B的变形)。1号线则是“2号线和3号线的累加”，也就是[△E(A的变形)+△E(B的变形)]，我们命名其为△E(Y)。从1号线我们可以得出结论：△E(Y)永远是正值——毕竟，从一个能量最优的全局最小点，无论向哪个方向偏移，都会造成能量升高。 图1中的4号线代表随着A和B的相互靠拢，A和B之间产生了相互作用，致使A和B之间的总势能降低（在不考虑A和B变形的前提下）。我们称它为△E(Z)。很显然，△E(Z)是负值（如果A和B之间以吸引力为主）或正值（如果A和B之间以排斥力为主）。如果△E(Z)>0，那么我们就根本没有考虑的必要了——因为他们根本不可能结合。因此，我们只考虑△E(Z)<0的情况。如果△E(Z)<0，那么，A和B之间能否结合，将取决于△E(Y)和|△E(Z)|谁更大。如果△E(Y)<|△E(Z)|，也就是图1中的5号线表示的那样，整个体系的总势能下降，致使整个体系稳定性增加，因此A和B可以结合为AB。相反，如果△E(Y)>|△E(Z)|，则证明A和B之间的吸引不足以克服A和B变形所带来的阻碍，A和B则无法结合为AB。 在mdbbs我同时发了这个帖子，有网友问： 不知道lz根据什么得出的 △E(interaction)=E(X)-E(A2)-E(B2)。 以下是Interaction energy（From Wikipedia） ‘ Supermolecular interaction energy’ A straightforward approach for evaluating the interaction energy is to calculate the difference between the energies of isolated objects and their assembly. In the case of two objects, A and B, the interaction energy can be written as: △Eint=E(A,B) -（E(A) +E(B) ） where E(A) and E(B) are the energies of the isolated objects (monomers), and E(A,B) the energy of their interacting assembly (dimer). The supermolecular approach has an important disadvantage in that the final interaction energy is usually much smaller than the total energies from which it is calculated, and therefore contains a much larger relative uncertainty. 我的解释如下： 维基百科似乎不是很科学吧？它也就比百度百科的“不严谨程度”稍好一点而已 维基百科的△Eint并非对应于我所说的△E(interaction) 维基百科的△Eint事实上对应于我所说的△E(reaction) 我在这个帖子中严格区分了interaction和reaction 在我的这个帖子中，interaction仅仅是A与B的相互作用——这个相互作用完全独立于A和B的变形所造成的能量变化。如果要考察A和B的由于变形所造成的能量变化，请看我的帖子中的△E(A的变形)和△E(B的变形)。 维基百科中，他所谓的interaction，并非只有A与B的相互作用，而是包含着A和B的变形所造成的能量变化（这种能量变化往往是不可忽视的）。所以，维基百科中所谓的△Eint远非真正的相互作用能。 笼统考虑，在某些时候，并不是科学态度（当然，我不是说全部）。 我们应该把各个因素分开考虑，这才有了我这里严格区分的△E(interaction)、△E(reaction)、 △E(A的变形)、△E(B的变形)、E(X)、E(A1)、E(B1)、E(A2)、E(B2)、E(Y)、E(Z)。我之所以分成这么多项目，就是为了将各个因素分开，而不是笼统的考虑它们 做科研应该有自己的思想，而不是人云亦云。别人说的话，是需要自己考虑一番才行的，而不是全盘接受。我在1楼发的帖子，很明显比维基百科上的说法要深入一些，毕竟维基百科上的定义比较笼统，而我将这个反应过程的两个因素（吸引和形变）分开来考虑。可能robustljb 先生发现维基百科上的符号△Eint和我的符号△E(interaction)在表面上是一致的，而没有深入考虑到维基百科上的符号△Eint和我的符号△E(reaction)才是本质上一致的 mdbbs的斑竹chenys的回答： 要考察与平衡结构有偏离的结构的能量变化时，都可以引入这个所谓的变形能，这不仅对超分子体系适用。 在一般的情况下interactin energy的确是按照维基那个定义处理的，其实这个定义更形象的说法应该是stablization energy。 可是stablization energy这个词在化学中有多种用法，比如 The amount by which the energy of a delocalized chemical structure is less than the theoretical energy of a structure with localized bonds. It is obtained by subtracting the experimental heat of formation of the compound (in kJ mol−1) from that calculated on the basis of a classical structure with localized bonds. interactin energy，deformation energy，stablization energy这些词要在论文中用还是给个定义比较好。 有网友又问： 以前只是想当然的认为△Eint=E(A,B)-E(A)-E(B) ，没有考虑吸引和形变这些因素，这次学习啦。不过还有两个问题不太明白 我觉得△E(reaction)中这个reaction使用欠妥，你把它认为一个反应还是别的什么呢？ 我的回答 我认为所有的相互作用都是广义上的反应。我给“反应”下的定义是：“在化学过程中，终态和始态相比，发生了变化”或“虽然终态和始态相同，但是化学过程中有体系的变化” 其实这个我所说的这个“反应”倒是跟“化学过程”类似了。是我不严谨，呵呵。谢谢您的提醒 有网友又问： 你在文字中说△E(reaction)=△E(interaction)+△E(A的变形)+△E(B的变形) 但是你的图像中却是△E(interaction)=△E(reaction)+△E(A的变形)+△E(B的变形) 我的回答 注意正负号 [ Last edited by yjcmwgk on 2010-6-16 at 21:17 ] |
回复此楼» 收录本帖的淘贴专辑推荐
 第一性原理 | 量子化学资料 | Gaussian | 我学习计算的一些帖子 |
| 量化计算有意义的讨论 | 计算化学 | 基础 | 气态分子化学动力学 |
| 自己感觉好的东西 | 高斯计算反应过渡态 | gaussian | Bao_S_J |
» 本帖已获得的红花（最新10朵）
coco207» 猜你喜欢
电催化或能源材料合作
已经有0人回复
-
k型热电偶的测温线出现两个地方短接，而只在一端有高温对测温有影响吗？
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有104人回复
-
[博士招聘]伍伦贡大学 University of Wollongong 王彩云老师课题组博士招聘
已经有0人回复
-
山东科技大学（青岛）化工学院招收还有调剂名额（仅限化学专业），速联系
已经有15人回复
-
华算科技诚聘理论计算工程师（兼职，有限元方向）
已经有0人回复
-
硅碳负极材料问题
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 能量变化解释化学反应的问题
已经有4人回复
- 量子化学参数△N解释
已经有4人回复
- 反应扩散和化学反应的区别和联系
已经有3人回复
- A和B分别生成C或D，如果判断C、D哪个先生成？
已经有4人回复
- 高分子近程相互作用与远程相互作用区别
已经有7人回复
- 修改意见反馈，关于表面相互作用
已经有5人回复
- 这几个化学反应能发生吗？生成什么呢
已经有4人回复
- 弱问：相互作用能负值越大应该描述成相互作用能越大还是越小
已经有15人回复
- 怎么能用量子化学计算含氮杂环芳香化合物与水分子的相互作用
已经有17人回复
- CP校正的必要性？
已经有8人回复
- 【求助】化学反应键能计算
已经有15人回复
- 【求助】间歇式和连续式反应的要求和区别
已经有13人回复
2楼2009-12-02 16:34:19
Joannaouc
木虫 (著名写手)
- 应助: 46 (小学生)
- 贵宾: 2.2
- 金币: 1686.5
- 散金: 894
- 红花: 10
- 帖子: 1262
- 在线: 266.1小时
- 虫号: 305787
- 注册: 2006-12-10
- 专业: 理论和计算化学
3楼2009-12-03 11:43:38
superfan
至尊木虫 (著名写手)
耕耘者
- 应助: 3 (幼儿园)
- 金币: 12221.2
- 红花: 7
- 帖子: 1816
- 在线: 699.2小时
- 虫号: 17588
- 注册: 2003-06-28
- 性别: GG
- 专业: 理论和计算化学
4楼2009-12-03 21:22:23
yrclq
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 249.2
- 散金: 157
- 帖子: 865
- 在线: 324.2小时
- 虫号: 529134
- 注册: 2008-03-20
- 性别: GG
- 专业: 催化化学
5楼2009-12-07 16:50:13
6楼2009-12-07 19:24:54
superfan
至尊木虫 (著名写手)
耕耘者
- 应助: 3 (幼儿园)
- 金币: 12221.2
- 红花: 7
- 帖子: 1816
- 在线: 699.2小时
- 虫号: 17588
- 注册: 2003-06-28
- 性别: GG
- 专业: 理论和计算化学
7楼2009-12-09 09:08:57
8楼2009-12-09 11:31:32
bianliang555551
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2162.5
- 散金: 3
- 红花: 4
- 帖子: 910
- 在线: 227.6小时
- 虫号: 425093
- 注册: 2007-07-22
- 性别: GG
- 专业: 理论和计算化学
9楼2009-12-09 18:31:29
10楼2009-12-10 11:50:11
10