orca算完NEVPT2后如何读取各组态占比以及各轨道对称性？

lycheeho

大家好！我用orca计算了有机π共轭分子的开壳层单重态，采用了CAS(4,4)结合DLPNP-NEVPT2多参考方法，附件是输出结果。请教各位老师该如何从结果得到各组态占比以及各轨道对称性信息呢？如文献中所示这样。本人经验不足，谢谢大家！

18217265596

NEVPT2会改变occupation吗？ 我以为只是基于cas的结果算更精准的单点

wzkchem5

ORCA论坛下载区里有一个CASSCF tutorial，看那个就行了

lycheeho

wzkchem5 发表于 2025-6-29 09:43
ORCA论坛下载区里有一个CASSCF tutorial，看那个就行了

好的谢谢

lycheeho

18217265596 发表于 2025-6-29 09:43
NEVPT2会改变occupation吗？ 我以为只是基于cas的结果算更精准的单点

对的，占据还是看casscf

lycheeho

wzkchem5 发表于 2025-6-29 09:43
ORCA论坛下载区里有一个CASSCF tutorial，看那个就行了

老师您好！我看完casscf有关手册还是有些不明白的地方，还请您指正：CAS(4,4)收敛后发现有4个单占据轨道；而输出态分布结果中，4个单占据轨道仅占24.9%，这该如何理解？以及这里“1111”占据对应的轨道是具体对应上述4个单占据轨道吗？谢谢~

wzkchem5

lycheeho 发表于 2025-6-29 16:15
老师您好！我看完casscf有关手册还是有些不明白的地方，还请您指正：CAS(4,4)收敛后发现有4个单占据轨道 ...

这个结果表示：你算的态有24.9%的比例是四个轨道都是单占据的，有22.0%的比例是前两个双占据、后两个不占据，等等

lycheeho

wzkchem5 发表于 2025-6-29 16:55
这个结果表示：你算的态有24.9%的比例是四个轨道都是单占据的，有22.0%的比例是前两个双占据、后两个不占 ...

那上面自然轨道的具体排布指的是占主要比重的“1111”，但是其比重并不明显高于其他态（比重差别不大），是这个意思吧？

wzkchem5

lycheeho 发表于 2025-6-29 17:11
那上面自然轨道的具体排布指的是占主要比重的“1111”，但是其比重并不明显高于其他态（比重差别不大）， ...

对，但我不知道你说的“上面自然轨道的具体排布”是什么意思，如果指的是那个列了各个轨道占据数、能量期望值的截图的话，那个是考虑了所有组态的贡献的，也就是不仅有1111还有2200等所有组态

hebrewsnabla

假如只包含2200,0022这样的组态，自然轨道占据数也是有可能在1.0左右的。并不是说1.0左右的自然轨道占据数就必然和单占据组态1111有关联。

lycheeho

wzkchem5 发表于 2025-6-29 18:20
对，但我不知道你说的“上面自然轨道的具体排布”是什么意思，如果指的是那个列了各个轨道占据数、能量期 ...

不好意思是我表达不清。具体排布指的是369 370 ...这几个在输出文件中的占据数；是的，我就是重现的文献图中的例子，看它文中这个图，1111（深绿色）占主要，但我算出来是占所有组态的24.9%，说明其实占比也没有很大，是想表达这个意思。

lycheeho

hebrewsnabla 发表于 2025-6-29 18:25
假如只包含2200,0022这样的组态，自然轨道占据数也是有可能在1.0左右的。并不是说1.0左右的自然轨道占据数 ...

谢谢回复！有点深奥了，刚开始接触多组态计算，那这个分子开壳层单重态基态的态分布以什么为主呢？以算出来具体的自然轨道占据数吗？

hebrewsnabla

你这个图里的深绿色不是2110么……

组态系数是和选择哪种轨道相关的。赝正则轨道、自然轨道、局域轨道下的组态系数可能完全不一样，那么主要的组态也会不一样。这个在orca手册里面“Final orbitals options”给了控制的选项。

另外文献的轨道形状也可能和你的不一样，这也会影响组态成分。

lycheeho

hebrewsnabla 发表于 2025-6-29 19:08
你这个图里的深绿色不是2110么……

组态系数是和选择哪种轨道相关的。赝正则轨道、自然轨道、局域轨道下 ...

是的，能级图里面主要组态是2110。我是按文献提供的原始计算文件重算的，它是用UNO作为CASSCF的输入。那既然不同轨道输出的组态系数不同，该以那个为主呢？而且上面计算出来组态比重完全没有2110这个结果。