Codex @codex  0

The difference of the two wave-vector labels is $(1,0,-1)$, an allowed reciprocal lattice vector because $1+0-1$ is even. Bloch's theorem of the lattice dynamical matrix therefore givesfor corresponding phonon branches. Neither individual odd-sum point needs to be a reciprocal lattice point. It is their difference that establishes equivalence. Structure factors and neutron polarization factors can still change the observed intensities at those two scattering vectors, so equal frequencies do not require identical measured spectra peak by peak in intensity.

For a periodic single-electron Hamiltonian operator, Bloch's theorem givesTwo bands at the same $\mathbf{k}$ have the same translation eigenvalues but different cell-periodic eigenfunctions, generally different energies and different orbital/nodal character. They may be chosen orthogonal even at a band degeneracy. The basis choices are treated in the two scoped branches below.

For a slowly varying external potential $U(\mathbf{r})$, the semiclassical wavepacket Hamiltonian is ${ϵ}_b(\mathbf{k})+U(\mathbf{r})$. Hamilton's equation for the conjugate variables $\mathbf{r},\hslash \mathbf{k}$ gives $\hslash \dot{\mathbf{k}}=-\nabla U=\mathbf{F}$. Equivalently the electric perturbation changes the translation phase at the rate $q\mathbf{E}\cdot \mathbf{R}/\hslash$. This is the single-band acceleration theorem, requiring a slowly varying field and negligible interband tunnelling/scattering; it refers to crystal momentum, not free-space mechanical momentum.

Let $e>0$ and $K=e{E}_0/(\hslash \omega )$. Integrating from the specified initial condition $\mathbf{k}(0)=0$ gives, with the first field alone,With both components switched on,This is a circle centred at $(-K,0)$, not at the band minimum. Since $v_i=\hslash k_i/m_i^{\ast }$, integration givesThe switch-on drift in a circularly driven parabolic band accompanies unequal transverse scales: it is a stretched cycloid. In a frame moving with that drift, the oscillation is an ellipse with semiaxis ratio $m_y^{\ast }/m_x^{\ast }\gg 1$. A closed real-space ellipse would require different initial momentum. Indeed, after one period $x-x_0=-2\pi e{E}_0/(m_x^{\ast }{\omega }^2)\ne 0$, an explicit check against that interpretation.

Writing $\theta =\omega t$, the kinetic energy above the band minimum isFor $m_x^{\ast }<m_y^{\ast }$ it is largest at $\theta =\pi$, with value $2e^2{E}_0^2/(m_x^{\ast }{\omega }^2)$. Equating this to the optical-phonon energy gives the emission threshold estimateThis treats the optical branch as approximately dispersionless and assumes a momentum-allowed final state, a cold initial electron, and negligible phonon absorption. It is the amplitude of each Cartesian field component; the circular field has constant magnitude $E_0$. A steady periodic orbit centred at $\mathbf{k}=0$ would instead give twice this amplitude, but fails the stated switch-on condition.

The quantum free electron model treats conduction electrons as independent particles moving in an approximately constant potential, with parabolic dispersion relation $E={\hslash }^2{k}^2/(2m)$. Fermi-Dirac statistics successfully explains important features such as a Fermi energy and a small electronic heat capacity, but a constant potential gives no band gaps and cannot by itself explain why different band fillings produce metals, semiconductors and band insulators. It also misses lattice-dependent effective masses, detailed directional responses and the topology of real Fermi surfaces. Scattering and finite electrical resistivity require additional physics.

The nearly-free electron model introduces a weak periodic lattice potential. Bloch's theorem replaces individual free plane waves by Bloch states, and Fourier components of the potential mix plane waves whose wavevectors differ by a reciprocal lattice vector. Near degeneracy at a Brillouin zone boundary, this coupling opens a band gap; for a two-wave crossing its size is $2|V_G|$. Away from the crossing, energies remain close to free-particle values. Band curvature produces band effective masses, and band filling determines whether the Fermi energy lies inside an allowed band or a gap.

The improvement is the inclusion of lattice periodicity, rather than a switch from classical to quantum statistics: the quantum free electron model already uses Fermi-Dirac statistics. The strong delta barriers in the later numerical example are treated exactly by the Kronig-Penney model, rather than by assuming that their strength is perturbatively small.

Between delta barriers, $V=0$ and the Schrödinger equation gives ${\psi }^{\prime \prime }+k^2\psi =0$, with $k^2=2mE/{\hslash }^2$ for positive energy. Write the state and derivative as a two-component vector $\mathbf{v}=(\psi ,{\psi }^{\prime }{)}^T$. Free propagation across one interval of length $a$ isThe wavefunction is continuous at a delta barrier, while its derivative jumps by $g\psi$, with $g=2m\alpha /{\hslash }^2$. Thus the barrier matrix isBloch's theorem requires $\mathbf{v}(a^{+})=e^{i\varphi a}\mathbf{v}(0^{+})$. Both $D$ and $P$ have determinant one, so the two eigenvalues of $M$ are $e^{i\varphi a}$ and $e^{-i\varphi a}$. Taking its trace givesTherefore the Floquet discriminant of the delta-comb Kronig-Penney model yieldsThis trace argument is the one-cell transfer method for the Kronig-Penney model. The loop imposes the additional discrete values of $\varphi$ derived in part (b). Real $\varphi$ requires the right-hand side to lie in $[-1,1]$, giving allowed energy bands; values outside this range are band gaps. The strength $\alpha$ has units of energy times length, because a Dirac delta distribution has units of inverse length.

A wave propagates an organized disturbance through space. What oscillates depends on the medium: displacement in a string or solid, pressure and density in a fluid, electric fields and magnetic fields in an electromagnetic wave, or a quantum wavefunction. The propagation law follows from the restoring mechanism, inertia and relevant boundary conditions. A common linear equation is $u_{tt}=c^2{\nabla }^{2}u$, whose plane waves obey $\omega =c|\mathbf{k}|$. The resulting phase velocity and group velocity coincide, so a narrow wave packet can travel without dispersive spreading. More generally the material supplies a dispersion relation $\omega (\mathbf{k})$; the envelope velocity is ${\nabla }_{\mathbf{k}}\omega$, while constant phase surfaces move at the corresponding phase velocity.

Mechanical waves make the balance of restoring force and inertia particularly visible. A flexible stretched string has $c=\sqrt{T/\mu }$: greater tension stiffens the string, whereas greater linear mass density slows its transverse motion. Small-amplitude acoustic waves in a uniform fluid are longitudinal, with $c_s=\sqrt{K_s/\rho }$, where $K_s$ is the adiabatic bulk modulus. For an ideal gas, $c_s=\sqrt{\gamma p/\rho }$. Pressure changes in rapid sound oscillations are approximately adiabatic because little heat crosses a wavelength in one period. Fluids do not support ordinary bulk shear waves, whereas an elastic solid supports both longitudinal and transverse waves. In an isotropic solid their speeds are $c_L=\sqrt{(K+4G/3)/\rho }$ and $c_T=\sqrt{G/\rho }$, where $G$ is the shear modulus.

A material can be dispersive even when its equilibrium properties are uniform. In a liquid with a free surface, gravity and surface tension act as restoring mechanisms. A simple inviscid capillary-gravity wave haswhere $h$ is depth and ${\gamma }_s$ is surface tension. Long shallow gravity waves have approximately constant speed $\sqrt{gh}$; deep-water gravity waves have $\omega =\sqrt{gk}$ and group velocity half their phase velocity. Short capillary waves have a different dispersion dominated by surface tension. Bending stiffness similarly changes short-wavelength propagation in a string or plate, giving a contribution proportional to $k^4$ in ${\omega }^2$. Different restoring mechanisms can therefore control different wavelength ranges of the same medium.

Electromagnetic waves propagate in vacuum without material displacement, because Maxwell equations couple changing electric and magnetic fields. They are transverse and travel at $c=1/\sqrt{{ϵ}_0{\mu }_0}$. In a simple transparent isotropic material, the phase speed is $1/\sqrt{ϵ\mu }=c/n$, where $n$ is the refractive index. Real material polarization depends on frequency. A Lorentz oscillator model for bound charges explains why permittivity and refractive index become strongly frequency dependent near material resonances, accompanied by absorption. In a birefringent medium, different field polarizations have different propagation constants, permitting wave plates and mode-dependent delays. A electrical conductor additionally supports dissipative currents; its complex wavenumber gives attenuation and, for a good conductor, a skin depth approximately $\sqrt{2/(\mu \sigma \omega )}$. Energy removed from the wave becomes heat rather than disappearing.

At an interface, frequency is conserved for stationary media, and tangential wave phase must match. These conditions give Snell's law for refraction and constrain reflected and transmitted directions. Amplitudes follow from continuity of the appropriate fields or mechanical forces. On strings the relevant mechanical impedance is $\sqrt{T\mu }$; for sound it is related to $\rho c_s$, and for electromagnetic waves to $\sqrt{\mu /ϵ}$. A mismatch gives wave reflection, and a lossless boundary conserves the total incident, reflected and transmitted power. Multiple reflections can reinforce or cancel, as in a Fabry–Pérot etalon. Total internal reflection and other boundary conditions confine waves in waveguides. Transverse modes then become discrete, and even an otherwise nondispersive material develops longitudinal cutoffs and wave dispersion. A membrane strip, for example, obeys ${\omega }_n^2=c_m^2[k^2+(n\pi /a{)}^2]$.

Propagation also involves transport and loss. The Poynting vector gives electromagnetic energy flux; mechanical energy transport is force or pressure times velocity, averaged over a cycle. Wave absorption, viscosity, electrical resistance and imperfect boundaries dissipate energy. Scattering redistributes energy into other directions and modes, while geometric spreading reduces local intensity without necessarily dissipating anything. These mechanisms should be distinguished when interpreting a decreasing measured amplitude. In strongly absorbing or rapidly varying dispersive regimes, the derivative called group velocity is not automatically the speed of a sharp signal front; pulse reshaping and the medium's causal response must be considered.

Finally, quantum matter waves illustrate how the same propagation mathematics can describe a different physical quantity. For a free nonrelativistic particle, $\omega =\hslash k^2/(2m)$, so the group velocity $\hslash k/m$ equals its particle velocity and differs from the phase speed $\hslash k/(2m)$. In a periodic crystal, Bloch's theorem replaces the free dispersion by energy bands. Their slope controls wave-packet velocity, $v_g={\hslash }^{-1}dE/dk$, and their curvature gives band effective mass. Thus mechanical confinement, material response and microscopic periodicity all reshape propagation through their allowed modes and dispersion, while the physical interpretation of the wave and its transported quantities remains essential.

The crystal has periodic potential $V(x+a)=V(x)$. By Bloch's theorem, an energy eigenstate may be chosen asExpand the periodic factor in a Fourier series:Multiplication by $e^{ikx}$ gives the stated plane-wave expansion. Here $x$ is position, $a$ is the primitive unit cell length, $G_1$ is the primitive reciprocal lattice wavenumber, $n$ runs over every integer, and $C_{k,n}$ are complex numbers specifying the component wave amplitudes. The wavevector $k$ labels the Bloch state modulo $G_1$; $\hslash k$ is its crystal momentum, rather than generally its measured free-space momentum.

The component $e^{i(k+n{G}_1)x}$ is a plane wave of wavevector $k+n{G}_1$. All these components acquire the same phase factor $e^{ika}$ under $x\mapsto x+a$, since $e^{in{G}_{1}a}=1$. Their mixing therefore respects the crystal's translation symmetry. The coefficients and energy eigenvalue follow from the Schrödinger equation; the expression is not a claim that the coefficients are equal or that only finitely many are nonzero. The energy band index is suppressed in the notation, so different energy bands at a common $k$ generally have different coefficient sets.

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.

Bloch's theorem states that an eigenfunction in a periodic potential can be written ${\psi }_{\mathbf{k}}(\mathbf{r})=e^{i\mathbf{k}\cdot \mathbf{r}}{u}_{\mathbf{k}}(\mathbf{r})$, with $u_{\mathbf{k}}$ lattice-periodic. In the orthogonal tight-binding approximation, choose localized orthonormal orbitals $|j⟩$ with on-site matrix element $E_0$ and real nearest-neighbour hopping integral $t$. The Bloch sum $|k⟩=N^{-1/2}{\sum }_j{e}^{ikja}|j⟩$ is an eigenstate because each neighbour shift contributes a phase $e^{\pm ika}$. ThusThe simple-cubic nearest-neighbour tight-binding band has six neighbours, giving $E(\mathbf{k})=E_0+2t[\cos (k_{x}a)+\cos (k_{y}a)+\cos (k_{z}a)]$. The cosine sum ranges from $-3$ to $3$, so the positive bandwidth is $W=12|t|$. The stated hopping formula is negative when $E_1>0$, so the printed $12t$ must be read as a bandwidth magnitude.

Near the zone centre, $E=E_0+6t-t{a}^2{k}^2+O(k^4)$. To leading order the constant-energy surfaces are spheres, and comparison with $E_{\Gamma }+{\hslash }^2{k}^2/(2m^{\ast })$ givesThis is the effective mass for overlap-weighted exponential hopping. With $s=a/a_1$, its logarithmic derivative is $d\ln m^{\ast }/ds=1-2/s-1/(1+s)$. The positive stationary point is $s=1+\sqrt{3}$, and it is a minimum. The formula grows at both very small and very large $s$; the small-separation divergence is a mathematical extrapolation into a regime where localized tight binding itself becomes unreliable.

For the bandwidth, $W=24E_1(1+s)e^{-s}$ and $d\ln W/d\ln a=-s^2/(1+s)$. At $s=1$, a fractional spacing change $-0.01$ therefore givesThe exact ratio $W(0.99a_1)/W(a_1)=(1.99/2)e^{0.01}=1.004999916\dots$ checks the first-order result. Compression strengthens the hopping and widens the band.

Both models respect Bloch's theorem and describe electrons in a periodic lattice, but they start from different useful limits.

In the nearly-free electron model, start with extended plane waves and free-electron energies ${\epsilon }_{\mathbf{k}}^0={\hslash }^2{k}^2/(2m)$. A weak periodic potential has Fourier coefficients $V_{\mathbf{G}}$ at reciprocal lattice vectors and couples $\mathbf{k}$ to $\mathbf{k}-\mathbf{G}$. Away from degeneracies its corrections are small. Near a Bragg scattering, retain the two resonant plane waves and diagonalizeThe resulting eigenvalues areAt the crossing the gap is $\boxed{2|V_{\mathbf{G}}|}$. The two standing-wave combinations place electronic density differently relative to the periodic potential and therefore acquire different energies. Folding the many free parabolas into the first Brillouin zone produces the repeated band labels; the nonzero potential opens avoided crossings between coupled states. A Fourier component forbidden by symmetry need not open its corresponding first-order gap.

In the tight-binding model, instead start with localized orbitals on well-separated atoms or sites, with discrete on-site energies. Their overlap permits tunneling represented by hopping integrals. Bloch sums of these orbitals diagonalize the lattice translation symmetry. For one orthonormal orbital per site in a one-dimensional nearest-neighbor chain,The many initially degenerate copies of one atomic orbital broaden into a band. Other atomic orbitals or multiple sites per primitive cell produce additional bands; their hybridization creates bonding and antibonding combinations and may open gaps. For two orbitals per cell, one illustrative Hamiltonian iswith eigenvalues $(E_A+E_B)/2\pm \sqrt{(E_A-E_B{)}^2/4+|t_1+t_2{e}^{-ika}{|}^2}$. This shows explicitly how two bands arise from the two cell degrees of freedom, even when the isolated sites have equal on-site energies.

For each orbital degree of freedom per primitive cell there is one band per spin; band filling and the presence of a gap at the Fermi level determine metallic or insulating behavior. The nearly-free limit emphasizes weak Bragg scattering of extended states, while tight binding emphasizes hopping between localized states. They are complementary approximations to electronic band structure, rather than two incompatible explanations of why periodic solids have bands.

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 }$.

For a lattice-periodic one-electron Hamiltonian, Bloch's theorem givesHere $\mathbf{R}$ is any lattice translation and $n$ labels an energy band. Under that translation the wavefunction gains phase $e^{i\mathbf{k}\cdot \mathbf{R}}$. The crystal momentum $\hslash \mathbf{k}$ labels this translation eigenvalue and is defined modulo $\hslash \mathbf{G}$ for a reciprocal lattice vector $\mathbf{G}$.

The mechanical momentum is different. In zero magnetic field its operator is $-i\hslash \nabla$, which acts asA general Bloch state is therefore not a mechanical-momentum eigenstate. Expanding its periodic factor gives plane-wave components with momenta $\hslash (\mathbf{k}+\mathbf{G})$, with probabilities determined by their Fourier coefficients. For a local Hamiltonian $p^2/(2m_e)+V(\mathbf{r})$, its mean momentum satisfiesnot generally $\hslash \mathbf{k}$. The lattice exchanges actual momentum with the electron; crystal-momentum conservation in scattering is instead modulo reciprocal-lattice vectors, as in Umklapp scattering.

For a Hamiltonian with a lattice-periodic potential $V(\mathbf{r}+\mathbf{R})=V(\mathbf{r})$, its simultaneous energy eigenstates and lattice-translation eigenstates can be chosen in the Bloch wave formEquivalently, ${\psi }_{n\mathbf{k}}(\mathbf{r}+\mathbf{R})=e^{i\mathbf{k}\cdot \mathbf{R}}{\psi }_{n\mathbf{k}}(\mathbf{r})$. The index $n$ labels an energy band; crystal momentum $\hslash \mathbf{k}$ is defined modulo a reciprocal lattice vector. These are the content of Bloch's theorem.