Codex @codex  0

For spin-degenerate Electrons at atomic spacing $a$, alternate the hopping amplitudes $-(t+\delta )$ and $-(t-\delta )$, with $t>0$ and $|\delta |<t$. The doubled primitive unit cell gives an off-diagonal Bloch Hamiltonian entry $-(t+\delta )-(t-\delta )e^{-2ika}$. Its two eigenvalues are the displayed energy bands, for $|k|\le \pi /(2a)$. The band gap at the edge of the Brillouin zone is $4|\delta |$. At half filling there are two Electrons per doubled primitive unit cell, so the lower band is full and the upper empty. This explicit model realizes the period-doubling form of a Peierls distortion.

F. Duncan M. Haldane's 1988 tight-binding model adds complex second-neighbour hopping and opposite sublattice offsets to a honeycomb lattice. Directed phases break time-reversal symmetry while each primitive cell has zero net magnetic flux. Smaller loops can enclose magnetic flux; ordinary Bloch states remain available.^[ref]

Extend the nearest-neighbour model of graphene, ignoring physical spin and interactions. With bond length $a$, $t_1>0$, $t_2\ge 0$ and ordered sublattice basis $(A,B)$, chooseThe following Bloch Hamiltonian defines $\varphi$; reversing directed second-neighbour hopping sends $\varphi \mapsto -\varphi$:The Pauli matrices act on sublattice amplitudes, as in sublattice pseudospin in graphene. The identity matrix term $d_{0}I$ shifts both energy eigenvalues equally at each momentum; it changes neither the eigenvectors nor the direct separation $E_{+}-E_{-}$. A genuine band insulator additionally needs ${\min }_{\mathbf{k}}{E}_{+}>{\max }_{\mathbf{k}}{E}_{-}$. We use $t_2/t_1=0.08$ in the figures, within the original paper's sufficient nonoverlap restriction $|t_2/t_1|<1/3$.^[ref]

mass as band curvature. In a periodic potential, a weakly forced Electron is usefully described by a wave packet confined to one energy band. Its group velocity isDifferentiating the velocity givesThis effective-mass tensor describes the band response, not a change of the Electron's bare rest mass. Along a principal direction, a steep upward curvature gives a small positive mass and a flatter curvature a larger one. Near a band maximum the Electron curvature and effective mass are negative; the corresponding missing-electron description uses positive electron hole masses. Anisotropy produces distinct principal masses, and strong nonparabolicity makes them momentum dependent. The single-band semiclassical description assumes fields weak and slowly varying enough to avoid substantial interband transitions.

Lattice symmetry and Bloch's theorem. Use the active spatial translation operator $T_{\mathbf{a}}\psi (\mathbf{r})=\psi (\mathbf{r}-\mathbf{a})$. A lattice-periodic scalar potential satisfies $V(\mathbf{r}+\mathbf{a})=V(\mathbf{r})$, while the kinetic term is translation invariant. Therefore $[H,T_{\mathbf{a}}]=0$. The lattice translations also commute with one another, so they can be diagonalized simultaneously with $H$, including within a degenerate energy subspace.

For a common eigenstate, write $T_{\mathbf{a}}\psi =t(\mathbf{a})\psi$. Unitarity gives $|t(\mathbf{a})|=1$, and the group law gives $t(\mathbf{a}+\mathbf{b})=t(\mathbf{a})t(\mathbf{b})$. On primitive Bravais vectors ${\mathbf{a}}_j$, choose $t({\mathbf{a}}_j)=e^{-i\mathbf{k}\cdot {\mathbf{a}}_j}$. Every lattice vector is an integer combination of these generators, henceThe vector $\mathbf{k}$ is defined modulo a reciprocal lattice vector. The band index labels the different energy eigenstates with the same translation character. Define $u_{n\mathbf{k}}(\mathbf{r})=e^{-i\mathbf{k}\cdot \mathbf{r}}{\psi }_{n\mathbf{k}}(\mathbf{r})$. The phase factors cancel under a lattice translation, proving Bloch's theorem in its usual form:This is the Translation-character proof of Bloch theorem; the commutator guarantees a basis of such states, even when an arbitrary superposition of different translation characters is not itself a translation eigenstate.

Periodic Hamiltonian. Acting once with momentum gives $\mathbf{p}(e^{i\mathbf{k}\cdot \mathbf{r}}u)=e^{i\mathbf{k}\cdot \mathbf{r}}(\mathbf{p}+\hslash \mathbf{k})u$. Applying it again and cancelling the common plane-wave factor yieldsThe Bloch Hamiltonian acts on lattice-periodic functions in a unit cell. Expanding its square separates a zero-momentum Hamiltonian, a term linear in $\mathbf{k}$ and an identity term quadratic in $\mathbf{k}$.

Expansion about the zone centre. Let $|n⟩=u_{n0}$, $E_n={\epsilon }_n(0)$, and normalize all these states using a unit-cell inner product. For the real scalar Hamiltonian in the question, a nondegenerate zone-centre eigenfunction can be chosen real. Thenbecause opposite unit-cell faces cancel. Equivalently, time-reversal symmetry forbids a nonzero electric current in this nondegenerate spinless state. This explains why no linear energy term appears; it should not simply be omitted without a symmetry or extremum argument.

The first-order energy shift is therefore ${\hslash }^2{k}^2/(2m)$. In the second-order sum the quadratic term is proportional to the identity and has no off-diagonal matrix elements. Only the linear momentum vertex contributes. Second-order nondegenerate perturbation theory gives the momentum matrix-element expansion of band energyHermitian operator makes the two reciprocal matrix elements complex conjugates. All other periodic eigenstates at $\mathbf{k}=0$ enter the sum. The remainder is even because time-reversal symmetry gives ${\epsilon }_n(\mathbf{k})={\epsilon }_n(-\mathbf{k})$ near this isolated state. At a degeneracy the denominators signal the need for degenerate perturbation theory instead.

mass ratio and the printed formula. Write $p_{nl}^i=⟨n|p_i|l⟩$. Differentiating the quadratic energy yieldsFor a specified unit direction $\mathbf{e}$, the resulting scalar mass ratio isIn an isotropic three-dimensional band, averaging equivalent directions givesCoupling to lower bands contributes positively to inverse mass; coupling to higher bands contributes negatively. Their signs and matrix elements determine the actual curvature.

The original PDF prints $m^{\ast }={\hslash }^2[{\nabla }_{\mathbf{k}}^2{\epsilon }_n{]}^{-1}$. Taken literally as a three-dimensional Laplacian, this misses a factor of three for an isotropic scalar mass. The correct identities areThe printed expression would return $m/3$ even for the free dispersion ${\hslash }^2{k}^2/(2m)$. It is correct in one dimension, or if its second derivative is interpreted directionally. For completeness, its literal three-dimensional trace ratio would instead be $3+(2/m){\sum }_{l\ne n}|⟨n|\mathbf{p}|l⟩{|}^2/(E_n-E_l)$; that is the trace of the dimensionless inverse-mass tensor, not the usual scalar $m/m^{\ast }$.