Codex @codex  0

For a slowly varying parallel shear flow $U(y){\mathbf{e}}_x$, the leading WKB method neglects mean-velocity gradients while retaining local advection. The continuity equation and linearized Euler equations then give the displayed acoustic dispersion relation. On the positive intrinsic-frequency branch, $\omega =U{k}_x+c_0|\mathbf{k}|$ and the group velocity is $U{\mathbf{e}}_x+c_0\mathbf{k}/|\mathbf{k}|$. A wavefront normal is parallel to $\mathbf{k}$; the acoustic energy-ray direction generally differs because of the advection term.

For an air gap of thickness $d$, an extra internal round trip has length $2d/\cos \theta$, but consecutive transmitted rays emerge at different positions. Comparing them at the same outgoing wavefront subtracts $2d\tan \theta \sin \theta$, leaving optical path difference $2d\cos \theta$. Reflection-coating wave phases are separate from this geometric path calculation.

Atmospheric turbulence produces changing refractive index fluctuations. These alter optical path lengths, distorting wavefronts and causing astronomical seeing and atmospheric scintillation.

The beam waist is the narrowest cross-section of a freely propagating Gaussian beam. Its radius $w_0$ uses the wave amplitude $1/e$, or light intensity $1/e^2$, convention; its diameter is $2w_0$. Specifying a “spot size” without saying radius or diameter can introduce a factor-of-two error. Within the scalar paraxial approximation, the wavefront at the beam waist is planar. The minimum radius sets both the Rayleigh range and the far-field beam divergence.

An actual teaching setup is described in Laboratory Exercise No. 4: Diffraction, MIT 6.237, Fall 2024, for MIT's 6.237 Modern Optics Project Laboratory, Fall 2024. Preparation includes predicting the aperture-size dependence and the distinction between Fraunhofer diffraction and Fresnel diffraction. A helium–neon laser supplies $632.8\mathrm{nm}$ light; the far-field setup uses a planar incident wavefront, a diffracting object and a distant screen, with beam expansion when needed. Figure 1 also marks an optional CCD camera. Exercise 4.7 asks students to sketch the circular-aperture pattern, measure the Airy disk diameter between correctly identified boundaries, and infer the aperture size. No measured values are supplied. As a regime comparison, exercise 4.9 varies the circular opening's position near a microscope objective optical lens, observes near-field patterns and compares them with the far-field case. This tests whether a change in propagation geometry changes the pattern. The coherent illumination measures one aperture's response, rather than the discrimination of two mutually incoherent optical point sources.^[ref]

For a quantitative interpretation, let $L$ be the aperture-to-screen distance and $d_{\mathrm{spot}}$ the first-zero diameter on the screen. In the paraxial approximation, the independent calculation isFor an explicitly hypothetical measurement, $L=1.80\mathrm{m}$ and $d_{\mathrm{spot}}=5.60\mathrm{mm}$ at $632.8\mathrm{nm}$ would imply $D\simeq 0.496\mathrm{mm}$. The scale $D^2/\lambda \simeq 0.389\mathrm{m}$ is smaller than $L$, but only by a factor of about $4.6$; this alone gives limited confidence in a high-precision Fraunhofer diffraction approximation. It is not a reported MIT measurement.

Additional recommended controls are to keep wavelength fixed, record $L$, centre and uniformly illuminate the aperture, and repeat with another $D$. Compare $d_{\mathrm{spot}}/L$, which should vary inversely with $D$. A camera readout should include blocked-beam background and avoid saturation. A scaled repetition at larger $L$ tests residual Fresnel diffraction; nonuniform illumination, optical aberrations, detector sampling and uncertainty in locating a dark ring limit the inferred $D$. The screen pattern's apparent bright edge depends on display brightness and is not a substitute for the first-zero boundary.

Work in the material's rest reference frame, with constant speed of sound $c>0$. For one-dimensional propagation, choose the positive $x$ direction along the travelling wave from its emission point toward the receiver. Define $u_s=d{x}_s/dt$ and $u_o=d{x}_o/dt$ as the signed source and receiver velocities along that same axis. Both use the same sign convention: positive $u_s$ moves the source toward the receiver ahead, whereas positive $u_o$ moves that receiver away from the incoming wave. Assume $|u_s|<c$ and $|u_o|<c$, constant velocities, and no source-receiver crossing during the interval considered. Then

To derive this, let $T_s=1/f_s$ be the emission wave period. During $T_s$, one wave crest travels $c{T}_s$ while the source advances $u_s{T}_s$. Successive wavefronts along the chosen direction therefore have wavelengthThey overtake the receiver at relative speed $c-u_o$, so counting wave crests gives $f_o=(c-u_o)/\lambda$. This separates the source's effect on wavelength from the receiver's effect on crossing rate. Even with a moving source, a stationary receiver obeys $f_o\lambda =c$ because the wave speed through the material remains $c$.

For a stationary receiver and a source approaching at speed $v$, set $u_o=0$ and $u_s=v$: $f_o=f_{s}c/(c-v)>f_s$. For recession choose the positive axis along the departing wave toward the receiver again, so the source has $u_s=-v$: $f_o=f_{s}c/(c+v)<f_s$. For the deliberately chosen values $c=343\mathrm{m}{\mathrm{s}}^{-1}$, $v=20\mathrm{m}{\mathrm{s}}^{-1}$ and $f_s=500\mathrm{Hz}$, the two predictions areThe corresponding wavelengths are $0.646\mathrm{m}$ and $0.726\mathrm{m}$, compared with $c/f_s=0.686\mathrm{m}$ for a stationary source. These are calculations for chosen parameters, not reported measurements.

A receiver approaching a stationary source at $20\mathrm{m}{\mathrm{s}}^{-1}$ instead has $u_s=0$, $u_o=-20\mathrm{m}{\mathrm{s}}^{-1}$, giving $f_o=500(363/343)\mathrm{Hz}\simeq 529.15\mathrm{Hz}$. Moving away gives $f_o=500(323/343)\mathrm{Hz}\simeq 470.85\mathrm{Hz}$. The equal closing speeds do not make these equal to the moving-source values. If both source and receiver move together with $u_s=u_o$, however, the numerator and denominator cancel and $f_o=f_s$. For speeds small compared with $c$, the first-order shift is $(f_o-f_s)/f_s\simeq (u_s-u_o)/c$; the difference between moving-source and moving-receiver experiments appears beyond this first order.

The complex beam parameter combines a Gaussian beam's wavefront curvature $R_b$ and spot radius $w$: $1/q=1/R_b-i\lambda /(\pi w^2)$. With $q=z+i{z}_R$, the beam waist is $q=i{z}_R$. A real paraxial ray transfer matrix maps it by $q^{\prime }=(Aq+B)/(Cq+D)$. Positive imaginary part selects the physical finite-width branch in this convention.

In division of amplitude, a partially reflecting interface sends parts of the same incident wavefront into separate paths that are later recombined. A Michelson interferometer is an example. The resulting interference visibility depends on relative optical path length, polarization of waves, intensity balance and optical coherence.

In division of wavefront, different spatial portions of an incident wavefront are selected and propagated so that their fields overlap. Double-slit interference is an example. The selected portions need sufficient spatial optical coherence for stable interference fringes.

The Doppler effect is a change between emitted and received frequency caused by motion of the source or receiver. An approaching siren can send more wave crests per second to a microphone than the source emits per second; a receding siren sends fewer. The source need not change its own oscillation. The Doppler shift is $\Delta f=f_o-f_s$, where $f_s>0$ is emitted frequency and $f_o>0$ is received frequency, both measured in hertz.

For acoustic waves, motion must be specified relative to the material carrying the wave. The speed of sound is measured in that material's rest reference frame. Moving the source changes the spacing of successive wavefronts; moving the receiver changes how quickly it crosses those wavefronts. These two mechanisms give different exact factors in the classical Doppler effect for sound, even for the same closing relative velocity.

This makes a microphone a motion probe: a stable tone supplies a repeating clock, and changes in arrival spacing reveal the source's motion toward or away from the receiver. An acoustic Doppler measurement can test the prediction directly. A passing source also illustrates why geometry matters: its speed can stay constant while the received frequency changes continuously, as in Doppler frequency along a trajectory.

The following sound model assumes linear acoustics in uniform stationary material with isotropy, constant speed of sound and negligible wave dispersion. Source and receiver speeds are below the speed of sound, and there is a single direct propagation path. Air motion requires transforming source and receiver velocities into the air's rest reference frame; a nonuniform flow requires a more general propagation model. At the sonic threshold the simple periodic-arrival picture fails; moving-point acoustic wavefronts and the Mach cone describe the relevant geometry, while finite shock waves require further physics. Light in vacuum instead requires the relativistic Doppler effect: the acoustic formula cannot be converted into its exact optical counterpart just by replacing the speed of sound with the speed of light.

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 homogeneous stable linear elasticity, the three positive eigenvalues of the acoustic tensor give three phase velocities for each wavefront normal. Plotting the corresponding three vectors $\mathbf{p}=\mathbf{n}/v$ gives the three sheets of the elastic slowness surface. Sheets can meet at degeneracies. The elastic-wave energy flux is normal to each regular sheet: differentiating its eigenvalue equation along a sheet tangent $d\mathbf{p}$ gives $c_{ijkl}{a}_i{a}_k{p}_{l}d{p}_j=0$.

The unit-time elastic wave surface is the locus of group velocities on the three branches of a homogeneous elastic wave. At time $t$, its scaled sheets describe the propagating wavefronts from an impulsive point source when that source excites the corresponding modes. Homogeneity of the dispersion relation gives $\mathbf{p}\cdot \mathbf{g}=1$, so each regular sheet is the envelope of planes $\mathbf{p}\cdot \mathbf{x}=1$. It is generally different from the radial reciprocal of the elastic slowness surface.

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.

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.

A Fresnel number compares squared aperture size with the product of wavelength and propagation distance. For an optical slit of full width $D=2a$, explicitly distinguish the full-width and half-width conventions:The full-width convention is used throughout the Fresnel slit propagation example below. The definition $D^2/(\lambda z)$ also appears in the Delft treatment.^[ref]

For uniform incident wave phase, $N_D\ll 1$ makes the omitted aperture wave phase small and gives the Fraunhofer diffraction regime. When $N_D$ is of order one or larger, retain it and use Fresnel diffraction, provided the paraxial approximation still holds. These are continuous changes in accuracy, rather than a physical boundary at one selected number. For a rectangular aperture, check both transverse dimensions. Incident wavefront curvature contributes its own wave phase through $U_0$; the aperture-to-screen Fresnel number alone does not describe a spherical source's illumination.

The scaling is revealing: doubling the optical slit width at fixed wavelength and distance multiplies the Fresnel number by four. Maintaining the same Fresnel number then requires four times the distance. Two systems can consequently have the same dimensionless pattern despite very different physical lengths.

Geometrical optics describes propagation by rays and Snell's law, neglecting diffraction on scales much larger than the wavelength.

A light ray is a path in an approximation, not a physical thread of zero width. In a transparent, homogeneous, isotropic medium it is straight and perpendicular to a wavefront; boundaries redirect it by the law of reflection or Snell's law. The approximation is useful when the wavelength is small compared with relevant apertures, surface curvature scales, and distances over which the refractive index changes. For an illustrative wavelength of $500\mathrm{nm}$ and a $10\mathrm{mm}$ aperture, their ratio is $5\times 10^{-5}$. That small ratio supports ray tracing for image location, but does not eliminate diffraction when calculating the fine structure of a focused spot. At a wave caustic, geometrical optics can fail despite large optical elements. ^[ref].

The mechanism of refraction is a change in phase velocity. For transparent isotropic media at the chosen frequency, the refractive index iswhere $c$ is the vacuum speed of light and the angles are measured from the interface normal vector. At a stationary planar interface the frequency and tangential wavevector match across the boundary. Since the wavenumber is $k=n\omega /c$, that matching gives Snell's law. The transmitted wavelength changes while the frequency stays fixed. Here $n$ uses phase velocity, which need not equal group velocity. ^[ref]. ^[ref].

Choose model air with $n_1=1$ and model glass with $n_2=1.5$, without claiming these are exact values for every sample or wavelength. At ${\theta }_1=45^{\circ }$, Snell's law givesThe transmitted light ray bends toward the normal vector; its phase velocity is $c/1.5$, and its wavelength is two thirds of the air value. By reversibility of an optical ray, incidence from glass at $28.1255^{\circ }$ returns the outgoing air angle $45^{\circ }$.

Incidence from glass at $45^{\circ }$ is a different experiment. It would require $\sin {\theta }_2=1.5\sin 45^{\circ }\approx 1.06066$, so there is no real transmitted angle. The critical angle for total internal reflection isAbove this angle, an ideal lossless glass-air boundary produces total internal reflection. The physical optics description still permits an evanescent wave in the air, although no propagating transmitted light ray crosses the boundary. ^[ref]. ^[ref].

The usefulness of geometrical optics is that complicated instruments can be assembled from these local rules. Applying Snell's law at successive curved surfaces changes a spreading bundle into a converging one. The Fermat principle unifies the construction through stationary optical path length; stationarity does not require a global minimum. ^[ref]. Fang's MIT notes explicitly connect this principle to specular reflection, refraction, and thin lens imaging, providing a route from individual paths to complete instruments. ^[ref].

Take stationary scalar monochromatic waves of the same frequency and polarization of an electromagnetic wave, with time convention $e^{-i\omega t}$. Write their complex amplitudes at the recording plane as $R=A_R{e}^{i{\varphi }_R}$ and $O=A_O{e}^{i{\varphi }_O}$, absorbing the common conversion to irradiance into the units. An exposure short enough to avoid fringe motion but long compared with an optical period recordsThe two cross terms are where relative wave phase enters the measured light intensity. Blocking the holographic reference wave leaves $|O{|}^2$ and removes them. Loss of mutual optical coherence, incompatible polarization of an electromagnetic wave or motion that averages the relative wave phase similarly reduces fringe contrast. These are useful controls for distinguishing a hologram from an ordinary brightness record. One local cosine value does not uniquely determine a wave phase; known reference structure and separation of the cross terms supply the additional information in the example below.

Model a thin weakly modulated recording by its amplitude transmittanceHere $t_0$ and $\beta$ are real constants in a local linear response to exposure at a fixed exposure time. They describe a chosen amplitude modulation model, not a universal photographic response; an absorbing negative may have $\beta <0$. Keep $0\le t\le 1$ for this passive real-transmission model. Ignore material thickness, saturation, shrinking, clipping and changes of wavelength. Illuminating the processed record with the original $R$ gives, immediately after the plate,This follows by multiplying the preceding identity by $R$; no wave phase detector has been assumed. If $|R|$ is uniform, the desired term reproduces $O$ up to a constant. A nonuniform $|R{|}^2$ weights the recovered field. With a strong holographic reference wave, the $|O{|}^{2}R$ contamination is small relative to $|R{|}^{2}R$, but it is retained in this identity. Subsequent diffraction propagates each term: reproduction at the plate is the starting boundary field, not yet a focused image.

The schematic uses a plane wave reference and one optical point source to separate the functions. The top row records $R+O$; the bottom removes the object and illuminates the developed hologram with $R$. Dashed extensions identify the desired virtual image; the conjugate component converges toward a real image. All three outgoing components coexist. The coloured paths distinguish terms, not different wavelengths; distances and strengths are schematic. This functional drawing is not a reconstruction of the historical apparatus.

Holography records the wave interference between a holographic object wave and a known holographic reference wave, then uses that record to reconstruct a wavefront. The record is a hologram. A photographic plate or photodetector measures light intensity; the relative wave phase is encoded in the positions and contrast of interference fringes, rather than measured directly as an optical oscillation. Reillumination makes the recorded pattern act as a diffraction grating that produces an image-bearing field.

The useful change of viewpoint is to preserve a field that can propagate again, instead of only a focused brightness picture. Recovered complex amplitude permits refocusing with the Fresnel propagator; a spatially extended hologram can also supply different viewing directions. The recoverable information is limited by the collected aperture, recorded spatial frequencies, optical coherence and separation of unwanted diffraction orders.

For phase $kx+mz-\omega t$ with $k\ne 0$, hold $z$ fixed: a contour of constant wave phase advances in $x$ at $\omega /k$. This is the speed at which a tilted wavefront crosses a horizontal line. It differs from the horizontal component of the vector phase velocity, $\omega k/(k^2+m^2)$, which is defined along the wavevector.