On 26 October 2012 07:55, Bo Qiu <span dir="ltr"><<a href="mailto:200210qb@gmail.com" target="_blank">200210qb@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Dr. Paulatto,<div><br></div><div>Thanks a lot for the explanations! I'm actually writing a small code to order phonon bands based on eigenvector continuities on a 3d grid. Thanks for pointing out there is already such option in matdyn.x, I'll definitely check it out. I have finished a preliminary test for my code on Al for a band crossing along Gamma-K direction -- it seems to be well captured, but band orderings are wrong at some random point where there is actually no band crossing. I'll try to figure it out and also referring to the matdyn.f90 to see what else I can do. </div>

<div>On the other hand, I think for electronic band structures one will have the same problem with band ordering -- though most of time people are only interested about those very few bands near Fermi level. I wonder in the case of electronic band structure whether continuity of wavefunction should be used to deal with band ordering with band crossing? Thanks a lot!</div>
<span class="HOEnZb"><font color="#888888">
<div><br></div></font></span><div><span class="HOEnZb"><font color="#888888"></font></span></div></blockquote></div><div><br></div><div>The bands.x code already does that, i.e. it sorts the bands using the overlap of wavefunctions at consequent k-points. Again, this is an option that works well for line plots, but I'm not sure it always make sense to use in a 3D dispersion. Or at least, it would need to be modified to sort according to neighboring points, not just according to the previous point in the input list. It is not very hard to do, but a bit boring as it requires a rethinking of the code structures.</div>
<div><br></div><div>bests</div><div><br></div>-- <br><div>Dr. Lorenzo Paulatto </div><div>IdR @ IMPMC -- CNRS & Université Paris 6</div><div>phone: +33 (0)1 44275 084 / skype: paulatz</div><div>www:   <a href="http://www-int.impmc.upmc.fr/~paulatto/" target="_blank">http://www-int.impmc.upmc.fr/~paulatto/</a></div>
<div>mail:  23-24/4é16 Boîte courrier 115, 4 place Jussieu 75252 Paris Cédex 05</div><br>
<div><br></div><div><br></div>