<div dir="ltr">Thank you very much</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Sep 13, 2023 at 8:32 PM Igor I Mazin <<a href="mailto:imazin2@gmu.edu">imazin2@gmu.edu</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">1. Undoped BaTiO3 is an insulator. An oxygen vacancy introduces electron <br>
doping. When you dope an insulator, you get a metal. In a metal, the gap <br>
is zero.<br>
2. U=15 eV is an insanely large and unphysical value. You are better off <br>
with no U at all. For small doping (small concentration of vacancies) <br>
Hubbard correlations are irrelevant.<br>
3. 2x2x1 with one vacancy is doping of 0.5 e per formula unit. This is <br>
unphysically large (like, by three orders of magnitude) and I would not <br>
trust a single number from this simulation. Check <br>
<a href="https://journals.aps.org/prresearch/abstract/10.1103/PhysRevResearch.5.033080" rel="noreferrer" target="_blank">https://journals.aps.org/prresearch/abstract/10.1103/PhysRevResearch.5.033080</a><br>
<br>
On 9/13/2023 13:19, Natalia Andreeva wrote:<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 can <br>
> 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 defect <br>
> 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>
Igor Mazin, Prof. of Advanced Studies<br>
Quantum Science and Engineering Center<br>
Department of Physics and Astronomy<br>
George Mason University<br>
phone 1-703-503-8152 (h)<br>
<a href="http://mason.gmu.edu/~imazin2" rel="noreferrer" target="_blank">http://mason.gmu.edu/~imazin2</a><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>