mopac2xyz：将MOPAC的几何优化过程转换为xyz轨迹文件的工具

mopac2xyz：将MOPAC的几何优化过程转换为xyz轨迹文件的工具
mopac2xyz: A tool to convert geometry optimization process of MOPAC program into xyz trajectory file

文/Sobereva @北京科音
First release: 2013-Nov-12    Last update: 2021-Sep-20

 

最近用MOPAC优化一些大体系，虽然有现成的程序诸如gabedit可以观看优化后的结构，但是没法像在GaussView里那样观看优化过程中的结构 变化。而看不到优化过程的话，特别是对于较大体系，就心里没谱。另外这些支持MOPAC的可视化程序的界面也着实不好。遂写了个小程序，把MOPAC的优化任务的输出文件转化为xyz格式的轨迹文件（介绍见谈谈记录化学体系结构的xyz文件http://sobereva.com/477），然后将xyz文件直接拖到VMD里就可以很方便地观看优化过程中结构怎么变化了。

mopac2xyz 1.0.3版下载地址：http://sobereva.com/soft/mopac2xyz_1.0.3.zip

使用说明：

mopac2xyz用来把MOPAC的EF方式的几何优化的结构变化过程提取出来构成多帧xyz文件。然后直接将生成的xyz文件拖进VMD就可以观看优化过程的结构变化了。

当前mopac2xyz版本已测试能够兼容MOPAC 2012/2016。不保证能兼容其它版本的MOPAC的输出文件。

压缩包内的test.mop是示例输入文件，test.out是相应的输出文件。mopac2xyz.exe是编译好的Windows版程序，无后缀的mopac2xyz文件是编译好的Linux版程序。

启动mopac2xyz后先输入MOPAC的优化任务的输出文件名。由于优化过程可能步数很多，为了减少输出文件的帧数，程序会让你指定每多少步输出一次到.xyz文件里，如果输入1，那么每一步的结构都输出。程序处理完之后，结构变化过程就都存到当前目录下的.xyz文件中了。

注意MOPAC的输入文件必须用笛卡尔坐标，并且必须写PRNT=2关键词。并且输入文件第三行的内容必须是All coordinates are Cartesian（是为了让mopac2xyz能够定位）。如果任务不是EF方法优化，而是比如L-BFGS方法优化（对较大体系默认如此），那么此程序提取出来的只有初始结构和优化后的结构，也有可能提取不成功。此时用户必须加上EF关键词才能提取优化过程的所有帧。建议用户总带着PRNT=2 EF关键词。
 
有问题请电子邮件联系我（程序启动时显示了）。