|
geom=connectivity对当前计算完全没有意义,建议删去。IOp(1/8=x)和scf=qc这两个关键词有些需要注意的地方,详见《常见的多余的和被滥用的Gaussian关键词》
21 IOp(1/8=x)、IOp(6/7=3)
IOp(6/7=3)相当于pop=full、IOp(1/8=x)相当于在opt里定义maxstep,明明有对应的关键词,却还老看到有人在输入文件里直接写这俩IOp。有对应关键词的场合不要写IOp,否则自己容易记错,别人看这输入文件的时候还得查IOp表。
顺带强调一下,优化极小点时用maxstep定义步长上限的时候一定要带着notrust,即opt(maxstep=x,notrust),否则步长上限会随着优化进行自发发生变化而导致设定没起到作用。然而洒家经常看到有人用maxstep却不写notrust。而对于优化过渡态,本身默认就是用了notrust的,因此opt=TS任务还写notrust那是完全多余的。
22 scf=xqc或scf=qc
在《解决SCF不收敛问题的方法》中介绍过。很多人看到不收敛上来就用这个,甚至将之当成平日计算默认使用的关键词,这样做是很糟糕的!用qc(二次收敛)做迭代比一般方式迭代每步多花好几倍时间,虽然收敛的可能性高了不少,但是远不如先尝试其它解决不收敛的办法。而且,如果迫不得已只能用qc方式才能收敛,也极有可能最后得到的波函数并不稳定,此时也暗示很可能当前体系结构、方法、自旋多重度、基组等是不合适的。
scrf里写solvent=water也是无意义的,scrf默认就是水做溶剂。
需要确认当前体系带三个正电荷且自旋多重度为2。
用了genecp关键词,所以结尾需要定义基组和赝势,给出这个部分的截图。
|
|