<div dir="ltr">Hi thanks for the reply. I need to compute a sort of finite temperature electron-phonon coupling, so I needed to go into the code and I wrote my own routine to do that, I am debugging the code and I noticed that the result at different but equivalent q points, let's say X=(1,0,0) and X=(1,1,0) in FCC are not the same as I should expect.<div>Then I noticed that the wave function (evq) in reciprocal space differs (that is expected) but I do not understand  the relation between the two, that was why I was asking the previous question. In the meantime I found a mistake in the way I was computing the potential derivative, that is probably the main source of my problem. I was wondering also what does the routine symdyn_munu_new exactly do? Why does the dynamical matrix have to be symmetrized in the basis of the modes ?</div><div><br></div><div><div style="color:rgb(80,0,80)">Thanks in advance,<br></div><div style="color:rgb(80,0,80)">Jacopo Simoni, Lawrence Berkeley National Lab.</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, 27 Sept 2021 at 00:20, Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com">paulatz@gmail.com</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 dir="auto"><div>Hello Jacopo,</div><div dir="auto">instead of answering your question, I may ask you what you actually want to do, because chances are that your problem has already been met by others.</div><div dir="auto"><br></div><div dir="auto">E.g. there are ways to eliminate this phase, or to neutralize it to compute the derivative w.r.t. the wave vector. But more often, when you have an observable quantity that comes from a sum over the k-points, the phases will cancel out in the total, if done correctly. It is actually a good way to check that your formulas are correct.</div><div dir="auto"><br></div><div dir="auto">Hth<br><br><div dir="auto">-- <br>Lorenzo Paulatto</div><br><div class="gmail_quote" dir="auto"><div dir="ltr" class="gmail_attr">On Sat, Sep 25, 2021, 21:28 Jacopo Simoni via users <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</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 dir="ltr">So this means the phase factor is of the form e^{iG(k)r} where G(k) is the translation vector in rec. space such that<div>G(k)+q+k=q'+k' (inside the first Brillouin zone)</div><div>the phase factor is therefore dependent on the k vector and in reciprocal space the wave functions of the two equivalent q points are shifted by the vector G(k) ?</div><div>Is there a variable in the code corresponding to this vector G or to the phase factor itself ?</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, 25 Sept 2021 at 03:42, Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com" rel="noreferrer" target="_blank">p.giannozzi@gmail.com</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 dir="ltr"><div dir="ltr">On Sat, Sep 25, 2021 at 2:00 AM Jacopo Simoni via users <<a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer" target="_blank">users@lists.quantum-espresso.org</a>> wrote:</div><div dir="ltr"><br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">This wave function appears different from the wave function at an equivalent q point, for instance if I look at evq at q=(0,0,1), this is different from evq at (1,1,0) that are equivalent by translation of a G vector (I am thinking here at a FCC periodic lattice). The two functions just differ by a phase factor or I am missing something ?</div></blockquote><div><br></div><div>They differ by a phase factor; moreover, in the presence of degenerate eigenvalues, you have no guarantee that the eigenvectors in the degenerate subspace are the same. 
Finally, the ordering of k+G components is not necessarily the same in the two cases<br></div><div><br></div><div>Paolo<br></div><div><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 dir="ltr"><div></div><div>Thanks in advance, <br></div><div>Jacopo Simoni, Lawrence Berkeley National Lab.<br></div></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div></div>
</blockquote></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" rel="noreferrer" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div></div></div>
</blockquote></div>