<div dir="ltr"><span class="">On Fri, Nov 18, 2016 at 8:44 AM, Antoine JAY <span dir="ltr"><<a href="mailto:antoine.jay@polytechnique.edu" target="_blank">antoine.jay@polytechnique.edu</a><wbr>></span> wrote:</span><br><span class=""></span><div class="gmail_quote"><div dir="ltr"><span class=""></span><br><div class="gmail_extra"><div class="gmail_quote"><span class=""><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote">I had trouble performing phonon DOS on a nq=5^3 grid:<br>
frc_blk','wrong total_weight'<br></blockquote><div><br></div></span><div>I confirm that this is true ...<br><br></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote"><span class="">
But the total_weight seems to be calculated through parameters defined in<br></span>
the q_gen subroutine (matdyn.f90). In this subroutine, nr1 n2 and nr3 are fixed:<span class=""><br>
INTEGER, PARAMETER:: nr1=4, nr2=4, nr3=4<br>
which give in my example, a tota_weight that is 124.1, far from 125 (5^3)<br>
<br>
When I changed in q_gen the 4^3 to 5^3 values, I obtained the good weight<br>
of 125 and was able to run properly the phonon DOS.<br></span></blockquote><div><br></div><div>... and that this is not true. I made the same change and the error is still there.<br> <br></div><span class=""><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote">
I think the nr1, nr2 and nr3 that are use in q_gen must be the same as the<br>
one ask by the user and not fixed parameters.<br></blockquote><div><br></div></span><div>I don't think so. There is no relation between the nr1,nr2,nr3 parameters that define the wavevector grid and those in the q_gen subroutine that define the largest supercell for the "mass approximation" (that you did not use anyway).<br><br></div><div>The problem of occasional "wrong total_weight" errors is well known. Like for several other known phonon problems, no solution is currently available.<span class="HOEnZb"><font color="#888888"><br><br></font></span></div><span class="HOEnZb"><font color="#888888"><div>Paolo<br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="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>
<br><br></div></font></span></div></div></div></div><font color="#888888">
</font></div>