<p dir="ltr">I do try and trap various odd cases during a normal scf run, and have never thought of just using mixer as a conversion so it might be broken by a trap :-(. I am pretty sure I put something in for the next version for some other reasons that probably breaks this; there might already be something there in the current version. (If there are problems with input files during an scf run the mixer switches to Pratt as a last resort. A full Pratt step could then be very bad.)</p>
<p dir="ltr">Copying the case.clmval to case.vrespval and using -vresp (I.e. x mixer_vresp) may be safest as this is forced to a full Pratt step with no renormalization. The copy case.vrespsum to case.clmsum or edit the aim.def file.</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 8:20 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">I'm pretty sure mixer works that way (unless Lauri has broken it<br>
recently ?):<br>
<br>
when you provide just a case.clmval file (no clmcor, no clmsum_old !!!)<br>
it cannot "mix" anything, just transfer clmval into clmsum.<br>
<br>
And to switch off renormalization, put "NO" for NORM in case.inm.<br>
<br>
<br>
On 02/19/2015 03:06 PM, Laurence Marks wrote:<br>
> A small comment on what Peter suggested. It would not be hard to modify<br>
> aim to use a clmval since the code has similarities to lapw5 where that<br>
> is already done - look for the switch in lapw5 then do something similar<br>
> in aim. (I am slightly surprised it is not already there.)<br>
><br>
> Using filtvec or lapw2 with energy limits you may then be able to come<br>
> close to what you want, but perhaps not exactly what you are looking for.<br>
><br>
> I am not certain about Peter's idea of fiddling the input to mixer to do<br>
> the normalization -- it might work or could just crash. You probably<br>
> need to create case.clmcor files with nothing in them by changing<br>
> case.inc so it has no electrons. Then force PRATT 1.0 to be used (but<br>
> this might be trapped). Some fiddling needed here.<br>
><br>
> You may also come close by using filtvec or limited energy ranges and<br>
> look at the :QTL numbers from lapw2 or qtl.<br>
><br>
> Or some combination of filtvec/energy ranges with lapw2, then lapw5 and<br>
> critc2.<br>
><br>
> Or switch to some orbital code, perhaps the best option if you really<br>
> need the numbers.<br>
><br>
> ___________________________<br>
> Professor Laurence Marks<br>
> Department of Materials Science and Engineering<br>
> Northwestern University<br>
> <a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a> <<a href="http://www.numis.northwestern.edu" target="_blank">http://www.numis.northwestern.edu</a>><br>
> <a href="http://MURI4D.numis.northwestern.edu" target="_blank">MURI4D.numis.northwestern.edu</a> <<a href="http://MURI4D.numis.northwestern.edu" target="_blank">http://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<br>
> nobody else has thought"<br>
> Albert Szent-Gyorgi<br>
><br>
> On Feb 19, 2015 7:05 AM, "Peter Blaha" <<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a><br>
> <mailto:<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a>>> wrote:<br>
><br>
> 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<br>
> inverse<br>
> > participation ratio. For that I need to get the Mulliken point<br>
> 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<br>
> someone<br>
> > advice me how?<br>
> > Also if possible I would like to get it inside atomic basins<br>
> 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>
> <mailto:<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:<br>
> <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><br>
> <mailto:<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> <mailto:<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:<br>
> <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>
> 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>