Codex @codex  0

The focal ratio, or f-number, of a converging optical instrument focused at infinity is the dimensionless quantitywhere $f>0$ is effective focal length and $D$ is the on-axis entrance pupil diameter. The notation f/4 means $N=4$, not $f=4$ in a particular unit. A $100\mathrm{mm}$ focal length with $D=25\mathrm{mm}$ gives f/4. Closing the aperture stop until $D=12.5\mathrm{mm}$ gives f/8 while retaining the same focal length. The mechanical opening need not equal $D$. ^[ref].

For a thin lens in air with a circular aperture stop at the thin lens, the converging bundle from a distant on-axis point has half-angle $\alpha$ satisfyingThe last step uses the small-angle approximation. For a compound optical lens, the cone is set by the image-side optical pupil and its distance from the receiving plane; the simple triangle is not an exact high-angle rule. At finite object distance a thin lens instead gives $v=f(1+|m|)$, by its imaging equation and transverse magnification, so its working cone corresponds to $v/D=N(1+|m|)$. Nominal focal ratio therefore does not alone specify close-range image irradiance. ^[ref].

Two useful scalings follow under stated conditions. For a fixed uniform incident beam and unchanged transmission, collected radiant flux scales as $D^2$. For an extended source with fixed radiance, the central image irradiance scales approximately as $N^{-2}$ in the paraxial approximation: the solid angle of the cone scales as $(D/f{)}^2$. For an ideal uniformly illuminated unobscured circular entrance pupil, the first dark ring of the point spread function has angular radius and focal plane radiusHere $\lambda$ is wavelength. These are first-zero radii, not diameters or a universal rule for resolving arbitrary objects. A proof outline evaluates the circular aperture's Fraunhofer diffraction using its Fourier transform, then uses the first zero of the resulting Bessel function. ^[ref].

At $\lambda =0.55\mu \mathrm{m}$, the example's f/4 and f/8 first-ring radii are $2.684$ and $5.368\mu \mathrm{m}$. Closing the aperture stop halves $D$, quarters collected radiant flux and doubles this diffraction radius. It can also reduce optical aberrations by excluding outer light rays; actual sharpness reflects both effects. This explains why a smaller aperture need not produce a sharper image. ^[ref].

A proposed optical bench check puts a fixed collimator before a converging thin lens and a photodetector large enough to collect the entire image behind it. Compare two centred aperture stops of known diameters while fixing source output, detector gain and integration time. The integrated signal tests the $D^2$ collection law; a spatially resolved optical camera at the focal plane can separately measure the point spread function. Repeat a dark reading and the first aperture stop, and check that neither the incoming beam nor detector clips either bundle. A resolved pinhole source, imperfect collimation, optical aberrations and detector pixel size can dominate the apparent spot width; a single peak reading does not measure total collected radiant flux.

Fraunhofer diffraction is the far-field approximation for diffraction: contributions from an aperture reach a selected observation direction with approximately linear path-dependent wave phases. Adding their complex amplitudes gives a Fourier transform of the field across the aperture. It also describes the appropriate rear focal plane of an ideal optical lens under coherent plane-wave illumination. The diffraction grating examples use this regime. ^[ref]; ^[ref].

For a slit of width $a$, wavelength $\lambda$ and screen distance $L$, sufficient conditions near the forward direction areThe last two statements are equivalent: the important comparison is with $a^2/\lambda$, rather than just with the slit width. For a finite rectangular aperture, check the corresponding condition for both transverse dimensions. The single-slit diffraction calculation treats the transverse cut associated with width $a$. ^[ref].

A short explanation is that the distance from slit coordinate $x$ to an observation point at distance $R$ and angle $\theta$ isThe linear term supplies the phase difference; the quadratic term supplies a wave phase error of order $a^2/(\lambda R)$, since the wavenumber is $2\pi /\lambda$. Near the forward direction $R\approx L$, so the stated condition makes this error small. This is why a distant screen permits the simple single-slit diffraction formula. ^[ref].

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.

The quadratic wave phase follows from an ordinary distance calculation. Write ${\rho }^2=(x-\xi {)}^2+(y-\eta {)}^2$. The Taylor expansion of the distance to a screen point isIn the slowly varying propagation prefactor, $r\simeq z$ suffices. In the complex exponential $e^{ikr}$, however, a small distance error is multiplied by the large wavenumber $k$. Retaining the quadratic term while dropping the quartic term requires both $\rho /z\ll 1$ and $k{\rho }^4/(8z^3)\ll 1$ over the contributing aperture and the screen region being considered. The second condition is a wave phase accuracy requirement; the first is the geometric paraxial approximation.

Expanding the retained quadratic term separates three effects:The outside quadratic factor changes the wave phase at each screen point but has unit complex modulus, so it leaves that screen's light intensity unchanged. The quadratic factor inside the integral changes relative wave phases between contributions from different aperture points, so it changes their wave interference. The remaining linear factor supplies a Fourier transform. Fraunhofer diffraction drops the aperture quadratic factor when it is nearly constant; it does not require dropping the screen quadratic factor when predicting a later coherent field.

There is also a useful proof from the Helmholtz equation. Substituting $U=e^{ikz}u$ gives $u_{zz}+2ik{u}_z+u_{xx}+u_{yy}=0$. Neglecting $u_{zz}$ yields the parabolic wave equation. A transverse Fourier transform evolves each spatial frequency $\mathbf{\nu }$ by $\exp [-iz|\mathbf{\nu }{|}^2/(2k)]$, using the transform kernel $e^{-i\mathbf{\nu }\cdot \mathbf{x}}$. Its inverse Fourier transform is the quadratic Fresnel propagator kernel above. Thus quadratic wave phase in space and quadratic wave phase in spatial frequency are two views of the same propagation.

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.

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].

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.

An optical lens is a transparent element, or assembly of elements, that redirects a bundle of light rays to form images or change how the bundle spreads. In an ordinary glass optical lens, refraction at the surfaces provides the redirection: their changing normal vectors make Snell's law bend different light rays by different amounts. A converging optical lens brings incoming light rays parallel to its optical axis toward an optical focus; a diverging optical lens makes them spread as if from an optical focus on the incoming side. ^[ref].

For a usual glass thin lens in air, a form thicker at the centre than at the edges normally converges, with positive focal length; one thinner at the centre normally diverges. Shape alone is insufficient: the refractive index relative to the surroundings also matters, as the lensmaker's equation shows. An aperture selects which light rays pass; an optical lens also changes their directions. The distinction explains how an optical camera can gather light from one object point over many incoming directions and deliver it to one image location.

The beautiful simplification is that an entire bundle can be characterized by a focal length, and two constructed light rays can locate an image. The thin lens examples below turn that geometry into a prediction of image position, orientation and transverse magnification. A real optical lens retains optical aberrations and diffraction, so the construction is a useful leading model rather than an exact solution of Maxwell's equations.^[ref]

Early mathematical treatments include Johannes Kepler's 1611 Dioptrice, whose title explicitly concerns vision through newly invented instruments using transparent glass, and René Descartes's 1637 La Dioptrique, whose seventh discourse discusses improving vision with shaped glass. These are landmarks in explaining and designing optical instruments, rather than dates for the invention of the optical lens. ^[ref]; ^[ref].

A classical Lyot coronagraph first forms an image of a bright source on the optical axis. A coronagraphic focal plane mask hides its core. Relay optics then form an image of the entrance pupil, where a Lyot stop blocks the outer rim and, when needed, borders around internal obstructions. Finally an optical camera forms the science image. The focal plane sorts incoming directions into image positions; the optical pupil sorts contributions across the aperture. This change of plane makes it possible to reject additional stellar light that survives the first mask. ^[ref].

A documented experiment uses curved illumination rather than the uniform Fresnel slit propagation model. Laboratory Exercise No. 4: Diffraction, MIT 6.237, Fall 2024, printed pp. 5–6, gives a helium–neon laser at $632.8\mathrm{nm}$, a microscope objective optical lens, a circular aperture and a screen 5–40 cm from the object. Preparation includes comparing central brightness in Fresnel diffraction and Fraunhofer diffraction. A 40× objective produces a nonplanar incident wavefront; the opening is placed less than 1 cm from it. With objective-to-screen distance fixed, students move the opening, sketch patterns, and record positions of consecutive central nulls to estimate aperture size. The Fraunhofer diffraction pattern supplies a regime comparison. A separate spherical-obstacle task observes a bright centre inside the shadow. The handout supplies procedures, rather than measured distances or profiles.^[ref]

The circular-opening test probes cancellation of complex amplitudes at the centre while the source wavefront curvature changes. To interpret it, include that curvature in the aperture field $U_0$ and use a circular aperture, rather than substituting the plane-illuminated Fresnel slit propagation profile. Useful additional controls are to retain the same laser and opening, centre the illumination, repeat position readings and check background with the beam blocked. If a camera replaces the visual screen readout, avoid saturation and calibrate its position scale. Alignment, imperfect aperture edges, illumination variation and limited contrast near a null constrain the inferred size. A bright central obstacle spot illustrates why the sharp shadows of geometrical optics can fail; it does not imply that every Fresnel diffraction pattern has a bright centre.

Off-axis holography tilts the holographic reference wave relative to the holographic object wave. Its interference fringes acquire a carrier spatial frequency; the two cross terms occupy different spatial frequency regions when the carrier is sufficiently large compared with the object's signal bandwidth. This permits filtering a single image term.

Gabor's lecture credits Emmett Leith and Juris Upatnieks with the later skew-reference method using helium–neon laser illumination, separating reconstructed images angularly.^[ref]

For an original periodic calculation, let $x$ run over a unit-length interval and chooseAll spatial frequencies here are cycles per unit length; the data are a chosen field at the recording plane, not measurements or a physical object profile. Direct expansion givesConsequently its Fourier series coefficients are $1.05$ at zero, $0.02$ at $\pm 8$, $0.2$ at $\pm 17$ and $0.1$ at $\pm 25$. The $R^{\ast }O$ component lies at $-17,-25$; the $R{O}^{\ast }$ component lies at $17,25$; the background lies at $0,\pm 8$. Filtering the negative-frequency cross term and multiplying by $R$ recovers $O$, including its relative wave phase.

The figure samples $N=1024$ equally spaced positions with the endpoint excluded, computes a normalized discrete Fourier transform using a Fast Fourier transform, retains bins $-30\le m\le -12$, then applies the inverse discrete Fourier transform and multiplies by $R$. The recovered real and imaginary parts coincide with the chosen $O$ to numerical precision. This exact periodic example avoids leakage by placing every frequency on an integer bin; finite nonperiodic records need allowance for broadening and an aperture.

For a continuous object field with Fourier transform confined to $|f|\le B$, the background $|O{|}^2$ occupies $|f|\le 2B$. A carrier $f_c>3B$ is a sufficient separation condition for the sideband centred at $-f_c$, and $f_c+B$ must remain below the Nyquist spatial frequency. Here $B=5$, $f_c=20$ and the Nyquist spatial frequency is $512$. This sufficient bound describes a whole bandwidth, rather than exploiting the gaps in this particular two-component field.

The same example also makes optical playback explicit. Choose $t_0=0.3$ and $\beta =0.01$. Since $|O|\le 0.3$, the bounds $0.49\le I\le 1.69$ imply $0.3049\le t\le 0.3169$, and $|\beta I|/t_0\le 0.0564$. The affine response is stipulated; its weak-modulation interpretation is approximate. In $Rt$, the desired field $0.01O$ has frequencies $3,-5$; the conjugate term $0.01R^2{O}^{\ast }$ has frequencies $37,45$; the background has frequencies $20,12,28$. Thus the same carrier that separates the recorded cross terms separates the playback components. Reillumination restores a field through diffraction; computational filtering extracts the corresponding encoded field through Fourier analysis.

An optical pinhole is a small aperture used to select a narrow part of a beam or approximate an optical point source. Its finite diameter and diffraction limit that approximation; an illuminated optical pinhole at a collimator's optical focus supplies a narrow bundle of nearly parallel light rays.

The detector is intended for arcsecond-scale astronomical seeing, rather than the telescope's diffraction limit. A circular aperture of this area has diameter $1.13\mathrm{m}$, giving a Rayleigh criterion angle about $0.11$ arcsec at the observing wavelength. One-arcsecond pixels conventionally sample an astronomical seeing profile of roughly two arcseconds across with two samples; a one-arcsecond profile would be undersampled. The pixel scale alone does not determine an exact astronomical seeing width.

Use the supplied surface brightness as radiative flux integrated over the photometric passband: its units contain no inverse hertz. Treat $h\nu$ as the representative photon energy, and assume unit throughput and quantum efficiency. A genuinely frequency-specific brightness would require a bandwidth. With angles measured consistently in arcseconds, an annulus has area $2\pi rdr$ square arcseconds, soThis is flux per collecting area; collected energy is $AtF$. There is no extra conversion to steradians when both brightness and angular area use square arcseconds.

For Poisson photon counting, a pixel with expected count $N$ has energy $e=Nh\nu$ and standard deviation ${\sigma }_e=h\nu \sqrt{N}$. Thus the local four-sigma rule gives $N=4\sqrt{N}$, or 16 photons at the boundary. With pixel area $\Delta \Omega =1{\mathrm{arcsec}}^2$,Consequently $r_m=(1-16/32.6991)/k=10.2138$ arcsec andIntegrating to the model's outer radius $R=1/k=20$ arcsec instead givesThe pixel cut retains 51.6% of this total. It is an imposed isophotal photometry rule, not a fundamental limit: averaging annuli or fitting a justified profile can recover signal below the per-pixel cut.

In practice, sky background photons add shot noise even after subtraction of the mean sky. read noise, dark current, imperfect flat-field correction, throughput calibration and atmospheric extinction also matter. At large radius the dominant uncertainty can be a tiny error in the sky level, scattered light or the wings of the point spread function, rather than source photon statistics. Mask foreground stars and background objects, compare blank-sky regions and dithered exposures, and check profile recovery by inserting artificial extended sources. astronomical seeing and finite pixels convolve the intrinsic profile; model that response when extrapolating the total.

The luminosity convergence of an outer surface-brightness profile follows directly: $I\propto r^{-p}$ contributes $dF\propto r^{1-p}dr$, hence $F(<R)\propto R^{2-p}$ at large radius for $p<2$. For either shallow index considered, ever larger annuli dominate: extending the law to infinity gives an unphysical divergent luminosity. Such a law must eventually steepen or truncate. Over a finite range, however, substantial light can lie below a surface-brightness selection threshold and be mistaken for sky. This biases galaxy luminosities, the luminosity function and the inferred diffuse light budget; distant diffuse galaxies become especially incomplete because bolometric cosmological surface-brightness dimming scales as $(1+z{)}^{-4}$. An assumed finite-aperture total is therefore a model-dependent quantity.

A thin lens has ray transfer matrixThe complex beam parameter transformation gives $q_2=q_1/(1-q_1/f)$, henceThe sign is also obtained by multiplying the incident wave amplitude by the thin lens phase $\exp [-ik{\rho }^2/(2f)]$.

Choose the thin lens plane as $z=0$ and put $a=\pi w_0^2/\lambda$, the incident Rayleigh range. Since the incident beam waist lies there, $q_1=-ia$. Immediately after the thin lens,Subsequent free propagation adds $z$, so write $q(z)=z-z_f-i{z}_R^{\prime }$. Comparing real and imaginary parts givesFor a positive focal length, the new beam waist lies downstream but closer than the geometrical focus. At fixed $f,\lambda$, as $w_0\to \infty$,This is the formal unbounded-aperture scalar limit; eventually finite aperture and the paraxial approximation limit the extrapolation.

Use time dependence $e^{-i\omega t}$, positive wavenumber $k$ and aperture normal $\mathbf{n}={\mathbf{e}}_x$. Write $F(y^{\prime },z^{\prime })=E_{\mathrm{inc},y}(0,y^{\prime },z^{\prime })$ on the aperture and extend $F$ by zero outside it. DefineThe outgoing branch has $\kappa \ge 0$ for propagating components and $\mathrm{Im}\kappa >0$ for evanescent components.

Fix the Green function normalization by $G_0(\mathbf{r})=e^{ik|\mathbf{r}|}/(4\pi |\mathbf{r}|)$, so $(\Delta +k^2)G_0=-\delta$. Its transverse Fourier transform satisfies $({\partial }_x^2+{\kappa }^2){\hat{G}}_0=-\delta (x)$. An outgoing solution on both sides is proportional to $e^{i\kappa |x|}$, and its derivative jump fixes ${\hat{G}}_0=i{e}^{i\kappa |x|}/(2\kappa )$. Thus the Weyl plane-wave representation for observation at $x>0$ isSince $\mathbf{n}\times {\mathbf{E}}_{\mathrm{inc}}=F{\mathbf{e}}_z$, the given aperture integral becomes a superposition of curls of ${\mathbf{e}}_z{e}^{i\mathbf{q}\cdot \mathbf{r}}$. Each curl supplies $i\mathbf{q}\times {\mathbf{e}}_z=i(\nu ,-\kappa ,0)$. Therefore, with the unit prefactor in that integral,Every vector amplitude is transverse to $\mathbf{q}$, as Maxwell's equations require. The polarizer specifies the incident aperture field; diffraction can regenerate an $x$ component outside the aperture. In contrast, its $z$ component remains zero in this representation. Evanescent plane waves are needed for the complete near field, not just the propagating disk ${\nu }^2+{\eta }^2\le k^2$.

The factor $1/2$ follows explicitly from the standard outgoing Green function for the three-dimensional Helmholtz equation and the printed unit curl prefactor. If the intended kernel is the doubled half-space kernel $2G_0$, or the aperture formula carries the usual compensating factor $2$, delete that $1/2$. The latter normalization makes the tangential output trace equal to $F$. Keeping the convention explicit avoids silently changing the supplied representation; the angular dependence and transversality of an electromagnetic plane wave are the same in either normalization. Reversing the specified normal also reverses the overall sign.

Let the two tangential aperture fields be $F_y(y,z)$ and $F_z(y,z)$, extended by zero off the aperture. The normal component does not enter $\mathbf{n}\times {\mathbf{E}}_{\mathrm{inc}}$, which equals $-F_z{\mathbf{e}}_y+F_y{\mathbf{e}}_z$. Apply the Weyl plane-wave representation and curl to both components. With the same kernel normalization as part (a),This is the electromagnetic aperture angular spectrum. The two electromagnetic polarization vectors are linearly independent and both orthogonal, bilinearly, to $\mathbf{q}$; they need not be orthogonal to each other. Arbitrary tangential electromagnetic polarization is therefore represented, with its outgoing longitudinal component fixed by $\kappa E_x+\nu E_y+\eta E_z=0$ rather than independently prescribed. Away from grazing components one could equally express the spectrum in a transverse-electric and transverse-magnetic basis. The same limiting outgoing branch handles grazing and evanescent components.

For each component the magnetic amplitude follows from Faraday's law as $\mathbf{B}=\mathbf{q}\times \mathbf{E}/\omega$. Together with ${\mathbf{q}}^2=k^2$ and transversality of an electromagnetic plane wave, this verifies the source-free Maxwell equations. If the doubled-kernel aperture convention is used, multiply the entire expression by two, exactly as in part (a).

For $p^2+q^2>1$, $m=i\sqrt{p^2+q^2-1}$ and the Weyl plane-wave representation contains evanescent waves. Their amplitudes decay exponentially with the distance from the scattering region. Excluding them removes high transverse spatial frequencies, including subwavelength information. Near-field measurements can in principle detect them, but back-propagating their decay exponentially amplifies measurement errors.

For propagating data, the Fourier transfer satisfies $|{\mathbf{k}}_s-{\mathbf{k}}_i|\le 2k$. One illumination still samples only the Ewald sphere, so this frequency bound does not imply access to every mode inside that ball. The lack of a volumetric spectrum is the separate nonuniqueness described in the previous part.

The Born approximation requires the incident field to remain a good approximation inside the scatterer. A useful sufficient smallness estimate for a weak, extended inclusion is small accumulated phase $kD|n-1|\ll 1$, alongside weak reflection and multiple scattering. Small pointwise contrast alone can fail for a thick inclusion or a resonance. The scalar Helmholtz and background-medium assumptions, and the extra linearization of $n^2-1$, also belong to the model's validity conditions.

Numerically, a finite aperture windows the measured fields and limits angular and Fourier resolution. Discrete sampling must satisfy the Nyquist–Shannon sampling theorem for the retained transverse bandwidth; spacings $\Delta x,\Delta y\lesssim \pi /K_{\max }$ avoid its simplest aliasing. Finite discrete Fourier transforms also introduce periodic continuation, truncation, interpolation and window errors. Noise, unresolved evanescent modes and near-grazing modes with $m\simeq 0$ limit reliable recovery. Refining the grid cannot by itself fill the missing three-dimensional Fourier data.

Fix a receiver position $\mathbf{\zeta }$ and define $g_{\mathbf{\zeta }}(\mathbf{z})=G_L(\mathbf{\zeta },\mathbf{z})$. Choose the two scalar fields in the Wigner distribution to be $u=v=g_{\mathbf{\zeta }}$. Fourier inversion in the separation variable givesIndeed, substituting the defining Wigner distribution makes the $\mathbf{p}$ integral a $(2\pi {)}^2$ Dirac delta function, canceling its normalization and setting the separation equal to $\mathbf{z}-{\mathbf{z}}^{\prime }$.

Average over the receiver plane, including the aperture:Then the receiver-averaged Wigner distribution represents the back-propagated field asThe fast decay of $f$ justifies its use as a test function. For oscillatory Green functions that are not integrable, the Wigner distribution and its inversion can be understood as tempered distributions, or derived with smooth cutoffs before taking their limits.

The printed $W_A$ expression can also be used literally, but its displayed ${\mathbf{z}}^{\prime }$ integral is over source-plane coordinates, not receiver positions. The missing receiver average and source weight can be encoded by choosing vector fields whose channels are the receivers. LetContract the receiver channels in the product of the vector fields, so thatWith exactly the printed definitionthe same Fourier inversion now yields the requested compact formFor a finite array, the receiver integral in the contraction is a weighted sum. For a continuous array, these are fields valued in the receiver Hilbert space. An uncontracted vector outer product instead gives a matrix-valued Wigner distribution, whose receiver trace must be taken. Thus the printed formula is usable with these choices and this contraction convention, but its phrase “over the plane of the receiver” does not describe the displayed integral.

In time reversal acoustics, the array records the incoming signal, reverses each recorded time trace, and re-emits it through the same medium. This is phase conjugation in the frequency domain. For the time-harmonic wave convention $e^{-i\omega t}$, reversing a real time trace replaces its positive-frequency wave amplitude by its complex conjugate. The medium must remain unchanged between recording and re-emission.

Let $f(\mathbf{z})={\Psi }_0(0,\mathbf{z})$ and $G_L(\mathbf{\zeta },\mathbf{z})=G(L,\mathbf{\zeta };0,\mathbf{z})$, with the chosen source and array normalizations incorporated into the Green function. Let $a(\mathbf{\zeta })\ge 0$ be the array's aperture weight, equal to the indicator of its receiving region for an ideal uniform array. The recorded field isBy wave reciprocity, back-propagation has the same Green function with the source and receiver exchanged. Therefore the physically re-emitted, back-propagated wave amplitude isIf $Hf=R$ and $P_A$ multiplies by $a$, then ${\Psi }_B=H^T{P}_A\overline{Hf}$, where $T$ denotes the transpose without conjugation. Taking a final complex conjugate instead defines the adjoint reconstruction $\overline{{\Psi }_B}=H^{\ast }{P}_{A}Hf$. This distinction prevents an erroneous conjugation in the time reversal operator.

For a localized Gaussian beam or acoustic point source in a homogeneous medium, a finite aperture admits a limited range of angles. The focal width is of order $\lambda L/A$ when $A$ denotes the aperture diameter. In a random medium, multiple scattering creates paths with a larger angular spread. Each reversed path retraces its route, and the paths interfere constructively at the source. This can produce a larger effective aperture in time reversal and a narrower focus, even though the unreversed field has a complicated speckle pattern.

This comparison concerns a homogeneous reference medium; a deterministic heterogeneous medium can also provide useful multipath propagation. Suitable scale limits or frequency and spatial averaging can make refocusing self-averaging. Such self-averaging is not automatic for every monochromatic source and every random realization. Wave absorption, changing medium parameters, unresolved paths or poor array coverage can spoil refocusing. With complete capture of the propagating modes and a lossless unitary operator $H$, ideal adjoint reconstruction is already exact in either medium. Random scattering can improve finite-aperture wave focusing through angular diversity. See the regime-dependent analysis in Statistical stability in time reversal.

The polarization generated by the relevant third-order nonlinear optical susceptibility is proportional to $E_1{E}_2{E}_3^{\ast }$. Multiplication of its temporal and spatial wave phase factors givesThe two opposing pump arrows in the original diagram mean ${\mathbf{k}}_2=-{\mathbf{k}}_1$. Thus the mixing source points along $-{\mathbf{k}}_3$, and its conjugate signal amplitude drives a backward field. This is counterpropagating-pump phase matching in four-wave mixing.

There is a phase matching qualification to the initial frequency statement. In the stipulated nondispersive isotropic medium, $|{\mathbf{k}}_j|=n{\omega }_j/c$. Exact equality of the backward free-wave vector to $-{\mathbf{k}}_3$ therefore requires ${\omega }_4={\omega }_3$, hence ${\omega }_3={\omega }_P$. For unequal frequencies the backward direction remains the geometry suggested by the source, but the longitudinal phase mismatch has magnitude $2n|{\omega }_P-{\omega }_3|/c$. Efficient generation over a finite length requires that mismatch times $L$ be small; a sinc function factor suppresses conversion otherwise. Nondispersion does not remove the magnitude mismatch. The remainder of the calculation is explicitly degenerate four-wave mixing, so this difficulty disappears.

At the pump frequency, for the strong field $E_1$, self-phase modulation and cross-phase modulation by $E_2$ dominate the weak signal's influence. They change its propagation wave phase through its own intensity and the other pump's intensity. With continuous uniform plane waves, these are nonlinear wave phase shifts; finite pulses acquire temporal chirp, and finite transverse beams can acquire a lens. No universal factor of two should be assigned between the two shifts without specifying their polarization contraction and the field-amplitude convention.

Write $u={\tilde{A}}_3$ and $v={\tilde{A}}_4$. Differentiating the exponential substitutions cancels the terms proportional to $2\kappa$:Consequently $u^{\prime }=i\kappa v^{\ast }$ and $v^{\prime }=-i\kappa u^{\ast }$. Conjugate the latter before eliminating a field: $(v^{\ast }{)}^{\prime }=i\kappa u$, so $u^{\prime \prime }=-{\kappa }^{2}u$. The backward four-wave-mixing boundary problem requires opposite-end boundary conditions: $u(0)=a$, where $a=A_3(0)$, and $v(L)=0$. The transfer matrix isAt $L$, $0=ia\sin \kappa L+b\cos \kappa L$, so $b=-ia\tan \kappa L$. Substitution givesThe printed target expression for $A_4$ has the opposite overall sign. This is present in the original PDF, not just the transcription. With its stated first-order equations, the derivative at the far face must be $v^{\prime }(L)=-i\kappa a^{\ast }/\cos \kappa L$. The boxed result has that derivative; the printed expression has its negative. Reversing the convention for the mixing coefficient could reverse the generated phase, but would also change the coupled equations. The result above satisfies the equations and both boundary conditions together.

At the entrance, $A_4(0)=i{a}^{\ast }\tan \kappa L$, so the returned amplitude is proportional to the complex conjugate of the signal. An ordinary ideal mirror gives $ra$, with a fixed reflection coefficient $r$, rather than $a^{\ast }$. The mixing pumps supply both the extra energy and a wave phase reference. The ideal conjugate reflectivity and signal transmission areThe constancy of $|u{|}^2+|v{|}^2$ along the medium is another direct check. The relation between the incoming signal and outgoing waves allows net output energy to come from the pumps. As $\kappa L\to \pi /2$, the undepleted-pump approximation fails: the predicted fields would extract an appreciable pump fraction. This is the ideal four-wave-mixing oscillation threshold; saturation, pump depletion and noise seeding replace the divergent linearized solution. It is not an unlimited source of reflected power. Amplified phase-conjugated propagation and parametric oscillation are also the effects identified in the primary analysis by Amnon Yariv and David M. Pepper.^[ref]

For the aberrating linear medium, assume its permittivity is real and stationary at this frequency. Taking the complex conjugate of the forward-envelope equation givesThis is the backward-envelope equation with carrier $e^{-i{k}_{3}z}$. Launching the conjugate boundary field therefore permits retracing of the distorted wave. In a transverse lossless phase screen model, let transmission be $t(x)=e^{i\psi (x)}$. A return proportional to $(tU{)}^{\ast }$ acquires $t$ on its second passage and becomes $t(tU{)}^{\ast }=U^{\ast }$. Thus the aberration phase cancels. A mirror instead gives $r{t}^{2}U$, adding the aberration twice; merely reversing the central propagation direction does not reverse the wavefront phase.

The one-dimensional printed equation describes longitudinal variation; genuine transverse aberrations require the corresponding transverse Helmholtz equation or paraxial operator, for which the same real-coefficient conjugation argument applies. Exact reconstruction requires adequate aperture, the correct conjugate amplitude and polarization, uniform pump-product phase, and an unchanged reciprocal lossless medium. Absorption, missing scattered modes or a changed aberrator prevent a universal reconstruction. This is phase conjugation and thin phase-plate cancellation by phase conjugation under specified conditions, rather than reversal of every possible propagation loss.