AN-04 proof edition · CC0; exact source credit and separately licensed prerequisites

A fundamental solution without ellipticity

Let P(D)P(D) be any nonzero constant-coefficient differential operator on Rn\mathbb R^n, with D=−i∂D=-i\partial. We construct a distribution EE satisfying P(D)E=δ0P(D)E=\delta_0, and then the fixed smooth local inverse needed by the hyperbolicity transport recursion. Neither ellipticity nor a continuous choice as PP varies is assumed.

Independent programme exposition and proof: GPT-6 Astra (OpenAI), Ultra, 5 October 2026; CC0. The mathematical source is Hörmander I, approved 2003 eBook, ISBN 978-3-642-61497-2, Theorems 7.3.10–7.3.12. Hörmander II, approved 2005 eBook, ISBN 978-3-540-26964-9, Theorem 10.2.1 gives a stronger regularity treatment. We need only the fundamental solution and smooth convolution bounds proved here. Exact earlier polynomial, Fourier, measure and finite-calculus proofs are identified in the proof map.

The exact map connects the complete polynomial proof, exponential series and circle calculus, Fourier inversion and measure and distribution foundations. The convex-support companion and exact figure supply the other U032 prerequisite and illustrate the rotating set used here.

F1. Avoiding polynomial zeros on a rotating compact set

Let Pm\mathcal P_m be the finite-dimensional complex vector space of polynomials in nn variables of degree at most mm, equipped with

∥Q∥Pm2=∑∣α∣≤m∣∂αQ(0)∣2.(FS1) \|Q\|_{\mathcal P_m}^2=\sum_{|\alpha|\le m}|\partial^\alpha Q(0)|^2. \tag{FS1}

This is a norm because the finite Taylor formula gives every coefficient. On a fixed compact set of ζ\zeta's, evaluation is a bounded linear functional of these coefficients, uniformly in ζ\zeta.

For each nonzero QQ, there is a point a∈Cna\in\mathbb C^n such that the polynomial q(z)=Q(za)q(z)=Q(za) is not identically zero. Here is a direct algebraic verification. A nonzero polynomial in one variable has at most its degree many zeros by successive division. In several variables, regard it as a polynomial in the last variable; by induction choose the earlier variables so one of its nonzero coefficient polynomials is nonzero. The resulting one-variable polynomial is nonzero, and a value outside its finite zero set gives aa with Q(a)≠0Q(a)\ne0. Hence q(1)≠0q(1)\ne0.

The polynomial qq has finitely many nonzero roots. Choose a radius r>0r>0, small enough that ∣ra∣<1/4|ra|<1/4, distinct from their absolute values. If a=0a=0, any sufficiently small radius is suitable and the orbit below is the single point zero. The compact orbit

{eiθra:0≤θ≤2π}(FS2) \{e^{i\theta}ra:0\le\theta\le2\pi\} \tag{FS2}

contains no zero of QQ. Its minimum ∣Q∣|Q| is positive. A sufficiently small neighborhood of this orbit is still contained in {∣ζ∣<1/2}\{|\zeta|<1/2\} and has a positive lower bound for ∣Q∣|Q|. Take a nonnegative smooth bump ψ\psi of integral one supported in a small ball about rara, whose entire rotation orbit is in that neighborhood, and define, for Lebesgue measure dAdA on Cn=R2n\mathbb C^n=\mathbb R^{2n},

ϕQ(ζ)=12π∫02πψ(eiθζ) dθ.(FS3) \phi_Q(\zeta)=\frac1{2\pi}\int_0^{2\pi}\psi(e^{i\theta}\zeta)\,d\theta. \tag{FS3}

It is smooth, nonnegative, rotation invariant, compactly supported in the ball of radius 1/21/2, and has integral one. Differentiation under the compact integral proves smoothness; a shift of the periodic angular variable proves invariance; the real rotation has Jacobian one, so the affine substitution theorem and Fubini prove the integral assertion. Its support consists only of rotations of points in supp⁡ψ\operatorname{supp}\psi, where ∣Q∣|Q| has the positive lower bound.

F2. One smoothly varying averaging kernel

Normalize polynomials by FS1. For each point QjQ_j of the unit sphere the preceding bump works for every polynomial in a sufficiently small neighborhood UjU_j of QjQ_j, with a uniform lower bound cj>0c_j>0 on its support. Indeed, evaluation is uniformly continuous in the coefficients on that compact support. Compactness of the coefficient sphere gives finitely many such neighborhoods.

We spell out the finite partition used here. Shrink finitely many coefficient balls so they still cover the sphere and their closed larger balls lie in the respective UjU_j. Smooth coefficient bumps equal to one on the smaller balls and zero outside the larger balls have a positive sum on the sphere. Divide each by that sum. Their restrictions λj\lambda_j form a smooth nonnegative partition with supports contained in UjU_j. Set

Φ(Q,ζ)=∑jλj(Q/∥Q∥Pm)ϕQj(ζ),Q≠0.(FS4) \Phi(Q,\zeta)=\sum_j\lambda_j(Q/\|Q\|_{\mathcal P_m})\phi_{Q_j}(\zeta), \qquad Q\ne0. \tag{FS4}

It is smooth jointly in Q≠0,ζQ\ne0,\zeta, nonnegative, rotation invariant in ζ\zeta, of integral one, and all its ζ\zeta-supports lie in a fixed compact set Z⊂{∣ζ∣<1/2}Z\subset\{|\zeta|<1/2\}. Finiteness gives constants c,C>0c,C>0 with

0≤Φ≤C,∣Q(ζ)∣≥c∥Q∥Pmon supp⁡ζΦ(Q,⋅).(FS5) 0\le\Phi\le C,\qquad |Q(\zeta)|\ge c\|Q\|_{\mathcal P_m} \quad\hbox{on }\operatorname{supp}_{\zeta}\Phi(Q,\cdot). \tag{FS5}

At a support point at least one relevant bump has that lower bound; taking the minimum of the finitely many constants proves FS5, including limits at support boundaries. The quotient Φ(Q,ζ)/Q(ζ)\Phi(Q,\zeta)/Q(\zeta), defined as zero near zeros of QQ, is smooth. Locally in Q≠0Q\ne0, FS5 separates the support by a positive margin from those zeros, which justifies this extension.

F3. The needed averaging identity, proved from exponentials

For a compactly supported smooth test φ\varphi, define for complex ww

φ^(w)=∫Rne−ix⋅wφ(x) dx.(FS6) \widehat\varphi(w)=\int_{\mathbb R^n}e^{-ix\cdot w}\varphi(x)\,dx. \tag{FS6}

The integral is over a compact real set. For fixed real ξ\xi and complex ζ\zeta, uniform convergence of the exponential series on that set shows

12π∫02πφ^(−ξ−eiθζ) dθ=φ^(−ξ).(FS7) \frac1{2\pi}\int_0^{2\pi} \widehat\varphi(-\xi-e^{i\theta}\zeta)\,d\theta =\widehat\varphi(-\xi). \tag{FS7}

In detail, expand exp⁡(ieiθx⋅ζ)\exp(i e^{i\theta}x\cdot\zeta). Its degree jj term contains eijθe^{ij\theta}, whose angular average is zero for j≥1j\ge1 by direct integration, and is one for j=0j=0. Uniform convergence permits integration term by term. The factor eix⋅ξφ(x)e^{ix\cdot\xi}\varphi(x) is integrable, and uniform convergence on the compact support permits the remaining exchange. The elementary exponential and circle integrals are proved in U001.

If Φ(Q,ζ)\Phi(Q,\zeta) is our rotation-invariant kernel, average the change of variables ζ↦eiθζ\zeta\mapsto e^{i\theta}\zeta in its integral against FS6. Fubini is legitimate on these compact sets, and FS7 gives

∫CnΦ(Q,ζ)φ^(−ξ−ζ) dA(ζ)=φ^(−ξ).(FS8) \int_{\mathbb C^n}\Phi(Q,\zeta) \widehat\varphi(-\xi-\zeta)\,dA(\zeta) =\widehat\varphi(-\xi). \tag{FS8}

Thus the precise mean-value assertion used below has its full proof. No general theorem about holomorphic extension or division is being substituted for it.

F4. Construction of the distribution and verification of its equation

For a fixed nonzero polynomial PP of degree at most mm, put Pξ(ζ)=P(ξ+ζ)P_\xi(\zeta)=P(\xi+\zeta). If dd is its actual degree and ∂αP\partial^\alpha P is a nonzero derivative of degree dd, that derivative is a fixed nonzero constant. Its term in FS1 gives

∥Pξ∥Pm≥cP>0for every real ξ.(FS9) \|P_\xi\|_{\mathcal P_m}\ge c_P>0\quad\hbox{for every real }\xi. \tag{FS9}

The statement also holds for a nonzero constant PP. Define the bilinear distribution candidate

E(φ)=(2π)−n∫Rn∫CnΦ(Pξ,ζ)P(ξ+ζ)φ^(−ξ−ζ) dA(ζ) dξ.(FS10) E(\varphi)=(2\pi)^{-n}\int_{\mathbb R^n}\int_{\mathbb C^n} \frac{\Phi(P_\xi,\zeta)}{P(\xi+\zeta)} \widehat\varphi(-\xi-\zeta)\,dA(\zeta)\,d\xi. \tag{FS10}

We first justify it. If supp⁡φ⊂A\operatorname{supp}\varphi\subset A, with AA compact, write FS6 at −ξ−ζ-\xi-\zeta as the ordinary Fourier integral with real phase x⋅ξx\cdot\xi and amplitude eix⋅ζφ(x)e^{ix\cdot\zeta}\varphi(x). Repeated integration by parts with 1−Δx1-\Delta_x, whose action on eix⋅ξe^{ix\cdot\xi} is multiplication by 1+∣ξ∣21+|\xi|^2, gives

∣φ^(−ξ−ζ)∣≤CA,N,Z(1+∣ξ∣2)−Nmax⁡∣α∣≤2Nsup⁡A∣∂αφ∣,ζ∈Z.(FS11) |\widehat\varphi(-\xi-\zeta)| \le C_{A,N,Z}(1+|\xi|^2)^{-N} \max_{|\alpha|\le2N}\sup_A|\partial^\alpha\varphi|, \qquad \zeta\in Z. \tag{FS11}

The exponential and all the finitely many derivatives are bounded uniformly on A×ZA\times Z. The support has finite volume, so the integral remainder is bounded as written. Choose 2N>n2N>n; the ξ\xi weight is integrable, as follows by summing dyadic annuli of volume at most C2jnC2^{jn}. FS5 and FS9 uniformly bound the quotient in FS10. The ζ\zeta-support has finite volume. Thus FS10 is absolutely convergent, linear and bounded by a finite smooth-test seminorm on each compact AA. It defines a distribution of finite order on every such compact.

The transpose of P(D)P(D) in the bilinear convention is P(−D)P(-D), without coefficient conjugation. Integration by parts in the compact test integral proves, also for complex frequencies,

P(−D)φ^(−ξ−ζ)=P(ξ+ζ)φ^(−ξ−ζ).(FS12) \widehat{P(-D)\varphi}(-\xi-\zeta) =P(\xi+\zeta)\widehat\varphi(-\xi-\zeta). \tag{FS12}

Using the absolutely convergent defining formula with this test, cancel the denominator. The resulting integral is still absolutely convergent by FS11 and boundedness of Φ\Phi. Equations FS8 and the earlier full Fourier inversion theorem now give

(P(D)E)(φ)=E(P(−D)φ)=(2π)−n∫Rnφ^(−ξ) dξ=φ(0).(FS13) (P(D)E)(\varphi)=E(P(-D)\varphi) =(2\pi)^{-n}\int_{\mathbb R^n}\widehat\varphi(-\xi)\,d\xi =\varphi(0). \tag{FS13}

Hence P(D)E=δ0P(D)E=\delta_0. This proves existence for every nonzero constant-coefficient operator. It asserts neither temperedness of this construction nor causal support. If there are zero variables, a nonzero constant operator is inverted simply by multiplication by its reciprocal.

F5. A fixed local right inverse, with every parameter derivative

Let V⋐UV\Subset U be relatively compact patches, and choose a fixed compact smooth χ\chi supported in UU, equal to one on a neighborhood of V‾\overline V. For g∈C∞(U)g\in C^\infty(U), extend χg\chi g by zero and set

Tg(x)=Ey((χg)(x−y)).(FS14) Tg(x)=E_y\bigl((\chi g)(x-y)\bigr). \tag{FS14}

For xx in any compact output set AA, the test in yy is supported in the fixed compact difference set A−supp⁡χA-\operatorname{supp}\chi. On a slightly enlarged such set, EE has some order kk. Difference quotients of the smooth compact input, with all derivatives through kk, converge uniformly by the integral Taylor formula. Passing them through the finite-order distribution bound proves differentiation under the pairing. Repeating it for output derivatives and passive parameters gives joint smoothness and

sup⁡A×L∣∂xα∂λβTg∣≤C∑∣γ∣≤ksup⁡supp⁡χ×L∣∂xα+γ∂λβ(χg)∣.(FS15) \sup_{A\times L}|\partial_x^\alpha\partial_\lambda^\beta Tg| \le C\sum_{|\gamma|\le k} \sup_{\operatorname{supp}\chi\times L} |\partial_x^{\alpha+\gamma}\partial_\lambda^\beta(\chi g)|. \tag{FS15}

Here LL is any fixed compact parameter set in its open parameter domain. The same Taylor argument on a slightly larger parameter compact proves continuity of these derivatives. The finite product rule replaces the right side by finitely many seminorms of gg, with constants fixed by χ,A,L\chi,A,L.

Constant-coefficient differentiation of FS14, and the definition of distributional derivatives, give P(D)Tg=χgP(D)Tg=\chi g, hence P(D)Tg=gP(D)Tg=g on VV. The output is smooth on all of Rn\mathbb R^n, so a finite transport recursion can always form its next right side on the original larger patch UU, apply the same cutoff and inverse, and retain the equation on the same VV. There is no accumulated loss of domain. For a smooth nonvanishing coefficient aa, applying this inverse to g/ag/a solves aP(D)w=gaP(D)w=g; no commutation with aa is needed.

These statements supply the general local inverse in U032, including its nonelliptic cone operator. The separate planar inverse for Dt−iDs=−2∂zˉD_t-iD_s=-2\partial_{\bar z}, with its precise sign, remains the previously proved function-normal-form provider.