差别

这里会显示出您选择的修订版和当前版本之间的差别。

--- adf:openshellftagmentanalysis_new2020 [2024/02/23 17:05] – [碎片轨道SFO的查看] liu.jun
+++ adf:openshellftagmentanalysis_new2020 [2024/02/23 17:37] (当前版本) – [碎片轨道SFO的查看] liu.jun
@@ 行 25: / 行 25: @@
 打开NOCV功能：
 {{ :adf:2020unrestrictedfrag03.png?600 |}}
-关闭对称性，因为NOCV分析不允许使用点群（因此，**对于有对称性的配合物，我们往往会把EDA和NOCV分开做**，两遍的结果是一样的，但是整体配合物的分子轨道编号、SFO编号，一个对应这里的Oh群，用Oh群对应的不可约表示编号，一个没有对称性，统一不可约表示符号为A）：
+关闭对称性，因为NOCV分析不允许使用点群（因此，**对于有对称性的配合物，我们往往会把EDA和NOCV分开做**，两遍的结果是一样的，但是整体配合物的分子轨道编号、SFO编号，都对应这里的Oh群的不可约表示编号，一个没有对称性，分子轨道和碎片轨道SFO统一的不可约表示符号都是A）：
 {{ :adf:2020unrestrictedfrag04.png?600 |}}
@@ 行 239: / 行 239: @@
 对与这些碎片轨道，尤其是**金属原子那个Region，对应着什么原子轨道，如何查看呢？**
-这里我们以采用了点群Oh的EDA为例，此时Sr碎片每个SFO也是按照不可约表示来分类列出的。每个原子轨道对应的p、d、f编号可以在原子对应这个片段作业的out文件里面去看：Output→Properties→SFO construction会有一个列表，列出每个碎片轨道对应的实际原子轨道编号。但是注意，这个编号是自然编号，例如1P，指第一个P轨道，如果没有冻芯近似（使用全电子基组），则是指2P轨道，类似地，1D是指3d轨道。
+===EDA中的不可约表示===
+对于有对称性的分子，EDA分析是支持点群设置的，整体、碎片沿用同样的点群（软件默认是这样设置的，但是保险起见，建议确认一样后，再去计算），计算完毕后，在配合物整体的*.out文件中，列出的配合物整体的分子轨道，以及碎片轨道（在*.out中称为SFO），都是按照这个点群去分类列出的。碎片轨道SFO可以在Output → Properties → SFO construction 中看到，注意编号。
+例如，Sr(CO)8采用Oh群计算，则包括如下不可约表示：
 <code>
+ A1.g
+ A2.g
+ E.g:1  E.g:2
+ T1.g:1  T1.g:2  T1.g:3
+ T2.g:1  T2.g:2  T2.g:3
+ A2.u
+ A1.u
+ E.u:1  E.u:2
+ T2.u:1  T2.u:2  T2.u:3
+ T1.u:1  T1.u:2  T1.u:3
+</code>
+主任务中Output → Properties → SFO construction给出:
+<code>
+......省略
                                        === A1.g ===
- Nr. of SFOs :   12
+......省略
- Cartesian basis functions that participate in this irrep (total number =    27) :
-     2     3     4     5     6     7     8     9    10
-    12    39    34    37    45    40    43    51    46
-    60    55    58    66    61    64
-    SFO  (index         Fragment          Generating    Expansion in Fragment Orbitals
+   A    20      20        --        15.710 au  CO            1.00     20 A1.g            1
-  indx  incl.CFs)   Occup   Orb.Energy   FragmentType  Coeff.   Orbital     on Fragment
+                             (     427.500 eV)
- --------------------------------------------------------------------------------------
+   B                      --        15.710 au
-       1      2.000    -590.007 au  Sr            1.00      1 S               1
+                             (     427.500 eV)
-                        (  -16054.900 eV)
+   A    21      21        --        53.281 au  CO            1.00     21 A1.g            1
-       2      2.000     -80.987 au  Sr            1.00      2 S               1
+                             (    1449.846 eV)
-                        (   -2203.763 eV)
+   B                      --        53.281 au
-       3      2.000     -12.958 au  Sr            1.00      3 S               1
+                             (    1449.846 eV)
-                        (    -352.617 eV)
+   A    22      22        --       126.617 au  CO            1.00     22 A1.g            1
-       4      2.000      -1.713 au  Sr            1.00      4 S               1
+                             (    3445.420 eV)
-                        (     -46.610 eV)
+   B                      --       126.617 au
-       5      2.000      -0.168 au  Sr            1.00      5 S               1
+                             (    3445.420 eV)
-                        (      -4.575 eV)
+   A    23      23       1.000    -589.802 au  Sr            1.00      1 A1.g            2
-       6       --         0.038 au  Sr            1.00      6 S               1
+                             (  -16049.317 eV)
-                        (       1.024 eV)
+   B                     1.000    -589.802 au
-       7       --         0.312 au  Sr            1.00      7 S               1
+                             (  -16049.326 eV)
-                        (       8.484 eV)
+   A    24      24       1.000     -80.767 au  Sr            1.00      2 A1.g            2
-       8       --         1.743 au  Sr            1.00      8 S               1
+                             (   -2197.781 eV)
-                        (      47.439 eV)
+   B                     1.000     -80.764 au
-       9       --        11.756 au  Sr            1.00      9 S               1
+.......省略
-                        (     319.901 eV)
-      10       --       441.428 au  Sr            1.00     10 S               1
-                        (   12011.878 eV)
-      11       --      3253.969 au  Sr            1.00     11 S               1
-                        (   88545.002 eV)
-      12       --     16322.070 au  Sr            1.00     12 S               1
-                        (  444146.115 eV)
 </code>
-可以看到，其中A1.g不可约表示中，对应的全是Sr的s轨道，例如6 A1.g就是6S。类似地，E.g:1里面的则是Dz2轨道：
+那么SFO编号为24的（这里有2列编号，如果不使用Frozen Core，则两列编号一样，否则会不一样，简单起见，建议用户如果搞不清楚，就不使用Frozen core），对应Sr这个碎片的2 A1.g这个轨道。而Sr这个2 A1.g轨道又是什么原子轨道呢？我们打开Sr这个碎片的能级图，鼠标放在Sr的分子轨道上（第二列），找到2 A1.g：
+{{ :adf:asfrg.png?650 }}
+可以看到成分实际上是Sr的2S轨道。
+===NOCV中的不可约表示===
+此时点群被关闭了，因此只有一个A不可约表示。Output → Properties → SFO construction 就是这种样子：
 <code>
-                                      === E.g:1 ===
+......省略
- Nr. of SFOs :    5
+                                       === A ===
- Cartesian basis functions that participate in this irrep (total number =    15) :
+......省略
-    34    37    45    40    43    51    46    49    60
+   A   415     415        --        53.549 au  CO            1.00     35 T1.u:3          1
-    58    66    61    64
+                             (    1457.137 eV)
+   B                      --        53.549 au
+                             (    1457.137 eV)
+   A   416     416        --       126.891 au  CO            1.00     36 T1.u:3          1
+                             (    3452.893 eV)
+   B                      --       126.891 au
+                             (    3452.893 eV)
+   A   417     417       1.000    -589.802 au  Sr            1.00      1 A1.g            2
+                             (  -16049.317 eV)
+   B                     1.000    -589.802 au
+                             (  -16049.326 eV)
+   A   418     418       1.000     -80.767 au  Sr            1.00      2 A1.g            2
+                             (   -2197.781 eV)
+   B                     1.000     -80.764 au
+                             (   -2197.699 eV)
+   A   419     419       1.000     -12.737 au  Sr            1.00      3 A1.g            2
+                             (    -346.591 eV)
+   B                     1.000     -12.736 au
+                             (    -346.559 eV)
+   A   420     420       1.000      -1.500 au  Sr            1.00      4 A1.g            2
-    SFO  (index         Fragment          Generating    Expansion in Fragment Orbitals
+.......省略
-  indx  incl.CFs)   Occup   Orb.Energy   FragmentType  Coeff.   Orbital     on Fragment
- --------------------------------------------------------------------------------------
-       1      2.000      -5.180 au  Sr            1.00      1 D:z2            1
-                        (    -140.960 eV)
-       2       --         0.046 au  Sr            1.00      2 D:z2            1
-                        (       1.247 eV)
-       3       --         0.789 au  Sr            1.00      3 D:z2            1
-                        (      21.479 eV)
-       4       --         5.701 au  Sr            1.00      4 D:z2            1
-                        (     155.138 eV)
-       5       --        36.611 au  Sr            1.00      5 D:z2            1
-                        (     996.241 eV)
 </code>
+这是所有SFO都堆在一起，不那么方便分析了。不过也可以看到418这个SFO，对应着Sr碎片的2 A1.g。由于此时碎片还是用了点群Oh的，因此打开Sr碎片作业的能级图，一样去找2 A1.g，一样可以看到是2S轨道。
+实际上EDA、NOCV计算可以沿用相同的碎片adf.rkf文件，而不需要重新去计算碎片，这样EDA和NOCV的一致性会严格得到保证，也节省了时间。既然是同一个adf.rkf文件，打开的能级图就是一样的了。

费米维基

用户工具

站点工具

差别

页面工具