用ORCA做结构优化及轨迹查看

用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分析等计算》一文。