<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:12.499999046325684px">David,</span><div style="font-family:arial,sans-serif;font-size:12.499999046325684px"><br></div><div style="font-family:arial,sans-serif;font-size:12.499999046325684px">
Thanks.  Yes, I later realized that with zero localization step the (U(dis)(k))^dagg.H(k).U(dis)(k)  essentially projects the bloch-states onto the trial basis within the chosen energy window (and also diagonalises the H in that subspace by construction), which as you mention can include unoccupied bands as well.  This is what I needed.</div>
<div style="font-family:arial,sans-serif;font-size:12.499999046325684px"><br></div><div style="font-family:arial,sans-serif;font-size:12.499999046325684px">Ganesh</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Fri, Jan 25, 2013 at 12:07 PM, David Vanderbilt <span dir="ltr"><<a href="mailto:dhv@physics.rutgers.edu" target="_blank">dhv@physics.rutgers.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Ganesh,<br>
<br>
Regarding this part<div class="im"><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Does anyone know if it is possible to transform the current<br>
h(r) matrix elements, which is in the WF-basis, to (say<br>
an atom centered) atomic-orbital basis within wannier90?<br>
The transformation might be lossy (i.e. it might be a projection)<br>
depending on the system of course, but would be good to compare<br>
with other TB-fitting codes which use such atomic-orbital basis<br>
to represent the Hamiltonian.<br>
</blockquote>
<br></div>
To do this, we normally choose trial functions that have the<br>
desired chemical symmetry, and then do zero iterations of the<br>
maxloc procedure.  That is, we use the WFs that are obtained<br>
just by projection on the trial functions.  We also include<br>
unoccupied bands if needed to "complete the shell" of WFs on<br>
each atom.  Because we skip the iterative localization, these<br>
may not be quite as localized as the MLWFs, but they may be<br>
good for your purposes.<br>
<br>
David<br>
<br>
------------------------------<u></u>------------------------------<u></u>------------<br>
Prof. David Vanderbilt               Phone: <a href="tel:%28732%29%20445-5500%20x2514" value="+17324455500" target="_blank">(732) 445-5500 x2514</a><br>
Department of Physics and Astronomy  Fax:   <a href="tel:%28732%29%20445-4400" value="+17324454400" target="_blank">(732) 445-4400</a><br>
Rutgers University                   Email:  <a href="mailto:dhv@physics.rutgers.edu" target="_blank">dhv@physics.rutgers.edu</a><br>
136 Frelinghuysen Road               <a href="http://www.physics.rutgers.edu/~dhv" target="_blank">http://www.physics.rutgers.<u></u>edu/~dhv</a><br>
Piscataway, NJ 08854-8019   USA<br>
------------------------------<u></u>------------------------------<u></u>------------<br>
</blockquote></div><br></div>