Page 1 of 1

How to get k-resolved Eliashberg spectrum

Posted: Tue Mar 27, 2018 2:47 am
by hellolori
Dear all,

As we know, Eliashberg spectrum reads,

Code: Select all

\alpha^{2}F(\omega) = 1/N_F/N_k/N_q \sum_{k,k',\nu} |g_{v}_{k,k'}|^2*\delta(\epsilon_k)*\delta(\epsilon_{k'})*\delta(\omega - \omega_{q,\nu})

I would like to know can I get k-resolved Eliashberg spectrum, \alpha^{2}F(k,k',\omega) from EPW directly by setting some tags? Where can I find the corresponding output?
Can any one help me?

Bests!

Re: How to get k-resolved Eliashberg spectrum

Posted: Tue Apr 03, 2018 2:21 pm
by roxana
Hi,

The anisotropic Eliashberg spectral function is not printed on file. You can make some small changes to the code if you would like to print it on file in one of the following two ways:

1) In eliashberg.f90 "CALL evaluate_a2fij" after "CALL evaluate_a2f_lambda"
You will need to add the lines to print to file in the subroutine evaluate_a2fij

or

2) In the subroutine evaluate_a2f_lambda calculate the anistropic Eliashberg spectral function as it is done in evaluate_a2fij and print it to file from there.

Best,
Roxana

Re: How to get k-resolved Eliashberg spectrum

Posted: Thu Apr 05, 2018 3:12 pm
by hellolori
roxana wrote:Hi,

The anisotropic Eliashberg spectral function is not printed on file. You can make some small changes to the code if you would like to print it on file in one of the following two ways:

1) In eliashberg.f90 "CALL evaluate_a2fij" after "CALL evaluate_a2f_lambda"
You will need to add the lines to print to file in the subroutine evaluate_a2fij

or

2) In the subroutine evaluate_a2f_lambda calculate the anistropic Eliashberg spectral function as it is done in evaluate_a2fij and print it to file from there.

Best,
Roxana


Thank you very much. What I have done is direct modify the subroutine
"evaluate_a2f_lambda" to print the k-resolved Eliashberg spectrum. After that, I test the validity of my method by summing up the output and the resultant Eliashberg spectrum is in lines with the "prefix.a2f" generated by EPW. I did not modify the subroutine "evaluate_a2fij" as you mentioned above. I wonder whether the above two methods are equivalent?

Re: How to get k-resolved Eliashberg spectrum

Posted: Fri Apr 06, 2018 11:17 am
by roxana
Hi,

Sounds good! Yes, the two methods are equivalent.

Best,
Roxana