<div dir="ltr"><div>I think 4x4x4 is not directly related to the q-point grid, but to the maximum size of the supercell (if you want to build a supercell with force constants computed in the unit cell. It was used once upon a time in GaAs/AlAs heterostructure, I think)<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Nov 15, 2016 at 9:00 AM, Antoine JAY <span dir="ltr"><<a href="mailto:antoine.jay@polytechnique.edu" target="_blank">antoine.jay@polytechnique.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear community and developers,<br>
<br>
Looking at the code of qe_6.0<br>
<br>
I am wondering why in PHonon/PH/matdyn.f90<br>
subroutine q_gen<br>
<br>
the values nr1, nr2 and nr4 are fixed to be 4:<br>
INTEGER, PARAMETER:: nr1=4, nr2=4, nr3=4<br>
<br>
whereas they should be fixed as nr1=nq1, nr2=nq2, and nr3=nq3.<br>
as done by the first call of read_ifc_param(nr1, nr2, nr3).<br>
<br>
Is there a physical meaning?<br>
Is this subroutine right if someone don't whant a 4^3 q-mesh?<br>
<br>
Antoine Jay<br>
<br>
______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br>
</blockquote></div><br><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>
</div>