计算化学公社

 找回密码 Forget password
 注册 Register

CCSD输入文件错误

查看数: 13425 | 评论数: 12 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2017-11-14 11:52

正文摘要:

我要用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 ...

回复 Reply

younggood 发表于 Post on 2017-11-15 09:01:29
sobereva 发表于 2017-11-14 21:19
谈谈赝势基组的选用
http://sobereva.com/373

谢谢老师
sobereva 发表于 Post on 2017-11-14 21:19:24
younggood 发表于 2017-11-14 17:49
非常感谢您,那么Rh的话用什么基组比较好呢

谈谈赝势基组的选用
http://sobereva.com/373
younggood 发表于 Post on 2017-11-14 17:49:09
本帖最后由 younggood 于 2017-11-14 17:54 编辑
liyuanhe211 发表于 2017-11-14 16:52
Aug-cc-pVTZ 中没有 Rh 的定义。
CCSD(T)配cc-pVDZ是搞笑,耗时又多,结果还不如双杂化/TZ or QZ级别
...

非常感谢您,那么Rh的话用什么基组比较好呢
liyuanhe211 发表于 Post on 2017-11-14 16:52:06
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还算的动
younggood 发表于 Post on 2017-11-14 16:36:28
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



beefly 发表于 Post on 2017-11-14 16:16:37
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更快更方便
younggood 发表于 Post on 2017-11-14 16:15:31
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
younggood 发表于 Post on 2017-11-14 16:12:21
liyuanhe211 发表于 2017-11-14 15:37
为什么Ni的基组不需要带F12呢?为什么不用aug-cc-pVTZ-F12

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

您好,我的文章里没有用过CCSD(T)-F12,一直用的M06算单点,审稿人让我补CC,然后我看到一篇发表的文献里用了CCSD(T)-F12(如截图),想模仿一下这么做
liyuanhe211 发表于 Post on 2017-11-14 15:37:43
younggood 发表于 2017-11-14 15:11
您好,打扰您一下,我不是很明白这个CCSD(T)方法,因为审稿人让我用这个方法再精确一下计算结果,然后我 ...

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

我不懂审稿人要干什么,而且你的文章里不是已经用CCSD(T)-F12算过了,为什么还被要求补CC?以及为什么之前算过CCSD(T)-F12现在又要拿高斯算?
sobereva 发表于 Post on 2017-11-14 12:57:24
用任何程序都不能脱离手册,这是基本原则,关键词不是想怎么写就能怎么写
支持耦合簇的F12计算的有Molpro、ORCA、Turbomole
beefly 发表于 Post on 2017-11-14 12:42:47
和高斯讲molpro的话?

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 07:22 , Processed in 0.183205 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list