Splitting of coarse el-ph matrix elements for polar material
Posted: Fri Sep 01, 2017 10:12 pm
Dear EPW Users,
I am attempting to add electrostatic screening (similar to the procedure in Nature Communications,8,15769,2017) for polar materials to EPW. My understanding is that according to PRL (115,176401,2015), EPW does the following:
So I would like to apply screening to both coarse and fine grid gL. However, the code does not seem to split the coarse grid g into gL and gS. It does add the gL for the fine grid in the rgd_blk_epw subroutine and ephwann_shuffle.f90. It also splits the dynamical matrix into L and S parts using the rgd_blk subroutine.
I have two questions:
1. Where does the code actually split the coarse matrix elements into the L and S components as suggested in (ii)?
2. Would it be sufficient to apply screening by dividing both the coarse and fine gL by epsilon(q), as suggested in the Nature Communications paper or are there other parts of the code that should be modified.
Thanks,
Vahid
Vahid Askarpour
Department of Physics and Atmospheric Science
Dalhousie University,
Halifax, NS, Canada
I am attempting to add electrostatic screening (similar to the procedure in Nature Communications,8,15769,2017) for polar materials to EPW. My understanding is that according to PRL (115,176401,2015), EPW does the following:
Code: Select all
(i) evaluate the complete matrix elements g on coarse Brillouin zone grids
(ii) subtract gL so as to obtain the short-ranged part of the matrix element, gS
(iii) apply Wannier-Fourier electron-phonon interpolation to the short-range matrix element, following Ref. [22]
(iv) add up the short-range part and the long-range part at arbitrary k and q points after interpolation.
So I would like to apply screening to both coarse and fine grid gL. However, the code does not seem to split the coarse grid g into gL and gS. It does add the gL for the fine grid in the rgd_blk_epw subroutine and ephwann_shuffle.f90. It also splits the dynamical matrix into L and S parts using the rgd_blk subroutine.
I have two questions:
1. Where does the code actually split the coarse matrix elements into the L and S components as suggested in (ii)?
2. Would it be sufficient to apply screening by dividing both the coarse and fine gL by epsilon(q), as suggested in the Nature Communications paper or are there other parts of the code that should be modified.
Thanks,
Vahid
Vahid Askarpour
Department of Physics and Atmospheric Science
Dalhousie University,
Halifax, NS, Canada