Phonopy 是一个由 python 实现的晶体声子分析程序
它是目前提供了 VASP 的 Wien2k 的接口用来计算原子受力
它的主要功能有: 计算声子色散谱; 计算声子态密度,包括分立态密度; 声子热力学性质,包括自由能,热容量,焓; Phonopy 通过力常数的方法计算声子谱
力常数由计算原子在超晶胞中被移动后的受力得到(Parlinsk-Li-Kawasoe 方法)
同样类型的程序还有 phon, fropho, phonon
其中 phonon 是商业软件,卖的很贵,fropho 和 phonopy 的代码其实都是来自于 phon, fropho 是为了代替 phon 而开发的,目的是为了使用 phon 更方便,phon 和 fropho 主要都是由 fortran 开发的,而现在 fropho 已经停止开发,由 python 开发的 phonopy 代替了 fropho, phonopy 在使用上更为方便,在计算量上更为减少
因为 phon, fropho 和 phonon 在移动原子位置时都是一次只移动一个原子的一个方向,而 phonopy 则可以一次移动一个原子的多个方向,所以和其它程序相比, phonopy 最多可以减少 2/3 的计算量
Phonopy 是使用 python 以及 C 等高级语言编写的用于计算声子谱的程序,它可以很方便地在现有的 unix 或linux 操作系统上进行安装
它利用 VASP 等第一性原理计算程序来计算有限位移下的原子受力,然后利用 phonopy 程序处理原子受力获得声子谱
Phonopy 使用 python 的 matplotlib 等工具包来绘图,因此在使用 Phonopy 绘制声子谱等图形时,要求系统有图形用户界面,对于远程机器,需要安装x-manager 来获得虚拟的图形用户界面,方法如下: 在自己的