On 9 November 2012 01:14, Yao Yao <span dir="ltr"><<a href="mailto:yao.yao@unsw.edu.au" target="_blank">yao.yao@unsw.edu.au</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">
Looks this is similar to the problem reported here<br>
<a href="http://www.democritos.it/pipermail/pw_forum/2009-June/013250.html" target="_blank">http://www.democritos.it/pipermail/pw_forum/2009-June/013250.html</a><br>
But seems that problem has been solved according to the following email. Is<br>
this another problem in q2r or matdyn?<br><br></blockquote><div></div></div><div><br></div>As I stated several times in the discussion that you are linking, those "weird jumps" are <b>an actual physical phenomenon</b>. If you remove them <b>you are doing it wrong</b>! The problem is not with the code, but your understanding of the problem. I do not want to be rude, this is just to <i>warn</i> you that you may be going in the wrong direction.<div>
<br></div><div>Anyway:</div><div>* if you have understood that LO-TO splitting is induced by the long range interaction of effective charges in polar insulators/semiconductors</div><div>* if you still think that removing this splitting is the right thing to do</div>
<div>* if you promise to <b>not</b> cite me as the source of this information</div><div><br></div><div>then, edit PH/matdyn.f90 and around line 466 replace</div><div><div>           ELSE</div><div>              lo_to_split=.TRUE.</div>
<div>           ENDIF</div><div>           !</div><div>           CALL nonanal (nat, nat_blk, itau_blk, epsil, qhat, zeu, omega, dyn)</div></div><div>with</div><div><div>           ELSE</div><div>              lo_to_split=.FALSE.</div>
<div>           ENDIF</div><div>           !</div><div>           !!! CALL nonanal (nat, nat_blk, itau_blk, epsil, qhat, zeu, omega, dyn)</div></div><div>i.e. set lo_to_split to .false. in every case and do not compute the non-analitical long-range term.</div>
<div><br></div><div>Then recompile the phonon code</div><div><br></div><div>best regards</div><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>
</div>