[Pw_forum] Order of the q in the star for each dyn files

[Samuel Poncé]

Dear all,

I'm a bit confused by the ordering of the q points that belong to the star
of a q-point.

I've made a small test (diamond) on a 6x6x6 q-mesh.
The second dynamical file produced by the phonon code (the first one is
Gamma and poses no problem here) corresponds to the irr. q-point:
-0.166666666666667       0.166666666666667      -0.166666666666667

The other 7 q-point that can be obtained by sym. operation are the
following (in that order):
  0.166666666666667       0.166666666666667       0.166666666666667
 -0.166666666666667      -0.166666666666667      -0.166666666666667
  0.166666666666667      -0.166666666666667       0.166666666666667
  0.166666666666667      -0.166666666666667      -0.166666666666667
  0.166666666666667       0.166666666666667      -0.166666666666667
 -0.166666666666667      -0.166666666666667       0.166666666666667
 -0.166666666666667       0.166666666666667       0.166666666666667

Those are also reported in the NAME.dyn_q02 file.

My problem is that when I use the PHonon/PH/q2qstar.x program on the same
dynamical file it gives me
     Number of q in the star =    4
     List of q in the star:
          1  -0.166666667   0.166666667  -0.166666667
          2   0.166666667   0.166666667   0.166666667
          3   0.166666667  -0.166666667  -0.166666667
          4  -0.166666667  -0.166666667   0.166666667
     In addition there is the -q list:
          1   0.166666667  -0.166666667   0.166666667
          2  -0.166666667  -0.166666667  -0.166666667
          3  -0.166666667   0.166666667   0.166666667
          4   0.166666667   0.166666667  -0.166666667

Which is NOT the same ordering.

To be more accurate, I'm interested by the variable 'sxq' in the code. I
have printed the value of this variable juste before writing of the
dynamical file and it corresponds well to the first set of 8 q-points and
then I have printed it in the q2qstar.x program and get
  sxq -0.166666667000000       0.166666667000000
-0.166666667000000
  0.166666667000000       0.166666667000000       0.166666667000000
  0.166666667000000      -0.166666667000000      -0.166666667000000
 -0.166666667000000      -0.166666667000000       0.166666667000000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
....

How can I obtain, in the q2qstar.x program, the same ordering as in
PHonon/PH/dynmatrix.f90 line 146 just at the exist of
  call star_q (xq, at, bg, nsym, s, invs, nq, sxq, isq, imq, .TRUE. )
write(*,*) 'sxq',sxq
?

Thank you a lot in advance,

Samuel Ponce,
Department of Materials, University of Oxford
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20150605/68511e2a/attachment.html>