How to get the el-ph coupling matrix in complex number form?

General discussion around the EPW software

Moderator: stiwari

Post Reply
edward
Posts: 1
Joined: Tue May 17, 2022 8:54 am
Affiliation: Fudan University

How to get the el-ph coupling matrix in complex number form?

Post by edward »

Dear Developers,

Along with the growth of the diversity of the electron-phonon topic greatly powered by your code, it's sometimes necessary to get the coupling matrix element in complex number form rather than |g|. I wonder if there is already a way to get it or maybe would you like to consider it to be a new feature in the next version of EPW code.

hlee
Posts: 415
Joined: Thu Aug 03, 2017 12:24 pm
Affiliation: The University of Texas at Austin

Re: How to get the el-ph coupling matrix in complex number form?

Post by hlee »

Dear edward:

For the moment, we don't have any plan to add it.
One of the reasons to print out |g| instead of g is to confirm the Wannier-interpolated electron-phonon matrix elements are comparable to the corresponding values from the direct DFPT calculations. As you know well, g is not a gauge-invariant quantity and if we print out g, it is not trivial to compare it with g from the phonon code.

I think that you can easily print out these values by modifying the subroutine of print_gkk in printing.f90.

Sincerely,

H. Lee

baiyunfei
Posts: 5
Joined: Thu Dec 22, 2022 2:12 am
Affiliation: Institute of Physics CAS

Re: How to get the el-ph coupling matrix in complex number form?

Post by baiyunfei »

As far as I know, EPW use wannier interpoltation to interpolate dvscf and wavefunction on a finer mesh ,then use them to construct electron-phonon matrix. However wannier interplotation is known to bring some gauge problems.
My question is : if I properly modify the ''print_gkk in printing.f90.'' ; set the coarse mesh of nscf calculation the same as fine mesh of wannier interpolation , and get the output of complex electron-phonon matrixs, will the el-ph matrixs' gauge be in consistent with nscf *SAVE file? if not ,what should I do to use wannier interpoltated dvscf and uninterpoltated nscf wave function to construct electron-phonon matrix?

baiyunfei
Posts: 5
Joined: Thu Dec 22, 2022 2:12 am
Affiliation: Institute of Physics CAS

Re: How to get the el-ph coupling matrix in complex number form?

Post by baiyunfei »

I am quite interested in complex value of electron-phonon matrix in order to do exciton-phonon interactions.I have read your discussion about it. As far as I know, EPW use wannier interpoltation to interpolate dvscf and wavefunction on a finer mesh ,then use them to construct electron-phonon matrix. However wannier interplotation is known to bring some gauge problems.
My question is : if I properly modify the ''print_gkk in printing.f90.'' ; set the coarse mesh of nscf calculation the same as fine mesh of wannier interpolation , and get the output of complex electron-phonon matrixs, will the el-ph matrixs' gauge be in consistent with nscf *SAVE file? if not ,what should I do to use wannier interpoltated dvscf and uninterpoltated nscf wave function to construct electron-phonon matrix?

Post Reply