计算化学公社

标题: molpro中casscf计算HCHO时的输入文件中OCC轨道是如何计算的呢？ [打印本页]

作者
Author: ttxx 时间: 2019-9-5 20:22
标题: molpro中casscf计算HCHO时的输入文件中OCC轨道是如何计算的呢？
本帖最后由 ttxx 于 2019-9-5 21:53 编辑

这是，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

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

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

作者
Author: yiwei 时间: 2019-10-24 15:23

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

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

作者
Author: lastzealot 时间: 2019-10-24 16:14

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

用hf算一下看输出文件

作者
Author: yiwei 时间: 2019-10-25 13:41

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

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

作者
Author: lastzealot 时间: 2019-10-25 14:50

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

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

作者
Author: yiwei 时间: 2019-10-25 19:42

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

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

作者
Author: lastzealot 时间: 2019-10-25 21:55
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

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

作者
Author: 橘垚垚 时间: 2022-5-17 22:08

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的空轨道

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