 Solution

For a lattice-periodic potential, Bloch's theorem permits energy eigenfunctions ${\psi }_{\mathbf{k}}(\mathbf{r})=e^{i\mathbf{k}\cdot \mathbf{r}}{u}_{\mathbf{k}}(\mathbf{r})$, where $u_{\mathbf{k}}(\mathbf{r}+\mathbf{R})=u_{\mathbf{k}}(\mathbf{r})$ for every lattice translation $\mathbf{R}$. The crystal momentum is $\hslash \mathbf{k}$, defined modulo a reciprocal lattice vector. A reciprocal vector $\mathbf{G}$ satisfies $\mathbf{G}\cdot \mathbf{R}\in 2\pi \mathbb{Z}$, so $e^{i\mathbf{G}\cdot \mathbf{r}}$ is lattice-periodic.

Expand the Bloch wave as plane waves in one momentum class, $\psi ={\sum }_{\mathbf{q}}{c}_{\mathbf{q}}{e}^{i\mathbf{q}\cdot \mathbf{r}}$ with $\mathbf{q}=\mathbf{k}+{\mathbf{G}}^{\prime }$. Substitution into the Schrödinger equation, using $U={\sum }_{\mathbf{G}}{U}_{\mathbf{G}}{e}^{i\mathbf{G}\cdot \mathbf{r}}$, gives a kinetic coefficient ${\hslash }^2{q}^2{c}_{\mathbf{q}}/(2m)$ and potential coefficient ${\sum }_{\mathbf{G}}{U}_{\mathbf{G}}{c}_{\mathbf{q}-\mathbf{G}}$. Independence of the plane waves therefore givesThis uses $\mathbf{q}$ for the paper's running plane-wave index, to distinguish it from the fixed reduced Bloch wavevector.

For $G=2\pi /a$, the sine potential has coefficients $U_G=-iV$, $U_{-G}=iV$. At the Brillouin zone boundary, the free states $e^{i\pi x/a}$ and $e^{-i\pi x/a}$ are degenerate. Keeping these two resonant waves in the nearly-free electron model givesThe determinant is $(\beta -E{)}^2-V^2$, so the two energies are $E_{\pm }=\beta \pm |V|$, separated by $2|V|$. With $V>0$, the lower combination has coefficients $(1,-i)/\sqrt{2}$ and density proportional to $1-\sin Gx$; the upper has $(1,i)/\sqrt{2}$ and density $1+\sin Gx$. The lower wave places less weight where the potential is positive. Coupling to nonresonant plane waves gives higher-order corrections beyond this two-wave approximation.

The first Brillouin zone of a primitive cubic lattice is the cube $|k_i|\le \pi /a$. Its boundary along $(100)$ is the face center $(\pi /a,0,0)$, while along $(111)$ it is the corner $(\pi /a,\pi /a,\pi /a)$. ThereforeA spin-degenerate band accommodates two electrons per cell. To be a band insulator, the filled first band must be separated from the second at every wavevector; a direct local splitting alone is insufficient. Under the paper's uniform-boundary-gap estimate, take the lower edge as $E^{(0)}-V$ and the upper edge as $E^{(0)}+V$. The highest lower edge is $3\beta -V$ and the lowest upper edge is $\beta +V$. Thus the indirect gap in a uniformly split zone-boundary model isEquality gives touching rather than a positive gap; for smaller $V$, overlapping energy bands prevent a band insulator. This final threshold is the stated uniform-splitting band-edge estimate. A quantitative three-dimensional potential mixes several plane waves at corners and requires its full band structure rather than assuming the two-wave sine calculation remains exact there.