phonopy计算声子谱时考虑LO-TO劈裂

经验分享

在使用VASP+加phonopy计算材料声子色散结构时，可以考虑激活非分析项矫正，这时需要增加计算的工作量以获得BORN文件。以此计算获得的声子色散则会考虑LO-TO的劈裂效果。

在phonopy的手册中提到了可以通过自洽计算获得BORN，手册中推荐的INCAR为

• PREC = Accurate

• IBRION = -1

• NELMIN = 5

• ENCUT = 500

• EDIFF = 1.000000e-08

• ISMEAR = 0

• SIGMA = 1.000000e-02

• IALGO = 38

• LREAL = .FALSE.

• LWAVE = .FALSE.

• LCHARG = .FALSE.

• LEPSILON = .TRUE.

这是在使用有限位移法也就是冻声子法时，增加的额外的计算任务。当使用DFPT方法计算原子受力时，可直接在INCAR中添加下面这条tag，然后继续计算

• LEPSILON = .TRUE.

在vasp的计算任务完成后，通过phonopy的代码可直接获得BORN文件，命令如下

• outcar-born >BORN #旧版

• phonopy-vasp-born>BORN #新版

在获得BORN文件后，在处理原子间力常数获得声子色散时的命令中加入--nac即可，如手册中的命令案例

• phonopy -p --nac band.conf

在phonopy手册中给出考虑和不考虑LO-TO劈裂的NaCl的声子谱的对比，分别如下：

不考虑LO-TO劈裂

​

考虑LO-TO劈裂

​

愿有所成

引喻失义 妄自菲薄