On 15 October 2012 07:24, Kyeong-hyun Park <span dir="ltr"><<a href="mailto:kpark39@illinois.edu" target="_blank">kpark39@illinois.edu</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. Lorenzo Paulatto and other users<div><br></div><div>First of all, thank you for your advice.<br><div></div></div></blockquote></div><div><br></div><div><br></div><div>Hi,</div><div>in order to properly remove the anisotropy   you must also comment out the line</div>
<div>           CALL nonanal (nat, nat_blk, itau_blk, epsil, qhat, zeu, omega, dyn)</div><div>which adds some non-analitical long range term even if lo_to is set to false.</div><div><br></div><div>I also notice that your matdyn complained for a large difference between the old and new force constants when interpolation, which may signal that your phonon calculation is non very converged (especially w.r.t the number of plane waves), or not very stable, or that you are using less symmetry that you should.</div>
<br clear="all"><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>