Codex @codex  0

A classical Thomson scattering frequency window for a harmonically bound Electron requires ${\omega }_0\ll \omega \ll {\omega }_C$, with ${\omega }_0=\sqrt{k/m}$ and recoil scale ${\omega }_C=m{c}^2/\hslash$. There is no parametrically separated window if ${\omega }_0\gtrsim {\omega }_C$, equivalently $k\gtrsim m^3{c}^4/{\hslash }^2$. This is a scale comparison, not a prediction of a sharp transition between scattering mechanisms.

At that binding scale, an oscillator quantum has energy comparable with the Electron rest energy. The nonrelativistic Lorentz oscillator model is then inadequate, so a physical cross-section requires a relativistic quantum description. Compton scattering itself does not switch on discontinuously at $\hslash \omega =m{c}^2$; its recoil corrections become substantial near that scale and tend to the Thomson scattering limit at low photon energy for a free target.

The dynamic electric polarizability relates the induced electric dipole moment at a driving frequency to the applied electric field at that frequency in linear response theory. A complex coefficient can encode both a dispersive phase and absorption. For a charged one-dimensional quantum harmonic oscillator with charge magnitude $e$, mass $m$, and natural frequency ${\omega }_0$, the undamped response isOne derivation uses the oscillator electric-dipole selection rule and first-order time-dependent perturbation theory. The same result follows from the classical Lorentz oscillator model: $m\ddot{z}+m{\omega }_0^{2}z=-eE$ and $d=-ez$ imply $\ddot{d}+{\omega }_0^{2}d=e^{2}E/m$. At exact resonance there is no bounded periodic response without damping; a finite switch-on produces secular growth rather than a finite steady-state value of the pole expression.

An optical oscillator strength is a dimensionless weight assigned to an optical resonance. In an independent-resonance Lorentz oscillator model using a reference electron mass $m_e$ and a specified participating number density $n$,Its normalization must be stated: assigning one unit of total strength per participating Electron means ${\sum }_j{f}_j=1$. A one-resonance approximation with $f_1=1$ represents the whole chosen electron population, rather than asserting that the first allowed interband transition carries all optical strength. The smallest band gap therefore need not equal the fitted resonance energy.

For $\omega \ne {\omega }_0$, isolate the terms proportional to $e^{-i\omega t}$ in the previous answer. Their coefficient isHence the driven electric dipole moment component isFor $t>0$, away from the instantaneous switch, the radiation electric field isThus the dynamic electric polarizability isBelow resonance the driven electric dipole moment is in phase with the electric field; above resonance its sign reverses. The same susceptibility follows from the classical Lorentz oscillator model without damping. At exact resonance the separation into bounded driven and transient terms fails: their apparent poles combine into the finite-time secular term identified in part (c). Without damping there is no bounded steady resonant response, and the first-order state expansion eventually ceases to apply.

Put ${\omega }_0^2=k/m$, where $k$ is the binding spring constant rather than an optical wavenumber. The steady complex amplitude solution of the forced Lorentz oscillator model isThe acceleration complex amplitude is $-{\omega }^2{\mathbf{x}}_0$. HenceThis is the scattering cross-section of a harmonically bound charge. The displacement and induced electric dipole moment are parallel to the incident linear polarization, and integrating their dipole radiation over all directions gives the total cross section above.

For $\omega \ll {\omega }_0$, binding dominates andrecovering Rayleigh scattering with static electric polarizability $q^2/k$. For $\omega \gg {\omega }_0$, inertia dominates: ${\mathbf{x}}_0\simeq -q{\mathbf{E}}_0/(m{\omega }^2)$, so its second time derivative becomes frequency independent andThis is a classical Thomson scattering regime only while $\hslash \omega \ll m{c}^2$. Exactly at $\omega ={\omega }_0$, an undamped driven oscillator has no finite steady-state amplitude; the divergent formula must not be used. Finite damping and radiation reaction regulate the physical resonant response. Coherent contributions from many charges in an extended medium can also require a structure calculation rather than multiplying this single-charge result by a particle count.

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.

Expand the cloud's potential energy about a stable equilibrium: $U(z)=U(0)+U^{\prime }(0)z+U^{\prime \prime }(0)z^2/2+\cdots$. Equilibrium gives $U^{\prime }(0)=0$, stability gives $U^{\prime \prime }(0)=k>0$, and reflection symmetry also removes odd terms. For small displacements, $U-U(0)\simeq k{z}^2/2$ and $F=-dU/dz\simeq -kz$.

Including the electric and damping forces, Newton's second law is $m\ddot{z}=qE(t)-b\dot{z}-kz$. Dividing by $m$ gives the Lorentz oscillator modelThe specified convention uses $2\gamma$ as the damping coefficient per unit mass.

For the steady response, represent the field by $E_0{e}^{i\omega t}$ and the displacement by $\tilde{z}{e}^{i\omega t}$, taking real parts for physical quantities. Substitution givesDefine the nonnegative amplitude $z_0=|\tilde{z}|$ and $\delta =atan2(2\gamma \omega ,{\omega }_0^2-{\omega }^2)$. ThenFor positive $q$, $\varphi =-\delta$. For electrons $q<0$, the phase of a negatively charged driven oscillator is $\varphi =\pi -\delta$: the low-frequency displacement is opposite to the field, rather than in phase with it. A signed amplitude could absorb this extra $\pi$, but the phase convention must then be stated. The induced dipole $p=qz$ has phase $-\delta$ for either charge sign. Transient free oscillations have been omitted because the question asks for steady state.

For the electron cloud, $q/m=-e/m_e$ andThe visible range corresponds to $\omega =2\pi c/\lambda$ between $2.35\times 10^{15}$ and $4.71\times 10^{15}\mathrm{rad}{\mathrm{s}}^{-1}$, below this resonance. The resonant wavelength is about $180\mathrm{nm}$, in the ultraviolet.

The acceleration amplitude is ${\omega }^2{z}_0$. At fixed incident field amplitude, the scattered intensity is therefore proportional toWell below resonance and outside strong damping corrections, $z_0$ is nearly frequency-independent and this becomes the Rayleigh scattering wavelength dependence $I_{\mathrm{sc}}\propto {\omega }^4\propto {\lambda }^{-4}$. The shortest visible wavelengths are consequently scattered much more strongly into directions away from the direct solar beam than the longest ones, producing a blue sky. The simple fourth-power approximation gives a factor 16 between 400 and 800 nm at equal incident intensity. Retaining the undamped denominator with the given resonance raises that factor to about 22.7, since the upper end is not extremely far below resonance. The exact ratio needs the unspecified damping and incident solar spectrum; the stronger scattering of the blue end is the robust conclusion here.

The model describes frequency-dependent scattering, not the detailed perceived color after solar-spectrum weighting and human visual response. It also explains why treating the electrons as free at visible frequencies would miss the bound-charge fourth-power limit.

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 the $e^{-i\omega t}$ convention, so passive photon absorption has positive imaginary permittivity. For a single Lorentz oscillator model resonance, write the photon energy as $E=\hslash \omega$ andHere $S$ has units of energy squared, $\Gamma =\hslash \gamma$ is a damping energy, and ${ϵ}_b$ is an approximately constant background. The corresponding real and imaginary parts areThus the real part starts above the background, rises on the low-energy side of the resonance, drops rapidly through ${ϵ}_r^{\prime }={ϵ}_b$ at $2\mathrm{eV}$, and approaches the background from below at larger energy. The imaginary part starts at zero, has a positive photon absorption peak near $2\mathrm{eV}$, and decays above it. For weak damping its peak is close to the resonance, but not exactly fixed there for arbitrary damping. At asymptotically large energy the corrections scale as ${ϵ}_r^{\prime }-{ϵ}_b\sim -S/E^2$ and ${ϵ}_r^{\prime \prime }\sim S\Gamma /E^3$.

The resonance position is fixed at 2 eV; peak height, width and background are not specified. The figure chooses an explicit finite damping and oscillator strength only to illustrate the shape over 0–5 eV. In particular, a negative real permittivity above resonance is possible for suitable strength and damping, but is not forced by the resonance energy alone. Setting damping exactly to zero would give a singular idealized response instead of a finite-width photon absorption peak.

Choose the $e^{-i\omega t}$ convention. The steady displacement in the Lorentz oscillator model has complex amplitudeThe induced electric dipole moment is $p_{\omega }=q{u}_{\omega }$, so the microscopic electric polarizability is $\alpha (\omega )=q^2/[m({\omega }_T^2-{\omega }^2-i\gamma \omega )]$. For $N$ identical oscillators per unit volume, $P_{\omega }=N{p}_{\omega }={ϵ}_0{\chi }_{\omega }{E}_{\omega }$. ConsequentlyHere the dimensionless quantity called polarizability in the paper is the electric susceptibility; the single-electron polarizability $\alpha$ has different units. The sign of the Electron charge disappears because both the driving force and dipole moment contain $q$.

Multiplying by the complex conjugate of the denominator givesFor positive damping and positive frequency, ${ϵ}_r^{\prime \prime }\ge 0$, as required for passive wave absorption in this Fourier convention. At zero frequency, ${ϵ}_r^{\prime }=1+{\omega }_p^2/{\omega }_T^2$ and ${ϵ}_r^{\prime \prime }=0$. Near a weakly damped resonance the imaginary part has a narrow positive peak, while the real part rises on the low-frequency side and then falls steeply through ${ϵ}_r^{\prime }=1$ at ${\omega }_T$. It can become negative above resonance for sufficiently strong oscillator strength and weak damping. At large frequency, ${ϵ}_r^{\prime }=1-{\omega }_p^2/{\omega }^2+O({\omega }^{-4})$ and ${ϵ}_r^{\prime \prime }\sim {\omega }_p^2\gamma /{\omega }^3$: both approach the vacuum values. The absorption maximum is close to, but need not be exactly at, ${\omega }_T$ for finite damping.

Take $q=-e$ as the signed electron charge. The Lorentz oscillator model is a driven damped bound charge:The first term is inertia, the second is a phenomenological frictional loss of coherent motion to other degrees of freedom, and the third is the linear restoring force from binding to the atom. The driving term is the electric force. The natural angular frequency ${\omega }_0$ represents the electronic transition; this classical oscillator describes its near-resonant linear response. The friction convention is $m\gamma \dot{u}$, so no factor two is hidden in $\gamma$.