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]

If the constant symmetric matrix $A$ satisfies $A\xi \cdot \xi \ge \theta |\xi {|}^2$, then every smooth function $u$ with compact support satisfiesIndeed, the Fourier transform of a derivative and the Plancherel theorem turn the squared norms into integrals with Fourier multipliers $|\xi {|}^4$ and $(A\xi \cdot \xi {)}^2$.

For a unit vector $v$, the directional Hilbert transform integrates along lines parallel to $v$. Its Fourier multiplier is $-i\mathrm{sgn}(\xi \cdot v)$ under the usual transform convention. Applied to a rectangle whose long axis is $v$, it reduces to the one-dimensional logarithmic tail of an interval indicator function, separately on each transverse coordinate.

For the upward derivative at a flat boundary, the Dirichlet-to-Neumann map for a Helmholtz half-space is the Fourier multiplier $i\beta (q)$. The outward normal of the upper half-space points downward and has the opposite sign. Its evanescent multiplier grows like $|q|$, explaining why fine roughness matters even when the height is small.

The Radon transform records line integrals of a planar density. Apply the transverse Hilbert transform and a derivative to each projection, then average the filtered values over all lines through the reconstruction point. With $Hh(\rho )={\pi }^{-1}\mathrm{PV}\int h(r)/(\rho -r)dr$, the filter has Fourier multiplier $|\kappa |$. The Fourier inversion theorem, or the jump of the complexified transport Cauchy kernel, proves the displayed formula. The factor is $1/(2\pi )$ when directions are integrated over a half-circle. The growing Fourier multiplier explains amplification of measurement noise and the need for regularization of an inverse problem.

For a polynomial with no sufficiently large real zero, a compact low-frequency cutoff function removes its remaining zeros. A reciprocal satisfying $|{\partial }^{\alpha }b|\lesssim ⟨\xi {⟩}^{-m-\delta |\alpha |}$, $m\ge 0$, $0<\delta \le 1$, defines a Fourier multiplier mapping $H^s$ to $H^{s+m}$. Its inverse-transform kernel is smooth off zero: repeated frequency integration by parts supplies arbitrary decay even after multiplying by powers of frequency for spatial derivatives. Consequently a compactly supported distribution supported away from the target contributes only a smooth function there. In a local regularity proof, the cutoff function commutator is placed in that separated region rather than estimated only by its differential order.

The two boundary values have jumpThe line-tail contributions cancel. Because the complexified transport Cauchy kernel is holomorphic on both sides of the unit circle and vanishes at spectral infinity, the Cauchy integral formula reconstructs it from this jump:There is no additional entire term: its inside and outside boundary values would agree, so it would define an entire function vanishing at infinity. Expand outside the circle as $\mu ={\mu }_1/\lambda +O({\lambda }^{-2})$. ThenIn the original coordinate $w=x_1+i{x}_2$, the transport equation is $(\lambda {\partial }_w+{\lambda }^{-1}{\partial }_{\bar{w}})\mu =f$. Its constant term at infinity is ${\partial }_w{\mu }_1=f$. Sincewe obtain the inverse Radon transform in filtered backprojection form:Opposite directions represent the same unoriented line. Reversing $\theta$ by $\pi$ reverses $\rho$; the Hilbert transform changes sign under this reversal, and its derivative changes it back. Thus the integrand agrees for opposite directions and the factor becomes $1/(2\pi )$ if the angular integral is over $[0,\pi )$. With the convention $\hat{h}(\kappa )=\int e^{-i\kappa \rho }h(\rho )d\rho$, the Hilbert transform has Fourier multiplier $-i\mathrm{sgn}\kappa$; the derivative makes the filter $|\kappa |$. This also identifies the formula with the standard ramp-filter reconstruction.

In ideal computed tomography, the unknown is an absorption coefficient $a(x)$. Along a straight ray the transmitted intensity satisfies $dI/d\tau =-aI$, soThus calibrated transmission measurements supply the Radon transform of the slice. The inversion above reconstructs the coefficient by filtered backprojection: filter each projection in its transverse variable and then accumulate it over the rays through each image point. This turns a collection of line-integral measurements into a spatial map, which is the key usefulness of the result.

The ideal derivation assumes known ray geometry, sufficiently complete angular sampling and a straight-ray attenuation model. The ramp Fourier multiplier grows at high frequency, so measurement noise is amplified. Practical reconstruction therefore uses discretization and regularization of an inverse problem, for example a window on the ramp filter. Limited angles or missing detector positions cannot be treated as if they were complete noiseless Radon transform data.

Use positive carrier wavenumber $k$ and retain waves traveling mainly in the positive $x$ direction. Factor out their rapid wave phase by writing $\psi =e^{ikx}E(x,z)$. Substitution into the Helmholtz equation is exact and givesFor a transverse Fourier mode $e^{i\nu z}$ at a small angle, $|\nu |/k\ll 1$ and its forward longitudinal wavenumber isThus the envelope varies along $x$ on scale $k/{\nu }^2$, much more slowly than the carrier wavelength, and $|E_{xx}|/|k{E}_x|=O({\nu }^2/k^2)$. Dropping this smaller term derives the parabolic wave equation,With transform $\hat{E}(x,\nu )=\int E(x,z)e^{-i\nu z}dz$, its free propagation isThis is a one-way paraxial approximation, not a replacement for backward waves or large-angle Helmholtz components. Its Fourier multiplier has unit modulus and therefore does not attenuate any admitted component.

Choose a spatial interval with periodic boundary conditions and a uniform grid, so no separate inflow boundary closure of a difference scheme is required. The discrete Fourier transform is a unitary operator in the discrete L2 norm. Every circulant stencil is diagonal in that transform, and a one-step Fourier multiplier is a contraction exactly when its modulus is at most one for every grid frequency. For a mesh-uniform statement, use all $\theta \in [-\pi ,\pi ]$, the limiting frequency range.

A Fourier mode $e^{im\theta }$ giveswith the coefficients from (a). Squaring moduli and putting $v=\cos \theta$ yieldsThis has the required nonnegative sign for every frequency exactly whenOne must also check that the implicit equation is solvable, rather than divide by a vanishing symbol. HereFor $\mu \ne 1/2$, a zero would have to occur at zero or $\pi$; the latter occurs only at $(1\pm \sqrt{3})/2$, both outside the displayed stable intervals. At $\mu =1/2$, $D=3/4+(1/4)\cos \theta \ge 1/2$. Thus there is no denominator zero in the stable range. For fixed such $\mu$, continuity on the compact frequency interval gives a uniformly bounded inverse, and $|G|\le 1$ proves $\Vert U^n{\Vert }_{2,h}\le \Vert U^0{\Vert }_{2,h}$ for all steps.

Outside these intervals the displayed difference is negative at every nonzero frequency where $D\ne 0$, so $|G|>1$. On periodic refining meshes there are frequencies approaching any selected such value. Their amplification is bounded away from one and is repeated $O(h^{-1})$ times over a fixed physical interval when $k=\mu h$, violating a mesh-uniform stability bound. At the two singular parameters the implicit operator itself fails on meshes containing the Nyquist mode. These observations prove necessity as well as sufficiency.

For the physical convention $k>0$, $h>0$, the result is $0<\mu \le 1$ or $\mu \ge 2$, with zero admitted as the identity limit. The negative intervals describe the algebraic contraction result for signed Courant numbers. In particular, the additional large positive interval must not be discarded by imposing the usual explicit-scheme Courant–Friedrichs–Lewy condition: this scheme is implicit. Stability at large fixed Courant number is not a guarantee of a small error constant.

Interpret $H_{-\infty }({\mathbb{T}}^n)={\bigcup }_{s\in \mathbb{R}}{H}^s({\mathbb{T}}^n)$, the periodic distributions of finite negative Sobolev space order. Smooth multiplication is bounded on every periodic Sobolev space. One direct proof uses the rapid decay of the Fourier coefficients of $V$: the inequalityreduces the weighted convolution formula for $\widehat{Vf}$ to an ${\ell }^1\ast {\ell }^2$ estimate. Thus $\Vert Vf{\Vert }_{H^s}\le C_{V,s}\Vert f{\Vert }_{H^s}$ for every real $s$.

Suppose initially $f\in H^s$. The equation gives $\Delta f=u-Vf\in H^s$, since a smooth function belongs to every periodic Sobolev space. For $k\ne 0$, the Fourier multiplier of $\Delta$ is $|k{|}^2$, andThe zero coefficient is already finite. Hence $f\in H^{s+2}$. Iterating the argument gives $f\in H^{s+2j}$ for every $j\ge 0$. Applying the Sobolev embedding theorem also to each distributional partial derivative produces a smooth representative. Therefore $\boxed{f\in C^{\infty }({\mathbb{T}}^n)}$. This is an explicit elliptic regularity bootstrap, valid for complex or real $V$ and without a sign assumption on it.

Use $D_j=-i{\partial }_{x_j}$ and $⟨\xi ⟩=(1+|\xi {|}^2{)}^{1/2}$, so the Fourier transform of $P(D)u$ is $P(\xi )\hat{u}(\xi )$. Write $P_N$ for the homogeneous degree-$N$ part. An elliptic differential operator has $P_N(\xi )\ne 0$ for every real $\xi \ne 0$. By continuity on the unit sphere, $m={\min }_{|\omega |=1}|P_N(\omega )|>0$. The lower-degree terms are bounded by $C|\xi {|}^{N-1}$ for $|\xi |\ge 1$, whencefor sufficiently large $|\xi |$. This proves the high-frequency lower bound for an elliptic polynomial. Degree zero just means a nonzero constant and has no derivative gain.

With the Fourier convention in Question 3, the Sobolev space isThe condition includes that the weighted transform is represented by an $L^2$ function. For open $X$, the Local Sobolev space consists of $u\in {\mathcal{D}}^{\prime }(X)$ for which $\chi u$, extended by zero, belongs to $H^s({\mathbb{R}}^n)$ for every $\chi \in C_c^{\infty }(X)$. We use the elementary Sobolev multiplication by a smooth cutoff fact for every real $s$. It follows from Fourier convolution with the rapidly decreasing $\hat{\chi }$, the weighted inequality $⟨\xi {⟩}^s\lesssim ⟨\xi -\eta {⟩}^{|s|}⟨\eta {⟩}^s$, and the $L^1\ast L^2\to L^2$ convolution bound.

For a compactly supported distribution, continuity on test functions supported in a fixed compact neighborhood gives finite order: for some integer $M$,Insert a compact smooth cutoff function equal to one near the support, times $e^{-ix\cdot \xi }$. The resulting Fourier transform is smooth and satisfies $|\hat{u}(\xi )|\le C^{\prime }⟨\xi {⟩}^M$. The distribution is also tempered by this same finite-order bound. ThusIndeed the square of the weighted bound is integrable exactly when $2(t+M)<-n$. This is negative Sobolev regularity of a compactly supported distribution.

We next prove local regularity by a high-frequency reciprocal parametrix kernel, including its off-diagonal smoothing property. Choose $\psi \in C_c^{\infty }({\mathbb{R}}^n)$ equal to one on a ball containing all real zeros of $P$, and put $b(\xi )=(1-\psi (\xi ))/P(\xi )$, defining it smoothly as zero in the inner ball. Differentiation of the reciprocal and the elliptic lower bound giveThe Fourier multiplier $E=b(D)$ maps $H^s$ to $H^{s+N}$ by its zeroth-order bound, andThe kernel $K={\mathcal{F}}^{-1}b$ is smooth away from the origin: for $x\ne 0$, repeatedly integrate by parts using $e^{ix\cdot \xi }=(i|x{|}^2{)}^{-1}x\cdot {\partial }_{\xi }{e}^{ix\cdot \xi }$. After sufficiently many integrations, the differentiated symbol is integrable. For any desired $x$ derivative, repeat the argument with the additional polynomial ${\xi }^{\beta }$. A large-radius cutoff function justifies each step and its removal uniformly on compact sets away from zero. Thus convolution by $K$ carries a compactly supported distribution to a smooth function at points separated from its support. Also $R$ carries a compactly supported distribution to a Schwartz function, because $\psi \hat{u}$ is smooth with compact support.

For a target compact set in $X$, choose $\chi \in C_c^{\infty }(X)$ equal to one on a neighborhood of it. Then $\chi u$ and $[P(D),\chi ]u$ have compact support, and the commutator is supported where derivatives of $\chi$ occur, away from the target. The parametrix identity givesThe first term is in $H^{s+N}$ since $\chi P(D)u\in H^s$. The other two terms are smooth near the target by the proved kernel property. Since the target was arbitrary,This proof does not discard the cutoff function commutator; it places its support away from the set where regularity is sought.

For the final polynomial, select a multi-index ${\alpha }_0$ with $|{\alpha }_0|=N$ and ${\partial }^{{\alpha }_0}Q$ a nonzero constant. The derivative-ratio hypothesis immediately gives $|Q(\xi )|\ge c|\xi {|}^{\delta N}$ at large frequency. In particular $Q$ has no real zero there. For $N>0$, its degree bound also forces $\delta \le 1$. Differentiating $1/Q$ gives products of ratios ${\partial }^{\beta }Q/Q$, whose total derivative order is $|\alpha |$. Hence the corresponding high-frequency reciprocal satisfiesIts multiplier maps $H^s$ to $H^{s+\delta N}$. Its kernel is still smooth off zero: each frequency integration by parts now lowers the order by $\delta$, so more iterations may be needed, but $\delta >0$ supplies arbitrarily much decay. The same separated-support commutator identity applies without a change. The derivative-ratio Sobolev gain for a polynomial operator is thereforeSimply estimating the commutator by its differential order would lose this sharper gain; the off-diagonal kernel argument uses the full derivative-ratio hypothesis.

We give a quantitative construction for the whole open set, including the case of infinitely many component intervals. Pointwise divergence at a single point is insufficient here. Write $z=e^{2\pi it}$ and use normalized measure on $\mathbb{T}$.

First let $K$ be a finite union of closed arcs with $|K|\le 3\epsilon$, where $\epsilon$ is small. Choose a smooth function nonnegative function $u$, equal to one on $K$, with $u\le 1$ and $\eta :=\int u\le 6\epsilon$. This follows by making the transition regions around the finitely many arc endpoints sufficiently short. Define the analytic Herglotz integralIts real part is $(\eta +P[u](z))/(2\eta )>0$, where $P[u]$ is the disk Poisson integral, and $G(0)=1$. Smoothness of $u$ makes $G$ continuous and smooth function up to the circle; on $K$ its boundary real part is at least $1/(2\eta )$. The branch $H=\log G$ therefore has $H(0)=0$ andThe complex logarithm has no branch problem because $G$ lies in the right half-plane. Its smooth function boundary values give absolutely summable Taylor coefficients. Choose an analytic polynomial $A$, with zero constant coefficient, uniformly within $1/4$ of $H$ on the closed disk. Put $P(t)=\mathrm{Im}A(e^{2\pi it})/\pi$. Then $\Vert P{\Vert }_{\infty }\le 1/2+1/(4\pi )<1$.

If $D=\mathrm{deg}A$ and $m>D$, the modulated trigonometric polynomial $Q(t)=e^{2\pi imt}P(t)$ has spectrum in $[m-D,m+D]$. A sharp Fourier partial sum at $m$ selects exactly the negative-frequency half of $P$, since $A(0)=0$:Consequently, for every $t\in K$,The carrier frequency $m$ can be arbitrarily large. This compact-set Fourier amplification lemma uses a bounded imaginary part to obtain a large one-sided Fourier multiplier; no assertion about pointwise divergence is used.

Now let $E$ be any open set of measure $\delta$. DefineIf this is zero take $f=0$. Otherwise choose a decreasing sequence ${\delta }_j\downarrow 0$, with ${\delta }_1=\delta$, so small for $j\ge 2$ thatFor $j=1$ this inequality also holds: $B(\delta )>4$ and $B(\delta )/2\ge B(\delta )/4+1$. Let $L_j$ be an increasing compact exhaustion of $E$.

Using compact batching of a small open set, we recursively choose finite closed-arc unions $K_j\subset E$, of measure at most $3{\delta }_j$, such that $L_j$ is covered by the interiors of $K_1,\dots ,K_j$ andHere is the measure-theoretic justification. At step $j$, the uncovered portion has measure at most ${\delta }_j$. Its intersection with $L_j$ after removing the previously covered interiors is compact, and can be covered by finitely many closed arcs lying in $E$ with total length at most ${\delta }_j$ plus an arbitrarily small excess. By inner regularity of Lebesgue measure, a further compact subset of the uncovered open portion captures all but ${\delta }_{j+1}$ of its measure; cover it similarly. The combined finite arc union has length at most $3{\delta }_j$. This both covers the outstanding part of $L_j$ and ensures the next remainder estimate. In particular $E\subset {\bigcup }_j{K}_j$.

Apply the finite-arc construction to $K_j$ to get $P_j$ of degree $D_j$. Choose carriers $m_j$ recursively so that the bands $[m_j-D_j,m_j+D_j]$ are positive and strictly ordered with no overlap. SetThe frequency-separated Fourier block series and the Weierstrass M-test give a continuous complex-valued function with $\Vert f{\Vert }_{\infty }\le {\sum }_j{2}^{-j}=1$. At the cutoff $m_j$, all earlier bands have been included in full and every later band is absent. Thus on $K_j$,Every point of $E$ lies in some $K_j$, soValues of $\psi$ outside $0<\delta \le 1$ may be defined arbitrarily. The construction handles an arbitrary open set, not just one point or a finite collection of intervals. It also prevents cancellation explicitly: earlier bands contribute a bounded complete sum, and later bands contribute nothing at the chosen cutoff.

For the second-order rough-surface scattered field, the Dirichlet boundary condition expanded at the mean plane isAt normal incidence, ${\psi }^{[0]}=e^{-ikz}-e^{ikz}$, so its second normal derivative vanishes at zero. Define the Dirichlet-to-Neumann map for a Helmholtz half-space through the Fourier multiplier $\mathcal{B}$:The first-order trace is $2ikh$, hence ${\partial }_z{\psi }_s^{[1]}(x,0)=-2k\mathcal{B}h(x)$. It follows thatwhere $h=O(\epsilon )$ with fixed regular profile. In integral notation the quadratic contribution isThe second-order field above the mean plane is $\mathcal{E}[2kh\mathcal{B}h]$ added to $-e^{ikz}+\mathcal{E}[2ikh]$. No local replacement of $\beta (q)$ by $k$ has been made; such a replacement would be an additional long-spatial-scale approximation.

The physical surface trace and reference-plane trace are distinct. At the actual rough boundary $z=h(x)$, the condition itself says ${\psi }_s(x,h(x))=-e^{-ikh(x)}$, soThe first boxed expression is the reference-plane trace needed in part (d), at $z=0$, using the perturbative continuation where that plane lies below the actual boundary; the second answers the literal “at the surface” wording if it means the physical boundary. Taylor-expanding the first expression and its normal derivatives from $z=0$ to $z=h$ reproduces the second, so there is no contradiction between them.

Let the stationary height have covariance function $C_h(\xi )=⟨h(x)h(x+\xi )⟩$, with $C_h(0)={\sigma }^2$. Define the power spectrum of surface height byFor real stationary heights this spectrum is even and nonnegative. The Fourier multiplier identity in part (c) givesThus coherent reflection from a stationary rough surface at normal incidence isRequire the corresponding weighted spectral moment to exist. If the stationary process has a spectral measure rather than a density, the same formula uses that measure with the matching normalization.

Splitting the propagating and evanescent parts makes the effect clear:through second order. The positive real correction reduces the magnitude of the initially negative unit coherent reflection at this order, as some reflection becomes diffuse. The evanescent part produces a coherent wave phase correction. In contrast, the first-order mean was exactly the flat reflected wave.

The height-correlation dependence of coherent reflection cannot generally be determined from $\sigma$ alone: the quadratic term weights the whole spectrum by $\beta (q)$, while ${\sigma }^2=(2\pi {)}^{-1}\int S_h(q)dq$ is unweighted. If the roughness varies only on scales much longer than the wavelength, so its spectrum is concentrated at $|q|\ll k$, then $\beta \simeq k$ andThis is a useful limiting formula, not the general second-order answer under only small-height assumptions. Also the mean field sampled at the moving physical boundary is $-1+k^2{\sigma }^2/2$ through second order, from part (c), and is a different observable.

For smooth periodic initial data, Fourier expansion reduces the equation to $\dot{\hat{u}}(m,t)=-(1+m^2)\hat{u}(m,t)$. Therefore define the Fourier multiplier solution operators byTheir Sobolev norm satisfies $\Vert S(t)u_0{\Vert }_{H^s}^2\le e^{-2t}\Vert u_0{\Vert }_{H^s}^2\le \Vert u_0{\Vert }_{H^s}^2$, so they extend to contractive linear operators on $H_{\mathrm{per}}^s$. Multiplication of the Fourier factors gives $S(t+s)=S(t)S(s)$ and $S(0)=I$. Finallyby dominated convergence theorem, dominated by four times the summable norm terms. This proves a strongly continuous contraction semigroup for every requested $s$; the semigroup property extends continuity from zero to other times. For smooth data, termwise derivatives converge and recover the original Cauchy problem for a partial differential equation uniquely.

For the heat equation, viewed as a second-order equation in all $(t,x)$ variables, the principal symbol is $|\xi {|}^2$. It vanishes on $dt$, so $t=0$ is a characteristic hypersurface andWriting $u_t={\Delta }_{x}u$ does not repair its hypotheses: the right side has second spatial derivatives, exceeding the first-order normal-form allowance. Indeed real analytic initial data need not produce a jointly real analytic solution through time zero. In one dimension, $u_0(x)=(1+x^2{)}^{-1}$ would forceThe time Taylor coefficients $(-1{)}^k(2k)!/k!$ have zero radius of convergence.

Nevertheless the forward heat equation has Hadamard well-posedness in standard $L^2$ or Sobolev space classes for $t\ge 0$; its Fourier multiplier $e^{-t|\xi {|}^2}$ gives stable smoothing. Backward evolution has the growing multiplier $e^{|t||\xi {|}^2}$ and is unstable. Neither forward stability nor backward instability is decided simply by whether the real analytic theorem applies.

Taking the Fourier transform of the equation in the sense of distributions givesThe Fourier multiplieris everywhere positive. Define $\hat{u}=\hat{f}/m$. Sincewe haveThus $u\in H^{s+4}$ exists and obeys the required estimate. If two solutions existed, their difference would have $m\hat{u}=0$; positivity of $m$ proves uniqueness.

A classical solution needs continuous derivatives through order four. Applying the Sobolev embedding theorem with $k=4$ to $u\in H^{s+4}$ shows that this is guaranteed whenor equivalently