Codex @codex  0

Equation and boundary conditions. The delta potential has strength $V_0$ with units of energy times length, because the Dirac delta function has units of inverse length. The time-independent Schrödinger equation isThe wavefunction is continuous at the origin. A finite jump in $\psi$ would produce a derivative of a Dirac delta function in ${\psi }^{\prime \prime }$, with no corresponding term in this potential. A jump in ${\psi }^{\prime }$ does produce a delta contribution and is allowed. Integrate through a symmetric interval:For continuous finite $\psi$, the right-hand side vanishes as $ϵ\to 0$. ConsequentlyThis derivative jump is the boundary condition associated with the attractive delta potential.

Unique bound state. For $E<0$, write $\kappa =\sqrt{-2mE}/\hslash >0$. Away from the origin, ${\psi }^{\prime \prime }={\kappa }^2\psi$. Normalizability excludes exponentials growing at either infinity, leavingContinuity gives $A=B=C$. The derivative jump then becomes $-2\kappa C=-2m{V}_{0}C/{\hslash }^2$. For a nonzero state,There is exactly one positive value of $\kappa$, so exactly one negative-energy bound state. An odd solution would have $\psi (0)=0$; the allowed decaying exponentials would then both vanish identically. Normalize the even solution usingChoosing the overall phase real and positive gives the bound state in a delta potential:It peaks at the origin, has a downward cusp there, and decays exponentially on both sides. Its value, but not its slope, is continuous.

Scattering amplitudes and fluxes. For $E>0$, put $k=\sqrt{2mE}/\hslash$ and retain $\kappa =m{V}_0/{\hslash }^2$ as a strength parameter. A unit-amplitude wave incident from the left hasContinuity gives $t=1+r$. The jump condition givesSolving,The probability current of a plane wave of amplitude $A$ is $\hslash k|A{|}^2/m$ in its propagation direction. Since the asymptotic wavenumber is the same on both sides, reflection probability and transmission probability are just $|r{|}^2$ and $|t{|}^2$:They satisfy $R+T=1$, expressing conserved probability current in this real potential.

The right-hand plot uses $\eta =\kappa /k=V_0\sqrt{m/(2E)}/\hslash$, proportional to $V_0$ at fixed incident energy. Then $R={\eta }^2/(1+{\eta }^2)$ rises monotonically from zero to one, while $T=1/(1+{\eta }^2)$ falls from one to zero; they are equal at $\eta =1$. Despite the absence of a classical forbidden barrier, an abrupt attractive interaction reflects a quantum wave. Increasing attraction strengthens the derivative mismatch and can approach complete reflection at fixed positive energy. At fixed $V_0$, increasing incident energy instead makes the well increasingly transparent. The negative-energy bound state and the positive-energy scattering states are different stationary solutions, not an absorptive loss channel.

Let $\rho ={\psi }^{\ast }\psi$ be the probability density. For a real potential energy $V(x)$, the Schrödinger equation and its complex conjugate areDifferentiate $\rho$ and substitute both equations. The potential terms cancel, leavingComparison with the probability continuity equation ${\partial }_t\rho +{\partial }_{x}J=0$ identifiesFor a plane wave of amplitude $A$ and positive wavenumber $k$, $J=\hslash k|A{|}^2/m$, the density times the velocity. Integrating the continuity equation over an interval makes its probability change equal to the incoming current minus the outgoing current.

From the second-edge amplitude in part (c),This uses ${\psi }_{\mathrm{R}}=F{e}^{i{k}_2(x-a)}$. If the transmitted wave is instead written $T{e}^{i{k}_{2}x}$, its coefficient is $\boxed{T/A=i\sqrt{k_0/k_2}{e}^{-i{k}_{2}a}}$; the extra factor only changes the quantum phase convention.

Using the probability current of a plane wave,Thus the flux transmission probability is one. The larger transmitted density $|F{|}^2>|A{|}^2$ compensates for the smaller transmitted velocity $\hslash k_2/m$; equal currents need not have equal wave amplitudes.

Both the incoming and outgoing spherical wavefronts contribute to the path excess. At aperture radius $u$, with axial distances $s$ before and $r$ after the aperture,The outer edge of the $n$th Fresnel half-period zone has $\Delta \ell =n\lambda /2$. HenceTherefore the outer radius is $\boxed{u_2=4.47\times 10^{-4}\mathrm{m}=0.447\mathrm{mm}}$. The paraxial approximation is excellent because $u_2\ll s,r$. Treating the illumination as a plane wave would incorrectly omit the source-side path curvature.

In Fraunhofer diffraction, the relative propagation phase across the aperture is linear in aperture coordinates. If $b$ is its largest transverse size, the neglected quadratic path term must be small: for an observer at distance $R$,For a finite-distance point source a distance $s$ away, also require $b^2/(\lambda s)\ll 1$ if the illumination is to be treated as a plane wave. Numerical factors depend on whether $b$ is a radius or full width and on the permitted phase error; the condition expresses a small Fresnel number.

One can place both source and screen sufficiently far away, or collimate a point source with a lens and observe the aperture in the back focal plane of a second lens. In the latter arrangement rays of each outgoing direction focus at one point, with $X\simeq f{\theta }_x$; the lens removes the need for a physically distant screen. A laser beam approximately uniform across the aperture can supply the incident plane wave. The sketches distinguish an actual far field from its focal-plane equivalent.

A localized quantum particle is described by a wave packet: a coherent quantum superposition of momentum components, rather than a single delocalized plane wave. The Born rule assigns the position probability density $|\psi (x,t){|}^2$. A useful normalized initial Gaussian wave packet isIts mean position and momentum are $x_0$ and $p_0$, and its standard deviations satisfy ${\sigma }_x(0)={\sigma }_0$, ${\sigma }_p=\hslash /(2{\sigma }_0)$. This unchirped initial state saturates the Heisenberg uncertainty relation.

For a free particle, the momentum-space Schrödinger equation multiplies each amplitude by $e^{-i{p}^{2}t/(2m\hslash )}$. The inverse Fourier transform is a Gaussian integral, giving the free Gaussian wave packet evolutionOne way to obtain this expression explicitly is to write the initial momentum amplitude as proportional to $\exp [-{\sigma }_0^2(p-p_0{)}^2/{\hslash }^2-ip{x}_0/\hslash ]$, multiply by the free evolution phase and complete the square in $p$. Taking the modulus squared removes the position-dependent phase and leavesThe expectation value follows uniform classical motion while the density spreads. A narrower initial packet needs a wider momentum distribution and therefore spreads more rapidly. At long times its width grows linearly with $|t|$; its variance grows quadratically, unlike the linear variance growth of ordinary diffusion. The probability continuity equation preserves total probability despite the decreasing peak height. The mean momentum and the mean energy $(p_0^2+{\sigma }_p^2)/(2m)$ remain constant, since the momentum populations do not change.

The same result has a wave interpretation. For a free nonrelativistic particle $\omega (k)=\hslash k^2/(2m)$, so the group velocity is $d\omega /dk=\hslash k/m$, equal to $p_0/m$ at the packet center. The curvature $d^2\omega /d{k}^2=\hslash /m$ makes different momentum components travel with different group velocities, producing dispersion. The phase velocity is $\omega /k=\hslash k/(2m)$ for this energy convention and does not describe the motion of the envelope. The evolving quadratic spatial phase is essential to the spreading; a density snapshot alone does not specify the full quantum state.

A freely spreading Gaussian wave packet remains a pure Gaussian state but develops position-momentum correlations. Its symmetrized covariance is $C_{xp}=⟨\{x-⟨x⟩,p-⟨p⟩\}⟩/2=t{\sigma }_p^2/m$. ConsequentlyEquality in the Robertson-Schrödinger uncertainty relation survives at every time, whereas the simpler product ${\sigma }_x{\sigma }_p$ equals $\hslash /2$ only initially. A chirped Gaussian wave packet can already have $C_{xp}(0)\ne 0$. More generally, $x(t)=x(0)+p(0)t/m$ gives ${\sigma }_x(t{)}^2={\sigma }_x(0{)}^2+2t{C}_{xp}(0)/m+t^2{\sigma }_p(0{)}^2/m^2$. A sufficiently negative initial covariance permits initial contraction before subsequent spreading.

A potential energy changes this behavior. A quadratic potential preserves the Gaussian form because the evolution of position and momentum is linear. For a quantum harmonic oscillator, the Gaussian covariance in a harmonic potential isA coherent state has ${\sigma }_x^2=\hslash /(2m\omega )$, ${\sigma }_p^2=m\hslash \omega /2$ and zero centered covariance. Its width is constant, and its center follows $x_c(t)=x_0\cos \omega t+p_0\sin \omega t/(m\omega )$. In oscillator notation its coherent state label rotates as $\alpha (t)=\alpha (0)e^{-i\omega t}$, apart from the common ground-energy phase of the state vector. A squeezed coherent state generally has a width that breathes at $2\omega$. Thus being Gaussian does not imply either free spreading or constant width; the Hamiltonian operator and initial covariance determine the result.

For a general nonquadratic potential, different parts of the packet experience different nonlinear forces, and an initially Gaussian profile need not remain Gaussian. Ehrenfest theorem gives the exact relations $d⟨x⟩/dt=⟨p⟩/m$ and $d⟨p⟩/dt=-⟨V^{\prime }(x)⟩$. Replacing the last mean by $-V^{\prime }(⟨x⟩)$ is exact for a quadratic potential and approximate for a packet narrow compared with the scale on which the force varies. This connects packet motion to the correspondence principle without treating a quantum state as a definite classical path. Scattering can split a packet into reflected and transmitted components, while confinement can produce oscillations and interference between discrete energy eigenstates. Their relative phases evolve even when their populations remain constant. Localization, spreading and coherent motion all follow from the same linear unitary time evolution.

Build a narrow-band wave packet from plane waves with wavenumbers close to $k_0$:The Taylor expansion of the dispersion relation is $\omega (k)\simeq {\omega }_0+{\omega }_0^{\prime }(k-k_0)$, where ${\omega }_0^{\prime }=(d\omega /dk{)}_{k_0}$. Put $\kappa =k-k_0$ to obtainThe carrier's wave phase travels at the phase velocity ${\omega }_0/k_0$, while the envelope has constant shape at constant $x-{\omega }_0^{\prime }t$. Its group velocity is thereforeThe approximation neglects higher derivatives of the dispersion relation. Those derivatives cause the wave packet to spread or distort over sufficiently long propagation times; the result describes the local envelope speed of a sufficiently narrow band.

In Fresnel diffraction, the propagation wave phase varies quadratically across the aperture and that curvature must be retained. For a transverse scale $w$ at propagation distance $L$, its size is of order $\pi w^2/(\lambda L)$. In Fraunhofer diffraction, that quadratic term is negligible, approximately $w^2/(\lambda L)\ll 1$, while the linear phase encoding the observation angle remains. A lens can also form the far-field angular pattern in its focal plane.

For a normally incident plane wave and aperture transmission $h(x,y)$, the far-field complex amplitude is its Fourier transform:The detected light intensity is proportional to $|\psi {|}^2$. For identical parallel slits, write the transmission as a sum of shifted copies of one slit. The Fourier shift theorem gives one single-slit transform multiplied by the sum of positional phase factors:The single-slit diffraction envelope sets the broad scale, while coherent interference between slit centres produces the finer array pattern. One sums complex amplitudes before squaring, rather than adding independent slit intensities.

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.

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.

The crystal has periodic potential $V(x+a)=V(x)$. By Bloch's theorem, an energy eigenstate may be chosen asExpand the periodic factor in a Fourier series:Multiplication by $e^{ikx}$ gives the stated plane-wave expansion. Here $x$ is position, $a$ is the primitive unit cell length, $G_1$ is the primitive reciprocal lattice wavenumber, $n$ runs over every integer, and $C_{k,n}$ are complex numbers specifying the component wave amplitudes. The wavevector $k$ labels the Bloch state modulo $G_1$; $\hslash k$ is its crystal momentum, rather than generally its measured free-space momentum.

The component $e^{i(k+n{G}_1)x}$ is a plane wave of wavevector $k+n{G}_1$. All these components acquire the same phase factor $e^{ika}$ under $x\mapsto x+a$, since $e^{in{G}_{1}a}=1$. Their mixing therefore respects the crystal's translation symmetry. The coefficients and energy eigenvalue follow from the Schrödinger equation; the expression is not a claim that the coefficients are equal or that only finitely many are nonzero. The energy band index is suppressed in the notation, so different energy bands at a common $k$ generally have different coefficient sets.

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.

Use Fresnel diffraction, since the detector moves along the axis through successive near-field interference extrema. Write $z$ for its distance from the opening and $\rho$ for the integration coordinate within the opening. For a uniform incident plane wave, the on-axis complex amplitude, in the paraxial approximation, is proportional toThus the light intensity relative to the incident value isThe Fresnel number $F=r^2/(\lambda z)$ is an odd integer at a maximum and an even integer at a minimum. It decreases as the detector moves away. Adjacent extrema therefore differ by one:ConsequentlyAs a check, $F=5$ at $z=1\mathrm{m}$ and $F=4$ at $z=1.25\mathrm{m}$, so these are indeed a maximum followed by the next minimum. Also $r/z\simeq 1.5\times 10^{-3}$, consistent with the paraxial approximation.

For the normally incident plane wave, the incident complex amplitude and phase are uniform across the aperture. Treat the field using scalar diffraction, with the aperture represented by its prescribed amplitude transmittance $h(x,y)$. Let $R$ be the distance from the aperture centre to the observation point and let ${\ell }_x,{\ell }_y$ be the transverse direction cosines of that point.

Across an aperture of characteristic size $a$, the Fraunhofer diffraction assumptions are that $a/R\ll 1$, so $1/r\simeq 1/R$ and the obliquity factor is nearly constant, and that the quadratic contribution to the propagation phase is negligible:The remaining path-length dependence is linear:A coherent monochromatic wave and an aperture large enough relative to wavelength for the scalar aperture model complete the usual physical assumptions. The complex amplitude becomeswhere $C$ includes the common propagation phase, distance factor and nearly constant obliquity factor. For the central part of a planar screen in the paraxial approximation, ${\ell }_x\simeq X/R$ and ${\ell }_y\simeq Y/R$. Small observation angles are useful for that screen-coordinate formula; a general far-field Fourier transform can instead use the exact direction cosines. A lens can also produce this pattern in its focal plane without requiring a physically very distant screen.

For a plane wave wavefunction ${\psi }_k(x)\propto e^{ikx}$, impose the cyclic or periodic boundary condition ${\psi }_k(x+L)={\psi }_k(x)$. It requiresAdjacent allowed wavenumbers on the full signed $k$ axis are spaced by $\Delta k=2\pi /L$. Ignoring spin, the density of states per unit $k$ is consequentlyBoth positive and negative $k$ are included in this density. It is not yet the density per unit magnitude $|k|$, which has an extra factor of two away from zero. In $d$ independent dimensions the same periodic boundary conditions give one spatial state per $k$-space cell of measure $(2\pi /L{)}^d$.

For a plane wave propagating along $\hat{\mathbf{k}}$ in vacuum, Maxwell equations give $\mathbf{B}=\hat{\mathbf{k}}\times \mathbf{E}/c$, with $\mathbf{E}\perp \hat{\mathbf{k}}$. HenceTaking $E=E_0\cos (\mathbf{k}\cdot \mathbf{r}-\omega t)$, the period average of the squared cosine is $1/2$. ThusThe fields are in phase, so the travelling-wave energy flux never reverses. The peak amplitude convention is essential; using a root-mean-square electric field would remove the factor $1/2$.

In charge- and current-free vacuum, Maxwell equations reduce toTaking a curl and using $\nabla \times \nabla \times \mathbf{E}=\nabla (\nabla \cdot \mathbf{E})-{\nabla }^2\mathbf{E}$ givesand the same wave equation for $\mathbf{B}$. For a plane wave $\mathbf{E}=\mathrm{Re}[{\mathbf{E}}_0{e}^{i(\mathbf{k}\cdot \mathbf{r}-\omega t)}]$, $\omega =ck$, $\mathbf{k}\cdot {\mathbf{E}}_0=\mathbf{k}\cdot {\mathbf{B}}_0=0$ and ${\mathbf{B}}_0=\hat{\mathbf{k}}\times {\mathbf{E}}_0/c$. Thus the fields are transverse, mutually perpendicular and in phase; polarization describes the electric-field trajectory at a fixed point.

The energy density and Poynting vector areFor a linearly polarized sinusoid with peak electric amplitude $E_0$, the mean intensity is ${ϵ}_{0}c{E}_0^2/2$. Electromagnetic momentum density is $\mathbf{S}/c^2$; normal absorption produces radiation pressure $⟨S⟩/c$, while ideal retroreflection doubles it. Linear superposition permits interference and standing waves, whose local electric and magnetic fields need not have the traveling-plane-wave equality point by point.

In microscopic SI form, Maxwell's equations couple charge density $\rho$ and current density $\mathbf{J}$ to electric and magnetic fields:Gauss's electric law expresses electric flux from charge and reduces to the inverse-square Coulomb field for a point charge. Gauss's magnetic law gives zero net magnetic flux through every closed surface; magnetic field lines do not begin or end on magnetic monopoles in this theory.

Faraday's induction law says changing magnetic flux produces a circulating electric field. For a fixed loop, its circulation is minus the rate of change of enclosed magnetic flux, explaining the sign of induction and Lenz's law. For example, changing a primary coil's electric current changes the flux through a stationary secondary coil and induces a voltage. A moving wire also has the magnetic part of the Lorentz force on its charges; the fixed-loop field equation and the full motional-emf calculation must be distinguished.

The Ampère-Maxwell law relates magnetic circulation to conduction current and displacement current. The latter is essential around a charging capacitor: a surface cutting the lead and another passing through the gap must predict the same circulation. Taking the divergence givesusing Gauss's electric law. Charge conservation therefore fixes the consistency of the time-dependent equations. In matter, the equivalent macroscopic equations use $\mathbf{D}$ and $\mathbf{H}$, free sources and constitutive relations. The Lorentz force density $\rho \mathbf{E}+\mathbf{J}\times \mathbf{B}$ describes exchange of momentum with charges.

In source-free vacuum, take the curl of Faraday's law and use $\nabla \times (\nabla \times \mathbf{E})=\nabla (\nabla \cdot \mathbf{E})-{\nabla }^2\mathbf{E}$. Both fields satisfyFor a plane wave, $\omega =ck$, $\mathbf{k}\cdot \mathbf{E}=\mathbf{k}\cdot \mathbf{B}=0$ and $\mathbf{B}=\hat{\mathbf{k}}\times \mathbf{E}/c$. The fields are transverse, mutually perpendicular and in phase. Their energy density and energy flux areThis is Poynting's theorem: field energy is transported and can do work on matter. Electromagnetic radiation is a propagating, energy-carrying solution of the coupled electric and magnetic field equations.

A small displacement $\mathbf{u}$ in a homogeneous isotropic elastic solid obeys linear elasticity. With Lamé constants $\lambda ,\mu$ and mass density $\rho$, momentum balance givesA longitudinal plane wave has displacement parallel to its propagation direction and changes local volume. A transverse wave has perpendicular displacement and produces shear. Substitution givesAn ideal fluid has no static shear modulus and cannot support this bulk transverse elastic wave. A thin rod with free lateral surfaces instead supports an extensional long-wavelength mode with speed $\sqrt{E/\rho }$; its transverse contraction changes the effective stiffness, so it should not be identified with the unbounded-solid longitudinal mode.

Stress and displacement boundary conditions produce reflection, transmission, mode conversion and guided waves. Viscoelasticity and scattering can attenuate waves; finite geometry can make guided modes dispersive even when the bulk elastic constants are frequency independent. Wave speed depends on both material stiffness and the allowed deformation.

For a homogeneous isotropic solid with small displacement $\mathbf{u}$, linear elasticity gives ${\sigma }_{ij}={\lambda }_L(\nabla \cdot \mathbf{u}){\delta }_{ij}+2\mu e_{ij}$. Momentum balance becomesA plane wave with displacement parallel to its wavevector is longitudinal, whereas perpendicular displacement is transverse. Writing $K={\lambda }_L+2\mu /3$ for the bulk modulus and $\mu$ for the shear modulus, their bulk speeds areThere are two independent transverse polarizations. The bulk modulus supplies resistance to compression; shear rigidity supplies transverse restoring stress. An ordinary fluid has no static shear rigidity and cannot sustain these propagating bulk shear waves. A thin rod with freely contracting sides instead has longitudinal speed $\sqrt{Y/\rho }$, because its lateral boundary conditions differ from a bulk plane wave. Reflection, refraction and mode conversion occur at interfaces; a free surface also supports Rayleigh waves.