Electron distribution using tight binding band structure?

1. Apr 6, 2017

Hi everyone,

I am just wondering how to calculate electron distribution using tight binding band structure for a system like graphene or any other solid.
So the goal is to get |\psi(r)|^2 which \psi is the band state and it is the linear combination of Bloch sum:
\psi={\sum_n,R}{c_nk}exp(ik.R){\phi(r-r_nR)}

c_nk is the eigenvector of the tight binding band structure that we already have them. n is the the orbital position in R unit cell and \phi is \pi orbital.