Gaussian中如何为不同原子指定不同基组[精].docVIP

Gaussian中如何为不同原子指定不同基组???? We may need a larger basis to describe the atoms involved in reaction, whil ??e describing ?? the rest of the system with smaller basis sets. ?? In this case, we need the Gen keyword ?? A sample input: ?? #gfinput iop(6/7=3) ?? #B3LYP/Gen Opt ?? …… ?? blank line ?? N1 N2 N3 0 ?? 6-311++G(d,p) ?? **** ?? N4 N5 N6 0 ?? 6-31G ?? **** ?? blank line 17:08?| 添加评论?| 固定链接?| 查看引用通告 (0)?| 写入日志 如何寻找transition state 如何寻找transition state??Answer: ?? A sample route section ?? #gfinput iop(6/7=3) ?? #B3LYP/6-31G(d) Opt(TS,Noeigen) ?? In order to increase the efficiency of the saddle point search,we could cal ??culate ?? the force constants by adding CalcFC keyword. ?? #gfinput iop(6/7=3) ?? #B3LYP/6-31G(d) Opt(TS,Noeigen,CalcFC) ?? We can also ask Gassian to automatically generate a gues structure for the ??reaction ?? by using keyword QST3 or QST2 ?? #gfinput iop(6/7=3) ?? #B3LYP/6-31G(d) Opt(QST3,Noeigen,CalcFC) ?? A+B--C Reactant //title section ?? 0 1 ?? structure of A+B ?? A+B--C Product ?? 0 1 ?? structure of C ?? A+B--C TS ?? 0 1 ?? guess structure for the TS ?? Note:the corresponding atoms need to appear in the same order within all th ??e molecule specifications. ??发信人: ghb (Never is a long time), 信区: Gaussian ??标 题: Re: 如何寻找transition state? ??发信站: BBS 大话西游站 (Sat Jan 5 11:10:02 2002) ??找TS好像也不是那么简单 ??我试了一下用QST2优化一个光环化反应的TS ??用PM3方法，竟然out文件有240M！ ??而且link died at L9999 ??ft死了 ??仔细想想，其实也就是对反应物分子和产物分子的Redandunt coordinate ??按照设定的path＝N做N等分，从这条路径找到一个近似的TS ??如果用Z-matrix一般都不会收敛 ??然后再逐步调整，根据Hessian判断是否到达了真正的TS ??我试过这样做，不知道是否有用，大家讨论讨论吧： ??用WinMopac2.0也可以做半经验的IRC计算，但不一定要从TS开始 ??可以选定Reaction Coordinate，让他从反应分子变到产物分子 ??例如对某个键断裂反应，选键长为反应坐标 ??可以写成这样的形式：1.5,1.6,1.7,1.725,1.75,…,3.4,3.6,4.0,… ??在可能的TS附近可以写多一点 ??然后开始IRC，这样能够得到每个反应坐标处的一个初始构型 ??当然也有能量，键级等等 ??Winmopac做这个非常方便，只要吧要做IRC的设为-1，不优化的设为0，优化的设为1 ??利用这些初始构型，固定每个反应坐标，用比较高的基组进行优化 ??比如B3LYP/6-31G* ??这样得到一系列构性，其基态能量最高者是否就是TS？ ??当然许多情况下还应该计算