<div dir="ltr">Hi Pietro,<br><br>Thank you for the reply.<div><br></div><div>This system actually has an issue with the symmetry while running an EPW calculation (which has been fixed temporarily for now in the code). I had the phonon done with splitting a job in each q-point generated by pw.x, and later when I ran an EPW, it crashed with error in the specific q-point which is not in the same order with ph.x. For many of the systems I checked, the way I run works. </div><div><br></div><div>I am not really worried, but trying to understand what's going on with this specific case.</div><div><br></div><div>Sincerely,</div><div>Hari</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 6, 2020 at 12:06 PM Pietro Delugas <<a href="mailto:pdelugas@sissa.it">pdelugas@sissa.it</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div lang="EN-US"><div class="gmail-m_4402175496600331969WordSection1"><p class="MsoNormal"><span lang="IT">Hi Hari <u></u><u></u></span></p><p class="MsoNormal"><span lang="IT"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="IT">They are generated by two different routines and probably with two different algorithms.  <u></u><u></u></span></p><p class="MsoNormal"><span lang="IT">Why are you worried about their ordering ? <u></u><u></u></span></p><p class="MsoNormal"><span lang="IT"><u></u> <u></u></span></p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt">they are two different things. <u></u><u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt"><u></u> <u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt">One is a MP mesh used to estimate integrals of periodic functions in the Brillouin Zone exploiting phase cancelations;   they are used all together to make an SCF calculation.<u></u><u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt"> <u></u><u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt"> The other one is a uniform grid of q points at each one of which to compute independently a dynamical matrix using all k points (plus others) of the previous list.  Then taken  all together the dynamical matrices  form the Fourier transform of the interatomic force constants in a 6X6X4 supercell.  <u></u><u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt"><u></u> <u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt">Regards - Pietro <u></u><u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt"><u></u> <u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt"><u></u> <u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt"><u></u> <u></u></span></span></p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt"><u></u> <u></u></span></span></p><p class="MsoNormal">Sent from <a href="https://go.microsoft.com/fwlink/?LinkId=550986" target="_blank">Mail</a> for Windows 10</p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt"><u></u> <u></u></span></span></p><div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(225,225,225);padding:3pt 0in 0in"><p class="MsoNormal" style="border:none;padding:0in"><b>From: </b><a href="mailto:hpaudya1@binghamton.edu" target="_blank">Hari Paudyal</a><br><b>Sent: </b>Thursday, August 6, 2020 4:00 PM<br><b>To: </b><a href="mailto:users@lists.quantum-espresso.org" target="_blank">Quantum Espresso users Forum</a><br><b>Subject: </b>[QE-users] difference in the order of k/q points with pw.x and ph.x</p></div><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt"><u></u> <u></u></span></span></p><div><p class="MsoNormal">Hi all,</p><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">I am doing a phonon calculation of a trigonal system with SG 164. I noticed that the order of the k (q) points printed from pw.x  and ph.x is different for the same k (q) grids.</p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">my inputs/outputs are as the following:</p></div><div><p class="MsoNormal"><b>pw.x input  Â </b> Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â <b>ph.x input</b></p></div><div><p class="MsoNormal">K_POINTS {automatic}  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ldisp = .true.<br>6 6 4 0 0 0  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  nq1 = 6, nq2 = 6, nq3 = 4</p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><b>pw.x output</b></p></div><div><p class="MsoNormal">  Â  Â number of k points= Â  Â 24 Â Marzari-Vanderbilt smearing, width (Ry)= Â 0.0200<br>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â cart. coord. in units 2pi/alat<br>  Â  Â  Â  k( Â  Â 1) = ( Â  0.0000000 Â  0.0000000 Â  0.0000000), wk = Â  0.0138889<br>  Â  Â  Â  k( Â  Â 2) = ( Â  0.0000000 Â  0.0000000 Â  0.1707981), wk = Â  0.0277778<br>  Â  Â  Â  k( Â  Â 3) = ( Â  0.0000000 Â  0.0000000 Â -0.3415962), wk = Â  0.0138889<br>  Â  Â  Â  k( Â  Â 4) = ( Â  0.0000000 Â  0.1924501 Â  0.0000000), wk = Â  0.0833333<br>  Â  Â  Â  k( Â  Â 5) = ( Â  0.0000000 Â  0.1924501 Â  0.1707981), wk = Â  0.0833333<br>  Â  Â  Â  k( Â  Â 6) = ( Â  0.0000000 Â  0.1924501 Â -0.3415962), wk = Â  0.0833333<br>  Â  Â  Â  k( Â  Â 7) = ( Â  0.0000000 Â  0.3849002 Â  0.0000000), wk = Â  0.0833333<br>  Â  Â  Â  k( Â  Â 8) = ( Â  0.0000000 Â  0.3849002 Â  0.1707981), wk = Â  0.0833333<br>  Â  Â  Â  k( Â  Â 9) = ( Â  0.0000000 Â  0.3849002 Â -0.3415962), wk = Â  0.0833333<br>  Â  Â  Â  k( Â  10) = ( Â  0.0000000 Â -0.5773503 Â  0.0000000), wk = Â  0.0416667<br>  Â  Â  Â  k( Â  11) = ( Â  0.0000000 Â -0.5773503 Â  0.1707981), wk = Â  0.0833333<br>  Â  Â  Â  k( Â  12) = ( Â  0.0000000 Â -0.5773503 Â -0.3415962), wk = Â  0.0416667<br>  Â  Â  Â  k( Â  13) = ( Â  0.1666667 Â  0.2886751 Â  0.0000000), wk = Â  0.0833333<br>  Â  Â  Â  k( Â  14) = ( Â  0.1666667 Â  0.2886751 Â  0.1707981), wk = Â  0.1666667<br>  Â  Â  Â  k( Â  15) = ( Â  0.1666667 Â  0.2886751 Â -0.3415962), wk = Â  0.0833333<br>  Â  Â  Â  k( Â  16) = ( Â  0.1666667 Â  0.4811252 Â  0.0000000), wk = Â  0.1666667<br>  Â  Â  Â  k( Â  17) = ( Â  0.1666667 Â  0.4811252 Â  0.1707981), wk = Â  0.1666667<br>  Â  Â  Â  k( Â  18) = ( Â  0.1666667 Â  0.4811252 Â -0.3415962), wk = Â  0.1666667<br>  Â  Â  Â  k( Â  19) = ( Â  0.3333333 Â  0.5773503 Â  0.0000000), wk = Â  0.0277778<br>  Â  Â  Â  k( Â  20) = ( Â  0.3333333 Â  0.5773503 Â  0.1707981), wk = Â  0.0555556<br>  Â  Â  Â  k( Â  21) = ( Â  0.3333333 Â  0.5773503 Â -0.3415962), wk = Â  0.0277778<br>  Â  Â  Â  k( Â  22) = ( Â  0.0000000 Â  0.1924501 Â -0.1707981), wk = Â  0.0833333<br>  Â  Â  Â  k( Â  23) = ( Â  0.0000000 Â  0.3849002 Â -0.1707981), wk = Â  0.0833333<br>  Â  Â  Â  k( Â  24) = ( Â -0.1666667 Â  0.4811252 Â -0.1707981), wk = Â  0.1666667</p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><b>ph.x output</b></p></div><div><p class="MsoNormal">  Â  Â Dynamical matrices for ( 6, 6, 4) Â uniform grid of q-points<br>  Â  Â ( Â 24 q-points):<br>  Â  Â  Â N Â  Â  Â  Â  xq(1) Â  Â  Â  Â  xq(2) Â  Â  Â  Â  xq(3) <br>  Â  Â  Â 1 Â  0.000000000 Â  0.000000000 Â  0.000000000<br>  Â  Â  Â 2 Â  0.000000000 Â  0.000000000 Â  0.170798095<br>  Â  Â  Â 3 Â  0.000000000 Â  0.000000000 Â -0.341596190<br>  Â  Â  Â 4 Â  0.000000000 Â  0.192450090 Â  0.000000000<br>  Â  Â  Â 5 Â  0.000000000 Â  0.192450090 Â  0.170798095<br>  Â  Â  Â 6 Â  0.000000000 Â  0.192450090 Â -0.341596190<br>  Â  Â  Â 7 Â  0.000000000 Â  0.192450090 Â -0.170798095<br>  Â  Â  Â 8 Â  0.000000000 Â  0.384900179 Â  0.000000000<br>  Â  Â  Â 9 Â  0.000000000 Â  0.384900179 Â  0.170798095<br>  Â  Â  10 Â  0.000000000 Â  0.384900179 Â -0.341596190<br>  Â  Â  11 Â  0.000000000 Â  0.384900179 Â -0.170798095<br>  Â  Â  12 Â  0.000000000 Â -0.577350269 Â  0.000000000<br>  Â  Â  13 Â  0.000000000 Â -0.577350269 Â  0.170798095<br>  Â  Â  14 Â  0.000000000 Â -0.577350269 Â -0.341596190<br>  Â  Â  15 Â  0.166666667 Â  0.288675135 Â  0.000000000<br>  Â  Â  16 Â  0.166666667 Â  0.288675135 Â  0.170798095<br>  Â  Â  17 Â  0.166666667 Â  0.288675135 Â -0.341596190<br>  Â  Â  18 Â  0.166666667 Â  0.481125224 Â  0.000000000<br>  Â  Â  19 Â  0.166666667 Â  0.481125224 Â  0.170798095<br>  Â  Â  20 Â  0.166666667 Â  0.481125224 Â -0.341596190<br>  Â  Â  21 Â  0.166666667 Â  0.481125224 Â -0.170798095<br>  Â  Â  22 Â  0.333333333 Â  0.577350269 Â  0.000000000<br>  Â  Â  23 Â  0.333333333 Â  0.577350269 Â  0.170798095<br>  Â  Â  24 Â  0.333333333 Â  0.577350269 Â -0.341596190</p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">For example, you can see the 7th point from pw.x is the 8th point from ph.x and so on.</p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Any suggestions and explanations are highly appreciated.</p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Sincerely,</p></div></div><p class="MsoNormal">Hari Paudyal</p><p class="MsoNormal"><span class="gmail-m_4402175496600331969DefaultFontHxMailStyle"><span style="font-size:12pt"><u></u> <u></u></span></span></p></div></div>_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>