<div dir="ltr">Many thanks. Yangzheng<br><div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Sep 4, 2014 at 12:43 PM, Paolo Giannozzi <span dir="ltr"><<a href="mailto:paolo.giannozzi@uniud.it" target="_blank">paolo.giannozzi@uniud.it</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">>From the Doc/release-notes file:<br>
===<br>
Fixed in 5.0.3 version:<br>
[...]<br>
   * PHonon: restart with convt=.true. was not working with PAW<br>
   * PHonon: in matdyn.x, displacement patterns were no longer correct<br>
     if eigenvectors of the dynamical matrix were also written to file<br>
   * PHonon: compatibility with files produced by previous releases<br>
     was broken if pools were used<br>
   * VERY NASTY PHonon bug: symmetry Sq=-q+G introduced in v.5.0.2<br>
     could lead to bad phonons in some cases and has been disabled<br>
   * PHonon: effective charges eu + PAW were not accurate with pools;<br>
     effective charges eu + US PP + spin-orbit nonmagnetic were not<br>
     working in parallel.<br>
<br>
I guess the "VERY NASTY PHonon bug" is the one you encountered.<br>
V.5.0.3 is very similar to v.5.0.2 so it shouldn't be difficult<br>
to adapt your codes built upon v.5.0.2<br>
<br>
P.<br>
<div><div class="h5"><br>
On Thu, 2014-09-04 at 11:57 -0400, Lin, Yangzheng wrote:<br>
> Dear all QE users,<br>
><br>
> Although phonon-5.1 is available now, I have to use phonon-5.0.2 for<br>
> some other codes I am using were built on phonon-5.0.2.<br>
><br>
> Recently I found some problem that phonon-5.0.2 did not give the same<br>
> phonon frequencies as phonon-5.1 for some q points of my system. The<br>
> difference is big. For example,<br>
><br>
> Phonon-5.1 gave<br>
>      q = (    0.000000000   0.000000000   0.053115094 )<br>
>  **************************************************************************<br>
>      freq (    1) =       3.205949 [THz] =     106.938946 [cm-1]<br>
> but Phonon-5.0.2 gave<br>
>      q = (    0.000000000   0.000000000   0.053115094 )<br>
>  **************************************************************************<br>
>      omega( 1) =      -7.856216 [THz] =    -262.055171 [cm-1]<br>
><br>
> I found some codes in phonon-5.0.2 were modified by phonon-5.1 in<br>
> subroutine smallg_q (xq, modenum, at, bg, nrot, s, ftau, sym, minus_q)<br>
><br>
>  phonon-5.0.2 used<br>
>      if (.not.minus_q) then<br>
>         raq = - raq<br>
>         minus_q = eqvect (raq, aq, zero, accep)<br>
>      endif<br>
><br>
>  phonon-5.1 modified it as<br>
>      if (sym(irot).and..not.minus_q) then<br>
>         raq = - raq<br>
>         minus_q = eqvect (raq, aq, zero, accep)<br>
>      endif<br>
><br>
> So I did this modification in phonon-5.0.2 and the modified<br>
> phonon-5.0.2 gave the same phonon frequencies as those from phonon-5.1<br>
> for the same q point.<br>
><br>
> Since I have to continue use phonon-5.0.2 now, my question is does<br>
> anyone of you know if there were some other problems found in<br>
> phonon-5.0.2? Thanks.<br>
><br>
> Sincerely,<br>
> Yangzheng Lin<br>
> Postdoctoral Associate<br>
> Geophysical Laboratory, Carnegie Institution for Science<br>
> 5251 Broad Branch Rd. NW, Washington, DC 20015, USA<br>
</div></div>> _______________________________________________<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>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
 Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
 Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
 Phone <a href="tel:%2B39-0432-558216" value="+390432558216">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222">+39-0432-558222</a><br>
<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>
</font></span></blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr"><div><div><div>Yangzheng Lin<br></div>Postdoctoral Associate<br></div>Geophysical Laboratory, Carnegie Institution for Science<br></div>5251 Broad Branch Rd. NW, Washington, DC 20015, USA<br>
</div>
</div></div></div>