<div dir="ltr"><div><div>The introduction of a compensating background is needed because the energy of a charged supercell is infinite. In the case of a charged defect, the supercell of charged defects + compensating background will converge to the isolated defect in the limit of large supercell, but the convergence will be slow. Various techniques (Makov-Payne, Martyna-Tuckerman, ...) can be used to obtain better results for reasonably-sized cells. Note that "tot_charge=+1" means "one electron missing from the system".<br><br></div>I have no reason to think that forces in charged supercells are wrong (not just approximated: wrong). Your results look strange, but it is hard to draw any conclusion from a comparison of two very different approaches and XC functionals. If you can provide your data, maybe we can figure out what happens and why.<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Apr 25, 2018 at 1:03 PM, Aleksandra Oranskaia <span dir="ltr"><<a href="mailto:aleksandra.oranskaia@kaust.edu.sa" target="_blank">aleksandra.oranskaia@kaust.edu.sa</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>Hi dear users and developers of QE,</div><div><br></div><div>If one is interested in studying charged state defects -- is it correct to run supercell optimizations with a compensating jellium background ? </div><div>Let’s say, if one is interested in +1 point defect state -- is it correct to add tot_charge = +1 ?</div><div><br></div><div>I quickly checked if the gradients for molecular cation optimization are reasonable with such an approach. For this I took small cation with known structure and:</div><div>- gaussian09-d.01 optimization with 6-11G**/B3LYP, +1 charge, 1 spin multiplicity gave perfectly matching structure with the real one </div><div>- QE-6.2.1 optimization with vdw-df2/USPP and tot_charge = +1 destroyed the cation, it dissociated; without “jellium” the optimized structure is very close to the real one.</div><div><br></div><div>I would greatly appreciate your advice on how to obtain realistic structures of the supercells with charged defects,</div><div><br></div><br><div>
<div style="word-wrap:break-word"><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">Thanks in advance,<br>Alex.<br>___<br>Aleksandra Oranskaia (M.Sc.)<br>ChemS PhD student, KAUST<br>Phone: +966 50 1335254<br><br><br></div><div style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><br></div><br class="m_6636881564231025686Apple-interchange-newline"></div><br class="m_6636881564231025686Apple-interchange-newline">
</div>
<br></div>
<br>
<div><hr></div><font face="Arial" size="1">This message and its contents, including attachments are intended solely for the original recipient. If you are not the intended recipient or have received this message in error, please notify me immediately and delete this message from your computer system. Any unauthorized use or distribution is prohibited. Please consider the environment before printing this email.</font><br>______________________________<wbr>_________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.<wbr>org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-<wbr>espresso.org/mailman/listinfo/<wbr>users</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div>