计算高斯出现NB too large for integers on this machine问题

Dragonmk3

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的晶胞大小

sobereva

Gaussian目前版本的PBC连谐振都做不了，还想做非谐振，根本不可能
而且6-31g这么垃圾的基组就算程序能给出结果又有何意义？哪怕你只是算非谐振校正量，这么烂的基组也绝对不可接受。太缺乏相关常识了
算晶胞尺寸直接对晶胞做优化就完了，根本不牵扯做振动分析

喵星大佬

优化为啥要算非偕振？

随便找个QE/VASP/CASTEP之类的算PBC不香么？金刚石一个元胞就两个原子

北大-陶豫

周期性体系的freq是和孤立体系很不一样的，涉及到声子什么的（随便找本固体物理书看看就知道了），高斯现在还不能实现，更不用说算非谐了。

Freeman

                               

高斯的PBC频率计算有bug，以下是我刚琢磨出的凑合的解决方法。
高斯只能算PBC的数值频率，不能算解析频率。如果你只写freq，它就会傻乎乎地当成解析频率来算（猜测高斯看到计算要求的是纯泛函，心想：这还不简单，上解析二阶导数！然后就崩了。这是第一个bug。），所以你要写freq=numer。数值频率慢如乌龟，几乎没什么用。
注意这样一来它可能仍然不能正常结束（第二个bug），也不能输出热力学数据，不过可以在报错退出前输出振动模式了。
有了振动模式，喂给GaussView，仍然不能显示振动的动画，vibration那个按钮是灰色的（第三个bug）。你要手动改一改输出文件，才能观看振动。以下是修改后的输出文件，能正常显示振动动画。你可以照着这个文件来改自己的输出文件。（本来想把修改前的输出文件放上来，可是太大了，不让上传）
2rev.log (230.71 KB, 下载次数 Times of downloads: 11)

2021-11-24 00:17 上传 Uploaded

点击下载
Click to download

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