It is currently Fri Dec 13, 2019 8:34 am

Tell a friend!


Post new topic Reply to topic  [ 2 posts ] 
How to output the el-ph strength (g2) of coarse grid? 
Author Message

Joined: Wed Apr 24, 2019 7:21 am
Posts: 2
University: TU Wien
Post How to output the el-ph strength (g2) of coarse grid?
Dear all,

How to output the el-ph coupling strength (g2) of coarse grids without wannier interpolation in EPW?

Referring to what has been done for “epmatf” or “epf17”, I have tried “g2_coarse=abs(sum(uf*epmatq))**2*inv_wq”, where “uf” is phonon eigenvector and “empatq” is coarse el-ph coupling matrix as said in EPW. However, the calculated g2_coarse is different from wannier interpolated values. Then what should I do to output correct values?


Wed Apr 24, 2019 8:27 am
Profile E-mail
Site Admin
User avatar

Joined: Wed Jan 13, 2016 7:25 pm
Posts: 575
University: Oxford
Post Re: How to output the el-ph strength (g2) of coarse grid?
Hello,

You can try to uncomment
!write(*,*)'epmatq(:,:,2,:,nqc)**2',SUM((REAL(REAL(epmatq(:,:,2,:,nqc))))**2)+&
! SUM((REAL(AIMAG(epmatq(:,:,2,:,nqc))))**2)

inside elphon_shuffle_wrap.f90

Then also print the same elements from the variable "epmatf" in ephwann_shuffle.f90

Make sure you print g and its c.c. to avoid phase issues. You also need to average on degenerate states (average all the bands, modes and k-points). Only the average on degenerate state is expected to be the same.

Best,
Samuel

_________________
Dr. Samuel Poncé
Department of Materials
University of Oxford
Parks Road
Oxford OX1 3PH, UK


Tue May 14, 2019 7:37 am
Profile E-mail
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 2 posts ] 


Who is online

Users browsing this forum: No registered users and 1 guest


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum

Search for:
Jump to:  
cron
Powered by phpBB © phpBB Group.
Designed by Vjacheslav Trushkin

(All content on this board is governed by and is the sole responsibility of the board administrator.)


Gratis forum Free forum hosting| gratis phpbb3 forum | phpbb3 styles