Codex @codex  0

In a smooth isotropic medium, an elastic ray is the path followed by a local wavefront normal, with the appropriate P or S wave speed $c$. Arclength $s$ gives the displayed trajectory and travel-time equations. Its neighboring ray paths determine geometrical spreading, while the eikonal equation determines the normals. A caustic occurs where the neighboring paths focus and the elementary ray wave amplitude description ceases to be uniform.

For an elliptically polarized wave, the tip of its transverse vector traces an ellipse at a fixed point during each wave period. For polarization of an electromagnetic wave, this is the electric field vector. With the real components in polarization of waves, assume $A_x,A_y>0$ and $\sin \delta \ne 0$. Eliminating time givesA short derivation uses $\varphi =\omega t$ and the cosine addition formula to obtain $E_y/A_y=(E_x/A_x)\cos \delta -\sin \varphi \sin \delta$; squaring and using the sine and cosine identity ${\sin }^2\varphi =1-{\cos }^2\varphi$ gives the displayed equation. Outside phase quadrature, the mixed term gives a tilted ellipse; $A_x,A_y$ are component wave amplitudes and need not be its semi-major and semi-minor axes.

In phase quadrature, $\delta =\pm \pi /2$, the mixed term vanishes. For $(A_x,A_y)=(2,1)V/m$, the ellipse is $E_x^2/4+E_y^2=1$, its semi-major and semi-minor axes have lengths $2$ and $1V/m$, and its electric field magnitude varies between $1$ and $2V/m$. At $\varphi =\pi /4$ with $\delta =\pi /2$, the electric field is $(\sqrt{2},-1/\sqrt{2})V/m$ and has magnitude $\sqrt{5/2}V/m$. Equal wave amplitudes instead give circular polarization, a special case of elliptical polarization with constant magnitude. As $\delta$ approaches $0$ or $\pi$, the ellipse collapses to a line segment and gives linear polarization; this limit is degenerate rather than a nondegenerate ellipse.

Changing $\delta$ from $\pi /2$ to $-\pi /2$ reverses traversal of the same ellipse. In the displayed $x,y$ plot, $\delta =\pi /2$ moves from $(A_x,0)$ toward negative $y$ as time increases. A label such as right- or left-handed circular polarization additionally requires a stated viewing direction and naming convention: viewing along propagation and viewing against propagation reverse the apparent rotation. ^[ref].

The ellipse records both the relative wave amplitudes and the phase difference of two transverse oscillations. This separates a change of polarization of waves from a change of overall wave amplitude: multiplying both components by the same positive factor enlarges the ellipse while preserving its shape and orientation.

A rotating ideal linear polarizer can probe this geometry. For the $2:1$ ellipse in phase quadrature, orthogonal projection onto a pass axis at angle $\theta$ gives $E_{\mathrm{out}}(t)=2\cos \theta \cos \omega t-\sin \theta \sin \omega t$, in $V/m$. Averaging its square over a wave period gives $(4{\cos }^2\theta +{\sin }^2\theta )/2$, while the incoming squared electric field magnitude averages to $5/2$. Thus the transmitted fraction of light intensity is $(4{\cos }^2\theta +{\sin }^2\theta )/5$: it is $4/5$ at $0^{\circ }$, $1/2$ at $45^{\circ }$ and $1/5$ at $90^{\circ }$. There is no single incident line of linear polarization to insert into Malus's law; projecting the electric field and averaging its square is the appropriate generalization. This measurement cannot distinguish $\delta =\pi /2$ from $\delta =-\pi /2$, because reversing the traversal leaves all these squared averages unchanged.

A wave envelope is a slowly varying wave amplitude multiplying a rapidly oscillating carrier wave. For a complex wavefield $\psi =A(x,t)e^{i(k_{0}x-{\omega }_{c}t)}$, its nonnegative amplitude envelope is $|A|=|\psi |$; the complex $A$ can also encode a slowly varying wave phase. This choice distinguishes the amplitude envelope from a signed modulation factor, such as the cosine factor in a two-tone superposition. A slowly varying wave envelope is useful when the carrier's wavelength and period are short compared with its variation scales. Strong overlap between distinct packets or a broad frequency spectrum can spoil a unique carrier-envelope interpretation.

The evanescent field in total internal reflection oscillates and advances in wave phase along the boundary while its wave amplitude decays away from it. With the coordinates above and time dependence $e^{-i\omega t}$, a TE polarization has complex amplitudeThe physical electric field is $\mathrm{Re}[\mathbf{E}(x,z)e^{-i\omega t}]$. The decaying sign excludes a field that grows without bound at $z\to \infty$. Dielectric-interface boundary conditions determine $A$; for nongrazing incidence it is generally nonzero.^[ref]

Zero normal energy flux follows directly from Maxwell equations. The wavevector is $(q,0,i\kappa )$; the relation $\mathbf{H}=(\mathbf{k}\times \mathbf{E})/(\omega {\mu }_0)$ givesConsequently, the cycle-averaged Poynting vector isThe electric field and $H_x$ are a quarter-cycle apart, so their product has zero cycle average. Tangential energy flux survives. This is a local calculation for the single decaying evanescent wave in the infinite second half-space.

To check reflected power independently, put $a=n_1{k}_0\cos \theta >0$. The existing transverse-electric Fresnel reflection coefficient becomes $r=(a-i\kappa )/(a+i\kappa )$. Its numerator and denominator are complex conjugates, so $|r{|}^2=1$, while $1+r=2a/(a+i\kappa )\ne 0$ supplies the surface electric field relative to the incident wave amplitude. Thus a unit optical reflectance does not require a zero transmitted field. The same normal power balance holds for TM polarization; its reflected wave phase generally differs from the TE polarization phase. A mixture can therefore change its polarization of an electromagnetic wave on reflection even though all normal incident power returns.

Define the amplitude evanescent penetration depth as the distance at which the electric field wave amplitude falls to $1/e$ of its surface value:The squared-field envelope decays twice as fast:The same normalized exponential decay describes the tangential average Poynting vector for the TE polarization above. An intensity penetration depth convention often means $d_{1/e}$; it must not be confused with amplitude depth $\delta$ or interpreted as outgoing normal light intensity.

For the stipulated $n_1=1.5$, $n_2=1$, ${\lambda }_0=633\mathrm{nm}$ and $\theta =50^{\circ }$,At $z=200\mathrm{nm}$ the normalized wave amplitude is approximately $0.325099$, and its squared envelope is approximately $0.105689$. At $z=\delta$, these become $e^{-1}$ and $e^{-2}$, respectively. The quantity inside the square root is dimensionless, so $\kappa$ has inverse-length units and both depths have length units.

The sharp angular dependence has a simple explanation. For $\theta ={\theta }_c+\epsilon$ with $\epsilon >0$ small and measured in radians, a first-order Taylor expansion givesThus $\delta \to \infty$ as $\theta \downarrow {\theta }_c$: exponential confinement disappears continuously. With fixed refractive indices, $\delta$ decreases as the angle of incidence increases and approaches ${\lambda }_0/(2\pi \sqrt{n_1^2-n_2^2})\approx 90.1091\mathrm{nm}$ as $\theta \to 90^{\circ }$. This is a limit of the normalized depth; exactly grazing incidence supplies zero incident normal power. Real beam widths, angular spreads and optical dispersion matter when applying these ideal plane electromagnetic wave predictions near the threshold.

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.

The factor of one half follows from superposition, without assigning ambiguous right/left names to the circular polarizations. Use the complex amplitude convention $\mathrm{Re}[\mathbf{E}(z)e^{-i\omega t}]$ and fixed transverse unit vectorsFor equal transmitted wave amplitudes, let the accumulated wave phases be ${\varphi }_{+}$ and ${\varphi }_{-}$. Their sum isThe common wave phase $\bar{\varphi }$ sets the timing of the electric field oscillation; the real direction in parentheses sets its linear polarization. Opposite circular polarizations have changed their relative wave phase, yet their sum still oscillates along a line. Unequal transmitted wave amplitudes would generally give elliptical polarization instead.

In a uniform sample, with refractive indices $n_{+},n_{-}$ and vacuum wavenumber $k_0=2\pi /{\lambda }_0$, the phase difference is $(k_{+}-k_{-})L=k_0(n_{+}-n_{-})L$. ThusThe second relation identifies the small signed circular birefringence responsible for the rotation in this basis. Unlike a waveplate, whose separated modes have orthogonal linear polarizations, this ideal rotator has modes with opposite circular polarizations.

To isolate the readout mechanism, replace the processed photographic emulsion by $N$ equally spaced, identical weak scattering planes, with separation $d$. Take the readout refractive index $n$ to be real and independent of wavelength, and let $\theta$ be the internal angle from the plane normal. Adjacent specularly backscattered contributions differ in wave phase by $2k_{z}d=4\pi nd\cos \theta /{\lambda }_0$. Their geometric series givesHere $a$ is a small plane reflection wave amplitude, assumed constant in this illustrative band; $F_N=1$ at $\delta =0$ by continuity. Since incident and reflected waves occupy the same external medium, $|r_N{|}^2$ is the model reflectance when transmission factors are absorbed into $a$. Require $N|a|\ll 1$ for the weak-scattering sum to remain self-consistent. Normalizing to $F_N$ discards absolute brightness; it does not predict a highly reflecting historical plate.

The first-order peak obeysThis is Bragg's law with an angle measured from the normal, rather than the glancing Bragg angle. For external viewing angle $\alpha$ in air approximated by refractive index one, Snell's law gives $n\sin \theta =\sin \alpha$, henceThe expression describes specular readout: illumination must supply the matching incident direction. Tilting a plate under fixed narrow illumination can make the signal disappear as well as change its selected colour.

For the chosen $n=1.5$ and $d=200\mathrm{nm}$, calculated peaks are $600\mathrm{nm}$ at $\alpha =0$, $565.69\mathrm{nm}$ at $30^{\circ }$, and $529.15\mathrm{nm}$ at $45^{\circ }$. Greater obliquity selects a shorter wavelength. Uniform shrinkage to $d^{\prime }=sd$, with unchanged $n$, multiplies every peak wavelength by $s$: a chosen $s=0.95$ moves the normal peak to $570\mathrm{nm}$. More generally, ${\lambda }_B^{\prime }/{\lambda }_B=(n^{\prime }/n)s$ at normal incidence. A change of thickness alone is insufficient if processing also changes the refractive index.

Finite $N$ gives a finite spectral response. The nearest zeros satisfy $\delta =\pm 2\pi /N$, giving ${\lambda }_{\pm }={\lambda }_B/(1\pm 1/N)$. At the normal $600\mathrm{nm}$ peak, $N=30$ gives neighbouring zeros at approximately $580.65$ and $620.69\mathrm{nm}$; $N=10$ gives $545.45$ and $666.67\mathrm{nm}$. Thus greater sampled depth narrows the central band in this model. The small-detuning distance from peak to a nearest zero is approximately ${\lambda }_B/N$; it is not the full width at half maximum.

Real processed photographic emulsions have a continuous, irregular depth response rather than identical planes. Absorption, wavelength-dependent refractive index, unequal scattering strengths, a nonideal exposure reflector, changes during photographic processing and multiple scattering of light can alter both peak positions and shapes. A transfer calculation with measured material parameters would be needed for quantitative reflectance. The present model explains the phase-selection mechanism and the tradeoff between depth and spectral width; it cannot establish faithful reproduction of an arbitrary incident optical spectrum. Its useful lesson is that a photograph can record colour through sub-wavelength structure, while its eventual appearance still depends on illumination, viewing geometry and processing.

Consider one selected polarization of waves, represented by complex fields $E_1(t)$ and $E_2(t)$ at two sampling positions. Suppose their second moments are finite and stationary: $J_j=⟨|E_j(t){|}^2⟩$ and the following uncentered wave-field correlation do not depend on the absolute time $t$. DefineThe star denotes the complex conjugate; brackets denote an ensemble average. A long detector average estimates this quantity only if it adequately samples the fluctuations. The superscript $1$ identifies a field correlation: it does not mean the autocorrelation of the measured light intensity. The one-position version is $g^{(1)}(\tau )=⟨E^{\ast }(t)E(t+\tau )⟩/⟨|E{|}^2⟩$, so $g^{(1)}(0)=1$ for every nonzero field covered by this definition; this alone says nothing about its correlation at nonzero delay. ^[ref].

The Cauchy-Schwarz inequality gives $0\le |g_{12}^{(1)}(\tau )|\le 1$. At a specified delay, magnitude one means $E_2(t+\tau )=a{E}_1(t)$ in mean square for a fixed nonzero complex $a$: the samples have a reproducible relation in both wave amplitude and wave phase. Magnitude zero means the averaged field cross term vanishes. It need not mean statistical independence. A deterministic phase difference changes the argument of $g^{(1)}$, not its complex modulus; a distorted but fixed wavefront can therefore remain fully coherent. ^[ref].

To see the observable, recombine the samples in the same detected spatial mode with matched wave polarizations. With a common positive conversion factor $C$, let $I_j=C{J}_j$ be the single-arm light intensities reaching that detector. Adding $E_1(t)+e^{i\varphi }{E}_2(t+\tau )$ givesThis follows simply by expanding the squared complex modulus and averaging, then taking the maximum and minimum as $\varphi$ spans a full cycle. The extra wave phase scan must leave $I_1,I_2$ and the magnitude of the correlation effectively unchanged. Thus interferometric visibility equals $|g^{(1)}|$ only for equal nonzero light intensities and ideal mode matching. For $I_2/I_1=1/4$, perfectly coherent beams have interferometric visibility $0.8$; a measured value $0.4$ would imply $|g^{(1)}|=0.5$ under these assumptions. Orthogonal linear polarizations, background light, unresolved interference fringes or vibrations can reduce measured contrast further. One absent arm gives zero contrast and makes the normalized two-sample correlation undefined.

These characteristic wave amplitudes split a one-dimensional elastic wave into forward and backward components. Their instantaneous net elastic-wave energy flux is ${\varphi }_{+}^2-{\varphi }_{-}^2$. For real-part harmonic phasors the mean flux is $(|{\varphi }_{+}{|}^2-|{\varphi }_{-}{|}^2)/2$. In a static graded medium with real positive seismic impedance, put $g=(\log Z{)}^{\prime }/2$; then ${\partial }_x{\varphi }_{\pm }\pm {\beta }^{-1}{\partial }_t{\varphi }_{\pm }=g{\varphi }_{\mp }$. Seismic impedance gradients couple the two directions while preserving net flux.

Consider an ideal coherent monochromatic wave from an unresolved source on the optical axis, a uniform circular entrance pupil of diameter $D$, and ideal relay optics without wavefront errors. In the scalar Fraunhofer diffraction approximation, a Fourier transform carries the pupil's complex amplitude into the image plane. Choose dimensionless pupil coordinates $\mathbf{x}$ in units of $D$ and angular coordinates $\mathbf{u}=D\mathbf{\theta }/\lambda$, where $\lambda$ is the wavelength. DefineThe coronagraphic focal plane mask has transmission $M=1-W$, with $W=1$ inside the occulting disk and $W=0$ outside. Thus the field immediately after the mask is $MA$. Outside the disk $M=1$, so the local wave intensity is still $|A{|}^2$: hiding the core has left those rings intact. Reimaging an unlimited pupil giveswhere $\ast$ is convolution, $L$ is the Lyot stop transmission, and the proportionality constant in the irradiance is omitted. Coordinate reversals from physical relay optics are absorbed into the coordinate choice. These equations use wave amplitudes throughout propagation; applying a Fourier transform to wave intensity would give a different and incorrect model. This is the pupil/image calculation described by Sivaramakrishnan et al.; the circular disk and numerical choices below are our illustrative specialization. ^[ref].

The subtraction explains the beauty of the arrangement: in much of the interior, the original pupil field and the convolution term nearly cancel. Near a boundary, that cancellation is incomplete. The Lyot stop can therefore block a region containing disproportionately much unwanted light. A hard-edged coronagraphic focal plane mask has a transform with extended tails, so the interior need not be perfectly dark. With an unlimited relay, $L=1$ gives $A_{\mathrm{final}}=MA$ exactly: no further ring suppression occurs. A stop covering the full original optical pupil instead uses $L=P$ and already clips field outside the optical pupil; an undersized Lyot stop clips more of the remaining rim.

The toy uses a coronagraphic focal plane mask of angular radius $3\lambda /D$ and a Lyot stop of diameter $0.8D$. The discrete Fourier transform grid has $1024\times 1024$ samples, with 128 across $D$ and angular sampling $0.125\lambda /D$; it approximates the continuous equations on a finite periodic domain. The plot averages wave intensity in annuli of width $0.25\lambda /D$. Over the annulus $4\le |\mathbf{u}|\le 10$, the mean wave intensity divided by the direct stellar peak is approximately $1.40\times 10^{-4}$ without a coronagraph, the same with only the focal mask and an unlimited relay, $4.44\times 10^{-5}$ with the mask and full pupil, and $1.66\times 10^{-6}$ with the mask and $0.8D$ Lyot stop. Thus the undersized opening reduces this particular annular mean about 27-fold relative to the full-pupil comparison. It does not suppress every direction by that factor, and these values do not predict a real instrument's detection limit.

The ideal toy excludes a central obstruction, support vanes, scattering of light, wavefront errors, atmospheric effects, a spread of wavelengths, pointing drift and detector noise. It also contains only the star: no simulated planet or solar corona. For a far off-axis source whose image misses the coronagraphic focal plane mask, the continuous geometric coronagraph throughput of the $0.8D$ opening tends to $0.8^2=0.64$ relative to the full disk. Rejecting stellar light therefore costs some desired light as well. Enlarging the coronagraphic focal plane mask narrows the characteristic scale of its pupil-plane transform by the scaling property of the Fourier transform, but also hides a larger neighborhood of the source.

Frequency modulation carries a message by changing the instantaneous frequency of a carrier wave. In the ideal model the wave amplitude stays constant: information changes the rate at which the wave phase advances. Write a real electrical signal and its chosen complex representation asThe wave envelope in this representation has modulus $|z(t)|=A$ throughout, although $s(t)$ continually passes through zero. A smaller instantaneous value of $s(t)$ therefore does not mean a smaller wave amplitude. The continuously unwrapped wave phase $\varphi$ is measured in radians; its derivative is an angular frequency, while division by $2\pi$ gives frequency in cycles per second, or hertz.

For a message $m(t)$ expressed as an electric potential difference, an ideal voltage-controlled oscillator has sensitivity $K_f$ in hertz per volt and obeysThe integral is essential: putting the message directly into the wave phase gives phase modulation, whose instantaneous frequency instead depends on the message's derivative. Amplitude modulation changes the multiplying wave amplitude. These are different ways to encode the same kind of message signal.

An original slow illustration uses $f_c=20\mathrm{Hz}$, message frequency $f_m=1\mathrm{Hz}$, frequency deviation $\Delta f=5\mathrm{Hz}$, and FM modulation index $\beta =5$. With zero initial wave phase,Here $t$ is in seconds. The instantaneous frequency stays positive, between $15$ and $25\mathrm{Hz}$; the wave phase always advances. The oscillations bunch up where $f_i$ is large and spread out where it is small, while the wave envelope remains at $A$.

This separation is useful in a receiver: changing a positive multiplying wave amplitude does not change the zero crossings of $\cos \varphi$, so their timing still carries the message. Noise that moves those crossings, a vanishing wave amplitude, and frequency response distortion can still damage it. The surprising beauty of FM is that a constant wave envelope can conceal a rich frequency spectrum: the components interfere to change timing instead of the envelope's height.

Fresnel diffraction describes diffraction when the quadratic propagation wave phase across an aperture must be retained. Its Fresnel propagator becomes the Fraunhofer diffraction description when that curvature is negligible.

A Fresnel diffraction pattern is a map of how propagation changes the wave interference of light from an aperture. It can resemble the opening, develop bright and dark structure near its edges, and eventually approach a Fraunhofer diffraction pattern. The controlling comparison involves wavelength, opening size and propagation distance together; a distance large compared with the opening need not give the far field.

For a monochromatic wave, use time dependence $e^{-i\omega t}$, wavenumber $k=2\pi /\lambda$ and complex amplitude $U$. Let $U_0(\xi ,\eta )$ be the transmitted field immediately behind an aperture in the plane $z=0$. In a homogeneous medium, the scalar paraxial approximation givesThis is the two-transverse-coordinate Fresnel propagator, with its carrier $e^{ikz}$ restored. It adds complex amplitudes before taking their squared complex modulus to obtain light intensity. It is a scalar paraxial approximation, rather than a full vector solution of Maxwell equations.^[ref]

The source in this constructed geometry supplies a plane wave; the observation screen is a plane perpendicular to the optical axis. The right panel isolates the aperture-dependent wave phase that distinguishes Fresnel diffraction from Fraunhofer diffraction. Its vertical scale is in radians, and the drawing's lengths are schematic.

Keeping this wave phase is useful when transporting a field between nearby optical components: both the shape and the complex amplitude matter for subsequent wave interference. A particularly striking feature is that propagation can create a dark centre even from a uniformly illuminated symmetric aperture; the incident wave amplitudes need not be unequal for their outgoing contributions to cancel.

A fundamental Gaussian beam is an ideal beam of light whose transverse wave amplitude has a Gaussian function profile. Here consider a circularly symmetric beam in vacuum, with propagation along $z$, two transverse coordinates $x,y$, and radius $r=\sqrt{x^2+y^2}$ from the optical axis. With the beam waist at $z=0$, total power $P>0$, and vacuum wavelength $\lambda$, its time-averaged light intensity, meaning irradiance in $\mathrm{W}{\mathrm{m}}^{-2}$, isHere $w_0$ is the minimum beam radius, and $z_R$ is the Rayleigh range. Far from the beam waist, $w(z)\sim \theta |z|$ with beam divergence half-angle $\theta =\lambda /(\pi w_0)$. These are the fundamental circular Gaussian beam laws within the scalar paraxial approximation, appropriate when $w_0\gg \lambda$ and the beam divergence is small. The underlying envelope solves the parabolic wave equation exactly; this does not make it an exact vector solution of Maxwell equations. ^[ref].

The radius convention matters. Writing $U$ for the complex wave amplitude envelope, the same $w(z)$ is the wave amplitude $1/e$ radius and the light intensity $1/e^2$ radius:At $r=w(z)$ these ratios are $e^{-1}\approx 0.367879$ and $e^{-2}\approx 0.135335$, respectively. The light intensity $1/e$ radius is instead $w(z)/\sqrt{2}$. None of these radii is a sharp edge: the ideal Gaussian function has nonzero tails at every finite radius. Kogelnik and Li explicitly use the wave amplitude $1/e$ radius in their fundamental-mode treatment. ^[ref].

For a numerical example, take vacuum wavelength $\lambda =632.8\mathrm{nm}$ and beam waist radius $w_0=0.500\mathrm{mm}$. Substitution gives Rayleigh range $z_R\approx 1.24114754\mathrm{m}$ and far-field beam divergence half-angle $\theta \approx 0.402852992\mathrm{mrad}$. At $z=z_R$, the beam radius is $w=\sqrt{2}{w}_0\approx 0.707107\mathrm{mm}$ and the on-axis irradiance is half its waist value. At the same wavelength, halving the beam waist radius to $0.250\mathrm{mm}$ quarters the Rayleigh range to $0.310286885\mathrm{m}$ and doubles the beam divergence half-angle to $0.805705984\mathrm{mrad}$. At fixed power, that smaller beam waist also quadruples the peak irradiance.

The contours in the left panel are $r=\pm w(z)$ in a plane through the optical axis; they mark an light intensity level, not an outer boundary. The right panel compares the exact model radii with the far-field asymptotes $w=\theta z$. The horizontal and vertical axes use different length units, so the drawn slopes are not physical angles.

The normalization expresses conserved power during free propagation. In transverse polar coordinates, an annulus has area $2\pi rdr$, soA proof outline is the substitution $s=2r^2/w(z{)}^2$ in this integral. Spreading reduces peak irradiance as $w^{-2}$ while increasing transverse area as $w^2$, leaving power unchanged. In particular, the disk of radius $w(z)$ contains $1-e^{-2}\approx 0.8646647168$, or $86.4665\%$, of the total power; the local light intensity at its rim is only $13.5335\%$ of the on-axis light intensity. These are different quantities.

The left panel compares normalized wave amplitude magnitude and irradiance; the right panel integrates the transverse area to show enclosed power. For an ideal absorbing circular aperture of radius $a$ centered on the optical axis, the same fraction gives the power passing the aperture at that plane. An aperture of radius $w$ clips about $13.5335\%$ of the incident power; radius $2w$ clips only $e^{-8}\approx 0.000335463$, or $0.0335463\%$. After clipping, extra diffraction means the transmitted field no longer obeys the unmodified ideal Gaussian beam law. Every finite aperture cuts some ideal tail. Kogelnik and Li distinguish their large-aperture beam analysis from their finite-aperture diffraction analysis. ^[ref].

The useful tradeoff is that a small bright spot cannot remain narrow indefinitely. An optical lens can create a smaller beam waist, but diffraction then makes the beam spread faster. The particularly simple relation $w_0\theta =\lambda /\pi$ ties localization directly to beam divergence. A proof outline for the propagation law starts with $U(x,y,0)=U_0{e}^{-(x^2+y^2)/w_0^2}$. The parabolic wave equation separates into two transverse factors, and a Gaussian integral in each coordinate givesHere $q$ is dimensionless and $k$ is the wavenumber; the carrier convention is $e^{i(kz-\omega t)}$, with angular frequency $\omega$. Taking the squared complex modulus gives the stated $w(z)$ and its $w^{-2}$ light intensity scaling. Each transverse factor contributes $q^{-1/2}$, explaining the $q^{-1}$ prefactor here and the different prefactor in the preserved Gaussian beam with one transverse coordinate.

Historically, Gaussian beam theory developed alongside early laser resonator theory and lens-guided transmission. In their 1966 review, Kogelnik and Li credit Boyd and Gordon, and Goubau and Schwering, with early beam treatments; they summarize earlier work rather than claim the whole construction as a new discovery. Their paper brings beam propagation, focusing and finite-aperture diffraction into a common practical account. ^[ref].

Change in wave amplitude as the same propagating energy is distributed across changing wavefront area. It is distinct from intrinsic material absorption even though both can reduce observed amplitude with distance.

A gravity-capillary wave has restoring contributions from both gravity and surface tension. Consider an initially resting, horizontally unbounded liquid layer of constant depth $h$, with uniform mass density $\rho$, incompressible flow, inviscid flow and irrotational flow. The upper gas exerts constant pressure and has negligible inertia. Assume a clean interface with constant surface tension $\gamma$: no surfactant stresses, imposed current or wind. The wave amplitude $a$ must satisfy $ka\ll 1$ and $a\ll h$, so linearized boundary conditions apply. Deep water means $kh\gg 1$, rather than a particular depth in metres. ^[ref]

Let $z=0$ be the undisturbed free surface, with liquid below it. Potential flow gives a velocity potential $\varphi$ with $\mathbf{u}=\nabla \varphi$, and incompressible flow gives the Laplace equation. The linearized kinematic boundary condition for a free-surface graph and the Unsteady Bernoulli equation, including gravity and the Young–Laplace equation, giveHere $\eta (x,t)$ is the interface displacement; the signed small-slope pressure jump is $p_{\mathrm{liquid}}-p_{\mathrm{gas}}=-\gamma {\eta }_{xx}$. A wave crest has ${\eta }_{xx}<0$, so its capillary pressure is positive. Choose a right-travelling Fourier mode of positive wavenumber $k$ and angular frequency $\omega$:The factor $e^{kz}$ decays into the deep liquid, satisfies the Laplace equation, and makes ${\varphi }_z={\eta }_t$ at $z=0$. Substituting in the dynamic boundary condition leavesThis derivation displays the two restoring contributions with the same sign. Dimensional analysis also checks the result: $\gamma /\rho$ has units ${\mathrm{m}}^3{\mathrm{s}}^{-2}$, so both $gk$ and $(\gamma /\rho )k^3$ have units ${\mathrm{s}}^{-2}$. ^[ref]

At finite depth, a flat impermeable bottom at $z=-h$ requires ${\varphi }_z=0$. Replacing $e^{kz}$ by the normalized vertical factor $\cosh (k(z+h))/\cosh (kh)$ makes the surface ratio ${\varphi }_z/\varphi =k\tanh (kh)$. Eliminating the mode amplitudes then givesThus the same model supplies its own depth check: use the deep-water formula only when $\tanh (kh)$ is sufficiently close to one for the required accuracy. For example, $kh=3$ reduces $\omega$ by about $0.248\%$ relative to deep water, while $kh=1$ reduces it by about $12.7\%$. These percentages are calculations from this factor. ^[ref]

In deep water the phase velocity magnitude follows a wave crest, while the group velocity magnitude follows the wave envelope of a narrow wave packet. Differentiating the positive-frequency dispersion relation givesFor $k{\ell }_c\ll 1$, the surface gravity wave limit has $c_g\simeq c_p/2$. For $k{\ell }_c\gg 1$, the deep-water capillary-wave dispersion relation has $c_g\simeq 3c_p/2$: a wave crest slips backward relative to the faster wave envelope, even while both advance. Neither speed is the oscillating fluid element's own velocity. ^[ref]

There is a particularly neat coincidence at the crossover. Differentiate $c_p^2=g/k+\gamma k/\rho$: its unique minimum is at $k_c=\sqrt{\rho g/\gamma }=1/{\ell }_c$, exactly where the restoring terms are equal. ThereforeFor the stipulated parameters in capillary length, the minimum is $0.231\mathrm{m}{\mathrm{s}}^{-1}$. At $\lambda =3\mathrm{mm}$, direct substitution gives $f=\omega /(2\pi )\simeq 131\mathrm{Hz}$, $c_p\simeq 0.394\mathrm{m}{\mathrm{s}}^{-1}$ and $c_g\simeq 0.580\mathrm{m}{\mathrm{s}}^{-1}$. This short wavelength is capillary dominated, but its group velocity is still slightly below the pure-capillary ratio $3c_p/2$ because gravity has not vanished.

The group delay is the frequency derivative of an unwrapped transfer wave phase, using the $e^{-i\omega t}$ convention. It predicts the leading delay of a narrow wave envelope when the transfer wave amplitude and the delay vary little across its band; it is not generally a signal-front arrival time.

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 solid-to-fluid interface can turn spherical source P waves into parallel vertical rays if its height $Z(R)$ satisfies the displayed travel-time identity, with $r=\sqrt{R^2+Z^2}$. The positive branch through $Z(0)=z_0$ is a hyperboloid when $\alpha >{\alpha }_f$. Differentiating the travel-time identity enforces Snell law for elastic and acoustic waves. Travel time above the interface is then $z_0/\alpha +(z-z_0)/{\alpha }_f$, independent of cylindrical radius. The leading wavefront is planar although its wave amplitude need not be uniform across that plane.

At a stationary, uniform planar boundary, incident and transmitted plane waves have the same variation of wave phase along the interface and in time. The figure marks places where equal-wave phase wavefronts meet the boundary. Although their perpendicular spacing differs, their intersections must advance along the boundary together. This is the local mechanism behind Snell's law.^[ref]

For a short derivation, let $x$ run along the interface and $z$ increase into the transmitted medium. Write each wave phase as ${\Phi }_j=k_{j,x}x+k_{j,z}z-{\omega }_{j}t+{\varphi }_j$. The electromagnetic boundary conditions must hold at every $x,t$ on $z=0$; their matching requires ${\omega }_i={\omega }_t$ and $k_{i,x}=k_{t,x}$. Constant offsets ${\varphi }_j$ and wave amplitudes are determined separately. Since $k_{j,x}=k_j\sin {\theta }_j$, this givesFor nonzero tangential wavenumber, the spacing of equal-wave phase intersections along the interface is ${\lambda }_i/\sin {\theta }_i={\lambda }_t/\sin {\theta }_t$. It is not the perpendicular wavelength. At normal incidence, the tangential wavenumber vanishes and the interface is uniform in wave phase.^[ref]