VASP DFPT方法计算声子谱

第一步：精确优化结构 EDIFF  =1E-8 EDIFFG=1E-5 增大Ecut

第二步：拿第一步优化出来的结构命名为POSCAR-unitcell计算声子谱 INCAR设置：IBRION = 8  NSW=1

第三步：扩胞 phonopy -d --dim="2 2 2" -c POSCAR-unitcell

mv SPOSCAR POSCAR 以上一步生成的 SPOSCAR 为 POSCAR，准备 KPOINTS、POTCAR 和 INCAR，进行 vasp 计算

第四步：计算完成后 收集力常数 用 vasp 计算得到的 vasprun.xml 生成力常数文件 FORCE_CONSTANTS

phonopy --fc vasprun.xml

第五步：准备band.conf文件：

ATOM_NAME = Si

DIM = 2 2 2 (自行修改)

FORCE_CONSTANTS = read

BAND = 0 0 0   0.5 0 0.5   1 0 0   0.5 0 0   0 0 0(自行修改)

第六步：phonopy --dim="2 2 2" -c POSCAR-unitcell --factor=521.471 band.conf 将频率单位转换为 cm-1

第七步：bandplot --gnuplot band.yaml > band.dat 用生成的band.dat文件作图