Codex @codex  0

In a one-dimensional two-wave nearly-free electron model, set ${\epsilon }_B={\hslash }^2{G}^2/(8m)$ and let the coupling magnitude be $U>0$. Near the Bragg plane,Expansion at $q=0$ gives $E_{\pm }={\epsilon }_B\pm U+[{\hslash }^2/(2m)](1\pm 2{\epsilon }_B/U)q^2+O(q^4)$, henceThe lower edge has negative electronic curvature, and the upper edge positive curvature. The last expression is a leading weak-potential approximation, not the exact curvature of the two-wave dispersion.

For free-electron dispersion proportional to $k^2$, states at $\mathbf{k}$ and $\mathbf{k}-\mathbf{G}$ are degenerate on this plane, obtained by equating their squared norms. A nonzero periodic-potential Fourier component can mix them and open an avoided crossing in the nearly-free electron model. First-zone boundaries are selected planes of this kind.

With spin degeneracy, each energy band contains two electrons per primitive unit cell. Divalency therefore permits a filled first band, but only if its maximum lies below the second band's minimum. With negligible ionic periodic potential, the circular Fermi surface and the rectangular zone have equal enclosed area, not equal shape. Electrons occupy some states outside the first zone while some higher-energy first-zone states remain empty. Folded bands overlap, and the material is a metal.

A periodic potential mixes plane waves differing by a reciprocal lattice vector. At a two-wave Bragg point, the nearly-free electron model gives a direct gap $2|V_{\mathbf{G}}|$. Increasing appropriate Fourier coefficients can remove the band overlap. With the first band filled, the system is then a band insulator, provided the temperature is small compared with the insulating band gap in energy units.

The rectangular free-electron band-overlap scale follows from two different points in reciprocal space. Put $C={\hslash }^2/(2m_e)$. Before gap opening, the highest first-zone energy is at the corner, while the lowest second-zone energy is at the midpoint of the nearer face:Thus the downward shift ${\delta }_v$ of the first-band maximum plus the upward shift ${\delta }_c$ of the second-band minimum must overcomeIn a rigid-splitting estimate, an opening $\Delta$ lowers the first band by $\Delta /2$ and raises the second by $\Delta /2$, giving the familiar threshold $\Delta >{\hslash }^2{\pi }^2/(2m_e{b}^2)$. Equality is only band touching, not a positive insulating gap. This expression is the overlap that must be removed, under that explicit band-shift approximation.

Positive direct gaps can coexist with band overlap: a local direct splitting and the global indirect band gap are not interchangeable. The shape and Fourier coefficients of the ionic potential are unspecified, so there is no universal exact critical direct-gap width determined by $a,b$ alone. For a concrete demonstration, a separable potential $2V_x\cos (2Kx)+2V_y\cos (2Ly)$ has axial direct gaps ${\Delta }_x=2|V_x|$ and ${\Delta }_y=2|V_y|$. To first order, the corner's lowest level is $C(K^2+L^2)-|V_x|-|V_y|$, and the two candidate second-band minima are $C{K}^2+|V_x|$ and $C{L}^2+|V_y|$. Therefore, within this model and neglecting higher-order band shifts, nonoverlap requiresEqual axial gaps would give the first-order estimate $\Delta >2C{L}^2/3$, rather than the rigid-splitting estimate. Such a large opening can exceed the quantitative range of weak-potential perturbation theory, but the example identifies the missing model assumption: corner levels feel both couplings. The exact general answer is the positive indirect-gap criterion above; a geometry-only direct-gap threshold would claim more than the supplied potential information supports.

Write $U(x)=-e\varphi (x)$ for the potential energy of an electron in an electrostatic potential $\varphi$, where $e>0$. With $G=2\pi /a$, the coupling in the nearly-free electron model is the Fourier coefficientIf a symbol $V(x)$ already denotes the electron potential energy, as in part (c), it replaces $U(x)$ without another factor $-e$.

At the Brillouin zone boundary the two resonant plane waves have unperturbed energy $E_{\mathrm{ZB}}$. Their Hamiltonian operator in the two-state basis has diagonal entries $E_{\mathrm{ZB}}$ and off-diagonal entries $V_G,V_G^{\ast }$; its eigenvalues are $E_{\mathrm{ZB}}\pm |V_G|$. ThusThis is the local band gap in the two-wave nearly-free electron model.

Define the dimensionless coupling $C=\pi |V_0|/E_{\mathrm{ZB}}$, so that $ϵ=C|\delta |$. Substitution into the energy density from part (b) gives, within the prescribed continuation of the two-wave nearly-free electron model,Its value at $\delta =0$ is understood by continuity. For $C|\delta |\ll 1$,Both electrons per doubled unit cell have been included. The fixed $2G$ harmonic affects the underlying energy bands beyond the two-wave approximation; the expression here is the energy density requested from part (b), rather than an exact solution for all harmonics of the full periodic potential.

The nearly-free electron model expands the Bloch state in plane waves,All basis functions have the same translation eigenvalues because $e^{i\mathbf{G}\cdot \mathbf{R}}=1$. A weak periodic potential mixes them most strongly near Bragg degeneracies and opens band gaps. Delocalized $s,p$ conduction states of simple metals, such as sodium and aluminium, are typical applications; the approximation need not describe tightly bound core states.

Normalize the plane waves in a periodic crystal of area or volume $\mathcal{V}$. The Fourier coefficients of the periodic potential obeyTranslation of the integral by any lattice vector $\mathbf{R}$ leaves it invariant but multiplies it by $e^{-i(\mathbf{k}-{\mathbf{k}}^{\prime })\cdot \mathbf{R}}$. A nonzero result therefore requires $\mathbf{k}-{\mathbf{k}}^{\prime }$ to be a reciprocal lattice vector.

Near an isolated Bragg scattering degeneracy, ordinary nondegenerate perturbation theory fails because the free energies coincide. Retain the resonant plane waves $|\mathbf{k}⟩$, $|\mathbf{k}+\mathbf{G}⟩$ and apply degenerate perturbation theory. Absorb the constant $V_0$ into the energy origin. Since a real potential has $V_{-\mathbf{G}}=V_{\mathbf{G}}^{\ast }$, the projected Schrödinger equation has coefficient matrixA nonzero eigenstate requiresOff-resonant plane waves give higher-order energy shifts omitted in this leading nearly-free electron model.

There are two electrons per primitive unit cell; a spin-degenerate energy band accommodates exactly that number. A semiconductor results when the first band is full and does not overlap the second. In the leading two-wave/four-wave approximation, compare the first-band maximum at $M$ with the second-band minimum at $X$:This is an indirect band gap; a nonzero local splitting at $X$ alone is insufficient. At the onset,With $a=0.3\mathrm{nm}$, ${\hslash }^2/(2m_e)=0.03810\mathrm{eV}\mathrm{n}{\mathrm{m}}^2$ gives $E_X\simeq 4.178\mathrm{eV}$, henceFor $V_{10}$ more negative than this value the estimated indirect band gap is positive; at equality the bands just touch. This is the onset estimate within the source's nearly-free electron model truncation, not an exact critical potential after coupling all plane waves.

The potential has period $2\pi$, so the first Brillouin zone is $-1/2\le k<1/2$. At its boundary the free-particle statesare degenerate, with energyThe relevant Fourier coefficient is supplied by $V_0\cos x$:The term $-V_0\cos 2x$ changes wavevector by two and has no matrix element within this degenerate pair. Thus degenerate perturbation theory gives the matrixand the two band-edge energiesThe nearly-free electron model therefore predicts the lowest band gap

Let $G_n=2\pi n/a$. The free states $|k⟩$ and $|k-G_n⟩$ are degenerate at the Bragg point $k=G_n/2=n\pi /a$. Sincedegenerate perturbation theory restricts the Hamiltonian near this crossing toWriting $k=n\pi /a+\kappa$ and diagonalizing gives the nearly-free electron dispersion near a one-dimensional band gapThus the two continuous bands avoid crossing, and at $\kappa =0$ their separation isThis is the nearly-free electron model. The dispersion relation is the relation $E=E(k)$ between energy and Bloch wavevector.

For the specified potential,Henceand all other nonconstant Fourier coefficients vanish. The gaps are thereforeTheir centre energies are respectivelyThe extended-zone sketch consists of the free-particle parabolas shifted upward by $V_0$, with avoided crossings of these two widths at the listed Bragg points; elsewhere they meet to this order because the corresponding Fourier coefficients vanish.

Write the lattice reciprocal vector as $G=2\pi /a$ in one dimension. The nearly-free electron model starts from ${\epsilon }_k={\hslash }^2{k}^2/(2m_e)$ and couples plane waves separated by G. Near a zone boundary, the two-wave Hamiltonian isAt $k=G/2$ the free waves are degenerate and the periodic potential opens a gap $2|V_G|$. A reduced-zone drawing folds k modulo G into $[-G/2,G/2]$, displaying multiple bands at each reduced wavevector. The original plot shows free folded parabolas, avoided crossings and an allowed vertical optical transition. ^[ref]

Including two spin states, filling a free Fermi sphere gives $n=2(4\pi k_F^3/3)/(2\pi {)}^3=k_F^3/(3{\pi }^2)$. A monovalent BCC cell contributes two electrons, soThe sphere therefore lies inside every first-zone face, since the supplied distance is the nearest boundary.

For the shortest BCC reciprocal vector, $G=2\sqrt{2}\pi /a$. In the zero-potential folded-band limit, the energy of a vertical transition from k along G to the neighbouring free branch is ${\hslash }^2(G^2-2Gk)/(2m_e)$. It decreases over the occupied interval and is smallest at $k=k_F$. Thus the monovalent BCC direct optical threshold isThe final state is empty because $G-k_F>k_F$. Finite weak $V_G$ changes the energies and makes optical matrix elements nonzero; an exactly free electron gas by itself has no interband absorption. The expression is a leading weak-potential threshold, not a claim that every peak equals the onset.

Using the supplied lattice constants and bare-electron free bands gives $E_F=(3.24,2.12,1.85)\mathrm{eV}$ and thresholds $(2.06,1.35,1.18)\mathrm{eV}$ for Na,K,Rb. The source plot's broad higher-frequency features shift toward lower energy in the same order as a increases, qualitatively consistent with the predicted $a^{-2}$ scale. Their peak shapes, heights and exact positions need matrix elements, band corrections and scattering; low-frequency intraband conductivity is a separate contribution. No measured trace is reproduced.

The free electron model replaces the detailed periodic potential by a constant background, treating conduction electrons as mobile particles with approximately parabolic dispersion $E={\hslash }^2{k}^2/(2m)$. It is useful to distinguish the classical Drude model from the quantum free electron gas. They share a simple transport picture, but differ fundamentally in occupation statistics.

With a momentum-relaxation time $\tau$, both lead in a simple isotropic one-band approximation toThis explains approximately Ohmic longitudinal conduction and a carrier-density scale from the Hall effect. The relaxation time is phenomenological: introducing it does not itself predict scattering from defects or phonons. A four-probe resistance measurement separates the sample's longitudinal response from current-contact voltage drops; field and current reversals isolate the transverse signal.

The classical model incorrectly assigns roughly $3N{k}_B/2$ of electronic heat capacity, far too large in a degenerate metal. Fermi-Dirac statistics fixes this: most occupied states cannot change occupation after a small thermal perturbation because nearby low-energy states are already filled. Only a fraction of order $k_{B}T/E_F$ near the Fermi energy participates. For $T\ll T_F$,The quantum model also gives a finite Fermi velocity even at zero temperature and Pauli paramagnetism. Under predominantly elastic scattering, its electronic heat and charge conductivities obey the Wiedemann-Franz law,This is a quantum success; classical equipartition does not give the correct small heat capacity and Lorenz number together.

Its limitations are equally instructive. Without the periodic potential it has no band gaps, so it cannot distinguish a filled-band insulator from a metal or explain semiconductor doping. The simplest single-electron Hall formula cannot explain positive Hall coefficients; band curvature and electron holes, or a multiband transport treatment, are needed. It neglects much Fermi-surface anisotropy and detailed optical interband absorption. Noninteracting electrons alone do not explain ferromagnetism or superconductivity, and lattice thermal behavior needs phonons. The nearly-free electron model, which allows reciprocal-lattice Fourier components of the potential to mix plane waves, provides a controlled first correction near Bragg planes.

Its main success is a simple quantitative baseline for degenerate conduction electrons; its main failure is treating the crystal and interactions as too featureless. A fit of conductivity alone is not a test of the whole model: compare heat capacity, Hall sign, optical response and temperature dependence using independently constrained parameters.

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.

Choose the Fourier series conventionFor a real potential energy, $V_{-1}=V_1^{\ast }$. The symbol $V(G_1)$ means the Fourier coefficient $V_1$, with units of energy, rather than the pointwise value of $V(x)$ at a position $G_1$.

Retain the two plane waves $|1⟩=L^{-1/2}{e}^{ikx}$ and $|2⟩=L^{-1/2}{e}^{i(k-G_1)x}$ in the nearly-free electron model. Their kinetic energies are ${ϵ}_1={\hslash }^2{k}^2/(2m)$ and ${ϵ}_2={\hslash }^2(k-G_1{)}^2/(2m)$. The matrix representation of a linear map for the Quantum Hamiltonian projected onto this two-dimensional basis isThus $\det (H-EI)=0$ givesSolving this quadratic equation yieldsChoosing the zero of potential energy so that its spatial average is $V_0=0$ reproduces the requested expression. If that average has not been removed, both energy bands are shifted by $V_0$; no band gap, group velocity or effective mass is changed by this constant shift.

The two kinetic energies coincide at $k=G_1/2$. Their coupling gives an avoided crossing with band gap $2|V_1|$. The result is exact within the projected two-wave Quantum Hamiltonian. It is an approximation to the full crystal problem, appropriate near this crossing for a weak periodic potential whose other Fourier coefficients do not cause strong mixing with remote plane waves. Small $|V_1|$ alone would not justify ignoring arbitrarily large other Fourier coefficients.

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.

The nearly-free electron model treats a weak periodic potential as a perturbation of plane waves. Near a Bragg plane, states with wavevectors $\mathbf{k}$ and $\mathbf{k}-\mathbf{G}$ are nearly degenerate and must be mixed together. Their two-state Hamiltonian isAt the unperturbed crossing, the band gap is $2|V_{\mathbf{G}}|$. The corresponding standing waves place different charge densities near the ions and between them, giving different potential energies. Far from a crossing, ordinary perturbation theory is adequate; near it, the degenerate two-wave calculation is essential. In three dimensions, gaps at individual zone boundaries do not guarantee a global insulating gap because bands from different directions can overlap.

The tight-binding approximation starts with localized atomic orbitals. Their overlap allows hopping, broadening each atomic energy into a band. For example, a single orbital on a simple cubic lattice with nearest-neighbor hopping gives $E(\mathbf{k})=E_0-2t[\cos (k_{x}a)+\cos (k_{y}a)+\cos (k_{z}a)]$. Small overlap means a narrow band; orbital symmetry and hopping signs determine its extrema and dispersion.

The nearly-free electron model starts instead with free-electron parabolas, folded into the Brillouin zone, and treats the periodic crystal potential as a perturbation. Near a Bragg plane, states differing by a reciprocal vector become degenerate and mix, opening a gap approximately $2|V_{\mathbf{G}}|$. Away from such degeneracies, the energy remains close to the free-electron dispersion. These are complementary limits of the same Bloch-band framework, not different counts of electronic states. Both can give metals or insulators depending on filling and band overlap, but their natural basis functions and useful expansion parameters differ.

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.

In the nearly-free electron model, the free parabola ${\epsilon }_k^0={\hslash }^2{k}^2/(2m)$ is coupled to copies shifted by reciprocal wavevectors $G=2\pi /a$. The first degeneracy is at $k=\pm \pi /a$, where the counterpropagating waves have equal kinetic energy. Retaining their coupling gives two branchesConsequently a nonzero Fourier component $V_G$ opens a gap $\boxed{\Delta E=2|V_G|}$ at the zone edge. The dispersion remains approximately free-electron-like away from the crossing, but the branches avoid crossing and their slopes flatten at the band edges. The horizontal axis can be labeled either wavevector $k$, with edges at $\pm \pi /a$, or physical momentum $p=\hslash k$, with edges at $\pm \pi \hslash /a$.

At the edge, an appropriate choice of spatial origin makes $V_G$ real. The eigenstates are equal-weight standing waves. In a normalization length $L$ containing an integer number of periods, representative real forms areThey come from the sum and difference of $e^{i\pi x/a}$ and $e^{-i\pi x/a}$, up to a global phase. Their densities have maxima and nodes displaced by $a/2$ relative to one another. Both wavefunctions change sign under translation by $a$, as required by their zone-edge Bloch phase, while their densities are lattice-periodic.

For the figure's choice $V(x)=2V_G\cos (2\pi x/a)$ with $V_G<0$, the potential minima occur at $x=na$. The cosine state puts density there and is lower, while the sine state puts a node there and is higher. Reversing the potential's phase relative to the coordinate origin exchanges that assignment. Since the question supplies only the spatial period, it fixes the zone-edge locations and qualitative standing-wave structure, but not a numerical gap or an origin-independent choice of which drawn sine/cosine is lower.

Why degenerate perturbation theory is needed. If several unperturbed states have equal or nearly equal energy, a perturbation can mix them at leading order. In that subspace, diagonalize the perturbation matrix rather than using nondegenerate formulas with small energy denominators. The resulting eigenvectors identify the combinations selected by the coupling, and the eigenvalues give the leading splitting. Other remote states can then be included perturbatively.

For the linear combination of atomic orbitals description of molecular hydrogen, equivalent atomic $1s$ states on the two nuclei supply the nearly degenerate starting functions. In an orthonormal two-state model, equal diagonal energy $E_a$ and real hopping $t$ give $E_a\pm t$ and symmetric/antisymmetric combinations. For actual overlapping atomic functions with overlap $S$ and Hamiltonian coupling $h$, one instead solves a generalized eigenvalue problem:The bonding state has enhanced density between the nuclei and the antibonding state a separating node. Two opposite-spin electrons occupying the bonding molecular orbital describe the simplest closed-shell hydrogen-molecule bond. Electron-electron interaction and nuclear repulsion are needed to calculate a realistic total binding energy; the two-state splitting alone is not the full molecular energy.

In the nearly-free electron model, the analogous degenerate states are plane waves $|\mathbf{k}⟩$ and $|\mathbf{k}-\mathbf{G}⟩$ at a Bragg scattering. A Fourier coefficient $U_{\mathbf{G}}$ of the periodic potential couples them. At exact degeneracy $E_B$, diagonalizing their two-by-two matrix givesThe associated standing waves place electronic density differently relative to the ionic potential. These are the same gap-opening combinations illustrated in Q1(b). The labels of the two starting basis states survive as mixed orbital or plane-wave character, rather than independent crossing eigenstates.

Origin of broad and narrow bands. Spatially extended s/p-derived orbitals have substantial overlap and therefore large hopping amplitudes and broad bands. More localized d orbitals overlap less strongly and produce narrower bands. A multiorbital tight-binding calculation assigns the on-site orbital energies and hopping matrices, forms Bloch sums, and diagonalizes the resulting matrix at each wavevector. Different hopping ranges and orbital symmetry determine the bandwidths and allowed interband coupling. The two-band description retains one broad combination and one nearly flat d-derived combination; a real transition metal contains several d orbitals and more complicated dispersions.

Hybridized eigenvalues and eigenstates. Let $|{\psi }_{\mathbf{k}}⟩=a|{\varphi }_{\mathbf{k}}⟩+b|{\chi }_{\mathbf{k}}⟩$, in an orthonormal basis. Hermitian operator fixes the reverse coupling to be $V^{\ast }$. Projection onto the two components givesSolving this quadratic yields the band hybridization dispersionFor nonzero $V$, one eigenvector representation is $(a,b)\propto (V,E-E_1)$, followed by normalization. Thus the wavefunctions exchange broad-band and d-band character through the avoided crossing; the levels are equal mixtures at the uncoupled crossing.

At $V=0$, the parabola crosses the horizontal level at $k=\pm k_0$, with $k_0=\sqrt{2m{E}_d}/\hslash$. For small nonzero $|V|$, the curves avoid crossing with minimum direct separation $2|V|$. The lower branch follows the parabola below the crossing and flattens toward $E_d$ above it; the upper does the opposite. The lower branch satisfies $E_{-}<E_d$ everywhere and tends to $E_d$ from below as $k\to \infty$, while the upper branch always lies above $E_d$. This is the perfectly flat-band continuum approximation, with zone-boundary intersections excluded as requested.

Fixed particle number and the larger Fermi sphere. Denote the uncoupled broad-band Fermi radius by $k_0$ to distinguish it from the new one $K$. With two spin states, the broad-band number density at $V=0$ is $k_0^3/(3{\pi }^2)$. The stipulated equal population in the flat sector makes the total density twice that value. At zero temperature and any nonzero hybridization, all occupied states fit into the monotonic lower branch below $E_d$, while the upper branch is empty. ThereforeThis is the stated factor increase. It uses fixed total particle number, equal spin counting, an exactly flat unhybridized band and enough available states below any Brillouin-zone restriction. In generic pairs of dispersive bands both hybridized pockets can remain occupied; then Fermi-volume conservation in two hybridized electron bands fixes a sum of volumes rather than this one radius. The exactly flat band is degenerate at $V=0$, so its initial occupation must be specified, as it is here.

Fermi-velocity suppression. Set $\Delta (k)=E_1(k)-E_d$. The exact derivative of the lower branch isAt $k=K$, $E_1(K)=2^{2/3}{E}_d$ and $\Delta (K)=(2^{2/3}-1)E_d>0$. For $|V|\ll \Delta (K)$, the square bracket divided by two is $|V{|}^2/\Delta (K{)}^2+O(|V{|}^4/\Delta (K{)}^4)$. Since group velocity is ${\hslash }^{-1}dE/dk$ and the broad-band reference velocity is $v_F^{(1)}=\hslash k_0/m$,so the flat-band hybridization Fermi-velocity suppression isThe reference velocity is that of the original broad-band Fermi sphere, not the broad-band slope at the larger $K$; omitting the radius change loses the factor $2^{1/3}$.

The low-temperature qualification becomes more stringent as $V\to 0$: the new chemical potential approaches $E_d$ with $E_d-\mu \simeq |V{|}^2/\Delta (K)$, and the relevant thermal energy must be small compared with this hybridization scale. At fixed nonzero temperature the singular flat-band redistribution is not automatically resolved for arbitrarily small coupling. A small but finite narrow-band dispersion can also change the occupation topology.