<div dir="ltr"><div><div><div><div><div>Dear all,<br><br></div>I have a question about the electron phonon coupling calculation in the subroutine "elphel" (in the file elphonon.f90).<br><br></div><div>This subroutine computes <n,k|dV|m,k+q>. |n,k> is the wavefunction for band "m" with wavevector "k".<br><br></div><div>First, the program computes <n,k|dV.  Then, it is multiplied by |m,k+q><br></div><div><br></div><div>If q=0, it only needs to read to evc as:<br>...<br>        IF (lgamma) THEN<br>           CALL get_buffer(evc, lrwfc, iuwfc, ikk)<br>        ELSE<br>           CALL get_buffer (evc, lrwfc, iuwfc, ikk)<br>           CALL get_buffer (evq, lrwfc, iuwfc, ikq)<br>        ENDIF<br>...<br></div>My question is in the q=0 case, array "evq" is used to multiply the computed <n,k|dV instead of evc. (line 401-406 in version 5.2.1)<br><br></div>Can anyone explain  why "evc" array is not used in q=0 case for these line (401-406)? <br><br></div>Thanks,<br></div>Fan<br></div>