<div dir="ltr"><div>Dear Eric,</div><div><br></div><div>just a couple of quick comments. <br></div><div>If you are already using a large supercell you can probably avoid supercells (all this would need to be checked to be sure).</div><div>Also keep in mind that from version 6.4.1 you can find in the QE distribution a code called HP (Hubbard parameters) that computes Hubbard U automatically using DFPT. More details here:</div><div>10.1103/PhysRevB.98.085127.</div><div>Hope this helps.</div><div>Best regards,</div><div><br></div><div>Matteo<br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mer 26 giu 2019 alle ore 21:51 Eric Glen Suter <<a href="mailto:esuter@uga.edu">esuter@uga.edu</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">




<div dir="ltr">
<div id="gmail-m_5236249242120536584divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p style="margin-top:0px;margin-bottom:0px">Hello all,</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">I'm trying to simulate transition metal dopants in particular charge states. <span style="font-size:12pt">I'm more or less familiar with the "tot_charge" tag and how it functions. For these transition metals, it seems
 the extra charge I introduce doesn't want to localize on my dopants. I've been looking into using the Hubbard U as a way to coerce this extra electron (or hole in some cases</span><span style="font-size:12pt">) to localize on my dopant. I've read up a bit
 on the literature about these types of corrections, but I'd like to know:</span></p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">1) Am I on the right track with this line of reasoning? I also take notice of the "starting_charge" tag for pw.x. I haven't tried it out, but is that a viable way of trying to coax the extra charge to stay on the dopant?</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">2) If Hubbard U really is the way to go, is there a way of tuning the U parameter that doesn't involve supercells? I'm already working with a pretty large cell and I don't have a ton of computational power at my disposal. </p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">I understand this is no trivial task, but if anyone has any insight on how to model a dopant in a particular charge state, I'd be grateful for a nudge in the right direction.</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Thank you in advance for any help you might have to offer,</p>
<p style="margin-top:0px;margin-bottom:0px">Eric Suter</p>
<p style="margin-top:0px;margin-bottom:0px">Dept. of Physics and Astronomy</p>
<p style="margin-top:0px;margin-bottom:0px">University of Georgia</p>
</div>
</div>

_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>