 Solution

Use one orthonormalized localized atomic orbital $|j⟩$ at each position $ja$. In the nearest-neighbour tight-binding model, translation invariance gives the on-site energy and real hopping integralThe Bloch state isSince $H|j⟩=E_0|j⟩+t|j-1⟩+t|j+1⟩$, shifting the two summation indices givesThe hopping integral is an energy, and its sign sets whether $k=0$ is a minimum or maximum. $E_0$ is the diagonal matrix element in the crystal Hamiltonian, which can differ from an isolated-atom energy.

The orthogonal tight-binding approximation matters. If raw neighbouring orbitals retain a nonzero orbital overlap $s$ in their normalization, projection instead gives a generalized eigenvalue problem and the nonorthogonal nearest-neighbour tight-binding band $E(k)=[E_0+2h\cos (ka)]/[1+2s\cos (ka)]$. The stated cosine band follows after orthonormalization and nearest-neighbour truncation, or after neglecting that overlap metric while retaining the hopping it motivates.