《gaussian基本概念和用法》由会员分享，可在线阅读，更多相关《gaussian基本概念和用法（50页珍藏版）》请在金锄头文库上搜索。
1、1 计计 算算 流流 程程上面的图应该是从 L301 说起的，即图中的“基组”，L301 的作用是产生基组信息。（在这之前还有 L1：处理计算执行路径，创建执行链接的列表，并初始化 scratch 文件；L101：读取标题和分子说明部分；L103：berny 优化到最小值；L202：重新定位坐标，计算对称性，检查变量）“基组”之后是“分子结构”即 L401，形成初始的分子轨道初猜。第一个框中有“分子轨道”和“能量”两项，我想这应该指的是 L502 即迭代求解 SCF 方程。下面红线圈的框中则是迭代求解的具体过程。每次迭代完毕都得到一个能量。即（SCF Done: E(RHF) = -xx.xx
2、xx A.U. after xx cycles）然后是 L103 判断力和位移是否满足收敛条件。如果满足（即四个 yes），则进行布居分析等，并完成计算。如果不满足，则继续调整分子结构（L401），再次进行迭代求解，反复这个过程，直到满足力和位移收敛条件为止。如果不是 HF 方法而是多体微扰，CI 等后-scf 方法的话，在每次迭代（L502）完成后还要多一步相关能的计算，以 MP2 为例，迭代完成后多出:L801(双电子积分变换的初始化),L906(半直接的 MP2)和 L1002(迭代求解 CPHF 方程)等计算。 对于楼主的图，我觉得力和位移判断为 NO 以后，箭头应该指向最头上的“基组
3、”，因为实际过程中，每次构型调整都会重新定位坐标和产生基组信息。常见问题分析1 检查是否有初始文件错误 在命令行中加入 %kJob L301 or %kJob L302 如果通过则一般初始文件 ok。常见初级错误：a. 自旋多重度错误; b. 变量赋值为整数; c. 变量没有赋值或多重赋值; d. 键角小于等于 0 度，大于等于 180度e. 分子描述后面没有空行; f. 二面角判断错误，造成两个原子距离过近; g. 分子描述一行内两次参考同一原子,或参考原子共线 2 SCF(自洽场)不收敛 则一般是 L502 错误,省却情况做 64 个 cycle 迭代(G03 缺省 128 cycles)
4、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.
5、 修改坐标，使之合理; 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) 注：N512 7 优化过渡态 opt=TS (给出过渡态);opt=qst2 (给出反应物和产物); opt=qst3 (给出反应物和产物和过渡态)a. 用
6、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,-1b. 改变计算方法 MP2=Direc
7、t 可以少占硬盘空间c. 限制最大硬盘 maxdisk=N GB，*MB,有些系统写 2GB 会出错，可以写 2000MB 9. FOPT 出错 原因是变量数与分子自由度数不相等。 可用 POPT 或直接用 OPT 10. 优化过渡态只能做一个 STEP 原因是负本征数目不对 添加 iop（1/11）=1 或者 noeigentest一些容易忽略的地方1。低频率振动模式时，力常数很小，必须使用 opttight 以确保适当的收敛和随后任务步骤中频率计算的可靠度。我曾经算过一个体系，第一个频率振动的力常数只有 1020，很不放心，加入 tight 后，优出来的结构却与前面大相径庭，而且还出现了一
8、个虚频。2。做 IRC 计算的时候，如果步长很小，必须要用 verytight。就曾有人步长设为 1，而不加 verytight，算出来的反应路径是 v 字形倒海鸥状的。3。使用 int 指定用于数值积分的积分网格来用于消虚频的时候，必须注意，在比较能量的时候，对所有的计算要使用相同的积分网格。低频振动模式就是力常数比较小，至于有多大算低频振动模式，我并不太清楚。 不过根据我的经验，50 的时候不用加tight， Convergence criterion not met.SCF Done: E(UB+HF-LYP) = -2003.48215555 A.U. after 201 cycles
9、Convg = 0.1927D+00 -V/T = 1.9356S*2 = 0.7500Annihilation of the first spin contaminant:S*2 before annihilation 0.7500, after 0.7500Convergence failure - run terminated.Error termination via Lnk1e in F:G98Wl502.exe.Job cpu time: 0 days 0 hours 25 minutes 31.0 seconds.File lengths (MBytes): RWF= 6 Int
10、= 0 D2E= 0 Chk= 1 Scr= 1请问是拿出错了？是不是初始构行不好？ 有时自旋多重度的设定也会影响收敛，你是不是换换多重度的数值试试输入文件为：%chk=qwq %mem=200mb # b3lyp/sto-3g opt scf(maxcyc=200) . 先试试 scf=qc 看看我认为你可以尝试把 SCF 的收敛标准降低，比如 5，然后在此结果的基础上，恢复 SCF 的默认收敛标准，继续计算。下面引用由下面引用由 zjx1103zjx1103 在在 2004/04/272004/04/27 09:48am09:48am 发表的内容：发表的内容： 输入文件为：%chk=qwq
11、 %mem=200mb # b3lyp/sto-3g opt scf(maxcyc=200) . 对 Cu 用 sto3g？ 可行吗？可信吗？下面引用由下面引用由 helpmehelpme 在在 2004/05/082004/05/08 09:01am09:01am 发表的内容：发表的内容： 先试试 scf=qc 看看 初始构型太差的缘故，建议调整构型。 我的输入%chk=al#hf/6-31+g(3df) opt=(qst3,bimolecular) Guess=Alter freq test输出Rotational constants (GHZ): 1.5082594 0.6204235
12、0.5203045我的输入应该怎样改阿？还望各位高手们不吝赐教GAUSSIAN 输出文件中出错信息大汇集输出文件中出错信息大汇集G98 的出错讯息都在用户手册中的Gaussian 98 的链接一节中就给出了，以前版面上也曾有人贴出来过，你仔细看看。计算千差万别，出错的原因也各有不同，一时是总结不了的。但大部分都是：1。输入有误，这没话说了。一般都是 L1，L101 出错2。基组有误。特别是刚开始用自定义基组的，或者不知道基组适合范围的，这是 L301 出错3。无法收敛、优化和频率计算出错，这就很多了，比如 L502，L716，L9999。这就只好具体问题具体分析了。Range of M.O.s
13、 used for correlation: 20 167NBasis= 167 NAE= 57 NBE= 57 NFC= 19 NFV= 0NROrb= 148 NOA= 38 NOB= 38 NVA= 110 NVB= 110Semi-Direct transformation.ModeAB= 4 MOrb= 38 LenV= 38809382LASXX= 49931981 LTotXX= 49931981 LenRXX= 101093547LTotAB= 51161566 MaxLAS= 65294640 LenRXY= 0NonZer= 151025528 LenScr= 265347
14、674 LnRSAI= 65294640LnScr1= 136751342 MaxDsk= -1 Total= 568487203SrtSym= T ITran= 4File extend in NtrExt1 failed; probably out of disk space.*呵呵，老大， LINK DIE 容易找但是，比如我上面的出错信息是什么意思啊？ 硬盘空间不够了，哥们。 不过你说的对，像 NtrExt1 之类的，我的确看不懂：（ 我在用 OPT 的一个参数 STEPSIZE 的时候，从来没有成功过，不知道是为什么？# UHF/6-31 OPT(Stepsize=30)-QPERR
15、 - A SYNTAX ERROR WAS DETECTED IN THE INPUT LINE.# UHF/6-31 OPT(STEPSIZE=30) Last state=“OPT2“TCursr= 9023 LCursr= 15Error termination via Lnk1e in D:G98Wl1.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这个贴子最后由 helpme 在 2004/03/29 04:58pm 第 1 次编辑我试了一下，OPT(STEPSIZE=30) 确实不可以用。 （网上有人说过）可以用 IOp(1/8)来实现。默认的 IOP(1/8=30)IOp(1/8)L103, L109, L112: MAXIMUM STEP SIZE A
