计算化学公社

标题: Gaussian计算出现Allocation failure for numerical quadrature问题 [打印本页]

作者
Author: gym 时间: 2018-1-29 11:26
标题: Gaussian计算出现Allocation failure for numerical quadrature问题
我的输入文件：
%chk=ZnPcPDI4_model.chk
%mem=32GB
%nprocs=24
# opt b3lyp/genecp integral=(finegrid,acc2e=10)

ZnPcPDI4_model

0 1
出现如下错误：
ECPInt:  NShTT= 1023165 NPrTT= 6494122 LenC2=  188380 LenP2D=  525616.
LDataN:  DoStor=T MaxTD1= 6 Len=  172
Calling FoFJK, ICntrl=    2127 FMM=T ISym2X=0 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0.
In Quad4, NSgAt2=    417 but IDerNu=1 NAtomA=    465 NAtAll=    401.
Allocation failure for numerical quadrature.
Error termination via Lnk1e in /share/home/majing/g16/l703.exe at Fri Jan 26 04:58:36 2018.
Job cpu time:    8 days 17 hours 36 minutes 27.1 seconds.
Elapsed time:    0 days 12 hours 41 minutes 40.0 seconds.
File lengths (MBytes):  RWF= 9360 Int=    0 D2E=    0 Chk= 1184 Scr=    32

请求各位大神，帮助解惑。

作者
Author: 我本是个娃娃 时间: 2018-1-29 11:40
泛函用MN15L比较好。【计算过渡金属配合物：M06、BP86、B3LYP、MN15L（哪个更好不一定，BP86有时很不错）】如果用G16，MN15L测试的结果比较好。

积分格点用int=finegrid。你的积分格点有点低。

PS：你是南*大学的吧？

作者
Author: gym 时间: 2018-1-29 12:04

我本是个娃娃发表于 2018-1-29 11:40
泛函用MN15L比较好。【计算过渡金属配合物：M06、BP86、B3LYP、MN15L（哪个更好不一定，BP86有时很不错）】 ...

你好，确实是**大学。不好意思，我输入文件粘贴的时候错了，是integral=(finegrid,acc2e=10)，还是出现这种错误，请问这种错误出现，是什么原因呢？

作者
Author: 我本是个娃娃 时间: 2018-1-29 12:09

gym 发表于 2018-1-29 12:04
你好，确实是**大学。不好意思，我输入文件粘贴的时候错了，是integral=(finegrid,acc2e=10)，还是出现 ...

用MN15L泛函试试

尝试用G16的默认格点设置。

作者
Author: gym 时间: 2018-1-29 12:14

我本是个娃娃发表于 2018-1-29 12:09
用MN15L泛函试试

尝试用G16的默认格点设置。

您好，请问您知道怎么降低g16的收敛标准吗？

作者
Author: zjxitcc 时间: 2018-1-29 12:37
SCF收敛标准:scf(conver=N)，默认是8，不建议低于7.
opt结构收敛标准:opt=xxx。其中opt=loose选项比较粗糙，适用于预优化，成功后再去掉loose优化至更严格的收敛标准。可以参看http://gaussian.com/opt/的说明。
g16的SCF收敛标准在conver=7以下时，它是不允许做opt的，会警告用户收敛限太低；实际上这也比较合理，conver小于等于6时，算单点勉强可以，做opt是过于粗糙的。
另外积分格点coarsegrid也过于粗糙，作为预优化是可以的，不适合做正式计算。
这错误比较少见，可以尝试换个泛函；如果实在想用这个泛函，可以先在其他泛函下优化好了，再改用这个泛函优化。
另外，如果是优化了几步之后出现这个问题，可以试试在一开始加上opt=calcfc；如果是第一个结构就出现这种问题，calcfc就可能没啥用了。

作者
Author: liyuanhe211 时间: 2018-1-29 13:23
检查一下计算的时候有没有那么多内存，例如你声明允许高斯用40G内存，但系统实际只剩30G这种情况

作者
Author: gym 时间: 2018-1-29 14:23

zjxitcc 发表于 2018-1-29 12:37
SCF收敛标准:scf(conver=N)，默认是8，不建议低于7.
opt结构收敛标准:opt=xxx。其中opt=loose选项比较粗糙 ...

您好，您说的g16的scf在conver=7以下，无法做opt，是在guassian16手册的哪里，能否发一下链接，谢谢

作者
Author: zjxitcc 时间: 2018-1-29 15:02

gym 发表于 2018-1-29 14:23
您好，您说的g16的scf在conver=7以下，无法做opt，是在guassian16手册的哪里，能否发一下链接，谢谢

不需要看手册，直接运行一个甲烷的结构优化，设置scf(conver=6)，报错之后打开log文件最后一百行就看到了

作者
Author: sobereva 时间: 2018-1-29 15:10
典型的乱用关键词(coarsegrid)自找麻烦

作者
Author: gym 时间: 2018-1-29 15:24

sobereva 发表于 2018-1-29 15:10
典型的乱用关键词(coarsegrid)自找麻烦

sob老师，我是用了finegrid，而coarsegrid是我问问题是打错了，那原因是什么，恳请老师解答？

作者
Author: sobereva 时间: 2018-1-29 16:03

gym 发表于 2018-1-29 15:24
sob老师，我是用了finegrid，而coarsegrid是我问问题是打错了，那原因是什么，恳请老师解答？

上传输入输出文件，否则不好说
如果确认输入文件肯定没问题，尝试G09也跑一下，没准是G16 bug

作者
Author: gym 时间: 2018-1-29 16:27

sobereva 发表于 2018-1-29 16:03
上传输入输出文件，否则不好说
如果确认输入文件肯定没问题，尝试G09也跑一下，没准是G16 bug

sob老师，我已经上传了文件，在g09上运行也会出现这样的情况。

作者
Author: sobereva 时间: 2018-1-29 16:45

gym 发表于 2018-1-29 16:27
sob老师，我已经上传了文件，在g09上运行也会出现这样的情况。

用ultrafine再试，还不行尝试换个基组再试，还不行给官方客服发邮件反馈bug
如置顶新人必读贴所示，上传较大文本文件前一定要压缩。输出文件已删除

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