<div dir="ltr"><div>I do not think the routine you mention should have any effect whatsover, unless you introduce a supercell (as in old-fashioned "mass approximation").<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Nov 15, 2016 at 1:49 PM, 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 Paolo Giannozzi and community,<br>
<br>
This 4*4*4 is important because it modifies the "total_weight" variable.<br>
<br>
Example:<br>
I did an interpolation  of a 5^3 q-mesh to calculate the phonon DOS.<br>
I get the error message:<br>
<br>
frc_blk, wrong total_weight<br>
<br>
because the condition<br>
ABS(total_weight-nr1*nr2*nr3).<wbr>GT.1.0d-8)   of matdyn.f90/subroutine frc_blk<br>
is not fulfilled.<br>
in frc_blk, nr1, nr2 and nr3 have been changed through read_ifc_param and<br>
are each 5 as asked in <a href="http://ph.in" rel="noreferrer" target="_blank">ph.in</a>.<br>
<br>
But the total_weight was 124.1 and not 125.<br>
This is due to the fact that total_weight is calculated through variables<br>
that are in the subroutine q_gen which uses fixed nr:<br>
INTEGER, PARAMETER:: nr1=4, nr2=4, nr3=4<br>
<br>
When I change in q_gen<br>
INTEGER, PARAMETER:: nr1=5, nr2=5, nr3=5<br>
<br>
the total_weight changes and is exactly 125<br>
and the condition was fulfilled so that the phonon DOS run without any error.<br>
<br>
What hypothesis did you do to keep these variables fixed?<br>
Am I doing something wrong?<br>
<br>
<br>
Antoine JAY<br>
<br>
<br>
<br>
<br>
<br>
<br>
I think 4x4x4 is not directly related to the q-point grid, but to the<br>
maximum size of the supercell (if you want to build a supercell with force<br>
constants computed in the unit cell. It was used once upon a time in<br>
GaAs/AlAs heterostructure, I think)<br>
<br>
<br>
Paolo<br>
<br>
On Tue, Nov 15, 2016 at 9:00 AM, Antoine JAY <<a href="mailto:antoine....@polytechnique.edu">antoine....@polytechnique.edu</a><wbr>><br>
wrote:<br>
<br>
> 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>