<div dir="ltr">Dear WIEN2k users,<br><br>I am calculating the band gap of bulk BaTiO3 using the LDA+U with PBE functional. For the values below, the band gap was 2.002 eV (which I can explain by the choice of the PBE functional). When I moved to a supercell with a defect (oxygen vacancy, supercell dimensions 2x2x1), the band gap became 0.0 eV.<br>I tried to increase the Hubbard parameter U. For example, for bulk calculations with U=1.121 Ry = 15.252 eV, the band gap was 3.391 eV. However, running with the same parameters for a supercell with a defect gave a calculated value of 0.889 eV.<br>I have the following questions:<br>1. Why does the band gap decrease so much when going from bulk to supercell with a defect?<br>2. Is it worth changing the PBE functional to another one if the band gap for bulk is less than the experimental one?<br>3. If I continue the calculations on PBE+U, how can I improve the values?<div><br></div><div>With Best Regards,</div><div>Natalia</div><div><br></div></div>