计算化学公社

标题: MRCC中采用ccsd(t)和ccsdt(q)两种方法无法计算出甲基的能量 [打印本页]

作者
Author: 王晓东 时间: 2021-1-25 17:34
标题: MRCC中采用ccsd(t)和ccsdt(q)两种方法无法计算出甲基的能量
本帖最后由王晓东于 2021-1-27 14:27 编辑

MRCC中采用ccsd(t)和ccsdt(q)两种方法无法计算出甲基的能量，输入文件如下：
gprint, orbitals
memory, 2000, m
angstrom

geometry={
C                0.00000000 -0.00000800 0.00000100
H                0.93276200 -0.53849100 -0.00000200
H                -0.93275800 -0.53849800 -0.00000200
H                -0.00000400 1.07703500 -0.00000300}

basis=cc-pvdz
rhf
basis=aug-cc-pvtz
rhf
uccsd(t)
uhf
mrcc,method=ccsd(t)
mrcc,method=ccsdt(q)
---
输出out文件末尾：
MRCC Input end
MRCC now being run, temporary output in /public3/home/sc55389/tmp/molpro.Ac6LzCpp5Q/mrcc_35726/mrcc.out

一直停在这，是算不动吗？望各位老师指正

作者
Author: wzkchem5 时间: 2021-1-25 17:35
你之前在这台机子上正常跑通过MRCC计算吗

作者
Author: 王晓东 时间: 2021-1-25 18:53

wzkchem5 发表于 2021-1-25 17:35
你之前在这台机子上正常跑通过MRCC计算吗

有成功过，比如计算甲烷的时候就成功了

作者
Author: sobereva 时间: 2021-1-25 20:12
这根本就没体现是高阶激发态
直接用MRCC计算，没必要用molpro调用。具体有什么报错自行看提示中相应的文件

作者
Author: zjxitcc 时间: 2021-1-25 20:40
我也想问哪里能体现出你在算“高阶激发态”。。

作者
Author: biogon 时间: 2021-1-26 13:36
这也没有体现用EOM或者LR，只用耦合簇不就是算基态

作者
Author: wzkchem5 时间: 2021-1-26 20:44

biogon 发表于 2021-1-26 13:36
这也没有体现用EOM或者LR，只用耦合簇不就是算基态

估计楼主指的是高阶激发，不是高阶激发态

作者
Author: biogon 时间: 2021-1-27 00:55
本帖最后由 biogon 于 2021-1-27 00:56 编辑

wzkchem5 发表于 2021-1-26 20:44
估计楼主指的是高阶激发，不是高阶激发态

在群里了说是要算基态的能量，对比CCSD(T和CCSDT(Q，不知道这种对比有何意义

作者
Author: wzkchem5 时间: 2021-1-27 09:33

biogon 发表于 2021-1-27 00:55
在群里了说是要算基态的能量，对比CCSD(T和CCSDT(Q，不知道这种对比有何意义

这么小的体系与其做CCSDT(Q)，不如做个selected CI，直接得到近似FCI结果

作者
Author: biogon 时间: 2021-1-27 09:38

wzkchem5 发表于 2021-1-27 09:33
这么小的体系与其做CCSDT(Q)，不如做个selected CI，直接得到近似FCI结果

这个体系似乎用FCI也能算的动？还有ICE-CI之类的也行吧

作者
Author: wzkchem5 时间: 2021-1-27 09:39

biogon 发表于 2021-1-27 09:38
这个体系似乎用FCI也能算的动？还有ICE-CI之类的也行吧

cc-pVDZ能算动，aug-cc-pVTZ不确定，可能也行。
ICE-CI就是属于selected CI的一种

作者
Author: wzkchem5 时间: 2021-1-27 09:42

biogon 发表于 2021-1-27 09:38
这个体系似乎用FCI也能算的动？还有ICE-CI之类的也行吧

其实说到这个，这种体系CCSD(T)的方法误差应该明显小于aug-cc-pVTZ的基组误差，与其做CCSDT(Q)/aug-cc-pVTZ，还不如做CCSD(T)/aug-cc-pVTZ、CCSD(T)/aug-cc-pVQZ外推到基组极限。

作者
Author: 王晓东 时间: 2021-1-27 13:30

wzkchem5 发表于 2021-1-27 09:42
其实说到这个，这种体系CCSD(T)的方法误差应该明显小于aug-cc-pVTZ的基组误差，与其做CCSDT(Q)/aug-cc-pV ...

基组外推也做了，现在就是要用ccsd(t)，ccsdt(q)这两种方法算，不知道输入文件应该怎么修改才能算出来。MOLPRO手册上说mrcc只能应用于闭壳层，是不是甲基就不能算了

作者
Author: 王晓东 时间: 2021-1-27 13:31

biogon 发表于 2021-1-27 00:55
在群里了说是要算基态的能量，对比CCSD(T和CCSDT(Q，不知道这种对比有何意义

不是要对比的

作者
Author: zjxitcc 时间: 2021-1-27 13:50

王晓东发表于 2021-1-27 13:30
基组外推也做了，现在就是要用ccsd(t)，ccsdt(q)这两种方法算，不知道输入文件应该怎么修改才能算出来。M ...

那你咋不看MRCC手册怎么写的

作者
Author: Liuai613 时间: 2021-8-13 22:43
您好，您的问题解决了吗？我最近正好遇到同样的情况，想要请教一下

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