第一步:精确优化结构 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文件作图
联系客服