应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 用cbs-qb3计算I时出错
101/1返回列表
查看: 2530  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

amin0241

金虫 (小有名气)

[求助] 用cbs-qb3计算I时出错 已有2人参与

我的输入信息如下：
%chk=hi
%mem=96mw
%nproc=2
#P MP2(direct)/GEN PSEUDO=READ

Title Card Required

0 1
I                  0.00000000    0.00000000    0.02990700
H                  0.00000000    0.00000000   -1.58505500




I 0
Lanl2dz
****
H 0
CBSB3 CBSExtrap=(MinPop)
****

I 0
Lanl2dz

出错信息如下：
Number           0
Base         20480
End          67072
End1         67072
Wr Pntr      20480
Rd Pntr      20480
Length       46592
Error termination in NtrErr:
NtrErr Called from FileIO.
回复此楼

» 收录本帖的淘帖专辑推荐

Gaussian

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-05-14 10:13:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhouwohua

银虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
PSEUDO=READ 是要读取chk文件里面的赝势信息
你的错误说明chk文件里面这个信息。
赞一下 回复此楼
学无止境；天天学习；努力向上。
2楼2014-05-14 11:04:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
第一个问题，分子描述和附加部分中间空太多行了，请严格按照输入格式写。
第二个问题，CBS-QB3应该不能这么用在赝势里，至少我尝试了几种写法都报语法错误了，具体怎么用，请翻阅资料或者google
赞一下 回复此楼
3楼2014-05-14 11:53:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amin0241

金虫 (小有名气)
引用回帖:
2楼: Originally posted by zhouwohua at 2014-05-14 11:04:49
PSEUDO=READ 是要读取chk文件里面的赝势信息
你的错误说明chk文件里面这个信息。

那应该用什么代替呢？

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
4楼2014-05-16 14:20:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amin0241

金虫 (小有名气)
引用回帖:
3楼: Originally posted by 枪下游魂 at 2014-05-14 11:53:14
第一个问题，分子描述和附加部分中间空太多行了，请严格按照输入格式写。
第二个问题，CBS-QB3应该不能这么用在赝势里，至少我尝试了几种写法都报语法错误了，具体怎么用，请翻阅资料或者google

如果对I直接用CBS也出错，基组不行，所以对它用赝势。第一个不是问题呢

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
5楼2014-05-16 14:24:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)
引用回帖:
5楼: Originally posted by amin0241 at 2014-05-16 14:24:31
如果对I直接用CBS也出错，基组不行，所以对它用赝势。第一个不是问题呢
...

你没懂我的意思：
1、我直接复制你给的输入文件去跑，高斯会报和你给出的一样的错误，原因就是你分子描述和附加输入部分中间空了太多行了，我删去多于的行数，不再报你贴出的这个错误，但是报赝势基组定义错误；
2、I 肯定用赝势来算，但是CBS方法在使用时是不需要输入基组的，而赝势定义时要分别定义不同的基组，所以我觉得不能直接在定义赝势的时候写CBS类方法，但是怎么用，我也没查到，所以你需要看看你的文献里是怎么说这个的
赞一下 回复此楼
6楼2014-05-16 14:44:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhouwohua

银虫 (小有名气)
引用回帖:
4楼: Originally posted by amin0241 at 2014-05-16 14:20:50
那应该用什么代替呢？
...

不好意思，之前没细看。以为是guess=read。我在高斯09说明说上看了一下关于CBS方法的说明：
这个方法关键词指定使用由 Petersson 等人开发的各种完备基组（Complete Basis Set；
CBS）方法，用于进行非常精确的能量计算[Nyden81, Petersson88, Petersson91, Petersson91a, Montgomery94,
Ochterski96, Montgomery99, Montgomery00]。 这 些 关 键 词 分 别 包 括 改 进 版 的 CBS-4 [Ochterski96,
Montgomery00]，CBS-Q//B3 [Montgomery99, Montgomery00]，和 CBS-APNO [Ochterski96]方法。使用这些关
键词时不需要指定基组。可以为 CBS-QB3 加上前缀 RO，请求 ROCBS-QB3 方法[Wood06]。
这些方法是复杂的能量计算，对指定的体系包括多个预先定义的计算。当指定了这些关
键词中的一个，即会自动进行多步的计算工作，并将最后计算的能量显示在输出中。
Opt=Maxcyc=n，QCISD=Maxcyc=n，或 CCSD=Maxcyc=n 关键词可以和这些关键词中
的任何一个结合使用，分别指定优化和 QCISD 或 CCSD循环的最大数量。
适用性
只能用于能量计算。CBS-4M和 CBS-QB3 只能用于第一和第二行原子；CBS-APNO 只
适用于第一行原子。ROCBS-QB3 [Wood06]用于使用 CBS 方法的限制性开壳层计算。


在不同温度下重新运行计算。下面的两步任务描述了在第二个（非常快的）CBS 计算中
使用不同温度的方法。这个任务在 298.15 K 和 300 K 的温度下计算 CBS-QB3 能量：
%Chk=cbs
# CBS-QB3 Test

CBS-QB3 on formaldehyde

0 1
分子说明

--Link1--
%Chk=cbs
%NoSave
# CBS-QB3(Restart,ReadIso) Geom=AllCheck Test

300.0 1.0
同位素说明

所以CBS-QB3不是放在基组那里的，直接摆在命令行的。
中文版g09说明书：
http://scc.ustc.edu.cn/zlsc/jsrj ... 410485230908642.pdf
60页，就是
赞一下 回复此楼
学无止境；天天学习；努力向上。
7楼2014-05-16 20:44:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amin0241

金虫 (小有名气)
引用回帖:
7楼: Originally posted by zhouwohua at 2014-05-16 20:44:19
不好意思，之前没细看。以为是guess=read。我在高斯09说明说上看了一下关于CBS方法的说明：
这个方法关键词指定使用由 Petersson 等人开发的各种完备基组（Complete Basis Set；
CBS）方法，用于进行非常精确的能 ...

这种写法计算H2一点问题没有，可是算I就不行了
赞一下 回复此楼
8楼2014-05-18 19:13:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amin0241

金虫 (小有名气)
引用回帖:
6楼: Originally posted by 枪下游魂 at 2014-05-16 14:44:05
你没懂我的意思：
1、我直接复制你给的输入文件去跑，高斯会报和你给出的一样的错误，原因就是你分子描述和附加输入部分中间空了太多行了，我删去多于的行数，不再报你贴出的这个错误，但是报赝势基组定义错误；
...

%chk=hi
%mem=96mw
%nproc=4
#P CCSD(T)/gen  pseudo=read

Title Card Required

0 1
I                  0.00000000    0.00000000    0.02990700
H                  0.00000000    0.00000000   -1.58505500

I 0
Lanl2dz
****
H  0
6-31+G(d')
****

I 0
Lanl2dz



--Link1--
%chk=hi
%mem=96mw
%nproc=4
#P MP4SDQ/gen pseudo=read Geom=AllCheck Guess=Read

I 0
Lanl2dz
****
H 0
CBSB4
****

I 0
Lanl2dz




--Link1--
%chk=hi
%mem=96mw
%nproc=4
#P Geom=AllCheck Guess=Read MP2/gen pseudo=read CBSExtrap=(NMin=10,MinPop)

I 0
Lanl2dz
****
H 0
CBSB3
****

I 0
Lanl2dz


我这样分布写时，前面计算都没有出问题，就是最后一步出错，出错信息：
  Complete Basis Set (CBS) Extrapolation:
M. R. Nyden and G. A. Petersson, JCP 75, 1843 (1981)
G. A. Petersson and M. A. Al-Laham, JCP 94, 6081 (1991)
G. A. Petersson, T. Tensfeldt, and J. A. Montgomery, JCP 94, 6091 (1991)
J. A. Montgomery, J. W. Ochterski, and G. A. Petersson, JCP 101, 5900 (1994)

Available Memory in CBS =  100663296
Minimum Number of PNO for Extrapolation =  10
Absolute Overlaps: IRadAn =     99590
LocTrn: ILocal=3 LocCor=F DoCore=T.
MOMinB can only project down in basis size.
Error termination via Lnk1e in d:\G09W\l904.exe at Fri Jan 18 16:18:35 2013.
Job cpu time:  0 days  0 hours  0 minutes 17.0 seconds.
File lengths (MBytes):  RWF=     14 Int=      0 D2E=      0 Chk=      1 Scr=      1
赞一下 回复此楼
9楼2014-05-18 19:17:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

syq.0522

银虫 (正式写手)
引用回帖:
9楼: Originally posted by amin0241 at 2014-05-18 19:17:53
%chk=hi
%mem=96mw
%nproc=4
#P CCSD(T)/gen  pseudo=read

Title Card Required

0 1
I                  0.00000000    0.00000000    0.02990700
H                  0.00000000    0.00000000   -1 ...

在B3lyp/6-311++G(2df,2pd)基础上，用CBS-QB3再次计算，可是CBS-QB3和B3LYP构型不一样。CBS-QB3不是相当于做单点么？
赞一下 回复此楼
为了生活而更好的生活！！
10楼2015-03-11 19:23:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 amin0241 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭