软件版本为AMS2019.3及之前的版本,AMS2020以后的版本请参考链接:电荷分解分析(CDA)
理论上,采用的是Mulliken布居方法。
ADF的任何一个片段分析计算,均默认给出CDA情况。片段分析的参数设置,参考:
点击ADF LOGO > Output > Properties > SFO Populations,即列出片段轨道的布居情况,例如:
=== AA === SFO contributions (%) per orbital (multiplication by the orbital occupation yields the SFO Gross Populations) Orb.: 24 25 26 27 28 29 30 31 32 33 34 35 36 37 occup: 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 0.00 0.00 0.00 0.00 CF+SFO ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ------ 12: 0.00 0.00 0.00 2.64 2.64 0.00 3.12 3.16 -0.00 0.00 0.00 -0.00 -0.11 -0.11 18: 0.00 0.00 0.00 2.64 2.64 0.00 3.16 3.12 0.00 0.00 0.00 -0.00 -0.11 -0.11 省略……………………… 137: 0.00 19.43 0.00 29.82 0.00 24.64 19.17 0.00 5.35 0.00 0.00 0.20 0.02 0.00 138: 0.00 0.99 0.00 0.69 0.00 0.09 0.65 0.00 1.04 0.00 0.00 22.13 46.41 0.00 Summation over all MOs, multiplied by occupation: Total SFO Gross Populations in this Irrep =========================================================================================== 2.00 2.00 2.00 0.01 0.03 0.00 0.00 -0.00 0.00 2.00 2.00 0.33 0.01 0.00 0.00 2.00 2.00 0.33 0.01 0.00 0.00 1.92 0.02 -0.00 1.84 0.01 0.00 1.08 -0.03 -0.00 0.00 0.00 0.01 0.01 2.00 2.00 1.99 1.99 1.51 0.02 0.03 0.01 0.01 -0.00 0.00 0.00 -0.00 -0.00 0.00 0.00 -0.00 -0.00 1.98 0.08 -0.01 -0.00 0.00 -0.00 0.00 0.00 省略…………………… 0.00 0.00 -0.00 -0.00 1.98 0.08 -0.01 -0.00 0.00 -0.00 0.00 0.00 0.00 0.00
注意:
那么第N个片段轨道分别是什么片段轨道呢?
点击Properties > SFO construction,可以看到所有的片段轨道分别是什么。同样地,也是按照不可约表示来列的,例如AA不可约表示的片段轨道:
Nr. of SFOs : 146 Cartesian basis functions that participate in this irrep (total number = 159) : 1 2 3 4 11 12 40 41 42 5 6 8 9 13 14 16 17 19 20 43 44 27 33 39 22 23 25 28 29 31 34 35 37 51 54 46 47 49 52 56 57 58 59 60 61 62 64 65 67 68 75 70 71 73 76 77 78 79 80 81 82 84 85 87 88 95 90 91 93 96 97 98 99 100 101 102 104 105 107 108 115 110 111 113 116 117 118 119 120 121 122 124 125 127 128 135 130 131 133 136 137 138 139 140 141 142 144 145 147 148 155 150 151 153 156 157 158 159 160 161 162 164 165 167 168 175 170 171 173 176 177 178 179 180 181 182 184 185 187 188 195 190 191 193 196 197 198 199 200 201 202 204 205 207 208 215 210 211 213 SFO (index Fragment Generating Expansion in Fragment Orbitals indx incl.CFs) Occup Orb.Energy FragmentType Coeff. Orbital on Fragment -------------------------------------------------------------------------------------- 1 1 2.000 -299.542 au Ni 1.00 1 S 1 ( -8150.961 eV) 2 2 2.000 -35.461 au Ni 1.00 2 S 1 ( -964.940 eV) 3 3 2.000 -3.742 au Ni 1.00 3 S 1 ( -101.818 eV) 4 4 -- -0.068 au Ni 1.00 4 S 1 ( -1.859 eV) 5 5 -- 0.131 au Ni 1.00 5 S 1 ( 3.557 eV) 6 6 -- 0.882 au Ni 1.00 6 S 1 ( 23.994 eV) 7 7 -- 5.326 au Ni 1.00 7 S 1 省略…………………… ( 70.293 eV) 0.71 7 PI:y 5 144 144 -- 4.431 au CO1 0.71 8 PI:x 5 ( 120.586 eV) 0.71 8 PI:y 5 145 145 -- 0.957 au CO1 1.00 1 DELTA:xy 5 ( 26.046 eV) 146 146 -- 1.584 au CO1 1.00 2 DELTA:xy 5 ( 43.100 eV)
可以看到AA不可约表示有146个片段轨道,编号以 (index incl.CFs)这一列为准。如果不使用冻芯近似,第一列编号和第二列编号是一致的,本例就是这种情况。
所以可以看到前面说的第4个片段轨道(也就是有0.01个电子的那个片段轨道)是:
4 4 -- -0.068 au Ni 1.00 4 S 1 ( -1.859 eV)
也就是Ni这个片段,片段轨道的能级是-0.068 au=-1.859 eV。可以看到这个轨道在片段孤立存在的时候,上面的电子个数(也就是Occup这一列数据)显示为“- -”也就是没有电子。所以形成配合物的时候,这个片段轨道得到了0.01电子。
如果做了冻芯近似,那么第二列数字就不是从1开始了,有多少个电子被冻结,就缺多少个编号。但那些被冻结的电子,在前面 Properties > SFO Populations也是列出来了的,只是电子的个数设置为0了。但如果记住Properties > SFO construction是以第二列数字为准,来数Properties > SFO Populations就不会数错了。