计算化学公社

标题: Molpro 计算NO的时候如何确定occ不可约表示的占据轨道数目及总对称性 [打印本页]

作者
Author: wulalalana 时间: 2020-8-31 21:17
标题: Molpro 计算NO的时候如何确定occ不可约表示的占据轨道数目及总对称性
本人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'
}

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

作者
Author: wulalalana 时间: 2020-9-1 11:02

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

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

作者
Author: wzkchem5 时间: 2020-9-1 16:23

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

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

作者
Author: wulalalana 时间: 2020-9-1 19:56

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

好的，谢谢

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