<div dir="ltr">Hi everyone,<div><br></div><div>I recently came across some phonon softening in a system calculated with the 5.0.2 version of QE. This came as a surprise because, in addition to having some physical reasons for not believing it, it was in contrast to previous calculations I had done some time back on the same system (and input files) with QE 4.3.2. The deviations with respect to QE version over a 3x3x3 grid of q vectors in this system amounted to a maximal error of 13.5 cm^-1 for a particular q vector. </div>
<div><br></div><div>This previously mentioned system was fairly large, so I tried to find simpler systems to run so that I could narrow down where the problem was coming from. I found that many simple systems (like cubic diamond) exhibited essentially zero error regardless of version. However, in guessing that the problem could possibly be related to systems with hexagonal symmetry, I chose to test the hexagonal diamond phase of Ge (4 atoms/cell) and found that over a 3x3x3 grid of q vectors, I got deviations on the order of 95 cm^-1 for some q vectors. I actually do not think the problem is limited to hexagonal symmetry systems however, as I've seen some deviations in other orthorhombic systems (although they have been more minor). </div>
<div><br></div><div style>Using the hexagonal diamond system, I tested a wide variety of QE releases and SVN versions. The problem occurs for the more recent releases of 5.0.2 and 5.0.3 (i.e., the patch doesn't fix the problem). The latest SVN version is consistent with 5.0.1 and 5.0.0 as well as older versions tested (4.3.2 and 4.2.1). I'm not sure the exact SVN version where the problem first occurred, but rev 9724 is the first commit which gives the "correct" phonon frequencies. </div>
<div style><br></div><div style>Rev 9724 differs from 9723 as follows:</div><div style><br></div><div style><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
login1$ svn -r 9723:9724 diff<br>Index: PHonon/PH/sgam_ph.f90<br>===================================================================<br>--- PHonon/PH/sgam_ph.f90<span class="" style="white-space:pre">    </span>(revision 9723)<br>
+++ PHonon/PH/sgam_ph.f90<span class="" style="white-space:pre">      </span>(revision 9724)<br>@@ -161,7 +161,8 @@<br>            sym(irot) = sym(irot) .and. (abs(raq(ipol)-aq(ipol)) < 1.0d-5)<br>         enddo<br>      endif<br>
-     if (.not.minus_q) then<br>+!     if (.not.minus_q) then<br>+     if (sym(irot).and..not.minus_q) then<br>         raq = - raq<br>         minus_q = eqvect (raq, aq, zero, accep)<br>      endif</blockquote><div><br></div>
<div style>Anyway, hope these tests are useful. Unless I've made a mistake somewhere, I think it suggests that it's extremely dangerous to do phonon calculations with 5.0.2 or 5.0.3 and you should either stick to 5.0.1 or svn revision 9724 or later. </div>
<div style><br></div><div style>Best,</div><div style>Brad Malone</div><div style>Harvard University</div><div><br></div></div><div><br></div><div><br></div></div>