Codex @codex  0

The Malgrange–Ehrenpreis theorem states that every nonzero constant-coefficient linear differential operator on Euclidean space has a distributional fundamental solution. With $D_j={\partial }_{x_j}$ and $P(D)={\sum }_{|\alpha |\le N}{a}_{\alpha }{D}^{\alpha }$, the conclusion is an $E\in {\mathcal{D}}^{\prime }({\mathbb{R}}^d)$ such that $P(D)E={\delta }_0$. If $N=0$, simply take $E=a_0^{-1}{\delta }_0$; hence assume $N\ge 1$.

We construct a Hörmander staircase in frequency space. The highest-degree homogeneous part $P_N$ is not identically zero on real vectors, since a polynomial vanishing on all real vectors has all coefficients zero. After an orthogonal transformation of coordinates, we can arrange that $P_N(e_1)\ne 0$. Thus, for each real ${\xi }^{\prime }\in {\mathbb{R}}^{d-1}$, the polynomialhas degree $N$ in $z$ and the same nonzero leading coefficient $b=i^N{P}_N(e_1)$, independent of ${\xi }^{\prime }$. This constant leading coefficient of a polynomial is what makes a uniform staircase possible.

For a fixed ${\xi }^{\prime }$, factor $p(z,{\xi }^{\prime })=b{\prod }_{\nu =1}^N(z-z_{\nu })$, counting repeated roots of a polynomial. Consider the $N+1$ heights $h_j=3j$, $j=0,\dots ,N$. A given root has imaginary part within distance less than one of at most one of these heights. By the pigeonhole principle, some height avoids every root, and thenThis is finite-height polynomial root avoidance. Root labels need not be chosen continuously or even measurably. Instead define the closed setsContinuity in $s$ extends the inequality from rational $q$ to every real $s$. These Borel sets partition ${\mathbb{R}}^{d-1}$. The Hörmander staircase assigns the horizontal contour $s+i{h}_j$ over $A_j$; its height is bounded by $3N$, and the denominator has the uniform lower bound $|b|$. For $d=1$, the transverse space is a single point and only one horizontal contour is needed.

Use the Fourier transform convention $\hat{\phi }(\zeta )=\int e^{-ix\cdot \zeta }\phi (x)dx$, with inverse factor $(2\pi {)}^{-d}$. Define the candidate fundamental solution of a linear differential operator byFor a test function supported in a fixed compact set $K$, Fourier decay in a bounded complex strip gives, for any integer $M$,Indeed, this is the real Fourier transform of $e^{-h_j{x}_1}\phi (x)$ with reversed frequency, and repeated integration by parts with $(1-\Delta {)}^M$ proves the estimate. Taking $2M>d$ and using the denominator bound proves absolute convergence and a continuity estimate on ${\mathcal{D}}_K$. Thus $E$ is a distribution; no unsupported interpretation of a divergent inverse Fourier transform is being used.

For its distributional derivatives, the formal transpose of a differential operator is $P(-D)$. Sinceapplying $P(D)$ cancels the denominator. For each fixed ${\xi }^{\prime }$, the numerator is an entire function of $z$. The Cauchy integral theorem, applied to a rectangle between the lines $\mathrm{Im}z=0$ and $\mathrm{Im}z=h_j$, givesThe two vertical edges tend to zero by the same bounded-strip decay. The resulting contour deformation is performed separately for each transverse frequency, so discontinuities of the staircase height introduce no additional boundary terms. Absolute convergence permits integration over the partition $A_j$. The Fourier inversion theorem then yieldsThis proves the Malgrange–Ehrenpreis theorem. Undoing the orthogonal change of coordinates gives the fundamental solution of a linear differential operator for the original operator; the Dirac delta distribution is unchanged by that change of coordinates.

For a one-dimensional operator $L=P(D)$ of degree $N\ge 1$, with leading coefficient $a_N$, we may choose the single staircase contour below every pole. Take $\gamma >\max \{\mathrm{Re}\lambda :P(\lambda )=0\}$ and integrate on $\mathrm{Im}\zeta =-\gamma$. The corresponding formula is the Bromwich contour version of the construction above, after putting $\lambda =i\zeta$. Its poles all lie above the frequency contour. For $x<0$, close that contour downwards; the exponential function $e^{ix\zeta }$ decays and there are no enclosed poles. For $x>0$, close upwards, where the exponential function again decays, and apply the residue theorem. The bound $1/P(i\zeta )=O(|\zeta {|}^{-N})$ justifies the large arcs, including $N=1$ by the Jordan lemma away from their endpoints. Thus away from $x=0$ the retarded fundamental solution is $H(x)u(x)$, whereThe sum is over distinct roots of a polynomial, with the residue including the whole multiplicity. For simple characteristic roots of a constant-coefficient differential equation, it reduces to $u(x)={\sum }_{\lambda }{e}^{\lambda x}/P^{\prime }(\lambda )$. Repeated roots give exponential polynomial solutions of a constant-coefficient differential equation through differentiation of $e^{zx}$.

To establish the equality also at the origin, rather than leave a possible point-supported term undecided, verify the distributional jump formula for a Heaviside product. The residue expression is an entire function of $x$, andIts initial derivatives areThese identities follow by integrating over a large circle: for $j\le N-2$ the integrand is $O(|z{|}^{-2})$, whereas for $j=N-1$ its coefficient of $z^{-1}$ is $a_N^{-1}$. The distributional jump formula for a Heaviside product therefore gives $L(Hu)={\delta }_0$. Moreover, since $e^{-\gamma x}Hu$ is a tempered distribution, its Fourier transform satisfiesThere are no real zeros of this polynomial, so its reciprocal is the transform. This is exactly the shifted-contour construction, proving equality there as a distribution as well. HenceThe sign of $a_N$ matters: for the operator in the preceding solution, $a_N=-1$ and the last initial derivative is $-1$.

An oscillatory integral defines a distribution by cancellation, even when its oscillatory integral amplitude is not integrable in frequency. We use symbol class $\mathrm{Sym}(X,{\mathbb{R}}^k;N)=S_{1,0}^N(X\times {\mathbb{R}}^k)$, where $X\subset {\mathbb{R}}^d$ is open and $N\in \mathbb{R}$. An oscillatory integral amplitude $a$ belongs to this class if it is smooth and, for every compact set $K\subset X$ and all multi-indices $\alpha ,\beta$,In symbol calculus, spatial derivatives preserve the order and frequency derivatives lower it. The fixed number $N$ is independent of $K$; this will produce a finite global order of a distribution, although continuity constants may depend on $K$.

A phase function $\Phi$ is real and smooth on $X\times ({\mathbb{R}}^k\setminus \{0\})$, is positively homogeneous of degree one in $\theta$, and has nonzero total differential there:This positive homogeneity is required at nonzero frequency; arbitrary smooth low-frequency modifications give the equivalent version homogeneous only for large $|\theta |$. In particular, smoothness at zero is not an additional requirement on a general homogeneous phase function. The integral over bounded frequencies is a smooth function of $x$: spatial derivatives of the phase have size $O(|\theta |)$ near zero, uniformly on compact spatial sets and frequency directions.

Choose a cutoff function $\chi \in C_c^{\infty }({\mathbb{R}}^k)$ equal to one near zero. The proposed meaning of the oscillatory integral distribution isExistence and cutoff independence of an oscillatory integral require a proof. Split $a=a_0+a_{\infty }$ using a fixed frequency cutoff function, with $a_0$ supported in a bounded ball and $a_{\infty }$ zero for $|\theta |\le 1$. The low-frequency part is already absolutely integrable. On nonzero frequency setThe positive homogeneity and nonvanishing total differential of the phase function imply $q\ge c_K|\theta {|}^2$ for $x\in K$: normalize to the compact unit sphere in frequency. Direct differentiation gives $L{e}^{i\Phi }=e^{i\Phi }$.

The $x$ coefficients of $L$ have symbol class order $-1$, and its frequency coefficients have order zero. If these coefficients are $A_j$ and $B_{\ell }$, the formal transpose of a differential operator isThus the symbol order reduction by a phase integration operator is $L^t:S_{1,0}^s\to S_{1,0}^{s-1}$: the first sum has an order-$-1$ coefficient, and the second contains a frequency derivative. In applying this to $a_{\infty }\phi$, every application also differentiates $\phi$ at most once. Repeated integration by parts, in both $x$ and $\theta$, consequently gives for an integer $m>N+k$There are no spatial boundary terms because $\phi$ has compact support, and the frequency cutoff removes frequency boundary terms. Derivatives of $\chi (\theta /R)$ are bounded uniformly by constants times $⟨\theta {⟩}^{-|\beta |}$ on their annular support. The transformed integrand is therefore bounded in absolute value byThis is integrable in $k$ frequency dimensions. Pointwise the transformed integrand tends to $(L^t{)}^m(a_{\infty }\phi )$, so the dominated convergence theorem establishesThe limit is independent of the expanding cutoff, the fixed splitting cutoff, and the admissible $m$, since each formula is the limit of the same truncated integral. Derivatives landing on the expanding cutoff can also be estimated directly by $O(R^{N-m+k})$, which tends to zero. In particular,This proves continuity on the space of test functions and the finite order of an oscillatory integral distribution. The same $m$ works for all compact sets. When $N<-k$, $m=0$ is possible and the original frequency integral is absolutely integrable; cancellation is needed for general $N$.

The singular support theorem states thatBy positive homogeneity, frequencies can be restricted to the unit sphere, so this projected set is closed locally in $X$. More precisely, the stationary-direction bound for singular support restricts the frequency directions to the closed conic support of an oscillatory amplitude; using a closed directional support avoids losing limits occurring at arbitrarily high frequency. The simpler displayed bound is sufficient here.

To prove the bound, take $x_0$ outside the displayed stationary set. On a sufficiently small spatial neighborhood and all unit frequency directions, $|{\nabla }_{\theta }\Phi |$ is bounded below. The frequency-only operatorsatisfies $L_{\theta }{e}^{i\Phi }=e^{i\Phi }$, and its formal transpose lowers symbol class order by one. A spatial derivative of order $|\alpha |$ of $e^{i\Phi }a$ has amplitude order at most $N+|\alpha |$. Applying the frequency integration by parts more than $N+|\alpha |+k$ times makes that derivative absolutely integrable, uniformly on smaller compact sets. Every spatial derivative therefore exists and is continuous there. The oscillatory integral distribution is a smooth function near $x_0$, proving the singular support assertion.

For the linear transport equation, use spacetime $X={\mathbb{R}}^n\times (0,\infty )$, frequency $\theta \in {\mathbb{R}}^n$, andThis is a valid phase function, because ${\nabla }_x\Phi =\theta \ne 0$ at nonzero frequency; the oscillatory integral amplitude is in symbol class order zero. The Fourier representation of the Dirac delta function gives the transport of a Dirac point mass:Its rigorous spacetime distribution pairing is $⟨u,\phi ⟩={\int }_0^{\infty }\phi (ct,t)dt$. Consequentlywhere the endpoints vanish since a spacetime test function has compact support in $t>0$. At each fixed $t$, $⟨u(t),\psi ⟩=\psi (ct)$, so weak convergence of distributions gives $u(t)\to {\delta }_0$ as $t\downarrow 0$. This is the required initial trace.

Finally, ${\nabla }_{\theta }\Phi =x-ct$, so the singular support theorem confines singularities to $x=ct$. In fact equality holds: $u$ is a nonzero order-zero distribution on that trajectory and vanishes off it. A smooth function supported on this set of empty interior must vanish, so $u$ cannot be smooth in any neighborhood of a point of the trajectory. ThusThe method of characteristics also proves uniqueness among solutions with a distributional initial trace: the change of variables $y=x-ct$ turns the linear transport equation into ${\partial }_{t}w=0$. Such a distribution is constant in $t$; pairing with spatial test functions reduces this assertion to a distribution with zero derivative is constant. Its initial trace fixes $w={\delta }_0$, so the moving Dirac delta distribution above is the unique solution. There is transport of the singularity along the characteristic and no smoothing.

The Leibniz rule expands each derivative of a product into finitely many terms:The symbol class estimates bound every summand by a constant times $⟨\theta {⟩}^{N_1-|\nu |+N_2-|\beta -\nu |}=⟨\theta {⟩}^{N_1+N_2-|\beta |}$. Hence symbol orders add under multiplication:

We now address the remaining unheaded requests. To define the oscillatory integral as a functional on the space of test functions, choose $\eta \in C_c^{\infty }({\mathbb{R}}^k)$ equal to one near zero and setThe limit is not an assertion of absolute convergence of the original frequency integral. Split off a bounded-frequency part, which is smooth in $x$. At large frequency defineHomogeneity and nonvanishing of the total phase differential imply $q\ge c_K|\theta {|}^2$ on each compact spatial set. Also $L{e}^{i\Phi }=e^{i\Phi }$. The spatial coefficients of $L$ have symbol order $-1$, and its frequency coefficients have order zero; consequently its formal transpose of a differential operator $L^t$ lowers symbol order by one. After $r>N+k$ integrations by parts, the high-frequency pairing has an absolutely integrable amplitude $(L^t{)}^r[(1-\chi (\theta ))a(x,\theta )f(x)]$, where $\chi$ is a fixed cutoff near zero. Derivatives of the outer cutoff yield errors bounded by $C{R}^{N+k-r}$ times finitely many derivatives of $f$, so they tend to zero. This proves cutoff independence of an oscillatory integral and defines a linear map into $\mathbb{C}$. Its distributional continuity is permitted as an assumption, and is also visible from this finite-derivative estimate.

The singular support of a distribution is the complement of the largest open subset on which it equals a smooth regular distribution. The stationary-direction bound for singular support needs conic support of an oscillatory amplitude. WriteIts radial lift is the closed conic enlargement of the ordinary support. The generally valid bound isIf the amplitude support is conic, this is exactly the printed bound. It is also the standard interpretation when support in the frequency directions is understood conically.

To prove the bound, take a point outside its right-hand side. Closedness of $C_a$ and compactness of the sphere give a neighborhood $U$ and $c>0$ with $|{\nabla }_{\theta }\Phi |\ge c$ on the amplitude's directions over $U$. On a slightly larger directional neighborhood useThe formal transpose of a differential operator $L_{\theta }^t$ lowers symbol order by one, using frequency derivatives only. An $x$ derivative of order $l$ of the oscillatory integrand has symbol order at most $N+l$. Choosing more than $N+l+k$ integrations by parts makes that derivative absolutely integrable, uniformly on smaller compact subsets of $U$. This works for every $l$, so $I_{\Phi }(a)$ is smooth on $U$.

The ordinary amplitude support can miss a singular-support limit phenomenon requires a qualification here: with unrestricted ordinary support the printed inclusion is false. An explicit counterexample uses $X={\mathbb{R}}^2$, $x=(s,t)$, one frequency variable and $\Phi (s,t,\theta )=s\theta$. Choose $t_j=1/j$, $w_j=1/(10j^2)$ and $\chi \in C_c^{\infty }((-1/4,1/4))$ with $\chi (0)=1$. Its translates ${\chi }_j(t)=\chi ((t-t_j)/w_j)$ have disjoint supports. Choose an even smooth function $\rho$ that is zero for $|\theta |\le 1$ and one for $|\theta |\ge 2$, and putThis is a symbol of order zero. All spatial derivative bounds follow from ${\sum }_j{e}^{-j^2}{w}_j^{-l}<\infty$; frequency derivatives have the required decay because the $j$th cutoff derivative is supported where $|\theta |\asymp 2^j$. Near every point $(s,0,\theta )$ with finite $\theta$, all large-$j$ terms vanish through the frequency cutoff and all remaining terms vanish in a small $t$ neighborhood. Thus no point $(0,0,\theta )$ lies in its ordinary support.

Nevertheless its oscillatory integral iswhere $G$ is smooth. Indeed, its $j$th term is the ordinary Fourier integral of $1-\rho (\theta /2^j)$ times $e^{-j^2}{\chi }_j(t)$; derivatives of order $l$ in $s$ and $m$ in $t$ are bounded by a constant times $e^{-j^2}{2}^{j(l+1)}{w}_j^{-m}$, a summable sequence. Since $B(t_j)\ne 0$, each $(0,t_j)$ is singular. Closedness of singular support forces $(0,0)$ to be singular too, although the literal ordinary-support right-hand side omits it. The closed conic support includes this limiting point and resolves the defect.

Finally use Fourier inversion distributionally. The constant amplitude gives $(2\pi {)}^{-n}\int e^{ix\cdot \theta }d\theta ={\delta }_0$, and differentiation of the exponential supplies a factor $i\theta$. Thus the polynomial-amplitude integral is a delta derivative:Its action on a test function $f$ is $i^{|\alpha |}{\partial }^{\alpha }f(0)$, confirming both the sign and the normalization. This is the delta derivatives from polynomial oscillatory amplitudes identity.

With $⟨\theta ⟩=(1+|\theta {|}^2{)}^{1/2}$, the symbol class $\mathrm{Sym}(X,{\mathbb{R}}^k;N)=S_{1,0}^N$ consists of complex-valued $a\in C^{\infty }(X\times {\mathbb{R}}^k)$ such that, for every compact $K⋐X$ and every pair of multi-indices $\alpha ,\beta$,An $x$ derivative preserves symbol order, while a $\theta$ derivative lowers it by one. The definition is local in $x$ and imposes estimates of every derivative order.

In the homogeneous convention, a phase function is real, smooth on $X\times ({\mathbb{R}}^k\setminus \{0\})$, positively homogeneous of degree one in $\theta$, and satisfiesOnly the large-frequency behavior is relevant to the regularization argument; any prescribed smooth modification near $\theta =0$ contributes an ordinary convergent integral. If the phase is presented on all of $X\times {\mathbb{R}}^k$, homogeneity can be required outside a bounded frequency region. Equivalently for this construction, its large-frequency symbol bounds and the following uniform estimate are sufficient:For a homogeneous phase function, this follows from compactness of $K\times S^{k-1}$. Nonvanishing of the total gradient is the condition; ${\nabla }_{\theta }\Phi$ may vanish. The displayed definition does not impose the stronger rank condition sometimes called a nondegenerate phase in Fourier-integral-operator theory.

We construct the oscillatory integral by integration by parts in both $x$ and $\theta$. For $|\theta |\ge 1$, setThen $L(e^{i\Phi })=e^{i\Phi }$. The coefficients of the $x$ derivatives have symbol order $-1$ and those of the $\theta$ derivatives have order zero. Use the bilinear formal transpose of a differential operator, not the Hermitian adjoint:The symbol order reduction by a phase integration operator follows from the product rule and the symbol class estimates: $L^t:S_{1,0}^s\to S_{1,0}^{s-1}$ locally in $x$. Differentiating the coefficients preserves precisely these orders because $q$ is elliptic of degree two on each compact $x$ set.

Let $\zeta \in C_c^{\infty }({\mathbb{R}}^k)$ equal one on $|\theta |\le 1$ and vanish on $|\theta |\ge 2$. The low-frequency part is an ordinary integral. For an integer $r>N+k$, define the high-frequency part byIf the homogeneous phase is only defined off zero, its value at the single point zero is immaterial to the low-frequency integral; its modulus remains one. In the high-frequency term all coefficients are used away from zero. For $\mathrm{supp}\phi \subset K$, repeated product rules giveSince $N-r<-k$, both terms converge absolutely. They are linear in the test function, with the continuity boundThus $I_{\Phi }(a)\in {\mathcal{D}}^{\prime }(X)$, with order of a distribution at most $r$ on each compact set.

To check that this construction is the intended oscillatory integral and does not depend on $r$, $\zeta$ or the large-frequency cutoff, let $\chi \in C_c^{\infty }({\mathbb{R}}^k)$ equal one near zero and formAfter $r$ applications of integration by parts, the term without a derivative on $\chi$ tends to the absolutely convergent high-frequency term by the dominated convergence theorem. Each extra term is supported in an annulus $|\theta |\asymp R$. A derivative of $\chi (\theta /R)$ contributes $R^{-1}$, so those terms have total absolute value at most $C_K{R}^{N-r+k}{\max }_{|\alpha |\le r}\Vert {\partial }^{\alpha }\phi {\Vert }_{\infty }\to 0$. The low-frequency part is unchanged for large $R$. Hence $J_R(\phi )$ tends to the displayed functional, independently of all cutoffs and of the permissible number of integrations. This proves the cutoff independence of an oscillatory integral as well as its continuity.

For the final distribution, the derivative convention $⟨{\delta }^{\prime },g⟩=-g^{\prime }(0)$ identifies it as $-x_2{\delta }^{\prime }(x_1)$. The delta derivatives from polynomial oscillatory amplitudes identity yieldsHere $k=1$, $a\in \mathrm{Sym}({\mathbb{R}}^2,\mathbb{R};1)$, and ${\nabla }_x\Phi =(\theta ,0)$ is nonzero for $\theta \ne 0$. Thus the phase is admissible even though ${\partial }_{\theta }\Phi =x_1$ vanishes on the distribution's support.

For a direct proof of the sign and normalization, put $g(s)={\int }_{\mathbb{R}}{x}_2\phi (s,x_2)d{x}_2\in C_c^{\infty }(\mathbb{R})$. The cutoff integral isThe Fourier transform $\hat{g}$ is rapidly decreasing, so the dominated convergence theorem and differentiated Fourier inversion giveThis is exactly the required pairing, so the oscillatory representation gives the prescribed distribution, not its negative or a multiple of it.

For a homogeneous phase function, let $q=|{\nabla }_x\Phi {|}^2+|\theta {|}^2|{\nabla }_{\theta }\Phi {|}^2$ and $L=(iq{)}^{-1}({\nabla }_x\Phi \cdot {\nabla }_x+|\theta {|}^2{\nabla }_{\theta }\Phi \cdot {\nabla }_{\theta })$. Then $L{e}^{i\Phi }=e^{i\Phi }$. The formal transpose of a differential operator $L^t$ lowers the symbol class order by one: the $x$ coefficients have order $-1$, while the order-zero frequency coefficients accompany a frequency derivative.