计算化学公社

标题: 新手求助,orca利用多参考计算能量 [打印本页]

作者
Author:
lm1997    时间: 2022-6-18 15:39
标题: 新手求助,orca利用多参考计算能量
本帖最后由 lm1997 于 2022-6-18 17:22 编辑

各位老师们好,小生最近学习使用orca用多参考态计算能量(某烷加氧过程),状况频出,特来请教。
我用casscf对反应分子进行了柔性扫描,原想利用mrci计算下能量一直遇到不收敛的情况,就改换成了NEVPT2,但是遇到计算很快的问题大概几十秒就完成了计算,我导估计每个点大概需要几个小时的计算,所以我觉得可能是我输入文件的错误,故向各位老师请教一下,希望各位老师批评指正。输入文件如下:

(, 下载次数 Times of downloads: 21)

(, 下载次数 Times of downloads: 35)

作者
Author:
zjxitcc    时间: 2022-6-18 15:51
体系多少原子、多少基函数
作者
Author:
lm1997    时间: 2022-6-18 15:52
zjxitcc 发表于 2022-6-18 15:51
体系多少原子、多少基函数

C6H3O3
作者
Author:
wzkchem5    时间: 2022-6-18 15:55
NEVPT2本来就是远快于MRCI的。如果是比CASSCF还快,那才是有问题。
如果还不确定,贴完整输出文件上来
作者
Author:
hebrewsnabla    时间: 2022-6-18 15:56
那很正常,就是这么快
作者
Author:
lm1997    时间: 2022-6-18 15:58
wzkchem5 发表于 2022-6-18 15:55
NEVPT2本来就是远快于MRCI的。如果是比CASSCF还快,那才是有问题。
如果还不确定,贴完整输出文件上来

老师已经上传
作者
Author:
lm1997    时间: 2022-6-18 15:59
hebrewsnabla 发表于 2022-6-18 15:56
那很正常,就是这么快

第一次做这个,我觉得这种方法起码都要几十分钟起步的,竟然是秒级别的,超出我的认知,谢谢老师
作者
Author:
wzkchem5    时间: 2022-6-18 16:10
lm1997 发表于 2022-6-18 08:58
老师已经上传

所有计算步骤都正常结束了,结果确实是可以用的。
不过既然SC-NEVPT2做得这么快的话,建议用更可靠的FIC-NEVPT2,SC-NEVPT2不满足轨道酉不变性,做势能面可能出现势能面不连续的问题
作者
Author:
lm1997    时间: 2022-6-18 16:21
wzkchem5 发表于 2022-6-18 16:10
所有计算步骤都正常结束了,结果确实是可以用的。
不过既然SC-NEVPT2做得这么快的话,建议用更可靠的FIC ...

好的老师,我去学习一下,十分感谢各位老师的帮助,解决小生心中的疑惑。
作者
Author:
zjxitcc    时间: 2022-6-18 16:30
计算时间没问题。记得检查CASSCF活性轨道(算前、算后)合理性,不要当成HF/MP2那样黑箱式地算下去。
作者
Author:
ionexchangeC    时间: 2022-6-18 16:35
mrci的不收敛是指casscf步不收敛还是后面的mrci步不收敛?
作者
Author:
lm1997    时间: 2022-6-18 17:22
ionexchangeC 发表于 2022-6-18 16:35
mrci的不收敛是指casscf步不收敛还是后面的mrci步不收敛?

输出文件 已经上传 您看一下
作者
Author:
ionexchangeC    时间: 2022-6-18 17:38
lm1997 发表于 2022-6-18 17:22
输出文件 已经上传 您看一下

为什么mrci计算时不用ORCA5?ORCA5的CASSCF能收敛。
另外,做NEVPT2计算,基组至少3zeta
作者
Author:
lm1997    时间: 2022-6-18 21:04
ionexchangeC 发表于 2022-6-18 17:38
为什么mrci计算时不用ORCA5?ORCA5的CASSCF能收敛。
另外,做NEVPT2计算,基组至少3zeta

是用的 orca5.01算的 大概我有三个结构只有一个结构能算出来,其他的算到中间就卡住了,需要算dz和tz之后外推到完全基组




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