Dear all,<br><br>I've read about sources in espresso/PH dir. Please correct me if my understanding about following quantities is incorrect.<br><br>1. elph_mat(nbnd, nbnd, ik, ipert): It is  <\psi_ikq | dVscf(q) / dR_ipert | \psi_ik> while ipert is along one of three Cartesian direction for a particular atom.
<br>2. dyn(iat, jat): The dynamical matrix which obeys -- dyn^{\diag} * amass * dyn = 1<br><br>Now, If I want get electron phonon matrix gkq( ibnd, jbnd, ik, imode) out, should I just do this?<br>---gkq(ibnd, jbnd, ik, imode) = \sum_{ipert} elph_mat(ibnd, jbnd, ik, ipert) * dyn(ipert, imode)
<br>The sqrt(mass) is hiden in dyn, so if I let gkq = gkq / sqrt( sqrt(w2(imode) ) ), will I get the final elph interaction matrix?<br><br>Regards,<br>Jin Zhang<br>-- <br>Dep. of Physics, Peking University<br>Beijing, P.R.China
<br>Phone: 86-10-62-753-190