Just now I did a test calculation for a zigzag graphene nanoribbon. The results obtained by from the LDA calculations are different from what you have stated, i.e., the Bloch state for the highest occupied eigenvalue of Z point (0, 0, 0.5) is localized at both edges of ribbon.