caspt2优化自由基出错

Will

这几天拿molpro的caspt2/vtz优化一个自由基老是会优化到直线状态，扩大活化空间又会出现"The problem occurs in muwvfc:spieig"的错误。看选的轨道是很正常的pai轨道，加上iprint,micro看收敛圈数也正常，请问到底是哪儿出了问题？现在知道的是这个自由基的基态和第一激发态的能量很接近。以下是未扩大活化空间的输入文件，谢谢！
***,occno
memory,5500,m;
gprint,civector;

geometry={
O                 -2.29939500    0.15648300    0.00000000
C                  0.00000000    0.35281800    0.00000000
C                  1.23286400   -0.07298500    0.00000000
O                  2.38335200   -0.20884800    0.00000000
N                 -1.15269100   -0.18001200    0.00000000
}

basis=VTZ

{hf
save,2100.2}

{multi
occ,15,4
closed,12,1
wf,35,2,1
start,2100.2
save,2140.2}

optg

put,molden,occno_caspt2.molden

xujian

感觉直线状态可能是存在虚频？

Y.zhang

楼主这个问题解决了吗？遇到同样的问题

zjxitcc

看到楼上有人提问，这个帖子被翻出来，我就顺便评论几句，楼主展示的文件内容并没有CASPT2部分（不知道是漏复制了rs2还是啥），当前的信息只能看作是CASSCF结构优化。我用MOKIT做了一个自动CASSCF单点计算，输入文件（OCCON.gjf）内容如下

1. %mem=48GB
   
2. %nprocshared=48
   
3. #p CASSCF/cc-pVTZ
   
4. 
   
5. mokit{CASSCF_prog=Molpro}
   
6. 
   
7. 0 2
   
8. O      -2.299395   0.156483   0.0
   
9. C       0.000000    0.352818   0.0
   
10. C       1.232864  -0.072985    0.0
    
11. O       2.383352  -0.208848    0.0
    
12. N      -1.152691  -0.180012    0.0
    

复制代码

约2 min算完，从输出文件中可以看到程序自动确定活性空间为CASSCF(7,7)，同时得到CASSCF自然轨道文件 OCCON_uhf_gvb12_CASSCF_NO.fch (979.69 KB, 下载次数 Times of downloads: 1)

2022-11-6 22:55 上传 Uploaded

点击下载
Click to download

执行

1. fch2com OCCON_uhf_gvb12_CASSCF_NO.fch

复制代码

生成文件Molpro输入文件OCCON_uhf_gvb12_CASSCF_NO.com和轨道文件occon_uhf_gvb1.a。打开该com文件，将{matrop及以下替换成结构优化关键词

1. {matrop;
   
2. read,mo,ORB,file=occon_uhf_gvb1.a;
   
3. save,mo,2140.2,ORBITALS}
   
4. {CASSCF;closed,14;occ,21;NoExtra}
   
5. optg

复制代码

提交给Molpro，轨道文件自动识别，计算也就几分钟吧，CASSCF(7,7)结构优化顺利完成，体系不是直线的，优化好的结构如下
OCCON_CAS77_opted.gjf (290 Bytes, 下载次数 Times of downloads: 5)

2022-11-6 22:55 上传 Uploaded

点击下载
Click to download