Hello everyone. I have a semiconductor, BeSiP2, which has a band gap slightly above 1eV in the GGA approximation (https://materialsproject.org/materials/mp-1009085/).<div><br></div><div>I calculate a vacancy, V_{P}^{+3}, first using diagonalization. Which produces a band gap around 1eV, indicating that the charge does not delocalize, and so OT should still be appropriate. When I switch to OT, however, there is a negative homo/lumo gap and I can see it in the eigenvalues:</div><div><br></div><div><p>0.17131297      0.17154145      0.18362819      0.18519005</p>
<p> Fermi Energy [eV] :    5.039278</p>
<p>  Lowest Eigenvalues of the unoccupied subspace spin            1</p>
<p> -----------------------------------------------------</p>
<p>  Reached convergence in            1  iterations</p>
<p>       0.18499549      0.22295666      0.23850766      0.25228060</p><p><br></p><p>This is bizarre to me. OT is usually very reliable for me. So this is a strange inconsistency, and maybe its due to how I'm running the calculation. For OT I use the full_all preconditioner, and the rest defaults. For diagonalization I the section at the bottom. Anyone have any guesses as to why there is this discrepancy? </p><p>-Nick</p><p>&SCF</p><p>    MAX_SCF 200</p><p>    added_mos 100</p><p>    eps_scf 1e-6</p><p>    max_diis 15</p><p>    &diagonalization</p><p>    &end</p><p>    &mixing</p><p>    &end</p><p>    &smear</p><p>        elec_temp 300</p><p>        method fermi_dirac </p><p>    &end</p><p>&end</p></div>