Problem with number of electrons

[eliephys78]

I am trying to calculate the electron phonon coefficient (lambda) in a 2D structure. I have been supplying the Fermi energy value in the epw input files using the flag : efermi_read=.true.
In the 'crystal.fmt' file, the number of electrons in my system is "22.000". However, when I run the calculations for lambda, I get :

"The Fermi level will be determined with -30.00000 electrons." Why is it -30 and not 22.. I have checked both the electronic and phonon dispersion bands and they look pretty much the same when compared to those obtained directly from QE. Would this "wrong" number of electrons affect lambda values?

Thanks

[eliephys78]

Please,

can anyone respond to my question whenever you guys are free.

Thanks

[carla.verdi]

Hi,

If you set efermi_read=.true. the Fermi level used in EPW is the one you provided. The message 'The Fermi level will be determined with -30.00000 electrons' is a bit misleading, as it is written also in the case efermi_read=.true.
It doesn't look like the wrong value of nelec should affect your result for lambda, since you read the Fermi level in. However, it is strange that you get a value of -30, and it might indicate there are some further problems in your EPW run. As a note, the value of nelec printed is not the one you have in crystal.fmt, but it is strictly speaking nelec - fac * nbndskip (where fac=1 if noncolin=T, otherwise it's 2).

Best,
Carla