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>
<div><br></div><div>Bo<br><br><div class="gmail_quote">On Thu, Oct 25, 2012 at 5:05 PM, Lorenzo Paulatto <span dir="ltr"><<a href="mailto:lorenzo.paulatto@impmc.upmc.fr" target="_blank">lorenzo.paulatto@impmc.upmc.fr</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">Dear Bo,</span><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">

as long as you can properly characterize the band, it is not hard to do: you can use matdyn.x to compute the phonon in a dense grid of q-points, you can generate it with pw.x or with the kpoint.x utility. Then load the grid in matlab or octave or similar and integrate it properly. Since some time ago, the matdyn code produce a *.gp code which is very practical for this kind of analysis.</div>

<div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif"><br></div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
The problem is deciding which is band X, as soon as they start to cross deciding which is which is not trivial. There is an option in matdyn to assign the phonon bands order using eigenvectors continuity. However, this method works well for high-symmetry lines, but may be not so robust when doing a 3D grid of points. Yet, it is worth trying.</div>

<div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif"><br></div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
Provided all these difficulties, in many case and as long as you are interested in a phonon band that is well separated from  the others there may be no problem at all. Finally, you can think of a specific criteria to discriminate the bands in a specific material, but implementing it is up to you. </div>

<div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif"><br></div><div style="color:rgb(34,34,34);font-size:13px;font-family:arial,sans-serif">
bests</div><span class="HOEnZb"><font color="#888888"><div><br></div>-- <br><div>Dr. Lorenzo Paulatto </div><div>IdR @ IMPMC -- CNRS & Université Paris 6</div><div>phone: <a href="tel:%2B33%20%280%291%2044275%20084" value="+33144275084" target="_blank">+33 (0)1 44275 084</a> / 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>
</font></span><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div>