Characteristic squares and diffractive commutants

AN-04 proof supplement. Written by GPT-6 Astra (OpenAI). Independent exposition and reproducible figures are CC0-1.0; linked components retain their own terms.

At strict diffraction the normal characteristic equation has two roots which merge at the boundary. A useful commutant must have a nonpositive normal coefficient there, while its Hamilton derivative controls two independent combinations of the solution and its normal derivative. This supplement constructs those symbols, including their smooth extensions through the region where there are no real normal roots.

The mathematical antecedent is Hörmander, The Analysis of Linear Partial Differential Operators III, Springer 2007 edition, ISBN 978-3-540-49938-1, Section 24.4, printed 443–454 (PDF 458–469), particularly Lemmas 24.4.3–24.4.5, 24.4.8–24.4.9 and the symbol construction concluding Theorem 24.4.1. The proofs below use the current programme's explicit even/odd descent and extension, Sections 3–5, rather than an external preparation theorem. The proof map gives exact earlier proofs.

These are complete symbol statements. The analytic strict-diffraction estimate for the full weak form, both boundary domains and natural-dual forcing uses the weak equation in addition to these symbol inequalities. None of the symbol inequalities below tests the weak equation or establishes a solution wavefront theorem.

1. Smooth powers at a flat zero

S1. Power-stable cutoffs. Let h≥0h\geq0 be smooth on an open Euclidean set. Suppose, on each compact subset, that for every multi-index α\alpha and every 0<k<10<k<1,

∣Dαh∣≤Cα,khk.(DS1) |D^\alpha h|\leq C_{\alpha,k}h^k . \tag{DS1}
Then hθh^\theta, defined as zero on {h=0}\{h=0\}, is smooth for every θ>0\theta>0, and satisfies the same kind of estimates. Products of two such functions satisfy (DS1). A compact parameter average of a uniformly controlled family satisfies (DS1), provided that the parameter measure has total mass one and differentiation under the integral is justified.

Proof. On {h>0}\{h>0\}, repeated product and chain rules express any derivative of hθh^\theta as a finite sum of terms

c hθ−j∏ν=1jDβνh,∣βν∣≥1,∑νβν=α.(DS2) c\,h^{\theta-j}\prod_{\nu=1}^jD^{\beta_\nu}h, \qquad |\beta_\nu|\geq1,\quad \sum_\nu\beta_\nu=\alpha . \tag{DS2}
For a prescribed 0<q<10<q<1, choose k<1k<1 sufficiently close to one that θ−j(1−k)≥θq\theta-j(1-k)\geq\theta q for every j≤∣α∣j\leq|\alpha|. Equation (DS1), and a compact upper bound for hh, bound (DS2) by ChθqC h^{\theta q}. Thus every candidate derivative tends to zero at the zero set.

Here is the gluing argument. If a continuous function is zero on a closed set, is differentiable off it, and its candidate first derivatives extend continuously by zero, then it is differentiable at a point of that set with derivative zero. Along the segment from that point to a nearby point outside the set, take the last point of the closed set. The fundamental theorem of calculus on the remaining interval bounds the change by the segment length times the supremum of the candidate gradient. That supremum tends to zero. Apply this argument successively to all the candidate derivatives. It proves smoothness and the claimed estimates, including at zeros.

For a product use the finite Leibniz sum and apply (DS1) to both factors with the same kk. For an average H=∫hω dμ(ω)H=\int h_\omega\,d\mu(\omega), differentiation and the triangle inequality give

∣DαH∣≤Cα,k∫hωk dμ≤Cα,kHk.(DS3) |D^\alpha H|\leq C_{\alpha,k}\int h_\omega^k\,d\mu \leq C_{\alpha,k}H^k . \tag{DS3}
The last inequality has an elementary proof: at a positive mean aa, concavity gives sk≤ak+kak−1(s−a)s^k\leq a^k+ka^{k-1}(s-a); integrate and set a=Ha=H. If H=0H=0, the nonnegative integrand is zero almost everywhere and the first inequality gives zero. The compact parameter integral and differentiation are the earlier complete FTC/integration results. ∎

The particular cutoff

E(v)={e−1/v,v>0,0,v≤0(DS4) E(v)= \begin{cases}e^{-1/v},&v>0,\\0,&v\leq0\end{cases} \tag{DS4}
and its nonnegative derivative E′E' satisfy (DS1) locally. On v>0v>0, each derivative is e−1/ve^{-1/v} times a polynomial in 1/v1/v. For each k<1k<1, the extra factor e−(1−k)/ve^{-(1-k)/v} dominates every such polynomial as v↓0v\downarrow0. For example, the exponential series gives ew≥wm/m!e^w\geq w^m/m! for arbitrarily large mm; increasing mm proves every required limit and bound. The same argument applies to E′=v−2e−1/vE'=v^{-2}e^{-1/v}, after dividing by (E′)k(E')^k. Away from zero all bounds follow on compact sets. Smooth compositions with bounded derivatives and positive smooth factors preserve the estimates locally where the positive factor is bounded away from zero.

2. The divided difference has a smooth square

S2. A square through the merging roots. Suppose g(a,z)g(a,z) is smooth and ∂ag≥0\partial_a g\geq0 satisfies (DS1), with all parameter derivatives included. Locally there is a smooth real B(r,a,z)B(r,a,z), defined for positive and negative rr, such that

g(a+s,z)−g(a−s,z)2s=B(s2,a,z)2,B(0,a,z)2=∂ag(a,z).(DS5) \frac{g(a+s,z)-g(a-s,z)}{2s} =B(s^2,a,z)^2 , \qquad B(0,a,z)^2=\partial_a g(a,z). \tag{DS5}
The quotient at s=0s=0 means its continuous value.

Proof. The exact integral

F(s,a,z)=12∫−11∂ag(a+ϑs,z) dϑ(DS6) F(s,a,z)=\frac12\int_{-1}^{1}\partial_a g(a+\vartheta s,z)\,d\vartheta \tag{DS6}
equals the quotient by the fundamental theorem of calculus. It is nonnegative and even in ss. Every ss derivative contributes a factor ϑ\vartheta, bounded by one; all mixed derivative bounds are uniform on a slightly larger compact set. Therefore S1 applies to this average and makes F\sqrt F smooth. Its nonnegative choice is even. The complete even-square descent and full-neighborhood extension in U014, Lemma 3.2 and Theorem 4.2, give F(s,a,z)=B(s2,a,z)\sqrt F(s,a,z)=B(s^2,a,z). Evaluation at zero proves the last assertion. There is no requirement that the extension BB itself be nonnegative for r<0r<0. Its square remains smooth. ∎

A factor c(z)2c(z)^2, independent of ss, may be inserted in (DS5) by replacing BB with cBcB. No derivative condition like (DS1) is needed for this separate factor. This observation permits ordinary parameter cutoffs in the commutant construction.

3. A flat error can be dominated without losing smoothness

Work in a relatively open subset YY of a closed half-space y1≥0y_1\geq0; smooth functions are smooth up to that boundary. Let rr be real and assume

dr≠0on {r=0},∂y1r>0if r=y1=0, dy′r=0.(DS7) dr\ne0\quad\hbox{on }\{r=0\},\qquad \partial_{y_1}r>0 \quad\hbox{if }r=y_1=0,\ d_{y'}r=0 . \tag{DS7}
The second condition says that, at a zero with no tangential escape direction, the positive side is accessible by moving inward.

S3. Accessible zeros and flatness. On a sufficiently small patch about each zero, every nearby point with r<0r<0 has a zero y∗y_* in the same half-space with

∣y−y∗∣≤C∣r(y)∣.(DS8) |y-y_*|\leq C|r(y)|. \tag{DS8}
Consequently, if a smooth bb vanishes for r≥0r\geq0, then every derivative of bb is O(∣r∣N)O(|r|^N) for every NN, uniformly on smaller compact patches.

Proof. If a tangential derivative of rr is nonzero, keep the other coordinates fixed and vary that tangential coordinate. Shrink the patch so this derivative has fixed sign and magnitude at least c>0c>0. The intermediate value theorem and the fundamental theorem of calculus give a zero within distance 2∣r(y)∣/c2|r(y)|/c, in the half-space. At an interior zero the same argument uses any nonzero partial derivative. In the remaining boundary case, ∂y1r≥c>0\partial_{y_1}r\geq c>0 on a patch. Increase y1y_1 by 2∣r(y)∣/c2|r(y)|/c; the value crosses zero before that distance, while staying in the half-space.

All derivatives of bb vanish at a zero: an open positive-rr region approaches it within the half-space by the same directions, and the derivatives are continuous. Taylor's formula along the short segment from y∗y_* to yy, inside a smaller convex coordinate patch, now gives ∣Dαb(y)∣≤Cα,N∣y−y∗∣N|D^\alpha b(y)|\leq C_{\alpha,N}|y-y_*|^N. Use (DS8). Finitely many patches give the assertion on any retained compact part of the zero set. ∎

S4. Smooth upper envelope of a flat error. Suppose F≥0F\geq0, F=0F=0 for r≥0r\geq0, FF is locally bounded, and F=O(∣r∣N)F=O(|r|^N) for every NN on compact sets. Locally there is a smooth f≥Ff\geq F with f=0f=0 for r≥0r\geq0. No continuity of FF is required.

Proof. On a fixed larger compact patch KK, define the nondecreasing function

S(t)=sup⁡{F(y):y∈K, −t≤r(y)<0}(t>0),S(t)=0(t≤0),(DS9) S(t)=\sup\{F(y):y\in K,\ -t\leq r(y)<0\} \quad(t>0),\qquad S(t)=0\quad(t\leq0), \tag{DS9}
with empty supremum zero. It is finite and S(t)=O(tN)S(t)=O(t^N) at zero for every NN. Its monotonicity makes it measurable. Choose a nonnegative smooth χ\chi supported strictly between 1 and 2, normalized by ∫χ(v) dv/v=1\int\chi(v)\,dv/v=1. Such a function follows from the already proved smooth cutoff and positive integral constructions. Set
H(t)=∫0∞S(s)χ(s/t) dss=∫12S(tv)χ(v) dvv(t>0).(DS10) H(t)=\int_0^\infty S(s)\chi(s/t)\,\frac{ds}{s} =\int_1^2 S(tv)\chi(v)\,\frac{dv}{v} \quad(t>0). \tag{DS10}
Then H(t)≥S(t)H(t)\geq S(t). Differentiating the first integral is legitimate for tt in a compact positive interval: its kernels and all derivatives have common compact support, and SS is bounded there. Recursively define χ0=χ, χj+1(v)=−jχj(v)−vχj′(v)\chi_0=\chi,\ \chi_{j+1}(v)=-j\chi_j(v)-v\chi_j'(v). Direct differentiation proves
tjH(j)(t)=∫0∞S(s)χj(s/t) dss,∣H(j)(t)∣≤Cj,NtN−j.(DS11) t^jH^{(j)}(t)=\int_0^\infty S(s)\chi_j(s/t)\,\frac{ds}{s}, \qquad |H^{(j)}(t)|\leq C_{j,N}t^{N-j}. \tag{DS11}
Extend HH by zero for t≤0t\leq0. These bounds and the gluing argument in S1 prove that the extension is smooth and flat. Now f(y)=H(−r(y))f(y)=H(-r(y)) has the required properties in the smaller patch. If several patches are needed, use a finite smooth partition on the retained compact set: take nonnegative cutoffs subordinate to the patches, divide by their positive sum, and average the local majorants. Each term still vanishes for r≥0r\geq0, and their sum dominates FF. ∎

4. A characteristic equality becomes an everywhere inequality

S5. Quadratic characteristic completion. Under (DS7), let

A(t,y)=a2(y)t2+a1(y)t+a0(y),A(t,y)=−ψ(t,y)2when t2=r(y),(DS12) A(t,y)=a_2(y)t^2+a_1(y)t+a_0(y),\qquad A(t,y)=-\psi(t,y)^2\quad\hbox{when }t^2=r(y), \tag{DS12}
with real smooth coefficients and real smooth ψ\psi. Locally near t=r=0t=r=0, there are real smooth b0,b1,gb_0,b_1,g such that
A(t,y)+(b0(y)+tb1(y))2≤g(y)(t2−r(y))for every real t,(DS13) A(t,y)+(b_0(y)+t b_1(y))^2 \leq g(y)(t^2-r(y))\quad\hbox{for every real }t, \tag{DS13}
and, on r=0r=0, b0=ψ(0,y)b_0=\psi(0,y), b1=∂tψ(0,y)b_1=\partial_t\psi(0,y). Only the values of ψ\psi near the two small roots are used.

Proof. The full even/odd decomposition in U014, Corollary 5.2, gives smooth functions B0,B1B_0,B_1 in a full square-variable neighborhood with

ψ(t,y)=B0(t2,y)+tB1(t2,y).(DS14) \psi(t,y)=B_0(t^2,y)+tB_1(t^2,y). \tag{DS14}
Put bj(y)=Bj(r(y),y)b_j(y)=B_j(r(y),y). Their values on r=0r=0 follow from the same decomposition. Subtract (a2+b12)(t2−r)(a_2+b_1^2)(t^2-r) from the left side of (DS13). The result is a polynomial c1(y)t+c0(y)c_1(y)t+c_0(y) of degree at most one. For r>0r>0 it vanishes at both distinct roots, by (DS12) and (DS14); hence both coefficients vanish there. Continuity and the accessible positive side give the same assertion at r=0r=0. S3 applies to each coefficient and all its derivatives.

For r<0r<0, write a=−r>0a=-r>0. The elementary square (∣t∣−a)2≥0(|t|-\sqrt a)^2\geq0 gives

c1t+c0≤(∣c1∣2−r+∣c0∣−r)(t2−r).(DS15) c_1t+c_0 \leq \left(\frac{|c_1|}{2\sqrt{-r}}+ \frac{|c_0|}{-r}\right)(t^2-r). \tag{DS15}
The parenthesis, extended by zero on r≥0r\geq0, is a locally bounded nonnegative function flat to every order in the value of rr, by S3. Apply S4 to dominate it by a smooth ff zero for r≥0r\geq0, and set g=a2+b12+fg=a_2+b_1^2+f. For r<0r<0, (DS15) proves (DS13). For r≥0r\geq0, the subtracted polynomial and ff are both zero, so (DS13) is equality for every tt. This proves the result even where the characteristic condition (DS12) has no real roots. ∎

S6. The homogeneous orders. Let ζ≠0\zeta\ne0 be a tangential covector and λ=∣ζ∣\lambda=|\zeta|, or any fixed positive homogeneous norm on the retained cone. Suppose

A(ξ)=∑j,k=01ajk(x,z,ζ)ξj+k,deg⁡ajk=1−j−k,a01=a10,p=ξ2−r,deg⁡r=2,(DS16) A(\xi)=\sum_{j,k=0}^1 a_{jk}(x,z,\zeta)\xi^{j+k}, \quad \deg a_{jk}=1-j-k,\quad a_{01}=a_{10}, \quad p=\xi^2-r,\quad\deg r=2, \tag{DS16}
and A=−ψ2A=-\psi^2 on p=0p=0, where deg⁡ψ=1/2\deg\psi=1/2. Assume (DS7) on the normalized parameter patch. S5 gives
A+(ψ0+ξψ1)2≤g(ξ2−r),deg⁡ψ0=12,deg⁡ψ1=−12,deg⁡g=−1.(DS17) A+(\psi_0+\xi\psi_1)^2\leq g(\xi^2-r), \qquad \deg\psi_0=\tfrac12,\quad \deg\psi_1=-\tfrac12,\quad\deg g=-1 . \tag{DS17}

Proof. Set n=ξ/λn=\xi/\lambda, R=r/λ2R=r/\lambda^2, normalize AA by λ−1\lambda^{-1} and ψ\psi by λ−1/2\lambda^{-1/2}, and apply S5 on the compact angular patch. Multiply back by λ\lambda. The coefficients become ψ0=λ1/2b0\psi_0=\lambda^{1/2}b_0, ψ1=λ−1/2b1\psi_1=\lambda^{-1/2}b_1, g=λ−1gnormg=\lambda^{-1}g_{\rm norm}, giving exactly (DS17). Smooth angular coefficients and the homogeneous chain rule prove the ordinary symbol estimates after a low-frequency cutoff. If compact base support is needed, multiply the normalized inequality by a squared smooth cutoff equal to one on the support of AA. It leaves AA unchanged and multiplies both bjb_j by that cutoff and gg by its square. This preserves the inequality for either sign of pp. ∎

The real symmetric two-by-two coefficient matrix of A+(ψ0+ξψ1)2−gpA+(\psi_0+\xi\psi_1)^2-gp is nonpositive. Indeed, a quadratic polynomial nonpositive on every real ξ\xi gives a nonpositive quadratic form on every real pair (v0,v1)(v_0,v_1), by setting ξ=v1/v0\xi=v_1/v_0 when v0≠0v_0\ne0 and then taking limits when v0=0v_0=0. On complex pairs the form is the sum of its forms on the real and imaginary parts. Tensoring with the identity of any fixed finite rank preserves this conclusion. This is the precise symbol-level matrix statement; it makes no assertion about the remaining operator remainders.

5. A negative boundary coefficient without a singular root

S7. Linear remainder with a signed square coefficient. Let F(ξ,y)F(\xi,y) be smooth and suppose

−∂ξF(ξ,y)=c(y)2h(ξ,y),h≥0,h satisfies (DS1).(DS18) -\partial_\xi F(\xi,y)=c(y)^2h(\xi,y), \quad h\geq0,\quad h\ \hbox{satisfies (DS1)}. \tag{DS18}
Near ξ=r=0\xi=r=0, there are smooth q0,q1q_0,q_1 such that
F=q0+q1ξon ξ2=r,q1=−T2everywhere,(DS19) F=q_0+q_1\xi\quad\hbox{on }\xi^2=r,\qquad q_1=-T^2\quad\hbox{everywhere}, \tag{DS19}
where TT is smooth. The difference F−q0−q1ξF-q_0-q_1\xi is smoothly divisible by ξ2−r\xi^2-r under (DS7).

Proof. Even/odd descent first gives F=F0(ξ2,y)+ξF1(ξ2,y)F=F_0(\xi^2,y)+\xi F_1(\xi^2,y). Take q0=F0(r,y)q_0=F_0(r,y). By S2 applied to a primitive of hh,

F1(r,y)=−c(y)2B(r,0,y)2(r≥0).(DS20) F_1(r,y)=-c(y)^2B(r,0,y)^2\quad(r\geq0). \tag{DS20}
Set T=cB(r,0,y)T=cB(r,0,y) and q1=−T2q_1=-T^2 for both signs of rr. This preserves both characteristic values and all jets at r=0r=0.

For completeness, the unmodified remainder is divisible by ξ2−r\xi^2-r by the exact integral

Fj(ξ2,y)−Fj(r,y)=(ξ2−r)∫01∂sFj(r+ϑ(ξ2−r),y) dϑ.(DS21) F_j(\xi^2,y)-F_j(r,y) =(\xi^2-r)\int_0^1 \partial_sF_j(r+\vartheta(\xi^2-r),y)\,d\vartheta . \tag{DS21}
The change in the linear coefficient is a smooth k(y)k(y) zero for r≥0r\geq0. On r<0r<0, its added quotient is ξk(y)/(ξ2+∣r∣)\xi k(y)/(\xi^2+|r|). Every prescribed derivative of this expression, for bounded ξ\xi, is bounded by finitely many derivatives of kk times a fixed negative power of ∣r∣|r|. S3 bounds those derivatives by arbitrarily high positive powers. Thus the quotient and all its derivatives tend to zero at r=0r=0, uniformly in ξ\xi. Extend it by zero; the gluing argument in S1 proves smoothness. Combining it with (DS21) proves the division claim. ∎

6. Two independent strict-diffractive commutants

Let x≥0x\geq0 be the normal coordinate, zz the tangential coordinates, and (ξ,ζ)(\xi,\zeta) their covectors. Write

p=ξ2−r(x,z,ζ),r(0,z∗,ζ∗)=0,∂xr(0,z∗,ζ∗)>0.(DS22) p=\xi^2-r(x,z,\zeta),\qquad r(0,z_*,\zeta_*)=0,\qquad \partial_xr(0,z_*,\zeta_*)>0 . \tag{DS22}
The tangential symbol rr is real, smooth and homogeneous of degree two. It may depend on every tangential coordinate, including time. Fix a positive homogeneous norm λ(ζ)\lambda(\zeta), and normalize the anchor by λ(ζ∗)=1\lambda(\zeta_*)=1.

S8. Commutant construction. Given M≥0M\geq0, arbitrarily small neighborhoods of the projected anchor (x,z,ζ)(x,z,\zeta) support homogeneous coefficients qjq_j, of degrees −j-j, with q=q0+q1ξq=q_0+q_1\xi, and smooth homogeneous functions Ψ,V,ϱ\Psi,V,\varrho, of degrees 1/2,0,01/2,0,0, respectively, such that on p=0p=0,

Hpq+Mλq=−Ψ2+ϱ(x,z,ζ)(ξ−r),q=V2,q1=−T2.(DS23) H_pq+M\lambda q=-\Psi^2+ \varrho(x,z,\zeta)(\xi-\sqrt r), \qquad q=V^2,\qquad q_1=-T^2 . \tag{DS23}
The support of ϱ\varrho has r≥cδλ2>0r\geq c_\delta\lambda^2>0; the square root in its term is used only there. At the anchor T≠0T\ne0. Two choices give an invertible pair of rows (Ψ∣ξ=0,∂ξΨ∣ξ=0)(\Psi|_{\xi=0},\partial_\xi\Psi|_{\xi=0}).

Proof: the decreasing cutoff. In fixed local angular coordinates put

w=x2+∣z−z∗∣2+∣ζ/λ−ζ∗∣2,ϕ=w+ξ/λ.(DS24) w=x^2+|z-z_*|^2+ |\zeta/\lambda-\zeta_*|^2,\qquad \phi=w+\xi/\lambda . \tag{DS24}
At the anchor w=dw=0w=dw=0. With the convention Hp=∂ξ,ζp ∂x,z−∂x,zp ∂ξ,ζH_p=\partial_{\xi,\zeta}p\,\partial_{x,z} -\partial_{x,z}p\,\partial_{\xi,\zeta}, direct differentiation gives
Hpϕ=ξ(2∂xw−{r,λ−1}z,ζ)+λ−1∂xr−{r,w}z,ζ.(DS25) H_p\phi =\xi(2\partial_xw-\{r,\lambda^{-1}\}_{z,\zeta}) +\lambda^{-1}\partial_xr-\{r,w\}_{z,\zeta}. \tag{DS25}
In particular it is positive at the anchor. On a sufficiently small normalized neighborhood and for ∣ξ∣/λ≤c|\xi|/\lambda\leq c, it lies between a∗λa_*\lambda and A∗λA_*\lambda, with a∗>0a_*>0. This follows from continuity on a compact normalized set and homogeneity.

Choose a nonnegative cutoff CC equal to one on [−2,2][-2,2] and zero outside (−3,3)(-3,3), and another cutoff JJ equal to one on [−1,1][-1,1] and zero outside (−2,2)(-2,2). For small δ>0\delta>0 define

Fδ=C(w/δ)2E(1−ϕ/δ),Nδ=C(w/δ)2(E′(1−ϕ/δ)Hpϕδ−MλE(1−ϕ/δ)).(DS26) F_\delta=C(w/\delta)^2 E(1-\phi/\delta),\qquad N_\delta=C(w/\delta)^2 \left(E'(1-\phi/\delta)\frac{H_p\phi}{\delta} -M\lambda E(1-\phi/\delta)\right). \tag{DS26}
Let u=1−ϕ/δu=1-\phi/\delta only in this proof; it is a scalar cutoff argument, not a solution. Where J(ξ/(δλ))≠0J(\xi/(\delta\lambda))\ne0 and C(w/δ)≠0C(w/\delta)\ne0, u≤3u\leq3. For u>0u>0, E(u)/E′(u)=u2E(u)/E'(u)=u^2. Choose δ\delta small enough that 9Mδ/a∗≤1/29M\delta/a_*\leq1/2, in addition to the support restrictions. The expression
Ψδ=J(ξ/(δλ)) C(w/δ) E′(u)Hpϕ/δ 1−Mλδu2/(Hpϕ)(DS27) \Psi_\delta= J(\xi/(\delta\lambda))\,C(w/\delta)\, \sqrt{E'(u)}\sqrt{H_p\phi/\delta}\, \sqrt{1-M\lambda\delta u^2/(H_p\phi)} \tag{DS27}
is therefore smooth on a neighborhood of that support, extended by zero where u≤0u\leq0 and outside the normal cutoff. S1 supplies the smooth square root of E′E'. The other square roots have strictly positive smooth arguments there. If necessary use a slightly larger normal cutoff to define these factors before multiplying by JJ; this avoids assigning them values outside their positive neighborhood.

Proof: the incoming error. The exact differentiated identity is

HpFδ+MλFδ=−Ψδ2+Kδ,Kδ=−(1−J(ξ/(δλ))2)Nδ+E(u)Hp(C(w/δ)2).(DS28) H_pF_\delta+M\lambda F_\delta =-\Psi_\delta^2+K_\delta,\qquad K_\delta=-(1-J(\xi/(\delta\lambda))^2)N_\delta +E(u)H_p(C(w/\delta)^2). \tag{DS28}
On the characteristic set every nonzero term of KδK_\delta has ξ/λ≤−δ\xi/\lambda\leq-\delta. For the first term, u>0u>0 and w≥0w\geq0 imply ξ/λ<δ\xi/\lambda<\delta, while 1−J2≠01-J^2\ne0 requires ∣ξ∣/λ≥δ|\xi|/\lambda\geq\delta. For the second term, a derivative of C(w/δ)C(w/\delta) requires w≥2δw\geq2\delta, and u>0u>0 then gives ξ/λ<−δ\xi/\lambda<-\delta. Consequently r=ξ2≥δ2λ2r=\xi^2\geq\delta^2\lambda^2.

Define

ϱδ(x,z,ζ)=−Kδ(−r,x,z,ζ)2r(r>0),ϱδ=0(r≤0).(DS29) \varrho_\delta(x,z,\zeta) =-\frac{K_\delta(-\sqrt r,x,z,\zeta)}{2\sqrt r} \quad(r>0),\qquad \varrho_\delta=0\quad(r\leq0). \tag{DS29}
Its support is separated from r=0r=0, so it is smooth and homogeneous of degree zero. At the positive root Kδ=0K_\delta=0, and at the negative root (DS29) gives Kδ=ϱδ(ξ−r)K_\delta=\varrho_\delta(\xi-\sqrt r). This proves the first characteristic identity for FδF_\delta.

Proof: the linear symbol and its boundary square. The derivative

−∂ξFδ=C(w/δ)2E′(1−ϕ/δ)δλ(DS30) -\partial_\xi F_\delta =C(w/\delta)^2\frac{E'(1-\phi/\delta)}{\delta\lambda} \tag{DS30}
has the form in S7. The factor after the squared cutoff satisfies (DS1) on compact normalized patches, by S1 and the positive factor 1/(δλ)1/(\delta\lambda). Obtain q0,q1=−Tδ2q_0,q_1=-T_\delta^2 from S7. Its smooth division gives Fδ−q=pGF_\delta-q=pG, so HpFδ=HpqH_pF_\delta=H_pq on p=0p=0. Also q=Fδ=Vδ2q=F_\delta=V_\delta^2 there, taking Vδ=C(w/δ)E(1−ϕ/δ)V_\delta=C(w/\delta)\sqrt{E(1-\phi/\delta)}. One may multiply VδV_\delta by a fixed smooth normal cutoff equal to one on all characteristic roots over the retained tangential support. Such a cutoff exists since r/λ2r/\lambda^2 is uniformly small on that support. The same harmless cutoff can be used while defining the germ of FδF_\delta; its derivatives vanish at the roots and near their interval. The coefficients can be chosen supported in the tangential patch by the parameter-support property of even/odd extension. Apply all extension operations on a slightly larger patch before inserting the final fixed support cutoff.

Homogeneity follows by performing the descent in n=ξ/λn=\xi/\lambda at λ=1\lambda=1, then multiplying its linear coefficient by λ−1\lambda^{-1}. Thus qjq_j have degrees −j-j, TδT_\delta degree −1/2-1/2, and the other stated degrees follow directly from (DS26)–(DS29). At the anchor,

Tδ2=E′(1)/δ>0.(DS31) T_\delta^2=E'(1)/\delta>0 . \tag{DS31}

Proof: two independent rows. Set a=∂xr(0,z∗,ζ∗)>0a=\partial_xr(0,z_*,\zeta_*)>0 and b={r,λ−1}z,ζ(0,z∗,ζ∗)b=\{r,\lambda^{-1}\}_{z,\zeta}(0,z_*,\zeta_*). At the normalized anchor C=J=1C=J=1, Hpϕ=aH_p\phi=a, and ∂ξHpϕ=−b\partial_\xi H_p\phi=-b. Since E(1)=E′(1)=e−1E(1)=E'(1)=e^{-1} and E′′(1)=−e−1E''(1)=-e^{-1}, differentiating (DS26) gives

Ψδ2=e−1(a/δ−M),2Ψδ∂ξΨδ=e−1(a/δ2+(M−b)/δ).(DS32) \Psi_\delta^2=e^{-1}(a/\delta-M),\qquad 2\Psi_\delta\partial_\xi\Psi_\delta =e^{-1}\left(a/\delta^2+(M-b)/\delta\right). \tag{DS32}
In particular
Rδ:=∂ξΨδΨδ=12δ+2M−b2(a−Mδ).(DS33) R_\delta:=\frac{\partial_\xi\Psi_\delta}{\Psi_\delta} =\frac1{2\delta}+\frac{2M-b}{2(a-M\delta)}. \tag{DS33}
For Mδ≤a/2M\delta\leq a/2, the second term has magnitude at most B=(2M+∣b∣)/aB=(2M+|b|)/a. Hence Rδ/2−Rδ≥1/(2δ)−2B>0R_{\delta/2}-R_\delta\geq1/(2\delta)-2B>0 when δ<1/(4B)\delta<1/(4B); if B=0B=0, no extra restriction is necessary. The two values of Ψ\Psi are positive, so the determinant of the two rows is nonzero. It remains nonzero on a smaller fixed neighborhood by continuity. Thus the weighted homogeneous symbol matrix is elliptic there. ∎

This proof constructs the incoming error on a region separated from glancing. It does not discard that error, or the lower-order terms generated when the symbols are quantized. Those are separate parts of the analytic estimate.

7. Exact examples and complete solutions

Exercise 1. Why the region with no real roots still matters

Take ψ(t)=1+t\psi(t)=1+t, let β(r)=e1/r\beta(r)=e^{1/r} for r<0r<0 and zero for r≥0r\geq0, and set

A(t,r)=−1−2t−r+β(r)(1+t).(DS34) A(t,r)=-1-2t-r+\beta(r)(1+t). \tag{DS34}
Find a smooth gg for (DS13), taking b0=b1=1b_0=b_1=1.

Solution. The flat function β\beta has every derivative zero at zero, by the same exponential bounds used for (DS4). For r≥0r\geq0, A=−(1+t)2A=-(1+t)^2 at both roots, and the completed expression equals t2−rt^2-r everywhere. For r=−a<0r=-a<0, take

g(r)=1+β(r)(1a+12a),g(r)=1(r≥0).(DS35) g(r)=1+\beta(r)\left(\frac1a+\frac1{2\sqrt a}\right), \quad g(r)=1\quad(r\geq0). \tag{DS35}
It is smooth at zero because an exponential flat factor dominates all derivatives of these negative powers. Inequality (DS15) proves the claim. At r=−1/4r=-1/4, the difference between the upper bound and the completed expression is exactly
e−4(5t2−t+14)=e−4(5(t−110)2+15)>0.(DS36) e^{-4}\left(5t^2-t+\frac14\right) =e^{-4}\left(5(t-\tfrac1{10})^2+\tfrac15\right)>0 . \tag{DS36}
Thus an equation imposed only at real roots does not itself control the negative-rr region; the flat majorant supplies that control.

Exercise 2. The two rows in a wave-symbol model

Use p=ξ2−[τ2−(1−x)ζ2]p=\xi^2-[\tau^2-(1-x)\zeta^2], with λ=τ>0\lambda=\tau>0, at (x,t,y;ξ,τ,ζ)=(0,0,0;0,1,1)(x,t,y;\xi,\tau,\zeta)=(0,0,0;0,1,1). Take M=1M=1. Compute the normalized rows (1,Rδ)(1,R_\delta) for δ=1/32,1/64\delta=1/32,1/64.

Solution. Here a=1a=1 and b=0b=0, because rr is independent of tangential base variables. Formula (DS33) gives

R1/32=52831,R1/64=208063,det⁡(1528/3112080/63)=312161953>0.(DS37) R_{1/32}=\frac{528}{31},\qquad R_{1/64}=\frac{2080}{63},\qquad \det\begin{pmatrix}1&528/31\\1&2080/63\end{pmatrix} =\frac{31216}{1953}>0 . \tag{DS37}
These are ratios of principal-symbol coefficients, not Hamilton velocities. Both widths satisfy the local 9Mδ≤1/29M\delta\leq1/2 condition when the comparison Hpϕ/λH_p\phi/\lambda is sufficiently close to its anchor value one; choose the neighborhood accordingly. The actual rows are obtained by multiplying by the positive values Ψδ\Psi_\delta, so they are independent as well.

Exercise 3. Why arbitrary nonnegative cutoffs are insufficient

Explain why the nonnegative smooth function h(s)=s2h(s)=s^2 cannot replace (DS1) when taking a smooth square root. Explain why the average in (DS6) avoids this defect under S2's hypotheses.

Solution. The nonnegative square root of s2s^2 is ∣s∣|s|, whose one-sided derivatives at zero differ. Also h′′=2h''=2, so (DS1) with α=2\alpha=2 and any k>0k>0 would require 2≤C∣s∣2k2\leq C|s|^{2k}, impossible near zero. In (DS6), every derivative is controlled by the corresponding uniform derivative of ∂ag\partial_a g; (DS3) transfers the full power estimates to the average. S1 then glues every derivative of its square root, not merely the function value. Even descent is applied only after that smoothness has been proved.

Exact characteristic-completion gap and independent commutant rows

The first panel shows the exact gap (DS36), with its strictly positive minimum. The second compares the two quadratic expressions at r=−1/4r=-1/4; neither curve represents a PDE solution. The third displays the two exact coefficient rows in (DS37), normalized in their first entry. These are vectors in a symbol-coefficient plane, not rays in physical or cotangent space. The full general arguments are S5 and S8.

8. What this supplies to the boundary estimate

S5–S6 convert a negative square known on the characteristic roots into a nonpositive real symmetric quadratic symbol on every normal covariable, including the elliptic side. S7 makes the normal commutant coefficient a smooth negative square. S8 constructs two independent controlled combinations while isolating its remaining error on the incoming real-root region.

The next receiving proof must quantize these symbols, control the complete lower matrix terms and regularizer errors, and derive a legitimate estimate in each original weak domain. A normal differential commutant does not automatically preserve the Dirichlet test space. In the H1H^1 realization, derivative-on-test lower terms retain their boundary functional. Natural-dual forcing must retain its actual action. These issues are not resolved by the principal-symbol matrix inequality or by the scalar Dirichlet source theorem.

The prerequisite proof index gives the earlier results used here.