<div dir="ltr"><div><div><div><div>Dear all, <br><br></div>I'm a bit confused by the ordering of the q points that belong to the star of a q-point.<br><br></div>I've made a small test (diamond) on a 6x6x6 q-mesh. <br></div>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:<br>-0.166666666666667       0.166666666666667      -0.166666666666667   <br><br></div>The other 7 q-point that can be obtained by sym. operation are the following (in that order): <br><div>  0.166666666666667       0.166666666666667       0.166666666666667     <br> -0.166666666666667      -0.166666666666667      -0.166666666666667     <br>  0.166666666666667      -0.166666666666667       0.166666666666667     <br>  0.166666666666667      -0.166666666666667      -0.166666666666667     <br>  0.166666666666667       0.166666666666667      -0.166666666666667     <br> -0.166666666666667      -0.166666666666667       0.166666666666667     <br> -0.166666666666667       0.166666666666667       0.166666666666667<br><br></div><div>Those are also reported in the NAME.dyn_q02 file.<br><br></div><div>My problem is that when I use the PHonon/PH/q2qstar.x program on the same dynamical file it gives me <br></div><div>     Number of q in the star =    4<br>     List of q in the star:<br>          1  -0.166666667   0.166666667  -0.166666667<br>          2   0.166666667   0.166666667   0.166666667<br>          3   0.166666667  -0.166666667  -0.166666667<br>          4  -0.166666667  -0.166666667   0.166666667<br>     In addition there is the -q list: <br>          1   0.166666667  -0.166666667   0.166666667<br>          2  -0.166666667  -0.166666667  -0.166666667<br>          3  -0.166666667   0.166666667   0.166666667<br>          4   0.166666667   0.166666667  -0.166666667<br><br></div><div>Which is NOT the same ordering. <br><br></div><div>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<br>  sxq -0.166666667000000       0.166666667000000      -0.166666667000000     <br>  0.166666667000000       0.166666667000000       0.166666667000000     <br>  0.166666667000000      -0.166666667000000      -0.166666667000000     <br> -0.166666667000000      -0.166666667000000       0.166666667000000     <br>  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000<br>  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000<br>....<br><br></div><div>How can I obtain, in the q2qstar.x program, the same ordering as in <br>PHonon/PH/dynmatrix.f90 line 146 just at the exist of<br>  call star_q (xq, at, bg, nsym, s, invs, nq, sxq, isq, imq, .TRUE. )<br></div><div>write(*,*) 'sxq',sxq <br>?<br><br></div><div>Thank you a lot in advance, <br><br></div><div>Samuel Ponce, <br>Department of Materials, University of Oxford<br><br><br><br><br></div></div>