<div dir="rtl"><div dir="ltr">Dear Giovanni,<br><br></div><div dir="ltr">I got it, Thank you for the clarification and have a nice week-end.<br><br><span class="gmail-HOEnZb gmail-adL"><font color="#888888"><pre style="margin:0em">Youssef Aharbil,
Laboratory of Physics and Chemistry of Materials
Faculty of sciences Ben msik, Casablanca
Morocco</pre></font></span><br></div><div class="gmail_extra"><br><div class="gmail_quote"><div dir="ltr">2016-12-03 22:26 GMT+00:00 Giovanni Cantele <span dir="ltr"><<a href="mailto:giovanni.cantele@fisica.unina.it" target="_blank">giovanni.cantele@fisica.unina.it</a>></span>:</div><blockquote class="gmail_quote" style="margin:0 .8ex;border-left:1px #ccc solid;border-right:1px #ccc solid;padding-left:1ex;padding-right:1ex"><div style="word-wrap:break-word">they are the same but printed in different units, Wannier90 in Ang^-1, QE in 2*pi/a units.<div><br></div><div>If you multiply the QE reciprocal vectors by 2*pi/a = 2*pi/alat = 2 * pi / (10.577485*0.529177) Ang^-1 = 1.1225262 Ang^-1 you exactly get the Wannier90 ones.</div><div><br></div><div>Giovanni</div><div><br><div><br><div><blockquote type="cite"><div><div class="h5"><div>On 3 Dec 2016, at 21:54, Youssef Aharbil <<a href="mailto:aharbil@gmail.com" target="_blank">aharbil@gmail.com</a>> wrote:</div><br class="m_-6837141083114159771Apple-interchange-newline"></div></div><div><div><div class="h5"><div dir="rtl"><div dir="ltr">Dear All,<br></div><div dir="ltr"><br></div><div dir="ltr">I have noticed that the reciprocal vectors computed via wannnier90 and quantum espresso (QE) are different***, which means that the Kpath and segment length across the BZ aren't the some for both code though handling the same crystal structure.<br><br><br><br></div><div dir="ltr">In order to investigate this issue, I have dug inside the source code of wannier90 and QE and compared the formula used in each code and finally realized the root cause, In fact QE divide the reciprocal vectors by an additional term, the real space vector denominator (den variable in recips.f90), which isn't the case for wannier90 (utility.F90).<br></div><div dir="ltr"><br>So, I am really confused, I can manually fix this, but I wanted to share with you.<br></div><div dir="ltr">Any suggestions or clues are most welcome.<br><br><br></div><div dir="ltr">***<br></div><div dir="ltr">Wannier90 header:<br><br>                              <wbr>      ------<br>                              <wbr>      SYSTEM<br>                              <wbr>      ------<br> <br>                              Lattice Vectors (Ang)<br>                    a_1     2.784694  -1.607742   4.581610<br>                    a_2     0.000000   3.215484   4.581610<br>                    a_3    -2.784694  -1.607742   4.581610<br> <br>                   Unit Cell Volume:     123.07314  (Ang^3)<br> <br>                        Reciprocal-Space Vectors (Ang^-1)<br>                    b_1     1.128164  -0.651347   0.457131<br>                    b_2     0.000000   1.302693   0.457131<br>                    b_3    -1.128164  -0.651347   0.457131<br><br><br></div><div dir="ltr">QE Header:<br><br>     celldm(1)=  10.577485  celldm(2)=   0.000000  celldm(3)=   0.000000<br>     celldm(4)=   0.504986  celldm(5)=   0.000000  celldm(6)=   0.000000<br><br>     crystal axes: (cart. coord. in units of alat)<br>               a(1) = (   0.497501  -0.287232   0.818530 )  <br>               a(2) = (   0.000000   0.574464   0.818530 )  <br>               a(3) = (  -0.497501  -0.287232   0.818530 )  <br><br>     reciprocal axes: (cart. coord. in units 2 pi/alat)<br>               b(1) = (  1.005024 -0.580251  0.407234 )  <br>               b(2) = (  0.000000  1.160501  0.407234 )  <br>               b(3) = ( -1.005024 -0.580251  0.407234 )  <br><br></div><div dir="ltr"><br><pre style="margin:0em">Youssef Aharbil,
Laboratory of Physics and Chemistry of Materials
Faculty of sciences Ben msik, Casablanca
Morocco</pre><br></div></div></div></div>
______________________________<wbr>_________________<br>Wannier mailing list<br><a href="mailto:Wannier@quantum-espresso.org" target="_blank">Wannier@quantum-espresso.org</a><br><a href="http://mailman.qe-forge.org/cgi-bin/mailman/listinfo/wannier" target="_blank">http://mailman.qe-forge.org/<wbr>cgi-bin/mailman/listinfo/<wbr>wannier</a><br></div></blockquote></div><br></div></div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="rtl">والله ولي التوفيق<br>يوسف أحربيل<br></div></div>
</div></div>