molpro中casscf计算HCHO时的输入文件中OCC轨道是如何计算的呢？

ttxx

这是，molpro 说明书上用casscf 计算 HCHO（C2v)） 的例子，为什么占据轨道写成， OCC, 7,2,3 ，而不是OCC,5,1,2 呢？还有这个链接里面https://www2.chemistry.msu.edu/c ... olpro_examples.html   是molpro计算 CH2（triplet） 的例子，输入文件中   wf 这一项中的 对称性，occ 轨道数，和closed 轨道数都是怎么得来的呢？[size=10.6667px]初学molpro，求各位老师指导，谢谢大家了！ @sobereva

lastzealot

HCHO是c2v对称性分子，有四个不可约表示，在molpro里其顺序为A1 B1 B2 A2 。
occ后面接的是各不可约表示对称性下的轨道数，比如A1下有7个轨道，B1下有2个轨道。没写就是没有对应轨道。这里面是既包括占据轨道，又包括非占据轨道。
closed对应不放入cas中的轨道数，双占据
occ-closed就是活性空间轨道数。
[总电子数-closed×2]就是活性空间电子数。

我记得是这样，manual上面都有。

yiwei

lastzealot 发表于 2019-10-19 08:50
HCHO是c2v对称性分子，有四个不可约表示，在molpro里其顺序为A1 B1 B2 A2 。
occ后面接的是各不可约表示对 ...

请问老师，可以通过什么看出来A1下7个轨道，B1下两个轨道呢？如果不设置occ的话有影响吗？

lastzealot

yiwei 发表于 2019-10-24 15:23
请问老师，可以通过什么看出来A1下7个轨道，B1下两个轨道呢？如果不设置occ的话有影响吗？

用hf算一下看输出文件

yiwei

lastzealot 发表于 2019-10-24 16:14
用hf算一下看输出文件

谢谢老师，请问是高斯输出文件这一部分吗？

lastzealot

yiwei 发表于 2019-10-25 13:41
谢谢老师，请问是高斯输出文件这一部分吗？

你这个是c1对称性 不是c2v
另外你不是用molpro吗

yiwei

lastzealot 发表于 2019-10-25 14:50
你这个是c1对称性 不是c2v
另外你不是用molpro吗

是的，我想用molpro算caspt2,之前听一个老师说高斯输入文件就能知道occ信息，我把我脑子里的两处信息给记混了……I am sorry...老师，麻烦您再帮我看一下，这是我molpro算的hf的输出文件，请问是看哪一部分呀？C1的应该怎么设置occ呢？

lastzealot

Full point group                 C1      NOp   1
对称性是c1

Alpha Orbitals:
       Occupied  (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
       Virtual   (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A)
Beta  Orbitals:
       Occupied  (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A)
       Virtual   (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A)
The electronic state is 3-A.
电子结构对称性为3A，就是有两个单电子的，每个单占据轨道不可与表示为A，A×A=A

Yin

请问楼主知道HCOH的OCC为什么写成7，2，3呢？这几个数字是如何确定的呢？

橘垚垚

Yin 发表于 2021-7-1 16:17
请问楼主知道HCOH的OCC为什么写成7，2，3呢？这几个数字是如何确定的呢？

occ后面的数字加起来超过了总的占据轨道数的话说明包含了一部分空轨道，因为这个是在做CASSCF的时候定义的，CASSCF的活性空间等于occ减去closed，所以该occ,7,2,3代表在occ5,1,2的基础上包括了2,1,1的空轨道