计算化学公社

标题: MOLPRO做CASSCF计算挑选活性轨道的问题 [打印本页]

作者
Author: 量化菜鸟 时间: 2021-6-30 10:25
标题: MOLPRO做CASSCF计算挑选活性轨道的问题
各位老师各位同学大家好，最近在看一篇文献，是关于双原子分子C2的，下图中是关于它的一些基本信息。导师让重复一下文献里的计算，计算的第一步是用MOLPRO 做CASSCF计算，之前对MOLPRO 有一些简单的了解，看了MOLPRO手册和快速入门会写一些形式比较简单的输入卡，但是我自己也清楚的知道我现在对MOLPRO的使用还停留在表面，还不是很会举一反三。先用分子文献里的基组用高斯看了一下C2分子的轨道，发现occ轨道为(SGG) (SGU) (SGG) (SGU) (PIU) (PIU)，与文献中给出的C2分子的基态电子组态一致，文献中把(SGG) (SGU) 这两个轨道冻结了（这个我比较好理解），剩下的(SGG) (SGU) (PIU) (PIU)与1PIG, 3SGU, 1DLTG，1DLTU共同组成活化空间（即下文中红色部分），用高斯查看得到的轨道信息如下：
Orbital symmetries:
   Occupied  (SGG) (SGU) (SGG) (SGU) (PIU) (PIU)
   Virtual (SGG) (SGU) (SGG) (PIG) (PIG) (PIU) (PIU) (SGG)
               (PIG) (PIG) (SGU) (DLTG) (DLTG) (SGG) (SGU) (PIU)
               (PIU) (DLTU) (DLTU) (PIG) (PIG) (SGU) (PIU) (PIU)
               (SGG) (PIG) (PIG) (SGU) (SGG) (DLTG) (DLTG) (PIU)
               (PIU) (PHIU) (PHIU) (DLTU) (DLTU) (SGU) (SGG)
               (PIG) (PIG) (DLTG) (DLTG) (PHIG) (PHIG) (SGU)
               (PIU) (PIU) (DLTU) (DLTU) (PIG) (PIG) (SGG) (SGU)
               (PIU) (PIU) (SGG) (SGU) (DLTG) (DLTG) (PIG) (PIG)
               (SGG) (PIU) (PIU) (PHIU) (PHIU) (DLTG) (DLTG)
               (DLTU) (DLTU) (PHIG) (PHIG) (PIG) (PIG) (SGU)
               (PIU) (PIU) (DLTU) (DLTU) (SGU) (SGG) (PIG) (PIG)
               (SGG) (SGU)
The electronic state is 1-SGG.

现在我有两个疑问：1，*** , C2
print,basis, orbitals
geometry={
c1
c2,c1,r
}

r=2.50 bohr
basis=avtz

hf

{casscf
occ,3,1,1,1,4,1,1,1
closed,1,0,0,0,1,0,0,0
wf,12,1,0
}
这是我根据文献中信息直接写的输入卡。可以直接这样写吗？
2，仔细查了MOLPRO手册，有rotate这个操作，但是对这个操作具体应该怎么用理解的不是很深刻，如果上文问题1中那样写的输入卡是不对的，那是不是应该用rotate命令交换轨道？针对这种情况应该怎么正确使用rotate这个命令呢？
说的有点啰嗦了，还请大家多多指教，谢谢。

作者
Author: evenn 时间: 2023-1-10 13:22
请问现在解决了吗？求教如何确定occ和closed后面的数字

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