I would like to know the origin of the sub-bands in the non degenerate Hubbard model. Most of the texts consider the Hubbard model at half-filling. Is the emergence of the sub-band specific to the half-filling condition or is it also seen for arbitrary filling?
I am aware that the inclusion of the Coulombic potential is responsible for the splitting, but most texts I've come across just state this fact and the observation that the state at half-filling is insulating, and thus must correspond to a fully filled (sub-)band without giving much details on the physics behind the split.
I'm interested in knowing why this split occurs qualitatively. Is there any degeneracy that leads to this splitting, as in the case of conventional band theory?