计算化学公社

标题: CCSD输入文件错误 [打印本页]

作者
Author: younggood 时间: 2017-11-14 11:52
标题: CCSD输入文件错误
我要用CCSD(T)方法计算一个单点，这是我的输入文件：
%nprocshared=14%mem=10GB
# UCCSD(T)-F12a/gen gfinput pseudo=read

retsRhON2phphCC

1 1

Rh 0
aug-cc-pVTZ
****
C H O N 0
cc-pVDZ-F12
****

Rh 0
aug-cc-pVTZ

可是为什么一交上就出错呢，这是错误信息：

---------------------------------------
# UCCSD(T)-F12a/gen gfinput pseudo=read
---------------------------------------
QPErr --- A syntax error was detected in the input line.
# UCCSD(T)-F12a/gen gfinput pseudo=read
'
Last state="GCL"
TCursr= 1177 LCursr= 2
Error termination via Lnk1e in i:\gs9\l1.exe at Tue Nov 14 11:50:33 2017.
Job cpu time: 0 days 0 hours 0 minutes 1.0 seconds.
File lengths (MBytes): RWF= 8193 Int= 0 D2E= 0 Chk= 1 Scr= 1

这个方法是我在文献中看到的，应该木有错啊，大神们，请告诉我是咋回事呢

作者
Author: beefly 时间: 2017-11-14 12:42
和高斯讲molpro的话？

作者
Author: sobereva 时间: 2017-11-14 12:57
用任何程序都不能脱离手册，这是基本原则，关键词不是想怎么写就能怎么写
支持耦合簇的F12计算的有Molpro、ORCA、Turbomole

作者
Author: younggood 时间: 2017-11-14 15:11

beefly 发表于 2017-11-14 12:42
和高斯讲molpro的话？

您好，打扰您一下，我不是很明白这个CCSD(T)方法，因为审稿人让我用这个方法再精确一下计算结果 (, 下载次数 Times of downloads: 43) ，然后我模仿了一篇文献（如上），这是审稿人的原话：The functionals and basis sets employed in theDFT calculations are standard methods and can be considered as suitablefor the investigated mechanisms. However, a validation against higher levelmethods (truncated CC, e.g.) would have added more confidence to the choice ofM06. 请问您，我的输入文件该怎么构建呢？

作者
Author: liyuanhe211 时间: 2017-11-14 15:37

younggood 发表于 2017-11-14 15:11
您好，打扰您一下，我不是很明白这个CCSD(T)方法，因为审稿人让我用这个方法再精确一下计算结果，然后我 ...

为什么Ni的基组不需要带F12呢？为什么不用aug-cc-pVTZ-F12

我不懂审稿人要干什么，而且你的文章里不是已经用CCSD(T)-F12算过了，为什么还被要求补CC？以及为什么之前算过CCSD(T)-F12现在又要拿高斯算？

作者
Author: younggood 时间: 2017-11-14 16:12

liyuanhe211 发表于 2017-11-14 15:37
为什么Ni的基组不需要带F12呢？为什么不用aug-cc-pVTZ-F12

我不懂审稿人要干什么，而且你的文章里不是 ...

您好，我的文章里没有用过CCSD(T)-F12，一直用的M06算单点，审稿人让我补CC，然后我看到一篇发表的文献里用了CCSD(T)-F12（如截图），想模仿一下这么做

作者
Author: younggood 时间: 2017-11-14 16:15

liyuanhe211 发表于 2017-11-14 15:37
为什么Ni的基组不需要带F12呢？为什么不用aug-cc-pVTZ-F12

我不懂审稿人要干什么，而且你的文章里不是 ...

我想用高斯算：
%nprocshared=10
%mem=10GB
# CCSD(T)/gen gfinput pseudo=read

retsRhON2phphCC

1 1

Rh 0
aug-cc-pVTZ
****
C H O N 0
cc-pVDZ
****

Rh 0
aug-cc-pVTZ

可是这么写错误信息：Rotational constants (GHZ):    0.1289258    0.0929036    0.0812487
General basis read from cards:  (5D, 7F)
Atomic number out of range in ccpVTZ
Error termination via Lnk1e in /opt//g09/l301.exe at Tue Nov 14 15:53:12 2017.
Job cpu time:    0 days  0 hours  0 minutes  2.1 seconds.
File lengths (MBytes):  RWF=    5 Int=    0 D2E=    0 Chk=    1 Scr=    1

作者
Author: beefly 时间: 2017-11-14 16:16
F12是加速基组收敛的方法，CCSD(T)-F12用vnz级别的基组近似得到v(n+1)z基组级别的结果，计算量是CCSD(T)/vnz的几倍，一般远小于CCSD(T)/v(n+1)z

但是这种方法需要专门的基组，如果是自定义基组，最好用orca的dlpno-ccsd(t)，比CCSD(T)-F12更快更方便

作者
Author: younggood 时间: 2017-11-14 16:36

beefly 发表于 2017-11-14 16:16
F12是加速基组收敛的方法，CCSD(T)-F12用vnz级别的基组近似得到v(n+1)z基组级别的结果，计算量是CCSD(T)/vn ...

您好，我可以这么构建我的输入文件吗？
%nprocshared=10
%mem=10GB
# CCSD(T)/gen gfinput pseudo=read

retsRhON2phphCC

1 1

Rh 0
aug-cc-pVTZ
****
C H O N 0
cc-pVDZ
****

Rh 0
aug-cc-pVTZ

作者
Author: liyuanhe211 时间: 2017-11-14 16:52

younggood 发表于 2017-11-14 16:15
我想用高斯算：
%nprocshared=10
%mem=10GB

Aug-cc-pVTZ 中没有 Rh 的定义。
CCSD(T)配cc-pVDZ是搞笑，耗时又多，结果还不如双杂化/TZ or QZ级别
我看你写的是“RhON2phphCC”，俩苯环、Rh、ON2这么些原子一般机器高斯CCSD(T)、像样的基组肯定是算不动的，20个原子左右的话用8L说的ORCA里DLPNO-CCSD(T)/TZ还算的动

作者
Author: younggood 时间: 2017-11-14 17:49
本帖最后由 younggood 于 2017-11-14 17:54 编辑

liyuanhe211 发表于 2017-11-14 16:52
Aug-cc-pVTZ 中没有 Rh 的定义。
CCSD(T)配cc-pVDZ是搞笑，耗时又多，结果还不如双杂化/TZ or QZ级别
...

非常感谢您，那么Rh的话用什么基组比较好呢

作者
Author: sobereva 时间: 2017-11-14 21:19

younggood 发表于 2017-11-14 17:49
非常感谢您，那么Rh的话用什么基组比较好呢

谈谈赝势基组的选用
http://sobereva.com/373

作者
Author: younggood 时间: 2017-11-15 09:01

sobereva 发表于 2017-11-14 21:19
谈谈赝势基组的选用
http://sobereva.com/373

谢谢老师

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