Dear All, I am suspicious that a bud exit in epsilon.f90 qe-6.0. The check of nbnd: ======================================================= IF ( REAL(nbnd,DP)*full_occ <= nelec ) CALL errore('epsilon', 'bad band number', 1) ======================================================= I think it should be: REAL(nbnd,DP)*full_occ <= nelec/2.0
Let me know if I missed something... Cheers Jia
_______________________________________________ Pw_forum mailing list Pw_forum@pwscf.org http://pwscf.org/mailman/listinfo/pw_forum