Past exam of the quantum condensed matter course of the University of Cambridge / 2006

Official past-paper archive: www.sp.phy.cam.ac.uk/drp2/home. The direct PDF link has not yet been verified.

A protein is synthesized as a polypeptide whose sequence constrains its possible interactions. Protein folding is the acquisition of a functional three-dimensional structure: local hydrogen bonds organize secondary structure, while packing, electrostatic interactions and burial of nonpolar groups through the hydrophobic effect stabilize the overall fold. Some proteins also form protein disulfide bonds. Folding reduces the chain's conformational entropy, but the relevant balance includes the entropy of the surrounding water; it is the total Gibbs free energy that is minimized at fixed pressure and temperature.

A folding funnel describes a biased but rugged energy landscape: increasingly native-like conformations become favorable, while barriers can trap intermediates. The chain therefore does not have to test every possible conformation in a uniform random search. Molecular chaperones can prevent protein aggregation and help escape kinetic traps; some proteins remain disordered or require partners, so a unique isolated compact fold is not universal.

The ribonuclease A refolding experiment connects this thermodynamic account to a measurement. Unfolding the enzyme with urea while reducing its disulfide bonds removes activity; suitable removal of the denaturant and restoration of disulfide bonds recover activity. Oxidizing while the chain remains unfolded instead produces incorrectly paired disulfides and little activity. Thus sequence, solvent conditions and the route past kinetic traps all matter; restoration of enzyme activity provides a functional test of folding.

An amphiphile contains a solvent-compatible head and a solvent-avoiding tail. In water, aggregation can lower the Gibbs free energy by reducing exposed hydrocarbon area through the hydrophobic effect. This competes with loss of translational entropy, head-group repulsion and the cost of deforming the tails. Above the critical micelle concentration, added amphiphiles mainly enter aggregates rather than continuing to increase the isolated-molecule concentration proportionally.

The amphiphile packing parameter $P=v/(a_0{\ell }_c)$ compares tail volume $v$, preferred head area $a_0$ and tail length ${\ell }_c$. Rough geometric ranges favor spherical micelles for $P<1/3$, cylindrical aggregates for $1/3<P<1/2$, and lipid bilayers for $1/2<P<1$. These are organizing estimates rather than universal sharp phase boundaries. A bilayer can close into a vesicle to eliminate exposed edges; concentration, salt and temperature can change the preferred aggregate.

A simple experimental test varies concentration and measures surface tension: its slope against log concentration changes near the critical micelle concentration, while light scattering detects the appearance of larger aggregates. These observations distinguish the concentration threshold from an assumption that every molecule forms an aggregate at arbitrarily low concentration.

The fluctuation-dissipation theorem relates spontaneous fluctuations in thermal equilibrium to the dissipative part of the linear response to a weak conjugate force. For an overdamped trapped coordinate, the Langevin equation isThe fluctuation-dissipation relation for a Langevin particle fixes the noise strength from the same drag coefficient $\gamma$ that dissipates motion. With a drive proportional to $e^{-i\omega t}$, the linear response function is $\chi (\omega )=1/(K-i\gamma \omega )$. The two-sided angular-frequency power spectral density then obeysIts integral $\int S_{xx}d\omega /(2\pi )=k_{B}T/K$ agrees with the equipartition theorem. Without a restoring force, the same noise-drag relation gives the diffusion-mobility relation $D=k_{B}T/\gamma$.

In thermal calibration of an optical trap, record a bead's spontaneous displacement fluctuations and compare their spectrum with its response to a weak calibrated oscillatory force. The overdamped harmonic-trap power spectrum and the dissipative response must be consistent with the same temperature and drag. The equilibrium and weak-drive conditions are essential: an actively driven bath need not obey this equilibrium relation.

A liquid crystal has some molecular order while retaining liquid-like motion. The distinction between orientational and positional order separates its main phases.

A nematic liquid crystal has long-range orientational order about a nematic director $\mathbf{n}$, with $\mathbf{n}$ and $-\mathbf{n}$ equivalent for apolar molecules, but no periodic density lattice. Its nematic order parameter is $S=⟨(3{\cos }^2\theta -1)/2⟩$, where $\theta$ measures molecular alignment with the director. A cholesteric liquid crystal is a chiral nematic whose director rotates in space, producing a helical pitch while remaining locally nematic.

A smectic phase adds a periodic density modulation in one direction, forming layers with liquid-like motion within each layer. In a smectic A phase, the director is perpendicular to the layers; in a smectic C phase, it is tilted relative to the layer normal. A columnar liquid crystal instead has columns with a two-dimensional positional lattice and comparatively fluid motion along them.

These orders have different experimental signatures. Birefringence under crossed polarizers reveals orientational anisotropy, while X-ray diffraction can detect a smectic layer spacing $d$ through a peak near $q=2\pi /d$. Optical anisotropy alone does not establish a layer lattice.

For equal-size molecules, let $f$ be the mixing Helmholtz free energy per lattice site. Flory–Huggins solution theory reduces to the regular solution model:The first two terms are minus the combinatorial entropy of mixing; the last is the excess cost of unlike contacts. A positive Flory interaction parameter penalizes such contacts. Additive terms linear in composition change the chosen pure-component reference but do not change stability.

A small composition fluctuation at fixed mean composition costs energy proportional to $f^{\prime \prime }$. Thus local stability requires $f^{\prime \prime }>0$, and the spinodal isInside this curve, $f^{\prime \prime }<0$ and arbitrarily small composition fluctuations grow. The minimum of ${\chi }_s$ is at $\varphi =1/2$, giving the critical value ${\chi }_c=2$. At the thermodynamic critical point, $f^{(2)}=f^{(3)}=0$ while $f^{(4)}>0$.

For $\chi =0$, $f$ is strictly convex with its minimum at $1/2$. At $\chi =2$, the central quadratic term vanishes but a positive quartic term remains. At $\chi =2.5$, there are two minima and a central maximum; the spinodal compositions are $(1\pm \sqrt{1-2/\chi })/2\simeq 0.2764,0.7236$. The common-tangent construction gives coexisting compositions ${\varphi }_b\simeq 0.1448$ and $1-{\varphi }_b$, satisfying $\ln [{\varphi }_b/(1-{\varphi }_b)]+\chi (1-2{\varphi }_b)=0$. Between a binodal endpoint and the spinodal, the homogeneous state is locally stable but metastable; global stability requires the convex envelope of $f$.

For two species of chains with the same polymerization degree $N$, the number of independently translating objects per occupied site is reduced by $N$. The symmetric Flory-Huggins polymer blend hasThe symmetry still places the thermodynamic critical point at ${\varphi }_c=1/2$, but nowThe entropy of mixing per monomer is much smaller than for simple liquids, so even a weak unfavorable interaction can cause phase separation for long chains. If the energetic interaction parameter is $\chi =A/T$ with $A>0$, this model gives $T_c=AN/2$; the conclusion relies on that specified temperature law.

With polymer degree $N_A$ and single-site solvent molecules, Flory–Huggins solution theory gives $f/(k_{B}T)=(\varphi /N_A)\ln \varphi +(1-\varphi )\ln (1-\varphi )+\chi \varphi (1-\varphi )$. Expanding the solvent entropy at low concentration gives the dilute-polymer free-energy expansionThe linear term sets a reference chemical potential. The quadratic term controls the leading effective two-segment excluded volume, $v_{\mathrm{eff}}=v_0(1-2\chi )$ for site volume $v_0$.

In a good solvent, $\chi <1/2$ and positive excluded volume swells a dilute chain: the Flory approximation gives $R\sim b{N}_A^{3/5}$ in three dimensions. In a poor solvent, $\chi >1/2$ favors segment contacts and a sufficiently long isolated chain forms a polymer globule, with compact scaling $R\sim b{N}_A^{1/3}$ up to density-dependent factors.

At the theta temperature, the leading two-segment interaction cancels, $\chi =1/2$ in this model, and the long-chain conformation is approximately an ideal chain, $R\sim b\sqrt{N_A}$. The positive cubic term remains. The critical point of a Flory–Huggins polymer solution at finite $N_A$ has ${\chi }_c=\frac{1}{2}(1+N_A^{-1/2}{)}^2$; its limiting value, rather than every finite-chain critical point, is the theta condition.

Write $z_c$ for the lattice coordination number and ${ϵ}_{AA},{ϵ}_{BB},{ϵ}_{AB}$ for contact energies. Counting the excess of unlike contacts relative to the unmixed reference gives the Flory interaction parameterUnlike contacts that are less favorable than the mean like contact give $\chi >0$. For temperature-independent contact energies, $\chi$ is proportional to $1/T$, so heating weakens this energetic preference relative to the entropy of mixing. Effective parameters in real solutions can also include entropic terms, often represented schematically by $A/T+B$; the pure contact-energy model does not explain every possible temperature dependence.

For smooth uncharged spheres, the basic attraction is the Van der Waals force, arising from correlated electric fluctuations. The Hamaker constant $A_H$ incorporates the particles and intervening medium. For equal spheres of radius $a$ separated by a surface gap $h\ll a$, the near-contact van der Waals attraction of spheres isHere positive force means repulsion, and the attractive case has $A_H>0$. At molecular separation, strong short-range repulsion prevents overlap and replaces the continuum divergence with a deep primary minimum. The energy approaches zero from below at large separation; an uncoated, uncharged attractive pair has no electrostatic barrier.

For illustration, $A_H=10^{-20}\mathrm{J}$, $a=100\mathrm{nm}$ and $h=1\mathrm{nm}$ give $|U|\simeq 8.3\times 10^{-20}\mathrm{J}\simeq 20k_{B}T$ at $300\mathrm{K}$. These are chosen scales, not universal material constants. Such a well can trap particles after a thermal encounter and cause particle aggregation, unless another mechanism such as steric stabilization intervenes. The neutral curve in the figure below shows the attraction regularized by an illustrative contact repulsion.

Use the colloidal probe technique: attach a sphere of known radius to the cantilever of an atomic force microscope and approach a chemically matched planar surface in the same liquid as the suspension. A laser reflected from the cantilever measures its deflection $d$; calibration of the spring constant $k_c$ converts this to force, $F=k_{c}d$. The surface separation is reconstructed from piezo motion after correcting for cantilever bending and defining the contact position. Slow approach reduces hydrodynamic drag; residual surface charge, roughness and hydration must be controlled when identifying a dispersion force.

The Derjaguin approximation accounts for the geometry. With interaction free energy per area $W(h)$ between flat plates, $F=2\pi R_{\mathrm{eff}}W(h)$, where $R_{\mathrm{eff}}=a$ for a sphere and plane but $a/2$ for equal spheres. Since $W_{\mathrm{vdW}}=-A_H/(12\pi h^2)$, the sphere-plane force is $-A_{H}a/(6h^2)$; divide it by two to infer the equal-sphere result in the same near-contact approximation. Attractive force gradients can produce a jump to contact, limiting the directly accessible approach curve.

An actual demonstration attached a silica sphere of radius $3.5\mu \mathrm{m}$ to an AFM force sensor and measured its interaction with a plane in sodium chloride solutions. Its double-layer contribution and short-distance deviations show why charge and surface conditions must be measured rather than simply presumed absent.^[ref]

Distinguish the fixed surface charge density ${\sigma }_s$ from the mobile-ion charge magnitude $q=Ze>0$. Let $n_0$ be the bulk number density of each species in a symmetric $Z:Z$ electrolyte, choose bulk potential zero, and write $\epsilon ={\epsilon }_0{\epsilon }_r$. This convention removes the ambiguity between the surface charge and ion valence in the wording.

Constancy of each ionic electrochemical potential gives the Boltzmann distribution in the local electric potential:Positive potential repels positive ions and attracts negative counterions. The coordinate $s$ measures distance from a locally planar surface; the same local relations hold for a spherical potential with its appropriate radial coordinate.

Multiplying each ion number density by its signed charge gives the local charge densityThe fixed surface charge is supplied separately as a boundary condition; it is not part of this mobile charge density in the liquid.

Substitution into Poisson's equation, ${\nabla }^2\psi =-\rho /\epsilon$, gives the Poisson-Boltzmann equationFor a plane, ${\nabla }^2\psi =d^2\psi /d{s}^2$; for a sphere, it is $r^{-2}d(r^2{\psi }^{\prime })/dr$. The bulk conditions are $\psi \to 0$ and $n_{\pm }\to n_0$. At a planar wall with negligible internal field and normal pointing into the liquid, Gauss's law imposes $-\epsilon {\psi }^{\prime }(0)={\sigma }_s$; alternatively, a prescribed surface potential gives $\psi (0)={\psi }_0$.

The Debye–Hückel approximation linearizes the hyperbolic sine, giving $({\nabla }^2-{\kappa }^2)\psi =0$ withThe planar Debye-Hückel double layer has $\psi (s)={\psi }_0{e}^{-\kappa s}$. For an isolated sphere of radius $a$, the Debye-Hückel potential of a charged sphere is $\psi (r)={\psi }_0(a/r)e^{-\kappa (r-a)}$ outside the sphere; the factor $a/r$ accounts for spherical spreading as well as screening.

The Debye–Hückel screening length is the characteristic distance over which the ion cloud suppresses a weak surface potential. Adding salt increases ionic strength and therefore increases $\kappa$, shortening the interaction range. Linearization requires $|q\psi |\ll k_{B}T$, particularly near the surface where the potential is largest. It means the ionic electrostatic energy is small compared with thermal energy; for monovalent ions near room temperature, $k_{B}T/e\simeq 26\mathrm{mV}$. This condition alone does not remove other limitations of point-ion mean-field theory at high concentration or strong correlations.

The force is determined by the total normal stress of the ionic liquid. The ions contribute osmotic pressure because confining an excess number density costs translational entropy. The osmotic stress in a symmetric electric double layer, relative to the bulk reservoir, isThe second term is the normal electrostatic stress. The Poisson-Boltzmann equation makes this combination independent of position between planar surfaces. At the midpoint of identical equally charged surfaces, symmetry gives $E=0$, soCalculating only the repulsion of bare surface charges would omit both the screening ion cloud and its entropic stress. Evaluating the total stress at the midpoint packages their combined force in a particularly simple form.

For equal spheres, assume $h\ll a$, $\kappa a\gg 1$, weak surface potential and weak overlap, with prescribed ${\psi }_0$. Then ${\psi }_{\mathrm{mid}}\simeq 2{\psi }_0{e}^{-\kappa h/2}$ gives $\Pi \simeq 2\epsilon {\kappa }^2{\psi }_0^2{e}^{-\kappa h}$. Integrating once gives the planar interaction energy per area, $W={\int }_h^{\infty }\Pi (s)ds\simeq 2\epsilon \kappa {\psi }_0^2{e}^{-\kappa h}$, and the Derjaguin approximation gives the repulsive sphere force. Combining it with the near-contact van der Waals attraction of spheres yields the DLVO theory approximationPositive force is repulsive. The exponential term can produce a repulsive barrier in the corresponding pair energy; the short-range attraction leads to a primary minimum, while a shallower secondary minimum may occur farther out. Adding salt shortens the repulsion and can destroy the barrier, allowing particle aggregation. Molecular contact requires a repulsive-core correction, and changing the charge-regulation boundary condition changes the electrostatic prefactor.

In an isolated energy band, semiclassical motion satisfies $\hslash \dot{k}=F$ and $v={\hslash }^{-1}dE/dk$. The band is periodic in crystal momentum, so a constant force sweeps the state through equivalent points of the Brillouin zone. Its velocity repeats rather than increasing without bound. For an electron in a constant electric field, $F=-e\mathcal{E}$, and the Bloch oscillation period isFor $E(k)=E_0-2J\cos ka$, $v=(2Ja/\hslash )\sin ka$ and integration gives a bounded position oscillation of amplitude $2J/|F|$ about its mean. A filled band has no corresponding net transport oscillation after summing uniformly over all $k$.

Observation requires scattering times longer than $T_B$ and negligible interband tunnelling. Ordinary crystals often scatter too rapidly before a full zone traversal; larger-period structures and clean atom-optical lattices help. In the cold-atom Bloch oscillation experiment, accelerating a standing-wave lattice supplied a controlled force, and releasing the atoms at different times revealed an oscillating momentum distribution.^[ref]

A field-effect transistor uses a gate voltage to control the conducting channel between source and drain. In a MOSFET, the gate is insulated from the semiconductor: an electric field changes the channel carrier population with little steady gate current. An n-channel device creates an semiconductor inversion layer when the gate sufficiently attracts electrons to a p-type surface. The drain supply provides the output power; the gate controls that flow.

The long-channel MOSFET current law follows from the gradual-channel approximation. With gate capacitance per area $C_{\mathrm{ox}}$, threshold voltage $V_T$, channel width $W$, length $L$, constant charge carrier mobility $\mu$ and $V_{mov}=V_{GS}-V_T>0$, the local sheet charge magnitude is $C_{\mathrm{ox}}(V_{mov}-V(x))$. Integrating the drift current along the channel givesAt $V_{DS}=V_{mov}$, the drain end reaches pinch-off. The ideal long-channel saturation current is $I_D=\mu C_{\mathrm{ox}}W{V}_{mov}^2/(2L)$; this assumes negligible channel-length modulation. A junction field-effect transistor instead changes the channel width through a reverse-biased depletion region.

Experimentally, measure $I_D$ against $V_{GS}$ at fixed drain bias to see gate control, and against $V_{DS}$ at several fixed gate biases to distinguish the approximately linear and saturation regimes. Simultaneously checking gate leakage tests the insulated-gate assumption.

A Wigner crystal is an ordered electron arrangement favored when repulsion overwhelms the kinetic cost of localization. In a three-dimensional electron gas with a neutralizing background, define a mean spacing $a=(3/4\pi n{)}^{1/3}$ and effective Bohr radius $a_B^{\ast }=4\pi \epsilon {\hslash }^2/(m^{\ast }{e}^2)$. The Coulomb-to-kinetic energy ratio of an electron gas scales asThus low density favors the crystal in the degenerate quantum model. Translational symmetry is broken, and excitations include collective lattice vibrations. The melting criterion depends on dimension, screening, disorder and quantum fluctuations, so this scaling argument does not supply a universal numerical transition density.

The electron crystallization on liquid helium experiment provides an important related observation: a dilute two-dimensional electron sheet above liquid helium crystallized on cooling. It was a classical electron system, with thermal rather than Fermi motion competing with Coulomb energy; its coupling parameter is $\Gamma =E_C/(k_{B}T)$. It therefore illustrates Coulomb ordering without being a direct measurement of the zero-temperature three-dimensional quantum transition.^[ref]

In the quantum Hall effect, a two-dimensional electron system at low temperature in a strong perpendicular magnetic field has transverse conductivity of quantized magnitude $|{\sigma }_{xy}|=\nu e^2/h$. On a plateau with negligible longitudinal conductivity, quantum Hall conductivity tensor inversion gives $|{\rho }_{xy}|=h/(\nu e^2)$ and ${\rho }_{xx}\to 0$; the Hall sign depends on the carrier and orientation convention. The filling factor is $\nu =n_{2D}h/(eB)$ when all spin-resolved states are counted.

For the integer quantum Hall effect, filled Landau levels and a mobility gap give integer $\nu$. Disorder-localized states allow the density or field to vary without changing the number of extended occupied channels, producing a quantum Hall plateau rather than quantization only at isolated fields. Quantum Hall edge states provide an equivalent transport picture. In a Hall bar, separate voltage leads measure transverse and longitudinal voltages while a known current flows; sweeping the field reveals both the plateau and the resistance minimum. The von Klitzing-Dorda-Pepper quantum Hall experiment established this integer quantization.

The fractional quantum Hall effect requires electron correlations within a partially filled Landau level. An incompressible correlated state can occur at fractions such as $\nu =1/3$; a noninteracting partially filled level does not account for the associated gap. The Laughlin wavefunction at $\nu =1/m$, odd $m$, contains a factor ${\prod }_{i<j}(z_i-z_j{)}^m$, suppressing close electron approaches while retaining fermionic antisymmetry. Its charged collective excitations can carry magnitude $e/m$; the constituent electrons still have charge $-e$. The Tsui-Störmer-Gossard fractional Hall experiment observed the fractional Hall feature together with a longitudinal-resistance minimum, providing the experimental distinction from an arbitrary rational filling of an ordinary conducting state.

An intrinsic semiconductor has negligible electrically active impurity doping, and thermal excitation across its band gap produces equal electron and hole densities $n=p=n_i$. In an extrinsic sample, semiconductor doping supplies donor impurities or acceptor impurities, so one carrier species dominates over a range of temperatures. Doping does not remove the possibility of low-temperature dopant freeze-out or high-temperature intrinsic behavior.

Assume parabolic isotropic band edges, spin degeneracy two, a single valley, thermal equilibrium and nondegenerate carriers. Write $E_F$ for the Fermi level, and measure conduction-band energy by $E=E_c+ϵ$. Counting momentum states gives the density of states per unit volumeThe low-occupation Maxwell–Boltzmann statistics givesThe substitution $x=ϵ/(k_{B}T)$ and ${\int }_0^{\infty }\sqrt{x}{e}^{-x}dx=\sqrt{\pi }/2$ give the coefficient. Similarly $p=N_v{e}^{(E_v-E_F)/(k_{B}T)}$ and $N_v=2[m_h^{\ast }{k}_{B}T/(2\pi {\hslash }^2){]}^{3/2}$. Multiplication eliminates $E_F$, proving the semiconductor mass-action lawHere the supplied effective masses are treated as the masses entering the effective density of states. A detailed multivalley silicon model would require the appropriate valley degeneracy and density-of-states masses consistently; adding an extra factor without changing the stated model would alter the requested expression.

A pentavalent arsenic impurity supplies a weakly bound excess electron. In the hydrogenic donor state approximation, dielectric screening reduces the Coulomb attraction and the band effective mass replaces the free mass. Hence the donor binding energy isThis is the effective-mass estimate, without a material-specific central-cell correction.

For the temperature plot use donor density $N_D=10^{15}{\mathrm{cm}}^{-3}$ and no acceptors. The charge neutrality with partial donor ionization model puts $E_D=E_c-E_b$ and includes the twofold spin degeneracy of an occupied donor:These equations give three regimes. At low temperature, $p$ is negligible and strong dopant freeze-out gives $n\simeq \sqrt{N_D{N}_c/2}{e}^{-E_b/(2k_{B}T)}$. At intermediate temperature, the donors are almost ionized and intrinsic holes remain rare, so $n\simeq N_D$. At high temperature, $n\simeq p\simeq n_i\propto T^{3/2}{e}^{-E_g/(2k_{B}T)}$. Thus the log-density plot against inverse temperature has a steep intrinsic segment on the left, a nearly horizontal extrinsic plateau, and a freeze-out decline on the right; the $T^{3/2}$ and $T^{3/4}$ factors give slight curvature.

To quantify “approximately constant”, require $0.9N_D\le n\le 1.1N_D$. At $300\mathrm{K}$, the specified model gives $N_c\simeq 4.12\times 10^{18}$, $N_v\simeq 8.61\times 10^{18}$ and $n_i\simeq 2.33\times 10^9{\mathrm{cm}}^{-3}$. The full neutrality equations give the ten-percent plateau intervalThe lower boundary is set by donor ionization and the upper boundary by intrinsic electron-hole production. The large ratio $N_c/N_D$ allows substantial donor ionization even when $k_{B}T<E_b$; using only $k_{B}T\gg E_b$ would incorrectly exclude much of the plateau. The quoted interval depends on the chosen tolerance and the constant-gap band model, rather than being an exact measured range for arsenic-doped silicon.

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.