Yambo Community Forum

Dear all,

I'm trying to apply GW approximation on DFT calculations obtained from quantum espresso (QE). The system that I'm working on is an insulator, but Yambo is considering it as a metallic system even with "occupations = fixed" option in QE input. Could any one please help me to fix this problem.

The attached is my input and r_setup files.

Thanks in advances,
Reza

Dear Reza,
most probably it happens because the energy gap as calculated by QE is very small (check it). Yambo recalculate occupations and see the system as metallic.
If this is the case you can generate the setup input with high verbosity: and reduce the default value for the occupation threshold: or just insert this variable by hand in your input file.

Best,
Daniele

Dear Daniele,

It didn't work. Thanks anyway for replying.

Best,
Reza

Dear Reza,
I will try to reproduce your problem.
Best,
Daniele

Dear Daniele,

Thank you so much for your kind.

Best,
Reza

Dear Reza,
I've run your qe input, everything converge smoothly, but when running the "nscf" calculation at the end QE says:
As you can see, the lowest unoccupied state is lower in energy than the highest occupied, indication that the system is "metallic", and yambo indeed see it as metallic.
You can have a look to this thread in QE forum.

If you are sure, for some reason (literature etc.) that your system should be semiconductor also at DFT /BLYP level check your input file (structure and other parameters etc.). Consider you are using a very low energy cutoff, where hgh pseudopotentials are quite hard and do require a large energy cutoff (even if this will not solve your problem). Consider also to update your yambo version, you are using a quite old one, and many bugs have been fixed in the new release.


Best,
Daniele