计算化学公社

标题: Gaussian 计算CCSD(T)出错 [打印本页]

作者
Author: hee_csu 时间: 2017-5-10 16:12
标题: Gaussian 计算CCSD(T)出错
请教大神:
我在采用CCSD(T)计算小分子能量时，出错，最后十行内容如下：

LnScr1=165548924416 LExtra=          0 Total= 798228789673
MaxDsk=       -1 SrtSym=          T ITran=          4
DoSDTr:  NPSUse= 16
JobTyp=1 Pass  1:  I= 1 to  54.
Complete sort for first half transformation.
First half transformation complete.
Erroneous write. Write 5046272 instead of 8384512.
fd = 4
orig len = 8384512 left = 8384512
g_write

输入文件如下：

%mem=48GB
%nprocs=16
%chk=INT1.chk
#p ccsd(t)/aug-cc-pvtz
gfinput iop(5/13=1)
scf=tight integral=ultrafine

INT1

0 2
C                -2.45215400 -0.08445700 0.44267400
C                -1.65044400 -1.27332400 0.86901900
C                -0.26299600 -1.42870600 0.19529500
C                0.80700700 -0.40510700 0.66723900
H                -2.55572500 0.76161300 1.11795300
H                -2.21577000 -2.19044400 0.64232800
H                -1.49225600 -1.25493900 1.95107100
H                0.12473600 -2.42611700 0.43348300
H                -0.37685000 -1.38442400 -0.89433500
C                -3.16843000 -0.01492200 -0.86825800
H                -2.95254100 0.93368700 -1.38212500
H                -2.83705300 -0.82030800 -1.53534600
C                -4.71962900 -0.09522100 -0.73767300
O                0.83551500 -0.42236500 2.08398800
H                0.60576300 0.49205200 2.35875700
C                2.20954200 -0.77271100 0.13397300
C                3.33086400 -0.47297200 0.91890800
C                2.41522000 -1.36905700 -1.11802800
C                4.61870200 -0.75427400 0.46501300
H                3.17716300 -0.03643900 1.89915500
C                3.70479900 -1.65270700 -1.57377900
H                1.57163100 -1.63344500 -1.74847600
C                4.81272400 -1.34438600 -0.78536700
H                5.47321900 -0.51647000 1.09358200
H                3.83883000 -2.12203700 -2.54513100
H                5.81612100 -1.56746000 -1.13825700
C                0.46936100 1.03013400 0.19750500
C                0.39758500 1.40175900 -1.15190600
C                -0.00476500 3.18782400 0.87273000
C                0.09963500 2.72251400 -1.46943900
H                0.58767400 0.67440000 -1.93366100
C                -0.11067000 3.64188700 -0.43776700
H                -0.14817400 3.86417400 1.71302700
H                0.03972500 3.03553500 -2.50846100
H                -0.34048900 4.68251100 -0.64389100
N                0.27755600 1.91531200 1.18241700
H                -5.19076900 0.05892900 -1.71695100
H                -5.07929900 0.70405800 -0.07861200
C                -5.17425900 -1.38185200 -0.20253300
N                -5.51507300 -2.41280800 0.20877500

请高手指点！非常感谢！

作者
Author: yflchx 时间: 2017-5-10 16:43
这个分子不小了，硬盘不够？

iop(5/13=1)永远不要用。

int=ultrafine画蛇添足。

实在不行，试试：CCSD(T)/aug-cc-pVTZ ≈ CCSD(T)/aug-cc-pVDZ + MP2/aug-cc-pVTZ – MP2/aug-cc-pVDZ，没办法的办法，论文里要说清楚。
搜一下，应该有类似用法的参考文献。

作者
Author: sobereva 时间: 2017-5-10 19:14
你这种体系用ccsd(t)/aug-cc-pvtz根本是不可能的，光是cc-pVTZ都算不动。以你的计算条件，体系无对称性的情况下，用cc-pVTZ撑死了能算20个原子的有机分子。
这些都毫无意义还添乱：
gfinput iop(5/13=1)
scf=tight integral=ultrafine
看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331
你这种体系用双杂化泛函算单点还差不多。

作者
Author: niobium 时间: 2017-5-10 22:24
硬要算的话搞个大内存机器用fulldirect

作者
Author: hee_csu 时间: 2017-5-15 13:32

yflchx 发表于 2017-5-10 16:43
这个分子不小了，硬盘不够？

iop(5/13=1)永远不要用。

非常感谢指导！

作者
Author: hee_csu 时间: 2017-5-15 13:32

sobereva 发表于 2017-5-10 19:14
你这种体系用ccsd(t)/aug-cc-pvtz根本是不可能的，光是cc-pVTZ都算不动。以你的计算条件，体系无对称性的情 ...

感谢大神的指点

作者
Author: hee_csu 时间: 2017-5-15 13:33

niobium 发表于 2017-5-10 22:24
硬要算的话搞个大内存机器用fulldirect

谢谢

作者
Author: hee_csu 时间: 2017-5-15 13:34

sobereva 发表于 2017-5-10 19:14
你这种体系用ccsd(t)/aug-cc-pvtz根本是不可能的，光是cc-pVTZ都算不动。以你的计算条件，体系无对称性的情 ...

这种体系，用M06\6-31G(d)优化，用M06/6-311++G(d,p)算单点能，是否可以接受？

作者
Author: milkxx 时间: 2017-5-15 17:48
如果想要高精度的话可以考虑 MP2 方法 + delta（CCSD(T)(小基组)-MP2(小基组)）的近似算法，该体系优化不添加弥散函数应该是可以接受的

作者
Author: sobereva 时间: 2017-5-16 15:48

hee_csu 发表于 2017-5-15 13:34
这种体系，用M06\6-31G(d)优化，用M06/6-311++G(d,p)算单点能，是否可以接受？

何不用双杂化，这个大小完全算得动。
如果非要用普通泛函，用def2-TZVP比你这个基组强得多。你这种体系完全不需要弥散函数。

至于前面说的借助MP2来近似得到CCSD(T)大基组下的结果，虽然可以试，但精度未必比用最好的双杂化好，另外额外需要算的 CCSD(T)/小基组对于这种大小的体系也依然是很耗时的。

作者
Author: yflchx 时间: 2017-5-16 20:11

hee_csu 发表于 2017-5-15 13:34
这种体系，用M06\6-31G(d)优化，用M06/6-311++G(d,p)算单点能，是否可以接受？

我印象中M06-2X更好一点。用M0x系列，int=ultrafine应该是标配。

++不需要，想用的话，一个+就可以。用M0x，Truhlar文章中常用的MG3S拿来算单点应该比较准确。

看看别人如何处理类似的反应，这样有助于形成你自己的判断。时间充裕的话，多尝试。

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