计算化学公社

标题: 计算高斯出现NB too large for integers on this machine问题 [打印本页]

作者
Author:
Dragonmk3    时间: 2021-9-15 10:25
标题: 计算高斯出现NB too large for integers on this machine问题
log如下:
**** Warning!!: The largest alpha MO coefficient is  0.22162853D+02

Symmetrizing basis deriv contribution to polar:
IMax=3 JMax=2 DiffMx= 0.00D+00
G2DrvN: will do     9 centers at a time, making    1 passes.
Calling FoFCou, ICntrl=  3107 FMM=T I1Cent=   0 AccDes= 0.00D+00.
NB too large for integers on this machine.
Error termination via Lnk1e in G:\G16W\l1110.exe at Wed Sep 15 01:52:57 2021.
Job cpu time:       0 days  2 hours 36 minutes 50.0 seconds.
Elapsed time:       0 days  2 hours 36 minutes 39.8 seconds.
File lengths (MBytes):  RWF=  48040 Int=      0 D2E=      0 Chk=   4554 Scr=      1



我的gjf如下:
%mem=20GB
%chk=G:\G16W\e.chk
%CPU=0,1,2,3,4,5,6,7
# opt freq=anharmonic b3lyp/6-31g

Title Card Required

0 1
C                  0.00000000    0.00000000    0.00000000
C                  0.89000000    0.89000000    0.89000000
C                  0.00000000    1.78000000    1.78000000
C                  2.67000000    0.89000000    2.67000000
C                  1.78000000    0.00000000    1.78000000
C                  0.89000000    2.67000000    2.67000000
C                  1.78000000    1.78000000    0.00000000
C                  2.67000000    2.67000000    0.89000000
Tv                 3.56000000    0.00000000    0.00000000
Tv                 0.00000000    3.56000000    0.00000000
Tv                 0.00000000    0.00000000    3.56000000


我的完整最近的完整log如下:
https://paste.ubuntu.com/p/SnsWvDR6fq/

我的需求如下:
计算金刚石结构的C-14的结构优化和非谐振的计算,通过这样来计算C14的晶胞大小
作者
Author:
sobereva    时间: 2021-9-15 15:28
Gaussian目前版本的PBC连谐振都做不了,还想做非谐振,根本不可能
而且6-31g这么垃圾的基组就算程序能给出结果又有何意义?哪怕你只是算非谐振校正量,这么烂的基组也绝对不可接受。太缺乏相关常识了
算晶胞尺寸直接对晶胞做优化就完了,根本不牵扯做振动分析

作者
Author:
喵星大佬    时间: 2021-9-15 23:09
优化为啥要算非偕振?

随便找个QE/VASP/CASTEP之类的算PBC不香么?金刚石一个元胞就两个原子
作者
Author:
北大-陶豫    时间: 2021-9-16 11:27
周期性体系的freq是和孤立体系很不一样的,涉及到声子什么的(随便找本固体物理书看看就知道了),高斯现在还不能实现,更不用说算非谐了。
作者
Author:
Freeman    时间: 2021-11-24 00:21
                               
高斯的PBC频率计算有bug,以下是我刚琢磨出的凑合的解决方法。
高斯只能算PBC的数值频率,不能算解析频率。如果你只写freq,它就会傻乎乎地当成解析频率来算(猜测高斯看到计算要求的是纯泛函,心想:这还不简单,上解析二阶导数!然后就崩了。这是第一个bug。),所以你要写freq=numer。数值频率慢如乌龟,几乎没什么用。
注意这样一来它可能仍然不能正常结束(第二个bug),也不能输出热力学数据,不过可以在报错退出前输出振动模式了。
有了振动模式,喂给GaussView,仍然不能显示振动的动画,vibration那个按钮是灰色的(第三个bug)。你要手动改一改输出文件,才能观看振动。以下是修改后的输出文件,能正常显示振动动画。你可以照着这个文件来改自己的输出文件。(本来想把修改前的输出文件放上来,可是太大了,不让上传)
(, 下载次数 Times of downloads: 11)

因为有这些bug,最好不要用高斯算PBC了。不过GaussView的PBC建模功能很给力。这里带点儿私货,推荐一下我弄的GaussView+PWscf的解决方案http://bbs.keinsci.com/thread-20166-1-1.html






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