Molpro 计算NO的时候如何确定occ不可约表示的占据轨道数目及总对称性

wulalalana

本人Molpro初学者：在用casscf计算NO的势能曲线，occ应该如何写？NO是开壳层分子，是不是不需要写closed？下面是我编辑的molpro输入文件！谢谢！
molpro输入文件：
r=[0.5,1,1.5,2,2.5,3,3.5,4,4.5] ang


i=1
geometry={N
          O,1,r(i)}


basis=av5z,F=avtz


hf                     


do i=1,#r               
{multi
closed,x,x,x,x
occ,   x,x,x,x
wf;state,2                 
canonical,ci}


{rs2,MIX=2
wf;state,2}               


e1_caspt2(i)=energu(1)      
e2_caspt2(i)=energu(2)      


e1_mscaspt2(i)=msenergy(1)  
e2_mscaspt2(i)=msenergy(2)  


{rs2,MIX=2,XMS=1
wf;state,2}               


e1_xmscaspt2(i)=msenergy(1)  
e2_xmscaspt2(i)=msenergy(2)  
enddo


{table,r,e1_caspt2,e2_caspt2,e1_mscaspt2,e2_mscaspt2,e1_xmscaspt2,e2_xmscaspt2
title,MS-MR-CASPT2 for NO
plot,file='no_mr_mscaspt2.plot'
}

wzkchem5

closed指的是活性空间外双占轨道的数目，和分子本身是否是开壳层无关。开壳层分子CAS一般会把所有单占轨道以及一部分双占轨道、空轨道纳入活性空间，其余的占据轨道仿照ROHF处理，即alpha、beta自旋用同一套轨道。

wulalalana

wzkchem5 发表于 2020-8-31 21:52
closed指的是活性空间外双占轨道的数目，和分子本身是否是开壳层无关。开壳层分子CAS一般会把所有单占轨道 ...

[img][/img]
closed 后这个是闭壳层轨道数，如果是开壳层分子也有这个吗

wzkchem5

wulalalana 发表于 2020-9-1 11:02
[/img]
closed 后这个是闭壳层轨道数，如果是开壳层分子也有这个吗

有的，“闭壳层”这个词修饰的是“轨道”，不是“分子”。在这个语境下闭壳层轨道和双占轨道是一个意思

wulalalana

wzkchem5 发表于 2020-9-1 16:23
有的，“闭壳层”这个词修饰的是“轨道”，不是“分子”。在这个语境下闭壳层轨道和双占轨道是一个意思

好的，谢谢