《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