icehoqion的博客

首先将可计算得到无虚频的vasp输入文件放置在文件夹中，包括POSCAR、INCAR、KPOINTS和POTCAR，以及phonopy处理需要的设置文件，我这里使用vaspkit 305功能产生的KPATH.phonopy文件当作设置文件（请根据你的使用习惯更改脚本）。计算主要流程为:以声子计算无虚频的结构文件为基础，更改其缩放系数，在不同体积下计算其原子间受力、声子结构和热性质。以提取出的e-V的数据为基础，使用phonopy-qha计算包括热膨胀系数、热容和格林艾森参数等在内的热学性质。

声子谱的计算主要目的是获得二阶力常数文件FORCE_CONSTANTS，而为了计算光学性质，同时需要计算波恩电荷，写出BORN文件。它可以从Phonopy或VASP计算红外(IR)强度，和计算Raman活度张量和远离共振近似的标量平均强度。通过irreps.yaml的振动信息，并参考D3点组的字符表，判断可能显示Raman活动的模式，并生成计算文件。按照软件案例，作者可能更喜欢hdf5格式的数据文件，生成force_constants.hdf5文件的命令如下。bands的值为需要计算的模式。

根据笔者和诸多计算材料的研究者的经验，部分结构的材料，尤其是根据实验结果测试获得的材料，往往在低温下使用简谐近似计算得到的声子谱存在一个很难去除虚频的声子谱，但是带有温度的考虑非谐效应的声子谱时则会得到一个无虚频的结果。（DFPT方法，如果使用有限位移法计算的声子谱，得到的力常数文件为FORCE_SETS，可通过DynaPhoPy作者所写的脚本转换为FORCE_CONSTANTS，脚本在DynaPhoPy软件安装包中的example文件夹内）。这里尤其要注意NSW的设置，尽量设置的大一点，依据笔者经验，

笔者在计算某面心立方晶体材料时，使用原胞扩包4*4*4进行声子谱计算，最初采用DFPT方法，对phonopy扩包后的SPOSCAR进行计算，在走完第一步离子步后出现报错。其中d*代表的是所有有限位移计算的文件夹。在处理声子谱数据时在band.conf中也要删掉FORCE_CONSTANTS= READ，并添加上下面内容。使用的INCAR如下，计算过程没有报错，并成功得到声子谱。在初学计算体系声子谱的时候会遇到诸多问题，这里简单列举其中一个以及解决办法。引喻失义 妄自菲薄。