高斯常见问题

《高斯常见问题》由会员分享，可在线阅读，更多相关《高斯常见问题（21页珍藏版）》请在金锄头文库上搜索。

1、、混合基组#p b31yp/GenECP opt Gen自己在输入文件末尾定义基组，ECP使用贋势 空行g03 opt空行0 1空行C H O 0碳 氢 氧 数字0GEN6-31G*基组名*结束符号Ru P Cl 0Ru P Cl 数字 0GENLanl2dz基组名*结束符号空行Ru P Cl 0Ru P Cl 数字 0ECPLanl2dz基组建议你查看GEN关键字。上述，如不加ECP亦可。二、出现错误提（一）出现名言提示，表明是正常结束的； l9999 终止，是因为循环次数不够。建议用 optcyc=100-999 增加循环；或者把输出结果做为输入结构再次优化！1、是我感觉你要是maxcyc

2、le=200还是999死的话，就减小一下步长试试，可能会有用的！ maxstep=12、该先看下收敛指标是不是很接近了，如果很接近就继续做（读chk）或取出最后一个坐 标继续算，如果相差很远，可以加非线形实验下，再一个就是加大循环，缩短步长。（二） L202是重新定位坐标，计算对称性，检查变量的，它出错应该就是输入结构有问题。三、过渡态寻找过渡态的时候,可以用 TS,QST2 和 QST3 方法,那一种方法的效果更好一些?” 取决于 计算前所拥有的信息量。当然，信息越多，结果越有把握。所以， QST3 最有把握。如果没 有任何初始信息，那就没有什么选择了。“我想不同的输入方法肯定不会产生一样的

3、结果， 那输入过程中限制它们间相对位置的条件是什么呢？” 如果输入不同会导致结果不同，那 么应当全面考虑。首先可以用化学知识筛选。如果此后仍有多种可能性，就应当通过计算比 较。我算过的用HF方法虚频超过1000，而用DFT只要100-300.所以说方法很重要。但是，只 要你找到的TS虚频与反应的方向一致并且IRC已经验证，这个TS就是可信的。虚频也不 能太小，比如小于 100.过渡态与你的产物相差甚远 说明在这条反应路径上还有过渡态你没有找到在优化过渡态的时候圈数设为100，步长 5，怎么输入opt (TS, calcFC, maxcyc=100, maxstep=5)是这样寻找过渡态不是一件

4、容易的事(对于我和大多数刚涉及量化的人来说)，因此我希望通过写 这个经验小结能对大家有些帮助。1. 首先遇到的问题是,用哪种方法来寻找过渡态？GAUSSIAN提供的方法是QSTN和TSN方法。两种方法各有优点和缺点。QSTN方法特 别 QST3 方法要求输入反应物,过渡态的猜测结构,产物这三者的结构。特别麻烦。但很管 用，一般不会出现不收敛的情况。对于TSN (对应关键词为OPT=TS)方法，只要求输入过 渡态的初始结构,但这个初始结构非常的关键,如果结构不好,则很容易出现不收敛的情况。 所以我建议，如果是刚开始做过度态的话，用QSTN方法是好的选择，等有了感觉之后， 再用 TSN 方法。2.

5、 怎么解决经常出现的错误？ 在找过度态的时候，经常碰到的一些问题就是(1)不收敛，(2)有一个错误的本征值(错误信息为：there is a wrong sign eigenvalue in hessian matrix.)(3)和 LINK9999 错误 导致退出。对于不收敛的情况，可以分为两类，比如提示信息里的CONVERGENCE FAILER提 醒收敛到了 10 (-5),而此时你设定的SCF循环次数也仅仅是64步，那么完全有希望通过 加大SCF循环次数来达到收敛的目的。倘若只收敛到10(-3)或10 (-2),此时加大循环 次数可能就没用了。结果还是CONVERGE FAILED此时

6、可采用SCF=QC,来达到强制收敛 的目的。因为SCF=QC (LINK508)的计算量比默认的L502要大，所以不到万不得以就不 用它了。出现第二个错误可以直接用关键词 OPT=NOEIGEN 来实现。 LINK9999 出错是因 为已经走完了默认的步数，但还未完成。系统会自动跳出。出现这种情况大多数就是因为优 化步数和SCF步数超过了默认值。可用OPT (MAXCYCLE=100)和SCF (MAXCYCLE=300) 来改错。3. 怎么样控制过渡态的优化，使得过渡态不至于收敛到其他的分子结构中去？我用GAUSS VIEW可以解决这个问题，当刚开始运行GAUSSIAN时，你用GVIEW去打

7、 开输出文件时，你可以看到你的过渡态的初始输入结构，当一个循环过后(从上一个LINK502 到下一个LINK502)，你再打开输出文件，你就可以清晰地看到优化一步后分子的构型，这 样就可以随时监控过度态分子的结构，倘若已经有收敛到其他分子构型的趋势时，你就可以 把它给KILL 了，而不至于需要等全部工作结束后，打开输出文件才知道已经不是想要的过 渡态了。如果收敛到其它的构型上去，可以考虑缩小OPT的步长.iop(1/8=2或3)即可。4. 还需要加其它的关键词吗？建议在OPT中加入CALCFCo这样可以加大找到过渡态的几率。本人深有体会！5. 如何控制寻找过渡态的步长？这个可以用IOP来实现。

8、具体相应的IOP是iop(1/8)=2,此语句说明在找过度态的时候， 以2A (A为基本单位长度)为单位来寻找过渡态。6. 我用高斯03做了两个稳定点A和B的过渡态TS,发现频率文件的振动方向还不足以证 明TS就是A和B的过渡态，所以，我又接着做了IRC,结果发现IRC正方向走出的结构与 A结构差别很小，而反方向走出的结构与B相差太大，是不是说明这个求出的TS不是AB 的过渡态？那就应该不是，可是你的这个ts也是反映路径当中的一个ts,他前后走出的东西也是im,你如 果想要得出AB之间的过渡态的话，最好就用qst2和qst3方法找.四、奥奥的输入文件%mem=8mw%rwf=1,2000mb,

9、2,2000mb,3,2000mb%chk=ts-1+#P B3LYP/6-311G(d) fopt=(ts,Z-matrix,calcfc) optcyc=90 iop(5/13=1) iop(9/6=150) iop(1/11=1) freq*1*1 2CS,1,r1X,1,1.,2,a1Si,1,r2,3,a1,2,180.,0Variables:r1=2.85080391r2=1.88789137a1=21.74742362iop(5/13=1) iop(9/6=150) iop(1/11=1)吉大 ADF 20:54:02 第一个是优化用的 第三个是找过渡态用的吉大 ADF 20:5

10、4:15第二也可以不加吉大 ADF 20:54:44比较没有第一个 你就是基本的优化都会出问题 ComputerGirl 20:54:55这几个字段是连在一起用的吗？还是做优化用一个，做过渡态用另一个吉大 ADF 20:55:18优化稳定体时 第一个就行吉大 ADF 20:55:29过渡态时 第一个和第三个吉大 ADF 20:55:33第二个 不用加 第二个 算能量时会用到Image point= NEB五、IRC计算程老师：IRC=(forward 或 backward stepsize=n) n 是个自然数，一般可以为 6. 一定要有你算过渡态的 CHK 文件啊！是从 CHK 文件读数据的

11、。不是有混合基组的话，如 果有多个虚频，说明未找到能力最高点。有机金属配合物很复杂，又使用混合基组，所以情 况很复杂，要靠 IRC 内禀反应坐标扫描了。在我看来，这样的计算就未必准。结果非常依 赖基组的大小。DFT那是半经验方法哦在断裂键和生成键上家极化函数好一些，也好 不到哪去。不必了。这种方法用频率验证意义不大程老师： 可以进行频率验证的，过渡态一般有且只有一个虚频(负频)，大小不一定，与要断裂键或 新生成键的键长有关系。程老师 14:14:01一般要IRC扫描确定过渡态。如果熟练的话，可以不用IRC扫描即可确定过渡态宏师姐 20:27:35 算完后要导出坐标，在过渡态的基础上算频率的 宏

12、师姐 20:28:36过渡态直接ts就可以了，可以和频率同时算。不需要再优化宏师姐 20:29:13 出现虚频是很正常的，要把他们消去。只有一个的虚频才算是过渡态 宏师姐 20:31:26 不用，可以再优化，有时候如果结构有对称性可能会出现虚频。或者你上网上查怎么消虚频， 然后把它消去，消去以后在优化。过渡态比较难找的宏师姐 20:32:25 我消虚频的办法是找到虚频对应的坐标把它加在新的优化好的坐标上，然后重新优化。宏师姐 20:33:31只要有一个虚频都可以做ire,也就是你只要找到过渡态就可以做ire宏师姐 20:43:44很小的虚频也不行的，过渡态只能有一个虚频。过多的虚频是要消去的，

13、那就表明你现在的 结构不是过渡态，那么你在做IRC是没有意义的小木虫上说对于ire不收敛你可以加iop (1/13=5)试试关于过渡态与反应物接近，还与是与产物接近，有一个它Hammonds Postulate,大意是：如果在反应过程中有两个状态(比如说一个过渡态和一个反应中间体)相继出现，并且这两 个状态的能量相差不大，那么这两个状态相互转变时的分子构型也不会有太大变化。这个假定的一个有用推论就是，过渡态的分子结构会与在能量上和它相近的组分相似。如果 一个反应是放热反应，那么反应物的能量高于产物的能量，但过渡态能量最高，在能量上过 渡态就与反应物相近，根据Hammonds Postulate

14、，它的分子结构也与反应物相近。这个过渡 态就是前过渡态。要是走构型就直接#b3lyp/6-311g* IRC(MAXPOINT=400,CALCFC,stepsize=5,MAXCYCLE=600,forward)SCF(MAXCYCLE=400)1请问，在IRC计算中出现L 5 0 8错误应该怎么办啊？谢谢sef eannot eonverge even using sef=qe,there is no good enough solution for this error but to use better initial guess for calculation.e.g. use gu

15、ess=read to read sef from similar molecule ( cation, anion, different spin states, neibor in perodie table, ete.); or ehange starting geometrylink 下面加上这个 geom=check guess=check频率在进行Gaussian计算时，经常会出现多个虚频个，每个都只有一30一20左右，但是如果结 构改动过大，则一不小心，虚频就会更大。消除虚频的主要方法是改变构型。其次在计算上 还可以尝试：nosymm；加大循环次数；提高收敛度；iop(l/8=l)等。通常对于势能面很平， 较小的虚频很难消除。为了解决以上的这种现象,按如下的步骤一步步调整：1) .在优化时采用Scf(tight)的选项，增加收敛的标准。再去计算频率。如果还有虚频， 参见下一步。2) . 对称性的影响，很多情况下的虚频是由于分子本身的对称性造成的。这样，在优化 时，如果必要，要将对称性降低，还有，输入文件有时是用内坐标。建议如果有虚频的话， 将内坐标改成直角坐标优化。但是有时这一方法也并没有多大作用。3) . 如果上述方法还有虚频，看一下虚频，找