求助用B2PLYPD3/def2TZVP计算单点时报错

yuanbingling

老师好，我用B2PLYPD3/def2TZVP计算单点输出文件末尾出现以下内容：
DoSCS=T DFT=T ScalE2(SS,OS)=  0.270000  0.270000
Fully direct method using O(OVN) memory.
Internal consistency failure #1 in GetIJB.
Error termination via Lnk1e in l906.exe
查询时说可能是程序BUG，但不知道怎么解决。gaussian只有D.01版本的。还请老师指点

输入文件如下：
%chk=TS2optsp.chk
%nprocshared=24
%mem=48GB
# B2PLYPD3/def2TZVP scrf=(smd,read)

Title Card Required

0 1
N                 -0.38856300    2.95830600   -0.78842000
H                 -0.13945300    2.54228400   -1.68463800
H                 -0.93672400    3.79297600   -0.99342100
C                 -1.18449400    2.00600900    0.01434000
H                 -0.56226300    1.12493600    0.19291000
H                 -1.38954400    2.46951900    0.98208000
C                 -2.48581300    1.60542900   -0.66808200
H                 -3.11383300    2.49337000   -0.79564300
H                 -2.25868000    1.22482500   -1.67818600
N                 -3.22294000    0.62825300    0.13474100
H                 -2.64259500   -0.20425100    0.22116000
C                 -4.47882400    0.23419700   -0.51451000
H                 -5.11672700    1.12233000   -0.58587600
H                 -4.30914500   -0.12110600   -1.54310900
C                 -5.20161700   -0.84717100    0.27571200
H                 -5.37397000   -0.49365800    1.30485800
H                 -4.56397400   -1.73554200    0.34711000
N                 -6.45311800   -1.23561800   -0.38276200
H                 -7.03452100   -0.40333400   -0.46123200
C                 -7.19822700   -2.23169600    0.39173400
H                 -7.41564900   -1.87664900    1.41425200
H                 -6.57796200   -3.12777100    0.49323800
C                 -8.50566100   -2.59665400   -0.29578100
H                 -8.29052900   -3.05410900   -1.26520500
H                 -9.07399500   -1.67496300   -0.49252100
N                 -9.26342300   -3.56466400    0.51303000
H                 -9.58375800   -3.11339400    1.36693100
H                -10.10264000   -3.83977100    0.00983900
N                  1.78304000    3.69192900    0.46709300
H                  0.77209300    3.34201900   -0.10467000
H                  1.63808500    3.51830200    1.46228000
C                  2.92068700    2.88982900   -0.05060500
H                  2.88961700    2.93592900   -1.14048400
H                  3.85878600    3.33448900    0.28502900
C                  2.80059500    1.45215900    0.43485400
H                  2.87891800    1.43176100    1.52679000
H                  1.80549200    1.05720800    0.17767600
N                  3.87930800    0.63385000   -0.11945400
H                  3.76258400    0.60570300   -1.13075600
C                  3.82641100   -0.74392900    0.38117200
H                  4.00772500   -0.72034500    1.46158200
H                  2.83316100   -1.19629800    0.23273500
C                  4.87741400   -1.61662300   -0.28854700
H                  5.86701600   -1.14451200   -0.17963100
H                  4.67097700   -1.67803300   -1.36286900
N                  4.85593800   -2.97417600    0.26256800
H                  5.06951800   -2.90853400    1.25603500
C                  5.85927000   -3.84348400   -0.35700700
H                  6.86981700   -3.40154400   -0.30755100
H                  5.61462200   -3.95790300   -1.41774700
C                  5.88679400   -5.21167900    0.30784900
H                  4.91669600   -5.70030300    0.18103500
H                  6.04174300   -5.07810800    1.38917800
N                  6.91524700   -6.06456000   -0.30909300
H                  7.83738500   -5.69316200   -0.09232800
H                  6.88066500   -6.99083100    0.10787300
C                  1.85660100    5.17286500    0.26369700
O                  1.20814500    5.82590200    1.10337700
O                  2.50358300    5.54781200   -0.72707600

eps=54.1915
epsinf=2.0805

我本是个娃娃

#p B2PLYP/def2TZVP scrf=(smd,read) em=gd3

yuanbingling

我本是个娃娃 发表于 2020-6-4 10:14
#p B2PLYP/def2TZVP scrf=(smd,read) em=gd3

好的，谢谢老师

snljty

B2PLYPD3手册上说是B2PLYP加上BJ阻尼的D3，我觉得可以写B2PLYP EmpiricalDispersion=GD3BJ。另外你用双杂化泛函结合隐式溶剂模型要考虑一下有没有必要算外迭代（虽然昂贵不少）。

yuanbingling

snljty 发表于 2020-6-4 11:27
B2PLYPD3应该是G16才引入的关键词，手册上说是B2PLYP加上BJ阻尼的D3，我觉得可以写B2PLYP EmpiricalDispers ...

试了 EmpiricalDispersion=GD3BJ和em=gd3 的关键词，但是还是出现了相同的问题，前面也用同样的关键词算了其他的结构是没有出现问题的，但是这个分子不知道为什么

niobium

这个报错可以通过减少核数解决。另外，用了fulldir的话尽量给多些内存，你48G内存恐怕要算10天8天的

sobereva

snljty 发表于 2020-6-4 11:27
B2PLYPD3应该是G16才引入的关键词，手册上说是B2PLYP加上BJ阻尼的D3，我觉得可以写B2PLYP EmpiricalDispers ...

B2PLYPD3关键词从G09 D.01就有，当前不是这个问题。一般我都这么写，省事。

尝试给大内存。实在算不了的话用ORCA，4核跑这个都没压力

Dolantin

niobium 发表于 2020-6-4 15:12
这个报错可以通过减少核数解决。另外，用了fulldir的话尽量给多些内存，你48G内存恐怕要算10天8天的

试了一下减少核数，确实可行，四个单点任务顺利结束，感谢！