常见错误分析

常见错误分析

1．检查是否有初始文件错误

在命令行中加入%kJob L301 or %kJob L302

如果通过则一般初始文件ok。

常见初级错误：

a. 自旋多重度错误

b. 变量赋值为整数

c. 变量没有赋值或多重赋值

d. 键角小于等于0度，大于等于180度

e. 分子描述后面没有空行

f. 二面角判断错误，造成两个原子距离过近

g. 分子描述一行内两次参考同一原子,或参考原子共线

2．SCF(自洽场)不收敛

则一般是L502错误

省却情况做64个cycle迭代(G03缺省128 cycles)

a. 修改坐标，使之合理

b. 改变初始猜 Guess=Huckel或其他的，看Guess关键词。

c. 增加叠代次数SCFCYC=N （对小分子作计算时最好不要增加，很可能结构不合理）

d. iop（5/13=1）这样忽略不收敛，继续往下做。

3．分子对称性改变

a. 修改坐标，强制高对称性或放松对称性

b. 给出精确的、对称性确定的角度和二面角。如CH4的角度给到109.47122

c. 放松对称性判据 Symm=loose

d. 不做对称性检查iop(2/16=1) （最好加这个选项）

iop(2/16=2) 则保持新的对称性来计算

4．Opt时收敛的问题

a. 修改坐标，使之合理

b. 增加叠代次数optcyc=N

5．优化过渡态，若势能面太平缓，则不好找到。

iop(1/8=10) 默认30（下一个结构和该结构的差别0.3牛筛某10。如果每一步都要用到小的步长，应该加opt(notrustupdate)

6．在CI（组态）方法中

如QCISD(T),CCSD(T),CID方法中，省却最大循环50，若出错（L913错误）

解决方法：#P QCISD(maxcyc=N) 注：N≤512

7．优化过渡态

opt=TS (给出过渡态)

opt=qst2 (给出反应物和产物)

opt=qst3 (给出反应物和产物和过渡态)

a. 用G03时的出错 opt=ts必须加FC (force constant)

写法：opt=(TS, calcFc)

or opt=(TS,calchffc)计算HF力常数，对QCISD，CCSD等方法用；

or opt=(TS,modRedundant) (最好写这个)

b. 如果计算采用QCISD计算（不好计算FC）

则写为QCISD opt=(TS, calcHFFC) (用HF计算FC)

8. 无法写大的Scratch文件RWF

a. 劈裂RWF文件%rwf=loc1,size1,loc2,size2,……..,locN,-1

b. 改变计算方法MP2=Direct可以少占硬盘空间

c. 限制最大硬盘maxdisk=N GB，****MB,有些系统写2GB会出错，可以写2000MB

9. FOPT出错原因是变量数与分子自由度数不相等。可用POPT 或直接用OPT

10. 优化过渡态只能做一个STEP 原因是负本征数目不对添加iop（1/11）=1或者noeigentest 一

“G98的出错讯息都在用户手册中的…Gaussian 98的链接‟一节中就给出了，常见错误有：1。输入有误，这没话说了。一般都是L1，L101出错

2。基组有误。特别是刚开始用自定义基组的，或者不知道基组适合范围的，这是L301出错

3。无法收敛、优化和频率计算出错，这就很多了，比如L502，L716，L9999。这就只好具体问题具体分析了。”

也就是说，当计算非正常退出时，首先查看文件结束处的说明，比如类似的句子：

Error termination via Lnk1e in /para/soft/g03/l502.exe at Mon Mar 29 16:31:18 2004.

这说明是L502错误，可以在g98 help中的Links，具体看看L502是干吗的，对错误有个初步了解，再查找相应的关键字

二

出错：

我在用OPT的一个参数STEPSIZE的时候，从来没有成功过，不知道是为什么？

# UHF/6-31 OPT(Stepsize=30)

---------------------------

QPERR --- A SYNTAX ERROR WAS DETECTED IN THE INPUT LINE.

# UHF/6-31 OPT(STEPSIZE=30) '

Last state="OPT2"

TCursr= 9023 LCursr= 15

Error termination via Lnk1e in D:\G98W\l1.exe.

Job cpu time: 0 days 0 hours 0 minutes 0.0 seconds.

File lengths (MBytes): RWF= 6 Int= 0 D2E= 0 Chk= 5 Scr= 1

解决：

STEPSIZE=30 是g98的一个bug。可以用iop(1/8=30)来实现。

三

出错信息：

输入：B3LYP/gen( SDD for U and 6-311+G(3df) for O ) opt 输出：

The electronic state of the initial guess is 3-SGG.

of initial guess= 2.0000

Requested convergence on RMS density matrix=1.00D-04 within 200 cycles. Requested convergence on MAX density matrix=1.00D-02.

Requested convergence on energy=1.00D-02.

No special actions if energy rises.

Defaulting to unpruned grid for atomic number 92.

Spurious integrated density or basis function:

NE= 48 NElCor= 0 El error=2.25D-04 rel=4.68D-06 Tolerance=1.00D-03

Shell 25 absolute error=1.91D-02 Tolerance=1.20D-02

Shell 12 signed error=9.34D-05 Tolerance=1.00D-01

Inaccurate quadrature in CalDSu.

Error termination via Lnk1e in /para/soft/g03/l502.exe at Mon Mar 29 16:31:18 2004. Job cpu time: 0 days 0 hours 1 minutes 6.2 seconds.

File lengths (MBytes): RWF= 24 Int= 0 D2E= 0 Chk= 1 Scr= 1 解决：

猜想是增加积分网格。

四

出错：

Problem detected with inexpensive integrals.

Switching to full accuracy and repeating last cycle.

Cycle 2 Pass 1 IDiag 1:

Warning! Spurious integrated density:

NE= 188 NElCor= 0 Integral= -76.06821 Tolerance=1.00D-03

Consistency failure #2 in CalDSu.

Error termination via Lnk1e in /home/gaussian/g98/l502.exe.

Job cpu time: 1 days 3 hours 21 minutes 30.2 seconds.