On 16 October 2012 05:09, 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<div><br></div><div>Thank you for your suggestion though, unfortunately commenting out the statement didn't resolve </div><div>my problem. </div></blockquote><div><br></div><div>Well, I think your first problem was solved. However I'm not sure if you  understand why there was an anisotropy in your frequencies, and the fact that the anisotropy was good, i.e. a real effect of your system.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>I thought I was pretty close, but now it seems that it's one of the input variables.</div>
<div><br></div></blockquote><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div></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></div></blockquote></div><div><br></div><div>I see that you use 5 inequivalent atomic type to model just 2 types of atoms, this means you are losing some symmetry. Symmetry helps a lot with phonons, both in terms of stability and computing time. Do you have any specific reason to do so?</div>
<div><br></div><div>I would say you probably have to increase the number of kpoint, and/or the wavefunction and/or charge density cutoff. When you are studying the phonon dispersion it is often necessary to test convergence of the phonon dispersion, i.e. compute at different values of cutoff and verify that it does not change much.</div>
<div><br></div><div>If you have someone who is more experienced with DFT calculations in your institution I suggest you contact him/her as it is difficult to understand exactly what's your problem via email. </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>