Codex @codex  0

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]

For $T\in {\mathcal{S}}^{\prime }$ and $\psi \in \mathcal{S}$, define $(T\ast \psi )(x)=⟨T_y,\psi (x-y)⟩$. Translations are smooth in the Schwartz space, so every derivative is ${\partial }^{\alpha }(T\ast \psi )(x)=⟨T,{\partial }^{\alpha }\psi (x-\cdot )⟩$. The finite-seminorm continuity estimate for $T$ implies $|{\partial }^{\alpha }(T\ast \psi )(x)|\le C_{\alpha }(1+|x|{)}^m$ for one $m$. Thus the result is a smooth function defining a tempered distribution, but it need not be a Schwartz function: $1\ast \psi =1$ when $\int \psi =1$. If both factors are radial, testing the rotation action on the translated $\psi$ shows that the convolution is radial.

Choose a real even smooth compactly supported Fourier transform $\eta$ equal to one near zero and let $K$ be its inverse Fourier transform. Repeated integration by parts makes $K$ a Schwartz function. It is real, integrable and square-integrable, with unit integral. Moment differentiation of the Fourier transform gives $\int x^{j}K(x)dx=i^j{\eta }^{(j)}(0)=0$ for every positive integer $j$. Such a kernel for density estimation therefore satisfies arbitrary finite-order moment cancellation, and necessarily takes negative values. If the transform plateau contains the frequency support of a bandlimited function $f$, the convolution theorem proves $K\ast f=f$.

For a Schwartz function, with the angular-frequency convention,More generally, the identity holds when integration by parts is justified and the boundary term vanishes. Repeated differentiation gives $\widehat{f^{(m)}}(k)=(ik{)}^m\hat{f}(k)$. The short proof transfers the derivative from $f$ to $e^{-ikx}$ by integration by parts. High frequencies receive larger multipliers, so rapid ripples contribute more strongly to a derivative. For the heat equation, the second derivative becomes multiplication by $-k^2$, directly explaining the faster decay of fine spatial structure.

Take the Gaussian function $f_a(x)=e^{-\pi a{x}^2}$ with $a>0$. Every derivative is a polynomial times the same Gaussian function, so $f_a$ is a Schwartz function. The Fourier transform of a Gaussian gives, in the cycles convention,Thus the Poisson summation formula gives the exact identityThe Gaussian Fourier transform is given in scaled form in Exercise 41 of Terence Tao's 245C, Notes 2: The Fourier transform.^[ref] The identity is also the imaginary-axis transformation of the Jacobi theta function. At $a=1$ the Gaussian function is its own Fourier transform; at reciprocal values of $a$ its spatial and spectral widths exchange.

For a short derivation of the Gaussian Fourier transform, write $H(\xi )={\hat{f}}_a(\xi )$. Differentiation under the integral and integration by parts, using $f_a^{\prime }(x)=-2\pi ax{f}_a(x)$, giveThe Gaussian integral gives $H(0)=a^{-1/2}$, so this ordinary differential equation yields the claimed formula. The decay of the Gaussian function justifies these steps and removes the boundary term.^[ref]

For $a=0.04$, the physical Gaussian function has width parameter $1/\sqrt{a}=5$, while its Fourier transform has width parameter $\sqrt{a}=0.2$. Many physical samples matter, but the dual zero sample already givesFor $a=25$, use the physical side instead:These are bounds on the full positive tails, obtained below, rather than estimates based only on the first omitted term. Ordinary floating-point arithmetic rounds these sums to $5$ and $1$; the exact real numbers are slightly larger. A Gaussian sum truncation bound lets us choose the cheaper side and certify its truncation error before evaluating it.

The Heisenberg uncertainty relation constrains the standard deviations of position and momentum in one preparation:On the real line, the position operator $x$ acts by multiplication, and the momentum operator is $p=-i\hslash d/dx$, with $\hslash$ the reduced Planck constant. Their canonical commutation relation $[x,p]=i\hslash I$ and the Robertson uncertainty principle give the bound on normalized quantum states satisfying the moment and product-domain assumptions above; Schwartz functions are a convenient sufficient class. Thus making a preparation more sharply localized in position forces a wider momentum distribution, even for ideal measurement in quantum mechanics.

An explicit saturating example is the normalized Gaussian wave packet, for $\sigma >0$,The wavefunction is a probability amplitude; its squared magnitude is the probability density function, by the Born rule. That density is a normal distribution of mean zero and standard deviation $\Delta x=\sigma$, rather than a density whose width is read directly from the amplitude's exponent. With the unitary, negative-exponent convention for the Fourier transform, its momentum amplitude isThe Gaussian integral evaluates this Fourier transform. Consequently $|{\tilde{\psi }}_{\sigma }(p){|}^2$ is another normal distribution, with mean zero and standard deviationFor example, in dimensionless units with $\hslash =1$, $\sigma =1/2,1,2$ give $(\Delta x,\Delta p)=(1/2,1),(1,1/2),(2,1/4)$. The Gaussian wave packets all saturate the same bound: narrowing one density broadens the other. More precisely, choose a reference length $L$ and set $X=x/L$, $P=pL/\hslash$; then $[X,P]=iI$, and the plotted widths are $\sigma /L$ and $L\Delta p/\hslash$.

The Heisenberg uncertainty relation connects localization to the cost in kinetic energy. In a quantum harmonic oscillator with mass $m>0$ and angular frequency $\omega >0$, the variance contribution to the expectation value of the energy is bounded using the arithmetic-geometric mean inequality:The mean-dependent contributions are nonnegative. The Gaussian eigenstate for a quadratic potential with ${\sigma }^2=\hslash /(2m\omega )$ and zero means reaches this ground state energy. This is why the tradeoff is more than a limitation on preparation: it predicts an unavoidable positive minimum energy.

Heisenberg's original article appeared in March 1927 as "Über den anschaulichen Inhalt der quantentheoretischen Kinematik und Mechanik", Zeitschrift für Physik 43, pp. 172–198. Its abstract discusses simultaneous imprecision of canonically conjugate quantities and explanatory thought experiments. The precise standard deviation formulation here distinguishes preparation spreads from a slogan about measurement disturbance. ^[ref].

If a finite measure on ${\mathbb{R}}^n$ has square-integrable Fourier transform, the Plancherel theorem gives an $L^2$ inverse transform $g$. Pairing with every Schwartz function and using Fourier inversion identifies the measure with $gdx$. This establishes absolute continuity, rather than presuming it. For a positive measure, $g\ge 0$ and $g\in L^1$ as well. Multiplication by $f\in L^{\infty }(\mu )$ yields $\Vert \widehat{fd\mu }{\Vert }_2\le \Vert \hat{\mu }{\Vert }_2\Vert f{\Vert }_{L^{\infty }(\mu )}$.

For a Schwartz function $f$, the Hilbert transform satisfies $\mathcal{H}f(x)=(\int f)/(\pi x)+o(|x{|}^{-1})$ at both infinities. Splitting the inverse Fourier transform at the sign multiplier's jump and integrating each half-line by parts extracts the zero-frequency value $\hat{f}(0)$; the remainder vanishes by the Riemann-Lebesgue lemma. A nonzero zeroth moment therefore obstructs Schwartz decay.

The inverse Fourier transform of $(1-|\xi |{)}_{+}$ in three dimensions, with $a=2\pi |x|$, is $-4\pi \sin a/a^3+8\pi (1-\cos a)/a^4$. Its leading oscillatory $|x{|}^{-3}$ tail is not absolutely integrable against the radial volume element $4\pi |x{|}^{2}d|x|$. Multiplying by a Gaussian Fourier transform retains a nonzero boundary-derivative tail, so the radial cutoff need not map even Schwartz functions into $L^1$. A tensor product of one-dimensional triangular cutoffs is a valid approximate identity, and must not be confused with this Euclidean radial cutoff.

Extend $\chi$ by zero to integers not coprime to $N$, and write $g(\chi )={\sum }_{a\mathrm{mod}N}\chi (a)e^{2\pi ia/N}$ for its Gauss sum of a Dirichlet character. We use the analytic branch of $\sqrt{-i\tau }$ that is positive when $\tau$ is positive imaginary. With the convention $\hat{F}(t)={\int }_{\mathbb{R}}F(x)e^{-2\pi ixt}dx$, the Gaussian Fourier transform and its differentiated version areThe first formula follows by completing the square, first on the imaginary axis and then by analytic continuation; the second follows by differentiating the Fourier transform in $t$. All the functions are Schwartz functions on the real axis.

For completeness, the finite Fourier transform of a primitive Dirichlet character satisfiesFor a unit $r$ this is substitution by its inverse. If $p\mid (r,N)$, primitivity supplies a unit $u\equiv 1(\mathrm{mod}N/p)$ with $\chi (u)\ne 1$: otherwise the character would factor through modulus $N/p$. Multiplication by $u$ fixes the additive exponential but multiplies the sum by $\chi (u)$, so the sum vanishes. This proves the identity also on nonunits. Finite Parseval identity now gives $|g(\chi ){|}^2=N$, since exactly $\phi (N)$ frequencies have nonzero character values. In particular this Gauss sum of a Dirichlet character is nonzero.

Apply Poisson summation to each residue class: ${\sum }_{j}F(a+Nj)=N^{-1}{\sum }_r{e}^{2\pi ira/N}\hat{F}(r/N)$. Multiply by $\chi (a)$ and sum over $a$. The formulas above give the two character theta transformationsandIn the odd case the prefactor originally reads $g(\chi )/(N^2\tau \sqrt{-i\tau })$; using $\tau =i(-i\tau )$ explains the factor $i$. Opposite parity makes the corresponding series identically zero by pairing $n,-n$. This is why the weight-three-halves series is needed for odd Dirichlet characters.

To deduce the products from these transformations, use two specific primitive Dirichlet characters. The even character ${\chi }_{12}$ takes values $1,-1,-1,1$ on $1,5,7,11$ modulo twelve. It is primitive: its values do not factor through six or four, hence through no proper divisor of twelve. Its Gauss sum of a Dirichlet character is $2\sqrt{3}=\sqrt{12}$. PutThe even transformation gives $A(-1/\tau )=\sqrt{-i\tau }A(\tau )$; also $A(\tau +1)=e^{\pi i/12}A(\tau )$, since each nonzero term has $n^2\equiv 1(\mathrm{mod}24)$. Thus $A^{24}$ transforms as a weight-twelve modular form under the generators $S,T$ of the modular group. It is holomorphic and starts with $q$, so it is a normalized cusp form. The proof in Question 1 shows $S_{12}=\mathbb{C}{\Delta }_0$, giving $A^{24}={\Delta }_0$. In particular $A$ has no interior zeros, by the valence formula for the modular group.

The Serre derivative argument in Question 2 gives $D{\Delta }_0=E_2{\Delta }_0$. Hence $DA/A=E_2/24$. The analytic product $P(\tau )=q^{1/24}{\prod }_{n\ge 1}(1-q^n)$ is nonzero and has the same logarithmic derivative, since $DP/P=1/24-\sum {\sigma }_1(n)q^n$. Therefore $A/P$ is constant on the connected complex upper half-plane; its limit at infinity is one. This proves $A=P$, without assuming either of the product identities being deduced. Pair the terms $n,-n$ in $A$, choosing representatives $n=6r+1$. Since ${\chi }_{12}(6r+1)=(-1{)}^r$, we getThis is the Euler pentagonal identity, derived by the character theta proof of the eta product.

For the second product use the odd primitive character ${\chi }_4(1)=1$, ${\chi }_4(3)=-1$, with Gauss sum of a Dirichlet character $2i$. SetThe odd transformation gives $B(-1/\tau )=(-i\tau {)}^{3/2}B(\tau )$ and the series gives $B(\tau +1)=e^{\pi i/4}B(\tau )$. Thus $B^8$ is another normalized weight-twelve cusp form, and $B^8={\Delta }_0=A^{24}$. Because $A$ is nonzero, $B/A^3$ has eighth power one, so it is constant; its leading coefficient fixes it to one. Pairing positive and negative odd integers now givesSince $B=A^3=q^{1/8}{\prod }_{n\ge 1}(1-q^n{)}^3$, the Jacobi cubic product identity follows:All pairings and differentiations take place in locally uniformly convergent theta functions or products on $|q|<1$.

Use the Fourier transform convention $\hat{h}(\xi )=\int h(x)e^{-2\pi ix\cdot \xi }dx$. The deleted-ball kernel need not be integrable at infinity, so first work with finite annuli. WriteWe will prove a uniform bound, convergence as $R\to \infty$, and then convergence as $\epsilon \downarrow 0$. This makes all Fourier-transform assertions precise for the annularly cancelling singular convolution kernel.

Fix $\xi \ne 0$ and put $y=\xi /(2|\xi {|}^2)$, so $e^{-2\pi iy\cdot \xi }=-1$. For $R>r\ge 4|y|$, let $\Omega =\{r<|x|<R\}$ and defineThe Hörmander integral kernel condition givesOn changing variables $u=x-y$, the phase factor shows $J_{r,R}=-{\int }_{\Omega -y}K(u)e^{-2\pi iu\cdot \xi }du$. The symmetric difference between $\Omega$ and $\Omega -y$ lies in the two shells whose radii are within $|y|$ of $r$ and $R$. On each such shell, the kernel's size bound gives, by polar integration,Thereforeand adding the two comparisons proves the oscillatory tail estimate for a singular kernel:The integral on the right is at most $B$. Moreover, for fixed $\xi$, it tends to zero as $r\to \infty$ because the translated difference is absolutely integrable outside $2|y|$. The shell term also tends to zero. Thus the far Fourier integral is Cauchy and converges as its outer radius tends to infinity.

Set $r_0=4|y|=2/|\xi |$. If $\epsilon <r_0$, annular cancellation allows subtraction of the constant phase on the near annulus:The same estimate applies to a shorter near annulus when $R<r_0$. The far annulus is bounded by $B/2+C_{d}A$, using the tail estimate; if $\epsilon \ge r_0$, apply that estimate with $r=\epsilon$. At $\xi =0$ the annular integral vanishes by cancellation. We have provedFor each fixed $\epsilon >0$, the size bound makes $K_{\epsilon }\in L^2$, and $K{1}_{\{\epsilon <|x|<R\}}\to K_{\epsilon }$ in $L^2$ as $R\to \infty$. By the Plancherel theorem, their Fourier transforms converge in $L^2$ to ${\hat{K}}_{\epsilon }$. The pointwise improper limits just established identify this transform, up to a null set, withFor example identification follows by selecting an almost-everywhere convergent subsequence from the $L^2$ convergence. HenceThe near-annulus integral with the constant phase subtracted converges absolutely at zero, since its bound is proportional to $|x{|}^{1-d}$. The far integral has already converged. Therefore $m_{\epsilon }(\xi )\to m(\xi )$ for every $\xi \ne 0$, with the same uniform bound; the value at zero may be set arbitrarily.

For $f\in L^2$, the convolution $K_{\epsilon }\ast f$ is absolutely convergent for each $x$ by Cauchy-Schwarz inequality. Its Fourier transform is $m_{\epsilon }\hat{f}$: first prove this for Schwartz functions and then approximate any $L^2$ input. The convolutions converge uniformly under this approximation by the same Cauchy-Schwarz bound, and the multiplier outputs converge in $L^2$ by their uniform operator bound, so these two definitions agree. Define the bounded Fourier multiplier operator $T$ by $\widehat{Tf}=m\hat{f}$. Dominated convergence theorem and the Plancherel theorem now give

For the final request, use a normalized bump function, $\int \varphi =1$, so that its dilates form an approximate identity. Multiplying Fourier transforms givesThese products are integrable by Cauchy-Schwarz inequality, since $T{\varphi }_{\epsilon },f\in L^2$, and also define $L^2$ functions because $m$ and $\hat{\varphi }$ are bounded. Thus the convolution identity is justified even though $T{\varphi }_{\epsilon }$ need not be in $L^1$:The allowed standard approximate identity results, applied to $Tf\in L^2$, prove convergence to $Tf$ in norm and almost everywhere. This is the mass normalization in mollified singular integrals; it does not require proving almost-everywhere convergence of the sharp truncations $T_{\epsilon }f$.

If “bump function” means merely smooth with compact support, the stated limit needs the unit-integral qualification. For arbitrary such $\varphi$, the limit is $(\int \varphi )Tf$ instead. This is also the almost-everywhere limit: at each Lebesgue point of $Tf$, subtract $(\int \varphi )Tf(x)$ inside the convolution and bound the remaining local average by the boundedness and compact support of $\varphi$. As a concrete counterexample without normalization, take $K(x)=1/(\pi x)$ in dimension one and $\varphi =2\rho$ for a unit-integral mollifier $\rho$. This kernel has the required size, odd annular cancellation, and translation-difference integral $\log (3)/\pi$ on $\{|x|>2|y|\}$ for every $y\ne 0$. Its limit operator is the nonzero Hilbert transform; by the Hilbert-transform Fourier multiplier, it is an $L^2$ isometry. For any nonzero $f\in L^2$ the mollified limit is $2Tf$, not $Tf$.

Use the Fourier transform convention $\hat{h}(\xi )={\int }_{\mathbb{R}}h(x)e^{-2\pi ix\xi }dx$. For a Schwartz function $h$, the Poisson summation formula isTo prove it, periodize $h$ as $P(x)={\sum }_{n\in \mathbb{Z}}h(x+n)$. Rapid decay of all derivatives makes this a smooth periodic function, with uniform termwise differentiability. Its $m$th Fourier coefficient iswhere absolute convergence justifies splitting the real-line integral into unit intervals. Integration by parts gives rapid decay of the coefficients of $P$, so its Fourier series converges absolutely and uniformly to $P$. Evaluate it at zero to obtain the formula.

For $\mathrm{Im}\tau >0$ and $z\in \mathbb{C}$, take the complex Gaussian function $h(x)=e^{\pi i\tau x^2+2\pi izx}$. It is a Schwartz function of the real variable $x$. Its Fourier transform isThe real Gaussian integral proves this first for $\tau =iy$, $y>0$, and real $z$. Normal convergence of the defining integrals makes both sides entire in $z$ and holomorphic in $\tau \in \mathbb{H}$, so the identity theorem extends it to the stated domain. Choose the square-root branch holomorphic for $-i\tau$ in the right half-plane and positive when $\tau =iy$.

Applying Poisson summation gives the transformation of the Jacobi theta function:The series on the right is ${\vartheta }_{00}(z/\tau ,-1/\tau )$. Rearranging therefore givesThe branch choice is part of the formula; at $\tau =i$ its square-root factor equals one.

Use the Fourier transform conventionThis convention also gives the sign of the inverse transform in 15A. The Parseval identity, in its bilinear form, isIn particular $\int |f{|}^2=(2\pi {)}^{-1}\int |\hat{f}{|}^2$. Initially take $f,g$ to be Schwartz functions; the formula then extends to $L^2$ with Fourier transforms interpreted by completion.

Here is a proof which also justifies the needed limit operations. For $\epsilon >0$, formAbsolute integrability allows Fubini's theorem to replace the transformed factors by their defining integrals. The elementary Gaussian integral givesso $I_{\epsilon }=\int f(x)\overline{(K_{\epsilon }\ast g)(x)}dx$, with convolution $(h\ast g)(x)=\int h(x-y)g(y)dy$. The kernel is nonnegative, has integral one, and its mass outside any fixed neighborhood of zero tends to zero. Since a Schwartz function is bounded and uniformly continuous, splitting the convolution difference into small and large $|x-y|$ proves $K_{\epsilon }\ast g\to g$ uniformly. As $f\in L^1$, the last integral consequently tends to $\int f\overline{g}$.

On the spectral side, repeated integration by parts makes the transforms rapidly decreasing; their product is integrable. The dominated convergence theorem therefore gives $I_{\epsilon }\to (2\pi {)}^{-1}\int \hat{f}\overline{\hat{g}}$. Equating the two limits proves the displayed identity. For general $L^2$ functions, approximate by Schwartz functions. The norm identity makes their transforms Cauchy in $L^2$ with the prescribed factor; the limit defines the $L^2$ Fourier transform independently of the approximation. Taking limits in the bilinear identity proves it for all $L^2$ pairs. Thus the proof includes the hypotheses and avoids unjustified exchanges of nonabsolutely convergent integrals.

Take the usual Heaviside step function, $H(x)=1$ for $x>0$ and $0$ for $x<0$. Its Fourier transform must be interpreted as a tempered distribution. For $a>0$, exponential damping gives an ordinary integral:Against any Schwartz function $\phi$, the first term is an approximate identity of total mass $\pi$: changing variables $k=au$ gives the limit $\pi \phi (0)$. For the second term, pair positive and negative $k$ near zero. The difference $\phi (k)-\phi (-k)=O(k)$ removes the singularity and gives the Cauchy principal value of $1/k$. The tails converge by dominated convergence. Since the damped functions also approach $H$ as tempered distributions, continuity of the Fourier transform yields the Fourier transform of the Heaviside step function:The printed plus sign is inconsistent with the printed $e^{-ikx}$ convention for this $H$. In fact $H^{\prime }=\delta$ requires $ik\hat{H}=1$; the printed expression would give $-1$. Its plus sign would instead be correct for $H(-x)$.

To find the Fourier transform of a principal-value reciprocal, exponential damping of its odd part givesThe integral vanishes at $k=0$ and its $k$-derivative is $a/(a^2+k^2)$, so it equals $\arctan (k/a)$. Taking the distributional convergence limit givesAn ordinary improper integral of $1/x$ across zero is not a substitute for the specified Cauchy principal value.

The distributional identity ${\partial }_{\bar{z}}(1/(\pi z))=\delta$ gives the Cauchy-Green operator solutionwhere $\tilde{f}$ is $f$ in the changed real coordinates. The singularity $1/(z-z^{\prime })$ is locally integrable. For Schwartz functions, splitting the integral into a region near the origin and a rapidly decaying far region proves $\mu =O(|z{|}^{-1})$. If two solutions obey this decay, their difference has zero Wirtinger derivative and is an entire holomorphic function tending to zero; the Liouville theorem proves uniqueness.

The change-of-variables Jacobian has absolute value $a_r/|b_r|$. Therefore an equivalent form of the complexified transport Cauchy kernel isIn the original complex coordinate $w=x_1+i{x}_2$, it isFor $|\lambda |<1$ the denominator can be rewritten to give the kernel $\lambda /[\pi ((w-w^{\prime })-{\lambda }^2(\bar{w}-{\bar{w}}^{\prime }))]$. For $|\lambda |>1$ it is ${\lambda }^{-1}/[\pi ((\bar{w}-{\bar{w}}^{\prime })-{\lambda }^{-2}(w-w^{\prime }))]$. Consequently the spectral solution is holomorphic separately inside and outside the unit circle, with $\mu =O(\lambda )$ at zero and $\mu =O({\lambda }^{-1})$ at infinity. These normalizations will determine its jump reconstruction.

For positive time, convolution with the heat kernel givesDifferentiation under the integral verifies the heat equation, and the approximate identity property gives the initial value. Taking absolute values under the integral and using $e^{-|x-y{|}^2/(4t)}\le 1$ givesThe printed inequality omits the modulus even though $g$ is complex-valued. A complex number has no compatible order, so the displayed absolute-value estimate is the meaningful correction.

Substitution of $u=e^{it}v$ reduces the forced partial differential equation to $(i-\Delta )v=f$. With the Fourier transform convention for which $\widehat{\Delta v}=-|\xi {|}^2\hat{v}$, setThe denominator never vanishes for real $\xi$. Every derivative of its reciprocal has at most polynomial growth, so multiplying a Schwartz function by it gives another Schwartz function. The inverse Fourier transform therefore supplies the required $v$.

For the prescribed initial value the full solution is $u(t,x)=e^{it}v(x)+K_t\ast (g-v)(x)$. Applying the heat kernel estimate to the Schwartz function $g-v$ proves the stronger uniform conclusion

Use $\hat{f}(\lambda )=\int f(t)e^{-i\lambda t}dt$ and interpret the regularity assumption to include $f\in L^2(\mathbb{R})$ with the continuous representative furnished by Fourier inversion. Write $F=\hat{f}$. The bandwidth assumption implies $F\in L^1$ as well as $L^2$, by the Cauchy-Schwarz inequality on its finite support, sois continuous and defined at every $t$. The integer samples are precisely the Fourier coefficients of $F$ with the sign convention appropriate to the basis $e^{-in\lambda }$:The completeness and Parseval identity for the exponential orthonormal basis on $[-\pi ,\pi ]$ give $F={\sum }_{n\in \mathbb{Z}}f(n)e^{-in\lambda }$ in $L^2$ andThe last equality is the Plancherel theorem. Apply the inverse integral to the finite symmetric Fourier partial sums. Sincewe obtain the Nyquist–Shannon sampling theorem in this normalization:where the quotient is one at $t=n$. This is a sampling expansion by periodic Fourier projection, not merely an almost-everywhere inversion. Indeed, the inverse integral of an error $G\in L^2[-\pi ,\pi ]$ is bounded, uniformly in $t$, by $\Vert G{\Vert }_2/\sqrt{2\pi }$. Hence the symmetric finite sampling sums converge uniformly on all of $\mathbb{R}$. Moreover, the Parseval identity applied to $e^{it\lambda }$ givesThe Cauchy-Schwarz inequality therefore makes the sampling series absolutely convergent at every $t$, with its absolute tail bounded uniformly by $({\sum }_{|n|>N}|f(n){|}^2{)}^{1/2}$. Thus the integer samples recover the function everywhere.