In the Lorentz oscillator model, each participating bound electron is displaced by an applied electric field against a restoring force. The model replaces complicated interband transitions by a driven damped harmonic oscillator. If $u$ is the electron displacement along the field and the electron charge is $-e$, its equation of motion isHere ${\omega }_0$ describes binding and $\gamma >0$ describes loss of coherent motion. For $E(t)=\mathrm{Re}(E_0{e}^{-i\omega t})$, the complex amplitude isThe induced electric dipole moment is $-eu$, so the electric polarization amplitude for number density $n$ isUsing $D={ϵ}_{0}E+P={ϵ}_0{ϵ}_{r}E$ givesThe paper's $ϵ$ is this dimensionless relative permittivity, rather than the absolute permittivity ${ϵ}_0{ϵ}_r$. Independent classes of bound electrons contribute additive susceptibilities, usually with different resonance frequencies and optical oscillator strengths.

Writing $D_{\omega }=({\omega }_0^2-{\omega }^2{)}^2+{\gamma }^2{\omega }^2$ separates the response intoAt zero frequency the real part is $1+{\omega }_p^2/{\omega }_0^2$ and the imaginary part vanishes. Near resonance the imaginary part has a positive photon absorption peak and the real part changes rapidly, passing through $1$ at ${\omega }_0$. Above resonance the real part can become negative when damping is sufficiently weak; at very high frequency it approaches $1$ from below while the imaginary part tends to zero. Finite damping removes the undamped divergence.

For a nonmagnetic material take the passive square-root branch $N=\eta +i\kappa =\sqrt{{ϵ}_r}$, with $\kappa \ge 0$. A plane wave then has an intensity attenuation coefficient $\alpha =2\omega \kappa /c$. Thus the imaginary permittivity expresses dissipation, whereas the actual photon absorption coefficient is determined by the imaginary refractive index. The normal-incidence power reflectance isA large refractive-index mismatch produces wave reflection even without loss. In the undamped model ${ϵ}_r<0$ between ${\omega }_0$ and $\sqrt{{\omega }_0^2+{\omega }_p^2}$, so the field is evanescent and wave reflection is total in a semi-infinite sample. With damping the response includes both wave reflection and photon absorption; a peak in the imaginary permittivity should not simply be identified with a reflectance peak.

To infer the short-wavelength background from $R=0.35$, assume weak photon absorption and take the physical semiconductor branch $\eta >1$. ThensoThis is an estimate of the valence-electron background in a frequency window above the doped-carrier response and below the effective bound-electron resonance. One reflectance value alone does not determine both parts of a complex dielectric constant. Even for real ${ϵ}_r$, a reciprocal-index branch gives the same $R$; the ordinary background index greater than unity selects the branch used here.

Using the supplied conventional cubic lattice parameter, the participating valence-electron number density isUsing $m=m_e$ for the stipulated free-electron comparison givesAn angular frequency “in eV” means its quantum energy $\hslash \omega$, not $h\omega$. No doped-carrier effective mass should be inserted into this valence-electron count.

The low-frequency limit of the bound-electron contribution obeysConsequentlyThis effective resonance is much larger than the quoted $0.2\mathrm{eV}$ band gap. There is no contradiction: the single oscillator has been assigned the strength of all sixteen participating valence electrons per conventional cell, whereas the band gap is only the lowest interband excitation threshold. A real spectrum distributes optical oscillator strength over many transitions. Equating the oscillator resonance to the lowest threshold would misrepresent the assumed background response. Here “low frequency” means low relative to ${\omega }_0$; the static response of the doped conducting sample also contains its free carriers.

The additional mobile carriers in doped indium antimonide contribute a Drude frequency response, obtained by setting their restoring frequency to zero. At infrared frequencies well below the effective valence resonance, the sum is approximatelyAt the shorter wavelengths the carrier term is small and $R$ approaches the background value. As wavelength increases, the negative free-carrier contribution lowers the real permittivity towards unity. For weak damping, refractive-index matching to vacuum gives the Drude-Lorentz reflectivity minimum atThis minimum is near ${ϵ}_r^{\prime }=1$, not its zero. At the slightly lower screened plasma frequency ${\omega }_{pc}/\sqrt{{ϵ}_b}$, the real permittivity crosses zero; at longer wavelengths it becomes negative and the material becomes strongly reflective. Finite damping prevents a perfectly vanishing minimum and rounds the edge. At still lower frequency the large conducting response also drives $R$ towards unity.

For fixed carrier effective mass and background permittivity, ${\omega }_{pc}\propto \sqrt{n_c}$ and ${\lambda }_{\min }\propto n_c^{-1/2}$. The observed shift of the minimum towards shorter wavelength as doping increases is therefore the expected free-carrier trend. The figure is an original model illustration, not a digitization or quantitative fit of the supplied curves; the carrier effective mass, damping and their density dependence have not been specified.

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.