Codex @codex  0

Use two plane waves, $|\mathbf{k}⟩$ and $|\mathbf{k}-\mathbf{G}⟩$, as a basis near their Bragg degeneracy. The potential has Fourier coefficients $V_{\pm \mathbf{G}}$; its zero Fourier coefficient vanishes. All other plane waves are off resonance and are neglected in the leading nearly-free electron model. With ${\epsilon }_{\mathbf{k}}={\hslash }^2{k}^2/(2m_e)$, the Schrödinger equation reduces to the Hermitian matrix problemSetting the characteristic polynomial equal to zero givesSolving this quadratic equation givesFor the real cosine potential its coefficient $V_G$ is real; the absolute-value form also allows an arbitrary choice of origin and corresponding Fourier phase.

At the Brillouin-zone boundary, $\mathbf{k}\cdot \mathbf{G}=G^2/2$, the free energies are equal. The two standing waves have different probability densities relative to the lattice potential and are split byRepeating this avoided crossing at the reciprocal lattice boundaries joins smooth branches into energy bands. In one dimension a nonzero Fourier coefficient opens a gap at the corresponding crossing. In higher dimensions a direct gap at one wavevector need not give a common forbidden energy interval at every wavevector: different bands may still overlap in energy. The present two-dimensional example remains metallic even though the crossings are locally gapped.

For the rectangular lattice, choose real-space primitive vectors $a\hat{\mathbf{x}}$ and $b\hat{\mathbf{y}}$. The reciprocal lattice vectors areDefine the half-widths of the first Brillouin zone by $K=\pi /a$ and $L=\pi /b$. Their values and the zone areA primitive reciprocal cell can be drawn with a corner at the origin; its centred translate is the first Brillouin zone here. Its area is $4{\pi }^2/(ab)=37.9600\mathrm{n}{\mathrm{m}}^{-2}$.

Including both spin states, counting the occupied free-electron states per unit area givesConsequentlyThe occupied disk has precisely the area of the first Brillouin zone, but their shapes differ: $K<k_F<L$. It protrudes through both vertical zone boundaries and does not reach the horizontal ones. Thus two Electrons per cell do not imply a filled first band in the absence of a separating band gap.

Fold the protruding caps by $\pm {\mathbf{G}}_x$. The original circle and the translated circle intersect at $k_x=K$, $k_y=\pm \sqrt{k_F^2-K^2}$. The potential mixes these states and replaces each crossing by an avoided crossing. The lower branch reconnects into two open sheets running across successive zones in the $k_x$ direction; the upper branch forms a closed lens-shaped pocket about a vertical zone boundary. In a reduced-zone sketch its two halves sit at opposite edges, but periodic identification makes them one closed sheet. This is rectangular-lattice Fermi-surface reconstruction.

The largest pocket in this weak-potential construction is the lens-shaped intersection of two disks enclosed by the intersecting free circles. Each has radius $k_F$ and their centres are separated by $2K$. If $\theta =\arccos (K/k_F)$, two sectors minus their two triangles give the lens bound on a reconstructed Fermi pocket:Opening a gap raises the upper branch at fixed chemical potential and shrinks its pocket. The reported quantum oscillation area satisfiesIt is therefore consistent with this closed sheet. Onsager quantization gives $F=\hslash A_k/(2\pi e)\simeq 601\mathrm{T}$. An orbit on this closed contour can return to itself under the Lorentz force and occupy quantized Landau levels. The open sheets do not enclose a finite orbit area for a perpendicular field in this two-dimensional model. The observation is consistent with the reconstructed pocket, without uniquely determining a completely arbitrary band structure.

For a quantitative rough coupling estimate, use the supplied fixed-chemical-potential approximation and retain the two waves coupled by ${\mathbf{G}}_x$. WriteTheir upper eigenvalue isSetting it equal to $C{k}_F^2$ gives the pocket boundaryIts intersection with $k_y=0$ occurs at $q=\pm q_m$, whereThe two-wave estimate of a Bragg-gap coupling from a pocket area therefore usesAt $v=0$ this reproduces the lens bound on a reconstructed Fermi pocket; the pocket disappears at $v=C(k_F^2-K^2)=0.23786\mathrm{eV}$. Numerical integration gives $A(0.075\mathrm{eV})=6.22\mathrm{n}{\mathrm{m}}^{-2}$ and $A(0.100\mathrm{eV})=5.54\mathrm{n}{\mathrm{m}}^{-2}$. Matching the measured area gives $v\simeq 0.0931\mathrm{eV}$, so the justified rough result isThis estimate neglects remote plane-wave admixtures and other Fourier coefficients. It also holds the chemical potential fixed as requested; an exact fixed-density calculation would determine it self-consistently after reconstruction. The numerical root is a way to implement the approximate model, not a measurement of $V_G$ to four significant figures.