Thanks, Nicola,<br><br>So, the code does not implement the additional operation as a default option then? Then, I might have to ask this way. If you know of the iron pnictide, do you expect the resulting MLWS has the crystal symmetry? Under what condition would the MLWS gain the symmetry of the original lattice? or, maybe the answer is a short no, it is just not guaranteed since mathematically a unitary transformation would not change the nature of the problem, as you mentioned just now. <br>
<br>Sincerely,<br>Jon<br><br><br><div class="gmail_quote">On Fri, May 8, 2009 at 5:36 PM, Nicola Marzari <span dir="ltr"><<a href="mailto:marzari@mit.edu">marzari@mit.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
<br>
Hi Jon,<br>
<br>
<br>
(as usual, do sign your emails, so we know a bit more about you).<br>
<br>
The problem is that the Wannier functions you obtain from minimization<br>
of the spread can be invariant (mathematically, or just numerically)<br>
with respect to certain unitary transformations - so the WFs do not<br>
have the symmetry that you expect.<br>
<br>
For this reason, we suggested in the 2002 PRB to perfrom an additional<br>
transformation (that would not change the spread) that would remove this<br>
degeneracy - for the case of d orbitals studied there is was the<br>
diagonalization of <r^2>.<br>
<br>
The choice of what you diagonalize depends on the symmetry you want to recover, and you might need to construct by hand (working on the code)<br>
the additional matrix elements that are needed.<br>
<br>
I never tried it, but if you were to give trial projections that<br>
already had the riht symmetry, and the right shape, maybe the code<br>
would not move away from that (even if there would be nothing<br>
opposing it - so just numerical build-up could slowly shift you<br>
away).<br>
<br>
nicola<br>
<br>
<br>
<br>
jon wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im">
Dear Prof. Marzari,<br>
<br>
Let me see whether I interpret your paper right. Basically, for a crystal with lower symmetry, MLWS could still have arbitrariness in choosing its form. Thus, introducing a new criterion by minimizing <r^2> could further reduce the arbitrariness, but not fully remove it. So, is this what you mean by "to recover symmetries that are not built in in the second-moment minimization"?<br>
<br>
I guess the then relevant question related to my problem would be: is minimization on <r^2> implemented in the wannier90 code? To help understand the issue better, I would ask two more questions: is this minimization equivalent to attempts to build in lattice symmetry further, or in other words, is the removal of arbitrariness equivalent in some sense to restoration of crystal symmetry for the wannier states? Is there logic behind if this equivalence holds?<br>
<br>
Thanks for your kind help,<br>
<br>
Sincerely,<br>
Jon<br>
<br></div><div><div></div><div class="h5">
On Thu, May 7, 2009 at 12:30 PM, Nicola Marzari <<a href="mailto:marzari@mit.edu" target="_blank">marzari@mit.edu</a> <mailto:<a href="mailto:marzari@mit.edu" target="_blank">marzari@mit.edu</a>>> wrote:<br>
<br>
Nicola Marzari wrote:<br>
<br>
<br>
<br>
Dear Jon,<br>
<br>
<br>
have a look at the Posternak...Marzari 2002 PRB. I presume that<br>
that refers to your problem - you need to diagonalize in the<br>
subspace of<br>
localized orbitals to recover symmetries that are not built in in<br>
the second-moment minimization.<br>
<br>
Let us know if that is the problem that affects you,<br>
<br>
nicola<br>
<br>
<br>
jon wrote:<br>
<br>
Dear Colleagues,<br>
<br>
I am working on DMFT applied to iron pnictide and need to<br>
obtain the impurity Green function, which is basically the<br>
bare Green function being summed over all the momentum<br>
points in 1BZ. The tight-binding hoppings are extracted from<br>
fitting to a band structure calculated through MLWS using<br>
wannier90. These hoppings have cubic symmetry(on 2D plane<br>
though). The impurity Green function turns out to be<br>
diagonal in the orbital indices. This is not generally the<br>
case unless the wannier states have been taken further<br>
actions to form cubic harmonics(I am not an LDA expert, thus<br>
please correct me if I am wrong here) before output.<br>
However, the tutorial or introduction does not explicitly<br>
mention this action and I am thus confused.<br>
<br>
Thanks,<br>
<br>
Sincerely,<br>
Jon<br>
<br>
<br>
------------------------------------------------------------------------<br>
<br>
_______________________________________________<br>
Wannier mailing list<br>
<a href="mailto:Wannier@quantum-espresso.org" target="_blank">Wannier@quantum-espresso.org</a><br></div></div>
<mailto:<a href="mailto:Wannier@quantum-espresso.org" target="_blank">Wannier@quantum-espresso.org</a>><div class="im"><br>
<a href="http://www.democritos.it/mailman/listinfo/wannier" target="_blank">http://www.democritos.it/mailman/listinfo/wannier</a><br>
<br>
<br>
<br>
<br>
<br>
-- ---------------------------------------------------------------------<br>
Prof Nicola Marzari Department of Materials Science and Engineering<br>
13-5066 MIT 77 Massachusetts Avenue Cambridge MA 02139-4307 USA<br></div>
tel 617.4522758 fax 2586534 <a href="mailto:marzari@mit.edu" target="_blank">marzari@mit.edu</a> <mailto:<a href="mailto:marzari@mit.edu" target="_blank">marzari@mit.edu</a>><br>
<a href="http://quasiamore.mit.edu" target="_blank">http://quasiamore.mit.edu</a><br>
<br>
<br>
</blockquote><div><div></div><div class="h5">
<br>
<br>
-- <br>
---------------------------------------------------------------------<br>
Prof Nicola Marzari Department of Materials Science and Engineering<br>
13-5066 MIT 77 Massachusetts Avenue Cambridge MA 02139-4307 USA<br>
tel 617.4522758 fax 2586534 <a href="mailto:marzari@mit.edu" target="_blank">marzari@mit.edu</a> <a href="http://quasiamore.mit.edu" target="_blank">http://quasiamore.mit.edu</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>department of physics and astronomy<br>louisiana state univeristy<br>