經(jīng)驗分享

在初學(xué)計算體系聲子譜的時候會遇到諸多問題，這里簡單列舉其中一個以及解決辦法。

筆者在計算某面心立方晶體材料時，使用原胞擴包4*4*4進行聲子譜計算，最初采用DFPT方法，對phonopy擴包后的SPOSCAR進行計算，在走完第一步離子步后出現(xiàn)報錯

Linear response G [H, r] |phi>, progress :

? Direction:? ?1

? ? ? ?N? ? ? ?E? ? ? ? ? ? ? ? ? ? ?dE? ? ? ? ? ? ?d eps? ? ? ?ncg? ? ?rms

?internal error in SETUP_DEG_CLUSTERS: NB_TOT exceeds NMAX_DEG

? ? increase NMAX_DEG to? ? ? ? ?233

通過對K點、INCAR某些設(shè)置等進行調(diào)整發(fā)現(xiàn)仍不能解決問題。

最后我們改用有限位移法進行聲子譜的計算。使用的INCAR如下，計算過程沒有報錯，并成功得到聲子譜。

?PREC = Accurate

?IBRION = -1

?ENCUT = 500

?EDIFF = 1.0e-08

?EDIFFG = -0.00001

?ISMEAR = 0

?SIGMA = 0.05

?IALGO = 38

?LREAL = .FALSE.

?LWAVE = .FALSE.

?LCHARG = .FALSE.

注意在后續(xù)處理數(shù)據(jù)時，使用如下命令得到FORCE_SETS文件

phonopy -f ./d*/vasprun.xml

其中d*代表的是所有有限位移計算的文件夾。在處理聲子譜數(shù)據(jù)時在band.conf中也要刪掉FORCE_CONSTANTS= READ，并添加上下面內(nèi)容

?FORCE_SETS = READ

?FULL_FORCE_CONSTANTS = .TRUE.

?FORCE_CONSTANTS= WRITE #生成FORCE_CONSTANTS

愿有所成

IEchoQ

引喻失義? ?妄自菲薄

IEchoQ