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

Global boundary operators, dual action and elliptic parametrices

This component retains AN03-U034 Sections 2.1–2.5, 3.5 and 4.1. Original author: Codex, September 2026, CC0. Current exact proof connections and clarifications: AN-04 course-writing task and OpenAI Codex, 5 October 2026, CC0. All original mathematical displays remain unchanged.

Exact scope and prerequisite proofs

The geometry companion supplies all conventions and (GL1)–(GL24), including the smooth-family definition at the new face. The local boundary calculus supplies every ordered near/far product and adjoint remainder. Its conormal-action companion supplies full conormal continuity, residual receiving estimates and the negative-order no-gain examples. The positive-order proof supplies (PS1)–(PS17) on both original half-space Sobolev scales at every real index, including rectangular systems.

The intrinsic conormal symbol proof supplies the full coordinate law and quotient; the conormal test proof supplies all Besov and ordinary operator bounds. The dual conormal and jet proof supplies the definitions, (GA1)–(GA2), (GD1)–(GD12), weak density, actual trace and corrected intrinsic derivatives. Thus the retained references to those equations below are exact earlier proofs. Formal adjoints use the fixed Hermitian metrics and half-density pairing. The supported space includes boundary distributions; the restricted space is its actual interior quotient. All orders are real.

For the global chart patching use the complete locally finite partition PS5; on a boundary manifold its same relative-half-ball construction gives precompact charts with locally finite closures. Use the complete ordinary operator calculus for interior charts, all-real Hilbert and quotient duality, and Fourier measure theory for the integrals and limits. The approved antecedent is Hörmander III, Sections 18.2–18.3.

The compressed wave-front definition and its tangential consequences are subsequent work; the local parametrix below proves the complete elliptic operator identity needed there. It does not by itself restore the higher-order Cauchy lesson.

2. Proper global operators and conormal smoothing

2.1. The global kernel and its actual pushforward

Choose a smooth function tt on the stretched square, positive off its new face and vanishing simply on that face. Near a boundary diagonal point it can be the original (xn+yn)/2(x_n+y_n)/2. Define the order-mm class by kernels represented as

k=t−1/2H,H∈Im(X×^X,Δ^;Ω1/2),H flat on both original side faces.(GC1) k=t^{-1/2}H,\qquad H\in I^m(X\widehat\times X,\widehat\Delta;\Omega^{1/2}), \qquad H\text{ flat on both original side faces}. \tag{GC1}

Here conormality at the new face has precisely the smooth family meaning specified in the geometry companion. A bundle kernel additionally has values in Hom⁡(Ey,Fx)\operatorname{Hom}(E_y,F_x). Multiplication by its smooth frame changes is included in the conormal topology.

There is a well-defined pushforward despite the separate factor t−1/2t^{-1/2}. Indeed, in the original normal chart, let a test half-density have coefficient p(x,y)p(x,y). Its pullback has coefficient t1/2p(β(x′,y′,t,r))t^{1/2}p(\beta(x',y',t,r)), since the absolute normal determinant is tt. Thus the actual pairing is

⟨β∗k,p⟩=∫0∞⟨H(x′,y′,t,r),p(x′,t(1+r/2),y′,t(1−r/2))⟩x′,y′,r dt.(GC2) \langle\beta_*k,p\rangle =\int_0^\infty \left\langle H(x',y',t,r), p(x',t(1+r/2),y',t(1-r/2))\right\rangle_{x',y',r}\,dt. \tag{GC2}

The conormal family is a smooth distribution-valued function of tt. On a compact interval its pairing with the smooth compact test family is bounded by finitely many test derivatives and symbol seminorms, by (C4). The integral is consequently well defined and continuous. To verify compactness here, cover the given compact test support by finitely many relatively compact product charts. Over each closed smaller chart the resolved normal coordinate lies in the compact interval [−2,2][-2,2], and t=(xn+yn)/2t=(x_n+y_n)/2 is bounded. The coordinate description (GL4) therefore makes the inverse image compact. Away from the corner the blowdown is a diffeomorphism. Those facts reduce the general pairing to finitely many such charts. Equation (GC2) is the complete meaning of pairing kk with the pulled-back test half-density. Neither square-root factor is separately discarded, and no pairing of an arbitrary distribution with a nonsmooth test coefficient is asserted.

If tˉ\bar t is another admissible defining function, Hadamard's formula gives tˉ=ct\bar t=c t, with c>0c>0 smooth near the new face. Away from that face both functions are positive. Hence

Hˉ=tˉ1/2k=c1/2H.(GC3) \bar H=\bar t^{1/2}k=c^{1/2}H. \tag{GC3}

Multiplication by c1/2c^{1/2}, and by its smooth inverse, preserves the conormal class and side flatness. The kernel, its pairing and its order are independent of this choice.

Every compactly localized kernel of (GC1) is the original half-space kernel of a lacunary order-mm symbol by (GL19)–(GL21). Away from the lifted diagonal the localized resolved coefficient is smooth, and the same inverse gives a residual symbol. The latter assertion also applies between distinct boundary charts: coordinates in the two charts may be used as the two independent tangential variables in the smooth kernel; there is no conormal singularity to align. Near an actual diagonal point use the same coordinate chart on both factors. Interior charts give the ordinary conormal pseudodifferential kernel by (C18).

The local theorem on smooth inputs therefore proves a continuous map

A:Cc∞(X;Ω1/2⊗E)⟶C∞(X;Ω1/2⊗F).(GC4) A:C_c^\infty(X;\Omega^{1/2}\otimes E) \longrightarrow C^\infty(X;\Omega^{1/2}\otimes F). \tag{GC4}

For fixed input support and a fixed output compact set, only finitely many product charts occur. Each local estimate uses finitely many derivatives of the input, and the sum of these estimates is finite. This proves the stated LF-to-Fréchet continuity. Formula (GL24) supplies its boundary value, and the original formula (5.2) supplies every normal jet. Conversely the local quantizations have exactly (GC1) by (GL13)–(GL18), so this kernel description and the patched local definition determine the same class. We denote it by Ψbm(X;E,F)\Psi_b^m(X;E,F).

To check coordinate invariance at all orders, the stretched coordinate change is (GL5), its conormal change is (GL10), and its complete amplitude and determinant are (C19). Reducing that entire amplitude by (C5) gives the transformed symbol, with the full order-(m−N)(m-N) remainder after NN terms. All frequency factors created by a base derivative pair with a frequency derivative and retain order mm. Thus the operator class, not only its leading term, is coordinate invariant. The local construction and the pushforward above agree in the interior and then as distributions by the same conormal family pairing (GC2).

2.2. The symbol quotient and actual realization of every symbol

For a compressed symbol aa, the original principal half-density of the kernel is the class of

a(x′,t,ξ′,ρ)t−1/2∣dx′ dt dξ′ dρ∣1/2.(GC5) a(x',t,\xi',\rho)t^{-1/2} |dx'\,dt\,d\xi'\,d\rho|^{1/2}. \tag{GC5}

The singular half-density multiplying aa is the invariant absolute symplectic half-density in (GL12). The complete determinant law (C20) therefore gives a scalar symbol, or a section of Hom⁡(E,F)\operatorname{Hom}(E,F), on T~∗X\widetilde T^*X, of ordinary order mm, modulo order m−1m-1. Both (GL16) and (GL17) remain its actual lower-order kernel remainder. Define Sm(T~∗X)S^m(\widetilde T^*X) using every local compact base set, all base derivatives, and the frequency bound ⟨(ξ′,ρ)⟩m−∣α∣\langle(\xi',\rho)\rangle^{m-|\alpha|}. The fibre transformation (GL10) and its inverse have smooth bounded coefficients on each compact chart; the chain rule proves that this definition is intrinsic. The leading symbol map has kernel exactly Ψbm−1\Psi_b^{m-1}, by the local normal-form kernel assertion (C21).

We prove surjectivity while retaining the original prescribed symbol. Choose a locally finite cover by precompact boundary or interior charts UjU_j, with locally finite closures. Choose a partition φj\varphi_j with support in these charts and ψj∈Cc∞(Uj)\psi_j\in C_c^\infty(U_j) equal to one on a neighborhood of supp⁡φj\operatorname{supp}\varphi_j. For a given symbol aa, compactly localize a full representative aja_j in the jj-th frame so that it equals aa on that neighborhood. In a boundary chart take exactly the lacunary modification

(aj)ρ(x,ξ)=∫aj(x,ξ′,ξn−v)ρ(v) dv,ρ^∈Cc∞((−1/2,1)),ρ^=1 near 0.(GC6) (a_j)_\rho(x,\xi)=\int a_j(x,\xi',\xi_n-v)\rho(v)\,dv, \quad \widehat\rho\in C_c^\infty((-1/2,1)),\quad \widehat\rho=1\text{ near }0. \tag{GC6}

Keep aj=(aj)ρ+[aj−(aj)ρ]a_j=(a_j)_\rho+[a_j-(a_j)_\rho]: the bracket is residual by the full Taylor integral (4.8), and every original symbol seminorm remains controlled. Quantize the first term and set

A=∑jφjT(aj)ρψj(GC7) A=\sum_j\varphi_j T_{(a_j)_\rho}\psi_j \tag{GC7}

in boundary charts, using ordinary quantization for interior charts. The sum is locally finite and has order mm. Its two support projections are proper: above a compact base set only finitely many chart closures occur, and both factors in each term have compact support in that chart. The leading term is ∑jφja=a\sum_j\varphi_j a=a. Multiplication on the input has its complete product remainder, of order m−1m-1, from (8.3); ψj=1\psi_j=1 near the working diagonal. The residual bracket in (GC6) has not become an equality between the original full symbol and its modification; it is explicitly the error in this realization. This gives both quotients

Ψbm/Ψbm−1≃Sm(T~∗X)/Sm−1(T~∗X),(GC8) \Psi_b^m/\Psi_b^{m-1}\simeq S^m(\widetilde T^*X)/S^{m-1}(\widetilde T^*X), \tag{GC8}

with the same isomorphism when the operator spaces are restricted to proper support. The bundle version follows in the same charts, with the actual frame transitions and (GC5).

2.3. Proper composition, adjoints and all asymptotic remainders

Let A∈Ψbm(X;F,G)A\in\Psi_b^m(X;F,G) and B∈Ψbm′(X;E,F)B\in\Psi_b^{m'}(X;E,F) be properly supported. For every compact output set the support of AA meets only a compact set of intermediate variables, and that compact set meets only a compact input set under BB. The reverse argument starts at a compact input set. This proves proper support of the composed operator and permits all intermediate partitions to be finite on the compact sets under consideration.

We give the localization argument, since a residual kernel at the new face need not be a smooth kernel on the original square. Localize near an output-input diagonal point. A sufficiently small common chart contains both variables. Split the intermediate variable into a slightly larger part of that chart and its complement. On the complementary part both kernels are separated from their actual diagonal, so their resolved coefficients are smooth with side flatness. On the chart part use the original local ordered product theorem (8.1)–(8.3). Any localization of either factor off its own diagonal gives a residual symbol by (GL19)–(GL21). That theorem then gives a residual product: its finite-order bound holds for every real order assigned to the residual factor, while its complete far term (8.2) is residual as well.

For clarity, this also handles two different boundary charts. For a localized rectangular kernel R(x,y)R(x,y) whose resolved coefficient is smooth, with output and input coordinates taken from their respective charts, the literal inverse is

AR(x,z)=xnR(x,(x′−z′,xn(1−zn))),aR(x,ξ)=∫e−iz⋅ξAR(x,z) dz.(GC9a) A_R(x,z)=x_nR\bigl(x,(x'-z',x_n(1-z_n))\bigr), \qquad a_R(x,\xi)=\int e^{-iz\cdot\xi}A_R(x,z)\,dz. \tag{GC9a}

This is the same (GL19) inverse written at fixed xn>0x_n>0, with all chart half-density factors included in RR. Resolved smoothness, side flatness, compact chart cutoffs and the full far estimates of (GL19)–(GL21) make aRa_R residual and lacunary. The two coordinate tuples both range in Euclidean half spaces, so the original local product theorem applies to these rectangular symbols with the intermediate coordinate tuple used as its common integration variable. If only one factor has a diagonal singularity, use its chart for the intermediate variable and its adjacent variable; the other is represented by (GC9a). If neither factor has a singularity, either choice works. Near a fixed output-input point off the diagonal, use disjoint neighborhoods of those two points and split the intermediate variable into their neighborhoods and the complement. At least one factor in each term is then separated from its diagonal, and the preceding argument proves a residual result. This proves the required smoothness off the lifted diagonal, its family regularity at the new face, and its side flatness; it does not replace a corner residual by an ordinary smoothing kernel.

In the common chart the complete product symbol is c=c1+c2c=c_1+c_2, where c1c_1 and c2c_2 are the actual (8.1) and (8.2), with the original factor order aa before bb. The full expansion is

c∼∑α1α!∂ξαa(x,ξ)Dx′α′Dvαnb(x′,vxn,ξ′,vξn)∣v=1,c−∑∣α∣<N(⋯ )∈Sm+m′−N.(GC9) c\sim\sum_\alpha\frac1{\alpha!}\partial_\xi^\alpha a(x,\xi) D_{x'}^{\alpha'}D_v^{\alpha_n} b(x',v x_n,\xi',v\xi_n)\big|_{v=1}, \qquad c-\sum_{|\alpha|<N}(\cdots)\in S^{m+m'-N}. \tag{GC9}

The exact far contribution c2∈S−∞c_2\in S^{-\infty} is part of cc. All its seminorms and the remainder seminorms are controlled by finitely many seminorms of the two original factors, by that local proof. Hence

AB∈Ψbm+m′,σ(AB)=σ(A)σ(B).(GC10) AB\in\Psi_b^{m+m'},\qquad \sigma(AB)=\sigma(A)\sigma(B). \tag{GC10}

For matrices this is E→σ(B)F→σ(A)GE\xrightarrow{\sigma(B)}F \xrightarrow{\sigma(A)}G; no order is exchanged. The same argument proves that the residual class is a two-sided ideal among properly supported operators of any finite order.

Transposition of the stretched square exchanges x,yx,y and sends rr to −r-r, keeping t=(xn+yn)/2t=(x_n+y_n)/2. Its kernel half-density law is

HA∗(x′,y′,t,r)=HA(y′,x′,t,−r)∗.(GC11) H_{A^*}(x',y',t,r)=H_A(y',x',t,-r)^*. \tag{GC11}

Conormality and side flatness are preserved. The local adjoint (7.3), including the residual adjoint of a−aρa-a_\rho, gives a†−a∗∈Sm−1a^\dagger-a^*\in S^{m-1}. Thus

A∗∈Ψbm(X;F,E),σ(A∗)=σ(A)∗,(λA+μB)∗=λ‾A∗+μ‾B∗.(GC12) A^*\in\Psi_b^m(X;F,E),\quad \sigma(A^*)=\sigma(A)^*,\quad (\lambda A+\mu B)^*=\overline\lambda A^*+\overline\mu B^*. \tag{GC12}

Proper support is preserved by this exchange. All adjoints here use the fixed metrics and the original half-density pairing.

The calculus admits full asymptotic sums. Here is the construction needed later for parametrices. In each chart let aj∈Sm−ja_j\in S^{m-j} be the actual localized symbols to be summed. Choose ff zero on the unit frequency ball and one outside twice that ball. Choose RjR_j increasing so that f(ξ/Rj)ajf(\xi/R_j)a_j has each of the first jj seminorms in Sm−j+1S^{m-j+1} at most 2−j2^{-j}, on the first jj compact base sets. This is possible because the support has ∣ξ∣≥Rj|\xi|\ge R_j, giving the extra factor Rj−1R_j^{-1}; frequency derivatives of the cutoff have the same bound after the product rule. The exact sum

asum=∑j=0∞f(ξ/Rj)aj,asum−∑j<Naj∈Sm−N(GC13) a_{\mathrm{sum}}=\sum_{j=0}^\infty f(\xi/R_j)a_j, \qquad a_{\mathrm{sum}}-\sum_{j<N}a_j\in S^{m-N} \tag{GC13}

converges in the stated seminorms. For the remainder, its infinite tail lies in Sm−NS^{m-N}; each of its finitely many low-frequency differences is residual. Lacunarize the full sum by (GC6), retaining the new residual difference, and patch as in (GC7). Each finite off-diagonal discrepancy is residual by the localization argument above. Consequently the global operator differs from every prescribed finite operator sum by exactly the asserted lower-order class. This proves asymptotic completeness, including on the proper-support subspace, without omitting any finite term or its remainder.

2.4. The supported and restricted Sobolev maps

For m≥0m\ge0 and every real ss, the actual local theorem (PS1)–(PS17) gives the supported bound with loss mm. Smooth coordinate changes and frame multiplications on compact sets are bounded on every real Sobolev scale: integer bounds follow by the chain and product rules, negative integers by their exact adjoints with the Jacobian, and intermediate orders by the full dyadic estimate in the conormal test proof. Explicitly choose integers a<s<ba<s<b and frequency blocks Πj\Pi_j, j≥0j\ge0. For the localized coordinate/frame map TT, its two integer bounds give

2ls−js∥ΠlTΠj∥2→2≤C{2−(j−l)(s−a),j≥l,2−(l−j)(b−s),j<l.(GB1) 2^{ls-js}\|\Pi_lT\Pi_j\|_{2\to2} \le C \begin{cases} 2^{-(j-l)(s-a)},&j\ge l,\\ 2^{-(l-j)(b-s)},&j<l. \end{cases} \tag{GB1}

The right side is a summable sequence in l−jl-j. Applying the triangle inequality to the sum of its shifts on ℓ2\ell^2 bounds the output dyadic square norm by the input norm. Plancherel identifies those norms with HsH^s. Finite Fourier sums converge in HsH^s; their images converge in the same space and in distributions to the actual coordinate map. This proves the bound for the original distributional map, including negative noninteger ss. The inverse coordinate map obeys the same argument. Boundary coordinate changes extend to open neighborhoods, preserve the closed half-space locally, and their adjoints retain the full Jacobian; the same bounds therefore pass to supported subspaces and to restriction quotients by taking the infimum over representatives. On a fixed compact input set, proper support and a partition reduce AA to finitely many of those bounds. A rectangular residual term has every order, so use its original order-zero bound and the continuous inclusion Hs⊂Hs−mH^s\subset H^{s-m}, valid for m≥0m\ge0. Thus, for each compact KK, there is a compact K′K' and a constant depending on finitely many localized full-symbol seminorms such that

supp⁡u⊂K⟹supp⁡Au⊂K′,∥Au∥H˙s−m(K′)≤C∥u∥H˙s(K).(GC14) \operatorname{supp}u\subset K\quad\Longrightarrow\quad \operatorname{supp}Au\subset K',\qquad \|Au\|_{\dot H^{s-m}(K')} \le C\|u\|_{\dot H^s(K)}. \tag{GC14}

The supported action is the original transpose action (9.1), so it includes distributions supported at the boundary. Both localizations and coordinates keep those terms; they are not zeroed by a chosen extension. For a nonproper operator the same proof gives local output bounds for compact input.

The corresponding restricted map follows from the exact local quotient map (PS16)–(PS17). Local supported representatives give the bounded map; the local action preserves the full ideal of distributions supported only on the boundary by (9.2), so different representatives have the same interior output. Taking the infimum over representatives in each fixed chart gives the restricted norm bound. Patch the actual quotient maps using the same finite cutoffs. This proves H‾comps→H‾locs−m\overline H^s_{\mathrm{comp}}\to\overline H^{s-m}_{\mathrm{loc}}, and the compact-output assertion under proper support. No gain of −m-m is asserted when m<0m<0; the original residual examples prohibit such a claim.

2.5. The supported conormal class and the residual receiving map

For real mm, let Am(X)\mathcal A^m(X) consist of distributions supported in the closed manifold which, in every boundary chart and its extension, are conormal of order mm to xn=0x_n=0. Explicitly put κ=−m−n/4\kappa=-m-n/4 and require

Pu∈B2,∞,locκfor every P∈Diff⁡b(X),supp⁡u⊂X.(GA1) Pu\in B^\kappa_{2,\infty,\mathrm{loc}} \quad\text{for every }P\in\operatorname{Diff}_b(X), \qquad \operatorname{supp}u\subset X. \tag{GA1}

The topology uses all these local seminorms. Multiplication, coordinate changes and finite frame transitions are continuous on those Besov spaces by (C2)–(C3), and (GL9) identifies the intrinsic tangent fields. Moving a coefficient through a word in tangent fields leaves only shorter words with smooth coefficients. Hence (GA1) is intrinsic and agrees exactly with (C15), with the original shift −m−n/4-m-n/4.

The coordinate seminorms can use D′D' and Xn=xnDnX_n=x_nD_n. The precise comparison with the original weighted derivatives is

xnkDnk=∏j=0k−1(Xn+ij),Xnk=∑j=0kckj xnjDnj,(GA2) x_n^kD_n^k=\prod_{j=0}^{k-1}(X_n+ij),\qquad X_n^k=\sum_{j=0}^k c_{kj}\,x_n^jD_n^j, \tag{GA2}

where the first polynomial identity defines an invertible triangular change of basis and the second is its inverse. The first follows by Dnxn=xnDn−iD_n x_n=x_nD_n-i and induction: multiplying xnkDnkx_n^kD_n^k by Xn+ikX_n+ik on the right gives xnk+1Dnk+1x_n^{k+1}D_n^{k+1}. These identities retain every lower term and its factor ii. Tangential derivatives commute with XnX_n. On compact sets all smooth-coefficient tangent words reduce to these generators. In particular (GA1) makes uu smooth in the interior by repeated ordinary derivatives there and the local Sobolev estimate.

If A∈ΨbdA\in\Psi_b^d is proper, its local conormal preservation theorem 11.2, with the actual index κ=−m−n/4\kappa=-m-n/4, gives

A:Am(X;E)⟶Am(X;F)(GA3) A:\mathcal A^m(X;E)\longrightarrow\mathcal A^m(X;F) \tag{GA3}

continuously, for every real m,dm,d. Finite localized product estimates prove the global continuity as in (GC14). An off-diagonal resolved term is residual and obeys the same theorem in the rectangular chart. The order dd does not change the conormal index: the local proof factorizes each full tangent derivative of the original operator through an even-order totally characteristic differential operator and applies the order-zero Besov bound to its two complete factors. In particular neither the original frequency nor a positive-order remainder is omitted.

There is a stronger receiving statement for a residual operator R∈Ψb−∞R\in\Psi_b^{-\infty}. Fix an input compact set and N>0N>0. For every tangent word PP, the composition PRPR is still residual, by the local product formulas and (GC10). Its order-zero bound therefore gives

∥PRu∥H−N(K′)≤CP,N,K∥u∥H˙−N(K),Ru∈AN−n/4(X).(GA4) \|PRu\|_{H^{-N}(K')}\le C_{P,N,K}\|u\|_{\dot H^{-N}(K)}, \qquad Ru\in\mathcal A^{N-n/4}(X). \tag{GA4}

The inclusion H−N⊂B2,∞−NH^{-N}\subset B^{-N}_{2,\infty} is immediate from the dyadic ℓ2\ell^2 and ℓ∞\ell^\infty norms. Any compactly supported distribution of order LL belongs to H−NH^{-N} for N>L+n/2N>L+n/2: its Fourier transform is bounded by C⟨ξ⟩LC\langle\xi\rangle^L, and the weighted square integral converges for exactly that strict inequality. The same estimate is uniform on a family with fixed support and bounded distribution-order seminorm. Thus

R:E′(X)⟶A(X):=⋃m∈RAm(X),(GA5) R:\mathcal E'(X)\longrightarrow \mathcal A(X):=\bigcup_{m\in\mathbb R}\mathcal A^m(X), \tag{GA5}

with the explicit fixed-Sobolev-source continuity in (GA4). This is a conormal receiving map, not an assertion that a residual operator produces a smooth function at the boundary.

3.5. Proper totally characteristic operators on the dual class

Let B∈Ψbd(X;E,F)B\in\Psi_b^d(X;E,F) be properly supported. Its formal adjoint B∗B^* preserves Am\mathcal A^m continuously for every mm, by (GC12) and (GA3). Given a compact output test v∈Amv\in\mathcal A^m, proper support places B∗vB^*v in a compact set depending only on the support of vv; its full conormal seminorms are bounded by finitely many input seminorms. Define the actual dual action by

(Bu)(v)=u(B∗v).(GD13) (Bu)(v)=u(B^*v). \tag{GD13}

For a smooth boundary test, (GD3) and the preceding estimate prove Bu∈A′Bu\in\mathcal A'. For smooth uu, the local adjoint identity (7.6) makes (GD13) the original kernel action. For general uu, (GD5) and weak continuity extend that equality; in the interior it matches the usual distributional operator. The coefficient order and bundle maps remain those of the original BB, with no scalar commutation of matrices.

The boundary jet formula extends as well. In a local half-space chart, let a∈Slada\in S^d_{\mathrm{la}}, and define the kk-th interior normal derivative of uu by iterating (GD9), then taking (GD8). For smooth uu, the exact formula is

(∇nint,kTau)∣xn=0=∑j=0k(kj)akj(x′,D′)(∇nint,ju∣xn=0),akj(x′,ξ′)=∑i=0j(ji)(Dxnk−jDξnia)(x′,0,ξ′,0).(GD14) \begin{aligned} \bigl(\nabla_n^{\mathrm{int},k}T_au\bigr)|_{x_n=0} &=\sum_{j=0}^k\binom{k}{j} a_{kj}(x',D') \bigl(\nabla_n^{\mathrm{int},j}u|_{x_n=0}\bigr),\\ a_{kj}(x',\xi') &=\sum_{i=0}^j\binom{j}{i} \bigl(D_{x_n}^{k-j}D_{\xi_n}^{i}a\bigr) (x',0,\xi',0). \end{aligned} \tag{GD14}

Every coefficient in the inner sum and the outer binomial factor is retained. Its ii-th summand has order d−id-i. The tangent boundary operator akj(x′,D′)a_{kj}(x',D') acts continuously on boundary distributions after compact localization; this is the ordinary local symbol action. The left and right sides of (GD14) are weakly continuous in uu by (GD8), (GD9), (GD13) and the boundary operator continuity. Smooth functions are weakly dense by (GD5), so (GD14) holds for every u∈A′u\in\mathcal A'. It is a statement about their actual boundary traces, not about an arbitrary supported representative's raw distributional normal derivatives.

4. Compressed wave fronts

4.1. Ellipticity and the characteristic set at every symbol order

Let B∈Ψbm(X;E,F)B\in\Psi_b^m(X;E,F) and let bb be its complete local symbol. At a nonzero compressed covector q=(x0,ζ0)q=(x_0,\zeta_0), call BB elliptic if the two bundle ranks agree and there are a base neighborhood UU, an open cone Γ\Gamma containing ζ0\zeta_0, constants c,R>0c,R>0, and local frames such that

b(x,ζ):Ex→Fx is invertible,∥b(x,ζ)−1∥≤c−1⟨ζ⟩−m(x∈U, ζ∈Γ, ∣ζ∣≥R).(GW1) b(x,\zeta):E_x\to F_x\text{ is invertible},\qquad \|b(x,\zeta)^{-1}\|\le c^{-1}\langle\zeta\rangle^{-m} \quad(x\in U,\ \zeta\in\Gamma,\ |\zeta|\ge R). \tag{GW1}

The set of covectors where this fails is Char⁡B\operatorname{Char}B. Changing the representative by Sm−1S^{m-1} does not change (GW1): b0−1(b−b0)b_0^{-1}(b-b_0) is O(⟨ζ⟩−1)O(\langle\zeta\rangle^{-1}), so for large ∣ζ∣|\zeta| its ordered Neumann series is invertible. Explicitly for a matrix EE with ∥E∥<1\|E\|<1, the partial sums of ∑ν≥0(−E)ν\sum_{\nu\ge0}(-E)^\nu are Cauchy in the finite-dimensional matrix norm, and multiplication by I+EI+E leaves the error (−E)N+1(-E)^{N+1}, whose norm tends to zero. The limit is the two-sided inverse. Smooth frame changes conjugate or left/right multiply by uniformly invertible matrices on a smaller compact chart. Thus the definition is intrinsic. The complement of the characteristic set is open and conic; the set itself is closed and conic in T~∗X∖0\widetilde T^*X\setminus0.

The inverse in (GW1) has the exact order −m-m on a smaller cone. Differentiate b−1b=IEb^{-1}b=I_E:

∂ζjb−1=−b−1(∂ζjb)b−1,∂xjb−1=−b−1(∂xjb)b−1.(GW2) \partial_{\zeta_j}b^{-1} =-b^{-1}(\partial_{\zeta_j}b)b^{-1},\qquad \partial_{x_j}b^{-1} =-b^{-1}(\partial_{x_j}b)b^{-1}. \tag{GW2}

Repeated product differentiation keeps the original matrix order. Every frequency derivative lowers the order by one; every base derivative leaves it unchanged. Multiplying by a conic cutoff gives a global chart symbol of order −m-m. This proves the symbol estimate needed for an actual microlocal parametrix.

Here is the complete construction. Choose nested cones Γ0⋐Γ1⋐Γ\Gamma_0\Subset\Gamma_1\Subset\Gamma and base neighborhoods U0⋐U1⋐UU_0\Subset U_1\Subset U, with a smooth large-frequency cutoff χ\chi supported in U1×Γ1U_1\times\Gamma_1 and equal to one on U0×Γ0U_0\times\Gamma_0 for ∣ζ∣|\zeta| large. Set c0=χb−1c_0=\chi b^{-1}, with the original map F→EF\to E. The lacunary realization (GC6)–(GC7) changes this by a residual symbol only. The ordered product gives

c0#b=χIE+e1,e1∈S−1,e1 residual away from the elliptic working cone.(GW3) c_0\# b=\chi I_E+e_1,\qquad e_1\in S^{-1}, \quad e_1\text{ residual away from the elliptic working cone}. \tag{GW3}

Suppose the sum c0+⋯+cN−1c_0+\cdots+c_{N-1} has error eN∈S−Ne_N\in S^{-N}, residual away from that cone. Its next correction is

cN=−χ1eNb−1∈S−m−N,(eN+cNb)=0where χ1=1 on the nonresidual support of eN.(GW4) c_N=-\chi_1e_N b^{-1}\in S^{-m-N}, \qquad (e_N+c_N b)=0 \quad\text{where }\chi_1=1\text{ on the nonresidual support of }e_N. \tag{GW4}

The actual product cN#b−cNbc_N\# b-c_Nb is one order lower by (GC9), so the new error is in S−N−1S^{-N-1}, still residual off the elliptic cone. One fixed larger cutoff suffices: take χ1=1\chi_1=1 near the closed base-angular support of χ\chi, inside the elliptic chart. The finite local product formula shows that away from this support every term of e1e_1 is zero to arbitrary symbol order, because every term contains a derivative of c0c_0; its full remainder can be assigned arbitrarily negative order. Inductively the same is true for eNe_N and cNc_N. In the transition region of χ1\chi_1 these errors are already residual, so its derivatives create only residual terms. Low frequencies are compact and hence residual. This justifies the stated cancellation with no division outside the elliptic region.

The cutoff χ1\chi_1 is supported inside the cone where b−1b^{-1} exists. At each stage the part not cancelled by χ1=1\chi_1=1 was already residual; include it in the final residual instead of dividing it by bb. Asymptotic summation (GC13) produces a single properly supported C∈Ψb−m(F,E)C\in\Psi_b^{-m}(F,E) with the global localized identity

CB=Q+R,Q=Op⁡(χIE)∈Ψb0(E,E),R∈Ψb−∞,(GW5) CB=Q+R, \qquad Q=\operatorname{Op}(\chi I_E)\in\Psi_b^0(E,E), \quad R\in\Psi_b^{-\infty}, \tag{GW5}

and QQ is elliptic on U0×Γ0U_0\times\Gamma_0. The low-frequency part and all coordinate-patching errors are residual, and every finite error is retained until its correction in (GW4). The same construction on bc=IFb c=I_F yields a right parametrix when needed; it is a separate ordered calculation, not inferred by commuting the matrices in (GW4).