使用openmx计算含有两个氧空位的SrTiO3体系没有带隙的问题

freeee

有大佬使用过openmx这个程序吗，我使用openmx计算了含有两个氧空位的SrTiO3体系发现计算出来的总DOS图是没有带隙的，但是我使用原胞进行计算，使用同样的参数计算出来的体系是有带隙的，有大佬知道是什么原因吗？第一张图是原胞的，第二张图是含有两个氧空位的超胞的。这是我的参数#
# File Name
#

System.Name                       openmx
level.of.stdout                   1     # default=1 (1-3)
level.of.fileout                  0     # default=1 (0-2)
DATA.PATH                         /home/user11/tools/openmx3.9/DFT_DATA19
HS.fileout                        off    # on|off, default=off

#
# Definition of Atomic Species
#

Species.Number       3
<Definition.of.Atomic.Species
  O     O6.0-s2p2d1       O_PBE19
  Ti    Ti7.0-s3p2d1      Ti_PBE19
  Sr    Sr10.0-s3p2d2     Sr_PBE19
Definition.of.Atomic.Species>#
# SCF or Electronic System
#

scf.XcType              GGA-PBE        # LDA|LSDA-CA|LSDA-PW|GGA-PBE
scf.SpinPolarization        off        # On|Off|NC
scf.ElectronicTemperature  300       # default=300 (K)
scf.energycutoff           300.0       # default=150 (Ry)
scf.maxIter                 200        # default=40
scf.EigenvalueSolver       Band        # DC|GDC|Cluster|Band
scf.Kgrid                  2 2 2       # means n1 x n2 x n3
scf.Mixing.Type            RMM-DIISK   # Simple|Rmm-Diis|Gr-Pulay
scf.Mixing.History         40
scf.Mixing.EveryPulay      1
scf.Mixing.StartPulay      20
scf.criterion             1.0e-6      # default=1.0e-6 (Hartree)

#
# MD or Geometry Optimization
#

MD.Type                     Nomd       # Nomd|Opt|NVE|NVT_VS|NVT_NH
MD.maxIter                    1        # default=1
MD.TimeStep                   1        # default=0.5 (fs)
MD.Opt.criterion         1.0e-4        # default=1.0e-4 (Hartree/bohr)

#
# DOS and PDOS
#

Dos.fileout                  off        # on|off, default=off
DosGauss.fileout             on
DosGauss.Num.Mesh            500
DosGauss.Width               0.02
Dos.Erange              -25.0  20.0    # default = -20 20
Dos.Kgrid                9 9 9      # default = Kgrid1 Kgrid2 Kgrid3