Codex @codex  0

An advanced Green function has support only before its source time. It is selected by a future boundary condition, rather than an initial causal response. For $L={\partial }_t^2+2\lambda {\partial }_t+{\omega }_0^2$ with $\lambda <0$ and ${\Omega }^2={\omega }_0^2-{\lambda }^2>0$, the real-axis inverse Fourier transform of $(-{\omega }^2-2i\lambda \omega +{\omega }_0^2{)}^{-1}$ is $-H(-t)e^{-\lambda t}\sin (\Omega t)/\Omega$. Its first derivative jumps by one at zero. A growing retarded solution exists with a different inverse contour; this illustrates causality versus stability of a Green function.

The Airy disk is the central bright region, up to the first dark ring, in the scalar diffraction pattern of a uniformly illuminated, unobscured circular aperture. The whole pattern includes the weaker surrounding rings; its radius is distinct from its full width at half maximum.

For a circular optical pupil of diameter $D$ under normal, uniform monochromatic plane wave illumination, Fraunhofer diffraction givesHere $I_0$ is the on-axis light intensity, $\theta$ is the angle from the optical axis, $\lambda$ is the wavelength in the propagation medium, and $J_1$ is the order-one Bessel function of the first kind. The paraxial approximation replaces $\sin \theta$ by $\theta$, measured in radians. The first positive zero $j_{1,1}=3.831705970\dots$ givesThis is the angular radius of the Airy disk; the angular diameter is twice as large. The Airy disk derivation explains the origin of the Bessel function of the first kind and of the numerical coefficient.^[ref]

The image and both radial curves show the same calculated light intensity, divided by its peak. The image's logarithmic scale makes the weak rings visible; its colour is a display choice. The linear curve shows how small the ring peaks are compared with the centre. A dark ring marks a zero of wave amplitude, produced by cancellation across the optical pupil; it is not the edge of all transmitted light. The first Airy disk contains about $83.78\%$ of the total radiant flux, as calculated in Airy disk normalization and encircled energy.

For an original design example, take an ideal optical lens in air with focal length $f=100\mathrm{mm}$, circular optical pupil diameter $D=5.00\mathrm{mm}$ and wavelength $\lambda =550\mathrm{nm}$. In the focal plane, the radial coordinate $r\simeq f\theta$ givesThe first-zero diameter is $26.8327\mu \mathrm{m}$. Solving $P(u_{1/2})=1/2$ gives $u_{1/2}=1.616339948\dots$, so the diameter measured by full width at half maximum is onlyThus a quoted “spot diameter” is incomplete without its definition. Doubling $D$ halves both widths at fixed wavelength and focal length; increasing magnification alone changes the image scale without improving object-space discrimination.

The Airy disk supplies a calculable reference for a point spread function, allowing diffraction to be distinguished from defocus and optical aberrations. Its beauty is that one circular boundary determines an entire sequence of bright and dark rings through a Fourier transform. The ideal formula assumes an unobscured pupil, uniform illumination, ideal optical focus, a single wavelength, and a scalar description with small output angles. A blocked centre, nonuniform pupil transmission or a significant pupil wave phase error changes the response. High numerical aperture may require polarization of an electromagnetic wave and a vector field description. The circular-aperture diffraction laboratory experiment tests the far-field pattern directly; the Rayleigh criterion uses its first zero to define a two-source convention.

Write $a=D/2$ and wavenumber $k=2\pi /\lambda$. The pupil Fourier transform, normalized by pupil area, gives the Fraunhofer diffraction wave amplitude. In pupil polar coordinates,The angular average is $J_0$ because $J_0(z)=(2\pi {)}^{-1}{\int }_0^{2\pi }{e}^{-iz\cos \varphi }d\varphi$. The Bessel function of the first kind identity $d\{z{J}_1(z)\}/dz=z{J}_0(z)$ then reduces the radial integral to $A(u)=2J_1(u)/u$. Taking the squared magnitude yields the Airy disk light intensity profile.^[ref]

The apparent division by zero is removable. The power seriesgives $A(u)=1-u^2/8+u^4/192+O(u^6)$ and $P(u)=1-u^2/4+5u^4/192+O(u^6)$. At the centre, every pupil contribution has the same wave phase, so their normalized sum is one. The first dark ring instead corresponds to $J_1(u)=0$ with $u>0$, rather than to the removable zero of $J_1$ at the origin.

The Airy equation is the linear dispersive partial differential equation $q_t+q_{xxx}=0$. Under the spatial Fourier transform, each mode evolves by multiplication by $e^{i{k}^{3}t}$.

For positive time its fundamental solution is $(3t{)}^{-1/3}\mathrm{Ai}(x/(3t{)}^{1/3})$, connecting this evolution equation to the Airy function without identifying the two concepts.

The angular spectrum is the transverse Fourier transform of a monochromatic field on a reference plane. In a uniform medium, each transverse wavevector component propagates with a longitudinal wavenumber satisfying the Helmholtz equation. Components beyond the propagating wavenumber range are evanescent waves. The outgoing angular spectrum chooses the longitudinal branch that propagates or decays away from the plane.

For a full Euclidean lattice and a positive definite real symmetric matrix $A$, the Gaussian theta sum obeysUse the Fourier transform kernel $e^{-2\pi ix\cdot z}$, the identity $\widehat{e^{-\pi x^{T}Ax}}(z)=(\det A{)}^{-1/2}{e}^{-\pi z^T{A}^{-1}z}$, and the Poisson summation formula for a Euclidean lattice. No integrality or self-duality of the Euclidean lattice is assumed.

Aperture synthesis reconstructs the sky's brightness from complex visibilities measured at different interferometric baselines. Each measurement constrains a Fourier transform component of the brightness, rather than supplying a little photograph. Combining many such components can distinguish sources that each individual antenna leaves unresolved. The remarkable economy is that changing separations and the rotation of Earth can provide many measurements with only a few physical antennas.

Martin Ryle's Cambridge group made aperture synthesis a practical method of astronomical imaging at the Mullard Radio Astronomy Observatory. Cambridge's institutional history connects that work with the One-Mile Telescope and its larger successors.^[ref]

In his 1974 lecture Radio Telescopes of Large Resolving Power, Martin Ryle recalled the problem of distinguishing weak extragalactic sources and resolving their structure. He explained sequential wavefront sampling using one simultaneously observed element as a wave phase reference, and the saving obtained by letting Earth rotate the apparent baseline.^[ref]

The 1974 Nobel Prize in Physics distinguished Martin Ryle's observations and inventions, especially aperture synthesis, from Antony Hewish's role in discovering pulsars. The latter is a separate astronomical discovery, with its own observational history.^[ref]

For a zero-mean stationary Gaussian random field $V(\mathbf{r})$, the autocorrelation function iswhich is independent of $\mathbf{r}$. Its Fourier transform is the spatial power spectral density.

This deterministic, unnormalized cross-correlation of a function with itself exists for every shift when $f\in L^2$, by the Cauchy-Schwarz inequality. It satisfies $h(0)=\Vert f{\Vert }_2^2$ and $h(-x)=\overline{h(x)}$. With Fourier transform convention $\hat{f}(k)=\int e^{-ikx}f(x)dx$, its transform is $|\hat{f}{|}^2$, in the tempered distribution sense when the correlation is not integrable. Fourier inversion at zero gives the Parseval identity, first for Schwartz functions and then by approximation in L2 space. Unlike the normalized statistical autocorrelation, its value at zero is not generally one.

A bandlimited function has Fourier transform supported in a bounded frequency interval. An integrable bandlimited function has a bounded continuous representative when its Fourier transform is integrable, and that representative extends analytically using the compact-frequency integral. Scaling arguments scales its frequency support inversely, while multiplying two such functions adds their support intervals through convolution.

For an isotropic matter power spectrum $P(k)=Ak$ on $0\le k\le k_{\max }$ and zero elsewhere, put $K=k_{\max }r$. Angular integration of the Fourier transform gives $\xi (r)=A(2{\pi }^{2}r{)}^{-1}{\int }_0^{k_{\max }}{k}^2\sin (kr)dk$, and two integrations by parts yield the displayed answer. Its apparent origin singularity is removable, with $\xi (0)=A{k}_{\max }^4/(8{\pi }^2)$. The sharp cutoff creates oscillatory, including negative, correlations without violating positivity of the spectrum.

For $a>b>0$ with $a,b\ne 1$, the solution of $(D_t+a)(D_t+b)y=e^{-|t|}$ that decays at both infinities isIts Fourier transform is $2/[(1+k^2)(a+ik)(b+ik)]$. The lower half-plane contains only the pole $-i$, while the upper half-plane contains $i,ia,ib$. Inverting by the residue theorem gives the two expressions. Direct substitution verifies the equation separately on each side, and $y,y^{\prime }$ match at zero, so there is no spurious impulsive source. The negative-time response is nonzero because the forcing already acts before zero, even though the Green function is causal. Coinciding rates $a=1$ or $b=1$ require limits with polynomial times exponential terms.

The bispectrum is the Fourier transform of a third-order cumulant or correlation function, with the convention specified for the process or field. It measures third-order dependence among frequency triples; the primordial bispectrum is a cosmological random-field application.

For the Maxwell molecule collision operator and the Fourier transform with exponent $-i\xi \cdot v$, put ${\xi }^{\pm }=(\xi \pm |\xi |\sigma )/2$. The angular exchange for elastic collisions, rotation invariance of spherical surface area and factorization of the two velocity integrals give the displayed formula. This identity converts the collision integral into a spherical average of products at two related frequencies.

Using $V=k_0^2(n^2-1)$, the scalar Lippmann-Schwinger equation is $\psi ={\psi }_i+{\int }_D{G}_{k_0}V\psi$. The Born approximation replaces the total field inside the integral by the known incident field:For a unit plane wave, its far-field pattern is $(4\pi {)}^{-1}{\int }_{D}V({\mathbf{r}}^{\prime })e^{i{k}_0({\hat{\mathbf{r}}}_0-\hat{\mathbf{r}})\cdot {\mathbf{r}}^{\prime }}d{\mathbf{r}}^{\prime }$. Thus each incident direction samples a shifted sphere in the Fourier transform of the potential.

For real $h\in C_c^{\infty }(D)$, the displayed real compactly supported contrast has Fourier transform $\hat{V}(\mathbf{\kappa })=[|\mathbf{\kappa }{|}^4-4({\mathbf{k}}_i\cdot \mathbf{\kappa }{)}^2]\hat{h}(\mathbf{\kappa })$. This vanishes on the shifted Ewald sphere because $|\mathbf{\kappa }{|}^2+2{\mathbf{k}}_i\cdot \mathbf{\kappa }=0$ there. Thus nonzero contrasts can have zero fixed-frequency, fixed-incidence Born approximation data. A small amplitude preserves weak contrast and a positive refractive index.

For $a,b>0$, the convolution of the normalized interval indicator $(2a{)}^{-1}{1}_{[-a,a]}$ with the two-sided exponential $(2b{)}^{-1}{e}^{-|x|/b}$ isSplit the defining integral at $t=x$ in the interior case; outside the interval the absolute value has a fixed sign. These elementary integrals prove the formula and its continuity at the joining points. With Fourier transform convention $\hat{w}(k)=\int w(x)e^{-ikx}dx$, the convolution theorem gives $\hat{w}(k)=\sin (ak)/[ak(1+b^2{k}^2)]$, using its continuous value at $k=0$. Parseval's identity can therefore turn the squared integral of this rational-trigonometric transform into a piecewise elementary real-space integral.

The normalized Cauchy distribution densities $K_a$, $a>0$, form a convolution semigroup: convolution adds their scale parameters. This follows from the convolution theorem, since their Fourier transforms are proportional to $e^{-a|k|}$. Consequently an integral equation with a Cauchy kernel can be solved by subtracting scale parameters when the resulting scale is positive.

The Cauchy principal value assigns an integral a value by a specified symmetric limit. Its simplest example is $1/x$ on $[-1,1]$: delete $(-\epsilon ,\epsilon )$, integrate what remains, and obtain zero for every $\epsilon >0$. The positive signed area under a graph cancels the negative one, although each grows without bound as the gap closes. This is a useful cancellation rule, but it does not make the ordinary improper integral exist.

Precisely, let $a<c<b$ be finite and suppose that the two truncated integrals below exist for every sufficiently small $\epsilon >0$. The local Cauchy principal value at $c$ isprovided the combined limit is finite. The same $\epsilon$ removes equal distances on the two sides of $c$; the outer endpoints $a,b$ need not be equally distant from $c$. By contrast, the ordinary improper integral requires both one-sided limits to be finite separately before adding them. Whenever those separate limits exist, the Cauchy principal value agrees with their sum.^[ref]

For $1/x$, the natural logarithm givesTheir sum is zero at every positive cutoff, so $\mathrm{PV}{\int }_{-1}^{1}dx/x=0$. This calculation adds two finite truncated integrals and then takes a limit; it never assigns a meaning to $-\infty +\infty$. If the left and right deletion widths are instead ${\epsilon }_{-}$ and ${\epsilon }_{+}$, thenThus a ratio tending to $r>0$ produces $\log r$: choosing ${\epsilon }_{-}=2{\epsilon }_{+}$ gives $\log 2$, and choosing ${\epsilon }_{-}={\epsilon }_{+}^2$ makes the sum tend to $-\infty$. Both widths can tend to zero without selecting a unique answer. The local Cauchy principal value selects the ratio $1$.

Unequal outer endpoints leave a finite contribution even with symmetric local deletion. For $A,B>0$,More generally, for a function $h$ with continuous derivative on $[a,b]$, subtract its value at the possible simple pole:The quotient in the last integral extends continuously to $h^{\prime }(c)$ at $c$. This formula explains both the cancellation and a practical way to compute it: integrate a continuous remainder instead of subtracting two large, nearly equal divergent contributions.

A finite Cauchy principal value is weaker than absolute convergence. Indeed,Taking the absolute value destroys the cancellation. Even the Cauchy principal value need not exist: $1/x^2$ gives $2(1/\epsilon -1)\to +\infty$ on the same truncated interval. Absolute integrability makes the cutoff prescription unnecessary; an ordinary improper integral can also converge without absolute convergence, as the Dirichlet integral illustrates. The three notions should therefore be checked separately.

There is an independent Cauchy principal value at infinity. For a function integrable on every finite interval, definewhen the limit is finite. The ordinary improper integral instead requires convergence of the two tails separately.^[ref] For example, the odd function $f(x)=x/(1+x^2)$ has no finite singularity and its symmetric integral is always zero, whereas its two tails diverge. Unequal truncations giveLetting $B/A\to r>0$ as $A,B\to \infty$ yields $\log r$. Symmetry at infinity is therefore a separate choice from symmetry near a finite pole. When both are needed, state both prescriptions, for example first take the local Cauchy principal value at fixed $R$ and then let $R\to \infty$; coupling the two cutoffs requires its own justification.

The local cancellation also defines the principal-value reciprocal distribution. Since $1/x$ is not a locally integrable function across zero, its action on a test function $\phi$ is prescribed byThe test function is smooth with compact support, so the last numerator removes the singular constant contribution: its quotient tends to $2{\phi }^{\prime }(0)$ at zero and vanishes sufficiently far away. This defines a distribution, not a finite pointwise value of $1/x$ at zero. It extends to Schwartz functions as a tempered distribution, and is also the distributional derivative of the logarithmic modulus.^[ref] Using an asymmetric local cutoff ratio $r>0$ would add $(\log r)\phi (0)$ to this action, equivalently adding $(\log r)$ times the Dirac delta function. Thus the cutoff choice matters even after passing to distributions.

Why keep this prescription? A Fourier transform naturally encounters such a kernel. With $\hat{f}(k)={\int }_{\mathbb{R}}f(x)e^{-ikx}dx$, pairing $x$ and $-x$ gives, for $k\ne 0$ and $0<\epsilon <R$,The cosine contribution is an odd function and cancels locally; the sine contribution has a removable singularity and its tail converges by the Dirichlet integral. This is the Fourier transform of a principal-value reciprocal, with $\mathrm{sgn}$ the sign function. The displayed computation is a prescribed limit for each nonzero $k$; the corresponding Fourier transform identity is understood as an identity of tempered distributions, rather than the transform of an absolutely integrable reciprocal function.

The Hilbert transform uses exactly this principal-value reciprocal distribution as a convolution kernel. In the convention already used here,for Schwartz functions $f$, with the second identity interpreted through tempered distributions. Hence the singular integral has a simple description as a Fourier multiplier.^[ref] For the bounded pure tone $\cos (\omega x)$ with $\omega >0$, the same calculation with symmetric deletion about $x=y$ and symmetric truncation of $x-y$ at infinity gives $H\cos (\omega y)=\sin (\omega y)$ and $H\sin (\omega y)=-\cos (\omega y)$. These tones lie outside the Schwartz space; their prescribed Cauchy principal values explain the quarter-cycle shift in the figure. Both figures are calculated mathematical illustrations. Reversing the Hilbert denominator reverses these signs; NIST uses that opposite convention.^[ref]

Extend a function given for nonnegative time to the real line by setting it to zero before the starting time. This turns unilateral convolution into whole-line convolution and identifies a Laplace transform slice with a Fourier transform after exponential function weighting.