<p dir="ltr">A small comment on what Peter suggested. It would not be hard to modify aim to use a clmval since the code has similarities to lapw5 where that is already done - look for the switch in lapw5 then do something similar in aim. (I am slightly surprised it is not already there.) </p>
<p dir="ltr">Using filtvec or lapw2 with energy limits you may then be able to come close to what you want, but perhaps not exactly what you are looking for.</p>
<p dir="ltr">I am not certain about Peter's idea of fiddling the input to mixer to do the normalization -- it might work or could just crash. You probably need to create case.clmcor files with nothing in them by changing case.inc so it has no electrons. Then force PRATT 1.0 to be used (but this might be trapped). Some fiddling needed here.</p>
<p dir="ltr">You may also come close by using filtvec or limited energy ranges and look at the :QTL numbers from lapw2 or qtl.</p>
<p dir="ltr">Or some combination of filtvec/energy ranges with lapw2, then lapw5 and critc2.</p>
<p dir="ltr">Or switch to some orbital code, perhaps the best option if you really need the numbers.</p>
<p dir="ltr">___________________________<br>
Professor Laurence Marks<br>
Department of Materials Science and Engineering<br>
Northwestern University<br>
<a href="http://www.numis.northwestern.edu">www.numis.northwestern.edu</a><br>
<a href="http://MURI4D.numis.northwestern.edu">MURI4D.numis.northwestern.edu</a><br>
Co-Editor, Acta Cryst A<br>
"Research is to see what everybody else has seen, and to think what nobody else has thought"<br>
Albert Szent-Gyorgi</p>
<div class="gmail_quote">On Feb 19, 2015 7:05 AM, "Peter Blaha" <<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">The Mullican charges are not possible within wien2k, since we do not<br>
have "overlapping atomic wavefunctions".<br>
<br>
AIM should be possible with some trick. You have to run in 2 steps, in<br>
the first one you provide the total density and determine the basin, in<br>
the INTEG step you provide a density of a single state. However, the<br>
corresponding case.clmval file cannot be just copied to clmsum, since<br>
its normalization is different. Instead, you have to put it through<br>
"mixer" (without core and clmsum_old densities and without<br>
renormalization) to create a "pseudo-clmsum file.<br>
<br>
In any case, you will never recover eg. a delocalized 4s electron in any<br>
basin ....<br>
<br>
<br>
On 02/19/2015 01:33 PM, Pavel Ondracka wrote:<br>
> Dear Wien2k mailing list,<br>
><br>
> I would like to estimate electron localization by computing the inverse<br>
> participation ratio. For that I need to get the Mulliken point charges<br>
> q_i(E) residing at an atomic site i for an eigenstate with energy<br>
> eigenvalue E.<br>
><br>
> For start I would like to get it inside the atomic spheres, can someone<br>
> advice me how?<br>
> Also if possible I would like to get it inside atomic basins determined<br>
> with AIM. Is this possible, or can I get only the total charge inside<br>
> the basin with AIM?<br>
><br>
> Best regards<br>
> Pavel Ondračka<br>
><br>
> _______________________________________________<br>
> Wien mailing list<br>
> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" 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" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
><br>
<br>
--<br>
<br>
                                       P.Blaha<br>
--------------------------------------------------------------------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300             FAX: +43-1-58801-165982<br>
Email: <a href="mailto:blaha@theochem.tuwien.ac.at">blaha@theochem.tuwien.ac.at</a>    WIEN2k: <a href="http://www.wien2k.at" target="_blank">http://www.wien2k.at</a><br>
WWW:   <a href="http://www.imc.tuwien.ac.at/staff/tc_group_e.php" target="_blank">http://www.imc.tuwien.ac.at/staff/tc_group_e.php</a><br>
--------------------------------------------------------------------------<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" 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" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>