EPW Forum

Hello,

I am trying to solve the Eliashberg equations in the normal state (for example, T > Tc), in which way the delta function should be zero. However, the program gives out errors while delta function is approching zero. Is there a good way to implement this?

Thanks.

Zhen

Dear Zhen,

I'm afraid the code was not intended to study the Eliashberg equations in the normal state and there are indeed some numerical divergence very close to Tc.

I'm not sure what would be the best way to treat this. You could try to modify the code to avoid very small values.

Best,
Samuel