Gaussian中有用的IOp一览

1、*IOp(1/7=N) 设定优化时的收敛限。力的RMS 收敛限设定为NIOp(2/15=x) 控制 L202 对对称性处理的细节。nosymm相当于1。 2是仍然将分子弄到标准朝向，但不利用对称性，3 是根本不调用处理对称性的模块。4 是用较松的标准判断对称性，然后将体系严格地对称化成这种对称性。10-6( 即Threshold 显示的值)，最大受力被设为*N*10-6 。 RMS 位移会设为4*N*10-6bohr，最大位移会被设为6*N*10-6 bohr。*IOp(1/8=N) 在优化的时候以N* bohr/弧度为最初置信半径，如果优化中不动态更新置信半径，则每一步步长都不会大于这个值(

2、并非是必须等于这个值)。减小其数值有助于解决收敛震荡问题，如果势能面平缓，则应当加大以加快优化速度。默认是30。等价于opt关键词中maxstep选项。IOp(1/9=x) 设定对置信半径的处理。默认对于寻找极小点会每步自动更新置信半径，对于过渡态寻找则不会。=1 代表不更新，=2 代表更新，对应opt 中NoTrustUpdate和 TrustUpdate关键词。当步长超过置信半径时处理方法用10和20 来选择，10 代表将位移向量乘以刻度因子使其模等于置信半径(对寻找过渡态是默认的)，20 代表在置信半径对应的超球面上寻找能量极小点(对寻找极小点是默认的)。*IOp(1/11=x) =1

3、即 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=no

4、crowd。 =2 是小于埃就报错，即geom=crowd。 =3 是即便原子间距离为0， L202 也不报错，此选项似乎目前不能用。IOp(2/14=x) 使用内坐标时，控制是否测试内坐标线性依赖性，默认不做测试。见手册geom=independen。tIOp(2/16=x) 如果优化过程中遇到点群改变怎么处理。1 是继续走但关了对称性，2 是继续走仍使用旧对称性，3 是继续走并使用新对称性，4 是停止任务。默认是 3。IOp(2/17=N) 决定对称性时，距离比较的容差为10-N。IOp(2/18=N) 决定对称性时，非距离比较的容差为10-N。 当 IOp(2/17=4,2/18=3)时

5、相当于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 可以控制输出的基组、密度拟合基组信息是

6、否先做归一化。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/33=1) 输出单电子积分，核哈密顿，势能，重叠矩阵(overlap 那项)。 =3用标准格式

7、输出双电子积分，=4则用 debug格式输出(对半经验无效)。=5同时输出单、双电子积分。IOp(3/53=N) 设定 ECP 积分精度为10-N， -1 则不截断。10-N*IOp(3/59=N) 与 3/32 相关，控制扔掉重叠矩阵本征向量的标准为IOp(3/60=x) 控制对广义收缩基组的正交化和简化处理。-1 是不做正交化和简化，默认是2，即做正交化并去掉系数小的GTF。 1 是正交化并只去除系数为0 的GTF， N 是做正交化并去除系数小于10-N 的 GTF。在Gaussian与其它量化程序计算结果相比较时此选项有用。IOp(3/74=x) 设用什么密度泛函，负数都是交换和相关泛函

8、组合好的，包括那些杂化泛函，如-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=MMMMMNNNNN)使 LDA 和高阶相关泛函按照NNNNN/10000 比MMMMM/1

9、0000 的比例组合。* 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 校正项。而默认情况下，只有所用泛函在定义中包括

10、了色散校正项才会加上。IOp(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-N。* IOp(5/7=-1) MCSCF 中不做迭代只做CI。 设为 N 时 SCF和 MCSC

11、F 最大迭代次数为 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矩阵重新计算能带结构。

12、IOp(5/16=x) For L502，设定 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 矩阵。* IOp(5/38=x) 控制

13、直接SCF 过程(每次迭代重算积分)中积分精度。默认=5，允许VarAcc(积分精度可变)， =1 不允许， =2 比当前收敛情况精度高三位小数，=3 在最终迭代前用当前收敛情况同样的精度，最终迭代时高两位小数。=4 首先收敛到10-5(这个过程电子积分精度为10-6)，然后再完全收敛到期望精度，=7 全精度双电子积分。* IOp(5/70=N) 设定解决SCF 收敛问题的分数占据方法(温度展宽)的温度为N。IOp(5/75=N)(5/78=N) 分别设定分数占据涉及的Alpha 电子数、Alpha 轨道数、Beta电子数、Beta轨道数。IOp(5/80)、 IOp(5/81)、 IOp(5

14、/82) 设定SCF中共轭梯度法的选项，分别设定最大步长、所用方法和最大步数、收敛标准。IOp(5/83=N) SCF 中 DIIS 最大子空间维度，默认是20。IOp(5/86=x) 设定能量升高时是否降低DIIS 子空间维度，何时开启能级移动及能级移动的量。IOp(5/102=N) 设 CAS-MP2 最大组态数是N，默认是1000。* IOp(5/103=N) 设置 IOp(5/33=3)输出的每个k 点的占据和非占据轨道数目都为N，默认是5。* IOp(6/7=3) 输出所有的分子轨道信息，即pop=full。* IOp(6/8=3) 输出密度矩阵。* IOp(6/17=x) 设定 L

15、602(诸如 prop关键词会涉及)计算电场及梯度、静电势属性时包含哪些成分的贡献。=0 计算全部贡献，=1 只计算原子核的贡献，=2 只计算电子的贡献，-N 只计算 N 号壳层的贡献。IOp(6/20=x) 设定拟合静电势电荷的方法和元素半径的定义。-1=读入自行设定的拟合点，1000=只拟合重原子电荷，10000只拟合活性原子电荷。* IOp(6/26=x) 设定L602、 L604模块中使用哪种密度，默认是1=总密度，2=Alpha密度(好像程序不认)、3=Beta密度、4=自旋密度。IOp(6/27=22) 不 仅 输 出 Mulliken 电 荷 ， 还 输 出 Mulliken 键

16、 级 矩 阵 ， 即pop=bonding，实际上只是把Gross 布居对角元设成了0 而已。此选项数值千位数若是 1，可以读入自定义权重而不是均分交叉项来做Mulliken 分析。* IOp(6/33=2) 输出拟合静电势时的点(单位为埃)和每个点的精确静电势.)。IOp(6/3539、 6/5355) 详细设定aim 模块(L609)的任务类型和参数* IOp(6/41=N) 拟合ESP电荷时每个原子用N 个同中心的点层为，默认N=4。* IOp(6/42=N) 拟合ESP电荷时每单位面积的点为N 个， 默认 N=1 。 17时每个原子大概用2500 个点，对于大分子可以稍微降低这个值以减

17、少计算量。IOp(6/75=x) 输出 CM2 电荷。 Lowdin 电荷、 Mayer 键级也一起输出。x 应设为计算波函数所用的理论方法和基组，可以为以下值：1 = HF/MIDI!(5D)2 = HF/MIDI!(6D)3 = HF/6-31G*4 = BPW91/MIDI!(5D)5 = BPW91/MIDI!(6D)6 = B3LYP/MIDI!(5D)7 = BPW91/6-31G*8 = HF/6-31+G*9 = HF/cc-pVDZ10 = BPW91/DZVP(6D)11 = AM112 = PM314 = ZINDO/S2IOp(6/79=1) 输出 Hirshfeld

18、电荷(第一列是Hirshfeld 电荷，后三列分别是原子的 x/y/z 偶极矩.)，G09 可以直接用pop=hirshfeld。 g03中必须用6d 基函数，否则结果是NaN， g09 没这个问题。=3 时还利用原子偶极、四极矩计算原子间相互作用能。* IOp(6/80=1) 输出 Lowdin 电荷、 Mayer 键级。* IOp(6/83=111) 输出每个原子的不同角动量原子轨道的布居数，以及每个分子轨道中各种角动量原子轨道的成分，分子轨道中d、 f 原子轨道的贡献和及具体组成也会被输出。* IOp(6/84=-1) 6/83默认只分析占据轨道，此指令让其分析所有轨道。设为N 是分析全

19、部占据轨道和N 个最低空轨道。* IOp(7/32=x) 将能量、核坐标、受力、力常数之类的信息输出到scratch/文件中，可以通过参数设定输出哪些内容，以及什么坐标下、什么格式。IOp(7/33=1) 做 freq 时，输出更多细节信息，如偶极矩导数、极化率导数、振动 极化率来自各振动模式的贡献等。* IOp(7/64=3) 在 freq任务中， 将振动频率做洛伦兹展宽成为红外光谱x-y数据点，可以导入诸如origin 等程序绘制红外光谱图。如果也做了Raman， 也会得到Raman光谱 x-y 数据点。IOp(7/90=N) 对大体系只分析最低N 个振动模式，此功能似乎目前不能用。IOp

20、(9/25=x) 控制每一次迭代、或是第一次和最后一次迭代，或是不输出电子对 儿对相关能的贡献。IOp(9/26=1) 归一化后HF 波函数(微扰波函数除外)，即令总波函数的模Norm(A)=1 ，可参考手册CID 关键词。默认只是对参考态HF 行列式波函数归一化 (Intermediate normalization)。IOp(9/30=x) 控制计算单粒子密度矩阵。IOp(9/40=x) 对于L906(半直接MP2) ， =1 时 MP2 使用 CASSCF 波函数作为参考态，相当于CASSCF 任务中同时写上MP2 关键词。对于L914( CIS、 TDDFT等任务)，输出每个激发态信息时将其中组合系数大于10-x 的组态输出出来，默认是 x=1 。IOp(9/44=3) 算基态密度以及来自基态的跃迁密度，还计算所有激发态间的跃迁密度。默认是2，只算算基态密度以及来自基态的跃迁密度。1 是计算每个激发态密度。这些结果会被填到rwf 文件633 位置内。IOp(9/46=N) Davidson 迭代中(在CIS、 TDDFT 等任务中都会用到)，最低N个态系数收敛就宣告迭代阶数，而不是等所有要算的激发态的系数都收敛才结 束