《Gaussian中有用的IOp一览》由会员分享，可在线阅读，更多相关《Gaussian中有用的IOp一览（7页珍藏版）》请在金锄头文库上搜索。
1、*IOp(1/7=N) 设定 优化时的 收敛 限。力 的 RMS收敛限 设定 为 N*10-6( 即 Threshold 显示的值 ) ，最大受力被设为 1.5*N*10-6 。RMS位移会设为 4*N*10-6bohr ， 最 大 位 移 会 被 设 为 6*N*10-6 bohr 。 *IOp(1/8=N) 在优化的时候以 N*0.01 bohr/ 弧度为最初置信半径，如果优化中不动态更新置信半径，则每一步步长都不会大于这个值（并非是必须等于这个值）。减小其数值有助于解决收敛震荡问题，如果势能面平缓，则应当加大以加快 优 化 速 度 。 默 认 是30 。 等 价 于opt关 键 词 中m
2、axstep选 项 。IOp(1/9=x)设定对置信半径的处理。默认对于寻找极小点会每步自动更新置信半径，对于过渡态寻找则不会。=1代表不更新，=2代表更新，对应opt中NoTrustUpdate和TrustUpdate关键词。当步长超过置信半径时处理方法用10和 20 来选择， 10 代表将位移向量乘以刻度因子使其模等于置信半径（对寻找过渡态是默认的）， 20 代表在置信半径对应的超球面上寻找能量极小点（对寻找极小点是默认的）。*IOp(1/11=x)=1 即noeigentest关键词， =2 为总是测试势能面曲率符号，错了就停止任务，即 eigentest关键词。默认对 LST/QST方
3、法找过渡态及寻找极小点不测试，对 Z矩阵或笛卡尔下 TS任务进行检测。*IOp(1/40=N)每 N 步重新精确计算一次 Hessian 矩阵，期间只是使用一阶导数更新之前的Hessian矩阵。IOp(1/111=N)设定温度为N/1000 K ，若 N 为负数则为 N/1000000 K 。IOp(1/111=N)设定压力为 N/1000 atm ，若 N 为负数则为 N/1000000 atm 。IOp(2/12=x)默认是 1，仅当距离为 0时才报错，即 geom=nocrowd。=2 是小于 0.5埃就报错，即 geom=crowd。=3 是即便原子间距离为0，L202 也不报错，此选
4、项似乎目前不能用。IOp(2/14=x)使用内坐标时，控制是否测试内坐标线性依赖性， 默认不做测试。见手册geom=independent。*IOp(2/15=x)控制 L202 对对称性处理的细节。 nosymm相当于 1。2 是仍然将分子弄到标准朝向，但不利用对称性，3 是根本不调用处理对称性的模块。4 是用较松的标准判断对称性，然后将体系严格地对称化成这种对称性。IOp(2/16=x)如果优化过程中遇到点群改变怎么处理。1 是继续走但关了对称性， 2 是继续走仍使用旧对称性，3 是继续走并使用新对称性，4 是停止任务。默认是3。IOp(2/17=N)决定对称性时，距离比较的容差为10-N
5、。IOp(2/18=N)决定对称性时，非距离比较的容差为10-N。当IOp(2/17=4,2/18=3)时相当于symm=loose。*IOp(3/10=x) 设定以何种基函数方式储存波函数信息。 1000 和 2000 控制是否将指数相同的 S 和 P 壳层组合成 SP壳层，默认是 2000（组合）， 1000 代表分裂开（SplitSP关键词）。 100000 将收缩基组完全去收缩化，即基函数数目将与高斯函数数目一致。IOp(3/18=1)输出赝势信息，和用了GFPrint关键词的效果一样。*IOp(3/24=1) 1、10、100 可以分别以老方式、Gen的输入格式、某种易读的方式输出基
6、组信息。 1000 输出壳层所属中心坐标。 x0000、x00000 可以控制输出的基组、密度拟合基组信息是否先做归一化。IOp(3/26=11)控制 L311、L314（ scf=conventional时会在自洽场迭代前调用）的双电子积分精度，默认是计算到10-10 精度。设为 11 可以做到尽可能精确。IOp(3/27=N)扔掉数值小于10-N以下的双电子积分信息，默认N=10。IOp(3/29=N)使单电子积分（L302）的精度做到10-N，默认是N=13。*IOp(3/32=2)避免去除线性相关的基函数。默认是检验重叠矩阵本征值，发现存在线性相关问题就去掉多余基函数。*IOp(3/3
7、3=1)输出单电子积分，核哈密顿，势能，重叠矩阵(overlap那项 ) 。=3 用标准格式输出双电子积分，=4 则用 debug 格式输出（对半经验无效）。 =5同时输出单、双电子积分。IOp(3/53=N)设 定ECP积分精度为 10-N，-1则不截断。*IOp(3/59=N)与 3/32相关，控制扔掉重叠矩阵本征向量的标准为10-N。IOp(3/60=x)控制对广义收缩基组的正交化和简化处理。 -1是不做正交化和简化，默认是2，即做正交化并去掉系数小的GTF。1 是正交化并只去除系数为 0的 GTF，N是做正交化并去除系数小于 10-N 的 GTF。在 Gaussian 与其它量化程序计
8、算结果相比较时此选项有用。IOp(3/74=x) 设用什么密度泛函， 负数都是交换和相关泛函组合好的， 包括那些杂化泛函，如 -5=B3LYP，-24=O3LYP，-33=X3LYP。 0-99 以内的都是代表相关泛函， 100 及以上的都是代表交换泛函，这两个可以任意组合起来，比如402=BLYP。*IOp(3/76=MMMMMNNNNN)将 HF 交换项和DFT 交换项按照NNNNN/10000比MMMMM/10000的比例组合。*IOp(3/77=MMMMMNNNNN)使LDA和高阶交换泛函按照NNNNN/10000比MMMMM/10000的比例组合。*IOp(3/78=MMMMMNNN
9、NN)使LDA和高阶相关泛函按照NNNNN/10000比MMMMM/10000的比例组合。*IOp(3/86=N)扔 掉 角动 量 量子数=N的基 函数 。*IOp(3/93=1) 设定核电荷用点电荷描述， 而 Gaussian 对标量相对论计算默认用的 是 s 型 高 斯 函 数 描 述 ， 为 了 与 文 献 比 较 有 时 需 要 设 定 此 项 。IOp(3/111=1)PBC计算一开始额外输出每个晶胞与中心的距离。IOp(3/116=x)具体控制 SCF任务的类型，比如限制性、非限制性、实数、复数等问题。IOp(3/118=N)给估计 的直 接积分任务 所需 内存量增加N words
10、 内 存。IOp(3/124=1)不管是什么交换相关泛函， 都加上 VDW校正项。而默认情况下，只有所用泛函在定义中包括了色散校正项才会加上。*IOp(3/125=x)给 MP2 计算的二阶微扰能 (E2) 中自旋平行项和自旋相反项贡献设定校正系数， -3 是让自旋相反项为 0，-2是让自旋平行项为 0， -1是让两项都为 0（即还原到 HF结果）。设为 MMMMMNNNNN是让自旋平行项与自旋相反项以MMMMM/10000比 NNNNN/10000的比例组合，利用这个可以实现比如 SCS-MP2，即设为1200003333。IOp(4/6=x)设定读入的初猜波函数是否做正交归一性检查、正交化
11、和投影。IOp(4/11=x)设 置 初 猜 类 型、 Harris初猜时所用的积分格点设定。*IOp(4/15=N)设 定初 猜时 的自 旋多 重度 为N 。IOp(4/39=N)设 定guess=mix时的角度为Pi/N， 默 认N=4 。*Iop(4/43=3)在 CASSCF任务中且组态数很多时， 使所有组态的信息都输出出来（默认是仅在组态数少的时候才输出）。IOp(4/68=N)设定分子力学优化微迭代的收敛标准为10-N。*IOp(5/7=-1) MCSCF 中不做迭代只做CI。设为 N时 SCF和 MCSCF最大迭代次数为N次。IOp(5/9=1)使用 SCF直接最小化 (L503
12、) 使密度矩阵收敛到 10-3时切换到普通SCF过程。*IOp(5/13=1)SCF不收敛时仍继续执行后续任务。默认是0（等价 2），不收敛就停。*IOp(5/14=x)For L502，=1 对于 RHF任务，迭代结束后以UHF波函数形式记录，对 UHF 任务起到跳过自旋湮灭的作用（注意得到的密度矩阵将完全错误）； =2 保存自旋湮灭后的密度矩阵，故波函数分析模块 L607 使用的波函数将会是自旋湮灭后的；10=在第一次迭代计算完双电子项后结束； =20 仅从储存的实空间 Fock矩阵重新计算能带结构。IOp(5/16=x)ForL502，设定SCF中对角化方法。IOp(5/17=x) For L510，设定MCSCF中对角化方法，1000=为多参考态MP2生成数据，10000=试图控制根翻转。*IOp(5/33=3)在每次迭代后输出核哈密顿矩阵、MO系数、本征值、密度 / 自旋密度矩阵、 Fock 矩阵等一大堆矩阵。 PBC 计算时在最后迭代后还能输出不同k点的晶体轨道能。IOp(5/37=N)每 N 次迭代完整构成一次 Fock 矩阵，而不是默认的每 20 次一次（期间通过增量方法构成Fock 矩阵）。 -1则不用增量方法构成Fock 矩阵
