高斯计算常见错误及解决方案.pdfVIP

⾼斯计算常见错误及解决⽅案

GAUSSION计算常见错误及解决⽅案

1.⾃旋多重度错误

2.变量赋值为整数

3.变量没有赋值

4.键⾓⼩于等于0度，⼤于等于180度

5.分⼦描述后⾯没有空⾏

6.⼆⾯⾓判断错误，造成两个原⼦距离过近

7.分⼦描述⼀⾏内两次参考同⼀原⼦,或参考原⼦共线

运⾏出错

1.⾃洽场不收敛SCF

a.修改坐标，使之合理

b.改变初始猜Guess

c.增加叠代次数SCFCYC=N

d.iop（5/13=1）

2.分⼦对称性改变

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

b.给出精确的、对称性确定的⾓度和⼆⾯⾓

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

d.不做对称性检查iop(2/16=1)

3.⽆法写⼤的Scratch⽂件RWF

a.劈裂RWF⽂件

%rwf=loc1,size1,loc2,size2,……..,locN,-1

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

c.限制最⼤硬盘maxdisk=NGB

4.FOPT出错原因是变量数与分⼦⾃由度数不相等。可⽤POPT或直接⽤OPT

5.优化过渡态只能做⼀个STEP原因是负本征数⽬不对添加iop（1/11）=1

6.组态相互作⽤计算中相关能叠代次数不够，增加叠代次数

QCISD（Maxcyc=N）

Default.Rou设置

在Scratch⽂件夹中的Default.Rou⽂件中设置G03程序运⾏的省缺参数：

-M-200MW

-P-4

-#-MaxDisk=10GB

-#-SCF=ConventionalorDirect

-#-MP2=NoDirectorDirect

-#-OPTCYC=200

-#-SCFCYC=200

-#-IOPs设置如iop（2/16=1）

Default.Rou设置中的冲突

Defaultroute:MaxDisk=2GBSCF=DirectMP2=DirectOPTCYC=200SCFcyc=100iop(2/16=1)iop(5/13=1)

#ccsd/6-31G**opt

L903/L905andL906canonlydoMP2.

问题在于，ＭP2=Direct！去掉这个设置，CCSD的作业就能进⾏了。

因此，建议在Default设置中只设置，内存，最⼤硬盘，等项

振动分析的关键词

振动分析的选项

拉曼光谱Raman

⾮谐性Anharmonic势能函数不是简谐的

内转动，受阻转⼦Hinderedrotor⼒常数过⼩的转动

振-转耦合VibRot

振动圆⼆⾊性VCD⽤于光活性分⼦

同位素效应，温度，压⼒ReadIsotopes不影响频率数值

算法选项：

解析频率Analytic

数值频率Numeric

⼆阶导数数值Enonly

FREQ的缺省设置是：有解析⽅式就计算解析频率，不计算Raman性质

MP3，MP4，CISD，CCSD，QCISD等⽅法没有解析频率，只能算数值频率。

⼀般各种理论⽅法下，计算的频率和实验测得频率之间有系统误差，需要⽤校正因⼦校正（ScaleFactor）

各种⽅法计算频率的Scale因⼦

⽅法：约化因⼦约化因⼦

（频率）（ZPE）

HF/3-21G0.90850.9409

HF/6-31G(d)0.89290.9135

MP2(Full)/6-31G(d)0.94270.9646

MP2(FC)/6-31G(d)0.94340.9676

BLYP/6-31G(d)0.99401.0119

B3LYP/6-31G(d)0.96130.9804

SVWN/6-31G(d)0.98331.0079

Freq=ReadIsotopes输⼊

#。。。Freq=ReadIsotopes。。。

Title

02

H

O10.98

298.151.50.8929温度，压⼒约化因⼦

H2⽤质量数是2（程序读取精确值）的同位素

O18⽤质量数是18（程序读取精确值）的同位素

空⾏

Freq=（ReadFC，ReadIsotopes）

同⼀⽅法基组下计算不同同位素的⼒常数矩阵（Hessian）相同，因此从前⾯频率计算的chechpoint⽂件⾥读取⼒常数，计算

同位素效应的频率和热⼒学性质，可以节省许多时间。注：同种元素的不同同位素