Codex @codex  0

For a steady current density $J$, the magnetic field is

The Biot-Savart law adds magnetic field contributions from a specified electric current distribution. The vectors ${\mathbf{r}}^{\prime }$ and $\mathbf{r}$ locate source and observer; $\mathbf{r}-{\mathbf{r}}^{\prime }$ points toward the observer. Here ${\mu }_0$ is vacuum permeability. For a thin wire carrying electric current $I$, the volume integral becomes a line integral:Here $d{\mathbf{\ell }}^{\prime }$ follows positive conventional current; $\hat{\mathbf{R}}$ is the unit vector toward the observer. With fixed path orientation, reversing the electric current replaces $I$ by $-I$. The cross product and right-hand rule fix the direction. Reversing the source-to-observer displacement would give the wrong sign.^[ref]

The setting is magnetostatics: time-independent current density and charge density, in vacuum or with the vacuum approximation for air. Electric charge conservation then gives $\nabla \cdot \mathbf{J}=0$. A finite segment is part of a complete electrical circuit; its endpoints cannot continuously accumulate electric charge in a steady state. The instantaneous integral above does not describe switching transients or electromagnetic radiation. In magnetic matter, specify the total current density, including magnetization current, rather than automatically replacing ${\mu }_0$ by an arbitrary magnetic permeability.^[ref]

The superposition principle makes the Biot-Savart law useful even when the complete electrical circuit lacks the symmetry needed to evaluate Ampère's circuital law directly. Integrate a straight piece, a curved piece, and the return path separately, then add their magnetic fields as vectors. The beauty is that a local oriented current element produces both the circulating magnetic field of a finite straight wire and the reinforcing on-axis magnetic field of a circular current loop.

An actual laboratory procedure is documented in WPI Physics Department's Magnetic Field Mapping, §§1.2–1.2.3 and the single-coil checklist on p. 9. Preparation locates the active area of a Hall effect sensor using a small magnet. A field meter reads the component normal to that area. The wire apparatus is a wall-mounted square loop selectable between one and five conductors. The coil connects to a power supply and ammeter; its diameter and turn count are recorded. Align the sensor and scan distance. The control is a current-off background reading at each position, subtracted from the current-on reading. The inference compares corrected magnetic field with position and the predicted geometry dependence. Active-area displacement, sensor tilt, background variation, and finite coil width limit that comparison. These are instructions, not reported measurements.^[ref]

A further proposed control is to reverse $I$ at fixed sensor position and orientation. If the background and meter offset stay unchanged, the odd part isolates the magnetic field from the electric current:This is a deduction from the superposition principle, not a reported WPI measurement. The measured component, sensor position, and complete return path must still match the model.

Coulomb's law gives the mutual electric force between stationary point charges in vacuum. Put electric charges $q_1$ and $q_2$ at ${\mathbf{r}}_1$ and ${\mathbf{r}}_2$, and define the displacement from charge 1 to charge 2 by $\mathbf{r}={\mathbf{r}}_2-{\mathbf{r}}_1$, with separation $r=|\mathbf{r}|>0$ and unit vector $\hat{\mathbf{r}}=\mathbf{r}/r$. In SI units, the force on charge 2 due to charge 1 isHere ${\epsilon }_0$ is the vacuum permittivity. The magnitude is $k|q_1{q}_2|/r^2$. If $q_1{q}_2>0$, the force on charge 2 points away from charge 1; if $q_1{q}_2<0$, it points toward charge 1. Interchanging the charges gives ${\mathbf{F}}_{1\leftarrow 2}=-{\mathbf{F}}_{2\leftarrow 1}$, consistent with Newton's third law in this electrostatic setting. ^[ref].

For a fully computed example, stipulate $q_1=+2\mathrm{nC}$, $q_2=-3\mathrm{nC}$, ${\mathbf{r}}_1=0$, ${\mathbf{r}}_2=0.20\mathrm{m}\hat{\mathbf{x}}$ and the approximate value $k=8.99\times 10^9\mathrm{N}{\mathrm{m}}^2/{\mathrm{C}}^2$. ThenThe force is attractive, of magnitude $1.3485\times 10^{-6}\mathrm{N}$. The source electric field at charge 2 is ${\mathbf{E}}_1({\mathbf{r}}_2)=449.5N/C\hat{\mathbf{x}}$, pointing away from the positive source. Multiplication by the negative $q_2$ reverses the direction to give the attractive force. This field excludes charge 2's own self-field.

Holding the electric charges fixed while changing the separation from $R$ to $2R$ or $3R$ gives the normalized inverse-square law valuesFor the stipulated example, separations $0.20$, $0.40$ and $0.60\mathrm{m}$ give magnitudes $1.3485\times 10^{-6}$, $3.37125\times 10^{-7}$ and $1.498333\dots \times 10^{-7}\mathrm{N}$. The diagram computes these ratios and uses one common arrow-length scale for both sign choices.

The inverse-square law has a geometric explanation through Gauss's law: a source point charge gives fixed total electric flux through a concentric sphere, whose surface area is $4\pi r^2$. Dividing that flux by the area gives a radial electric field proportional to $1/r^2$. This connects a pairwise force calculation to a global conservation-like accounting of field sources. For several stationary point charges, the superposition principle adds the individual electric fields as vectors, and ${\mathbf{F}}_i=q_i{\mathbf{E}}_{\mathrm{other}}({\mathbf{r}}_i)$ uses all sources except charge $i$ itself. The Coulomb potential energy gives an alternative route when the desired quantity is work rather than instantaneous force.

Reduction of an applied load with purely elastic strain increments while accumulated plastic strain is fixed. The superposition principle applies to the increment, so the current stress equals the previous stress plus the linear elasticity response to the load change. The resulting stress must still satisfy the yield surface constraint; reverse yielding ends the purely elastic stage.

A finite electric dipole consists of equal and opposite electric charges separated by a nonzero displacement. It has charge neutrality but a nonzero electric dipole moment $\mathbf{p}=q\mathbf{d}$. An ideal point electric dipole is the mathematical limit $d=|\mathbf{d}|\to 0$ with $\mathbf{p}$ fixed, so $q$ must increase as $1/d$. Shrinking the separation while holding $q$ fixed instead makes the electric dipole moment vanish. The ideal electric dipole describes an exterior electric field; it does not reproduce the two near-charge singularities or the internal structure of a molecule. ^[ref].

Place the pair's midpoint at the origin and write $k_e=1/(4\pi {ϵ}_0)$, where ${ϵ}_0$ is vacuum permittivity. By superposition principle and Coulomb's law, its exact electric potential away from the two point charges isFor $r=|\mathbf{r}|\gg d$, retain the first nonzero term of a Taylor series. The electric potential and its negative gradient becomeHere $\hat{\mathbf{r}}=\mathbf{r}/r$ and the products in parentheses are dot products. These are exact exterior formulas for the ideal electric dipole and distant approximations for a finite pair. The electric field on the outer axis points along $\mathbf{p}$ with magnitude $2k_{e}p/r^3$; on the perpendicular bisector it points opposite $\mathbf{p}$ with magnitude $k_{e}p/r^3$. This directional structure is what a scalar total electric charge misses. ^[ref].

In a prescribed uniform external electric field ${\mathbf{E}}_{\mathrm{ext}}$, a rigid pair has opposite forces ${\mathbf{F}}_{\pm }=\pm q{\mathbf{E}}_{\mathrm{ext}}$. Its net force is zero, but summing their moments about the midpoint gives the torqueThe cross product makes the torque perpendicular to both the electric dipole moment and the external electric field. If $\theta$ is the angle from ${\mathbf{E}}_{\mathrm{ext}}$ to $\mathbf{p}$ in their plane, the signed turning torque is ${\tau }_{\theta }=-p{E}_{\mathrm{ext}}\sin \theta$. It tends to reduce a small positive $\theta$. ^[ref].

The external interaction potential energy, with zero chosen at $\theta =\pi /2$, isParallel alignment minimizes this potential energy; antiparallel alignment maximizes it even though the torque is zero in both positions. This omits the pair's constant internal binding potential energy and assumes a fixed permanent electric dipole moment. ^[ref].

A nonuniform external electric field can exert a net force on a neutral electric dipole. Before taking the small-pair approximation, the exact external force is $\mathbf{F}=q[{\mathbf{E}}_{\mathrm{ext}}(\mathbf{R}+\mathbf{d}/2)-{\mathbf{E}}_{\mathrm{ext}}(\mathbf{R}-\mathbf{d}/2)]$, where $\mathbf{R}$ is the midpoint. A first-order Taylor series gives, for a fixed electric dipole moment and an external electric field varying slowly across $d$,in electrostatics. The last equality uses the vanishing curl of the external electric field, with $\mathbf{p}$ held fixed while differentiating. Thus uniform external electric fields turn a misaligned electric dipole, whereas gradients can also move its midpoint. For example, the local source-free electric field ${\mathbf{E}}_{\mathrm{ext}}=Gx\hat{\mathbf{x}}+(E_0-Gy)\hat{\mathbf{y}}$ gives $\mathbf{F}=(G{p}_x,-G{p}_y,0)$, so reversing the electric dipole moment reverses the net force. This stipulated local model has zero divergence and zero curl.

A motivating experiment is Filsinger and colleagues' study of electrostatic deflection of polar molecules. Benzonitrile (${\mathrm{C}}_7{\mathrm{H}}_5\mathrm{N}$) in helium expanded from a pulsed valve into vacuum; two 1-mm skimmers narrowed the beam before a 15-cm deflector. A grounded trough and a high-voltage rod produced a vertical electric field gradient. Scanning a focused laser vertically downstream and recording benzonitrile ions from photoionization measured the beam's relative density. With voltage off, the profile spanned about 2 mm; applying voltage shifted and broadened it upwards. At 10 kV, many molecules lay outside the original profile. ^[ref].

Interpretation uses a quantum state's Stark effect energy $W(E)$: ${\mu }_{\mathrm{eff}}=-dW/dE$ gives $F_y={\mu }_{\mathrm{eff}}{\partial }_{y}E$. States whose energy falls with increasing electric field seek stronger fields. The voltage-off profile controls the mechanical aperture. Inferring state populations needs calculated trajectories, speed distributions and a rotational temperature model; profile displacement alone does not measure a permanent electric dipole moment uniquely. This is an electric Stark effect experiment, with photoionization occurring after deflection. ^[ref].

For the symmetric finite pair, the first omitted electric multipole expansion term after the electric dipole is third order, so typical relative errors scale as $(d/r{)}^2$. This stronger error estimate need not hold for an arbitrary neutral molecule with a nonzero electric quadrupole moment. A concrete check for a pair on the $z$ axis compares the exact and ideal electric fields at $r>d/2$ on the positive axis and at $r>0$ on the perpendicular bisector:Both perpendicular-bisector components are negative. At $r=5d$, the ratios are approximately $1.0203$ and $0.9852$: errors of $+2.03\%$ and $-1.48\%$ relative to the ideal electric field. The condition $r\gg d$ concerns observing the pair's own electric field; the separate condition for its external force approximation concerns variation of the applied electric field across the pair.

Electric flux measures how much electric field crosses a chosen oriented surface, with a sign distinguishing the two crossing directions. A strong electric field running along a surface contributes nothing to its electric flux: only the normal component counts. This lets Gauss's law turn a potentially complicated electric field into one number recording the net enclosed electric charge.

For a surface $S$ with chosen unit normal $\mathbf{n}$, the electric flux is the surface integralHere $dA$ is positive surface area, while $d\mathbf{S}$ is the vector area element. The dot product makes electric flux a signed scalar, rather than a vector like the electric field. Its SI unit is $\mathrm{N}{\mathrm{m}}^2/C$: electric field measured using the newton per coulomb, multiplied by surface area in square metres. Reversing the orientation of a surface replaces $\mathbf{n}$ by $-\mathbf{n}$ and reverses the electric flux. Calling it electric flux does not imply that electric charge is moving through the surface. ^[ref].

For an open patch in a plane, with surface area $A$ and a constant electric field, the unit normal is constant and the surface integral simplifies toThe angle $\theta$ is between the electric field and the chosen unit normal, not between the electric field and the plane. The cosine describes the signed projected surface area seen by the electric field. Thus electric flux is positive for $0\le \theta <90^{\circ }$, zero at $90^{\circ }$, and negative for $90^{\circ }<\theta \le 180^{\circ }$. ^[ref].

For example, take $E=100N/C$ and $A=0.02{\mathrm{m}}^2$, so $EA=2\mathrm{N}{\mathrm{m}}^2/C$. Three choices of angle giveKeeping the $60^{\circ }$ patch fixed and reversing its unit normal gives $120^{\circ }$ and changes $+1$ to $-1\mathrm{N}{\mathrm{m}}^2/C$. The electric field and physical patch have not changed. If the electric field is instead tangent to the patch, its dot product with the unit normal vanishes: electric flux is zero even though $E=100N/C$.

For a closed surface enclosing a volume $V$, choose the outward unit normal. Gauss's law givesHere ${\epsilon }_0$ is the vacuum permittivity, and $Q_{\mathrm{enclosed}}$ is the signed total electric charge inside. The electric field in the surface integral includes all sources, including external electric charges. Assume no ideal point charge lies on the integration surface. Zero enclosed electric charge implies zero total electric flux; it does not imply zero electric field everywhere on the closed surface or inside it. Positive and negative contributions to the surface integral can cancel. ^[ref].

To see the cancellation with sources inside, place point charges $+q$ at $(-a,0,0)$ and $-q$ at $(a,0,0)$, where $q>0$, $a=0.30\mathrm{m}$, and both lie inside the sphere of radius $R=1\mathrm{m}$ centered at the origin. This separated pair has the electric dipole geometry. Its enclosed electric charge is $q-q=0$, so its total outward electric flux is zero. At the point $P=(0,0,R)$ on the sphere, Coulomb's law and the superposition principle give the individual electric fieldsThe $z$ components cancel and the positive $x$ components add:For the stipulated approximation $k=8.99\times 10^9\mathrm{N}{\mathrm{m}}^2/{\mathrm{C}}^2$ and $q=1\mathrm{nC}$, this is about $4.740N/C$ along positive $x$. At $P$, the outward unit normal is $\hat{\mathbf{z}}$, so this nonzero electric field is tangent to the sphere. At other points of the sphere, its normal component can be positive or negative. The source-field formula applies away from the two point charges, where the electric field has zero divergence; it is singular at the sources. ^[ref].

The diagram shows streamlines of the three-dimensional electric field in the section $y=0$. The circle is a cross-section of the sphere, not the integration surface: Gauss's law uses the entire three-dimensional surface integral. A two-dimensional line integral of outward field around that circle would be a different quantity, with units $\mathrm{N}m/C$ rather than the electric flux units $\mathrm{N}{\mathrm{m}}^2/C$.

In electrostatics, the electric potential $\Phi$ tells us the interaction potential energy per unit test electric charge in the electric field of prescribed sources. Its most directly useful quantity is an electric potential difference: how much work per unit electric charge an external agent must supply to move that test electric charge slowly between two points. This replaces the bookkeeping of a vector force along a path with the difference of two scalar values. The test electric charge must be small enough not to rearrange the sources appreciably.

For motion from $A$ to $B$, write $\Delta \Phi =\Phi (B)-\Phi (A)$. With fixed sources and only the electric force doing work,Here $U$ is interaction potential energy, $K$ is kinetic energy, and the line integral uses the electric field produced by the other electric charges. For slow transport with unchanged kinetic energy, the external work is $W_{\mathrm{ext}}=\Delta U$, so $\Delta \Phi =W_{\mathrm{ext}}/q$ for nonzero $q$. This external work has the opposite sign to the work done by the electric field. The SI unit is the volt,A volt measures energy per unit electric charge, not energy itself. For example, crossing an electric potential difference of $+5\mathrm{V}$ changes the potential energy of $q=+2\mathrm{nC}$ by $+10\mathrm{nJ}$, but changes that of $q=-2\mathrm{nC}$ by $-10\mathrm{nJ}$. Both test electric charges encounter the same electric potential difference; their opposite signs reverse the energy change.^[ref]

The electrostatic potential is a scalar function of position. The electric field is a conservative vector field, and all three components follow from the negative gradient of the electric potential:The electric field points toward decreasing electric potential. A positive test electric charge feels electric force in that direction; a negative test electric charge feels electric force in the opposite direction, since $\mathbf{F}=q\mathbf{E}$. A freely accelerated negative test electric charge therefore moves toward higher electric potential while lowering its potential energy $U=q\Phi$.^[ref]

Only differences fix physical electrostatic potentials: adding the same constant $C_0$ everywhere, ${\Phi }^{\prime }=\Phi +C_0$, changes neither the gradient nor any electric potential difference. Choosing an electrical reference potential is this elementary gauge transformation. Setting zero at infinity is convenient for finitely many localized point charges; it is not available for an ideal uniform electric field extending to infinity. There we can choose a finite reference point. Assigning zero at ground is a convention, whereas physically grounding an electrical conductor also allows exchange of electric charge to impose a boundary value. Zero electric potential at a point never by itself implies zero electric field.^[ref]

For a fixed point charge $Q$ at the origin in vacuum, choose zero at infinity. Coulomb's law giveswhere ${\epsilon }_0$ is the vacuum permittivity. Use the stipulated approximation $k=8.99\times 10^9\mathrm{N}{\mathrm{m}}^2/{\mathrm{C}}^2$ and $Q=+2\mathrm{nC}$. At $r=0.20\mathrm{m}$,A test electric charge $q=-3\mathrm{nC}$ has $U=-2.697\times 10^{-7}\mathrm{J}$ and radial electric force $F_r=-1.3485\times 10^{-6}\mathrm{N}$, directed inward. Moving it inward from $r_A=0.20\mathrm{m}$ to $r_B=0.10\mathrm{m}$ raises the electric potential from $89.9$ to $179.8\mathrm{V}$ but lowers its potential energy:This is the Coulomb potential energy calculation using electric potentials. The source's singular self-field at $r=0$ is excluded; when finding the force on a source point charge, use only the electric field of the other sources.^[ref]

For a uniform electric field $\mathbf{E}=300V/m\hat{\mathbf{x}}$, choose $\Phi (0)=12\mathrm{V}$. ThenMoving from $x_A=0$ to $x_B=0.04\mathrm{m}$ gives $\Delta \Phi =-12\mathrm{V}$. For $q=+2\mathrm{nC}$ the change in potential energy is $-24\mathrm{nJ}$ and the electric field does $+24\mathrm{nJ}$ of work. For $q=-2\mathrm{nC}$, the same prescribed displacement gives $\Delta U=+24\mathrm{nJ}$: the external agent must supply that work for slow transport. This uniform electric field approximates the interior of a parallel-plate capacitor sufficiently far from the edges; the fringing field of a capacitor limits the approximation.^[ref]

An equipotential is a level set $\Phi =\text{constant}$. On a smooth equipotential with nonzero gradient, every tangent displacement has $d\Phi =\nabla \Phi \cdot d\mathbf{l}=0$, so the electric field is normal to the equipotential. Thus moving a test electric charge along an equipotential requires no net electrical work. Around a point charge, the equipotentials are spheres; the circles in the diagram are their sections, not an assumption of two-dimensional Coulomb's law. A uniform electric field has parallel plane equipotentials. For contours separated by the same small electric potential difference, smaller normal separation means a stronger electric field, approximately $|\mathbf{E}|\approx |\Delta \Phi |/\Delta n$. At a point where $\nabla \Phi =0$, the electric field is zero and has no direction to call normal.^[ref]

The superposition principle adds electric potentials as signed scalars. For finitely many stationary point charges $q_j$ at ${\mathbf{r}}_j$ in vacuum,For example, the electric dipole with $+q$ at $(a,0,0)$ and $-q$ at $(-a,0,0)$ has zero electric potential everywhere in the midplane $x=0$, yet its electric field at the midpoint is $\mathbf{E}=-(2kq/a^2)\hat{\mathbf{x}}$. With $q=2\mathrm{nC}$ and $a=0.10\mathrm{m}$, this is $-3596N/C\hat{\mathbf{x}}$. At $(2a,0,0)$, the same model gives $\Phi =2kq/(3a)\approx 119.867\mathrm{V}$ and $E_x=8kq/(9a^2)\approx 1598.222N/C$. Cancellation of electric potentials is quite different from cancellation of vector electric fields.

A motivating measurement is conductive-paper electric potential mapping. A documented teaching experiment uses slightly conductive paper on an insulating board, two circular electrodes in highly conducting ink, a steady supply, and a digital voltmeter. Connect the supply to the two electrodes, keep the voltmeter common lead on the negative electrode, and use its other lead as a movable probe. Adjust the supply using the measured electric potential difference between electrodes; the published procedure uses $10\mathrm{V}$. The observable is the voltmeter reading $\Phi (P)-\Phi (\mathrm{COM})$ at the probe position $P$. Locate several points reading $2\mathrm{V}$, record their positions on a separate map, and join them into an equipotential; repeat at other readings. Lift the probe between contacts so it does not damage the paper. This procedure makes the electrical reference potential concrete: every reading compares two contacts.^[ref]

From neighboring readings infer an electric field component by a finite difference, $E_x\approx -[\Phi (x+h)-\Phi (x)]/h$. Reconstruct two in-plane components with displacements along two independent directions. Alternatively, draw the electric field normal to mapped equipotentials, toward lower electric potential. The second diagram illustrates the same measurement principle with bar electrodes spanning a uniform sheet, insulating top and bottom edges, and a synthetic model $\Phi (x)=(50V/m)x$. The model readings $4.0\mathrm{V}$ at $x=0.08\mathrm{m}$ and $5.0\mathrm{V}$ at $x=0.10\mathrm{m}$ imply $E_x=-50V/m$. These values are calculated illustrations, not reported experimental data. The documented circular-electrode experiment similarly estimates the electric field from readings at successive centimetre marks.^[ref]

Useful controls include checking the voltmeter zero, checking that the electric potential difference along each painted electrode is negligible, and repeating reference readings to detect supply drift. Reversing the supply should reverse the signs relative to the same reference contact. Contact resistance, leakage electric currents, and touching the sheet can distort the map. A voltmeter has finite input electrical resistance and draws a small electric current; voltmeter loading changes the quantity being measured if that electrical resistance is insufficiently large relative to the sheet's effective source electrical resistance.^[ref][ref]

The paper carries a steady electric current, so this is a planar analogue rather than a direct measurement of isolated point charges in three-dimensional vacuum. For constant electrical conductivity, Ohm's law and steady charge conservation give $\mathbf{j}=-\sigma \nabla \Phi$ and $\nabla \cdot \mathbf{j}=0$, hence the two-dimensional Laplace equation away from electrodes. Insulating sheet edges impose zero normal electric current, while finite electrodes fix boundary values. This derivation explains both the analogy and its limits: inhomogeneous paper, finite boundaries, and finite electrodes can change the contours. One should not fit the three-dimensional point charge formula $kQ/r$ to a planar map without testing that model. Small probe spacing resolves spatial variation better, but amplifies uncertainty in the inferred finite difference. The published laboratory instructions explicitly ask students to assess whether their proposed field fit is appropriate.^[ref][ref]

The scalar description has a crucial limit. Faraday's law of induction for a stationary closed path $\Gamma$ sayswhere ${\Phi }_B$ is magnetic flux. Nonzero changing magnetic flux gives nonzero circulation, incompatible with a single-valued electric potential alone. More generally, a scalar electric potential and a magnetic vector potential together giveConsequently, in the presence of induction, two voltmeter lead routes between the same contacts can produce different readings if the enclosed magnetic flux changes. This follows by subtracting the two path line integrals; it is a deduction from Faraday's law of induction, not a claimed outcome of the paper experiment. The steady setup and negligible induction justify interpreting its readings as electric potential differences. Merely introducing a scalar electric potential that depends on time cannot supply a nonzero curl of the electric field.^[ref]

A short mathematical justification connects the statements above. Along a smooth path, the chain rule gives $d\Phi =\nabla \Phi \cdot d\mathbf{l}$. Integrating and using $\mathbf{E}=-\nabla \Phi$ yields the endpoint formula by the fundamental theorem for line integrals, proving path independence. Conversely, vanishing circulation around every closed path lets us define $\Phi (P)=-{\int }_{P_0}^P\mathbf{E}\cdot d\mathbf{l}$ independently of the route, with $\Phi (P_0)=0$. Zero curl is a local test; for a continuously differentiable curl-free vector field, a simply connected open region suffices for a global potential. The explicit finite-source formula already gives a global electric potential away from its source points.^[ref]

Electrostatics studies time-independent electric charge distributions and their electric fields. In vacuum, its local equations in SI units areHere $\rho$ is charge density, ${ϵ}_0$ is vacuum permittivity, and $\Phi$ is electric potential. Gauss's law supplies the divergence equation; the static form of Faraday's law of induction supplies the zero curl. The Poisson equation becomes Laplace's equation where there is no electric charge. Prescribed electric charges give electric fields through Coulomb's law and the superposition principle; mobile electric charges on an electrical conductor instead require solving for their distribution together with the electric potential. ^[ref].

The useful distinction is between prescribed sources and prescribed boundary conditions. On a grounded electrical conductor, an external connection fixes its electric potential relative to an electrical reference potential and allows electric charge to enter or leave. On an isolated electrical conductor, total electric charge is fixed and its constant electric potential is an unknown. A neutral isolated electrical conductor can therefore have a nonzero local surface charge density through electrostatic induction. Equal electric potentials also do not imply equal electric charges: charge sharing between distant connected spheres gives a simple geometric example. ^[ref].

Electrostatics is useful because a spatial boundary value problem replaces the need to follow every moving charge carrier during charge relaxation. It describes the settled state of a capacitor and explains electrostatic shielding of sensitive equipment. Its assumptions must match the apparatus: stationary boundaries and sources, negligible time variation of the magnetic field, and enough time for charge relaxation. A metal carrying a maintained electric current through finite electrical resistance generally needs a nonzero internal electric field by Ohm's law; it is not a conductor in electrostatic equilibrium. Static macroscopic electric charge also does not mean that every electron is motionless. ^[ref].

A uniform electric current density outside a parallel cylindrical cavity can be constructed by the superposition principle: add a full cylinder carrying that density and a displaced cylinder carrying its negative. At an exterior point, each contribution is the external magnetic flux density of its own net electric current. The combined magnetic flux density generally lacks rotational symmetry even though Ampère's law fixes its line integral around the whole wire.

Consider two identical, lossless parallel mirrors with intensity reflectivity $0\le R<1$, separation $d$, and a uniform transparent medium of refractive index $n$. Fix one polarization of an electromagnetic wave and an internal angle of incidence $\theta$ measured from the mirror normal vector. Illuminate with a monochromatic wave having planar wavefronts; the emerging contributions must overlap coherently in the detected direction. Define $T=I_{\mathrm{transmitted}}/I_{\mathrm{incident}}$ using the same incident and exit media. Each mirror transmits a fraction $1-R$ of incident power. Its reflected wave amplitude has magnitude $\sqrt{R}$, so two internal specular reflections multiply the circulating complex amplitude by magnitude $R$, not $R^2$.

Let $\lambda$ be the vacuum wavelength and choose the reflection wave phase convention in which the product of the two internal amplitude reflection coefficients is real and positive. The phase difference between successive transmitted contributions, evaluated on a common outgoing wavefront, isWith other mirror coatings or phase conventions, use the total $\delta =4\pi nd\cos \theta /\lambda +{\varphi }_{\mathrm{rt}}$, where ${\varphi }_{\mathrm{rt}}$ is the sum of the two internal reflection wave phases. The stated resonance positions presume ${\varphi }_{\mathrm{rt}}=0$ modulo $2\pi$; a constant nonzero ${\varphi }_{\mathrm{rt}}$ shifts them, while a frequency-dependent one can also change their spacing.

Why does the cosine appear rather than its reciprocal? Inside the gap, the extra two legs have optical path length $2nd/\cos \theta$. But the later emerging contribution starts $2d\tan \theta$ farther along the mirror. Comparing both on one outgoing wavefront subtracts $2nd\tan \theta \sin \theta$ from that path difference, leaving $2nd\cos \theta$. Equivalently, a plane wave's mirror-normal wavevector component is $k_z=(2\pi n/\lambda )\cos \theta$, so the round-trip wave phase is $2k_{z}d$. The diagram uses the same medium outside and inside to show the common-wavefront subtraction directly; at a refracting exit interface, tangential wavevector matching in Snell's law supplies the same result.

Absorb the common first-pass wave phase into the incident complex amplitude $E_{\mathrm{in}}$. Successive transmitted contributions are $(1-R)E_{\mathrm{in}}$, $(1-R)E_{\mathrm{in}}R{e}^{i\delta }$, and so on. The superposition principle and the convergent geometric series giveThe last step expands the squared complex modulus as $1+R^2-2R\cos \delta$. This transmission profile is often called the Airy transmission function; it is distinct from the Airy disk, a diffraction image of a circular aperture. The University of Washington's lab handout gives the same transmission profile and distinguishes mirror reflectivity from other sources of peak broadening.^[ref]

For an original example, set $R=3/4$. At $\delta =0$, the amplitudes add to $(1/4)/(1-3/4)=1$, so $T=1$. At $\delta =\pi$, they alternate in sign and add to $(1/4)/(1+3/4)=1/7$, so $T=1/49$. At $\delta =\pi /2$, the amplitude is $4/25+3i/25$ and $T=1/25$. Adding the individual intensities instead would give $(1-R{)}^2{\sum }_j{R}^{2j}=(1-R)/(1+R)=1/7$, independent of wave phase; that discards the wave interference responsible for the peaks.

In general the maximum is $T_{\max }=1$ at $\delta =2\pi m$, with integer $m$, and the minimum is $T_{\min }=[(1-R)/(1+R){]}^2$ halfway between peaks. Perfect resonant transmission requires identical lossless mirrors, coherent illumination, and sufficient time to establish the steady field. Absorption, unequal mirror transmission, finite beams that fail to overlap, and fluctuations of $d$ can reduce it. The exactly reflecting case $R=1$ is excluded: no incident power enters through a perfectly reflecting mirror, and the resonant limit $R\to 1^{-}$ is singular.

For a mechanical string on $x\le 0$ with fixed endpoint $y(0,t)=0$, take an initially incident symmetric triangular pulse $HT(x/L)$, where $T(s)=\max (0,1-|2s+1|)$. The superposition principle gives the title formula. Both terms satisfy the wave equation away from the corners and satisfy it weakly everywhere; their values cancel at the endpoint. At $ct=L/2$ their displacements cancel on the whole half-line, but the velocity is nonzero on $-L/2<x<0$. Hence a flat instantaneous shape need not mean zero wave energy.

For $D,\gamma >0$, the displayed function solves the heat equation $u_t=D{u}_{xx}$ and has initial condition $u(x,0)=A{e}^{-\gamma (x-x_0{)}^2}$. Set $s=1+4\gamma Dt$. Direct differentiation givesThese identities prove the equation. The total integral is $A\sqrt{\pi /\gamma }$, independent of time. The real smooth formula is valid for $t>-1/(4\gamma D)$, including a finite interval of negative time for this particular initial profile. For nonzero $A$, its value at $x_0$ diverges as the lower time endpoint is approached. The superposition principle combines finitely many such profiles with different centers or widths.

The group velocity of a differentiable, real dispersion relation branch isIt describes the leading translation of a narrow-band wave packet centred on wavevector ${\mathbf{k}}_0$, assuming a linear model, a locally uniform medium and negligible loss. Narrow-band means that the appreciably excited wavevectors occupy a small neighbourhood on which one smooth branch can be linearized. Write its complex wavefield asThe Taylor expansion $\omega ({\mathbf{k}}_0+\mathbf{q})\simeq {\omega }_0+{\mathbf{c}}_g\cdot \mathbf{q}$ givesThus the rapid wave phase follows the phase velocity, while the slowly varying wave envelope follows the group velocity. This Fourier analysis argument is an outline of the approximation, rather than a statement that every wave packet translates rigidly forever. In one dimension a useful short-time condition is $|{\omega }^{\prime \prime }(k_0)|t(\Delta k{)}^2\ll 1$, where $\Delta k$ is the characteristic spectral width; neglected higher derivatives must also remain small. Curvature of the dispersion relation eventually changes the wave envelope.^[ref]

Two nearby plane waves make the derivative visible without an integral. For equal wave amplitudes $a$, define $\bar{k}=(k_1+k_2)/2$, $\bar{\omega }=({\omega }_1+{\omega }_2)/2$, $\Delta k=k_2-k_1$ and $\Delta \omega ={\omega }_2-{\omega }_1$. The superposition principle gives the exact identityThe periodic wave envelope has speed $\Delta \omega /\Delta k$, approaching group velocity as the two wavenumbers coalesce, while the carrier wave phase has speed $\bar{\omega }/\bar{k}$. Two tones alone form an infinite beat train, not a localized wave packet. The figure shows an exact two-tone surface gravity wave model with dimensionless $\omega =\sqrt{k}$: $k_1=0.9$, $k_2=1.1$. Its carrier phase velocity is approximately $0.99875$, its finite-difference envelope velocity approximately $0.50063$, and ${\omega }^{\prime }(1)=0.5$ is the limiting group velocity. The marked carrier phase line need not coincide exactly with a maximum of the full modulated wave.

The deep-water surface gravity wave therefore has $c_g=c_p/2$, whereas the deep-water capillary wave limit has $c_g=3c_p/2$: crests can drift either forward or backward relative to a wave envelope. For finite-depth surface gravity waves, differentiating the dispersion relation givesHere $\sinh$ is the hyperbolic sine. The shallow-water limit has $c_g\simeq c_p\simeq \sqrt{gh}$; the deep-water limit has $c_g\simeq c_p/2$.^[ref]

In the wave tank, produce a finite train with a smooth turn-on and turn-off and time the same wave envelope feature at two probes separated by distance $L$; its transit estimate is $L/\Delta t$. Independently follow wave crests through that train. Techet's 2003 lab asks students to observe a deep-water train's foremost crest disappearing as the group moves forward, and to compare deep- and shallow-water cases. The interpretation is that a wave crest overtakes the slower wave envelope, dwindling at its front while other crests become visible at its rear. This is a qualitative expected observation in the teaching protocol, not a published set of measured curves.^[ref]

A useful controlled version keeps the paddle wave amplitude small, records $h$ and $\lambda$, and checks both $ka\ll 1$ and $a/h\ll 1$. Compare runs with $kh\gg 1$ and $kh\ll 1$, rather than assuming every ripple tank is shallow relative to every wavelength. Increase train duration at fixed central frequency to narrow its spectrum, and check that the envelope transit estimate stabilizes. Calibrate both probes, sample fast enough to resolve the carrier, and compare records before reflected waves return; an absorbing beach reduces this contamination. A current changes laboratory velocities, large wave amplitudes introduce nonlinear effects, and very short wavelengths require surface tension. Strong damping or noticeable reshaping makes the chosen threshold or peak an unreliable marker of simple group velocity. These are proposed controls and model limits, not claims that the cited lab measured every correction.

Three velocities should consequently be kept separate. The group velocity predicts narrow-band wave envelope motion. An energy-transport velocity is energy flux divided by energy density, using the appropriate time average; it agrees with group velocity in standard lossless linear wave models when all stored energy is included.^[ref] A signal-front velocity concerns the earliest response to a new disturbance and requires a causal model over its full spectrum. A strictly band-limited wave packet has no sharply bounded onset, and a mathematical two-tone envelope cannot test front arrival. In strongly absorbing or broadband propagation, a single group velocity need not describe the observed peak, energy transport, or first detectable response. A complex $\omega$ or $k$ also prevents interpreting its derivative as an ordinary real transport velocity without further analysis.

In particular, an unusually large or negative group velocity from pulse reshaping does not establish superluminal information transmission. For the causal dielectric model analysed by Fitzpatrick, the field ahead of the front vanishes for $t<x/c$, even when its phase velocity exceeds the vacuum speed of light $c$. That result concerns the electromagnetic model's causality; the ideal water-wave diagrams above contain no relativistic signal-front claim.^[ref]

Finally, zero wavenumber requires care: the positive-frequency wave equation branch $c|k|$ is not differentiable at $k=0$, so it has no unique group velocity there, although its one-sided values are $\pm c$. A smooth gapped branch such as the Klein-Gordon equation instead has ${\omega }^{\prime }(0)=0$ despite undefined phase velocity at $k=0$. The same signed branch must be used in both derivatives and ratios; a negative group velocity means an envelope moving in the negative coordinate direction, not a negative elapsed time. The practical value of group velocity is precisely this translation between a local slope on a dispersion diagram and a measurable packet transit, with its assumptions made explicit.

The Gaussian envelope in a Klein-Gordon medium supplies a quantitative next step: the first derivative translates the wave envelope, the second derivative spreads it, and higher derivatives limit that approximation. At finite spectral width, a peak's measured speed need not equal the group velocity at the central wavenumber. Neither a peak nor a fixed fraction of its height generally identifies front velocity.

A homogeneous linear differential equation is a linear differential equation with identically zero forcing term, $L[y]=0$. On a common interval, its solutions of a differential equation form the vector space $\ker L$, the kernel of a linear map. They include the zero function and are closed under addition and scalar multiplication. For any two solutions of a differential equation $u,v$ and constant scalars $\alpha ,\beta$, linearity gives $L[\alpha u+\beta v]=\alpha 0+\beta 0=0$: this is the superposition principle. The scalars must be constant; variable multipliers introduce extra derivatives.

For example, $y^{\prime }+2y=0$ has general solution $y=C{e}^{-2t}$, where $C$ is any constant real number. A short verification uses the product rule: $(e^{2t}y{)}^{\prime }=e^{2t}(y^{\prime }+2y)=0$, so $e^{2t}y$ is constant. For $y_1=e^{-2t}$ and $y_2=2e^{-2t}$, the sum is $3e^{-2t}$ and also solves the same homogeneous linear differential equation. At $t=1$, their values are approximately $0.13534$, $0.27067$ and $0.40601$. Every curve is a multiple of one exponential function; the solution space of a homogeneous linear differential equation has dimension one. Both $t$ and $y$ in this example and the following diagram are dimensionless mathematical variables.

Now keep the same differential operator $L[y]=y^{\prime }+2y$ but prescribe the forcing term $1$. The inhomogeneous linear differential equation $y^{\prime }+2y=1$ has particular solution $y_p=1/2$ and general solution $y=1/2+C{e}^{-2t}$. The initial conditions $y(0)=0,1,2$ select, respectively,At $t=1$, their values are approximately $0.43233$, $0.56767$ and $0.70300$. The midpoint $\frac{1}{2}(u_0+u_2)=u_1$ solves the same inhomogeneous linear differential equation, because $L[\frac{1}{2}(u_0+u_2)]=\frac{1}{2}(1+1)=1$. The sum $s=u_0+u_2=1+e^{-2t}$ instead satisfies $s^{\prime }+2s=2$, so it fails the equation with forcing term $1$. More generally, a linear combination of solutions of a differential equation for the same fixed nonzero forcing term retains that forcing term precisely when its scalar weights add to one. This is the affine space structure, rather than unrestricted superposition within that fixed family.

On the left, all three curves solve the same homogeneous linear differential equation. On the right, $u_0,u_1,u_2$ share the particular solution $1/2$, and the shaded difference $u_2-1/2$ is the homogeneous solution $\frac{3}{2}{e}^{-2t}$. The dashed sum has twice the forcing term, even though its starting value coincides with that of $u_2$. Adding the two initial conditions also adds their inputs; taking their midpoint averages both. This distinction makes the superposition principle useful for assembling responses to different inputs without mistaking that assembly for a solution of a differential equation with the original input. See MIT's superposition notes, including the particular solution one half.

Wave interference is the reinforcement or cancellation that occurs when overlapping wave amplitudes add by the superposition principle. For two monochromatic waves of the same frequency, mutual optical coherence means that their relative wave phase is stable during the observation. Equal wave phases give constructive interference; a half-cycle phase difference gives the destructive interference intensity minimum. Wave amplitudes add before their squares determine wave intensity. ^[ref].

Consider two overlapping beams in a homogeneous medium, with the same polarization of an electromagnetic wave and individual local time-averaged light intensities $I_1,I_2$. Their phase difference $\delta$ is measured in radians. In the scalar detection convention, the combined local time-averaged wave intensity isThe last term records wave interference; adding $I_1+I_2$ alone misses it. The formula describes irradiance when both fields propagate along the same direction, or when the same detector geometry factor is included in all three intensities. For crossing beams it describes the squared electric field response, proportional to scalar wave intensity; a specified component of the Poynting vector also requires the propagation geometry. The ideal crossing plane waves below make that distinction explicit.

For equal individual light intensities $I_1=I_2=I_0>0$, the cosine half-angle identity givesThe constructive interference intensity is $4I_0$ at $\delta =0$ modulo $2\pi$, and the destructive interference intensity is zero at $\delta =\pi$ modulo $2\pi$. At phase quadrature, $\delta =\pi /2$, it is $2I_0$. With $I_1=1$ and $I_2=0.25$ in the same arbitrary wave intensity unit, the maximum is $2.25$, the minimum is $0.25$, and the phase quadrature value is $1.25$: unequal wave amplitudes cannot cancel completely.

Both curves use the same wave intensity unit and a horizontal angle in radians. The dashed levels are the averages over a full phase difference cycle. Equal beams have interferometric visibility one; the unequal example has interferometric visibility $0.8$.

A spatially varying phase difference turns the same law into interference fringes. Take two ideal plane electromagnetic waves with common wavelength $\lambda$, equal individual light intensities $I_0$, and propagation directions at angles $+\theta$ and $-\theta$ from the $z$ axis in the $xz$ plane, with $0<\theta <\pi /2$. Both have linear polarization along $y$, perpendicular to both propagation directions. Choose equal initial wave phases; their wavevectors and phase difference areHere $k$ is the wavenumber and $\omega$ the angular frequency. The common $z$ and time terms cancel, so the interference fringes stay at fixed $x$ for every $z$. Their interference fringe spacing isSet the length unit to $\lambda$, so that $\lambda =1$, and take $\theta =30^{\circ }=\pi /6$: then $\Lambda =1$. The constructive interference intensity bands lie at integer $x/\lambda$; the destructive interference intensity bands lie at half-integer $x/\lambda$.

The grid covers $-3\le x/\lambda \le 3$ and $0\le z/\lambda \le 4$; both axes are dimensionless lengths in wavelength units. The colour scale is $I/I_0$ from zero to four, using the same single-beam normalization as the equal-beam curve above. This is an ideal plane wave superposition, with no slits or point sources and no near-field approximation. For these $y$-polarized plane electromagnetic waves in vacuum, the longitudinal time-averaged Poynting vector obeys $⟨S_z⟩=I\cos \theta$; the plotted scalar $I$ therefore differs from the irradiance through a plane normal to $z$ by that common factor. This relation follows by adding the magnetic fields as well as the electric fields.

The bright peaks do not create energy. Averaging over one interference fringe spacing, or equivalently over one full phase difference cycle, givesbecause the integral of the cosine is zero, with sine as its antiderivative. Equal beams thus average to $2I_0$, despite bright bands at $4I_0$. For the crossing-wave geometry, the average longitudinal energy flux is correspondingly $2I_0\cos \theta$. Bright and dark regions redistribute the incoming energy; conservation of energy concerns the full energy flux, rather than the height of one peak.

Wave interference is useful because phase differences convert small changes in optical path length into measurable changes of wave intensity, as in a Mach-Zehnder interferometer. Its beauty is that the same superposition principle produces both complete cancellation and fourfold brightening; adding more coherent contributions also produces the narrow peaks of a diffraction grating. Feynman's treatment connects this many-source pattern directly to addition of complex amplitudes. ^[ref].

A short proof outline uses local real fields $E_j(t)=A_j\cos (\omega t+{\varphi }_j)$ with parallel linear polarizations. Write $I=C⟨(E_1+E_2{)}^2{⟩}_t$ and $I_j=C{A}_j^2/2$ for the same positive conversion factor $C$. Averaging over many wave periods while the phase difference stays fixed uses $⟨{\cos }^2(\omega t+{\varphi }_j){⟩}_t=1/2$ and $⟨\cos (\omega t+{\varphi }_1)\cos (\omega t+{\varphi }_2){⟩}_t=\cos \delta /2$, giving the stated law. Equivalently, add the complex amplitudes and take the squared complex modulus. ^[ref].

For different fixed linear polarizations in a common transverse plane, the cross term acquires the dot product of their real unit vectors, equal to $\cos \alpha$ for relative angle $\alpha$. Orthogonal linear polarizations remove the cross term in polarization-insensitive scalar detection; different orientations generally only reduce it. A common polarizer can select overlapping components and recover wave interference. Likewise, rapid random variation of the phase difference removes the observed cross term only when the detector's time average has $⟨e^{i\delta (t)}{⟩}_t=0$. Slowly varying wave phases can instead move resolved interference fringes, and biased or incompletely sampled wave phases can leave residual contrast. Finite optical coherence therefore depends on the observation time as well as the fields.

A linear filter obeys the superposition principle: filtering $\alpha x+\beta z$ gives $\alpha T[x]+\beta T[z]$. Many familiar examples are also linear time-invariant systems, represented by an impulse response and convolution. Linearity alone does not imply time invariance: multiplying $x[n]$ by the index $n$ is linear but changes its behavior when the same input is shifted.

The moving average is a finite impulse response example implemented by arithmetic; a band-pass filter can instead be built from a series RLC circuit. Their frequency responses explain which oscillations pass through. Choosing a filter means choosing how much amplitude, timing and detail may change, not simply removing every visible fluctuation.

A linear partial differential equation is a partial differential equation in which the unknown function and its partial derivatives occur linearly, with coefficients independent of the unknown. Its homogeneous solutions obey the superposition principle.

For example, the heat equation permits adding independently decaying Fourier modes. For an equation $Lu=f$ with a fixed nonzero PDE source term, adding two solutions gives $L(u+v)=2f$, so the sum generally solves a different problem. Likewise, adding solutions adds their initial conditions and boundary conditions.

A linear time-invariant system combines two properties. Linearity means that the output from a linear combination of inputs is the same linear combination of their outputs. Time invariance means that shifting an input by an integer delay shifts its output by exactly that delay. The operation has no preference for a particular time origin.

For discrete time, write the unit impulse as $\delta [n]=1$ at $n=0$ and zero otherwise, a Kronecker delta. Its response $h[n]=T[\delta ][n]$ is the impulse response. A finite input can be decomposed into shifted impulses, givingThe superposition principle supplies the sum and time invariance supplies the shifted copies of $h$. Infinite sums need convergence assumptions; for a bounded input and an impulse response $h$ with absolute summability, the discrete convolution converges absolutely. For a physical system with stored initial energy, this construction describes the response with zero initial state; a separate free response accounts for nonzero initial conditions.

The gain is reuse. Measure or calculate one impulse response, and the convolution predicts the response to every allowed signal without solving a new dynamical problem each time. Alternatively, the frequency response predicts how each oscillatory component changes. ^[ref].

For the wave equation ${\psi }_{tt}=c^2{\psi }_{xx}$ on $0<x<L$, impose homogeneous Neumann boundary conditions ${\psi }_x(0,t)={\psi }_x(L,t)=0$, where $c,L>0$. Separation of variables gives $X^{\prime \prime }+\lambda X=0$ and $T^{\prime \prime }+c^2\lambda T=0$. Integration by parts gives ${\int }_0^L(X^{\prime }{)}^{2}dx=\lambda {\int }_0^L{X}^{2}dx$, so every nontrivial spatial eigenfunction has $\lambda \ge 0$.