<p dir="ltr">APW+lo methods have a step in the gradient of the density at the RMT. To avoid this use a lapw basis set: to reduce it increase RKMAX.</p>
<p dir="ltr">---<br>
Professor Laurence Marks<br>
"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Gyorgi<br>
<a href="http://www.numis.northwestern.edu">http://www.numis.northwestern.edu</a><br>
Corrosion in 4D <a href="http://MURI4D.numis.northwestern.edu">http://MURI4D.numis.northwestern.edu</a><br>
Partner of the CFW 100% gender equity project, <a href="http://www.cfw.org/100-percent">www.cfw.org/100-percent</a><br>
Co-Editor, Acta Cryst A<br>
    </p>
<div class="gmail_extra"><br><div class="gmail_quote">On Nov 21, 2016 07:39, "John Rundgren" <<a href="mailto:jru@kth.se">jru@kth.se</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div dir="ltr" style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Dear Peter Blaha and Gavin Abo,<br>
<br>
Non-overlapping muffin-tin spheres are used by WIEN2k and my LEED program eeasisss (Elastic Electron-Atom Scattering in Solids and Surface Slabs). But RMT(setrmt_lapw) is not automatically the best choice for RMT(eeasisss). LEED touching radii of atoms depend
 on exchange-correlation interaction between crystal electron gas (the WIEN2k electrons) and an incident LEED electron (energy 20-500 eV).<br>
<br>
This is a N+1 electron scattering situation, where "N" signifies the WIEN2k electrons and "1" an alien LEED electron.<br>
<br>
W2k can be reconciled with LEED using an atomic sphere approximation (ASA) extending into the Fourier expansion realm of W2k. A while ago you (P.B. and G.A.) suggested an ASA routine, in which I now use Poisson differentiation of vcoul_ASA in order to obtain
 clmsum_ASA. I consider the case LM=(0,0), sufficient for current LEED.<br>
<br>
The considered structure is a supercell = a surface slab 15 layers thick, where layers 1-2 and 14-15 are C-O and O-C, respectively. Mirror symmetry about layer 8. At the C-O layers vcoul_ASA(0,0) is continuous across the RMT radius, but clmsum_ASA(0,0) versus
 radius shows a step of the order of 10%.<br>
<br>
Is the step k-point dependent? It does not seem so. With 16 and 48 k-points the clmsum_ASA(0,0) steps are preserved within 6 digits.<br>
<br>
I shall be glad to supply the code. When the described numerical error is fixed, WIEN2k and eeasisss can be re-run self-consistently within the model of non-overlapping muffin-tin atoms.<br>
<br>
Regards,<br>
John Rundgren</p>
<p>KTH<br>
<br>
</p>
<p><br>
</p>
<p><br>
</p>
</div>

</blockquote></div></div>