两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 |
adf:fragmentanalysis [2019/12/07 01:36] – [片段轨道布居] liu.jun | adf:fragmentanalysis [2020/11/21 10:48] (当前版本) – 移除 liu.jun |
---|
======分子间相互作用能、分子之间相互作用能、键能、键解离能、结合能计算、键能分解EDA、片段轨道布居====== | |
使用该功能发表的代表性工作: | |
* [[http://www.fermitech.com.cn/adf/highlight-006|N-卤-鸟嘌呤四股结构的共振增强卤键]] | |
* [[http://www.fermitech.com.cn/adf/highlight-040|应变促进叠氮-炔环加成反应的高度加速(Nature Comm.,2014)]] | |
* [[http://www.fermitech.com.cn/adf/highlight-019|染料敏化太阳能电池中钌染料的效率(2014)]] | |
* [[http://www.fermitech.com.cn/adf/highlight-048|元结异构体(Agostic Isomers)-不同的颜色和结构(Angew. Chem. Int. Ed.,2013)]] | |
* [[http://www.fermitech.com.cn/adf/highlight-044|f区元素成键相互作用:新奇的铀化合物(Nature Chem.,2013;Angew. Chem. Int. Ed.,2011;JACS,2013)]] | |
* [[http://www.fermitech.com.cn/adf/highlight-051|十二个单电子配体的钼配合物(Angew. Chem. Int. Ed.,2008)]] | |
| |
氢键的片段分析、电子转移、电荷转移的分析,参考:[[adf:hbondanalysis|]]。关于键能分解的理论,参考:[[adf:eda]] | |
| |
片段是一个比较广义的概念,可以是分子的部分和部分之间相互作用,也可以是分子和分子之间的相互作用,方法原理都是一样。如果是分子内部一个区域与另一个区域相互作用,那么将这两个区域当成两个片段即可,如果是分子和分子间相互作用,那么将每个分子当成一个片段就可以了。如果要考虑相对论的片段分析,那么设置参数的时候,勾选 Main > Relativity > Scalar就可以了,目前ADF不支持Spin-Orbit的片段功能。 | |
| |
下面以分子内部区域与区域之间的相互作用为例子,来演示计算的过程。 | |
| |
1,创建分子(具体创建过程,可以参考:[[adf:1000atomsball]]、[[https://www.jianguoyun.com/p/Dfq5zjUQmZ2ZBhjprSc|建模:ADF模块分子的基本建模功能演示(视频)]]): | |
| |
首先应该进行结构优化,之后才进行片段分析。结构优化如何使用?请参考:[[adf:geoopt|]]、[[adf:geooptforbigsystem]] | |
| |
假定已经优化完成,如下演示的是片段分析的参数设置: | |
| |
{{ :adf:fragment01.png?650 }} | |
| |
关于参数设置详细介绍,请参考:[[adf:parameters]] | |
| |
{{ :adf:fragment02.png?650 }} | |
| |
分区操作,参考:[[adf:creatregions]]。手动勾选Multilevel > Fragment > Use fragment。 | |
| |
{{:adf:fragment03.png|}} | |
| |
| |
保存任务,将生成总共3个任务,其中一个是主任务,另外两个任务是两个碎片的任务。提交并行任务,尤其是这种碎片分析相关的任务,参考:[[adf:maintance|]] | |
| |
计算完成,查看结果(如何查看,参考:[[adf:showresults|]]): | |
====片段分析==== | |
SCM > level | |
| |
{{ :adf:fragment04.png?650 }} | |
| |
第二列轨道是分子轨道(底部有分子的名称,例如本例保存任务名字叫AT,因此底部分子也叫做AT),其他列为片段轨道,底部有碎片的名称。该图定性显示整个二聚体的分子轨道与两个碎片的分子轨道之间的关系。红线越粗,表示该碎片轨道贡献越大。鼠标放置在能级上,将显示定量贡献百分比。 | |
====EDA==== | |
| |
点击SCM的LOGO > Output > Properties > Bonding Energy Decomposition查看文本格式的输出文件中能量分解部分: | |
| |
{{ :adf:fragment06.png?650 }} | |
| |
其中键能,就是Total Bonding Energy。这个数字,在logfile末尾也有。 | |
| |
上图中进行了EDA(Energy Decomposition Analysis)分析: | |
- 理论参考:[[adf:eda]] | |
- 本例使用了色散修正泛函,因此还包含第四个部分:色散能。 | |
| |
====键能、键解离能、结合能计算==== | |
这个如果是两个片段的话,那么这个能量就是两个片段之间的结合能。 | |
| |
键解离能有两种描述方式,一种就是直接采用键能(一般为负值);另一种需要在键能的基础上加上准备能(聚合体中的片段们,分别进行能量最小化,之后的能量降低量之和,即准备能,为正值) | |
| |
====片段轨道布居==== | |
| |
=== AA === | |
| |
| |
SCM LOGO > Output > Properties > %SFO per orbital,会显示类似如下内容: | |
<code> | |
SFO contributions (%) per orbital | |
(multiplication by the orbital occupation yields the SFO Gross Populations) | |
| |
Orb.: 28 29 30 31 32 33 34 35 36 37 38 39 40 41 | |
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 ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- | |
------ | |
30: 0.01 1.34 0.00 0.40 0.02 0.00 0.05 0.13 0.29 0.07 0.00 0.00 0.00 0.01 | |
33: 0.31 3.25 0.01 0.22 0.08 0.00 0.06 0.08 0.06 0.02 0.00 0.00 0.00 0.00 | |
34: 0.00 11.58 0.33 2.55 0.01 0.01 0.00 0.14 0.15 0.13 0.00 0.00 0.00 0.00 | |
35: 5.48 0.91 81.44 9.52 0.65 0.04 0.16 0.02 0.09 0.02 0.00 0.00 0.00 0.00 | |
36: 0.00 0.10 0.03 0.03 0.21 93.18 5.39 0.29 0.17 0.07 0.00 0.00 0.00 0.00 | |
37: 0.01 0.28 0.11 0.25 0.70 2.87 54.47 27.53 4.35 5.78 -0.01 0.00 0.03 -0.01 | |
38: 0.02 0.29 0.07 0.22 0.55 0.61 8.62 0.20 3.74 81.10 -0.01 0.00 0.05 -0.01 | |
39: 0.00 0.07 0.00 0.00 0.00 0.00 0.01 0.06 0.06 0.06 0.24 83.34 0.90 8.57 | |
40: 0.00 0.16 0.00 0.01 0.00 0.00 0.02 0.12 0.13 0.16 0.09 16.04 8.82 53.46 | |
41: -0.01 0.03 0.00 0.00 0.00 0.00 0.02 0.03 0.00 0.01 0.09 0.03 0.11 2.80 | |
43: -0.02 0.16 0.00 0.07 -0.01 0.00 0.00 0.02 0.15 0.10 0.03 0.12 -0.09 10.63 | |
175: 0.03 0.00 0.77 0.22 0.01 0.58 1.21 0.38 0.00 0.62 -0.01 0.00 0.00 0.00 | |
177: 93.46 0.21 4.97 1.09 0.01 0.00 0.01 0.00 0.00 0.02 0.00 0.00 0.00 0.00 | |
178: 0.34 79.99 0.90 0.22 0.52 0.18 0.32 0.00 0.01 0.04 0.00 0.00 0.01 0.01 | |
179: 0.08 0.00 10.67 83.62 0.56 0.12 0.18 0.18 0.06 0.20 0.00 0.00 0.00 0.01 | |
180: 0.02 0.16 0.29 1.03 96.42 0.43 0.56 0.11 0.00 0.18 0.00 0.00 0.00 0.00 | |
181: 0.01 0.02 0.13 0.03 0.00 0.20 19.59 66.40 11.29 1.25 -0.01 0.00 0.01 0.00 | |
182: 0.01 0.07 0.03 -0.01 0.00 0.76 6.94 2.60 78.53 7.86 0.00 0.00 0.01 0.03 | |
183: 0.00 0.00 0.00 0.01 0.00 0.02 0.03 0.05 0.00 0.09 93.75 0.04 0.32 0.13 | |
184: 0.02 0.00 0.00 0.08 0.00 0.18 0.31 0.45 0.01 0.64 4.40 0.03 29.87 15.14 | |
185: 0.01 0.00 -0.02 0.03 0.00 0.05 0.02 0.15 0.01 0.24 0.73 0.29 59.33 0.23 | |
186: 0.00 0.00 0.00 0.01 0.00 0.05 0.10 0.07 0.00 0.11 0.05 0.00 0.32 7.04 | |
Summation over all MOs, multiplied by occupation: Total SFO Gross Populations in this Irrep | |
=========================================================================================== | |
| |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 2.00 1.99 2.00 2.00 1.99 | |
2.00 1.99 2.00 2.00 2.00 1.99 2.00 2.00 2.00 2.00 2.00 2.00 | |
1.97 1.95 0.01 0.02 0.00 0.01 0.02 0.00 0.00 0.00 0.00 0.01 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.99 1.98 2.00 2.00 | |
2.00 2.00 1.99 2.00 2.00 2.00 1.99 1.99 2.00 2.00 1.99 2.00 | |
1.99 1.96 0.01 0.05 0.02 0.01 0.00 0.01 0.00 0.01 0.01 0.00 | |
0.00 0.00 0.00 0.00 0.01 0.00 0.00 0.01 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 | |
| |
</code> | |
| |
注意, | |
* 如果体系有对称性,则该数据是按照不可约表示分别列出的。不可约表示的通俗理解,参考:[[adf:symmetryandrepresentation]] | |
* SFO contributions (%) per orbital显示的是整个二聚体的分子轨道(也就是 Orb.那一行,其数字是轨道的序号, occup是该轨道的电子占据数)与碎片轨道之间的组合关系 | |
* Total SFO Gross Populations in this Irrep显示的是,每个碎片轨道上,分布了多少个电子,这里面可以看到存在微弱的电子转移,例如第24号SFO失去了0.01电子,因此只有1.99电子了。最前面的19个轨道是冻芯轨道,没有列出电子占据数,显示为0。 | |
| |
24SFO是哪个轨道呢?可以查询。点击Properties > SFO construction,看到AA不可约表示中的SFO列表(同样地,需要看对应的不可约表示,例如这里需要看AA): | |
<code bash> | |
| |
=== AA === | |
Nr. of SFOs : 279 | |
Cartesian basis functions that participate in this irrep (total number = 323) : | |
1 20 39 58 77 96 115 134 153 172 | |
221 240 259 278 297 316 335 354 373 2 | |
3 4 5 6 8 9 11 12 19 14 | |
15 17 21 22 23 24 25 27 28 30 | |
31 38 33 34 36 40 41 42 43 44 | |
46 47 49 50 57 52 53 55 59 60 | |
61 62 63 65 66 68 69 76 71 72 | |
74 78 79 80 81 82 84 85 87 88 | |
95 90 91 93 97 98 99 100 101 103 | |
104 106 107 114 109 110 112 116 117 118 | |
119 120 122 123 125 126 133 128 129 131 | |
135 136 137 138 139 141 142 144 145 152 | |
147 148 150 154 155 156 157 158 160 161 | |
163 164 171 166 167 169 173 174 175 176 | |
177 179 180 182 183 190 185 186 188 191 | |
192 193 194 195 197 198 199 200 201 203 | |
204 205 206 207 209 210 211 212 213 215 | |
216 217 218 219 222 223 224 225 226 228 | |
229 231 232 239 234 235 237 241 242 243 | |
244 245 247 248 250 251 258 253 254 256 | |
260 261 262 263 264 266 267 269 270 277 | |
272 273 275 279 280 281 282 283 285 286 | |
288 289 296 291 292 294 298 299 300 301 | |
302 304 305 307 308 315 310 311 313 317 | |
318 319 320 321 323 324 326 327 334 329 | |
330 332 336 337 338 339 340 342 343 345 | |
346 353 348 349 351 355 356 357 358 359 | |
361 362 364 365 372 367 368 370 374 375 | |
376 377 378 380 381 383 384 391 386 387 | |
389 392 393 394 395 396 398 399 400 401 | |
402 404 405 406 407 408 410 422 411 423 | |
412 424 415 427 413 414 425 426 416 417 | |
418 419 420 | |
| |
SFO (index Fragment Generating Expansion in Fragment Orbitals | |
indx incl.CFs) Occup Orb.Energy FragmentType Coeff. Orbital on Fragment | |
-------------------------------------------------------------------------------------- | |
1 20 2.000 -0.960 au Region_1 1.00 1 AA 1 | |
( -26.120 eV) | |
2 21 2.000 -0.902 au Region_1 1.00 2 AA 1 | |
( -24.534 eV) | |
3 22 2.000 -0.865 au Region_1 1.00 3 AA 1 | |
( -23.545 eV) | |
4 23 2.000 -0.820 au Region_1 1.00 4 AA 1 | |
( -22.309 eV) | |
5 24 2.000 -0.786 au Region_1 1.00 5 AA 1 | |
( -21.386 eV) | |
6 25 2.000 -0.702 au Region_1 1.00 6 AA 1 | |
( -19.099 eV) | |
7 26 2.000 -0.613 au Region_1 1.00 7 AA 1 | |
( -16.688 eV) | |
8 27 2.000 -0.596 au Region_1 1.00 8 AA 1 | |
</code> | |
可以看到,24 SFO实际上就是Region_1的5 AA轨道,在Region_1碎片单独存在的时候,能级为-0.786 au = -21.386 eV | |
====注意==== | |
实际上,本例是一个较为复杂的体系,复杂之处,在于其中的Ni片段是开壳层的。本例中,没有对各个片段进行单独指定其电子占据方式。实际计算过程中,如果某个片段是开壳层的,或者非满占据的,片段计算完毕之后,都要检查片段的占据方式是不是正确。请参考[[adf:openshellfragmentanalysis]] | |
| |
其他详细的功能,参考:[[adf:hbondanalysis|]] | |