做GIC计算遇到L716报错问题“CrVal0: Division by zero detected!”

gingaorb

各位小哥哥，各位大佬们，我最近又遇到了一个比较特殊的报错问题，想请教一下大家呀~（这个问题我在Gaussian 常见报错及解决方法；新手求助报错时的注意事项 - 量子化学 (Quantum Chemistry) - 计算化学公社 (keinsci.com) 这个帖子里居然都没搜到）


问题是这样的，我想用GIC方法计算在气相中一个甲酸分子同时剥离出两个氢原子的能量曲线，参考了站长大大的文章：
http://sobereva.com/82 谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/474 详谈使用Gaussian做势能面扫描


以下是我输入文件的参数设置：
%chk=HCOOH-GIC.chk
%nprocshared=96
%mem=40GB
# opt=(modredundant, GIC, calcall) UPBE1PBE/def2tzvp empiricaldispersion=gd3bj pop=(always,nbo) guess=(mix,always) nosymm


Title Card Required


0 1
C                  0.00000000    0.41961800    0.00000000
O                  1.15928400    0.11679000    0.00000000
O                 -1.02924700   -0.44523700    0.00000000
H                 -0.39358100    1.44748200    0.00000000
H                 -0.64671500   -1.33761700    0.00000000
[size=18.6667px]

RCH(NSteps=50,StepSize=0.05)=R(1,4)
ROH=R(3,5)
F(Frozen)=RCH-ROH

[size=18.6667px]

然后，算到第41个点的时候，遇到报错：
Cartesian Forces:  Max     0.000675086 RMS     0.000304526
FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4.
CrVal0: IPrim =       6297
CrVal0: Division by zero detected!
Error termination via Lnk1e in /public/home/cjq/job_example/software/g16/l716.exe at Sun Apr 27 03:25:13 2025.
Job cpu time:       5 days  2 hours 33 minutes 57.3 seconds.
Elapsed time:       0 days  3 hours 52 minutes 40.1 seconds.
File lengths (MBytes):  RWF=     92 Int=      0 D2E=      0 Chk=      6 Scr=      2

第一次遇到这种特殊的情况诶~
可真是令人头疼啊~有没有有经验的小哥哥和大佬能给予我这个小白白一点点专业的指导和建议鸭？


希望大家能够不吝赐教啊！非常感谢大家的支持和帮助~