DFPT算声子

DFPT（vasp5.1以后可用）

（DFPT应该可以不扩胞，但学术之友给的示例扩胞了或许会更准确？）

• 自洽的INCAR POSCAR KPOINTS POTCAR

• 改INCAR中IBRION=8 #7/8即算DFPT得到Hessian矩阵，原子多的话增加SOC算贼慢

• 故我体系计算声子未添加SOC

phonopy --fc vasprun.xml #生成力常数文件FORCE_CONSTRAINTS

• 编辑band.conf文件，给出高对称点路径信息

ATOM_NAME = XXXX
DIM = m n l #建立超胞的大小，像我dfpt没扩胞就是1 1 1
BAND = 0.5 0.0 0.0  0.0 0.0 0.0  0.333333 0.333333 0.0  0.5 0.0 0.0  #高对称点路径，每组高对称点之间以两个空格隔开
FORCE_CONSTANTS = READ #表示读取力常数文件FORCE_CONSTANTS，我之前不知道在哪儿看有写write的我试了后就失败了 嗐

*phonopy软件默认在两个高对称点之间打点51个，且在PBAND.dat中每组高对称点数据以一个空行分割，据此即可得到完整的声子谱图像

• 生成band.yaml文件

phonopy --dim="m n l" -c POSCAR-unitcell band.conf #后加-f -s可保存为pdf

• 得到声子谱数字文件PBAND.dat文件，用gnuplot画图

phonopy-bandplot  --gnuplot> PBAND.dat

or phonopy-bandplot --gnuplot> phonon.out

本文参考贺勇前辈文章：https://yh-phys.github.io/2019/10/11/vasp-phonon/