这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
adf:nocv-disdinguishbonding [2019/12/07 12:47] – liu.jun | adf:nocv-disdinguishbonding [2020/11/22 10:58] (当前版本) – 移除 liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ======ETS-NOCV功能案例:二聚体能量分解(EDA)中轨道作用、电子在片段轨道间转移的分解====== | ||
- | =====理论参考===== | ||
- | * [[adf:eda]] | ||
- | * [[theoryofets-nocv]] | ||
- | |||
- | 该功能不局限于二聚体。实际上N(N> | ||
- | |||
- | {{ : | ||
- | |||
- | 两个分子形成了两个氢键,那么两个氢键,哪个氢键的轨道作用强呢?各自多大呢?可以通过ETS-NOCV得到定量的数值。 | ||
- | |||
- | =====参数设置===== | ||
- | |||
- | 将两个分子分为两个Region(如何分区,参考费米维基:[[adf: | ||
- | |||
- | {{ : | ||
- | |||
- | 并设置参数: | ||
- | |||
- | {{ : | ||
- | |||
- | ADF2017去掉了上图中Fragment Analysis这个选项,实际上这个选项的作用仅仅是自动勾选ADFinput > Multilevel > Fragment > Use fragment而已,因此此处勾选Single Point,手动勾选ADFinput > Multilevel > Fragment > Use fragment也一样。 | ||
- | |||
- | 其中NOCV的设置,选择Close-shell,因为这里两个分子都是闭壳层,alpha电子和beta电子的行为是一样的。如果其中有一个片段是开壳层,就一定要选择open-shell了,将分别分析α电子、β电子的轨道相互作用。 | ||
- | |||
- | {{ : | ||
- | |||
- | File > Save as保存任务,保存的时候,会生成三个任务,其中两个是单独计算两个片段孤立存在于真空的情况,第三个是计算两个片段放置在一起的情况。本例下载的文件中,任务名分别为02NOCV.A、02NOCV.T和02NOCV。在Jobs窗口选中02NOCV一行,点击Job > Run,则自动依次计算三个任务。或者在02NOCV的ADFinput窗口直接File > run也可以。如果是提交到集群计算,就要参考费米维基:[[adf: | ||
- | =====结果查看===== | ||
- | ====数据结果==== | ||
- | 在Out窗口,Properties > Bonding Energy Decomposition,可以看到: | ||
- | <code bash> | ||
- | Total Bonding Energy: | ||
- | </ | ||
- | |||
- | 也就是两个氢键总的键能为: | ||
- | <code bash> | ||
- | -0.032273260024145 Hatree = -0.8782 eV = -20.25 kcal/mol = -84.73 kJ/mol | ||
- | </ | ||
- | |||
- | 也可以看到里面列出了其中轨道相互作用能: | ||
- | <code bash> | ||
- | Total Orbital Interactions: | ||
- | </ | ||
- | |||
- | 那么两个氢键,每一组的轨道相互作用能分别是多少呢?轨道相互作用分别是谁贡献出来的呢? | ||
- | |||
- | 我们看ETS-NOCV的结果。在Out窗口 > Properties > ETS-NOCV可以看到: | ||
- | <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) | ||
- | 1 -12.54703 | ||
- | | ||
- | | ||
- | | ||
- | | ||
- | Total sum: -22.2200497016945 | ||
- | </ | ||
- | |||
- | |||
- | 轨道相互作用贡献(Orbital Interaction Energy Contributions)列表里面,贡献最大的两组轨道相互作用,一个是-0.01999 a.u. = -12.54703 kcal/ | ||
- | |||
- | 这两组轨道相互作用,NOCV orbital的本征值为±0.27662、±0.16882。 | ||
- | |||
- | 这些数字在上面都有。 | ||
- | |||
- | 还有类似如下数据结果: | ||
- | |||
- | <code bash> | ||
- | SFO decomposition of Delta rho k (major contributions): | ||
- | |||
- | | ||
- | | ||
- | | ||
- | |||
- | | ||
- | | ||
- | | ||
- | 220 SFO contribution: | ||
- | | ||
- | 224 SFO contribution: | ||
- | 221 SFO contribution: | ||
- | </ | ||
- | <color blue> | ||
- | |||
- | * 这是两组较强的轨道相互作用中的第一组,另一组在该数据的下方 | ||
- | * 19、18号SFO分别失去0.05453、0.03629个电子,220、224、221 SFO得到0.03846、0.01358、0.01191个电子 | ||
- | * <color blue> | ||
- | ====图形化显示==== | ||
- | 这两组轨道相互作用,具体是哪个键?我们可以在View里面去查看。 | ||
- | ====第一组:本征值为±0.27662==== | ||
- | ===成键前的样子=== | ||
- | 点击ADF LOGO > View > Add > Isosurface: With Phase在窗口下方点击Select Field > NOCV Orbitals,选择-0.2766这一个轨道,本征值为负的这个,表示轨道相互作用前,二者的样貌主要是关于N-H...N氢键: | ||
- | |||
- | {{ : | ||
- | |||
- | 并将轨道设置为透明显示(参考费米维基:[[adf: | ||
- | ===成键后的样子=== | ||
- | 选择NOCV Orbital 0.2766:就是右边碎片分子的占据轨道和左边碎片分子的空轨道形成的氢键的样子: | ||
- | |||
- | {{ : | ||
- | |||
- | ===成键导致的电子转移=== | ||
- | |||
- | 窗口底部的Select Field区域选择NOCV def Density,并选择本征值为-0.2766这个,即显示第一组轨道相互作用,导致的电子转移: | ||
- | |||
- | {{ : | ||
- | |||
- | 其中红色的区域,表示失去了电子,蓝色区域表示得到了电子,因为电子得失非常微弱,所以等值面的值需要改小才能看得见,这里改成了0.003。这个图,是文献中使用最多的图。其中红色的区域,代表这个区域失去了电子,蓝色区域值,代表这个区域得到了电子。N与N-H之间,出现了电子的富集,代表着氢键的形成。 | ||
- | |||
- | ====第二组:本征值为±0.16882==== | ||
- | ===成键前的样子=== | ||
- | Select Field > NOCV Orbitals,选择-0.16882这一个轨道,本征值为负的这个,表示轨道相互作用前,二者的样貌主要是关于O...H-N氢键: | ||
- | |||
- | {{ : | ||
- | |||
- | 可以看到是左边碎片O原子的孤对电子(占据轨道)与右边碎片分子H原子的空轨道在发生相互作用,作用能就是前面的-5.29874 kcal/mol。 | ||
- | |||
- | ===成键后的样子=== | ||
- | |||
- | {{ : | ||
- | |||
- | 选择NOCV Orbital +0.15558:就是左边碎片分子的占据轨道和右边碎片分子的空轨道形成的氢键的样子: | ||
- | |||
- | {{ : | ||
- | |||
- | ===成键导致的电子转移=== | ||
- | 窗口底部的Select Field区域选择NOCV def Density,并选择本征值为-0.15558这个,即显示第一组轨道相互作用,导致的电子转移: | ||
- | |||
- | {{ : | ||
- | |||
- | 其中红色的区域,代表这个区域失去了电子,蓝色区域值,代表这个区域得到了电子。N-H与O之间,出现了电子的富集,代表着氢键的形成。 | ||
- | =====注意一:===== | ||
- | 从上面可以看到: | ||
- | - 这种分法并不是严格的,比如第一组里面实际上包含了少量其他无关成分,第二组也包含了少量的无关成分,但定性的大小关系上可以给我们明确而有力的参考。 | ||
- | - 还有其他对轨道作用贡献较小的NOCV,我们可以自己去看看,会发现剩下的NOCV主要是分子内部的极化,也就是分子内部电子的转移。具体可以从NOCV def Density查看成键前后,电子的从什么轨道上转移到什么轨道上。 | ||
- | - NOCV def Densities,是“负本征值NOCV orbital的模方” - “正本征值NOCV orbital的模方”,不过在理解上要注意,负本征值NOCV Orbital实际上一部分是占据电子的,另一部分是不占据电子的。正本征值NOCV Orbital则是形成化学键,整个轨道上都占了电子的。当然程序自动求得NOCV def Density,已经考虑到了这一点。因此表示电子的流向,负值表示失去电子,正值表示得到电子,与“SFO contribution”列表中的贡献值(也就是具体每个SFO得失电子的个数)对应。真正能够可靠地、直观地、明确无误地表征成键作用的,是NOCV def Density。这是一个描述成键导致电子转移的分布图。所以文献中,一般并不讨论NOCV Orbital,而只讨论NOCV def Density。 | ||
- | =====注意二:===== | ||
- | 出现“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。 |