Codex @codex  0

The Poynting vector is the electromagnetic energy flux density,In vacuum it equals $\mathbf{E}\times \mathbf{B}/{\mu }_0$. For a linear nondispersive medium with time-independent constitutive coefficients, Maxwell's equations and the identity $\nabla \cdot (\mathbf{E}\times \mathbf{H})=\mathbf{H}\cdot \nabla \times \mathbf{E}-\mathbf{E}\cdot \nabla \times \mathbf{H}$ give Poynting's theorem,The source term is work done by the fields on free charges; for an Ohmic conductor it is the positive local Joule heating rate $\sigma E^2$. Integrating over a volume states that loss of stored field energy plus outward energy flow equals energy transferred to matter. In dispersive media, such as a plasma, material motion contributes to stored energy and the simple nondispersive expression cannot be used indiscriminately.

For a vacuum plane wave, $B=E/c$ and $\mathbf{E}$, $\mathbf{H}$ and the propagation direction form a right-handed triad. With real peak field $E_0$, the time-averaged wave intensity isMore generally, complex peak amplitudes give $⟨\mathbf{S}⟩=\frac{1}{2}\mathrm{Re}({\mathbf{E}}_0\times {\mathbf{H}}_0^{\ast })$. Along a transmission line, the cross-section integral of the longitudinal Poynting vector is the transmitted electrical power $VI$. Energy flows through the electromagnetic fields in the surrounding dielectric, even though the conductors guide the currents.

Magnetic field and force. A steady electric current produces a magnetic field. In vacuum the Biot-Savart law givesThe Lorentz force on the charge carriers transfers force to the material carrying them. The resulting force on a current element of another circuit is $d{\mathbf{F}}_2=I_{2}d{\mathbf{\ell }}_2\times {\mathbf{B}}_1$, or, for a continuous current distribution, the magnetic force per unit volume is $\mathbf{f}=\mathbf{J}\times \mathbf{B}$.

For two long parallel straight wires a distance $r$ apart, Ampère's law gives $B_1={\mu }_0{I}_1/(2\pi r)$. ThereforeThe force is attractive for parallel currents and repulsive for opposite currents, by the cross product rule. Its units are $\mathrm{N}{\mathrm{m}}^{-1}$. This approximation neglects wire thickness and end effects; a complete circuit includes its return current. Element-by-element forces need not form an isolated Newton's-third-law pair, but complete stationary circuits have balanced mutual forces when the entire system is included.

Energy description. For two linear coupled circuits the magnetic energy isHere $q$ is a mechanical coordinate and $M$ the signed mutual inductance. At prescribed currents, the generalized force is $I_1{I}_2\partial M/\partial q$ if the self-inductances do not vary. This positive derivative is obtained from magnetic coenergy, or by including the electrical sources' work; differentiating $W_B$ with a minus sign at fixed currents would omit that source work. At prescribed flux linkages the force instead equals the negative derivative of stored magnetic energy with those linkages fixed. The geometry thus determines both inductance and mechanical force.

Integral conservation. Electric charge cannot be created by ordinary transport. For a fixed volume $V$, outward current reduces the enclosed charge:Here ${\rho }_e$ has units $\mathrm{C}{\mathrm{m}}^{-3}$ and $\mathbf{J}$ has units $\mathrm{A}{\mathrm{m}}^{-2}$. Applying the divergence theorem to an arbitrary volume gives the local continuity equationIt is local charge conservation: a net outflow must be accompanied by depletion. In a steady state, $\nabla \cdot \mathbf{J}=0$. Integrating over a small circuit junction then gives Kirchhoff's current law. During a transient the junction can accumulate charge and the incoming and outgoing conduction currents need not yet match.

Consistency with Maxwell's equations. In vacuum, take the divergence of the Ampère-Maxwell equation:Gauss's law, $\nabla \cdot \mathbf{E}={\rho }_e/{ϵ}_0$, reproduces the continuity equation. The displacement current term makes these equations consistent for time-dependent charge distributions, such as a charging capacitor. It is not a stream of free charges across the insulating gap. In matter, $\nabla \cdot \mathbf{D}={\rho }_f$ and $\nabla \times \mathbf{H}={\mathbf{J}}_f+{\partial }_t\mathbf{D}$ give the analogous conservation law for free charge; total charge also includes bound charge.

From molecular dipoles to macroscopic fields. A dielectric contains charges which remain bound to its atoms or molecules. An applied electric field displaces positive and negative charge centres and can orient permanent electric dipoles. The electric polarization $\mathbf{P}$ is electric dipole moment per unit volume, in $\mathrm{C}{\mathrm{m}}^{-2}$. Its associated bound charge isat the boundary of the polarized material, where $\hat{\mathbf{n}}$ points outward. A uniform electric polarization has no volume bound charge, but generally has surface bound charge. The electric field remains the force per unit test charge; it includes the fields of both free and bound charges.

Introducing the electric displacement field $\mathbf{D}={ϵ}_0\mathbf{E}+\mathbf{P}$ packages the free-charge source separately:For a homogeneous linear isotropic dielectric,The electric susceptibility ${\chi }_e$ is dimensionless; the permittivity $ϵ$ has units $\mathrm{F}{\mathrm{m}}^{-1}$. This is a constitutive approximation. An anisotropic dielectric requires a permittivity tensor, while a nonlinear dielectric need not have electric polarization proportional to electric field. The field in these equations is the macroscopic field, not necessarily the local field at an individual molecule.

Interfaces. Let $\hat{\mathbf{n}}$ point from medium 1 into medium 2. Integrating Gauss's law over a thin pillbox, and the electrostatic curl equation around a thin loop, gives the electrostatic interface boundary conditionsThus normal electric displacement is continuous when there is no free surface charge, whereas tangential electric field is continuous even if the permittivity changes. Normal electric field generally jumps because the interface carries bound charge. The electric potential is continuous for an ordinary interface without an ideal dipole sheet. These conditions, together with Laplace's equation in each homogeneous charge-free region, determine the field.

Energy and electrical sources. Reversibly assembling free charge requires incremental work $\delta W=\int \mathbf{E}\cdot \delta \mathbf{D}dV$. For a linear dielectric whose response can be reversibly integrated from zero field, the electrostatic energy isThis includes the material polarization response, rather than simply using ${ϵ}_0{E}^2/2$ everywhere. In a capacitor, $W=Q^2/(2C)=C{V}^2/2$. Inserting a dielectric increases capacitance. At fixed charge the stored energy decreases. At fixed voltage the stored field energy increases, but the voltage source supplies $VdQ$; the appropriate mechanical energy including that source changes by $-V^{2}dC/2$. Both cases predict attraction of dielectric into the capacitor. A force calculation must state which electrical quantity is fixed. A lossy or hysteretic dielectric also dissipates energy, so the reversible quadratic expression does not describe its full history.

Shape and depolarization. Suppose the dielectric is surrounded by vacuum and the distant imposed electric field is ${\mathbf{E}}_0$. Bound surface charge produces a depolarization field. For a uniformly polarized ellipsoid, along a principal axis,Here $N_d$ is the depolarization factor of a dielectric ellipsoid. This self-consistent field explains why the same material responds differently in different shapes.

An infinite slab with electric field normal to its faces has $N_d=1$, giving $E_{\mathrm{in}}=E_0/{ϵ}_r$; its opposite surface bound charges reduce the internal field. For electric field parallel to the faces, $N_d=0$, giving $E_{\mathrm{in}}=E_0$. End effects are neglected.

A very long circular rod in an axial electric field has negligible depolarization away from its ends, so $E_{\mathrm{in}}\simeq E_0$. In a transverse electric field the infinite-cylinder solution has $N_d=1/2$, giving $E_{\mathrm{in}}=2E_0/({ϵ}_r+1)$. To see the latter directly, take the applied field along $x$. The interior and exterior electric potentials can be written ${\varphi }_{\mathrm{in}}=-Br\cos \theta$ and ${\varphi }_{\mathrm{out}}=-E_{0}r\cos \theta +A\cos \theta /r$. Continuity of potential and normal electric displacement at $r=a$ gives $B=2E_0/({ϵ}_r+1)$.

For a sphere of radius $a$, $N_d=1/3$ andFor example, write ${\varphi }_{\mathrm{in}}=-Br\cos \theta$ and ${\varphi }_{\mathrm{out}}=-E_{0}r\cos \theta +C\cos \theta /r^2$. Continuity of potential gives $B=E_0-C/a^3$, and continuity of normal electric displacement gives $ϵB={ϵ}_0(E_0+2C/a^3)$. Solving yields $B=3E_0/({ϵ}_r+2)$ and $C=a^3{E}_0({ϵ}_r-1)/({ϵ}_r+2)$. The $1/r^2$ potential is that of the electric dipole moment above. For ${ϵ}_r=1$ all these induced effects vanish; as ${ϵ}_r$ becomes very large the normal-field sphere and slab approach the corresponding conductor field exclusion. The ideal infinitely long rod in an exactly axial imposed field has a different, singular end-free limit.

Definition and equations. A normal mode is a motion in which every component of a linear system has the same time dependence with fixed amplitude ratios and phases. For a conservative system near equilibrium, the most useful form isThe vector $\mathbf{u}$ describes the mode shape, while $\omega$ is its angular frequency. Linear combinations of modes describe arbitrary sufficiently small motions when the modes form a complete basis.

Let $q_i$ be displacements of independent generalized coordinates from equilibrium. Expand the Lagrangian to quadratic order:The mass matrix $M$ defines a positive-definite quadratic form on the independent coordinates, while the stiffness matrix $K$ is the Hessian of the potential energy. Substituting a mode gives the generalized eigenvalue problem for small oscillationsA practical calculation identifies the equilibrium, chooses independent coordinates, computes the quadratic kinetic and potential energies, solves for the eigenvalues and then obtains their eigenvectors. One can equivalently diagonalize the real symmetric matrix $M^{-1/2}K{M}^{-1/2}$. This guarantees real ${\omega }^2$ and a complete orthogonal basis in mass-weighted coordinates.

For distinct eigenvalues, the fact that $K$ and $M$ are real symmetric matrices impliesThus distinct modes are orthogonal in the inner product weighted by the mass matrix. Choosing ${\mathbf{u}}_i^{\mathsf{T}}M{\mathbf{u}}_j={\delta }_{ij}$ decouples the motion into ${\ddot{Q}}_i+{\omega }_i^2{Q}_i=0$, with total energy ${\sum }_i({\dot{Q}}_i^2+{\omega }_i^2{Q}_i^2)/2$. A positive-definite matrix $K$ gives positive eigenvalues and oscillatory stability. Zero modes can describe rigid translations or rotations without restoring force; a negative eigenvalue gives exponential instability rather than an oscillatory mode.

Symmetry. A symmetry preserving both quadratic energies maps modes into modes of the same frequency. Coordinates adapted to symmetry divide the eigenvalue problem into smaller independent blocks. For two identical subsystems, exchanging them gives symmetric and antisymmetric modes. The degeneracy of normal modes means that independent modes share one frequency. Higher-dimensional group representations of a symmetry can enforce such degeneracies; symmetry does not require every pair of modes to be degenerate, and accidental degeneracies can also occur. Within a degenerate eigenspace, orthogonal linear combinations are equally valid mode shapes.

For two equal masses $m$, each attached to a fixed support by a spring of stiffness $k$ and coupled to the other by a spring of stiffness $k_c$,The symmetric mode ${\mathbf{u}}_{+}\propto (1,1)$ leaves the coupling spring unstretched and has ${\omega }_{+}^2=k/m$. The antisymmetric mode ${\mathbf{u}}_{-}\propto (1,-1)$ stretches it and has ${\omega }_{-}^2=(k+2k_c)/m$. Superposing the two gives beats when their frequencies are close.

Continuous systems. The same idea replaces matrices by differential operators and mode vectors by functions satisfying boundary conditions. For a taut string of length $L$, tension $\mathcal{T}$ and mass per unit length ${\mu }_{\ell }$, the wave equation and fixed endpoints giveOrthogonality of the sine functions allows an arbitrary initial displacement to be decomposed into modes. Changing the endpoints changes the spectrum: boundary conditions are part of the physical mode problem.

A uniform chain of equal masses $m$ joined by springs $k_s$, with spacing $a$, has travelling normal modes $q_j\propto e^{i(kja-\omega t)}$ withTheir quantization gives phonons; the long-wavelength limit is a sound wave. Rectangular electromagnetic cavity modes satisfy Maxwell's equations and conductor boundary conditions, with $\omega =c|\mathbf{k}|$ in vacuum. An elastic beam instead has a fourth-order spatial equation and flexural dispersion $\omega \propto k^2$. These examples show how restoring forces and geometry determine the frequency spectrum.

Real systems. Small-amplitude motion is often approximately linear even when the full equations are nonlinear. Large amplitudes produce frequency shifts, mode coupling and transfer of energy between modes. Friction broadens resonances and makes free oscillations decay. A damping matrix proportional to $M$ or $K$ can retain the undamped mode shapes; general damping need not, and the appropriate eigenfrequencies and modes can be complex. External periodic forcing predominantly excites modes near resonance, with amplitude set by coupling to the forcing and dissipation. Imperfections split ideal degeneracies, while radiation or coupling to a surrounding fluid can remove energy from the chosen subsystem. Normal modes remain useful when linearization, the boundary conditions and the treatment of damping accurately describe the motion being measured.

Equilibrium radiation and cavity frequencies. A blackbody absorbs all incident radiation at every frequency and direction. In thermal equilibrium its emission spectrum depends only on its temperature. Radiation inside an isothermal absorbing cavity, sampled through a small aperture that hardly disturbs equilibrium, realizes black-body radiation to a good approximation. Microscopic absorption and emission establish the equilibrium distribution; a perfectly reflecting wall is the convenient idealization used to count the modes once that distribution is specified.

In a rectangular electromagnetic cavity, the wave equation and perfect conductor wall boundary conditions select component wavenumbers $k_x=\pi l/L_x$, $k_y=\pi m/L_y$, $k_z=\pi n/L_z$. Combining the vacuum dispersion relation $\omega =c|\mathbf{k}|$ and $\nu =\omega /(2\pi )$ givesFor example, electric-field components can have the spatial dependenceFor three positive indices, Gauss's law requires $k_x{A}_x+k_y{A}_y+k_z{A}_z=0$, leaving two independent polarizations of an electromagnetic wave.

Index qualification. The original paper describes the indices as positive integers. That describes the three-nonzero-index family but omits valid electromagnetic modes with exactly one zero index, which have one polarization of an electromagnetic wave. For instance, with $n=0$, $l,m\ge 1$, only $E_z$ above remains; it is normal to the top and bottom faces and vanishes at the side walls. Two zero indices give no nontrivial cavity field. The omitted one-zero-index family does not affect the leading three-dimensional mode count, but is essential to the thin-cavity limit below. Retaining the paper's strict positivity restriction for every mode would instead give a thickness gap and would not yield its requested two-dimensional power law.

Three-dimensional mode count. A lattice cell in positive $\mathbf{k}$ space has volume ${\pi }^3/V$, where $V=L_x{L}_y{L}_z$. In a continuum approximation, a shell of radius $k$ and thickness $dk$ contains its positive octant, multiplied by two polarizations:Using $k=2\pi \nu /c$ givesFor a finite cavity the exact density of states is discrete; this formula is its bulk continuum limit, appropriate when many modes contribute and wall corrections are small.

Mean photon number. A mode is a quantum harmonic oscillator with thermal excitation energies $sh\nu$, $s=0,1,2,\dots$. The zero-point energy does not enter its thermal occupation. Put $q=e^{-h\nu /(k_{B}T)}$. Its partition function and mean occupation areThis is the Bose-Einstein distribution with zero photon chemical potential: equilibrium absorption and emission do not conserve photon number.

Thermal energy. Multiply the density of states by the mean excitation energy $h\nu ⟨s⟩$ and integrate:The substitution $x=h\nu /(k_{B}T)$ isolates all the dependence on cavity dimensions and temperature:Expanding $1/(e^x-1)={\sum }_{r=1}^{\infty }{e}^{-rx}$ gives an integral ${\sum }_{r}6/r^4=6\zeta (4)={\pi }^4/15$. HenceThis is the thermal energy only, excluding the zero-point contribution. The three-dimensional continuum regime requires $k_{B}T$ large compared with the geometric level spacings in all three directions.

Spectral sketch. The spectral energy density per unit ordinary frequency is Planck's law,At small frequency it rises as $8\pi k_{B}T{\nu }^2/c^3$, while at large frequency it decays exponentially. Raising temperature increases the spectrum at every fixed positive frequency, moves the peak to larger frequency in proportion to $T$, and increases the integrated area in proportion to $T^4$. A plot per unit wavelength has a different peak position because the integration variable changes.

Thin-cavity radiation. Let $A=L_x{L}_y$. The relevant low-temperature window isNonzero thickness indices are then thermally suppressed, while many in-plane modes are populated. The $n=0$ family has one polarization of an electromagnetic wave. Its allowed in-plane wavenumbers form a lattice of cell area ${\pi }^2/A$; a quarter of a circular shell therefore containsThus the dimensional crossover of thermal cavity radiation givesThe integral is $2\zeta (3)$ by the same geometric-series argument. A model with a different number of independent polarizations of an electromagnetic wave changes the prefactor but not the $T^3$ scaling. For a cavity with fixed finite $L_x,L_y$, cooling still further eventually freezes its lowest in-plane mode as well: the strict $T\to 0$ thermal energy is exponentially suppressed. The requested area law is the intermediate quasi-two-dimensional regime, not an arbitrarily low-temperature continuum limit.

The closed surface integral sums the outward electric flux: $\mathbf{E}$ is the electric field, $d\mathbf{S}$ its outward vector area, $Q$ the total enclosed charge and ${ϵ}_0$ the vacuum permittivity. Gauss's law relates that flux to enclosed charge for any closed surface. High symmetry can make the field's magnitude constant on a suitable Gaussian surface and fix its direction; without symmetry, flux alone does not determine the local field.

For an isolated conducting sphere, charges reside at its surface and the exterior field is radial. ThusThe conductor has zero interior field. Integrating exterior field energy givesLet the initial charge per drop be $q=4\pi {ϵ}_{0}a{V}_1$. After disconnecting the charging supply, charge is conserved. Volume conservation gives final radius $A=2^{1/3}a$, total charge $2q$ andInitial mutual energy is negligible at the stipulated large separation. HenceThe electrostatic energy increase on equal-drop coalescence is positive. It can be supplied by the decrease in surface energy during coalescence and/or by external mechanical work bringing like charges together. Remaining energy can become motion and heat. If the charge is too large, surface-energy release alone need not suffice to make merging spontaneous; conservation does not assert that it will. Keeping the drops connected to a voltage supply would instead be a different charge-balance problem.

In SI units, with microscopic charge density ${\rho }_e$ and current density $\mathbf{J}$, Maxwell's equations can be writtenThe first relates electric flux to charge. The second states that the net magnetic flux through a closed surface is zero. Faraday's law relates electric circulation to changing magnetic flux. The Ampère-Maxwell equation includes both conduction current and displacement current; the latter is essential to time-dependent fields and local charge conservation.

For fixed contours and spanning surfaces, the integral curl equations are $\oint \mathbf{E}\cdot d\mathbf{l}=-d{\Phi }_B/dt$ and $\oint \mathbf{B}\cdot d\mathbf{l}={\mu }_0{I}_{\mathrm{enc}}+{\mu }_0{\epsilon }_{0}d{\Phi }_E/dt$. For a moving material contour the first must include the motional force term; its total emf is $\oint (\mathbf{E}+{\mathbf{v}}_{\mathrm{loop}}\times \mathbf{B})\cdot d\mathbf{l}$. This prevents confusing a changing circuit area with a partial time derivative of a static magnetic field.

Taking the divergence of the Ampère-Maxwell equation and using electric Gauss's law gives ${\partial }_t{\rho }_e+\nabla \cdot \mathbf{J}=0$. In a source-free vacuum, combining the two curl equations gives ${\nabla }^2\mathbf{E}-{\mu }_0{\epsilon }_0{\partial }_t^2\mathbf{E}=0$, and similarly for $\mathbf{B}$, predicting electromagnetic waves of speed $c=1/\sqrt{{\mu }_0{\epsilon }_0}$. The Lorentz force law then specifies how these fields act on a charge. Macroscopic matter descriptions may instead use $\mathbf{D}$ and $\mathbf{H}$ with free sources and constitutive relations; those should not be silently mixed with the microscopic-source equations above.

An electric dipole is a localized pair of opposite electric charges, or more generally a neutral electric charge distribution with nonzero first electric moment. For charges $-q$ and $+q$ separated by displacement $\mathbf{d}$ directed from negative to positive, the electric dipole moment isA neutral molecule can have a permanent electric dipole moment when its positive and negative charge centres do not coincide. An applied electric field can also produce an induced electric dipole moment by distorting an initially nonpolar electric charge distribution.

At distance $r\gg |\mathbf{d}|$ in vacuum, expand the difference of the two point-charge potentials to first order in the separation. The leading electric potential isThe monopole term cancels because the net electric charge is zero. Taking $\mathbf{E}=-\nabla \Phi$ gives the electric dipole fieldOn the dipole axis its magnitude is $2p/(4\pi {\epsilon }_0{r}^3)$ and its direction is along $\mathbf{p}$; in the equatorial plane the magnitude is $p/(4\pi {\epsilon }_0{r}^3)$ and the field points opposite to $\mathbf{p}$. The $r^{-3}$ decay is faster than the $r^{-2}$ field of an isolated point charge.

In a uniform external electric field, the opposite forces on the two charges cancel in their resultant but form a couple. Its torque isFor a fixed permanent electric dipole moment, the interaction potential energy isDifferentiating with respect to orientation recovers the torque. Parallel alignment minimizes $U$ and is stable against a small rotation; antiparallel alignment maximizes $U$ and is unstable. Rotating a fixed dipole slowly from parallel to antiparallel requires external work $2pE$.

A nonuniform electrostatic electric field can also exert a net force. Expanding the two charge forces gives $\mathbf{F}\simeq (\mathbf{p}\cdot \nabla )\mathbf{E}$, which equals $\nabla (\mathbf{p}\cdot \mathbf{E})$ for a fixed dipole in a curl-free electrostatic field. Thus a suitably aligned dipole is attracted toward stronger field. An induced electric dipole moment requires care in computing energy: when $\mathbf{p}=\alpha \mathbf{E}$ for constant isotropic electric polarizability, building the dipole costs energy and the effective interaction is $U=-\alpha E^2/2$, rather than substituting directly into the fixed-dipole formula. These field, torque and energy effects explain molecular orientation and electric polarization.

For steady electric currents in vacuum, Ampère's circuital law isChoose an oriented closed contour and a spanning surface; the positive enclosed electric current crosses that surface according to the right-hand rule for the contour. First use the source symmetry to establish the direction of $\mathbf{B}$ and which positions have the same magnitude. Then choose the contour so its useful segments are parallel to $\mathbf{B}$ with constant magnitude, and the remaining segments are perpendicular to $\mathbf{B}$ or have vanishing field contribution. The line integral reduces to the unknown magnitude times a known length or sum of lengths.

For example, cylindrical symmetry around an infinite wire allows a circular contour, giving $B(2\pi r)={\mu }_{0}I$. A uniform current sheet instead calls for a rectangle straddling it. The law is always valid in this magnetostatic setting, but by itself does not determine the local magnetic field from an arbitrary enclosed current without these extra symmetry arguments. Here ${\mu }_0$ is the vacuum permeability.

For an infinite straight wire in vacuum, Ampère's circuital law around a circle of radius $R$ gives $B(2\pi R)={\mu }_{0}I$. ThereforeHere ${\mu }_0$ is the vacuum permeability and the azimuthal direction follows the right-hand rule around the electric current. The magnetic field is measured in $\mathrm{T}$.

A scanning tunneling microscope brings a tip and surface, both electrical conductors, close enough to leave only a small vacuum gap. An applied electric potential difference makes occupied electron states on one side able to feed unoccupied states on the other. Although the gap is a classically forbidden potential barrier, quantum tunnelling through its finite width produces an electric current.

For an approximate rectangular potential barrier of gap width $d$, the transmission probability, and thus the electric current at fixed small bias, varies approximately as $e^{-2\kappa d}$. This strong distance dependence makes the device sensitive to changes in separation on the scale of a chemical atom. Feedback control moves the tip vertically to keep the electric current constant during a scan; the required height changes form an image of the surface, with contrast also depending on its electronic states. The finite gap permits transmitted probability current, whereas the semi-infinite step in part (b) does not.

The vacuum wave equation gives ${\omega }^2=c^2[k_g^2+(\pi /a{)}^2]$. The cutoff frequency is $f_c=c/(2a)=1.499\mathrm{GHz}$. In a lossless empty guide, energy travels at group velocity,The larger phase velocity $\omega /k_g=c^2/v_E$ is not the speed of energy transport.

Different measurements reveal different aspects of electronic band structure; a convincing determination combines their constraints.

Angle-resolved photoemission. In angle-resolved photoemission spectroscopy, photons of known energy eject occupied electrons. Measuring their kinetic energy and emission angle reconstructs binding energy and the component of crystal momentum parallel to the surface:Here $\Phi$ is the sample work function when kinetic energy is referenced to the nearby vacuum; an analyzer energy scale requires its own work-function calibration. Momentum is conserved modulo a surface reciprocal lattice vector, neglecting photon momentum. Tracking the peaks maps occupied dispersion relations and the Fermi surface. The measured intensity also contains the matrix element, the Fermi-Dirac distribution and the spectral function, so an absent bright feature does not establish an absent band. For a three-dimensional solid, the perpendicular momentum requires a final-state model. Inverse photoemission spectroscopy probes unoccupied states by injecting electrons and detecting emitted photons.

Optical spectroscopy. Photon absorption and ellipsometry measure interband transitions, optical selection rules and the complex dielectric response. Because optical photon momentum is small on the scale of the Brillouin zone, direct transitions approximately preserve crystal momentum. Their onset constrains a direct band gap; a threshold involving a phonon can reveal an indirect band gap. Peak positions constrain energy differences between bands, but the response also depends on the joint density of states and optical matrix elements, so it does not uniquely determine every $E_n(\mathbf{k})$.

Quantum oscillations and cyclotron resonance. In clean samples at low temperature, the de Haas-van Alphen effect measures oscillations in magnetization and the Shubnikov-de Haas effect measures oscillations in electrical resistivity as the inverse magnetic field is varied. Their period in inverse magnetic field is $1/F$, where Onsager quantization relates $F$ to an extremal Fermi-surface area perpendicular to the field:Rotating the magnetic field maps the geometry of extremal orbits. The thermal damping determines the cyclotron effective mass,with the appropriate magnitude for a hole orbit. Cyclotron resonance supplies the corresponding dynamical mass through ${\omega }_c=eB/m_c$. These are detailed probes near the Fermi level, rather than measurements of all occupied bands. Diffraction supplies the reciprocal lattice needed to label the measured momenta, but diffraction by itself is not a measurement of electronic dispersion.

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.

Write $S=S_0+S_I$, where $S_0$ contains the quadratic kinetic term and mass term, and $S_I=-\lambda \int d^{D}x\varphi (x{)}^3/3!$. ThenLet $Z_0=\int D\varphi e^{i{S}_0}$ and let $⟨\cdot {⟩}_0$ denote the unit-normalized free Gaussian vacuum expectation value. The unnormalized numerator ${\mathcal{N}}_N$ is consequentlyThe Minkowski path integral is understood with the usual scalar Feynman propagator pole prescription, as a formal perturbation series. A real cubic potential energy is unbounded below; no nonperturbative stable vacuum is being assumed here.

Set $\Delta (x-z)=⟨T\varphi (x)\varphi (z){⟩}_0$, $K_i(x)=\Delta (y_i-x)$, and $t=\Delta (0)$, with a common regularization understood for coincident points. The zero-order one-point term has one field and vanishes. At first order, Wick's theorem gives three choices for the contraction with the external field, leaving a tadpole diagram:The second-order term contains seven fields, so it vanishes by the vanishing odd correlators of a free Gaussian field. HenceThe coefficient is $3/3!=1/2$, or symmetry factor $S=2$ when each diagram is weighted by $1/S$. No spacetime integral has been evaluated. The one-point tadpole diagram is shown in the first panel of the following figure; its other panels organize solution.

Let ${\int }_{\mathbf{p}}=\int d^{3}p/(2\pi {)}^3$. The conjugate momentum fields are ${\pi }_{\varphi }={\dot{\varphi }}^{†}$ and ${\pi }_{{\varphi }^{†}}=\dot{\varphi }$. Inserting the given mode expansion of a free field into $i\int d^{3}x({\dot{\varphi }}^{†}\varphi -{\varphi }^{†}\dot{\varphi })$ produces spatial Dirac delta factors. The $b^{†}{c}^{†}$ and $cb$ terms vanish after integration: they contain ${\omega }_{\mathbf{p}}-{\omega }_{\mathbf{q}}$ on ${\delta }^3(\mathbf{p}+\mathbf{q})$, and ${\omega }_{-\mathbf{p}}={\omega }_{\mathbf{p}}$. The number-conserving terms instead add their two frequency factors, cancelling the mode normalization. The un-subtracted ordered expression isThe last term is an infinite vacuum constant, not a finite physical vacuum charge. Choosing the Fock vacuum to have zero Noether charge, equivalently using normal ordering, givesA state containing $n_c$ quanta created by $c^{†}$ and $n_b$ by $b^{†}$ has Noether charge $n_c-n_b$. Thus the two families describe opposite unit Noether charges; with this source's current the $c^{†}$ quantum has positive Noether charge and the $b^{†}$ quantum negative Noether charge. This is an internal Noether charge; identifying it with electric charge requires a specified electromagnetic coupling. The vacuum subtraction of a complex scalar charge leaves all commutators unchanged.^[ref]

Scope of the expansion: a fixed-frequency expansion with time-independent oscillator coefficients is the free-field or asymptotic-field form. For arbitrary interacting $f$, it can be used as canonical oscillator coordinates at a chosen time, rather than as the exact interacting time dependence. For example, $f(s)=m^{2}s+\lambda s^2$ gives $(\square +m^2)\varphi +2\lambda {\varphi }^{\ast }\varphi \varphi =0$, which a generic free plane-wave superposition does not satisfy. The Noether charge and its generator action remain valid in the interacting theory; the symmetry does not depend on using a free Hamiltonian.