Wave smoothing uniform in the coefficients

A wave may be singular outside a forward cone while remaining smooth well inside it. We prove that the interior smoothness persists uniformly for a smooth family of operators approaching a constant spatial metric, including all coefficient-parameter derivatives and points on a Dirichlet boundary. The proof uses the full boundary observation estimate. Its observation operators act only where the solution is zero, including their residual kernels.

Independent exposition, proofs, exercises and figure: CC0-1.0. The antecedent is Lars Hörmander, The Analysis of Linear Partial Differential Operators III, Springer 2007 eBook, ISBN 978-3-540-49938-1, Lemmas 17.5.12–17.5.14 and the application in §24.7, printed pages 56–61 and 467–469. The proof map identifies the complete earlier programme arguments. Required human foundation proofs remain external.

1. The family and the conclusion

H0. A precise receiving statement. Write X=(r,y)∈R+nX=(r,y)\in\mathbb R^n_+, n≥1n\ge1, with r≥0r\ge0, and let aa range near a compact set K0K_0 in the face aq=0a_q=0 of a finite-dimensional parameter half-space. A smooth map c(a)∈∂R+nc(a)\in\partial\mathbb R^n_+ specifies the spatial center. Suppose

La=∂t2+Pa(X,DX),Pa=∑i,j=1nAaij(X)DiDj+∑i=1nbai(X)Di+ca0(X),D=−i∂,(WS1) L_a=\partial_t^2+P_a(X,D_X),\qquad P_a=\sum_{i,j=1}^n A_a^{ij}(X)D_iD_j+ \sum_{i=1}^n b_a^i(X)D_i+c_a^0(X),\quad D=-i\partial, \tag{WS1}

where AaA_a is real, symmetric and positive definite. All coefficients extend smoothly to an open parameter neighborhood of K0K_0, with uniform bounds on compact sets. The lower coefficients may be complex. At every a∈K0a\in K_0, the matrix Aa(X)A_a(X) is independent of XX. Let dad_a be the spatial length distance for the positive metric Aa−1A_a^{-1}, with paths in the closed half-space.

For aq>0a_q>0, let uau_a be jointly smooth in (a,t,X)(a,t,X), up to the spatial boundary, for ∣t∣<3|t|<3. Assume

Laua=0,ua(t,0,y)=0,supp⁡ua(0,⋅)∪supp⁡∂tua(0,⋅)⊂{X:da(X,c(a))≤2/5}.(WS2) L_au_a=0,\qquad u_a(t,0,y)=0,\qquad \operatorname{supp}u_a(0,\cdot)\cup \operatorname{supp}\partial_tu_a(0,\cdot) \subset\{X:d_a(X,c(a))\le2/5\}. \tag{WS2}

For every parameter multiindex α\alpha, assume the family ∂aαua(t,⋅)\partial_a^\alpha u_a(t,\cdot) is bounded in extendible spatial distributions, locally uniformly in t∈(−3,3)t\in(-3,3) and in aa near K0K_0. More explicitly, for each compact spatial test set and compact time interval there are a finite derivative seminorm pNαp_{N_\alpha} on tests smooth and flat at r=0r=0, and Cα<∞C_\alpha<\infty, such that

∣⟨∂aαua(t,⋅),ϕ⟩∣≤CαpNα(ϕ).(WS3) |\langle\partial_a^\alpha u_a(t,\cdot),\phi\rangle| \le C_\alpha p_{N_\alpha}(\phi). \tag{WS3}

A uniformly bounded set of distributions of a fixed finite order, as in the antecedent, has this property. The orders may increase with α\alpha.

Then there is one relative parameter neighborhood K′K' of K0K_0, independent of derivative order, such that for every α,j,β\alpha,j,\beta there is a finite constant with

sup⁡a∈K′, aq>0, 1≤t≤2da(X,c(a))≤2/5∣∂aα∂tj∂Xβua(t,X)∣≤Cαjβ.(WS4) \sup_{\substack{a\in K',\ a_q>0,\ 1\le t\le2\\ d_a(X,c(a))\le2/5}} |\partial_a^\alpha\partial_t^j\partial_X^\beta u_a(t,X)| \le C_{\alpha j\beta}. \tag{WS4}

The constants need not be bounded as the derivative order grows. No bound on the number of boundary reflections, no finite-contact hypothesis and no selfadjointness of the full operator is added. The proof below gives the assertion on a slightly larger spatial-time neighborhood of the displayed target. It does not assert smoothness on the entire surrounding backward cone.

2. Fixed geometric margins and finite propagation

F0. Reduce a reference spatial metric to the identity. Fix a0∈K0a^0\in K_0. Translate by c(a)c(a), retaining the parameter dependence of this translation. Choose a fixed linear map BB preserving the spatial half-space with

BTAa0−1B=I,B−1Aa0B−T=I.(WS5) B^T A_{a^0}^{-1}B=I, \qquad B^{-1}A_{a^0}B^{-T}=I . \tag{WS5}

To construct it, take an Aa0−1A_{a^0}^{-1}-orthonormal basis of the boundary tangent plane and its inward unit metric normal; use the latter as the first column. The columns form a basis, the first column has positive normal component, and the other columns have zero normal component. This proves both identities and the half-space property. The corresponding distance becomes Euclidean at a0a^0. Smooth change of variables preserves (WS3): transpose onto the test, retain its Jacobian, and use the finite chain rule. We henceforth use these coordinates and center zero.

On the fixed spatial ball ∣X∣≤10|X|\le10, shrink the parameter neighborhood so that

Aa≤c∗2I,c∗=101/100,(WS6) A_a\le c_*^2 I,\qquad c_*=101/100, \tag{WS6}

as quadratic forms, and retain a uniform positive lower bound. Each fixed coefficient derivative is bounded there. A path of metric length at most 2/52/5 cannot leave the unit ball: its portion reaching that sphere would have length at least 1/c∗>2/51/c_*>2/5. On paths remaining in this ball, metric length is at least Euclidean displacement divided by c∗c_*. Taking approximating paths in the distance infimum therefore gives

da(X,0)≤2/5⟹∣X∣≤(2/5)c∗=101/250.(WS7) d_a(X,0)\le2/5\quad\Longrightarrow\quad |X|\le(2/5)c_*=101/250. \tag{WS7}

For the closed endpoint, first use paths of length 2/5+ε<1/c∗2/5+\varepsilon<1/c_*, then let ε↓0\varepsilon\downarrow0. This argument also excludes a putative shortcut that travels outside the region of coefficient control.

F1. The common zero region comes from energy. In the working spatial ball ∣X∣<5|X|<5, finite propagation for (WS1) gives

ua(t,X)=0if∣X∣>c∗(2/5+∣t∣),∣t∣<3.(WS8) u_a(t,X)=0\quad\hbox{if}\quad |X|>c_*(2/5+|t|),\qquad |t|<3. \tag{WS8}

Here is the full domain-of-dependence use. Fix such a point and choose c′>c∗c'>c_* close enough that ∣X∣>2c∗/5+c′∣t∣|X|>2c_*/5+c'|t|. Its backward spatial cone to the plane t=0t=0, with radius c′c' times the remaining time, has a base disjoint from both initial supports in (WS2). Intersect the cone with the half-space. The conormal of its shrinking spherical side is c′ dt+d∣X−X∗∣c'\,dt+d|X-X_*|; its principal square is at least (c′)2−c∗2>0(c')^2-c_*^2>0, with future orientation. The artificial-face positivity DC:T008–T014 and the shrinking-ball proof DC:T032 show that weighted energy inside this cone is bounded by the initial energy, the forcing and the prescribed boundary value, all zero. In particular the solution vanishes at its tip, by taking nested cones with tips slightly beyond the chosen point. The physical boundary contributes the nonnegative zero-Dirichlet flux; no condition is imposed on an artificial timelike side. Corners and the radial center are handled by the monotone cutoff argument preceding DC:T014. The full complex lower terms are absorbed by its exponential weight.

All these cones lie in ∣X∣<9|X|<9 when ∣X∗∣<5|X_*|<5, ∣t∣<3|t|<3 and c′c' is close to c∗c_*; thus (WS6) suffices and no behavior at spatial infinity is assumed. For negative time use the past-directed inward multiplier of DC:T015, or repeat this argument after changing tt to −t-t. This proves (WS8). A bounded family of the already smooth solutions needs no limiting argument here.

3. Normalize before differentiating the parameter

N0. A time-preserving common normal form. Apply PBC:N0–G0 to LaL_a, using rr as the normal variable and (t,y)(t,y) as tangential variables. Since (WS1) has no mixed time-spatial principal coefficients, the flow that removes the normal cross terms leaves tt unchanged. Write its spatial part as

Φa(r,z)=(r,Ya(r,z)),Ya(0,z)=z,∂rYaj=(Aa1j/Aa11)(r,Ya),j>1.(WS9) \Phi_a(r,z)=(r,Y_a(r,z)),\quad Y_a(0,z)=z, \qquad \partial_rY_a^j= (A_a^{1j}/A_a^{11})(r,Y_a),\quad j>1. \tag{WS9}

The reference cross coefficients vanish. On the ball of radius ten they are uniformly small for aa close to a0a^0. The Picard construction and differentiated flow equations in PBC:N0 therefore work for 0≤r≤50\le r\le5, with a larger tangential buffer, not merely for an unspecified infinitesimal collar. For instance, a bound 1/5001/500 on the Euclidean norm of the right-hand side gives displacement at most 1/1001/100 over this interval. Further shrinking gives common inverse domains and

∣Φa(q)−q∣≤1/100(WS10) |\Phi_a(q)-q|\le1/100 \tag{WS10}

on the compact working boxes. In dimension n=1n=1 this flow is the identity. All fixed-order derivatives of the maps and their inverses are uniformly bounded.

Retain the full density factor WaW_a, the complex scalar normal gauge SaS_a, and the positive leading coefficient λa\lambda_a from PBC:D0–G0. With Fau=u(t,Φa(q))F_au=u(t,\Phi_a(q)), put

va=Sa−1WaFaua,Pa=Sa−1λa−1WaFaLaFa−1Wa−1Sa=Dr2−Ra(r,t,z,Dt,Dz).(WS11) v_a=S_a^{-1}W_a F_au_a, \qquad \mathcal P_a=S_a^{-1}\lambda_a^{-1} W_a F_a L_a F_a^{-1}W_a^{-1}S_a =D_r^2-R_a(r,t,z,D_t,D_z). \tag{WS11}

The operator RaR_a is a tangential differential operator of order at most two, with all lower terms retained and real quadratic principal symbol. The preceding coordinate, density and gauge formulas are exact full-operator formulas. Thus Pava=0\mathcal P_av_a=0, and va∣r=0=0v_a|_{r=0}=0. There is no hidden source in (WS11).

The transformed target in (WS4) has ∣q∣≤207/500|q|\le207/500 by (WS7)–(WS10). It lies in the interior, relative to the half-space, of the fixed compact set

M={99/100≤t≤201/100, ∣q∣≤21/50, r≥0}.(WS12) M=\{99/100\le t\le201/100,\ |q|\le21/50,\ r\ge0\}. \tag{WS12}

By (WS8) and (WS10), vav_a vanishes on the common open set

O={(t,q):∣q∣>∣t∣+1/2}(WS13) O=\{(t,q): |q|>|t|+1/2\} \tag{WS13}

within the working box. Indeed its possible support there satisfies

∣q∣≤c∗(2/5+∣t∣)+1/100≤∣t∣+111/250<∣t∣+1/2(∣t∣<3).(WS14) |q|\le c_*(2/5+|t|)+1/100 \le |t|+111/250<|t|+1/2\qquad(|t|<3). \tag{WS14}

The last gap is 7/1257/125. Multiplication by Sa−1WaS_a^{-1}W_a changes neither zero sets nor supports.

N1. Parameter derivatives retain an admissible weak bound. Set wα,a=∂aαvaw_{\alpha,a}=\partial_a^\alpha v_a. Differentiating (WS11) gives the exact triangular system

Pawα,a=fα,a,fα,a=∑0<β≤α(αβ)(∂aβRa)wα−β,a,f0,a=0.(WS15) \mathcal P_aw_{\alpha,a}=f_{\alpha,a},\qquad f_{\alpha,a}=\sum_{0<\beta\le\alpha} {\alpha\choose\beta}(\partial_a^\beta R_a) w_{\alpha-\beta,a},\qquad f_{0,a}=0. \tag{WS15}

Every ∂aβRa\partial_a^\beta R_a is tangential. In particular, it is a full boundary differential operator of order at most two. Ordinary Dr2D_r^2 does not appear in this sum. Differentiation also gives wα,a∣r=0=0w_{\alpha,a}|_{r=0}=0 and wα,a=0w_{\alpha,a}=0 on OO.

For completeness, each wα,aw_{\alpha,a} has a uniform local negative restricted Sobolev bound. Apply the finite parameter chain rule to Sa−1Waua(t,Φa(q))S_a^{-1}W_a u_a(t,\Phi_a(q)), including the preliminary translation. This gives finitely many coefficient multiples of spatial derivatives of ∂aγua\partial_a^\gamma u_a, with ∣γ∣≤∣α∣|\gamma|\le|\alpha|, pulled back by the same map. Transpose onto a flat test. The map preserves the boundary, its Jacobian and all relevant derivatives are uniformly bounded, and differentiation preserves flatness. Equation (WS3) therefore bounds the pairing by one finite seminorm on that test. Integrate over a compact time interval to obtain the same bound for spacetime tests. NGR:B1–B2, including its Fourier Sobolev inequality and Hahn–Banach restriction argument, now gives, for each fixed cutoff β0\beta_0,

∥β0wα,a∥H‾−Nα≤Cα.(WS16) \|\beta_0 w_{\alpha,a}\|_{\overline H^{-N_\alpha}} \le C_\alpha . \tag{WS16}

For the finitely many indices ∣α∣≤k|\alpha|\le k, take their largest order and enlarge their constants. No invariance of the original distribution family under coordinate changes, and no common order for infinitely many α\alpha, has been assumed.

4. Observation means actual input support

O0. An input-supported version of the chart estimate. In CRO:H0–H1 suppose the observation cone satisfies Γ2⊂Tb∗O∖0\Gamma_2\subset T_b^*O\setminus0, where OO is open in the base. The scalar observation test AoA_o can be chosen with its entire input projection compactly contained in OO. Consequently Aow=0A_o w=0 whenever w=0w=0 on OO.

To prove this refinement, make the following choices in that proof. For every nested test Ei,nE_{i,n} whose full microsupport is contained in Vi⋐Γ2V_i\Subset\Gamma_2, choose a right base cutoff χi∈Cc∞(O)\chi_i\in C_c^\infty(O) equal to one near the base projection of Vi‾\overline{V_i}, and a proper realization with complete input projection in its support. This retains ellipticity on UiU_i and the required full microsupport. Keep these exact tests as the observation list in CRO:T3. All comparisons with the local observation operators retain their actual full residuals and the separate source tests of CRO:R0–R1. Thus the defect-set argument CRO:D0–P3, including the parameter limit, is unchanged.

At its conclusion only finitely many Ei,nE_{i,n} occur. In the BSA:E0 comparison choose the common elliptic test QQ with a right cutoff in OO, equal to one near the union of their compact microsupports. The exact remainders Ei−BiQE_i-B_iQ have input projections in the union of the inputs of EiE_i and QQ. That union is compact in OO. In BSA:E1 choose the residual right cutoff χ\chi in OO, equal to one on this union. The scalar assembly is Ao=Q+USA_o=Q+US. Its relocation UU acts on the output alone, so it does not change these input supports. This proves the asserted support, with every residual included. It is stronger than a microsupport statement: a residual kernel could otherwise read the solution far outside OO.

O1. Full localization and tangential sources. Fix a compact cutoff β0\beta_0, equal to one on a neighborhood of an inner working box containing all rays and all actual input projections used below. Choose the coarse order in (WS16) first. Realize the chart estimates with their complete input projections in that inner box, as in PBC:T0. Coefficients can be smoothly extended to a larger half-space model while keeping the normal coefficient equal to one. For any resulting full source test TT,

TPa(β0wα,a)=Tfα,a,Ao(β0wα,a)=0.(WS17) T\mathcal P_a(\beta_0w_{\alpha,a})=T f_{\alpha,a}, \qquad A_o(\beta_0w_{\alpha,a})=0. \tag{WS17}

The first identity holds because TT has no input where a derivative of β0\beta_0 is nonzero. The second follows from O0. Targets with inputs in the same box similarly act on wα,aw_{\alpha,a} itself. The estimate is applied to smooth compact zero-Dirichlet inputs in the fixed Banach space H‾−N\overline H^{-N}. No density assertion for a larger graph domain is needed.

Moreover, for each summand of (WS15),

T(∂aβRa)∈Ψbord⁡T+2,WF⁡b′(T∂aβRa)⊂WF⁡b′T.(WS18) T(\partial_a^\beta R_a)\in\Psi_b^{\operatorname{ord}T+2}, \qquad \operatorname{WF}'_b\bigl(T\partial_a^\beta R_a\bigr) \subset\operatorname{WF}'_bT . \tag{WS18}

This is an actual full composition, including all residual terms. Tangential differentiation is a polynomial boundary symbol; the full composition and proper-support theorem GB:C6 applies. On the compact parameter set these operators form a bounded family in every fixed full-kernel seminorm. The differential factor does not enlarge the input projection. These facts are precisely what permits the source norm in (WS17) to use estimates already proved for lower parameter derivatives.

5. A compact family of backward rays

Which covectors can return from a forward cone?

G3. The exact flat reachability criterion. Use the flat symbol p=∣ξ∣2−τ2p=|\xi|^2-\tau^2, with spatial point q=(r,y)q=(r,y), r≥0r\ge0. Fix a real time shift aa. Let WaW_a be the characteristic covectors obtained by following every reflected ray backward from the open cone t0>a+∣q0∣t_0>a+|q_0|. Only nonzero characteristic covectors are under discussion: ∣ξ∣=∣τ∣≠0|\xi|=|\tau|\ne0. Then

(t,q,τ,ξ)∈Wa⟺Ha(t,q,τ,ξ):=τ⟨q,ξ⟩+(t−a)τ2>0.(WS33) (t,q,\tau,\xi)\in W_a \quad\Longleftrightarrow\quad H_a(t,q,\tau,\xi):= \tau\langle q,\xi\rangle+(t-a)\tau^2>0. \tag{WS33}
The assertion includes tangent rays and the compressed identification of the two normal roots at the boundary.

First work in the full spatial space. A future vertex (t0,q0)(t_0,q_0) lies on the characteristic line through the given covector precisely when

q0=q+(t−t0)ξ/τ,t0>t.(WS34) q_0=q+(t-t_0)\xi/\tau, \qquad t_0>t. \tag{WS34}
Put T=t−aT=t-a, T0=t0−aT_0=t_0-a, and b=τq+Tξb=\tau q+T\xi. The forward-cone condition is
∣q0∣<T0,T0>0⟺∣b−T0ξ∣2<T02τ2⟺∣b∣2<2T0⟨b,ξ⟩.(WS35) |q_0|<T_0, \quad T_0>0 \quad\Longleftrightarrow\quad |b-T_0\xi|^2<T_0^2\tau^2 \quad\Longleftrightarrow\quad |b|^2<2T_0\langle b,\xi\rangle . \tag{WS35}
Here the middle equivalence retains T0>0T_0>0, and uses ∣ξ∣2=τ2|\xi|^2=\tau^2. Such a T0T_0 exists if and only if ⟨b,ξ⟩>0\langle b,\xi\rangle>0: necessity follows directly from the strict inequality; for sufficiency choose T0>max⁡(0,T,∣b∣2/(2⟨b,ξ⟩))T_0>\max(0,T,|b|^2/(2\langle b,\xi\rangle)). This also guarantees the required future orientation. Since ⟨b,ξ⟩=Ha\langle b,\xi\rangle=H_a, the criterion is proved in full space. It implies t>a−∣q∣t>a-|q| by Cauchy--Schwarz. Every characteristic covector over t>a+∣q∣t>a+|q| satisfies it strictly. Equality Ha=0H_a=0 does not give a vertex in the open forward cone.

For the half-space, unfold the normal coordinate of a ray across r=0r=0. Simultaneously reflecting q1,ξ1q_1,\xi_1 preserves both the characteristic condition and HaH_a. A full-space future vertex with negative normal coordinate becomes a permitted half-space vertex on reflection, with the same norm. Conversely, unfolding any physical reflected segment produces exactly (WS34), possibly with that reflected vertex. At a hit r=0r=0, the normal term in ⟨q,ξ⟩\langle q,\xi\rangle is zero, so (WS33) descends to the compressed root identification. A tangent ray has zero normal velocity and remains in the boundary: its zero normal acceleration and the integrated normal equations force this continuation. Thus unfolding and the criterion apply to all flat generalized rays.

An open cone with genuine observation access

G4. The three base regions and the whole characteristic access. In Q={∣t∣<3,r≥0}Q=\{|t|<3,r\ge0\}, put a=1/2a=1/2 and

Ω1={t>a+∣q∣}∩Q,Ω2={∣q∣>∣t∣+a}∩Q,Ω3={t>a−∣q∣}∩Q.(WS36) \Omega_1=\{t>a+|q|\}\cap Q,\qquad \Omega_2=\{|q|>|t|+a\}\cap Q,\qquad \Omega_3=\{t>a-|q|\}\cap Q. \tag{WS36}
All cotangent cones below omit the zero section and use compressed covectors at the face. With Σ={p=0}\Sigma=\{p=0\}, define
Γj=Tb∗Ωj∖0 (j=1,2),Γ=((Tb∗Ω3∖0)∖Σ) ∪ Wa.(WS37) \Gamma_j=T_b^*\Omega_j\setminus0\ (j=1,2),\qquad \Gamma=\bigl((T_b^*\Omega_3\setminus0)\setminus\Sigma\bigr) \ \cup\ W_a. \tag{WS37}
This agrees with adjoining the whole cone over Ω1\Omega_1, since all its characteristic points belong to WaW_a. In the interior, the complement of Γ\Gamma relative to Tb∗Ω3∖0T_b^*\Omega_3\setminus0 is Σ∩{Ha≤0}\Sigma\cap\{H_a\le0\}, a closed set. At a boundary point, lift to its ordinary quadratic characteristic roots. Their compact normalized lifts vary continuously as an unordered pair; at a repeated root they coincide. A sequence of excluded compressed covectors therefore has a convergent excluded root lift, and at the face (WS33) has the same value on both representatives. This proves relative closedness there as well. Thus Γ\Gamma is an open conic set, rather than merely a relatively open part of the characteristic set.

Direct inequalities give Ω2⊂Ω3\Omega_2\subset\Omega_3 and Ω3∩{t≤0}=Ω2∩{t≤0}\Omega_3\cap\{t\le0\}=\Omega_2\cap\{t\le0\}. A characteristic covector over Ω2\Omega_2 need not lie in WaW_a. Accordingly the observation cone used for the access theorem is precisely Γ2obs=Γ∩(Tb∗Ω2∖0)\Gamma_2^{\rm obs}=\Gamma\cap(T_b^*\Omega_2\setminus0), which is open, lies in Γ\Gamma, and still reads only the zero region.

Every flat ray through a point (t0,q0)∈Ω1(t_0,q_0)\in\Omega_1 stays in Γ\Gamma as it is followed backward to s=−1s=-1. Indeed physical time is a global ray parameter since τ≠0\tau\ne0; HaH_a is constant between hits, continuous at each hit and unchanged on tangent segments. It remains positive, so the ray stays over Ω3\Omega_3. The explicit unfolded straight line and reflection give

∣q(s)∣≥t0−s−∣q0∣>a−s=∣s∣+a(−3<s≤0).(WS38) |q(s)|\ge t_0-s-|q_0|>a-s=|s|+a \qquad(-3<s\le0). \tag{WS38}
Hence the ray reaches Γ2obs\Gamma_2^{\rm obs} already at s=−1s=-1. No finite exit inside Γ\Gamma can precede that time: the displayed flat formula supplies its continuation. This verifies the maximal-ray access condition, including every reflected and tangent continuation. Only the cone Γ\Gamma is used; no smoothness on all of Ω3\Omega_3 is inferred. The later compact construction uses smaller neighborhoods of these rays, with the explicit positive margins in G0--G2.

At q=(1/2,0), the blue region consists exactly of characteristic directions that can return from the open forward cone.

V1. Read this covector slice. The horizontal coordinate is θ∈[−π,π]\theta\in[-\pi,\pi], with τ=1\tau=1 and ξ=(cos⁡θ,sin⁡θ)\xi=(\cos\theta,\sin\theta); the fixed spatial point is q=(1/2,0)q=(1/2,0). The vertical coordinate is T=t−aT=t-a. The boundary curve is exactly T=−cos⁡θ/2T=-\cos\theta/2; the blue region above it satisfies Ha>0H_a>0. The curve itself is excluded. This is a two-dimensional covector-direction slice, not a physical ray diagram. Above T=1/2T=1/2 every displayed direction is admitted; below T=−1/2T=-1/2 none is. The generator samples only this analytic curve for drawing; (WS33)--(WS35) prove the exact strict condition.

G0. Exact reference geometry, including reflection. At a0a^0 the principal wave operator is flat. In physical time, an unfolded characteristic from (t0,r0,y0)(t_0,r_0,y_0) with spatial velocity e=(er,ey)e=(e_r,e_y), ∣e∣=1|e|=1, is

r(s)=∣r0+(s−t0)er∣,y(s)=y0+(s−t0)ey.(WS19) r(s)=|r_0+(s-t_0)e_r|,\qquad y(s)=y_0+(s-t_0)e_y . \tag{WS19}

At an interior starting point either signed normal velocity is allowed; at r0=0r_0=0 both root representatives give the same reflected path with the corresponding tangential velocity. If er=0e_r=0 and r0=0r_0=0, (WS19) is the gliding path. The flat glancing equation has no other continuation: its normal acceleration and all contact derivatives vanish, and the integrated generalized-flow equations give r=0r=0, while the tangential motion is the displayed line. This is also the zero-comparison case PRS:S6.

Reflection preserves Euclidean norm. The reverse triangle inequality therefore gives, for s≤t0s\le t_0,

∣q(s)∣≥t0−s−∣q(t0)∣.(WS20) |q(s)|\ge t_0-s-|q(t_0)| . \tag{WS20}

For every target in MM, at s=−11/10s=-11/10 the right side exceeds ∣s∣+1/2|s|+1/2 by at least

99/100−21/50−1/2=7/100.(WS21) 99/100-21/50-1/2=7/100. \tag{WS21}

The same margin holds at s=−1s=-1. The entire segment from the target to s=−6/5s=-6/5 lies in ∣q∣≤21/50+201/100+6/5=363/100<4|q|\le21/50+201/100+6/5=363/100<4. These are bounds for all directions, including tangency, not samples of selected rays.

Let KMK_M be the complete unit compressed cotangent sphere over MM. Let K0K^0 be its union with every characteristic backward segment from its characteristic part down to t=−1t=-1. This set is compact. One proof uses the explicit formula (WS19) and compact unit directions, with the two reflected representatives identified by compression. Equivalently PRS:S2–S7 gives closedness of limits of all the bounded arcs. Noncharacteristic points of KMK_M contribute no arcs.

G1. Fix a parameter set before choosing derivatives. Choose an ambient open compressed cone Γ0\Gamma^0 in a slightly larger working box, containing K0K^0 and all the preceding segments down to −6/5-6/5. The reference paths from K0K^0 to t=−11/10t=-11/10 have positive support and observation margins by (WS21). The continuous flat reflected flow, or the compactness and exit-buffer proof PRS:S3–S10, gives an open neighborhood of K0K^0 from which every characteristic reaches OO backward while staying in Γ0\Gamma^0. Take two compactly nested smaller neighborhoods and call the inner one T1T_1. The closure of T1T_1 is a fixed compact normalized target set for CRO:H1.

There is a compact parameter neighborhood Z0Z_0 of a0a^0, chosen now and independently of any derivative order, such that all backward characteristic segments from KMK_M, down to t=−1t=-1, stay in T1T_1 for a∈Z0a\in Z_0, and their endpoints lie in OO with a uniform positive margin. To verify this stronger segment assertion, suppose otherwise and choose convergent parameters and initial states. Normalize characteristic covectors by ∣τ∣=1|\tau|=1. Spatial positive definiteness bounds every spatial component, and the physical-time component of the normalized Hamilton field is bounded away from zero. Thus physical time is a valid common parameter. On the larger compact box, PRS:S2–S7 supplies compactness and continuation; a first attempted exit is ruled out by its positive buffer. A subsequence would limit to a flat segment violating either the chosen neighborhood or (WS21), a contradiction. This also supplies the slight backward extension needed to reach −11/10-11/10 and proves the ray-access hypothesis used for T1T_1.

Shrink Z0Z_0 once more to the geometric parameter range in CRO:H1 for target T1‾\overline{T_1}, ambient cone Γ0\Gamma^0 and observation in OO. Equations (WS16)–(WS17), with zero source for vav_a, then give

sup⁡a∈Z0, aq>0sup⁡b∥Ca,bva∥2<∞ifWF⁡b′Ca,b⋐T1,(WS22) \sup_{a\in Z_0,\ a_q>0}\sup_b\|C_{a,b}v_a\|_2<\infty \quad\text{if}\quad \operatorname{WF}'_b C_{a,b}\Subset T_1, \tag{WS22}

for every prescribed bounded full family of finite order with the stated input margins. The finite norm bound depends on that family. In wave spacetime the tangential dimension is n≥1n\ge1, even for one spatial dimension, so the scalar assembly used here has its required positive-dimensional reserve. This does not invoke a scalar assembly theorem on an arbitrary compact boundary manifold.

6. Parameter differentiation on nested cones

G2. Nested neighborhoods for a fixed coefficient value. Fix a∗∈Z0a^*\in Z_0. Let K∗K^* consist of KMK_M and all backward Pa∗\mathcal P_{a^*} segments from its characteristic part to t=−1t=-1. The preceding compactness argument proves that K∗K^* is compact and lies in T1T_1; its section at t=−1t=-1 is compact in OO. It is closed under all further backward continuations until that plane. Indeed concatenating any such continuation with the segment from KMK_M is another permitted generalized characteristic; the endpoint law and compact continuation are the ones proved in GSP:C11 and PRS:S7. No unique continuation at infinite contact is selected.

Given a precompact open TνT_\nu containing K∗K^*, the compact section at t=−1t=-1 has a positive buffer in Tν∩Tb∗OT_\nu\cap T_b^*O. Uniform velocity bounds permit every endpoint continuation down to −1−ε-1-\varepsilon in that buffer, for some ε>0\varepsilon>0. Before t=−1t=-1, all the continuations lie in K∗K^*. The limit-arc and exit-buffer argument therefore gives a neighborhood Tν+1T_{\nu+1} with

K∗⊂Tν+1⋐Tν,every maximal ray in Tν meeting Tν+1 meets Tν∩Tb∗O.(WS23) K^*\subset T_{\nu+1}\Subset T_\nu, \quad\text{every maximal ray in }T_\nu \text{ meeting }T_{\nu+1}\text{ meets }T_\nu\cap T_b^*O. \tag{WS23}

Here is the contradiction step explicitly. If arbitrarily near initial states could leave TνT_\nu before the observation plane, take a subsequence converging to a state of K∗K^*. The bounded time interval, the nonzero physical-time velocity and PRS:S3–S7 give a limiting arc starting there. It lies in K∗K^* to t=−1t=-1, then in the endpoint buffer to t=−1−εt=-1-\varepsilon. Cover this compact limiting arc by finitely many smaller buffers inside TνT_\nu; uniform convergence and the local lifetime bound rule out the proposed exits. The same argument rules out missing OO at the final plane. Take Tν+1T_{\nu+1} sufficiently small that its starting times are above −1−ε/2-1-\varepsilon/2. Noncharacteristic starting states require no access assertion. This proves (WS23).

Choose the closure of Tν+1T_{\nu+1} inside the open access neighborhood just obtained. Apply CRO:H1 at the reference value a∗a^*, with that larger access neighborhood as its open target region, compact target Tν+1‾\overline{T_{\nu+1}}, ambient cone TνT_\nu, and observation Tν∩Tb∗OT_\nu\cap T_b^*O. Thus its compact target is strictly inside its open target region, as required. Use the observation refinement O0. The parameter radius is independent of the order or full target family. Only finitely many such applications will be used for any specified parameter derivative.

I0. The finite parameter induction. We claim that for every integer k≥0k\ge0, all ∣α∣≤k|\alpha|\le k and every bounded full target family with microsupport compact in Tk+1T_{k+1},

sup⁡a in a neighborhood of a∗a∈Z0, aq>0sup⁡b∥Ca,bwα,a∥2<∞.(WS24) \sup_{\substack{a\text{ in a neighborhood of }a^*\\a\in Z_0,\ a_q>0}} \sup_b\|C_{a,b}w_{\alpha,a}\|_2<\infty. \tag{WS24}

The neighborhood may depend on kk, but can be chosen independently of target order. For k=0k=0 this is (WS22). Assume it for k−1k-1, and take ∣α∣=k|\alpha|=k. Apply the chart estimate from TkT_k to Tk+1T_{k+1}, in a parameter neighborhood also contained in that of the induction hypothesis. Its coarse term is bounded by (WS16); its observation is zero by (WS17). Its source is the finite sum Af(∂aβRa)wα−β,aA_f(\partial_a^\beta R_a)w_{\alpha-\beta,a}. Each composed family in (WS18) has microsupport compact in TkT_k, is uniformly bounded in a finite full order, and retains the input margin. The induction hypothesis bounds its actual norm, including its residual. There are finitely many multiindices at order kk, so one parameter neighborhood and one finite coarse order work for all of them. Lower derivatives retain their bounds on the smaller cone. This proves (WS24).

The same reasoning applies to arbitrary bounded target families, rather than merely to one fixed test: that is part of CRO:H1. It is needed here because the coefficient derivatives in (WS15) vary with aa. No countable density assertion about all residual kernels is used.

7. Recover ordinary derivatives and keep one parameter neighborhood

J0. Tangential differential tests and normal recovery. Fix a finite parameter order kk. The full compressed sphere over MM lies in K∗⊂Tk+1K^*\subset T_{k+1}. Compactness of that sphere gives a base neighborhood of MM whose entire compressed sphere is still in Tk+1T_{k+1}. Choose finitely many nested cutoffs there. For each tangential multiindex γ\gamma, the actual operator χDt,zγ\chi D_{t,z}^\gamma is a compact boundary differential test of finite order, with microsupport in that sphere neighborhood. Thus (WS24) bounds every tangential derivative of every wα,aw_{\alpha,a}, ∣α∣≤k|\alpha|\le k, in local L2L^2. These are polynomial differential symbols, not an unjustified assertion that an arbitrary tangential pseudodifferential symbol belongs to the full boundary class at large normal frequency.

The full normal-recovery proof MH:R1–R2 now supplies the first two ordinary normal derivatives. For any integer s≥2s\ge2, use its indices

(r1,q1)=(0,s),(r2,q2)=(2,s−2),(a,b)=(2,s−2).(WS25) (r_1,q_1)=(0,s),\qquad (r_2,q_2)=(2,s-2), \qquad (a,b)=(2,s-2). \tag{WS25}

The seed is wα∈Lr2Ht,zsw_\alpha\in L^2_rH^s_{t,z}. Its source fαf_\alpha belongs to Lr2Ht,zs−2L^2_rH^{s-2}_{t,z} by (WS15) and the just-proved tangential bounds for lower parameter derivatives. All three recovery inequalities are satisfied, with total order ss. The receiving proof uses finitely many nested cutoffs and continuous coefficient multipliers, so its norm constants are uniform for the compact coefficient family. The integer-jet equivalence MH:N2 yields Drwα∈L2Hs−1D_rw_\alpha\in L^2H^{s-1} and Dr2wα∈L2Hs−2D_r^2w_\alpha\in L^2H^{s-2}, with uniform bounds. It also treats the interior parts of these patches.

Higher ordinary normal derivatives follow without a new microlocal theorem. Since Dr2=−∂r2D_r^2=-\partial_r^2, differentiate the exact equation to obtain

∂rj+2wα=−∑h=0j(jh)(∂rhRa)∂rj−hwα−∂rjfα.(WS26) \partial_r^{j+2}w_\alpha =-\sum_{h=0}^j {j\choose h} (\partial_r^h R_a)\partial_r^{j-h}w_\alpha -\partial_r^j f_\alpha . \tag{WS26}

For the finite collection ∣α∣≤k|\alpha|\le k, suppose all normal derivatives up to order j+1j+1, with arbitrary fixed tangential derivatives, have been bounded. Every term on the right of (WS26) involves at most jj normal derivatives of this same collection, because (WS15) is tangential. Smooth coefficient derivatives are uniformly bounded. Hence the next normal derivative is bounded as well. Starting with orders zero, one and two proves any requested finite ordinary jet. Only finitely many nested cutoffs are used for that jet, all around the same compact MM.

Finally the Fourier Sobolev estimate from NGR:B2 converts sufficiently high local ordinary Sobolev norms into uniform pointwise bounds, including the boundary. In the restricted norm choose a whole-space extension within twice the infimum; the whole-space estimate bounds its continuous derivatives. Their restrictions agree with the smooth interior derivatives, hence with the boundary limits. Taking the infimum proves the asserted local bound independently of the extension. This proves all the required mixed (a,t,q)(a,t,q) derivatives on MM.

J1. A finite cover for each derivative order. For a specified mixed derivative, J0 supplies a parameter neighborhood of each a∗∈Z0a^*\in Z_0 and a bound there. These open neighborhoods cover the same compact Z0Z_0. Choose a finite subcover and take the maximum of its finitely many constants. This gives a bound on all of Z0Z_0, for that derivative. Repeat for each derivative order. The set Z0Z_0 is never replaced by an intersection of the shrinking neighborhoods in (WS24).

Return through the exact inverse of (WS11), the fixed linear map and the parameter-dependent translation. The chain rule uses only finitely many of the derivatives just bounded, and all its coefficient and inverse-Jacobian derivatives are uniformly bounded. By (WS7), (WS10) and (WS12), the target in (WS4) is always mapped into the interior of MM. Thus the original derivatives have the same type of uniform bound.

Repeat the initial geometric construction at finitely many a0a^0 covering K0K_0. Choose a relative neighborhood K′K' of K0K_0 with closure in the union of the resulting compact parameter neighborhoods. These geometric choices preceded every derivative order. The maximum of the finitely many patch constants gives (WS4) on this one K′K'. This completes the coefficient-uniform wave smoothing statement, including the case n=1n=1.

8. Three examples with complete solutions

X1. A reflected ray and its observation margin. Take n=2n=2, target t0=3/2t_0=3/2, r0=1/5r_0=1/5, y0=1/10y_0=1/10, and unfolded unit velocity e=(3/5,4/5)e=(3/5,4/5). Find the boundary hit and the point at observation time −11/10-11/10.

Solution. The normal affine coordinate is 1/5+3(t−3/2)/51/5+3(t-3/2)/5, so its zero is t=7/6t=7/6. At that hit the tangential coordinate is −1/6-1/6. The compressed curve is continuous, with normal momenta interchanged and tangential velocity unchanged. At t=−11/10t=-11/10, the time difference is −13/5-13/5, and

r=34/25,y=−99/50,∣q∣2=577/100.(WS27) r=34/25,\qquad y=-99/50,\qquad |q|^2=577/100. \tag{WS27}

This is strictly outside the observation threshold ∣t∣+1/2=8/5|t|+1/2=8/5, whose square is 64/2564/25. The target has ∣q0∣=5/10<21/50|q_0|=\sqrt5/10<21/50. The uniform margin 7/1007/100 in (WS21) applies to every target and direction in the stated compact set; this particular ray has a larger margin.

X2. Why differentiating after normalization matters. Consider

La=∂t2−∂r2−2a∂r∂y−(1+a2)∂y2.(WS28) L_a=\partial_t^2-\partial_r^2-2a\partial_r\partial_y -(1+a^2)\partial_y^2. \tag{WS28}

Find a parameter-dependent coordinate change that removes all parameter dependence from the equation. Compare its parameter derivative with the derivative in the original coordinates.

Solution. Its spatial principal matrix is

Aa=(1aa1+a2)=BaBaTA_a=\left(\begin{smallmatrix}1&a\\a&1+a^2\end{smallmatrix}\right)=B_aB_a^T,

where Ba=(10a1)B_a=\left(\begin{smallmatrix}1&0\\a&1\end{smallmatrix}\right).

It is positive definite, with determinant one. Put y=z+ary=z+ar. The Jacobian is one, and the transformed full operator is

Dr2−Dt2+Dz2=Dr2−R,R=Dt2−Dz2.(WS29) D_r^2-D_t^2+D_z^2=D_r^2-R, \qquad R=D_t^2-D_z^2. \tag{WS29}

There is no density or gauge correction in this constant-coefficient example. If va(t,r,z)=ua(t,r,z+ar)v_a(t,r,z)=u_a(t,r,z+ar), then

∂ava=(∂aua+r∂yua)(t,r,z+ar).(WS30) \partial_a v_a=(\partial_a u_a+r\partial_yu_a)(t,r,z+ar). \tag{WS30}

The normalized equation for this derivative has zero source because ∂aR=0\partial_aR=0. Direct differentiation of the original equation instead produces the mixed normal-tangential term 2∂r∂yua+2a∂y2ua2\partial_r\partial_yu_a+2a\partial_y^2u_a. Formula (WS30) accounts for exactly the additional coordinate derivative. For a fixed flat wave ww, ua(t,r,y)=w(t,r,y−ar)u_a(t,r,y)=w(t,r,y-ar) gives va=wv_a=w, so the cancellation is immediate and exact. The metric distance from zero is r2+(y−ar)2\sqrt{r^2+(y-ar)^2}, since the shear is a half-space isometry for Aa−1A_a^{-1}. Both initial and target metric balls therefore become the same Euclidean balls.

X3. One spatial dimension and the differentiated source sign. Let Lc=∂t2−c2∂r2L_c=\partial_t^2-c^2\partial_r^2, with 99/100≤c≤101/10099/100\le c\le101/100. For smooth odd gg supported in [−2/5,2/5][-2/5,2/5], set

uc(t,r)=12{g(r+ct)+g(r−ct)},r≥0.(WS31) u_c(t,r)=\tfrac12\{g(r+ct)+g(r-ct)\},\qquad r\ge0. \tag{WS31}

Check the boundary condition, the empty interior-cone region and the normalized parameter equation.

Solution. Oddness gives uc(t,0)=0u_c(t,0)=0, direct differentiation gives Lcuc=0L_cu_c=0, and the initial velocity is zero. If 1≤t≤21\le t\le2, 0≤r≤2/50\le r\le2/5, both arguments have absolute value at least 99/100−2/5=59/100>2/599/100-2/5=59/100>2/5, so the function and all its c,t,rc,t,r derivatives vanish there. For this example the support and target are specified in fixed Euclidean coordinates; its initial metric radius is at most 2/(5c)2/(5c). The general metric-radius theorem above is not needed for this direct calculation.

Dividing by c2c^2 gives Pc=Dr2−c−2Dt2\mathcal P_c=D_r^2-c^{-2}D_t^2, so Rc=c−2Dt2R_c=c^{-2}D_t^2. Thus

Pc∂cuc=−2c−3Dt2uc.(WS32) \mathcal P_c\partial_cu_c=-2c^{-3}D_t^2u_c. \tag{WS32}

The sign agrees with (WS15). Physical time is the sole tangential coordinate, so the chart still has positive tangential dimension. This explains why the observation assembly used for the general proof also includes one spatial dimension.

9. Read the geometry in the figure

Left: spatial radius of the exact reflected ray in X1, the common support bound and the exterior observation boundary; the target box is marked. Right: its normal-coordinate projection, with the reflection at t=7/6.

V0. Exact interpretation. The left curve is the radius of the full two-dimensional ray in (WS19) with the data of X1. Reflection leaves that radius unchanged. The orange support envelope is c∗(2/5+∣t∣)+1/100c_*(2/5+|t|)+1/100, including the coordinate displacement bound; the dashed observation boundary is ∣t∣+1/2|t|+1/2. The blue rectangle is the radial projection of MM, not a set through which all rays pass. The right panel is only the normal projection r(t)r(t); its omitted tangential coordinate is exactly y(t)=1/10+4(t−3/2)/5y(t)=1/10+4(t-3/2)/5. The boundary hit and the observation sample are computed in (WS27). Analytic curves are sampled for rendering; the estimates and reflection law are proved above.

The wave smoothing theorem has the stated parameter and spatial-derivative scope. Whole-manifold scalar testing and general Airy mapping have separate geometric and analytic hypotheses. External human proof dependencies remain explicit.