AN-04 proof edition · CC0; exact source credit and separately licensed prerequisites
A fundamental solution without ellipticity
Let P(D) be any nonzero constant-coefficient differential operator on Rn, with D=−i∂. We construct a distribution E satisfying P(D)E=δ0, and then the fixed smooth local inverse needed by the hyperbolicity transport recursion. Neither ellipticity nor a continuous choice as P 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 be the finite-dimensional complex vector space of polynomials in n variables of degree at most m, equipped with
∥Q∥Pm2=∣α∣≤m∑∣∂αQ(0)∣2.(FS1)
This is a norm because the finite Taylor formula gives every coefficient. On a fixed compact set of ζ's, evaluation is a bounded linear functional of these coefficients, uniformly in ζ.
For each nonzero Q, there is a point a∈Cn such that the polynomial 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 a with Q(a)=0. Hence q(1)=0.
The polynomial q has finitely many nonzero roots. Choose a radius r>0, small enough that ∣ra∣<1/4, distinct from their absolute values. If a=0, any sufficiently small radius is suitable and the orbit below is the single point zero. The compact orbit
{eiθra:0≤θ≤2π}(FS2)
contains no zero of Q. Its minimum ∣Q∣ is positive. A sufficiently small neighborhood of this orbit is still contained in {∣ζ∣<1/2} and has a positive lower bound for ∣Q∣. Take a nonnegative smooth bump ψ of integral one supported in a small ball about ra, whose entire rotation orbit is in that neighborhood, and define, for Lebesgue measure dA on Cn=R2n,
ϕQ(ζ)=2π1∫02πψ(eiθζ)dθ.(FS3)
It is smooth, nonnegative, rotation invariant, compactly supported in the ball of radius 1/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ψ, where ∣Q∣ has the positive lower bound.
F2. One smoothly varying averaging kernel
Normalize polynomials by FS1. For each point Qj of the unit sphere the preceding bump works for every polynomial in a sufficiently small neighborhood Uj of Qj, with a uniform lower bound cj>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 Uj. 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 form a smooth nonnegative partition with supports contained in Uj. Set
Φ(Q,ζ)=j∑λj(Q/∥Q∥Pm)ϕQj(ζ),Q=0.(FS4)
It is smooth jointly in Q=0,ζ, nonnegative, rotation invariant in ζ, of integral one, and all its ζ-supports lie in a fixed compact set Z⊂{∣ζ∣<1/2}. Finiteness gives constants c,C>0 with
0≤Φ≤C,∣Q(ζ)∣≥c∥Q∥Pmon suppζΦ(Q,⋅).(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(ζ), defined as zero near zeros of Q, is smooth. Locally in Q=0, 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 φ, define for complex w
φ(w)=∫Rne−ix⋅wφ(x)dx.(FS6)
The integral is over a compact real set. For fixed real ξ and complex ζ, uniform convergence of the exponential series on that set shows
2π1∫02πφ(−ξ−eiθζ)dθ=φ(−ξ).(FS7)
In detail, expand exp(ieiθx⋅ζ). Its degree j term contains eijθ, whose angular average is zero for j≥1 by direct integration, and is one for j=0. Uniform convergence permits integration term by term. The factor eix⋅ξφ(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,ζ) is our rotation-invariant kernel, average the change of variables ζ↦eiθζ in its integral against FS6. Fubini is legitimate on these compact sets, and FS7 gives
∫CnΦ(Q,ζ)φ(−ξ−ζ)dA(ζ)=φ(−ξ).(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 P of degree at most m, put Pξ(ζ)=P(ξ+ζ). If d is its actual degree and ∂αP is a nonzero derivative of degree d, that derivative is a fixed nonzero constant. Its term in FS1 gives
∥Pξ∥Pm≥cP>0for every real ξ.(FS9)
The statement also holds for a nonzero constant P. Define the bilinear distribution candidate
E(φ)=(2π)−n∫Rn∫CnP(ξ+ζ)Φ(Pξ,ζ)φ(−ξ−ζ)dA(ζ)dξ.(FS10)
We first justify it. If suppφ⊂A, with A compact, write FS6 at −ξ−ζ as the ordinary Fourier integral with real phase x⋅ξ and amplitude eix⋅ζφ(x). Repeated integration by parts with 1−Δx, whose action on eix⋅ξ is multiplication by 1+∣ξ∣2, gives
∣φ(−ξ−ζ)∣≤CA,N,Z(1+∣ξ∣2)−N∣α∣≤2NmaxAsup∣∂αφ∣,ζ∈Z.(FS11)
The exponential and all the finitely many derivatives are bounded uniformly on A×Z. The support has finite volume, so the integral remainder is bounded as written. Choose 2N>n; the ξ weight is integrable, as follows by summing dyadic annuli of volume at most C2jn. FS5 and FS9 uniformly bound the quotient in FS10. The ζ-support has finite volume. Thus FS10 is absolutely convergent, linear and bounded by a finite smooth-test seminorm on each compact A. It defines a distribution of finite order on every such compact.
The transpose of P(D) in the bilinear convention is P(−D), without coefficient conjugation. Integration by parts in the compact test integral proves, also for complex frequencies,
P(−D)φ(−ξ−ζ)=P(ξ+ζ)φ(−ξ−ζ).(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 Φ. Equations FS8 and the earlier full Fourier inversion theorem now give
(P(D)E)(φ)=E(P(−D)φ)=(2π)−n∫Rnφ(−ξ)dξ=φ(0).(FS13)
Hence P(D)E=δ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⋐U be relatively compact patches, and choose a fixed compact smooth χ supported in U, equal to one on a neighborhood of V. For g∈C∞(U), extend χg by zero and set
Tg(x)=Ey((χg)(x−y)).(FS14)
For x in any compact output set A, the test in y is supported in the fixed compact difference set A−suppχ. On a slightly enlarged such set, E has some order k. Difference quotients of the smooth compact input, with all derivatives through k, 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
A×Lsup∣∂xα∂λβTg∣≤C∣γ∣≤k∑suppχ×Lsup∣∂xα+γ∂λβ(χg)∣.(FS15)
Here L 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 g, with constants fixed by χ,A,L.
Constant-coefficient differentiation of FS14, and the definition of distributional derivatives, give P(D)Tg=χg, hence P(D)Tg=g on V. The output is smooth on all of Rn, so a finite transport recursion can always form its next right side on the original larger patch U, apply the same cutoff and inverse, and retain the equation on the same V. There is no accumulated loss of domain. For a smooth nonvanishing coefficient a, applying this inverse to g/a solves aP(D)w=g; no commutation with a 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ˉ, with its precise sign, remains the previously proved function-normal-form provider.