<div dir="ltr"><div>Typically the hole (+1) gets delocalized between A and B molecules because of the self interaction error of DFT. If you are lucky, then it might get localized either on A or on B. You can find out where the hole gets localized by plotting the charge density difference between the neutral and the charged system.<br>
<br></div>Prasenjit<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 23 February 2013 21:06, Mahmoud Payami Shabestari <span dir="ltr"><<a href="mailto:mpayami@aeoi.org.ir" target="_blank">mpayami@aeoi.org.ir</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><font style="font-family:Tahoma;font-size:10pt">
<div>Dear QE Developers and Users,</div>
<div> </div>
<div>I have constructed a supercell containing two molecules "A" 
and "B" well separated by, say, 15 Bohrs distance. Now I give 
a total charge +1 for the cell. How the code decides to assign the 
charge +1 to molecule "A" or "B"?</div>
<div>Any comments is highly appreciated.</div>
<div> </div>
<div>Best regards,</div>
<div>
                    
Mahmoud Payami <br>
 </div>
</font></div>
<div> </div>
<div>Condensed Matter Section<br>
Physics Group, AEOI, Tehran-Iran<br>
<br>
Email: <a href="mailto:mpayami@aeoi.org.ir" target="_blank">mpayami@aeoi.org.ir</a></div>
<div>Phone: +98 (0) 21 82064393</div>
<div>Fax:     +98 (0) 21 88221105<br>
---------------------------------------------</div>

<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br>PRASENJIT GHOSH,<br>IISER Pune,<br>First floor, Central Tower, Sai Trinity Building<br>
Garware Circle, Sutarwadi, Pashan<br>Pune, Maharashtra 411021, India<br><br>Phone: +91 (20) 2590 8203<br>Fax: +91 (20) 2589 9790
</div>