这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版上一修订版两侧同时换到之后的修订记录 | ||
adf:openshellftagmentanalysis_new2020 [2024/02/23 17:02] – [参数设置] liu.jun | adf:openshellftagmentanalysis_new2020 [2024/02/23 17:37] – [碎片轨道SFO的查看] liu.jun | ||
---|---|---|---|
行 25: | 行 25: | ||
打开NOCV功能: | 打开NOCV功能: | ||
{{ : | {{ : | ||
- | 关闭对称性,因为NOCV分析不允许使用点群(因此,**对于有对称性的配合物,我们往往会把EDA和NOCV分开做**,两遍的结果是一样的,但是整体配合物的分子轨道编号、SFO编号,一个对应这里的Oh群,用Oh群对应的不可约表示编号,一个没有对称性,统一不可约表示符号为A): | + | 关闭对称性,因为NOCV分析不允许使用点群(因此,**对于有对称性的配合物,我们往往会把EDA和NOCV分开做**,两遍的结果是一样的,但是整体配合物的分子轨道编号、SFO编号,都对应这里的Oh群的不可约表示编号,一个没有对称性,分子轨道和碎片轨道SFO统一的不可约表示符号都是A): |
{{ : | {{ : | ||
行 239: | 行 239: | ||
对与这些碎片轨道,尤其是**金属原子那个Region,对应着什么原子轨道,如何查看呢?** | 对与这些碎片轨道,尤其是**金属原子那个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群计算,则包括如下不可约表示: | ||
< | < | ||
+ | A1.g | ||
+ | A2.g | ||
+ | | ||
+ | | ||
+ | | ||
+ | A2.u | ||
+ | A1.u | ||
+ | | ||
+ | | ||
+ | | ||
+ | </ | ||
+ | 主任务中Output → Properties → SFO construction给出: | ||
+ | < | ||
+ | ......省略 | ||
=== A1.g === | === A1.g === | ||
- | Nr. of SFOs : 12 | + | ......省略 |
- | | + | |
- | 1 | + | |
- | | + | |
- | | + | |
- | SFO (index | + | A 20 20 |
- | indx incl.CFs) | + | |
- | -------------------------------------------------------------------------------------- | + | B |
- | | + | |
- | ( -16054.900 eV) | + | A 21 21 |
- | | + | |
- | ( -2203.763 eV) | + | B |
- | 3 | + | |
- | ( -352.617 eV) | + | A 22 22 |
- | 4 | + | |
- | ( | + | B |
- | | + | |
- | ( -4.575 eV) | + | A 23 23 |
- | 6 | + | |
- | ( 1.024 eV) | + | B |
- | 7 | + | |
- | ( 8.484 eV) | + | A 24 24 |
- | 8 | + | |
- | ( 47.439 eV) | + | B |
- | 9 | + | .......省略 |
- | ( 319.901 eV) | + | |
- | | + | |
- | ( 12011.878 eV) | + | |
- | | + | |
- | ( 88545.002 eV) | + | |
- | | + | |
- | | + | |
</ | </ | ||
- | 可以看到,其中A1.g不可约表示中,对应的全是Sr的s轨道。类似地,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: | ||
+ | |||
+ | 可以看到成分实际上是Sr的2S轨道。 | ||
+ | |||
+ | ===NOCV中的不可约表示=== | ||
+ | 此时点群被关闭了,因此只有一个A不可约表示。Output → Properties → SFO construction 就是这种样子: | ||
< | < | ||
- | === E.g:1 === | + | ......省略 |
- | Nr. of SFOs : | + | === A === |
- | Cartesian basis functions that participate in this irrep (total number = 15) : | + | ......省略 |
- | 39 34 37 45 40 43 51 46 49 60 | + | |
- | 55 58 66 61 64 | + | |
+ | B -- 53.549 au | ||
+ | ( | ||
+ | A | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | B | ||
+ | | ||
+ | | ||
- | SFO (index | + | .......省略 |
- | indx incl.CFs) | + | |
- | | + | |
- | | + | |
- | ( -140.960 eV) | + | |
- | | + | |
- | ( 1.247 eV) | + | |
- | | + | |
- | ( 21.479 eV) | + | |
- | | + | |
- | ( | + | |
- | | + | |
- | ( | + | |
</ | </ | ||
+ | |||
+ | 这是所有SFO都堆在一起,不那么方便分析了。不过也可以看到418这个SFO,对应着Sr碎片的2 A1.g。由于此时碎片还是用了点群Oh的,因此打开Sr碎片作业的能级图,一样去找2 A1.g,一样可以看到是2S轨道。 | ||
+ | |||
+ | 实际上EDA、NOCV计算可以沿用相同的碎片adf.rkf文件,而不需要重新去计算碎片,这样EDA和NOCV的一致性会严格得到保证,也节省了时间。 |