用VESTA观看cp2k计算的振动模式的脚本

zzssxx

sob老师的MfakeG（http://sobereva.com/656）可以将cp2k的振动模式转化为GaussView可看的模式，但对于比较大的周期性体系，虽然可以用GaussView观看振动模式，但作图效果可能不够理想。

https://github.com/Stanford-MCTG/VASP-plot-modes提供了一个将VASP的振动模式转化为.vesta格式的脚本，我把这个脚本和@北大-陶豫 师兄的从cp2k拉曼计算的.out和.mol文件中读取振动信息的脚本结合起来，用于将cp2k的振动模式转化为.vesta文件，从而实现比较好的作图效果。


输入文件为：cp2k振动分析产生.out和.mol文件，以及vesta模板文件template.vesta（将结构用VESTA打开，保存为.vesta文件即可，注意原子顺序需要和cp2k计算的结构的原子顺序一致）
输出文件为：mode_第几个模式_频率.vseta文件，将该文件用VESTA打开，振动方向用箭头显示，箭头的长短表示振幅，不过不能动
脚本中保留了读取红外、拉曼活性的部分，这样方便通过对mask_mode函数进行修改，从而选择特定模式进行观看

效果如下, 这里只用了CsPbI3的原胞作为示意，实际计算需要建立超胞。

北大-陶豫

很棒的脚本！感觉这个脚本值一个威望

轻舟~

想问一下这个脚本只能在LInux系统上用吗？我用vccode运行这个脚本，显示'"ls *.mol"' 不是内部或外部命令，也不是可运行的程序或批处理文件“，Windows上是用不了嘛？