目录

Unrestricted二重态片段的能量分解(EDA)与ETS-NOCV分析:CH3-CH3

前言

本例以两个甲基形成乙烷为例。

参数设置

将乙烷分为两个区(如何分区,参考如何创建分区):

设置基组、泛函等,并勾选Unrestricted、周期性设置为None(在ADF的官方GUITutorial中,使用晶格常数很大的一维周期性Chain,那样也是可以的,但本文中设置为None的方式更好)

在Properties > PEDA-NOCV中打勾,表示要进行NOCV的计算:

在Multilevel > Fragment中勾选Use Fragments,表示要使用片段分析功能,并分别设置每个片段的自旋极化(未配对电子个数):

中Details > Numerical Quality中设置K点(对于NOCV的计算,只需要计算Gamma点即可):

保存任务。保存的时候,会生成两个片段的计算任务,生成片段的计算任务的时候,会提示:

表示软件发现我们设置了片段的自旋极化不为0,那么对片段的计算,也自动打开了Unrestriced。也就是说:对该计算,整体计算用了Unrestriced、片段的计算也用的是Unrestricted。

提交作业。

结果查看

键能能量分解(EDA)

点击ADF LOGO > Output > PEDA Energy Terms:

   P E D A   E n e r g y   T e r m s
  -------------------------------------------------------------------
  PBE                 au           eV     kcal/mol       kJ/mol
  -------------------------------------------------------------------
   E_int    |   -0.19052       -5.184      -119.56       -500.2
   E_Pauli  |    0.25880        7.042       162.40        679.5
   E_elstat |   -0.18875       -5.136      -118.44       -495.6
   E_orb    |   -0.26058       -7.091      -163.52       -684.2
            |
   E^0      |    0.07006        1.906        43.96        183.9
   T^0      |    0.67563       18.385       423.97       1773.9
   Elst^0   |   -0.43940      -11.957      -275.73      -1153.6
   XC^0     |   -0.16618       -4.522      -104.28       -436.3
            |
   E_orb    |   -0.26058       -7.091      -163.52       -684.2
   T_orb    |   -0.54525      -14.837      -342.15      -1431.6
   Elst_orb |    0.24049        6.544       150.91        631.4
   XC_orb   |    0.04418        1.202        27.73        116.0
  -------------------------------------------------------------------

表示EDA分析中,两个甲基之间的轨道相互作用能是-0.26058 au = -7.091 eV = -163.52 kcal/mol = -684.2 KJ/mol 其他各项能量的含义,具体参考pEDA能量各项的含义

NOCV分析:查看“轨道相互作用能”主要是什么轨道贡献出来的

点击ADF LOGO > Output > PEDA-NOCV 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
 
 
    ===K-Point  1   Spin 1
 
  NOCV      eigenvalue         E_tot         T         V
  -------------------------------------------------------------------
    1  |  -0.4736  +0.4736    -317.3    -496.2    +178.9
    2  |  -0.0585  +0.0585      -7.4     +65.6     -73.0
    3  |  -0.0548  +0.0548      -6.8     +56.1     -62.9
    4  |  -0.0496  +0.0496      -3.0    -154.5    +151.5
    5  |  -0.0456  +0.0456      -2.5    -147.6    +145.2
    6  |  -0.0417  +0.0417      -4.0     +14.9     -18.9
                     E_res=     -1.1 kJ/mol
                     E_tot=   -342.0 kJ/mol
 
 
    ===K-Point  1   Spin 2
 
  NOCV      eigenvalue         E_tot         T         V
  -------------------------------------------------------------------
    1  |  -0.4736  +0.4736    -317.3    -496.2    +178.9
    2  |  -0.0585  +0.0585      -7.4     +65.6     -73.0
    3  |  -0.0548  +0.0548      -6.8     +56.1     -62.9
    4  |  -0.0496  +0.0496      -3.0    -154.5    +151.5
    5  |  -0.0456  +0.0456      -2.5    -147.6    +145.2
    6  |  -0.0417  +0.0417      -4.0     +14.9     -18.9
                     E_res=     -1.1 kJ/mol
                     E_tot=   -342.0 kJ/mol
 
 =====================================================================

表示EDA分析中,两个甲基之间的轨道相互作用能-684.2 KJ/mol,主要是由第一组NOCV造成的(总共列出了6组)。其中第一组的alpha自旋这部分(也就是上表中的Spin 1),贡献了-317.3KJ/Mol;beta自旋这部分(也就是上表中的Spin 2),也贡献了-317.3KJ/Mol;第一组的本征值为-0.4736和+0.4736。其中本征值的正负的含义,参考:ETS-NOCV理论。下文中也略作说明。

NOCV分析:查看对“轨道相互作用能”做出主要贡献的NOCV

ADF LOGO > View > Fields → Grid → Fine, Add → Isosurface: With Phase:

分别选择Spin为A(A表示alpha电子、B表示beta电子),本征值分别为-0.4736、+0.4736的两个NOCV,显示如下:

其中,将NOCV设置为透明的显示,可以参考:如何设置View中各种空间分布图的透明度

其中Spin为A、本征值为-0.4736的这个,列出的是左边片段、右边片段,其中一个的占据轨道和另一个的空轨道。成键的时候,是这两个轨道在产生“轨道相互作用”(ADF模块使用Unrestricted计算<chem>CH3</chem>,可以发现AA不可约表示的4号轨道Spin A是占据的、Spin B是空的,其中空的B这个轨道比占据的A轨道略“胖”。分别查看两个片段的电子占据方式,也可以得到这个信息)。那么很显然是左边的甲基<chem>SP3</chem>杂化形成的一个未配对的alpha电子(按CTRL -缩小原子,可以看的更清楚一些,如下图所示)和右边甲基的<chem>SP3</chem>杂化形成的alpha空轨道在发生“轨道相互作用”;Spin B的情况类似,不过是反过来的。

其中Spin为A、本征值为0.4736的这个,表示两个甲基成键之后,左边甲基的那个占据轨道电子和右边甲基的空轨道发生“轨道相互作用”之后,变成了什么形状。Spin B也可以类似看到,和Spin A的朝向相反。

在View窗口可以分别显示所有的轨道的空间分布图。