<div dir="ltr"><div><div><div><div><div><div><div><div>Dear espresso users and developers,<br><br></div>I am trying to read through the electron-phonon coupling code and have a few questions. <br><br>1. From phq_readin.f90<br>  SELECT CASE( trim( electron_phonon ) )<br>  CASE( 'simple' )<br>     elph=.true.<br>     elph_mat=.false.<br>     elph_simple=.true. <br>  CASE( 'Wannier' )<br>     elph=.true.<br>     elph_mat=.true.<br>     elph_simple=.false.<br>     auxdvscf=trim(fildvscf)<br>  CASE( 'interpolated' )<br>     elph=.true.<br>     elph_mat=.false.<br>     elph_simple=.false.<br>  CASE DEFAULT<br>     elph=.false.<br>     elph_mat=.false.<br>     elph_simple=.false.<br>  END SELECT<br><br></div>What does the case electron_phonon = 'Wannier' correspond to? I see from phonon.f90 that the code would call ep_matrix_element_wannier.f90 if elph_mat = .true. but this does not seem to be explained anywhere in the documentation ( <a href="http://www.quantum-espresso.org/wp-content/uploads/Doc/INPUT_PH.html#idp74656">http://www.quantum-espresso.org/wp-content/uploads/Doc/INPUT_PH.html#idp74656</a> )? Is this similar to an external program EPW which performs Wannier interpolation from coarse to fine k-point grids. <br><br></div>2. From phonon.f90<br><br>           IF ( elph_mat ) then<br>              call ep_matrix_element_wannier()<br>           ELSE<br>              CALL elphon()<br>           END IF<br><br></div>elph_mat = .false. for both electron_phonon = 'simple' and electron_phonon = 'interpolated'. So I guess elphon.f90 is invoked in both simple and interpolated cases. Also in elphon.f90, there is no mention of the variable 'elph_simple'. It looks to me that both simple and interpolated cases would lead to the same code being run. Where exactly is there a difference between these two cases? <br><br></div>Thanks so much in advance!<br></div>Sridhar<br></div>Graduate Student<br></div>Purdue University<br><div><div><div><div><div><div><br><br></div></div></div></div></div></div></div>