<div dir="ltr">N.B., there can also be a discontinuity in the charge (small) due to the tails of the core states which can be eliminated by doing "touch .lcore".</div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 21, 2016 at 8:36 AM, Laurence Marks <span dir="ltr"><<a href="mailto:L-marks@northwestern.edu" target="_blank">L-marks@northwestern.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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" target="_blank">http://www.numis.northwestern.<wbr>edu</a><br>
Corrosion in 4D <a href="http://MURI4D.numis.northwestern.edu" target="_blank">http://MURI4D.numis.<wbr>northwestern.edu</a><br>
Partner of the CFW 100% gender equity project, <a href="http://www.cfw.org/100-percent" target="_blank">www.cfw.org/100-percent</a><br>
Co-Editor, Acta Cryst A<br>
    </p><div class="HOEnZb"><div class="h5">
<div class="gmail_extra"><br><div class="gmail_quote">On Nov 21, 2016 07:39, "John Rundgren" <<a href="mailto:jru@kth.se" target="_blank">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>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="font-size:12.8px">Professor Laurence Marks</span><br></div><div dir="ltr"><span style="font-size:12.8px">"Research is to see what everybody else has seen, and to think what nobody else has thought", </span><span style="font-size:12.8px">Albert Szent-Gyorgi</span><br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a> ; <span style="font-size:12.8px">Corrosion in 4D: </span><a href="http://MURI4D.numis.northwestern.edu" style="font-size:12.8px" target="_blank">MURI4D.numis.northwestern.edu</a><div><span style="font-size:12.8px">Partner of the CFW 100% program for gender equity, </span><a href="http://www.cfw.org/100-percent" style="font-size:12.8px" target="_blank">www.cfw.org/100-percent</a></div><div>Co-Editor, Acta Cryst A</div></div></div></div></div></div></div></div></div></div></div>
</div>