Thanks for the suggestion. Yes, OpenMX prints the energies in reduced-kpoints sets. I think that idea not 'expanding' seems physically more meaningful. However, the module 'bandgap.py' asks to return the 'eigenvalues' of the full grid. There may have