20240703
结构优化
在优化中出现l9999错误,可能结构在优化过程中出现了震荡,无法达到收敛限,可以尝试一下方法:
从实用性出发,对于大多数情况,对于初学者笔者推荐按照以下顺序去尝试解决震荡问题:
1 对于DFT,特别是明尼苏达系列泛函(如M062X、M06、M06L),首先加上int=ultrafine再试(从G16开始此为默认,不再需要尝试此项)。如果用了SMD,改用IEFPCM
2 如果体系小,尝试opt=calcall。如果体系大算得慢,此关键词会令计算量进一步激增,实在没辙再考虑这个,或者尝试opt=recalc=3或5
3 尝试opt=gdiis
4 尝试opt(gdiis,maxstep=x,notrust),x为3~5
5 如果上述做法都不灵,干脆改成loose收敛限,但如果之后做振动分析容易有虚频
如改为:
# opt=gdiis b3lyp/3-21+g* scrf=(iefpcm,solvent=water) em=gd3bj
采用em=gd3bj加色散校正。
可以将失败的log文件导入gaussianview,勾选read intermediate geometries(optimizations)