计算化学公社

标题: caspt2优化自由基出错 [打印本页]
作者
Author: Will    时间: 2019-7-18 16:47
标题: caspt2优化自由基出错
这几天拿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

作者
Author: xujian    时间: 2019-7-19 09:19
感觉直线状态可能是存在虚频？
作者
Author: Y.zhang    时间: 2022-11-6 21:27
楼主这个问题解决了吗？遇到同样的问题
作者
Author: zjxitcc    时间: 2022-11-6 22:56
本帖最后由 zjxitcc 于 2022-11-6 22:58 编辑

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

  5. mokit{CASSCF_prog=Molpro}

  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自然轨道文件 (, 下载次数 Times of downloads: 1)
执行
  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)结构优化顺利完成，体系不是直线的，优化好的结构如下
(, 下载次数 Times of downloads: 6)





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3