用ORCA做结构优化的关键词非常简单,与Gaussian类似,直接写上opt即可,同样也可以使用opt freq的组合,在优化完结构之后进行频率计算。示例如下:
!BLYP D3 def2-SVP def2/J opt freq tightSCF
%maxcore 1000
%pal nprocs 24
end
*xyz 0 1
C -1.2471894 -1.1718212 -0.6961388
C -1.2471894 -1.1718212 0.6961388
N -0.2589510 -1.7235771 1.4144796
C 0.7315327 -2.2652221 0.6967288
C 0.7315327 -2.2652221 -0.6967288
N -0.2589510 -1.7235771 -1.4144796
H -2.0634363 -0.7223199 -1.2472797
H -2.0634363 -0.7223199 1.2472797
H 1.5488004 -2.7128282 1.2475604
H 1.5488004 -2.7128282 -1.2475604
C -0.3380031 2.0800608 1.1300452
C 0.8540254 1.3593471 1.1306308
N 1.4701787 0.9907598 0.0000000
C 0.8540254 1.3593471 -1.1306308
C -0.3380031 2.0800608 -1.1300452
N -0.9523059 2.4528836 0.0000000
H -0.8103758 2.3643033 2.0618643
H 1.3208583 1.0670610 2.0623986
H 1.3208583 1.0670610 -2.0623986
H -0.8103758 2.3643033 -2.0618643
*
ORCA的默认结构优化收敛标准如下:
在手册上说必须所有条件都满足时,结构优化才算收敛。但在实际计算中,如果能量没有收敛,程序也认为收敛了,如上图的情形。程序中也给出了说明。
ORCA没有官方的可视化程序可以观看优化轨迹。卢天曾开发过OfakeG程序,可以将ORCA的结构优化和频率计算文件转化成一个伪装的Gaussian输出文件,这样便可借助GaussView观看优化轨迹和振动模式。其操作非常简单,只需将ORCA的输出文件拖入OfakeG打开后的界面,即可生成一个fake的Gaussian输出文件,具体介绍可参看卢老师的帖子:
http://sobereva.com/498
ORCA的结构优化功能并没有高斯强大,但是ORCA在RI近似方面却比高斯做得要好,对于大体系,在ORCA中使用纯泛函加RI近似,计算速度是非常快的。可以将ORCA的能量和力传给Gaussian进行结构优化,具体可参看《使用external关键词将其他程序与高斯对接进行结构优化、IRC分析等计算》一文。