MRCC中采用ccsd(t)和ccsdt(q)两种方法无法计算出甲基的能量

王晓东 · 发表于 Post on 2021-1-25 17:34:05

本帖最后由王晓东于 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

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

wzkchem5 · 发表于 Post on 2021-1-25 17:35:50

你之前在这台机子上正常跑通过MRCC计算吗

王晓东 · 发表于 Post on 2021-1-25 18:53:25

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

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

sobereva · 发表于 Post on 2021-1-25 20:12:28

这根本就没体现是高阶激发态
直接用MRCC计算，没必要用molpro调用。具体有什么报错自行看提示中相应的文件

zjxitcc · 发表于 Post on 2021-1-25 20:40:14

我也想问哪里能体现出你在算“高阶激发态”。。

biogon · 发表于 Post on 2021-1-26 13:36:34

这也没有体现用EOM或者LR，只用耦合簇不就是算基态

wzkchem5 · 发表于 Post on 2021-1-26 20:44:32

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

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

biogon · 发表于 Post on 2021-1-27 00:55:11

本帖最后由 biogon 于 2021-1-27 00:56 编辑

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

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

wzkchem5 · 发表于 Post on 2021-1-27 09:33:09

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

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

biogon · 发表于 Post on 2021-1-27 09:38:09

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

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

wzkchem5 · 发表于 Post on 2021-1-27 09:39:57

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

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

wzkchem5 · 发表于 Post on 2021-1-27 09:42:04

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外推到基组极限。

王晓东 · 发表于 Post on 2021-1-27 13:30:20

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

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

王晓东 · 发表于 Post on 2021-1-27 13:31:32

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

不是要对比的

zjxitcc · 发表于 Post on 2021-1-27 13:50:16

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[MRCC/CFOUR] MRCC中采用ccsd(t)和ccsdt(q)两种方法无法计算出甲基的能量

浏览过的版块