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

愿有所成

IEchoQ

引喻失义   妄自菲薄

IEchoQ