内容发布更新时间 : 2025/2/6 22:58:33星期一 下面是文章的全部内容请认真阅读。

排错过程：以SB-3QL体系为例，计算薄膜材料的自旋轨道耦合能带图

SYSTEM=Bi2Se3relax ENCUT=340 ISTART=0 ICHARG=2 ISMEAR=0 SIGMA=0.1 #ELECTRONIC PARAMETERS NELM=200 NELMIN=2 EDIFF=1E-6 IALGO=48 #Atom location potimization IBRION=2 NSW=200 EDIFFG=-0.01 ISIF=4 PREC=Accurate LREAL=.TRUE. NPAR=1 AMIN=0.01 LCHARG=.F. LWAVE=.F. 结构优化出现错误： Internal内部的、内在的 ERROR RSPHER:running out of buffer缓冲 0 0 13 1 0 nonlr.F:Out of buffer RSPHER 将NPAR=1修改成4（或者2），问题得以解决。 SYSTEM=Bi2Se3relax ENCUT=340 ISTART=0 ICHARG=2 ISMEAR=0 SIGMA=0.1 #ELECTRONIC PARAMETERS NELM=200 NELMIN=2 EDIFF=1E-6 ALGO=Fast #Atom location potimization IBRION=2 NSW=200 EDIFFG=-0.01 ISIF=4 SYSTEM=Bi2Se3scf ENCUT=340 ISTART=0 ICHARG=2 ISMEAR=0 SIGMA=0.1 #ELECTRONIC PARAMETERS NELM=200 NELMIN=2 ALGO=Fast EDIFF=1E-6 LELF=.F. PREC=A LREAL=.False. LCHARG=.T. LWAVE=.F. SYSTEM=Bi2Se3bands ENCUT=340 ISTART=0 ICHARG=11 ISMEAR=0 SIGMA=0.01 #ELECTRONIC PARAMETERS NELM=200 NELMIN=2 EDIFF=1E-6 LELF=.F. PREC=A ALGO=Fast LREAL=.False. NPAR=1 LORBIT=11 SYSTEM=Bi2Se3scf ENCUT=340 ISTART=0 ICHARG=2 ISMEAR=0 SIGMA=0.1 #ELECTRONIC PARAMETERS NELM=200 NELMIN=2 IALGO=48 EDIFF=1E-6 LELF=.F. PREC=A LREAL=.TRUE. LCHARG=.T. LWAVE=.F. NPAR=1 LORBIT=11 AMIN=0.01 ISPIN=2 LSORBIT=.TRUE. SYSTEM=Bi2Se3bands ENCUT=340 ISTART=0 ICHARG=11 ISMEAR=0 SIGMA=0.01 #ELECTRONIC PARAMETERS NELM=200 NELMIN=2 EDIFF=1E-6 LELF=.F. PREC=A IALGO=48 LREAL=.True. NPAR=1 LORBIT=11 AMIN=0.01 LCHARG=.F. LWAVE=.F. LSORBIN=.TRUE. ISYM=0 GGA_COMPAT=.FALSE. PREC=Accurate LREAL=.FALSE. NPAR=1 AMIN=0.01 LCHARG=.F. LWAVE=.F. 结构优化出现错误： NPAR=1 LORBIT=11 AMIN=0.01 ISPIN=2 LSORBIT=.TRUE. AMIN=0.01 LCHARG=.F. LWAVE=.F. LSORBIN=.TRUE. ISYM=0 GGA_COMPAT=.FALSE. WARNING: Sub-Space-Matrix is not hermitian in DAV 4-4.681828688433112E-002 通过将默认AMIX=0.4，修改成AMIX=0.2（或0.3），问题得以解决。 但是，最终得到的能带结构图，和文献并不相符。 方案一：我们修改原来的情况1-7中参数，并验证POTCAR的影响，进行计算（3SB 1-7---12-17（d--INCAR LMAXMIX=4））： 方案二：算Bi2Se3和Bi2Te3薄膜1-6QL，S1-S6，T1-T6 (1-6QL-SB-TB) 方案三：针对KPOINTS，取K点6-15，做一个测试，计算其是否对能带结构有影响。（3SB-1） 方案四：我们基于已经得出正确结果的Bi2Se3块体结果，算Bi2Te3块体结果，并算一系列块体的结果。(Bulk) 方案一： SYSTEM=Bi2Se3relax ENCUT=340 ISTART=0 ICHARG=2 ISMEAR=0 SIGMA=0.1 #ELECTRONIC PARAMETERS NELM=200 NELMIN=2 EDIFF=1E-6 #ALGO=Fast IALGO=48 #Atom location potimization IBRION=2 NSW=200 EDIFFG=-0.01 ISIF=4 PREC=Accurate LREAL=.TRUE. NPAR=4 AMIN=0.01 LCHARG=.F. LWAVE=.F. AMIX=0.2 参数设置： (一) IALGO=48 LREAL=.TRUE. 661 结构 a=4.186 c=74.380（50?） Se1 (0.333 0.667 0.312) Se2 (0.000 0.000 0.688) SYSTEM=Bi2Se3scf ENCUT=340 ISTART=0 ICHARG=2 ISMEAR=0 SIGMA=0.1 #ELECTRONIC PARAMETERS NELM=200 NELMIN=2 #ALGO=Fast IALGO=48 EDIFF=1E-6 LELF=.F. PREC=A LREAL=.TRUE. LCHARG=.T. LWAVE=.F. NPAR=2 LORBIT=11 AMIN=0.01 ISPIN=2 LSORBIT=.TRUE. AMIX=0.2 SYSTEM=Bi2Se3bands ENCUT=340 ISTART=1 ICHARG=11 ISMEAR=0 SIGMA=0.1 #ELECTRONIC PARAMETERS NELM=200 NELMIN=2 EDIFF=1E-6 LELF=.F. PREC=A #ALGO=Fast IALGO=48 LREAL=.True. NPAR=2 LORBIT=11 AMIN=0.01 LCHARG=.F. LWAVE=.F. LSORBIN=.TRUE. ISYM=0 GGA_COMPAT=.FALSE. AMIX=0.2 能带 ISMEAR=0 SIGMA=0.1 Scf: ISPIN=2 LSORBIT=.TRUE. Bands： LSORBIN=.TRUE. ISYM=0 GGA_COMPAT=.FALSE. （二） IALG=Fast LREAL=.FALSE. ISMEAR=0 SIGMA=0.1 Scf: ISPIN=2 LSORBIT=.TRUE. Bands： ISYM=0 LSORBIN=.TRUE. GGA_COMPAT=.FALSE. (三) IALGO=48 LREAL=.FALSE. Scf: ISPIN=2 ISMEAR=-5 LSORBIT=.TRUE. Bands： ISYM=0 LSORBIN=.TRUE. GGA_COMPAT=.FALSE. (四) IALGO=48 LREAL=.FALSE. Scf: ISMEAR=0 SIGMA=0.1 LSORBIT=.TRUE. Bands： ISYM=0 LSORBIN=.TRUE. GGA_COMPAT=.FALSE. (五) IALGO=48 LREAL=.FALSE. Scf: ISPIN=2 LSORBIT=.TRUE. LORBMOM=.TRUE. ISMEAR=0 SIGMA=0.1 Bands： ISYM=0 LSORBIN=.TRUE. LORBMOM=.TRUE. GGA_COMPAT=.FALSE. (六) IALGO=48 LREAL=.FALSE. Scf: ISPIN=2 LSORBIT=.TRUE. ISMEAR=0 SIGMA=0.1 Bands： LSORBIN=.TRUE. GGA_COMPAT=.FALSE. 771 a=4.187 c=74.316 同 881 a=4.192 c=74.161 同 991 a=4.192 c=74.160 同 (高能区有尖峰，不同K点能带图都一样) 10101 a=4.192 c= 74.161 同 11111 10101 a=4.187 c= 74.335 同 a=5.194 c= 48.301 Se1 (0.333 0.667 0.321) Se2 (0.000 0.000 0.679) 真空层 20? 未算完 10101 a=4.495 c= 48.361 Se1 (0.333 0.667 0.328) Se2 (0.000 0.000 0.672) 真空层 20? 10101 a=4.192 c= 84.959 Se1 (0.333 0.667 0.327) Se2 (0.000 0.000 0.459) 真空层 50?

10101 a=4.192 c= 74.161 Se1 (0.333 0.667 0.312) Se2 (0.000 0.000 0.688) 真空层 50? 10101 a=4.192 c= 74.161 Se1 (0.333 0.667 0.312) Se2 (0.000 0.000 0.688) 真空层 50? (七) IALGO=48 LREAL=.FALSE. 10101 a=4.191 c= 74.161