Gaussian 计算CCSD(T)出错

hee_csu

请教大神:
我在采用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



请高手指点！非常感谢！

yflchx

这个分子不小了，硬盘不够？

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

int=ultrafine画蛇添足。

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

sobereva

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

niobium

硬要算的话搞个大内存机器用fulldirect

hee_csu

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

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

非常感谢指导！

hee_csu

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

感谢大神的指点

hee_csu

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

谢谢

hee_csu

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

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

milkxx

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

sobereva

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

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

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

yflchx

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拿来算单点应该比较准确。

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