Gaussian中有用的IOp一览.doc

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 关 键 词 中 maxstep 选 项 IOp(1/9=x) 设定对置信半径的处理默认对于寻找极小点会每步自动更新置信半径，对于过渡态寻找则不会=1代表不更新，=2代表更新，对应opt中NoTrustUpdate和TrustUpdate关键词当步长超过置信半径时处理方法用10和 20 来选择， 10 代表将位移向量乘以刻度因子使其模等于置信半径（对寻找过渡态是默认的）， 20 代表在置信半径对应的超球面上寻找能量极小点（对寻找极小点是默认的）。

IOp(1/11=x)=1 即noeigentest关键词， =2 为总是测试势能面曲率符号，错了就停止任务，即 eigentest 关键词默认对 LST/QST方法找过渡态及寻找极小点不测试，对 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=nocrowd2 是小于 0.5埃就报错，即 geom=crowd3 是即便原子间距离为0，L202 也不报错，此选项似乎目前不能用IOp(2/14=x)使用内坐标时，控制是否测试内坐标线性依赖性， 默认不做测试见手册geom=independentIOp(2/15=x)控制 L202 对对称性处理的细节 nosymm相当于 12 是仍然将分子弄到标准朝向，但不利用对称性，3 是根本不调用处理对称性的模块。

4 是用较松的标准判断对称性，然后将体系严格地对称化成这种对称性IOp(2/16=x) 如果优化过程中遇到点群改变怎么处理 1 是继续走但关了对称性， 2 是继续走仍使用旧对称性， 3 是继续走并使用新对称性， 4 是停止任务默 认 是 3 IOp(2/17=N) 决定对称性时，距离比较的容差为 10^-NIOp(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的输入格式、某种易读的方式输出基组信息 1000 输出壳层所属中心坐标 x0000、x00000 可以控制输出的基组、密度拟合基组信息是否先做归一化。

IOp(3/26=11) 控制 L311、L314（ scf=conventional 时会在自洽场迭代前调用）的双电子积分精度，默认是计算到 10^-10 精度设为 11 可以做到尽可能精确IOp(3/27=N)扔掉数值小于10^-N以下的双电子积分信息，默认N=10IOp(3/29=N)使单电子积分（L302）的精度做到10^-N，默认是N=13IOp(3/32=2) 避免去除线性相关的基函数默认是检验重叠矩阵本征值，发现存 在 线 性 相 关 问 题 就 去 掉 多 余 基 函 数 IOp(3/33=1) 输出单电子积分，核哈密顿，势能，重叠矩阵 (overlap 那项 ) 3 用标准格式输出双电子积分，=4 则用 debug 格式输出（对半经验无效） =5同时输出单、双电子积分IOp(3/53=N)设 定ECP积分精度为 10^-N，-1则不截断IOp(3/59=N)与 3/32相关，控制扔掉重叠矩阵本征向量的标准为10^-NIOp(3/60=x)控制对广义收缩基组的正交化和简化处理 -1是不做正交化和简化，默认是2，即做正交化并去掉系数小的GTF1 是正交化并只去除系数为 0的 GTF，N是做正交化并去除系数小于 10^-N 的 GTF。

在 Gaussian 与其它量化程序计算结果相比较时此选项有用IOp(3/74=x) 设用什么密度泛函， 负数都是交换和相关泛函组合好的， 包括那些杂化泛函，如 -5=B3LYP，-24=O3LYP，-33=X3LYP 0-99 以内的都是代表相关泛函， 100 及以上的都是代表交换泛函，这两个可以任意组合起来，比如402=BLYPIOp(3/76=MMMMMNNNNN)将 HF 交换项和DFT 交换项按照NNNNN/10000比MMMMM/10000的比例组合IOp(3/77=MMMMMNNNNN)使LDA和高阶交换泛函按照NNNNN/10000比MMMMM/10000的比例组合IOp(3/78=MMMMMNNNNN)使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 内 存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，即设为1200003333IOp(4/6=x)设定读入的初猜波函数是否做正交归一性检查、正交化和投影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^-NIOp(5/7=-1) MCSCF 中不做迭代只做CI设为 N时 SCF和 MCSCF最大迭代次数为N次IOp(5/9=1)使用 SCF直接最小化 (L503) 使密度矩阵收敛到 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 矩阵。