
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none"><!--P{margin-top:0;margin-bottom:0;} --></style>

</head>

<body 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>

</body>

</html>