<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Apr 14, 2016 at 3:49 PM, yelena <span dir="ltr"><<a href="mailto:yelena@ipb.ac.rs" target="_blank">yelena@ipb.ac.rs</a>></span> wrote:<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I keep having negative phonon frequency in gamma point, -0.00007, nowhere else (...) but it keeps interfering with my lambda.x calculation </blockquote><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">and all i get is NAN for alpha2F and so on. So is there any trick to avoid it?</blockquote><div><br></div><div>modify the lambda.f90 file and other tools so that they discard small or negative values for frequency; or, manually set them to a small positive number<br>  <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"> This is clearly numerical issue, because other modes are positive and phonon dispersion looks as it should be and there is no any instability.</blockquote><div><br></div><div>it is a known problem (violation of the Acoustic Sum Rule in approximate calculations) for which no solution is known, other than imposing the Acoustc Sum Rule at the endo of the calculation<br><br></div><div>Paolo<br clear="all"></div></div><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div></div>