<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div>Dear Lorenzo, <br><br></div>The  q2qstar should give the same order but it does not. <br><br></div>Actually there is more than the order. The signs can be different. <br><br></div>For example I get <br></div>0.5 0.5 0.5 in the sxq variable of the phonon code and -0.5 -0.5 -0.5 in the sxq variable of the q2qstar code. <br></div><br></div>I know that it is the same point in diamond but this can leads to problems, especially when the crystal does not respect time reversal (maybe in that case the code produce the same result, I did not test). <br></div>Additionally, it does not really make sense that they produce different results. <br><br></div>My <a href="http://scf.in">scf.in</a> input is<br> &control<br>    calculation     = 'scf'<br>    prefix          = 'diam'<br>    restart_mode    = 'from_scratch'<br>    wf_collect      = .false.<br>    pseudo_dir      = '../../pp'<br>    outdir          = './'<br>    tprnfor         = .true.<br>    tstress         = .true.<br> /<br> &system<br>    ibrav           = 2<br>    celldm(1)       = 6.64245<br>    nat             = 2<br>    ntyp            = 1<br>    ecutwfc         = 40<br>    nbnd            = 4<br> /<br> &electrons<br>    diagonalization = 'david'<br>    mixing_beta     = 0.7<br>    conv_thr        = 1.0d-10<br> /<br>ATOMIC_SPECIES<br>  C  12.01078  C.UPF<br>ATOMIC_POSITIONS alat<br>  C   0.00  0.00  0.00<br>  C   0.25  0.25  0.25<br>K_POINTS automatic<br>4 4 4 1 1 1<br><br></div>My <a href="http://ph.in">ph.in</a> input is<br>  prefix   = 'diam',<br>  epsil    = .true.,<br>  fildyn   = 'diam.dyn',<br>  ldisp    = .true.<br>  fildvscf = 'dvscf'<br>  nq1=6,<br>  nq2=6,<br>  nq3=6,<br>  tr2_ph   =  1.0d-12,<br> /<br><br></div>Best Regards, <br><br></div>Samuel Ponce, <br>Department of Materials, University of Oxford<div><div><div><br><div><br><div><div><div><div><div><br> <br></div></div></div></div></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-06-05 21:57 GMT+01:00 Lorenzo Paulatto <span dir="ltr"><<a href="mailto:lorenzo.paulatto@impmc.upmc.fr" target="_blank">lorenzo.paulatto@impmc.upmc.fr</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><br>
<br>
On 05/06/15 22:26, Samuel Poncé wrote:<br>
> Dear all,<br>
><br>
> I'm a bit confused by the ordering of the q points that belong to the<br>
> star of a q-point.<br>
<br>
</span>Dear Samuel,<br>
the order is not important, if for some reason you need a specific order<br>
your code is going to break sooner or later.<br>
<br>
Anyway, q2qstar uses the same subroutines as phonon, so in principle it<br>
should give the same order. If you have a valid reason for requiring the<br>
same order and if you provide enough input and example files to<br>
reproduce your case, I may have a look.<br>
<br>
regards<br>
> --<br>
> Dr. Lorenzo Paulatto<br>
> IdR @ IMPMC -- CNRS & Université Paris 6<br>
> <a href="tel:%2B33%20%280%291%2044%20275%20084" value="+33144275084">+33 (0)1 44 275 084</a> / skype: paulatz<br>
> <a href="http://www-int.impmc.upmc.fr/~paulatto/" target="_blank">http://www-int.impmc.upmc.fr/~paulatto/</a><br>
> 23-24/4é16 Boîte courrier 115, 4 place Jussieu 75252 Paris Cédex 05<br>
_______________________________________________<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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br></div>