<div>Dear all QE users,</div>
<div> </div>
<div>Although phonon-5.1 is available now, I have to use phonon-5.0.2 for some other codes I am using were built on phonon-5.0.2.</div>
<div> </div>
<div>Recently I found some problem that phonon-5.0.2 did not give the same phonon frequencies as phonon-5.1 for some q points of my system. The difference is big. For example,</div>
<div> </div>
<div>Phonon-5.1 gave</div>
<div>     q = (    0.000000000   0.000000000   0.053115094 )<br> **************************************************************************<br>     freq (    1) =       3.205949 [THz] =     106.938946 [cm-1]</div>
<p>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]</p>

<p>I found some codes in phonon-5.0.2 were modified by phonon-5.1 in subroutine smallg_q (xq, modenum, at, bg, nrot, s, ftau, sym, minus_q)</p>
<p> phonon-5.0.2 used<br>     if (.not.minus_q) then<br>        raq = - raq<br>        minus_q = eqvect (raq, aq, zero, accep)<br>     endif</p>
<p> 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</p>
<p>So I did this modification in phonon-5.0.2 and the modified phonon-5.0.2 gave the same phonon frequencies as those from phonon-5.1 for the same q point.</p>
<p>Since I have to continue use phonon-5.0.2 now, my question is does anyone of you know if there were some other problems found in phonon-5.0.2? Thanks.</p>
<div dir="ltr">
<div>
<div>
<div>Sincerely,</div>
<div>Yangzheng Lin<br>Postdoctoral Associate<br>Geophysical Laboratory, Carnegie Institution for Science<br>5251 Broad Branch Rd. NW, Washington, DC 20015, USA</div></div></div></div>