前面我们学习了使用Phonopy，通过有限位移法和密度泛函微扰法(DFPT)计算得到声子谱以及对声子谱进行非解析项校正。那么如何得到声子相关的其他性质，比如自由能、热容、熵等热学性质呢？本文将简单介绍在DFPT计算声子谱的基础上，如何通过Phonopy后处理得到自由能、热容和熵等热学性质相关数据并绘图。

操作步骤

(1) 首先对结构扩胞，采用DFPT法进行声子谱计算，得到的声子谱没有虚频；

(2) 获得FORCE_CONSTANTS文件基础上，准备mesh.conf文件，文件内容如下；

       ATOM_NAME = Al N

       DIM = 3 3 3   #扩胞倍数

       MP = 10 10 10

       FORCE_CONSTANTS = READ

       TPROP = T

       TMIN = 200   #起始温度

       TMAX = 2000  #终止温度

       TSTEP = 50   #温度梯度为50

(3) 输入如下命令，即可获得thermal_properties.yaml文件；

       phonopy -t mesh.conf -c  POSCAR-unitcell

(4) 为了绘图方便，接着输入如下命令，将thermal_properties.yaml文件转换成thermal.dat文件；

       phonopy-propplot --gnuplot thermal_properties.yaml > thermal.dat

(5) 用Notepad++打开thermal.dat文件即可发现，第一列数据为温度，单位K，温度梯度为50，第二列数据为自由能，单位KJ/mol，第三列数据为熵，单位J/K/mol，第四列数据为热容，单位为J/K/mol；

(6) 将获得的thermal.dat文件导入Origin进行绘图即可。

结果如图所示：