计算化学公社

标题: G2算到最后一步出错 [打印本页]

作者
Author:
mfdsrax2    时间: 2022-2-1 15:00
标题: G2算到最后一步出错
G2算了很久,终于到最后一步
#N Geom=AllCheck Guess=TCheck SCRF=Check MP2/6-311+G(3df,2p)
但是刚一开始就出错

**** Warning!!: The largest alpha MO coefficient is  0.12072438D+03

Fully direct method using O(OVN) memory.
Internal consistency failure #1 in GetIJB.
Error termination via Lnk1e in *** at ***.


使用的软件是 Gaussian 09, Revision D.01
看论坛上说好像是一个bug,但是算了好几个不同的分子都有类似的问题,这是怎么回事



作者
Author:
alonewolfyang    时间: 2022-2-1 18:40
分子太大了?
作者
Author:
mfdsrax2    时间: 2022-2-1 19:09
alonewolfyang 发表于 2022-2-1 18:40
分子太大了?

24核48GB内存算了16天的样子,这个有关系吗
作者
Author:
alonewolfyang    时间: 2022-2-1 19:33
mfdsrax2 发表于 2022-2-1 19:09
24核48GB内存算了16天的样子,这个有关系吗

老师说过,MP2/6-311+G(3df,2p)这一步很耗内存,况且考虑其他步骤计算级别,一般也不超过30个原子
作者
Author:
mfdsrax2    时间: 2022-2-1 21:03
alonewolfyang 发表于 2022-2-1 19:33
老师说过,MP2/6-311+G(3df,2p)这一步很耗内存,况且考虑其他步骤计算级别,一般也不超过30个原子

很专业啊,但是这个课题已经基本算完了,就差这两个了,也没时间重算了

分子有18个CNO原子,含两个硝基,按理说不算大,不知道为何算的很慢,我再加内存试试
作者
Author:
sobereva    时间: 2022-2-1 21:26
这年头还用极其过时的G2很容易被comment,又贵精度又落伍
做后HF计算,内存能给多大给多大,绝对不要吝啬,否则轻则多花耗时,重则可能根本没法算
如果你的这个问题始终无法解决,尝试G16,或者干脆换其它方法重算所有的


作者
Author:
mfdsrax2    时间: 2022-2-1 21:40
sobereva 发表于 2022-2-1 21:26
这年头还用极其过时的G2很容易被comment,又贵精度又落伍
做后HF计算,内存能给多大给多大,绝对不要吝啬 ...

我知道G2已经很low了,有什么好的替换方案吗?

带MP2的看起来算得很快,这个精度影响大吗?
作者
Author:
sobereva    时间: 2022-2-1 22:23
mfdsrax2 发表于 2022-2-1 21:40
我知道G2已经很low了,有什么好的替换方案吗?

带MP2的看起来算得很快,这个精度影响大吗?

4L贴的我讲的北京科音基础量化班的ppt里已经充分体现了
用G4(MP2)

作者
Author:
alonewolfyang    时间: 2022-2-4 15:44
本帖最后由 alonewolfyang 于 2022-2-4 15:46 编辑
mfdsrax2 发表于 2022-2-1 21:40
我知道G2已经很low了,有什么好的替换方案吗?

带MP2的看起来算得很快,这个精度影响大吗?

你做含能生成焓计算,别被所谓某些含能主流计算方法误导,用Sob老师给的定义法计算不香吗,即便你用G2算完了,你还得按照他们的设计等键反应再反求你化合物的生成焓,不嫌麻烦吗,虽然你主合成,计算辅助,所以应该用简便快速有效的方法算出来,有个数就行了,个人见解,不喜勿喷
作者
Author:
mfdsrax2    时间: 2022-2-4 16:25
alonewolfyang 发表于 2022-2-4 15:44
你做含能生成焓计算,别被所谓某些含能主流计算方法误导,用Sob老师给的定义法计算不香吗,即便你用G2算 ...

这个是直接拿整个分子计算的,不用设计等键反应,所以说算的很慢
作者
Author:
alonewolfyang    时间: 2022-2-4 16:56
mfdsrax2 发表于 2022-2-4 16:25
这个是直接拿整个分子计算的,不用设计等键反应,所以说算的很慢

我说的不是你这里的意思
作者
Author:
mfdsrax2    时间: 2022-2-4 17:45
alonewolfyang 发表于 2022-2-4 16:56
我说的不是你这里的意思

定义法是怎么算的,能发一下吗
作者
Author:
alonewolfyang    时间: 2022-2-7 00:03
mfdsrax2 发表于 2022-2-4 17:45
定义法是怎么算的,能发一下吗

我意思是不论你怎么求,都是算的绝对焓,而你目的应该是计算焓变




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3