gaussian基本概念和用法

1 计 算 流 程

上面的图应该是从L301说起的，即图中的“基组”，L301的作用是产生基组信息。（在这之前还有L1：处理计算执行路径，创建执行链接的列表，并初始化scratch文件；L101：读取标题和分子说明部分；L103：berny优化到最小值；L202：重新定位坐标，计算对称性，检查变量）“基组”之后是“分子结构”即L401，形成初始的分子轨道初猜。第一个框中有“分子轨道”和“能量”两项，我想这应该指的是L502即迭代求解SCF方程。下面红线圈的框中则是迭代求解的具体过程。每次迭代完毕都得到一个能量。即（SCF Done: E(RHF) = -xx.xxxx A.U. after xx cycles）

然后是L103判断力和位移是否满足收敛条件。如果满足（即四个yes），则进行布居分析等，并完成计算。如果不满足，则继续调整分子结构（L401），再次进行迭代求解，反复这个过程，直到满足力和位移收敛条件为止。如果不是HF方法而是多体微扰，CI等后-scf方法的话，在每次迭代（L502）完成后还要多一步相关能的计算，以MP2为例，迭代完成后多出:L801(双电子积分变换的初始化),L906(半直接的MP2)和L1002(迭代求解CPHF方程)等计算。

对于楼主的图，我觉得力和位移判断为NO以后，箭头应该指向最头上的“基组”，因为实际过程中，每次构型调整都会重新定位坐标和产生基组信息。

常见问题分析

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 （对小分子作

1

计算时最好不要增加，很可能结构不合理）; 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

一些容易忽略的地方

1。低频率振动模式时，力常数很小，必须使用opt＝tight以确保适当的收敛和随后任务步骤中频率计算的可靠度。我曾经算过一个体系，第一个频率振动的力常数只有10～20，很不放心，加入tight后，优出来的结构却与前面大相径庭，而且还出现了一个虚频。 2。做IRC计算的时候，如果步长很小，必须要用verytight。就曾有人步长设为1，而不加verytight，算出来的反应路径是v字形倒海鸥状的。

3。使用int指定用于数值积分的积分网格来用于消虚频的时候，必须注意，在比较能量的时候，对所有的计算要使用相同的积分网格。低频振动模式就是力常数比较小，至于有多大算低频振动模式，我并不太清楚。 不过根据我的经验，>50的时候不用加tight，<20的时候就要小心了。这也跟你的体系有关，我做的是易变的自由基体系，这方面要多留个心眼了。不知你算的是什么分子？我在优化构型经常遇到需消虚频.请指教

不过象我上次提出的有关中间体的虚频问题，虚频是只有10-20啊，可是就是消不掉啊。不过它是一个整体晃动的情况。我也是用了opt＝tight试过的

Gaussian计算出错和对应的解决方法

supi指出：“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)

2

Last state=\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.

File lengths (MBytes): RWF= 118 Int= 0 D2E= 0 Chk= 17 Scr= 1 解决：(from ccl)

As far as I know it's a dft-related error. Gaussian checks the charge resulting from the integration of the density to the number of electrons, and in case of a large difference prints this error message. Try a denser grid with the \

我是对Rh，P，Cl用lanl2dz+P，C，H, O 用6－31G**,这是计算频率时出错的 (Enter /home/mqx/g98/linda-exe/l1002.exel) Minotr: Closed-shell wavefunction.

3

Direct CPHF calculation.

Solving linear equations simultaneously.

Requested convergence is 1.0D-08 RMS, and 1.0D-07 maximum. Secondary convergence is 1.0D-12 RMS, and 1.0D-12 maximum. Differentiating once with respect to electric field. with respect to dipole field.

Differentiating once with respect to nuclear coordinates. NewPWx=T KeepS1=F KeepF1=F KeepIn=F MapXYZ=F. MDV= 26214400 Using IRadAn= 2.

Integrals replicated using symmetry in FoFDir. MinBra= 0 MaxBra= 3 Meth= 1.

IRaf= 0 NMat= 57 IRICut= 57 DoRegI=T DoRafI=T ISym2E= 2 JSym2E=2. Raff turned off since only 62.75% of shell-pairs survive. Defaulting to unpruned grid for atomic number 45. Defaulting to unpruned grid for atomic number 45.

There are 114 degrees of freedom in the 1st order CPHF. 114 vectors were produced by pass 0.

AX will form 57 AO Fock derivatives at one time. 114 vectors were produced by pass 1. 114 vectors were produced by pass 2. 111 vectors were produced by pass 3. 111 vectors were produced by pass 4. 94 vectors were produced by pass 5.

Tuple checksum error for IPart= 1 tuple: Prism -- DenA L= 12973128 Ck= 89DC7555 CkI= 89DC75DF. Error termination in TplChk.

Error termination via Lnk1e. 我现在明白我上面贴的Gaussian的问题是由linda引起的，我在不使用linda时不会错，还有我将linda重新装一下也可以算。

END TIME: Fri Mar 26 20:49:43 EST 2004

大家帮忙看看这个错误信息应该怎么修改！！输入文件为： %chk=ta %mem=200mb

# rb3lyp/gen opt=(ts,calcfc,noeigentest) freq scf=(maxcyc=500) pseudo=read dd 0 1

W -0.1010.233-0.265 C -0.3910.1311.601 C 1.7330.283-0.672 C -1.5091.491-1.251 C -0.657-1.786-0.652 H-0.019-0.3852.347

4

H-0.6840.6602.148 H2.569-0.0900.215 H2.1430.730-1.590

H -1.2352.548-0.937 H-2.5781.265-1.309 H-1.4111.403-2.116

H -0.250-1.849-1.692 H -1.640-2.053-0.851 H -0.001-2.534-0.193 O-2.822-0.4480.330 O-2.527-0.2721.427 C H O 0 6-31+G(d) **** W 0 lanl2dz **** W 0 lanl2dz 错误信息为：

Bad length for file.

FileIO: IOper= 1 IFilNo(1)= -582 Len= 0 IPos= 0 Q= 66388072

dumping /fiocom/, unit = 1 NFiles = 28 SizExt = 524288 WInBlk = 512 defal = T LstWrd = 1581568 FType=2 FMxFil=10000

Number 0 501 502 503 511 512 551 552 Base 1057792 143360 149504 201216 259584 376320 374784 373760 End 1581568 144360 171554 201328 272305 1056320 374809 373776 End1 1581568 144384 172032 201728 272384 1056768 375296 374272 Wr Pntr 375808 143360 169504 201216 259584 1056320 374784 373760 Rd Pntr 375808 143407 169754 201216 259584 376320 374809 373760 Number 559 562 575 579 583 598 665 670 Base 1056768 272384 201728 374272 1057280 172032 199168 375808 End 1056769 373741 259348 374281 1057289 172034 200206 375810 End1 1057280 373760 259584 374784 1057792 172544 200704 376320 Wr Pntr 1056769 272384 201728 374272 1057289 172032 199168 375810 Rd Pntr 1056768 373741 201728 374272 1057280 172032 199168 375808 Number 674 698 989 991 992 993 994 995 Base 200704 375296 144384 147456 146944 142848 20480 142336

5

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库gaussian基本概念和用法在线全文阅读。