經(jīng)驗分享

在使用VASP+加phonopy計算材料聲子色散結(jié)構(gòu)時，可以考慮激活非分析項矯正，這時需要增加計算的工作量以獲得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，然后繼續(xù)計算

• LEPSILON = .TRUE.

在vasp的計算任務完成后，通過phonopy的代碼可直接獲得BORN文件，命令如下

• outcar-born >BORN #舊版

• phonopy-vasp-born>BORN #新版

在獲得BORN文件后，在處理原子間力常數(shù)獲得聲子色散時的命令中加入--nac即可，如手冊中的命令案例

• phonopy -p --nac band.conf

在phonopy手冊中給出考慮和不考慮LO-TO劈裂的NaCl的聲子譜的對比，分別如下：

不考慮LO-TO劈裂

考慮LO-TO劈裂

愿有所成

IEchoQ

引喻失義 ? 妄自菲薄