Dear Dr. Lorenzo Paulatto<div><br></div><div>Thank you for your suggestion though, unfortunately commenting out the statement didn't resolve </div><div>my problem. I thought I was pretty close, but now it seems that it's one of the input variables.</div>
<div><br></div><div>Based on what did you realize that my matdyn complained for non-converged or non-stable system?</div><div>Would lowering either convergence threshold (conv_thr) or threshold for selfconsistency (th2_ph) make</div>
<div>my system more stable? For now I have 'conv_thr=1.0d-8' and 'th2_ph=1.0d-14'</div><div><br></div><div>If it was a matter of symmetry, shouldn't I look at the 'K_POINTS' now?</div><div><br>
</div><div>I am sorry for asking each step of this process, but I have been stuck at this point for a while and lost.</div><div><br></div><div>Thank you.</div><div><br></div><div>Best Regards,</div><div><br></div><div><br>
</div><div><br></div><div><br><div class="gmail_quote">On Mon, Oct 15, 2012 at 5:32 AM, 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"><div class="im">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><div class="im"><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><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>
<div class="HOEnZb"><div class="h5">
<br clear="all"><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>
</div></div><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>