求助molpro的算cas的out文件中的natural orbital怎么与CI vector对应看

woshiwhm

求助molpro的算cas的out文件中的natural orbital怎么与CI vector对应看

各位老师好，我计算的是六个原子为CHHHIF;casscf(22,15)   frozen：3；closed：11；occ：26；
想知道第二个组态222222222202000     -0.1845906中的最后一个双占是什么，请问对应NATURAL ORBITALS的*.1呢，

NATURAL ORBITALS                                       
================                                       

   Orb     Occ        Energy       Coefficients

                                   1 1s      1 1s      1 1s      1 1s      1 2px     1 2py     1 2pz     1 2px     1 2py     1 2pz  
                                   1 2px     1 2py     1 2pz     1 3d0     1 3d2-    1 3d1+    1 3d2+    1 3d1-    1 3d0     1 3d2-
                                   1 3d1+    1 3d2+    1 3d1-    2 1s      2 1s      2 1s      2 2px     2 2py     2 2pz     2 2px  
                                   2 2py     2 2pz     3 1s      3 1s      3 1s      3 2px     3 2py     3 2pz     3 2px     3 2py  
                                   3 2pz     4 1s      4 1s      4 1s      4 2px     4 2py     4 2pz     4 2px     4 2py     4 2pz  
                                   5 1s      5 1s      5 1s      5 1s      5 1s      5 2px     5 2py     5 2pz     5 2px     5 2py  
                                   5 2pz     5 2px     5 2py     5 2pz     5 2px     5 2py     5 2pz     5 3d0     5 3d2-    5 3d1+
                                   5 3d2+    5 3d1-    5 3d0     5 3d2-    5 3d1+    5 3d2+    5 3d1-    5 3d0     5 3d2-    5 3d1+
                                   5 3d2+    5 3d1-    6 1s      6 1s      6 1s      6 1s      6 2px     6 2py     6 2pz     6 2px  
                                   6 2py     6 2pz     6 2px     6 2py     6 2pz     6 3d0     6 3d2-    6 3d1+    6 3d2+    6 3d1-
                                   6 3d0     6 3d2-    6 3d1+    6 3d2+    6 3d1-

   4.1  2.00000    -5.358271     0.000000  0.000000 -0.000002  0.000004  0.000000 -0.000006  0.000000 -0.000080 -0.001744 -0.000004
                                 0.000034  0.000750  0.000003  0.000000  0.000001 -0.000003  0.000021 -0.000062  0.000000 -0.000034
                                 0.000000 -0.000729 -0.000002 -0.000283 -0.004051  0.001573  0.000000  0.000005 -0.000002 -0.000026
                                -0.001369 -0.000662  0.000152  0.002187 -0.000849 -0.000001  0.000008  0.000001 -0.000399 -0.000790
                                 0.000356  0.000130  0.001865 -0.000724  0.000002  0.000008  0.000001  0.000291 -0.000758  0.000304
                                 0.000005  0.000010 -0.000001  0.000003 -0.000005  0.045847  0.998880  0.000486 -0.000754 -0.016437
                                -0.000008 -0.000028 -0.000618  0.000001  0.000006  0.000135  0.000002 -0.000003 -0.000001  0.000022
                                -0.000006  0.000478  0.000000  0.000001  0.000002  0.000011  0.000043  0.000000 -0.000001  0.000029
                                -0.000013  0.000624  0.000000 -0.000009  0.000000  0.000000  0.000051  0.001120  0.000002  0.000005
                                 0.000103  0.000001  0.000027  0.000587 -0.000001  0.000000  0.000000 -0.000001  0.000000 -0.000020
                                 0.000000  0.000000  0.000007  0.000001  0.000161
***
26.1  0.01794     1.348119    -0.019712 -0.156324 -0.012180  0.111953  0.000018  0.000537 -0.245272  0.000005 -0.000094  0.054510
                                -0.000006 -0.000358  0.129429 -0.024950  0.000000  0.000005 -0.000001  0.000091  0.008322  0.000001
                                -0.000002  0.000018 -0.000047  0.006515  0.007713  0.007293  0.000000 -0.000421 -0.000654 -0.000002
                                 0.008053  0.037644  0.006560  0.007545  0.007520  0.000360  0.000208 -0.000648 -0.006969 -0.004060
                                 0.037591  0.006565  0.007573  0.007522 -0.000360  0.000208 -0.000648  0.006973 -0.004065  0.037591
                                -0.016993  0.110484 -0.044134  0.016884 -0.257848  0.000000  0.000012 -0.004942  0.000002  0.000174
                                -0.028453 -0.000003 -0.000172 -0.171143  0.000001 -0.000210  0.028654 -0.105170  0.000000 -0.000003
                                 0.000000 -0.000074  0.160042  0.000000  0.000004  0.000001  0.000085  0.082389  0.000000  0.000005
                                -0.000004  0.000018  0.127292  0.676889 -0.591960 -0.120105  0.000023  0.000965  1.181540 -0.000025
                                -0.000993 -1.248009 -0.000006 -0.000215 -0.317275  0.020127  0.000000  0.000001  0.000000  0.000021
                                 0.034649  0.000000  0.000001  0.000000  0.000067

CI vector
=========

222222222220000      0.9634624
222222222202000     -0.1845906
222222222020200     -0.0595373
222222220220020     -0.0594988  

qinjiu

如果你是Cs点群，有两个对称性A`和A``，   n.1的意思就是第n个A`， n.2的意思就是第n个A``。
如果是C2v点群，有四个对称性A1,B1,B2,A2, n.1就是第n个A1, n.2就是第n个B2 ……以此类推
frozen3就是冻结1.1，2.1,3.1的轨道双占据不优化，能量从上一个参考轨道读取
closed11就是前11个.1的自然轨道占据数为2
occ26，closed11, 所以你的CI vector 有15个数字，分别代表从 12.1~26.1的占据数，这四个结果代表了不同的组态成分的比重
建议好好看一下说明书

woshiwhm

qinjiu 发表于 2021-5-25 17:15
如果你是Cs点群，有两个对称性A`和A``，   n.1的意思就是第n个A`， n.2的意思就是第n个A``。
如果是C2v点 ...

感谢您的解答，我已经了解这些设置以及对应的含义了，想再问下您想把第一个组态：222222222220000的轨道形状画出来应该怎么操作。我只会打开molden文件主功能0看轨道形状。

sobereva

woshiwhm 发表于 2021-5-27 09:35
感谢您的解答，我已经了解这些设置以及对应的含义了，想再问下您想把第一个组态：222222222220000的轨道 ...

所有组态都是共享同一套轨道，不是每个组态有各自的轨道