计算化学公社

标题: 求助molpro的算cas的out文件中的natural orbital怎么与CI vector对应看 [打印本页]

作者
Author: woshiwhm 时间: 2021-5-10 13:50
标题: 求助molpro的算cas的out文件中的natural orbital怎么与CI vector对应看
求助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

作者
Author: qinjiu 时间: 2021-5-25 17:15
如果你是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的占据数，这四个结果代表了不同的组态成分的比重
建议好好看一下说明书

作者
Author: woshiwhm 时间: 2021-5-27 09:35

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

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

作者
Author: sobereva 时间: 2021-6-1 15:39

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)