modes of phonons in EPW

Post here questions linked with issue while running the EPW code

Moderator: stiwari

Post Reply
gururaj

modes of phonons in EPW

Post by gururaj »

dear sir ,
i have run EPW calculation by considering a2f=.true , and i have got phonon LINEWIDTH , a2f(elaishberg) etc ..as output files .
i am trying to calculate the resistivity from elaishberg function . the question is how can i get resistivity offered by different types of modes such as (LO,TO,LA,TA,ZA,ZO)..??
In a2f output log file , i am having 12 columns and what all the columns represent??

thank you
sponce
Site Admin
Posts: 616
Joined: Wed Jan 13, 2016 7:25 pm
Affiliation: EPFL

Re: modes of phonons in EPW

Post by sponce »

Dear gururaj,

All the column are different q smearing factors (written at the bottom).

The mode decomposition can be obtained with "iverbosity 3" in the input file.

Once you have the eliashberg spectral function, you have to make the integral yourself to get the resistivity (using matlab or anything you like).

Best,

Samuel
Prof. Samuel Poncé
Chercheur qualifié F.R.S.-FNRS / Professeur UCLouvain
Institute of Condensed Matter and Nanosciences
UCLouvain, Belgium
Web: https://www.samuelponce.com
gururaj

Re: modes of phonons in EPW

Post by gururaj »

dear sir ,
as u said , for resistivity calculation , i am using the resistivity file given by EPW for how to compute resistivity in FCC LEAD example (http://epw.org.uk/Documentation/FCCLead).
but the paper i was referring to has resistivity calculated for all modes(LA,LO,TA,TO,ZA,ZO) and summed all to get total resistivity .
my question is how to get A2F elaishberg function values for all modes ??
i am getting phonon linewidth and lamda values for all modes.
please help me in resolving this.
thank you
sponce
Site Admin
Posts: 616
Joined: Wed Jan 13, 2016 7:25 pm
Affiliation: EPFL

Re: modes of phonons in EPW

Post by sponce »

In the case of a2F, this is not implemented.

However, it is implemented for the electronic self-energy when you use 'iverbosity 3'.

You can take a look at selfen_elec.f90 to see how it is done and then do something similar to a2f.90

Best,

Samuel
Prof. Samuel Poncé
Chercheur qualifié F.R.S.-FNRS / Professeur UCLouvain
Institute of Condensed Matter and Nanosciences
UCLouvain, Belgium
Web: https://www.samuelponce.com
Post Reply