<div dir="ltr">Thanks a lot</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Sep 13, 2023 at 8:56 PM Peter Blaha <<a href="mailto:peter.blaha@tuwien.ac.at" target="_blank">peter.blaha@tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">In addition to what Igor said:<br>
<br>
If you want a reasonable gap for the bulk, you can use   mBJ+U (with a <br>
small U of a few eV).<br>
<br>
Still, for your cell, you won't get an insulator (also doped Si is <br>
"metallic" and conducting, otherwise a computer would not work ...).<br>
<br>
PS: Draw the DOS and compare with the undoped DOS. Maybe you understand <br>
then better.<br>
<br>
<br>
Am 13.09.2023 um 19:19 schrieb Natalia Andreeva:<br>
> Dear WIEN2k users,<br>
><br>
> I am calculating the band gap of bulk BaTiO3 using the LDA+U with PBE <br>
> functional. For the values below, the band gap was 2.002 eV (which I <br>
> can explain by the choice of the PBE functional). When I moved to a <br>
> supercell with a defect (oxygen vacancy, supercell dimensions 2x2x1), <br>
> the band gap became 0.0 eV.<br>
> I tried to increase the Hubbard parameter U. For example, for bulk <br>
> calculations with U=1.121 Ry = 15.252 eV, the band gap was 3.391 eV. <br>
> However, running with the same parameters for a supercell with a <br>
> 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 <br>
> supercell with a defect?<br>
> 2. Is it worth changing the PBE functional to another one if the band <br>
> 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?<br>
><br>
> With Best Regards,<br>
> Natalia<br>
><br>
><br>
> _______________________________________________<br>
> Wien mailing list<br>
> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
<br>
-- <br>
-----------------------------------------------------------------------<br>
Peter Blaha,  Inst. f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-158801165300<br>
Email: <a href="mailto:peter.blaha@tuwien.ac.at" target="_blank">peter.blaha@tuwien.ac.at</a><br>
WWW:   <a href="http://www.imc.tuwien.ac.at" rel="noreferrer" target="_blank">http://www.imc.tuwien.ac.at</a>      WIEN2k: <a href="http://www.wien2k.at" rel="noreferrer" target="_blank">http://www.wien2k.at</a><br>
-------------------------------------------------------------------------<br>
<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div><br clear="all"><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature">Best Regards,<br>Natalia Andreeva.<br></div>