error in running pp.py

Running EPW
1

Dear all,

I am trying to calculate the superconducting properties of fcc Pb by solving the
isotropic Migdal-Eliashberg equations.

This is the error which I face when I run pp.py.
[code:3lb6kl35]python pp.py
File "pp.py", line 17
print ‘{0}.dyn{1} to {0}.dyn_q{1}.xml’.format(prefix,idyn)
^
SyntaxError: invalid syntax
[/code:3lb6kl35]

I cannot figure out where it is going wrong. I just follow the examples of Pb and made no modifications.

Best regards,
Mahdi

2

Dear Mahdi,

Can you try running the Pb test from QE/test-suite/epw_soc

This should be working.
Let us know if this solves your problem.

Best wishes,
Samuel

3

?Dear Samuel,

Thank you so much. It worked. But I am still having problems with "epw_super" (MgB2) example.
when I run pp.py, I get the following error:

[code:m50t73sd]
cp: cannot stat MgB2.fc: No such file or directory
[/code:m50t73sd]

as I know, MgB2.fc is the result of q2r.x, So, should I perform q2r calculation for this example?

the next problem happens when I run "epw.x". The Temperature has been set as follow In the input files:
[code:m50t73sd]
nstemp = 1
tempsmin = 15.00
tempsmax = 20.00
[/code:m50t73sd]

but the calculations just provides outputs for temperature 15.00, for example MgB2.imag_iso_015.00.
I got confused, what happen for the output files of other temperatures (16.00, 17,00, 18.00, 19.00 and 20.00)?

[quote="sponce":m50t73sd]Dear Mahdi,

Can you try running the Pb test from QE/test-suite/epw_soc

This should be working.
Let us know if this solves your problem.

Best wishes,
Samuel[/quote:m50t73sd]

4

Hi Mahdi,

By setting nstemp = 1, you are telling the code to calculate the superconducting properties at a single temperature, which is tempsmin = 15. If you set nstemp =2, then you will calculate at 15 and 20 K if tempsmin = 15.00 and tempsmax = 20.00. If you want more temperatures than you need to increase nstemp. To understand the meaning of these variables better, you can also take a look at the description of the input variables provided on the EPW website.

Best,
Roxana

5

Hi?I also meet with this problem when I running pp.py.
In my case, it is because I calculate phonon eigenvalue in different q point separately,
So the prefix.dyn file got four more line to describe crystal vector. Deleting that four line in every prefix.dyn or add four f.readline() in pp.py solved my problem
Hope this will help : )

6

Thank you very much.