Photon gas 2026-10-13
A photon gas is a collection of photons described statistically as radiation. Here the equilibrium model is homogeneous, isotropic electromagnetic radiation in vacuum at temperature : each photon has energy , with momentum and frequency , and there are two independent transverse states of electromagnetic polarization. An absorbing cavity with walls at a common temperature provides a way to establish this thermal equilibrium. The bulk formulas assume cavity dimensions much larger than the typical wavelength , so the electromagnetic mode density in a box can replace individual cavity modes by a continuum. [ref]
The crucial constraint is whether photon emission and photon absorption, or other number-changing processes, can equilibrate the photon population. When they can, photon number is not an independently conserved quantity and the equilibrium photon chemical potential is . For each polarization and frequency, the mean occupation number then follows the Bose-Einstein distribution:
If number-changing processes are ineffective but number-conserving interactions redistribute energy, a constrained Bose-Einstein distribution can instead have with chemical potential fixed by the number constraint. Thus thermal equilibrium with respect to energy exchange alone does not establish . The zero value can also persist in an already prepared Planck photon distribution during suitable subsequent evolution; continual photon emission is not required to preserve it. [ref]
Combining with the two-polarization density of states, , gives the Planck photon distribution and the spectral photon energy density:
Here is the spectral radiance in Planck's law; is stored energy density per unit frequency, not outgoing power per area. Their different units are connected by the full-sphere solid angle and the speed of light . The Boltzmann constant is , and is the Planck constant. [ref]
The equilibrium equation of state is particularly simple:
The radiation constant is the same constant called in black-body radiation; is the Stefan-Boltzmann constant. An outline of the energy density calculation substitutes and uses . The radiation pressure follows because the angular average of the squared directional component of photon momentum is for isotropy; the thermodynamic derivation of blackbody radiation pressure supplies an independent route. These formulas describe the radiation contribution, not the cavity walls. The isotropic massless relation survives a different isotropic spectrum, whereas the coefficient in requires the zero-chemical potential Planck photon distribution.
For a computed example with stipulated , the photon energy density is and the radiation pressure is . The photon number density in the International System of Units is
Here is the Riemann zeta function. The mean photon energy differs from the energy at the peak of a spectral curve: counting photons and weighting them by energy emphasize different frequencies. The calculations use the exact International System of Units values of the Planck constant, Boltzmann constant and speed of light. [ref]
Figure 1.
Analytic photon gas spectra: counting photons and weighting their energies
. Original equilibrium model, not measured data. The left panel evaluates at stipulated per gigahertz. The right panel uses and plots the normalized spectral fractions and , each with total integral one on . Their different shapes show why a photon-counting spectrum and an energy density spectrum have different peaks. Details.
A motivating observation is the cosmic microwave background spectrum measured by the Far Infrared Absolute Spectrophotometer (FIRAS) aboard COBE. Its apparatus used a polarizing Michelson interferometer, with one input viewing a roughly seven-degree sky patch and the other an internal blackbody. Scanning mirrors changed the optical path length difference; heat-sensitive detectors recorded the varying power, and Fourier transform spectroscopy recovered the spectral difference. Operation outside the terrestrial atmosphere removed atmospheric emission and absorption, while shielding and cooling the optics in a helium cryostat suppressed the instrument's own radiation. [ref]
Absolute calibration came from moving an external blackbody into the sky aperture, filling the same beam used for the sky. This was an iron-loaded epoxy absorber shaped as a re-entrant cone. Repeated reflections made it a close approximation to a blackbody, while thermometers supplied its temperature. Mather and colleagues explicitly evaluated reflections, temperature gradients, radiation leaking through the material and leakage around its edge. These are limits on how faithfully the calibrator realizes thermal equilibrium radiation, rather than details that can be eliminated merely by calling it black. [ref]
During sky observations the internal blackbody temperature was adjusted to nearly null the spectral difference, reducing gain and drift errors. Fixsen and colleagues' 1996 analysis used calibrated spectra over wavenumbers to , about to . They separated the sky's spatially uniform component, its dipole and Galactic emission, and varied sky cuts and foreground templates as controls. The resulting cosmic microwave background spectrum agreed with a single Planck function, with reported RMS deviations below 50 parts per million of its peak brightness. This is a residual statistic over the analyzed spectrum, not a pointwise relative error bound at every frequency. [ref]
The observation motivates using a photon gas as a physical model rather than merely a formal counting exercise. Its interpretation has a clear boundary: agreement with Planck's law constrains departures from an equilibrium spectral shape within the calibrated band and foreground assumptions; it does not directly measure radiation pressure or entropy density of a photon gas. Nor does it establish that today's cosmic microwave background is continually re-equilibrating with matter. A previously thermal spectrum can retain its Planck function shape as expansion of the universe lowers every photon frequency in the same proportion. The usefulness of the model is that one measured spectral temperature, together with its equilibrium assumptions, predicts energy density, radiation pressure, photon number density and entropy density of a photon gas consistently.
Planck's law 2026-10-13
Planck's law gives the spectral radiance of an ideal blackbody in thermal equilibrium at absolute temperature . It predicts how much power is carried in each small frequency or wavelength interval, per projected area and per solid angle. The whole curve is determined by : no material-specific parameter enters this ideal Planck function.
Writing for frequency in hertz and for wavelength in metres, the two forms of the Planck function are
Here is the Planck constant, the Boltzmann constant, and the speed of light in vacuum. Their exact values in the International System of Units are , , and . [ref]
The units matter: is measured in , whereas is measured in . Both are spectral radiance, rather than radiant exitance or energy density. To report per micrometre, multiply its per-metre value by ; a micrometre-wide interval is metres wide.
For a computed example, at and , the exponent is . Substitution into Planck's law gives , or per micrometre. At the same wavelength, the illustrative temperatures and give and in the latter units. These are ideal predictions for stipulated temperatures.
Figure 1. . The vertical axis is spectral radiance per micrometre, including the conversion from per metre. Increasing temperature raises the curve at every fixed wavelength and moves its marked peak toward shorter wavelengths. The peaks are approximately , , and for , , and , respectively; only to of each spectrum is shown. Details.
The key change of variables is . Corresponding positive intervals carry the same power, so and
Equal-width wavelength intervals have unequal widths in frequency; the factor compensates for this. Merely replacing the horizontal-axis labels of by does not produce . This explains the frequency and wavelength peaks of Planck's law.
A short derivation uses the Boltzmann distribution for a single electromagnetic mode with thermal energies , . Put . The geometric series gives and mean thermal energy . The electromagnetic mode density in a box, including both independent polarization of an electromagnetic wave states of a photon, is per volume per unit frequency. Multiplication gives spectral energy density ; isotropy gives . This counts thermal excitation energy, excluding zero-point energy. [ref]
Planck's law connects quantum mechanics to a measurable radiation curve. Its small- limit is the Rayleigh-Jeans law, ; its large- limit is the Wien approximation, . The suppression at high frequency makes the total energy density finite. Integrating the curve gives the Stefan–Boltzmann law, while locating its peak gives Wien's displacement law. Thus one formula predicts both the spectral shape and the total emission.
Point emitter 2026-10-18
An ideal point emitter emits radiant flux from one location. It approximates a physical source when its emitting dimensions are small relative to the distances relevant to the calculation. With isotropy, its radiant intensity is for total radiant flux . An extended source generally requires integrating contributions from different emitting positions.
Pressure Created 2026-09-24 Updated 2026-10-14
Pressure measures the compressive normal force per unit surface area. For a uniformly loaded surface of surface area , the magnitude of the perpendicular force is ; locally, . Its SI unit is the pascal, . A force parallel to the surface contributes shear stress, rather than pressure. In a gas, the steady macroscopic force comes from many microscopic transfers of momentum; ideal gas pressure from wall collisions makes this connection explicit. [ref]
For a computed example, a uniform pressure of on a piston of surface area produces a perpendicular force of . If the opposite face is exposed to , the net force due to the two pressures is , directed toward the lower-pressure side. The pressure on either face and their difference answer different questions; gauge pressure records a difference from a specified reference, often atmospheric pressure.
A scalar pressure describes a gas in local thermodynamic equilibrium with isotropic random particle velocities: every wall orientation has the same normal force per unit surface area. Without this isotropy, normal stresses can depend on direction and a pressure tensor is needed. Thus random motion alone is insufficient; a preferred direction matters.
In thermodynamics, pressure and volume are thermodynamic conjugate variables. For a simple compressible system, reversible expansion delivers work , and at fixed entropy and particle number ,
where is internal energy. A large pressure therefore means both a large wall force per unit surface area and a large amount of expansion work per small increase in volume.
Radiant intensity 2026-10-14
Radiant intensity describes how a source distributes its radiant flux over directions. A small cone of solid angle carries approximately when varies little across it; integrating over all emitting directions recovers the total radiant flux. For a source with isotropy emitting over a full sphere, and a cone carries . A planar source satisfying Lambert's cosine law instead has direction-dependent radiant intensity, even though its radiance is direction-independent.
For radiation in vacuum with isotropy, energy is distributed equally per solid angle: write for energy density per unit solid angle. Let be the angle between a ray's path and the wall normal. Its normal velocity determines the arrival rate, while normal momentum per energy is . Their product is . For a perfectly reflecting cavity wall,
The leading two is the normal momentum reversal; the incoming hemisphere contains half of the field's directional energy density. This is also the normal component of the bulk pressure tensor, . Isotropy makes all three diagonal components equal, while their sum is because the squared direction components sum to one. The factor one third therefore depends on three-dimensional angular geometry, not on the shape of the blackbody radiation spectrum.
Figure 1. . Original angular model of a reflecting cavity wall. The left panel shows representative incoming paths, not a literal two-dimensional photon gas. The right panel uses the correct three-dimensional solid angle weight, normalized on an incoming hemisphere: . The normal pressure contribution adds , giving total weight . Details.
A perfectly absorbing wall exposed to this incoming hemisphere receives only from absorption. If it is in thermal equilibrium with the isotropic radiation, its outgoing emission supplies another of recoil, recovering . A cold absorber before equilibrium is a different boundary condition. Similarly, the one-way incident energy flux is , whereas the net energy flux across a plane in the undisturbed isotropic field is zero. Dividing by misses the extra normal-momentum projection and does not give the cavity pressure.
For a stipulated isotropic radiation energy density , the bulk pressure is regardless of its spectral shape. In thermal equilibrium, blackbody radiation supplies with radiation constant , so ; Stefan–Boltzmann law defines . The equilibrium photon gas derives this same equation of state from its free energy.[ref]
Refraction Created 2026-10-13 Updated 2026-10-14
Refraction changes the direction of a transmitted wave when it crosses between media with different phase velocities. For light at a planar interface, Snell's law relates the angles measured from the normal vector through the refractive indices. At normal incidence the direction remains unchanged even though the wavelength and phase velocity can change. The worked example under geometrical optics distinguishes transmission from total internal reflection.
For light, refraction explains why a submerged object looks nearer the surface and why an optical lens can form an image. A wavefront crossing a boundary must retain a consistent wave phase along that boundary: its normal spacing changes with phase velocity, and oblique light rays change direction. Normal incidence separates these two effects especially clearly.
Take stationary, transparent, homogeneous media with isotropy, linearity of their optical response, and real positive refractive indices, at one chosen frequency . With the vacuum speed of light and the vacuum wavelength, medium has
The frequency is unchanged at a stationary boundary; the phase velocity, wavelength and wavenumber can change. Here is a phase refractive index, not the group refractive index. If there is optical dispersion, evaluate at this same frequency on both sides.[ref]
For an original idealized air–water example, stipulate and , with vacuum wavelength . Then , the water phase velocity is , and its wavelength is . At an air angle of incidence , Snell's law gives a water angle of refraction . These stipulated values illustrate the relations; actual refractive indices depend on material conditions and frequency. By reversibility of an optical ray, the same path emerging from water has outgoing air angle .
Refractive index Created 2026-09-29 Updated 2026-10-18
The refractive index compares the phase velocity of light in a medium with the vacuum speed of light . For a transparent, uniform medium with isotropy and a positive-index propagation branch,
Here is the angular frequency, is the wavenumber, and is the speed of a surface of constant wave phase. The refractive index is a dimensionless quantity; it is a property of a specified material response at a specified frequency, rather than a universal number attached to the word “glass”. More general absorbing or negative-index media require the conventions discussed in refractive index and electromagnetic response. [ref]
For an original model example, choose at vacuum wavelength . The phase velocity is , the frequency is , and the medium wavelength is . These are stipulated model values, not measured properties of glass or water. A value at this single frequency does not determine a pulse's group velocity.
The usefulness of refractive index is that one measured response connects several different observations: Snell's law predicts refraction, optical path length predicts accumulated wave phase, and changes with frequency describe optical dispersion. Even at normal incidence, where a light ray does not bend, a changed refractive index changes wavelength and wave phase. The following comparison isolates that effect; the oblique geometry is developed in wavefront.
Figure 1.
Chosen refractive index : shorter wavelength, slower wavefronts, unchanged interface frequency
. The top panel compares equal normalized electric field components at one time on either side of a stationary interface at normal incidence. The lower panel tracks constant-wave phase crests in space and time. Curves represent incident and transmitted components, with their amplitudes normalized separately; reflected light and interface amplitude matching are omitted. The chosen single-frequency model contains no measured material data. Details.
For a uniform linear medium with isotropy and local constitutive response, Maxwell equations give . Here is permittivity and is magnetic permeability. Since , the relative permittivity and the relative magnetic permeability obey
If both responses are real and positive with negligible loss, and the ordinary positive-index branch is intended, then . If additionally , then . 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 and explicitly fix the phase convention for complex amplitudes. For with , its wave amplitude has factor . A forward attenuating branch has , while its wave phase is controlled by . The complex relation 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.
Relative refractive index 2026-10-18
The relative refractive index of medium two with respect to medium one is , using their absolute refractive indices at the same frequency. For transparent media with isotropy and positive refractive indices, Snell's law gives , with the angles measured from the interface normal. Thus a refraction measurement in air directly determines an index relative to air. Converting it to an absolute refractive index requires the refractive index of that air; replacing it by one is an approximation.
Consider a mechanically free isotropic model of a saturated ferromagnet, at fixed temperature and fixed magnitude of magnetization. Let be a unit vector, where is the saturation magnetization. The magnetostrictive strain is an infinitesimal strain tensor, measured relative to a reference with randomly oriented magnetic domains. This reference convention matters: it sets the orientational average of the preferred strain to zero. Define , the saturation magnetostriction, as the longitudinal relative length change in this model.
Rotational symmetry permits a symmetric second-rank tensor made from of the form , where is the Kronecker delta. If we restrict attention to a change of shape with zero first-order volume change, its trace gives . Requiring the longitudinal component to equal gives . Hence
The tensor product is symmetric and unchanged under . It uses two copies of the magnetic direction, so the third-rank isotropic tensor obstruction in isotropic linear magnetostriction vanishes does not apply. More generally, a small-magnetization expansion can contain ; isotropy alone does not require its trace to vanish. The displayed saturated formula isolates the volume-preserving orientational contribution, rather than asserting that every material has zero volume magnetostriction.
Along a measuring unit vector , with , the fractional length change is
For a stipulated and reference length , alignment gives and a perpendicular orientation gives . Rotating from perpendicular to parallel therefore changes the length by , even though the saturated longitudinal value relative to the random reference is . A negative reverses these elongation and contraction signs.
The eigenvalues of this infinitesimal strain tensor are , so to first order. This is not exact finite incompressibility: using the illustrative deformation gradient gives
A finite-deformation law would need higher-order corrections if exact volume conservation were required.
Solid angle Created 2026-09-29 Updated 2026-10-18
A solid angle measures how much of the surrounding directions a region occupies, as seen from a specified point . Draw a sphere of radius centered at and mark where those directions meet its surface. If the marked patch has surface area , then
A Euclidean similarity scales both and by the same factor, so the chosen radius does not affect the answer. On the unit two-sphere, the numerical value is simply the patch's surface area. Unlike an angle drawn in one plane, a solid angle measures coverage in two angular directions.
For example, the full sphere has surface area and hence solid angle steradians. A hemisphere occupies steradians. Three mutually perpendicular coordinate planes divide the directions into eight congruent pieces, so each octant occupies steradians. A solid angle of one steradian therefore covers of all directions, rather than one full turn.
This makes a geometric question answer a collection question. For an ideal point emitter with isotropy, equal solid angles carry equal radiant flux. A detector accepting directions of total solid angle intercepts the fraction of that emitted radiant flux, before losses. With direction-dependent radiant intensity , the accepted radiant flux is instead . The solid angle of a disk below gives an exact example where the familiar inverse-square law is only a distant-source approximation.
Solid angle of a disk 2026-10-18
Consider a flat circular disk of radius , with its center a perpendicular distance from the viewpoint . The line from to its center is normal to the disk. Its rim defines a right circular cone with and . The solid angle of a circular cone immediately gives
For , this is the same half-opening example: about steradians. Both the radius and distance enter; doubling both leaves the solid angle unchanged.
A useful local version applies to a small surface patch of area , at distance from . If is the angle between its unit normal and the line to , its unsigned geometric contribution is
The numerator is the projected area. For a visible patch, perpendicular projection followed by radial projection onto a small piece of a sphere gives this formula. Integrating counts the directions once provided no direction encounters multiple counted patches. For an opaque detector, use its unoccluded face. At grazing view, the projected area tends to zero.
For this on-axis disk, a ring of radius has and . Its polar coordinates gives , so the same answer follows directly from the surface integral:
The changing distance and obliquity are both already accounted for; an extra detector-wide cosine factor would count the projection twice.
Figure 1.
The exact solid angle of an on-axis circular detector approaches a hemisphere at close range
. Left: a meridian section, drawn with . Right: the geometric acceptance fraction for a point emitter with isotropy, compared with the distant-source approximation . These are calculated geometric curves, not measured detector efficiencies. The dashed approximation is displayed outside its accurate range to show how it overestimates acceptance. Details.
For , the Taylor series of gives
Thus the accepted fraction of an ideal point emitter's isotropic output is approximately , where . This is the geometry behind the inverse-square law for a small detector. At fixed , however, taking gives , so the accepted fraction approaches rather than diverging.
For a chosen numerical example, let and . The exact solid angle is approximately steradians, giving an acceptance fraction . An ideal point emitter with isotropy and total radiant flux would therefore send about onto the disk. The distant-source approximation gives , about high. This calculation assumes unobstructed propagation and counts incident radiant flux; photon absorption, detector response, and any optical attenuator require separate factors.
Stefan-Boltzmann constant Created 2026-10-06 Updated 2026-10-14
The Stefan-Boltzmann constant converts the fourth power of absolute temperature into blackbody radiant exitance:
Here is the Boltzmann constant, the Planck constant and the speed of light. In the current International System of Units, these three constants have exact defining values, so is exact as an expression in them and . The displayed finite decimal is an approximation, not an exact decimal or a measured uncertainty. The 2022 recommended-constants table lists an ellipsis and an exact uncertainty for .[ref]
A dimensional check uses , and . Consequently , exactly the units needed in the Stefan–Boltzmann law.
For an equilibrium photon gas, the energy density is with radiation constant . The radiation pressure is . These describe the radiation filling an equilibrium cavity; has units , whereas the outward radiant exitance through a small opening is . The factor comes from averaging propagation directions and their projected area over a radiation field with isotropy, not from reducing the opening's emissivity.
Surface tension Created 2026-09-24 Updated 2026-10-14
Surface tension is the tensile force per unit length across a line drawn in a liquid interface. For a clean interface with isotropy, it has the same value in every tangential direction: each side pulls the other along the surface, perpendicular to the line. A straight segment of length therefore transmits a force of magnitude from one interface. Here we assume constant , with fixed temperature and composition. The same coefficient measures the reversible work required to create surface area:
The equality of units follows from ; the equality of the mechanical coefficients follows from the moving-rim argument below. [ref] [ref]
For a concrete example, take an ideal flat soap film spanning a rectangular wire frame with a straight sliding rim of length . Pull the rim outwards in the plane of the film, perpendicular to the rim. Each of the film's two liquid–gas interfaces pulls the rim back towards the film with force , so the outward holding force is . The two faces are two interfaces, each with its own surface tension. With the assumed toy value and , the required force is ; one interface would contribute . This value of is a stipulated model parameter, rather than a measured value for a particular soap solution.
Figure 1.
Two film faces double the force on a sliding rim
. The plan view shows the outward pull and the inward surface tension force on a sliding rim; the two interfaces overlap in this view. The graph compares one and two interfaces at the same assumed constant , with rim length in millimetres and force in millinewtons. Details.
Why does an interface resist being enlarged? Its molecules have a different environment from those in the bulk liquid; changing that environment has an energy cost. In the constant-surface tension model, reducing surface area lowers surface energy. This connects an easily measured rim force to the shape of a droplet, the Young–Laplace equation for its pressure, and the restoring mechanism of capillary waves. It is useful because one coefficient links forces, energy and geometry. [ref]
The moving-rim argument also explains the factor of two without memorizing a special formula. An outward displacement creates surface area on each face, hence . Slow pulling without viscous dissipation requires work , giving . For one interface, the same argument gives . A surface tension gradient instead produces Marangoni stress; the constant- assumption excludes that additional effect.
Total internal reflection Created 2026-09-29 Updated 2026-10-18
Total internal reflection occurs when a wave travels from higher to lower refractive index above the critical angle for total internal reflection, leaving no propagating transmitted plane electromagnetic wave normal to the interface.
At an infinite planar interface, total internal reflection returns all incident normal power in the reflected plane electromagnetic wave. Assume homogeneous dielectrics with isotropy, linearity, no dielectric loss, magnetic permeability , and positive real refractive indices . Above the critical angle for total internal reflection, both TE polarization and TM polarization have unit optical reflectance despite the external evanescent wave.[ref]
This makes total internal reflection useful for guiding light in an optical fiber: suitable internal incidence angles keep energy inside successive boundaries. Its more surprising feature is that perfect optical reflectance coexists with an external electric field. The evanescent field in total internal reflection can interact with a nearby object, as the frustrated total internal reflection experiment demonstrates.
For an original idealized example, choose , at vacuum wavelength . At angle of incidence , Snell law for electromagnetic waves gives an angle of refraction , so a transmitted light ray still emerges. At , its direction becomes tangent to the boundary. At , the proposed sine of a real angle of refraction would be : the transmitted wavevector must instead have an imaginary normal component. These are predictions for the stipulated media, not measured properties of a particular block.
Figure 1.
Constructed ray geometry just below, at, and above the critical angle for refractive indices 1.5 and 1; arrows indicate directions rather than relative powers
. Details.
A real semicircular-block experiment isolates this threshold. The University of Maryland laboratory instructions place a laser, motor-driven rotation stage and photodetector together; an aperture limits the detector voltage. The semicircular block is centered on the rotation axis, with the photodetector aligned along its flat face. The instructions identify a strong grazing-beam signal and a weaker signal just before the critical angle for total internal reflection. The readout is detector voltage versus stage angle. Retracing the incident laser beam identifies normal incidence for alignment.[ref]
The semicircular shape supplies a geometric control: a light ray passing radially through the curved face enters at normal incidence, avoiding additional entrance refraction. The weaker signal near the flat face locates a limiting angle of refraction; a vanishing remote-detector signal alone would not establish that the external electric field vanishes. Beam width and alignment limit angular resolution; no measured threshold or decay depth is supplied.[ref]
Volume magnetostriction 2026-10-18
Volume magnetostriction is the relative volume change associated with a magnetic-state change. To first order it is the trace of the magnetostrictive strain. Zero trace is an additional restriction, not a consequence of isotropy alone.
Wavefront Created 2026-10-05 Updated 2026-10-14
A wavefront connects points with the same wave phase at one specified time. It is a surface, rather than the trajectory of a particular piece of light. For a scalar monochromatic wave written as , a wavefront is a connected level set of the phase function ; different surfaces can share the same wave phase modulo . A plane wave has planar wavefronts. An outward wave from an ideal point source in a uniform medium has spheres as wavefronts, while an optical lens can turn nearly planar wavefronts into converging ones. The wavefront describes wave phase, not wave amplitude: drawing it does not say how bright the beam is. [ref]
The relation to light rays becomes precise in a stationary medium with isotropy, negligible loss and a real positive scalar refractive index . In the geometrical optics limit, spatial changes in the refractive index and wave amplitude are slow on the local wavelength scale. For a single smooth propagation branch, write
Here is the angular frequency, is the vacuum speed of light, has units of length, and the last relation is the optical eikonal equation. The local wavevector is . Its direction is normal to a smooth wavefront, since a gradient is normal to its level sets. A forward light ray has unit tangent vector , so it follows this normal under the stated assumptions. Along that light ray, : the change in is the optical path length. This is the short geometric argument linking wavefronts, light rays and optical path length. [ref]
Figure 1.
Wavefronts and normal light rays at an interface: an original plane wave model of refraction
. Details.
The diagram shows incident and transmitted branches separately, with the reflected branch omitted. Its stipulated refractive indices are and , and its incident angle is to the interface normal. At a stationary interface, the frequency and tangential wavevector agree on the two sides: Snell's law gives , hence . The wavefronts turn while their normal spacing decreases from to , where is the vacuum wavelength. More explicitly, with normal to the interface and tangential, take for and for . These functions agree at ; the drawn wavefronts are at with . This is an exact piecewise plane wave phase construction for the stipulated transparent media, not measured material data. [ref]
The normal-to-wavefront rule requires its assumptions. In a medium with anisotropy, the wavevector remains normal to a smooth wavefront, but the group velocity and energy flux need not follow it. [ref] Near a wave caustic, where multiple light rays meet, one smooth phase function may no longer describe the field; diffraction and the full wave amplitude matter. A wavefront is also not generally a surface of equal pulse arrival time. In a medium with optical dispersion, a monochromatic wave phase moves with the phase velocity, while a sufficiently narrow wave packet follows the group velocity. Even identical optical path lengths at one frequency can give different group delays. The usual equal-travel-time picture needs a common launch reference and a nondispersive ray model; an arbitrary broadband disturbance has no single monochromatic wavefront.[ref]