The line element of a Riemannian metric is its local expression for squared infinitesimal length. Along a parametrized curve $x(t)$, the arc length is $\int \sqrt{g_{ij}(x(t)){\dot{x}}^i{\dot{x}}^j}dt$. For a circular cylinder of radius $R$ parametrized by angle $\theta$ and axial position $z$, the induced metric gives $d{s}^2=R^{2}d{\theta }^2+d{z}^2$. Unrolling with $u=R\theta$ yields the Euclidean form $d{u}^2+d{z}^2$, explaining the straight unrolled descriptions of geodesics on a circular cylinder.

Consider a planar surface with outward normal along $z$. For an auxiliary equilibrium electron gas just outside it, let $n_v$ be the total number density, counting both signs of $v_z$. Its Maxwell-Boltzmann velocity distribution has normal-component densityAn electron with $v_z>0$ crosses more area per unit time when it moves faster. The outgoing number flux therefore uses $v_z{f}_z$, rather than just the probability of $v_z>0$:This integral contributes only $T^{1/2}$. Holding a classical bulk number density fixed and imposing a normal kinetic energy threshold $W$ would instead give $\Gamma =n\sqrt{k_{B}T/(2\pi m_e)}{e}^{-W/(k_{B}T)}$; that argument alone does not produce $T^2$.

For the ideal surface model, measure external electron energies from the metal's chemical potential $\mu$. A freely propagating external state has $E-\mu =W+p^2/(2m_e)$. When $W\gg k_{B}T$, its Fermi-Dirac distribution reduces to the Maxwell–Boltzmann statistics tail. Counting two spin states per phase space cell of volume $h^3$ givesEach of the three Gaussian integrals supplies $T^{1/2}$. Combining this $T^{3/2}$ state count with the $T^{1/2}$ crossing factor givesThe auxiliary atmosphere supplies an equilibrium state-counting argument; an operating collector removes the outgoing electrons. It need not maintain a real equilibrium gas throughout the gap. Unit transmission and the ideal free electron state count are substantive assumptions, so this calculation does not determine an arbitrary real surface's measured $A$.

In particular, electrons inside a metal generally form a degenerate Fermi gas. A classical distribution for the escaping tail does not imply a classical distribution for the whole interior. Richardson himself emphasized this distinction in his lecture.^[ref]

Choose a fixed positive reference $A_{\mathrm{ref}}$ with the same units as $A$. Taking the natural logarithm of a dimensionless ratio givesThus the Richardson plot is a straight line when $A$ and $W$ are constant, with slope $-W/k_B$ and intercept $\ln (A/A_{\mathrm{ref}})$. The familiar notation $\ln (J/T^2)$ silently chooses units; the reference makes that choice explicit. Changing $A_{\mathrm{ref}}$ shifts the intercept and leaves the inferred work function unchanged. If the horizontal axis instead uses $1000\mathrm{K}/T$, its slope is $-W/(1000\mathrm{K}{k}_B)$.

For the chosen $2.0\mathrm{eV}$ example, choose $A_{\mathrm{ref}}=120\mathrm{A}\mathrm{c}{\mathrm{m}}^{-2}{\mathrm{K}}^{-2}$. Then $y(1000\mathrm{K})=-23.2090$ and $y(1200\mathrm{K})=-19.3409$. Their slope is approximately $-23209\mathrm{K}$, giving $W=2.0\mathrm{eV}$. Two temperatures determine a line algebraically; several temperatures allow curvature and drift to be detected. These values are calculated, rather than experimental observations.

A straight trace alone does not establish a unique surface work function. Adsorption, changing surface coverage, or a mixture of regions with different barriers can change the inferred slope. For independent emitting patches with area fractions $f_i$, constant $A_i$ and barriers $W_i$,Lower barriers contribute disproportionately. Even with fixed patches, the slope of the Richardson plot gives a temperature-dependent weighted mean of the $W_i$, rather than their simple area average. If $A$ or $W$ themselves vary with temperature, the apparent slope-derived barrier becomesCheck collection, temperature calibration, and surface stability before assigning curvature to a new emission mechanism.

The coefficient $A$ has dimensions of electric current density divided by squared temperature. For ideal nonrelativistic particles with the free electron mass, two spin states, a planar uniform surface, and unit escape transmission, the thermionic tail-flux derivation giveswhere $e$ is the positive elementary charge and $h$ is the Planck constant. The two area units differ by a factor of $10^4$. A fitted $A$ need not equal $A_0$: reflection at the surface, the electronic band structure, and uncertain emitting area affect the inferred coefficient. A constant escape transmission $\tau$ would give $A=\tau A_0$ in this simple model; an energy-dependent transmission generally requires an integral instead.

The name joins Owen Willans Richardson with Saul Dushman. Dushman's 1923 paper develops a thermodynamic treatment of electron emission analogous to evaporation and evaluates the prefactor using quantum ideas.^[ref]

Consider the original pair of polynomial equationsTreat them as polynomials in $y$ with coefficients in $\mathbb{R}[x]$. Their Sylvester matrix and resultant areEvery common solution must therefore have $x=-1$ or $x=1$. Back-substitution into $F=0$ gives $y=x^2=1$, and both pairs also satisfy $G=0$. The complete real solution set isHere elimination and direct substitution agree: subtracting $F$ from $G$ gives $2x^2-2$, exactly the same condition.

In general, first solve the one-variable resultant equation, then recover and check the other coordinate in both original polynomial equations. A real coordinate candidate can correspond only to nonreal values of the eliminated variable. Also, vanishing leading coefficients require the separate check in degree drop in a specialized resultant. The general elimination-and-back-substitution method is described in the lecture's §3.1.^[ref]

Write $f(x)=a_m{\prod }_{i=1}^m(x-{\alpha }_i)$ and $g(x)=b_n{\prod }_{j=1}^n(x-{\beta }_j)$ in an algebraic closure, counting multiplicity of a root. The coefficient-row convention above givesThis makes the shared-root of a polynomial test transparent: a factor vanishes precisely when one ${\alpha }_i$ equals one ${\beta }_j$. Swapping the two polynomials reverses the sign of all $mn$ differences, giving $\mathrm{Res}(g,f)=(-1{)}^{mn}\mathrm{Res}(f,g)$.^[ref]

The value measures an algebraic product, not a normalized geometric distance. Replacing $f$ by $cf$ for nonzero $c\in K$ leaves its roots of a polynomial unchanged but multiplies the resultant by $c^n$. Therefore a small numerical resultant alone does not quantify how close two sets of roots of a polynomial are. Exact zero testing and numerical sensitivity are different questions.

A multiple root of a polynomial is also a root of a polynomial of its formal derivative. In characteristic zero, for $\mathrm{deg}f=m\ge 2$, the associated discriminant isFor example, $h(x)=x^2+s$ gives $\mathrm{Res}(h,2x)=4s$ and $\mathrm{Disc}(h)=-4s$. The two roots of a polynomial merge precisely at $s=0$.^[ref]

The respiratory control ratio divides the oxygen consumption rate during ADP-stimulated respiration by that after supplied ADP has been used. It compares reaction rates, not ATP formed per oxygen.^[ref]

A large ratio means strong stimulation under the comparison's conditions; it is not a universal measure of all mitochondrial functions.

For a constructed example, choose constant suspension volume $V=2\mathrm{mL}$ and oxygen concentration declines of $0.030$ and $0.006\mathrm{mmol}{\mathrm{L}}^{-1}\mathrm{mi}{\mathrm{n}}^{-1}$ during the two comparison windows. Assume no gas exchange or dilution and a correctly calibrated readout. A derivative converts each concentration slope into an amount of substance consumption rate:The common volume cancels in the ratio, so dividing the positive slope magnitudes gives the same answer. These are chosen inputs, not the 1955 observations. If an uncorrected influx of $3\mathrm{nmol}\mathrm{mi}{\mathrm{n}}^{-1}$ of oxygen occurred in both windows, the apparent rates would instead be $57$ and $9$, giving $19/3\approx 6.33$. The same absolute error can inflate a ratio because it removes a larger fraction of its smaller denominator.

In Respiratory enzymes in oxidative phosphorylation. I. Kinetics of oxygen utilization (1955), Britton Chance and G. R. Williams examined isolated liver mitochondria in isotonic medium with substrate and inorganic phosphate. A one-millilitre optical cuvette held the suspension. A vibrating platinum electrode lay outside the light path; calibrated electric current followed dissolved oxygen, with rapid changes recorded by a galvanometer oscillograph. Air-saturated and oxygen-depleted readings checked sensitivity. Their Figure 2 shows ADP additions accelerating oxygen consumption, followed by return to slow respiration; repeated additions restored acceleration. Figure 4B compares delayed substrate addition after ADP incubation with the reverse order, supporting little unrelated ADP loss in that preparation. These comparisons support coupling to ADP utilization. The authors discuss dilution, sensitivity changes and atmospheric oxygen entry, which especially biases slow respiration. The readout did not directly measure the membrane potential or prove the later chemiosmosis mechanism.^[ref]

A control can test either the instrument or the biology. Calibrating electric current against oxygen checks the former; comparing responses before and after an ADP addition checks the latter. A valid calibration cannot by itself establish ATP formation, and a repeatable biological response cannot by itself fix an instrument's sensitivity.

Near a renormalization-group fixed point, expand a difference between two nearby models in eigenvectors of the stability matrix. An irrelevant operator direction has $y_i<0$, so its dimensionless coefficient shrinks under increasing length factor $b$. A relevant operator direction grows and must be tuned to study the same critical behavior. This is a mechanism for losing microscopic detail, not an assertion that every parameter disappears. Critical temperature, individual amplitudes and regular contributions to free energy can remain material dependent.

For the solvable example, write $A=ϵ/g_0-1$. Its exact solution gives $g-ϵ=-ϵA{b}^{-ϵ}+O(b^{-2ϵ})$. At the thermal stopping scale $b_{\ast }=|t_0{|}^{-1/y_t}$ this is of order $|t_0{|}^{ϵ/y_t}$. If the stopped correlation length is smooth in $g$, the leading correction to its power law has that order. Thus two microscopic couplings can yield the same limiting critical exponent, yet approach it at different rates.

This mechanism makes universality useful: one tractable representative can predict critical exponents and normalized scaling forms for many systems in its universality class. A liquid-gas critical point and a scalar magnetic transition can share a long-distance description even though their microscopic constituents differ. The scope must be checked: spatial dimension, order parameter symmetry, interaction range, and extra soft variables can change the renormalization-group fixed point.

In his lecture, Kenneth Geddes Wilson connects limiting renormalization-group fixed points with common critical behavior and states the limits of his elementary truncation.^[ref]

Add a thermal scaling coordinate $t$, with its critical value subtracted, satisfyingThis is a relevant direction of a fixed point. A thermal scaling coordinate agrees to leading order with a nonzero multiple of reduced temperature; in a general model it can mix several bare parameters. An ordering conjugate field is held at its critical value. The combined stability matrix at $(g,t)=(ϵ,0)$ has eigenvalues $-ϵ$ and $y_t$. Attraction along $g$ therefore does not make the critical point attractive against thermal detuning. The critical surface in this example is $t=0$.

Let $X=\xi /a$ denote the correlation length in microscopic lattice units. A renormalization-group transformation of length factor $b$ givesIndeed, a tail with exponential envelope $\exp (-r/\xi )$ becomes $\exp [-r^{\prime }/(\xi /b)]$ under $r^{\prime }=r/b$. Choose the stopping scale $b_{\ast }=|t_0{|}^{-1/y_t}$ so that the transformed thermal coordinate is $+1$ or $-1$. If the stopped correlation length tends to a finite nonzero value $C_{\pm }$ as the irrelevant coupling approaches $g_{\ast }$, thenThe limit requires $|t_0|\ll 1$ and enough room between $a$, $\xi$, and the sample size. Singular dependence on a dangerously irrelevant coupling or a marginal operator needs separate analysis. On an ordered branch, use the connected correlation function in a pure thermodynamic phase.

For a completely specified illustrative scaling model choose $y_t=3/2$ and $X=C_{+}{t}^{-2/3}$ for $t>0$. Then $t=10^{-3}$ gives $X=100C_{+}$, whereas $t=10^{-6}$ gives $X=10^4{C}_{+}$. With $b=4$, $t^{\prime }=8t$ and $X(t^{\prime })=X(t)/4$. Models with $C_{+}=1$ and $C_{+}=2$ have different lengths but identical curves after plotting $X/C_{+}$ against $t$. This collapse checks the chosen scaling law; agreement of real data would test that law, rather than follow from this example. Neither $ϵ=1/5$ nor $y_t=3/2$ is being fitted to a three-dimensional Ising model.

Take the following bounded analytic example for a dimensionless positive coupling $g$, with fixed $ϵ>0$ and $g(0)=g_0>0$:Increasing $\ell$ means increasing observation length $L$, hence motion toward the infrared fixed point. This convention matters: if the inverse-length scale is $\mu ={\mu }_0{e}^{-\ell }$, the same trajectory has renormalization-group beta function ${\beta }_{\mu }(g)=\mu dg/d\mu =-ϵg+g^2$. The two derivatives have opposite signs.

To solve the equation, set $v=1/g$. It obeys $dv/d\ell =-ϵv+1$, a linear ordinary differential equation. ThereforeFor $\ell \ge 0$ its denominator is positive, and $g$ stays between $g_0$ and $ϵ$. Thus the trajectory is bounded on this domain. The renormalization-group fixed points are $g=0$ and $g_{\ast }=ϵ$. The derivative of the length-flow right-hand side is $ϵ-2g$: zero is repulsive from positive $g$, whereas $g_{\ast }$ has eigenvalue $-ϵ$ and attracts positive trajectories. Linearization gives $g-g_{\ast }\propto e^{-ϵ\ell }=b^{-ϵ}$ near $g_{\ast }$. The inverse-length beta function has positive derivative at the same infrared-attractive point.

For the chosen $ϵ=1/5$, take two starting models with $g_0=1/20$ and $g_0=3/5$. At $b=1024$, $e^{-ϵ\ell }=1/4$, so their couplings are respectively $4/35$ and $6/25$. Both move toward $1/5$, from opposite sides. They are still visibly different: slow decay of an irrelevant operator can leave substantial corrections at accessible lengths. The case $g_0=0$ stays exactly zero and is outside this positive-coupling basin.

Use a dimensionless statistical Hamiltonian $\mathcal{H}=H/(k_{B}T)$, so that the partition function is $Z={\sum }_s{e}^{-\mathcal{H}(s)}$. Choose a normalized blocking kernel $B(S\mid s)$ with ${\sum }_{S}B(S\mid s)=1$. For the majority example it is one when every block has its prescribed majority, and zero otherwise. Define the coarse statistical Hamiltonian by the displayed sum. Then an elementary interchange of finite sums provesThe probability of a retained pattern is therefore its original marginal probability. For any observable $A(S)$ of the retained variables,This is the sense in which elimination preserves predictions for retained information. It does not reconstruct every microscopic observable from a block sign.

Next express positions in units of the enlarged spacing: $x^{\prime }=x/b$, where $b>1$. In a continuum version, rescale the retained scalar field as well to restore a chosen gradient energy coefficient; any change in the integration measure contributes to the field-independent term. The result is a new set of dimensionless couplings. Repeating the procedure defines a renormalization-group flow. Additional interactions and the additive constant generally appear, even if the starting statistical Hamiltonian was simple. Keeping only a few terms is an approximation; the equality of partition functions above applies to the complete coarse weight. See also the Wilsonian coarse-grained statistical Hamiltonian for integration over fast Fourier modes.

Consider Ising spins $s_i=\pm 1$ on a square lattice of spacing $a$. Divide it into disjoint $3\times 3$ blocks $B$ and retain the majority sign $S_B$. Each block contains nine spins, so a tie is impossible. For example, five positive and four negative spins produce $S_B=+1$. All $2^9$ configurations are assigned a retained sign; the individual pattern within a block is lost. A coarse-grained variable has fewer possible values than the microscopic configuration it summarizes.

The spacing between block centres is $3a$. In coordinates $x^{\prime }=x/3$ it is again $a$, so the retained pattern can be compared with the original lattice at the same numerical cutoff. This coordinate change does not physically shrink the sample. A majority picture alone does not specify the new interactions: the blocking kernel and original Boltzmann weights determine their probabilities. The figure is a constructed configuration, not a simulated equilibrium state or a claim that majority blocking closes the nearest-neighbour Ising model.

Renormalization of spontaneously broken gauge theories means that the Higgs mechanism can generate gauge boson masses while a fixed finite set of renormalized couplings, masses and field normalizations still suffices to absorb ultraviolet divergences. The symmetry constrains the counterterms even when the chosen vacuum expectation value hides it. This makes radiative corrections to the electroweak interaction calculable after inputs have been fixed by renormalization conditions.

A useful distinction is between choosing variables that display the physical particles and choosing variables that display good power counting. Unitary gauge removes the would-be Goldstone bosons from the explicit fields. R-xi gauge keeps them and the required Faddeev-Popov ghosts in intermediate Feynman diagrams, exposing simpler large-momentum behavior. The physical scattering amplitude must agree after all contributions and counterterms have been included consistently.

In Renormalizable Lagrangians for massive Yang-Mills fields (1971), Gerard 't Hooft constructed spontaneously broken gauge theories with massive gauge bosons. His worked SU(2) gauge theory with an adjoint Higgs field leaves a massless photon; it is a demonstration model, not the full chiral Standard Model. The gauge fixing, Faddeev-Popov ghosts and Ward identities connect renormalizability with physical unitarity.^[ref]

Regularization and renormalization of gauge fields (1972), by Gerard 't Hooft and Martinus J. G. Veltman, developed dimensional regularization with recursive subtraction of ultraviolet subdivergences. It continued Feynman integrals in the dimension, keeping the relevant Ward identities in non-anomalous cases; its final section examined quantum anomalies. The procedure does not make an inconsistent chiral gauge theory consistent: gauge anomaly cancellation remains essential.^[ref]

The organizing proof outline is to use good large-momentum propagators, subtract every divergent subgraph, and constrain the remaining local divergences by the Ward identities. A finite family of allowed counterterms must close under this operation at every loop order. The same identities must remove unphysical intermediate contributions from physical amplitudes. Power counting alone establishes neither closure nor these cancellations. Modern BRST symmetry packages these restrictions; the scalar subtraction below illustrates only locality.

A Carnot refrigerator runs a reversible Carnot cycle in the direction that removes heat $Q_c$ from a cold reservoir at $T_c$ and delivers $Q_h$ to a hot reservoir at $T_h>T_c$. It requires supplied work $W=Q_h-Q_c$. Zero total entropy production gives $Q_c/T_c=Q_h/T_h$.

Its coefficient of performance, defined as useful cooling divided by supplied work, is thereforeAll temperatures here are absolute. The ratio can exceed one because supplied work moves heat rather than being the only energy delivered to the cold-side process. For reservoirs at $250\mathrm{K}$ and $300\mathrm{K}$, the reversible value is five. A real refrigerator produces entropy and requires more work for the same cooling between those reservoirs. The ideal result is a limit, not a claim of finite cooling power without thermal gradients.

An Abbé refractometer places a small sample between an illuminating prism and a measuring prism, then locates the boundary between bright and dark parts of an angular field. That limiting ray determines the sample's refractive index through Snell's law and the known measuring-prism geometry.

The prism surfaces must be clean and the sample must make good optical contact: a trapped air layer introduces a different interface. Because refractive index varies with wavelength and temperature, a useful reading specifies both; an instrument using white illumination needs compensation for the colored boundary caused by dispersion. A reference liquid checks the reading before an unknown sample is measured.

A prism minimum-deviation index measurement infers a transparent sample's relative refractive index from two angles: its known apex angle $A$ and the smallest deviation ${\delta }_{\min }$ of a selected spectral line. A spectrometer combines a source and slit, a collimator, a rotating optical prism table, and a viewing optical telescope with an angular readout. The collimator supplies a nearly parallel beam; the optical telescope locates its emergent direction. This turns an otherwise difficult measurement of the phase velocity of light into a measurement of angles.

For the symmetric path through a uniform optical prism, with the same external medium at both faces and positive refractive indices,The angles $i,e,r_1,r_2$ are measured from face normals. This is the minimum deviation relation, not a formula for an arbitrary prism orientation. An absolute refractive index additionally requires $n_{\mathrm{ext}}$.

For example, stipulate $A=60^{\circ }$ and $n_{\mathrm{rel}}=3/2$. Then $i=e=\arcsin (3/4)\simeq 48.5904^{\circ }$ and ${\delta }_{\min }\simeq 37.1808^{\circ }$. Substitution recovers $n_{\mathrm{rel}}=1.5$. These synthetic values illustrate the geometry; they are not readings for the prism in the laboratory protocol below.

A Michigan State University teaching experiment uses a mercury lamp, slit, collimator, prism table and optical telescope, with a subsidiary angular scale read to one arcminute. First focus the optical telescope on a distant object; with the prism removed, align its crosshair to a narrow vertical slit image, focus the collimator, and establish the direct-beam zero. Insert and centre the prism. Select the green mercury spectral line, tabulated at $546.1\mathrm{nm}$, then rotate the prism while tracking its image with the optical telescope. Locate the turning point where continued prism rotation reverses the image motion, and read the deviation. Re-establish the minimum separately for each other spectral line. The protocol assumes $A=60^{\circ }$; that assumed apex angle is not an independently measured result. ^[ref]

For an uncertainty-aware implementation, repeat the turning-point search from both prism-rotation directions, recheck the direct-beam reference, and record individual readings before averaging. Independently determine $A$ or include its calibration uncertainty; record the chosen spectral line, sample identity, temperature and surrounding-air conditions. These are proposed controls deduced from the inference, not claimed additional instructions or results of that manual. Finite slit width, unresolved spectral lines, imperfect collimation, table tilt, angular calibration and sample nonuniformity can limit the inferred refractive index. One arcminute readout resolution alone does not establish one arcminute total measurement uncertainty. A single selected wavelength gives one index, not a measured optical dispersion curve.

A short geometric derivation uses $r_1+r_2=A$ and $\delta =i+e-A$. Reversibility exchanges entrance and exit. At the symmetric stationary path, $i=e$ and $r_1=r_2$; applying Snell's law gives the displayed index relation. The plotted calculation checks that this stationary point is a minimum on the chosen transmitting branch. To assess local sensitivity, differentiate the inferred index using angles in radians:This first-order relation exposes both angular contributions; it is not an instruction to add signed systematic errors as independent random errors. In the synthetic example, a deviation error of $+1$ arcminute at fixed $A$ changes the inferred relative refractive index by approximately $+1.92\times 10^{-4}$.

For a uniform linear medium with isotropy and local constitutive response, Maxwell equations give $k^2={\omega }^2\epsilon (\omega )\mu (\omega )$. Here $\epsilon$ is permittivity and $\mu$ is magnetic permeability. Since $c^{-2}={\epsilon }_0{\mu }_0$, the relative permittivity ${\epsilon }_r=\epsilon /{\epsilon }_0$ and the relative magnetic permeability ${\mu }_r=\mu /{\mu }_0$ obeyIf both responses are real and positive with negligible loss, and the ordinary positive-index branch is intended, then $n=+\sqrt{{\epsilon }_r{\mu }_r}$. If additionally ${\mu }_r\simeq 1$, then $n\simeq \sqrt{{\epsilon }_r}$. All these responses must be evaluated at the electromagnetic wave's frequency; a static relative permittivity cannot simply be substituted at an optical frequency. ^[ref]

With loss, use a complex index $\tilde{n}=n^{\prime }+i\kappa$ and explicitly fix the phase convention for complex amplitudes. For $E\propto \exp [i(kz-\omega t)]$ with $k=\tilde{n}\omega /c$, its wave amplitude has factor $\exp (-\kappa \omega z/c)$. A forward attenuating branch has $\kappa \ge 0$, while its wave phase is controlled by $n^{\prime }$. The complex relation ${\tilde{n}}^2={\epsilon }_r{\mu }_r$ alone leaves a choice of sign: the physical solution also needs causality and the appropriate radiation or decay condition. In negative-index media the physically selected real part can be negative. The positive square root is therefore not a universal rule for complex or negative-index media. Direction and polarization of an electromagnetic wave also matter in a plane wave in an anisotropic dielectric, for which one scalar index need not describe every propagation branch.

For a linear stationary interface, phase matching at a planar wave interface preserves frequency: the incident and transmitted electric fields must oscillate together at each boundary point. Their phase velocities can differ, so their wavelengths differ. Writing $\omega =2\pi f$, the definitions giveThis explains the compressed wavefront spacing in the figure without changing the time between crests at the interface. Moving interfaces or time-dependent material responses need not preserve frequency. ^[ref]

The notation $n({\lambda }_0)$ labels refractive index by vacuum wavelength; it does not mean the wavelength inside the material. Dependence on ${\lambda }_0$, equivalently on frequency, is optical dispersion. The field driving the material induces oscillating electric charges; their emitted electromagnetic fields combine with the applied field to change the resulting wave phase. This collective response motivates refractive index, rather than a picture of independently delayed little pieces of light. ^[ref]

To distinguish three speeds, track a wavefront for phase velocity, a narrow wave envelope for group velocity, and the earliest newly introduced disturbance for front velocity. For a transparent differentiable propagation branch,The denominator is the group refractive index. It describes envelope transit when loss and reshaping are small. It need not equal $n$, and neither phase velocity nor a reshaped pulse peak generally measures signal arrival. The front velocity article explains the separate relativistic causality constraint and why a detector threshold is not an exact causal-front measurement.