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">
<div><div>Meanwhile, I will try different 'K_POINTS' setup in '<a href="http://scf.in" target="_blank">scf.in</a>' that have been set as 'automatic 4 4 4 0 0 0'</div>
<div>I thought '4 4 4 0 0 0' would be okay for most of structures even though it will take longer time due to not pre-defined symmetry points.</div><div>I will give it a shot with '6 6 6 0 0 0'. Please let me know if you have any better suggestion.</div>

<div>For the structure like rhombohedron, should it be some thing like '4 4 6 0 0 0' ?</div><div></div></div></blockquote></div><div><br></div><div>Dear  Kyeong-hyun Park,</div><div><br></div><div>1. K points have nothing to do with your problem. Phonon q-points and electronic k-points have nothing to do with each other except for having the same unit of measure (unless you are doing electron-phonon interactions)</div>
<div><br></div><div>2. phonon frequency at the Gamma point are computed by matdyn with a formula that is 100% direction independent:</div><div>    D(q=0) =  \sum_R F_{0,R} </div><div>    even if you put random data in the force-constants file the result will NOT depend on the direction.</div>
<div><br></div><div>3. the only direction dependent term computed by matdyn is the LO-TO splitting, you think you have disable it in the code, but maybe you have not (the code is unnecessarily complicated)</div><div><br></div>
<div>4. I do not have a handy example to check out how to enable/disable LO-TO splitting, but if you want to send me the *dyn files I'll have a look</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>