这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
adf:nocv-disdinguishbonding2020 [2020/11/22 11:11] – [参数设置] liu.jun | adf:nocv-disdinguishbonding2020 [2023/11/27 17:35] (当前版本) – [EDA-NOCV功能案例(分子间作用):二聚体能量分解(EDA)中轨道作用、电子在片段轨道间转移的分解] liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ======EDA-NOCV功能案例(闭壳层体系):二聚体能量分解(EDA)中轨道作用、电子在片段轨道间转移的分解====== | + | ======EDA-NOCV化学键分析(分子间作用):二聚体能量分解(EDA)中轨道作用、电子在片段轨道间转移的分解====== |
=====理论参考===== | =====理论参考===== | ||
行 12: | 行 12: | ||
=====参数设置===== | =====参数设置===== | ||
+ | **总的来说,对待片段,实际上也就是对片段进行的一个单点计算,因此电荷、自旋极化、收敛情况,都是同样需要注意的。** | ||
将两个分子分为两个Region(如何分区,参考费米维基:[[adf: | 将两个分子分为两个Region(如何分区,参考费米维基:[[adf: | ||
行 31: | 行 32: | ||
因为NOCV的计算,对称性被系统自动设置为NoSymm,分区的计算参数,保存任务时,会自动沿用主任务的参数。但实际上分区如果有对称性,最好是单独将各个分区的任务,对称性设置修改为Auto,这样生成的分区轨道也是按不可约表示来排列,经常会更方便我们的分析。 | 因为NOCV的计算,对称性被系统自动设置为NoSymm,分区的计算参数,保存任务时,会自动沿用主任务的参数。但实际上分区如果有对称性,最好是单独将各个分区的任务,对称性设置修改为Auto,这样生成的分区轨道也是按不可约表示来排列,经常会更方便我们的分析。 | ||
=====结果查看===== | =====结果查看===== | ||
- | ====数据结果==== | + | ====能量分解EDA==== |
在Out窗口,Properties > Bonding Energy Decomposition,可以看到: | 在Out窗口,Properties > Bonding Energy Decomposition,可以看到: | ||
<code bash> | <code bash> | ||
- | Total Bonding Energy: | + | Total Bonding Energy: |
</ | </ | ||
也就是两个氢键总的键能为: | 也就是两个氢键总的键能为: | ||
<code bash> | <code bash> | ||
- | -0.032273260024145 | + | -0.032081195090827 |
</ | </ | ||
- | 也可以看到里面列出了其中轨道相互作用能: | + | 也可以看到里面列出了其中**轨道相互作用能**: |
<code bash> | <code bash> | ||
- | Total Orbital Interactions: | + | Total Orbital Interactions: |
+ | </ | ||
+ | |||
+ | 泡利排斥能: | ||
+ | < | ||
+ | Total Pauli Repulsion: | ||
+ | </ | ||
+ | 静电作用能: | ||
+ | < | ||
+ | Electrostatic Interaction: | ||
</ | </ | ||
那么两个氢键,每一组的轨道相互作用能分别是多少呢?轨道相互作用分别是谁贡献出来的呢? | 那么两个氢键,每一组的轨道相互作用能分别是多少呢?轨道相互作用分别是谁贡献出来的呢? | ||
+ | ====轨道相互作用按共价作用拆分-NOCV分析==== | ||
我们看ETS-NOCV的结果。在Out窗口 > Properties > ETS-NOCV可以看到: | 我们看ETS-NOCV的结果。在Out窗口 > Properties > ETS-NOCV可以看到: | ||
<code bash> | <code bash> | ||
- | 1. Eigenvalue pairs from diagonalization of DeltaP expressed in Lowdin: | ||
- | 1 -0.27662 | ||
- | 2 -0.16882 | ||
- | 3 -0.09052 | ||
- | 4 -0.08262 | ||
- | 5 -0.05873 | ||
- | 2. Orbital Interaction Energy Contributions from each NOCV pair (in a.u.) | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | | ||
3. Orbital Interaction Energy Contributions from each NOCV pair (in kcal/mol) | 3. Orbital Interaction Energy Contributions from each NOCV pair (in kcal/mol) | ||
- | 1 -12.54703 | + | 1 -12.54553 |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | Total sum: -22.2200497016945 | + | Total sum: -22.0624598644894 |
</ | </ | ||
- | 轨道相互作用贡献(Orbital Interaction Energy Contributions)列表里面,贡献最大的两组轨道相互作用,一个是-0.01999 a.u. = -12.54703 kcal/ | + | 轨道相互作用贡献(Orbital Interaction Energy Contributions)列表里面,贡献最大的两组轨道相互作用,一个是-12.54553 kcal/ |
- | 这两组轨道相互作用,NOCV orbital的本征值为±0.27662、±0.16882。 | + | 这两组轨道相互作用,NOCV orbital的本征值分别为±0.27627、±0.16837: |
- | + | < | |
- | 这些数字在上面都有。 | + | 1. Eigenvalue pairs from diagonalization of DeltaP expressed in Lowdin: |
- | + | | |
- | 还有类似如下数据结果: | + | |
+ | 3 -0.09017 | ||
+ | 4 -0.08319 | ||
+ | 5 -0.05765 | ||
+ | </ | ||
+ | ====每组共价作用的电子转移情况==== | ||
+ | 每组轨道作用,其电子转移情况。例如第一组: | ||
<code bash> | <code bash> | ||
- | SFO decomposition of Delta rho k (major contributions): | + | SFO decomposition of Delta rho k (major contributions): |
| | ||
行 88: | 行 92: | ||
| | ||
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
</ | </ | ||
<color blue> | <color blue> | ||
- | * 这是两组较强的轨道相互作用中的第一组,另一组在该数据的下方 | + | * 另二组在该数据的下方 |
- | * 19、18号SFO分别失去0.05453、0.03629个电子,220、224、221 SFO得到0.03846、0.01358、0.01191个电子 | + | * 19、18号SFO分别失去0.05254、0.03552个电子,315、319、316 SFO得到0.03645、0.01359、0.01103个电子 |
- | * <color blue>SFO(symmetrized Fragment Orbital)即碎片轨道。用户可以在SCM > View > Add Isosurface: | + | * <color blue>SFO,symmetrized Fragment Orbital,即碎片轨道。用户可以在SCM > View > Add Isosurface: |
- | ====图形化显示==== | + | < |
- | 这两组轨道相互作用,具体是哪个键?我们可以在View里面去查看。 | + | |
- | ====第一组:本征值为±0.27662==== | + | indx incl.CFs) |
- | ===成键前的样子=== | + | |
- | 点击ADF LOGO > View > Add > Isosurface: With Phase在窗口下方点击Select Field > NOCV Orbitals,选择-0.2766这一个轨道,本征值为负的这个,表示轨道相互作用前,二者的样貌主要是关于N-H...N氢键: | + | |
- | + | ( -26.083 eV) | |
- | {{ :adf:nocv-hbonds05.png?600 |}} | + | |
- | + | | |
- | 并将轨道设置为透明显示(参考费米维基:[[adf: | + | 3 22 2.000 |
- | ===成键后的样子=== | + | |
- | 选择NOCV Orbital | + | 4 23 2.000 -0.819 au A |
- | + | ( -22.276 eV) | |
- | {{ :adf:nocv-hbonds06.png?600 |}} | + | |
- | + | ( | |
- | ===成键导致的电子转移=== | + | |
- | + | ( | |
- | 窗口底部的Select Field区域选择NOCV def Density,并选择本征值为-0.2766这个,即显示第一组轨道相互作用,导致的电子转移: | + | |
- | + | ( | |
- | {{ :adf:nocv-hbonds07.png?600 |}} | + | |
- | + | ( | |
- | 其中红色的区域,表示失去了电子,蓝色区域表示得到了电子,因为电子得失非常微弱,所以等值面的值需要改小才能看得见,这里改成了0.003。这个图,是文献中使用最多的图。其中红色的区域,代表这个区域失去了电子,蓝色区域值,代表这个区域得到了电子。N与N-H之间,出现了电子的富集,代表着氢键的形成。 | + | |
- | + | ( | |
- | ====第二组:本征值为±0.16882==== | + | 10 29 2.000 -0.535 au A |
- | ===成键前的样子=== | + | ( |
- | Select Field > NOCV Orbitals,选择-0.16882这一个轨道,本征值为负的这个,表示轨道相互作用前,二者的样貌主要是关于O...H-N氢键: | + | 11 30 2.000 -0.508 au A |
- | + | ( | |
- | {{ : | + | |
- | + | | |
- | 可以看到是左边碎片O原子的孤对电子(占据轨道)与右边碎片分子H原子的空轨道在发生相互作用,作用能就是前面的-5.29874 kcal/ | + | |
- | + | | |
- | ===成键后的样子=== | + | |
- | + | | |
- | 选择NOCV Orbital 0.16882:就是左边碎片分子的占据轨道和右边碎片分子的空轨道形成的氢键的样子: | + | |
- | + | | |
- | {{ : | + | |
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | ( -6.702 eV) | ||
+ | 19 38 2.000 -0.214 au A | ||
+ | ( -5.837 eV) | ||
+ | 20 39 | ||
+ | </ | ||
+ | 例如19号SFO就是来自分区A,在A分区的能级排序中,能级编号为19 AA,能量为-0.214 au = -5.837 eV。 | ||
+ | ====轨道作用导致的电子转移的图形化显示==== | ||
+ | 我们可以在View里面去查看。SCM > View > Add > Isosurface: With Phase在窗口下方点击Select Field > NOCV def Density,并选择本征值为-0.2763这个,即显示第一组轨道相互作用,导致的电子转移: | ||
- | ===成键导致的电子转移=== | + | {{ : |
- | 窗口底部的Select Field区域选择NOCV def Density,并选择本征值为-0.15558这个,即显示第一组轨道相互作用,导致的电子转移: | + | |
- | {{ : | + | 其中红色的区域,表示失去了电子,蓝色区域表示得到了电子,因为电子得失非常微弱,所以等值面的值需要改小才能看得见,这里改成了0.003。这个图,是文献中使用最多的图。N与N-H之间,出现了电子的富集,N-H上的N也得到了电子,而另一方的N原子失去了电子,这些都是典型的氢键特征。 |
- | 其中红色的区域,代表这个区域失去了电子,蓝色区域值,代表这个区域得到了电子。N-H与O之间,出现了电子的富集,代表着氢键的形成。 | + | 类似可以看到第二组轨道作用,是主要集中在第二个氢键上。 |
=====注意一:===== | =====注意一:===== | ||
从上面可以看到: | 从上面可以看到: | ||
- | - 这种分法和两个氢键分配并不是严格对应的,比如第一组里面实际上包含了少量其他无关成分,第二组也包含了少量的无关成分,但定性的大小关系上可以给我们明确而有力的参考。 | + | - 两组轨道作用,分别主要对应两个氢键,但并没有严格区分开,比如第一组里面实际上包含了少量其他无关成分,第二组也包含了少量的无关成分,但定性上可以给我们明确而有力的参考。 |
- | - 还有其他对轨道作用贡献较小的NOCV,我们可以自己去看看,会发现剩下的NOCV主要是分子内部的极化,也就是分子内部电子的转移。具体可以从NOCV def Density查看成键前后,电子的从什么轨道上转移到什么轨道上。 | + | - 还有其他对轨道作用贡献较小的NOCV,我们可以自己去看看,会发现剩下的NOCV还包括分子内部的极化,也就是分子内部电子的转移。 |
- | - NOCV def Densities,是“正本征值NOCV orbital的模方” - “负本征值NOCV orbital的模方”,不过在理解上要注意,负本征值NOCV Orbital实际上一部分是占据电子的,另一部分是不占据电子的。正本征值NOCV Orbital则是形成化学键,整个轨道上都占了电子的。当然程序自动求得NOCV def Density,已经考虑到了这一点。因此NOCV def Densities表示电子的流向,负值表示失去电子,正值表示得到电子,与“SFO contribution”列表中的贡献值(也就是具体每个SFO得失电子的个数)对应。真正能够可靠地、直观地、明确无误地表征成键作用的,是NOCV def Density。这是一个描述成键导致电子转移的分布图。所以文献中,一般并不讨论NOCV Orbital,而只讨论NOCV def Density。 | + | |
=====注意二:===== | =====注意二:===== | ||
出现“Not all NOCV orbitals pair! The results might be useless”这样的提示,或者NOCV轨道正负本征值没有配对,则表示NOCV方法对这个体系分析失败了,结果无效。 | 出现“Not all NOCV orbitals pair! The results might be useless”这样的提示,或者NOCV轨道正负本征值没有配对,则表示NOCV方法对这个体系分析失败了,结果无效。 | ||
- | =====注意三:===== | ||
- | AMS2018以后,在SCM LOGO > View中,可以显示SFO的空间分布情况。因此,结合该组NOCV的“SFO contribution”列表中负“contribution”的情况,可以具体看到电子从哪个SFO流到哪个SFO去了。失去电子的SFO为Donor orbital,得到电子的为Acceptor orbital。 |