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

Again use $D_j=-i{\partial }_j$ and distribution dual pairings with complex bilinearity. The Malgrange–Ehrenpreis theorem asserts that every nonzero constant-coefficient differential operator defined by a polynomial has a fundamental solution of a linear differential operator:If $P=a\ne 0$ is constant, take $E={\delta }_0/a$. For positive degree $m$, write $P_m$ for the highest homogeneous part. A nonzero polynomial cannot vanish on all of ${\mathbb{R}}^n$, so choose a real unit vector $e$ with $P_m(e)\ne 0$. A change of coordinates by an orthogonal matrix makes it the last coordinate direction. In these coordinates,Crucially the leading coefficient $a$ is constant and never vanishes as ${\xi }^{\prime }\in {\mathbb{R}}^{n-1}$ varies.

The finite-height polynomial root avoidance construction gives the needed Hörmander staircase. Set $h_{\ell }=3\ell$ for $0\le \ell \le m$. For each fixed ${\xi }^{\prime }$, the fundamental theorem of algebra supplies $m$ roots $z_1,\dots ,z_m$, with multiplicity. A root's imaginary part can be within distance strictly less than one of at most one candidate height, since those heights are separated by three. The pigeonhole principle therefore leaves at least one height with $|h_{\ell }-\mathrm{Im}{z}_j|\ge 1$ for all roots. At that height, for every real $s$,No continuous labelling of the roots is required. To select the heights measurably, defineEach $B_{\ell }$ is closed, since the coefficients are continuous in ${\xi }^{\prime }$. Continuity in $s$ makes the rational test equivalent to the bound for all $s\in \mathbb{R}$. The $A_{\ell }$ are therefore disjoint Borel sets covering ${\mathbb{R}}^{n-1}$. The union of horizontal fibresis a Hörmander staircase, with bounded imaginary height and a uniform nonzero denominator. For $n=1$, the transverse space has one point, and this is simply the choice of one horizontal line.

For example, $P({\xi }_1,z)={\xi }_1^2+z^2$ has roots $z=\pm i{\xi }_1$. One may select $h=0$ for $|{\xi }_1|\ge 1$, and $h=3$ for $|{\xi }_1|<1$. Then the first region has $|P|\ge 1$, while in the second each root is at least two units from the chosen height. The figure shows the height projection of these horizontal fibres, rather than additional connecting contours.

Define a linear functional on test functions byTo establish that this is a distribution, fix a compact set containing the support of $\phi$, inside a radius-$R$ ball. The Fourier decay in a bounded complex strip estimate follows by applying $(1-\Delta {)}^L$ to $e^{-h_{\ell }{x}_n}\phi (x)$:The finite bound on the heights absorbs all exponential factors and powers of $h_{\ell }$ into $C_{R,L}$. Taking $2L>n$ and using $|P|\ge |a|$ proves absolute convergence and the required finite-order continuity on each fixed support set. Hence $E\in {\mathcal{D}}^{\prime }$.

Now the formal transpose identity and the Fourier derivative identity giveThe polynomial denominator cancels in the pairing for $P(D)E$. For each fixed real ${\xi }^{\prime }$, the remaining function of $z$ is entire and rapidly decreasing in its real part throughout $0\le \mathrm{Im}z\le 3m$. The Cauchy integral theorem therefore allows each horizontal fibre to be shifted to the real axis; its truncated vertical end integrals tend to zero. Absolute convergence justifies Fubini's theorem and the finite sum over the measurable partition. ThusThe last equality is Fourier inversion. Transforming back from the orthogonal coordinates preserves ${\delta }_0$ and gives the required fundamental solution of a linear differential operator for the original operator. The staircase is allowed to jump: after cancellation the shift is performed on each one-dimensional fibre, so no assertion that $\Sigma$ is a smooth contour is used. This construction proves a distributional fundamental solution, without claiming a tempered-growth estimate.

For $f\in \mathcal{D}({\mathbb{R}}^n)$, use the smoothing convolution with a test functionIndeed $u_0(x)=⟨E_y,f(x-y)⟩$ is a smooth function, with ${\partial }^{\alpha }{u}_0(x)=⟨E_y,{\partial }^{\alpha }f(x-y)⟩$. On each compact set of $x$ values, all translated test functions have supports in one fixed compact set, so differentiation is justified by the continuity of the distribution. Compact support of $f$ makes the convolution well defined even when $E$ is not tempered. It does not imply compact support of $u_0$.

Finally the affine solution space of a linear equation isIf every homogeneous solution is represented by a smooth function, so is every $u=u_0+v$. Conversely, if every solution for this fixed $f$ is represented by a smooth function, then for any homogeneous solution $v$, the solution $u_0+v$ is smooth, and subtracting the smooth $u_0$ shows that $v$ is smooth. Hence All inhomogeneous solutions are smooth functions exactly when all homogeneous solutions are smooth functions. This is a global statement on ${\mathbb{R}}^n$ for the stated data; no additional local regularity theorem is being assumed.

The Malgrange–Ehrenpreis theorem states that every nonzero constant-coefficient linear partial differential operator $P(D)$ on ${\mathbb{R}}^n$ has a fundamental solution of a linear differential operator: there is an $E\in {\mathcal{D}}^{\prime }({\mathbb{R}}^n)$ such that $P(D)E={\delta }_0$.

Write $D=-i\partial$. After an orthogonal change of coordinates and multiplication by a nonzero constant, its polynomial symbol may be written as a monic polynomial in the last frequency,For each real ${\mu }^{\prime }$, this polynomial has $M$ complex roots counted with multiplicity. Among a fixed finite collection of horizontal lines at bounded heights, one can choose a line that stays a positive distance from all those roots. Continuity of the roots preserves the choice on a neighborhood $N({\mu }^{\prime })$. Take a countable locally finite cover by such neighborhoods, refine it to a measurable disjoint partition ${\mathbb{R}}^{n-1}={⨆}_j{\Delta }_j$, and let $c_j$ be the chosen height on ${\Delta }_j$. The resulting Hörmander staircasehas bounded heights and may be chosen so that $|P({\xi }^{\prime },s+i{c}_j)|\ge 1$ on each step.

For a test function $\phi$, defineThe Paley–Wiener–Schwartz theorem gives rapid decay in the real frequency directions and at most a fixed exponential factor in the bounded imaginary direction. Together with $|P|\ge 1$, this proves that the integral defines a continuous distribution. Applying $P(D)$ cancels the denominator. The remaining integrand is entire in $z$, so the Cauchy integral theorem shifts every horizontal contour to the real axis; the partition then recombines into ${\mathbb{R}}^{n-1}$. The Fourier inversion theorem giveswhich proves the theorem.

No. Non-entireness rules out a compactly supported solution, but the Malgrange–Ehrenpreis theorem gives a distributional fundamental solution of a linear differential operator $E$ with $P(D)E=\delta$. Since $v$ has compact support, the convolution $u=E\ast v$ is defined and satisfiesOne may choose a tempered fundamental solution for a constant-coefficient operator, so $u\in {\mathcal{S}}^{\prime }(\mathbb{R})$.

Use $D=-i\partial$, so $P(D)e^{i\lambda \cdot x}=P(\lambda )e^{i\lambda \cdot x}$. A concrete distributional division construction isThe locally integrable family, initially defined for sufficiently large $\mathrm{Re}s$, has a meromorphic continuation; $\mathrm{FP}$ denotes its finite part at zero. Multiplication before continuation givesThe right side is holomorphic at $s=0$ with value $1$, so comparison of constant Laurent coefficients yields $PT=1$. Thereforesatisfies $P(D)E={\delta }_0$. This finite-part formula explicitly realizes the Malgrange–Ehrenpreis theorem.