这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版上一修订版两侧同时换到之后的修订记录 | ||
adf:peda-nocv-polymer [2019/12/09 19:53] – liu.jun | adf:peda-nocv-polymer [2022/02/24 20:37] – [ETS-NOCV计算:聚合物吸附的化学键分析] liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ======ETS-NOCV计算:聚合物的片段相互作用、电子转移====== | + | ======ETS-NOCV计算:聚合物吸附、成键的化学键分析====== |
=====预备知识===== | =====预备知识===== | ||
**由于pEDA分析支持正常k空间布点,而NOCV则只支持Gamma点,因此建议二者分别计算。本文使用AMS2019.301完成。** | **由于pEDA分析支持正常k空间布点,而NOCV则只支持Gamma点,因此建议二者分别计算。本文使用AMS2019.301完成。** | ||
- | NOCV orbital、NOCV density、NOCV Def(ference) density的化学直观含义,参考:[[adf: | + | NOCV orbital、NOCV density、NOCV Def(ference) density的化学直观含义,参考:[[adf: |
=====模型===== | =====模型===== | ||
行 52: | 行 52: | ||
设置完成后,File > Save As,保存任务。 | 设置完成后,File > Save As,保存任务。 | ||
+ | =====提交任务===== | ||
- | **3,结果查看:** | + | 提交任务的方式,参考[[adf: |
- | **1)查看NOCV能量项:** | + | ====查看NOCV轨道==== |
- | 在ADFinput窗口点击SCM LOGO > Output,点击Properties > PEDA-NOCV Energy Terms | + | 因为NOCV本身是为了进一步分析EDA的结果,希望看到EDA中轨道相互作用(也就是共价作用)主要由哪些轨道贡献出来,因此我们首先看看有多少贡献比较大的NOCV。 |
- | 各项意义,参考预备知识。 | + | SCM - Output - Properties - PEDA Energy Terms: |
+ | < | ||
+ | P E D A N O C V E n e r g y T e r m s | ||
+ | ------------------------------------------------------------------- | ||
+ | Energies in kJ/mol | ||
- | **2)查看NOCV轨道:** | ||
- | 在ADFinput或者Output窗口点击: SCM LOGO > View > Fields → Grid → Fine, Add → Isosurface: With Phase | + | ===K-Point |
- | {{: | + | NOCV eigenvalue |
+ | ------------------------------------------------------------------- | ||
+ | 1 | -0.1419 +0.1419 | ||
+ | 2 | ||
+ | 3 | -0.0472 | ||
+ | 4 | -0.0284 | ||
+ | 5 | -0.0249 | ||
+ | 6 | -0.0218 | ||
+ | 7 | -0.0176 | ||
+ | 8 | -0.0144 | ||
+ | 9 | -0.0126 | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | </ | ||
+ | 可以看到贡献主要是第一组本征值为±0.1419的NOCV贡献出来(注意,本征值必须严格±成对,否则NOCV数据不可用),我们需要记住这个本征值,然后来查看该NOCV的具体情况。 | ||
+ | SCM → View → Add → Isosurface: With Phase → Select Field → NOCV Orbitals。 | ||
- | {{:adf:polymernocv07.png|}} | + | 注意: |
+ | 负本征值NOCV orbital则是成键前的样子,也就是一方碎片贡献占据轨道,另一方贡献空轨道(等值面调整到0.01看的更清楚一些,透明设置参考:[[adf:transparrency]]): | ||
- | **3)查看某NOCV引起的形变密度(电子从红色区域流向蓝色区域):** | + | {{ : |
- | {{: | + | 可以看到,贡献占据轨道的是聚合物分子的O原子的孤对电子(占据轨道),与水分子的空轨道(H原子上)作用。 |
- | {{: | + | 正本征值的NOCV orbital是成键之后的键轨道: |
- | 设置NOCV轨道和形变密度透明度、背景色、导出图片,参考[[adf:pEDA-NOCV-surface]] | + | {{ :adf:1d-peda09.png? |
- | + | ||
- | 注意: | + | |
- | * SCM LOGO > View> | + | |
- | * SCM LOGO > View> | + | |
- | * SCM LOGO > View> | + | |
- | * 所有的NOCV Def Densities加起来,就是形成吸附的过程中总的电子转移 | + | |
- | =====电子的转移===== | + | 形成了氢键。 |
+ | ====查看该组成键引起的密度变化的贡献==== | ||
- | NOCV是成对出现的(alpha、beta),一般研究第i个NOCV带来的电子转移,要把第i个alpha的NOCV的NOCV Def Densities和第i个beta的NOCV的NOCV Def Densities加起来:SCM LOGO > View, Fields > Grid > Fine, Fields > calculate,之后窗口下方 | + | SCM → View → Add → Isosurface: With Phase |
+ | {{ : | ||
- | {{: | + | 电子从红色区域流向蓝色区域。本例中,这是一个典型的氢键电子转移,从聚合物O原子转移电子到H...O之间,H原子本身也失去电子,与H直接连接的O原子也得到了电子。电子转移很微弱,因此等值面调整到0.0003看起来才比较清晰。 |
+ | ====导出图片==== | ||
- | - 第一个是编号,此处为C-1,后面会用到 | + | 参考:[[adf: |
- | - 第二个是选择要加的第1个量,选择例如NOCV Def Densities列表中第1个 | + | |
- | - 第三个是选择符号,默认为-,这里改为+ | + | |
- | - 第四个是要加的第二个量,选择例如NOCV Def Densities列表中第2个(通常第2n个和第2n+1个(n> | + | |
- | 之后,Add> | + | 延伸阅读:[[adf:theoryofets-nocv]]、[[adf: |
- | 根据“NOCV对”(包含alpha和beta)对ΔE< | ||
- | 如果将所有“NOCV对”的Def Densities加和,则得到形成该键的总的电荷转移。但一般而言,实际上主要贡献的“NOCV对”只有1~2对。因此一般看前1对,2对的Def Densities加和就能看出规律来了。 |