DFPT计算报错internal error in SETUP_DEG_CLUSTERS:NB_NOTexceeds解决办法

经验分享

在初学计算体系声子谱的时候会遇到诸多问题，这里简单列举其中一个以及解决办法。

笔者在计算某面心立方晶体材料时，使用原胞扩包4*4*4进行声子谱计算，最初采用DFPT方法，对phonopy扩包后的SPOSCAR进行计算，在走完第一步离子步后出现报错

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某些设置等进行调整发现仍不能解决问题。

最后我们改用有限位移法进行声子谱的计算。使用的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.

注意在后续处理数据时，使用如下命令得到FORCE_SETS文件

 phonopy -f ./d*/vasprun.xml

其中d*代表的是所有有限位移计算的文件夹。在处理声子谱数据时在band.conf中也要删掉FORCE_CONSTANTS= READ，并添加上下面内容

 FORCE_SETS = READ FULL_FORCE_CONSTANTS = .TRUE. FORCE_CONSTANTS= WRITE #生成FORCE_CONSTANTS

愿有所成

引喻失义   妄自菲薄