Contents

Scaled Weyl parametrices and the surviving differential degree

Written and dedicated to the public domain by Codex, September 2026 (CC0).

The index can be read from two ordered errors of a cutoff inverse. Scaling the original matrix symbol makes the analytic remainders small enough to expose a finite differential coefficient, but the cutoff part of each error must stay visible. We establish Fredholmness at every positive scale, retain both matrix orders, prove a uniform trace bound only for genuine high-order remainders, and then show exactly why the lower differential degrees cancel in the difference of traces. The surviving coefficient is still a finite Weyl expression; its exterior and boundary evaluation is a further calculation.

The prerequisites are Weyl products for a varying metric for the ordered CkC_k coefficients and remainders, Metric operator bounds for boundedness, Symbols, operators and Sobolev scales for the operator product, Traces that survive passage to cohomology for the powers-of-errors index identity, and Weyl kernels, operator traces, and a finite trace-class test for the exact matrix trace and finite weighted-derivative test. All objects below retain their original x,ξx,\xi coordinates and the original isotropic metric.

1. Original metric, symbol and cutoff

Let n,ν≥1n,\nu\ge1, z=(x,ξ)∈ℝ2nz=(x,\xi)\in\mathbb R^{2n}, and e=|dx|2+|dξ|2e=|dx|^2+|d\xi|^2. Preserve the original metric and its weight: h(z)=(1+|x|2+|ξ|2)−1,gz=h(z)e,a∈S(1,g;End⁡ℂν).(IP1) h(z)=(1+|x|^2+|\xi|^2)^{-1},\qquad g_z=h(z)e,\qquad a\in S(1,g;\operatorname{End}\mathbb C^\nu). \tag{IP1} Assume that a(z)a(z) is invertible outside an open ball BB, with sup⁡z∉B∥a(z)−1∥<∞\sup_{z\notin B}\|a(z)^{-1}\|<\infty. Choose a real scalar ψ∈C∞(ℝ2n)\psi\in C^\infty(\mathbb R^{2n}) which vanishes near the noninvertible region, whose support lies where aa is invertible, and which equals one outside a compact set. Put b(z)={ψ(z)a(z)−1,z∈supp⁡ψ,0,where ψ=0.ba=ab=ψIν.(IP2) b(z)= \begin{cases}\psi(z)a(z)^{-1},&z\in\operatorname{supp}\psi,\\ 0,&\text{where }\psi=0.\end{cases} \qquad ba=ab=\psi I_\nu . \tag{IP2} The two pointwise equalities keep their displayed matrix order. On the invertibility region, differentiating a−1a=Iνa^{-1}a=I_\nu gives ∂ja−1=−a−1(∂ja)a−1\partial_j a^{-1}=-a^{-1}(\partial_j a)a^{-1}. Induction over a multiindex expresses every higher derivative as a finite sum of ordered words alternating bounded inverse factors with derivatives of aa, whose derivative orders add to the total order. The defining estimate for (IP1) makes a derivative of total order kk at most Ckhk/2C_k h^{k/2}. Derivatives of ψ\psi are supported in a fixed compact region, on which hh has a positive lower bound. Consequently b∈S(1,g;End⁡ℂν).(IP3) b\in S(1,g;\operatorname{End}\mathbb C^\nu). \tag{IP3} This is the actual cutoff inverse; no global inverse of aa is assumed.

2. The scaled metric with uniform structural constants

For 0<ε≤10<\varepsilon\le1, define all four scaled objects: aε(z)=a(εz),bε(z)=b(εz),ψε(z)=ψ(εz),hε(z)=ε2h(εz)=1ε−2+|z|2,gε,z=hε(z)e.(IP4) a_\varepsilon(z)=a(\varepsilon z),\quad b_\varepsilon(z)=b(\varepsilon z),\quad \psi_\varepsilon(z)=\psi(\varepsilon z),\quad h_\varepsilon(z)=\varepsilon^2h(\varepsilon z) ={1\over\varepsilon^{-2}+|z|^2},\quad g_{\varepsilon,z}=h_\varepsilon(z)e . \tag{IP4} The Euclidean symplectic dual of gεg_\varepsilon is gεσ=hε−1eg_\varepsilon^\sigma=h_\varepsilon^{-1}e; its Planck weight is exactly hε≤1h_\varepsilon\le1. Write Lε(z)=ε−2+|z|2L_\varepsilon(z)=\varepsilon^{-2}+|z|^2. If gε,z(v)≤δ2g_{\varepsilon,z}(v)\le\delta^2 with 0<δ<10<\delta<1, then |v|≤δLε(z)|v|\le\delta\sqrt{L_\varepsilon(z)}. The triangle inequality in ℝ2n+1\mathbb R^{2n+1}, applied to (ε−1,z)(\varepsilon^{-1},z), yields (1−δ)2Lε(z)≤Lε(z+v)≤(1+δ)2Lε(z).(IP5) (1-\delta)^2L_\varepsilon(z) \le L_\varepsilon(z+v) \le(1+\delta)^2L_\varepsilon(z). \tag{IP5} This proves slow variation uniformly in ε\varepsilon. For any z,wz,w, |z|≤|w|+|z−w||z|\le|w|+|z-w| gives Lε(z)≤2Lε(w)+2|z−w|2L_\varepsilon(z)\le2L_\varepsilon(w)+2|z-w|^2, and the same inequality holds with z,wz,w exchanged. Since Lε(z),Lε(w)≥1L_\varepsilon(z),L_\varepsilon(w)\ge1, both ratios of Lε(z)L_\varepsilon(z) and Lε(w)L_\varepsilon(w) are bounded by a fixed multiple of 1+Lε(w)|z−w|2=1+gε,wσ(z−w)1+L_\varepsilon(w)|z-w|^2=1+g_{\varepsilon,w}^\sigma(z-w). Thus gεg_\varepsilon and each fixed power of hεh_\varepsilon have uniform symplectic-temperateness and weight constants. These are the exact hypotheses needed for the one-metric specialization of W31–W33 and the boundedness estimate B26.

For a multiindex α\alpha of length kk, the chain rule and (IP1),(IP3) give ∥∂zαaε(z)∥+∥∂zαbε(z)∥≤Cαεk(1+ε2|z|2)−k/2=Cαhε(z)k/2.(IP6) \|\partial_z^\alpha a_\varepsilon(z)\| +\|\partial_z^\alpha b_\varepsilon(z)\| \le C_\alpha\varepsilon^k(1+\varepsilon^2|z|^2)^{-k/2} =C_\alpha h_\varepsilon(z)^{k/2}. \tag{IP6} The constants do not depend on 0<ε≤10<\varepsilon\le1. The same argument treats ψε\psi_\varepsilon, so these three families are uniformly bounded in their respective S(1,gε)S(1,g_\varepsilon) seminorms.

3. Both ordered Weyl errors and their cutoff parts

The one-metric Weyl product W31–W33 with N=1N=1, applied in each factor order, gives d1,ε=bε#aε−ψεIν∈S(hε,gε),d2,ε=aε#bε−ψεIν∈S(hε,gε),(IP7) d_{1,\varepsilon} =b_\varepsilon\#a_\varepsilon-\psi_\varepsilon I_\nu \in S(h_\varepsilon,g_\varepsilon),\qquad d_{2,\varepsilon} =a_\varepsilon\#b_\varepsilon-\psi_\varepsilon I_\nu \in S(h_\varepsilon,g_\varepsilon), \tag{IP7} with all indicated symbol seminorms uniform in ε\varepsilon. Here C0(bε,aε)=bεaε=ψεIνC_0(b_\varepsilon,a_\varepsilon)=b_\varepsilon a_\varepsilon=\psi_\varepsilon I_\nu, and C0(aε,bε)=aεbε=ψεIνC_0(a_\varepsilon,b_\varepsilon)=a_\varepsilon b_\varepsilon=\psi_\varepsilon I_\nu; no interchange of derivative factors in the two noncommutative remainders is made.

Let Aε=aεwA_\varepsilon=a_\varepsilon^w and Bε=bεwB_\varepsilon=b_\varepsilon^w on L2(ℝn;ℂν)L^2(\mathbb R^n;\mathbb C^\nu). B26 makes both bounded. Their symbols have bounded derivatives of every Euclidean order, because hε≤1h_\varepsilon\le1, so they also belong to the Euclidean order-zero calculus. For compact smooth approximants, Fubini gives the Weyl operator product identity. Both W31 and the Euclidean operator-composition theorem E19 pass that identity through bounded symbol families converging locally smoothly; uniqueness of the Weyl kernel identifies the two products. Hence the actual operator errors are R1,ε=I−BεAε=(r1,ε)w,R2,ε=I−AεBε=(r2,ε)w,rj,ε=(1−ψε)Iν−dj,ε.(IP8) R_{1,\varepsilon}=I-B_\varepsilon A_\varepsilon =(r_{1,\varepsilon})^w,\qquad R_{2,\varepsilon}=I-A_\varepsilon B_\varepsilon =(r_{2,\varepsilon})^w, \quad r_{j,\varepsilon} =(1-\psi_\varepsilon)I_\nu-d_{j,\varepsilon}. \tag{IP8} In particular, the uniform assertion is exactly rj,ε−(1−ψε)Iν∈S(hε,gε)uniformly for 0<ε≤1.(IP9) r_{j,\varepsilon}-(1-\psi_\varepsilon)I_\nu \in S(h_\varepsilon,g_\varepsilon) \quad\text{uniformly for }0<\varepsilon\le1 . \tag{IP9} It is generally false to claim that the whole rj,εr_{j,\varepsilon} has a uniform S(hε,gε)S(h_\varepsilon,g_\varepsilon) seminorm. If ψ(0)=0\psi(0)=0, then 1−ψε(0)=11-\psi_\varepsilon(0)=1 while hε(0)=ε2h_\varepsilon(0)=\varepsilon^2.

For each fixed ε>0\varepsilon>0, 1−ψε1-\psi_\varepsilon is smooth and supported in |z|≤C/ε|z|\le C/\varepsilon. On that support hε(z)≥cε>0h_\varepsilon(z)\ge c_\varepsilon>0, and its derivatives obey (IP6). It therefore belongs to S(hεM,gε)S(h_\varepsilon^M,g_\varepsilon) for every fixed integer MM, with constants allowed to depend on ε\varepsilon. In particular (IP7)–(IP8) imply r1,ε,r2,ε∈S(hε,gε)for each fixed ε>0,(IP10) r_{1,\varepsilon},r_{2,\varepsilon} \in S(h_\varepsilon,g_\varepsilon) \quad\text{for each fixed }\varepsilon>0, \tag{IP10} without converting (IP10) into a false uniform claim.

4. Trace-class powers with the full finite threshold

For a fixed ε>0\varepsilon>0, the weight hεh_\varepsilon is admissible for the product theorem. Repeated W31 products and the operator identity in Section 3 therefore give, in each original factor order, Rj,εN=(rj,ε#N)w,rj,ε#N∈S(hεN,gε)(N∈ℕ).(IP11) R_{j,\varepsilon}^{\,N} =\bigl(r_{j,\varepsilon}^{\# N}\bigr)^w,\qquad r_{j,\varepsilon}^{\# N} \in S(h_\varepsilon^N,g_\varepsilon) \quad(N\in\mathbb N). \tag{IP11} There is no uniform-in-ε\varepsilon seminorm claim for the full power. To apply the exact CT1 criterion, fix N≥n+1N\ge n+1, put q=rj,ε#Nq=r_{j,\varepsilon}^{\# N}, and take polynomial weight degree d=|α|+|β|d=|\alpha|+|\beta| and derivative degree k=|α′|+|β′|k=|\alpha'|+|\beta'| with d+k≤n+1d+k\le n+1. Because hε(z)≍ε⟨z⟩−2h_\varepsilon(z)\asymp_\varepsilon\langle z\rangle^{-2} for this fixed ε\varepsilon, its symbol estimate gives |xαξβ∂xα′∂ξβ′q(z)|≤Cε,α,β,α′,β′⟨z⟩d−2N−k.(IP12) |x^\alpha\xi^\beta \partial_x^{\alpha'}\partial_\xi^{\beta'}q(z)| \le C_{\varepsilon,\alpha,\beta,\alpha',\beta'} \langle z\rangle^{d-2N-k}. \tag{IP12} The square of the right side is integrable in 2n2n dimensions: the smallest decay exponent 2N+k−d2N+k-d occurs at k=0,d=n+1k=0,d=n+1, where it is 2N−(n+1)≥n+1>n2N-(n+1)\ge n+1>n. Every term of the original four-multiindex CT1 sum is finite. CT1–CT2 thus prove that both R1,εNR_{1,\varepsilon}^N and R2,εNR_{2,\varepsilon}^N are trace class. The exact CI12/CT15 trace identity also applies because q∈L1q\in L^1: Tr⁡Rj,εN=(2π)−n∫ℝ2ntr⁡ℂν(rj,ε#N(x,ξ))dxdξ.(IP13) \operatorname{Tr}R_{j,\varepsilon}^{\,N} =(2\pi)^{-n}\int_{\mathbb R^{2n}} \operatorname{tr}_{\mathbb C^\nu} \bigl(r_{j,\varepsilon}^{\# N}(x,\xi)\bigr) \,dx\,d\xi . \tag{IP13} The matrix trace, both phase coordinates, and the inverse Fourier factor remain explicit.

5. Fredholmness, power index and scaling continuity

Apply the proved powers-of-errors identity T28 to the bounded pair T=Aε,S=BεT=A_\varepsilon,S=B_\varepsilon and the two trace-class powers in Section 4. It proves, for every 0<ε≤10<\varepsilon\le1, that AεA_\varepsilon is Fredholm and ind⁡Aε=Tr⁡R1,εN−Tr⁡R2,εN=(2π)−n∫ℝ2ntr⁡ℂν(r1,ε#N−r2,ε#N)(z)dz,N≥n+1.(IP14) \operatorname{ind}A_\varepsilon =\operatorname{Tr}R_{1,\varepsilon}^{\,N} -\operatorname{Tr}R_{2,\varepsilon}^{\,N} =(2\pi)^{-n}\int_{\mathbb R^{2n}} \operatorname{tr}_{\mathbb C^\nu} \bigl(r_{1,\varepsilon}^{\# N} -r_{2,\varepsilon}^{\# N}\bigr)(z)\,dz , \qquad N\ge n+1. \tag{IP14} This does not assume the two errors commute. At ε=1\varepsilon=1, A1=awA_1=a^w, so (IP14) already establishes Fredholmness of the original operator.

For completeness, index constancy along the entire positive scaling interval follows from actual operator-norm continuity, not an asserted limit at ε=0\varepsilon=0. Fix ε0>0\varepsilon_0>0 and restrict ε\varepsilon to a compact interval J⊂(0,1]J\subset(0,1] around it. The weights hεh_\varepsilon and hε0h_{\varepsilon_0} are uniformly comparable there. Differentiating an order-kk spatial symbol derivative gives ∂ε∂zαa(εz)=kεk−1(∂αa)(εz)+εkz⋅∇(∂αa)(εz).(IP15) \partial_\varepsilon\partial_z^\alpha a(\varepsilon z) =k\varepsilon^{k-1}(\partial^\alpha a)(\varepsilon z) +\varepsilon^k z\cdot\nabla(\partial^\alpha a)(\varepsilon z). \tag{IP15} The first term is bounded by CJhε0k/2C_Jh_{\varepsilon_0}^{k/2}; the second has the same bound because ε|z|/1+ε2|z|2≤1\varepsilon|z|/\sqrt{1+\varepsilon^2|z|^2}\le1 and ε−1\varepsilon^{-1} is bounded on JJ. The fundamental theorem of calculus in every finite S(1,gε0)S(1,g_{\varepsilon_0}) seminorm and B26 imply Aε→Aε0A_\varepsilon\to A_{\varepsilon_0} in operator norm. The Fredholm index is locally constant under that topology, hence ind⁡aεw=ind⁡awfor every 0<ε≤1.(IP16) \operatorname{ind}a_\varepsilon^w =\operatorname{ind}a^w \quad\text{for every }0<\varepsilon\le1. \tag{IP16} Nothing here states norm continuity at ε=0\varepsilon=0, where the rescaled symbol need not have such a limit.

6. The exact metric-weight integral and its limit

The uniformly controlled symbolic remainder in (IP9) involves powers of hεh_\varepsilon. Its phase-space integral has a separate exact scaling, for every real M>nM>n: ∫ℝ2nhε(z)Mdz=ε2M−2n∫ℝ2n(1+|w|2)−Mdw=ε2M−2nπnΓ(M−n)Γ(M).(IP17) \begin{aligned} \int_{\mathbb R^{2n}}h_\varepsilon(z)^M\,dz &=\varepsilon^{2M-2n} \int_{\mathbb R^{2n}}(1+|w|^2)^{-M}\,dw\\ &=\varepsilon^{2M-2n} \pi^n{\Gamma(M-n)\over\Gamma(M)}. \end{aligned} \tag{IP17} The first equality is the exact change w=εzw=\varepsilon z, whose 2n2n-dimensional Jacobian is ε−2n\varepsilon^{-2n}. For the second, insert (1+|w|2)−M=Γ(M)−1∫0∞tM−1e−te−t|w|2dt(1+|w|^2)^{-M}=\Gamma(M)^{-1} \int_0^\infty t^{M-1}e^{-t}e^{-t|w|^2}\,dt. Tonelli applies to the nonnegative integrand, the 2n2n-dimensional Gaussian integral is (π/t)n(\pi/t)^n, and the remaining integral is Γ(M−n)\Gamma(M-n).

Equation (IP17) is not a trace-norm bound for the full powers in (IP14). Those powers contain the compact cutoff contribution 1−ψε1-\psi_\varepsilon, whose S(hε)S(h_\varepsilon) seminorm is not uniform. The later differential-degree extraction must keep that contribution and prove its cancellations, as well as any separate trace estimate for uniformly controlled symbolic remainders. The top-degree alternating form, boundary Stokes sign and final integral of target 051 remain to be derived.

7. Uniform trace bounds for genuine symbolic remainders

There is nevertheless an exact trace estimate for a remainder that is uniform in its full weight. For either order of the original matrix factors, W32–W33 give, for an integer M≥n+1M\ge n+1, q1,ε(M)=bε#aε−∑k=0M−1Ck(bε,aε),q2,ε(M)=aε#bε−∑k=0M−1Ck(aε,bε),qj,ε(M)∈S(hεM,gε)with seminorms uniform in ε.(IP18) \begin{split} q_{1,\varepsilon}^{(M)} &=b_\varepsilon\#a_\varepsilon -\sum_{k=0}^{M-1}C_k(b_\varepsilon,a_\varepsilon),\\ q_{2,\varepsilon}^{(M)} &=a_\varepsilon\#b_\varepsilon -\sum_{k=0}^{M-1}C_k(a_\varepsilon,b_\varepsilon),\\ q_{j,\varepsilon}^{(M)} &\in S(h_\varepsilon^M,g_\varepsilon) \quad\text{with seminorms uniform in }\varepsilon. \end{split} \tag{IP18} The coefficients CkC_k are the actual ordered Weyl differential coefficients W29; in the matrix case the order of factors in each coefficient stays as written. By the CT1 argument used in (IP12), each qj,ε(M)wq_{j,\varepsilon}^{(M)w} is trace class for fixed ε\varepsilon. Its zeroth symbol seminorm and the trace identity then give the uniform trace bound |Tr(qj,ε(M))w|≤(2π)−nνCM∫ℝ2nhε(z)Mdz≤CM′ε2M−2n.(IP19) \left|\operatorname{Tr}(q_{j,\varepsilon}^{(M)})^w\right| \le(2\pi)^{-n}\nu C_M \int_{\mathbb R^{2n}}h_\varepsilon(z)^M\,dz \le C'_M\varepsilon^{2M-2n}. \tag{IP19} No trace-norm estimate is inferred from the final inequality. Repeated W31 products show the same conclusion when finitely many uniformly bounded S(1,gε)S(1,g_\varepsilon) factors stand before or after this S(hεM,gε)S(h_\varepsilon^M,g_\varepsilon) remainder, in any fixed displayed order.

A separate proof of the uniform trace-norm bound

The uniform symbol hypothesis gives a stronger conclusion than (IP19). If the original matrix symbols qε(z)q_\varepsilon(z) are uniformly bounded in S(hεM,gε)S(h_\varepsilon^M,g_\varepsilon), with an integer M≥n+1M\geq n+1, then ∥qεw∥1≤Cε2M−2n,0<ε≤1.(IP22) \|q_\varepsilon^w\|_1\leq C\varepsilon^{2M-2n}, \qquad 0<\varepsilon\leq1. \tag{IP22} The constant uses finitely many of those original symbol seminorms. This statement applies to the genuine remainders and every displayed finite ordered product containing one such factor. It does not supply a uniform estimate for the full errors (IP8), since their cutoff term fails the required weight hypothesis.

Put k=n+1k=n+1 and retain the full auxiliary oscillator Hε=I+ε2∑j=1n(xj2+Dj2),Hεhγ=(1+nε2+2ε2|γ|)hγ,∥(Hε−k/2⊗Iν)∥22=ν∑r=0∞(r+n−1n−1)(1+nε2+2ε2r)−k≤Cn,kνε−2n.(IP23) \begin{gathered} H_\varepsilon=I+\varepsilon^2\sum_{j=1}^n(x_j^2+D_j^2),\\ H_\varepsilon h_\gamma =(1+n\varepsilon^2+2\varepsilon^2|\gamma|)h_\gamma,\\ \|(H_\varepsilon^{-k/2}\otimes I_\nu)\|_2^2 =\nu\sum_{r=0}^\infty\binom{r+n-1}{n-1} (1+n\varepsilon^2+2\varepsilon^2r)^{-k} \leq C_{n,k}\nu\varepsilon^{-2n}. \end{gathered} \tag{IP23} The basis and every constant in its eigenvalue come from CT3–CT5, without shifting away II or nε2n\varepsilon^2. More precisely, Hε=(1−ε2)I+ε2HH_\varepsilon=(1-\varepsilon^2)I+\varepsilon^2H on the original D(H)D(H); for each fixed positive ε\varepsilon, its diagonal eigenvalues are comparable to those of HH, so this is exactly its selfadjoint domain. Write λε,γ=1+nε2+2ε2|γ|\lambda_{\varepsilon,\gamma}=1+n\varepsilon^2+2\varepsilon^2|\gamma|. Every real spectral power used here has the explicit domain and action D(Hεs)={f:∑γλε,γ2s|fγ|2<∞},Hεsf=∑γλε,γsfγhγ,s>0,Hε−s:Hν→D(Hεs),HεsHε−sf=f,Hε−sHεsf=f(f∈D(Hεs)).(IP23a) \begin{gathered} D(H_\varepsilon^s)=\left\{f:\sum_\gamma \lambda_{\varepsilon,\gamma}^{2s}|f_\gamma|^2<\infty\right\}, \qquad H_\varepsilon^sf=\sum_\gamma \lambda_{\varepsilon,\gamma}^sf_\gamma h_\gamma, \quad s>0,\\ H_\varepsilon^{-s}:H_\nu\longrightarrow D(H_\varepsilon^s), \qquad H_\varepsilon^sH_\varepsilon^{-s}f=f,\quad H_\varepsilon^{-s}H_\varepsilon^sf=f\ (f\in D(H_\varepsilon^s)). \end{gathered} \tag{IP23a} For vector-valued functions, each coefficient is in ℂν\mathbb C^\nu and the sums use its full squared norm. The diagonal coefficient calculation proves both identities and closedness of the graph directly. To prove the bound (IP23), the terms with r<ε−2r<\varepsilon^{-2} have total multiplicity at most Cnε−2nC_n\varepsilon^{-2n}. On each dyadic block 2jε−2≤r<2j+1ε−22^j\varepsilon^{-2}\leq r<2^{j+1}\varepsilon^{-2}, the multiplicity sum is at most Cn2jnε−2nC_n2^{jn}\varepsilon^{-2n}, and the displayed inverse eigenvalue power is at most Ck2−jkC_k2^{-jk}. Their full sum is bounded by Cε−2n∑j≥02j(n−k)C\varepsilon^{-2n}\sum_{j\geq0}2^{j(n-k)}, which converges because k=n+1k=n+1. Endpoint rounding changes only the finite dimensional constants.

The graph estimate needed for the other factor is ∥Hεk/2f∥2≤Cn,k∑|α|+|β|≤kε|α|+|β|∥xαDβf∥2.(IP24) \|H_\varepsilon^{k/2}f\|_2 \leq C_{n,k}\sum_{|\alpha|+|\beta|\leq k} \varepsilon^{|\alpha|+|\beta|} \|x^\alpha D^\beta f\|_2. \tag{IP24} First, for each Hermite coefficient, the full eigenvalue satisfies (1+nε2+2ε2|γ|)k≤Cn,k(1+ε2k∑j(γj+k)!/γj!)(1+n\varepsilon^2+2\varepsilon^2|\gamma|)^k \leq C_{n,k}(1+\varepsilon^{2k}\sum_j(\gamma_j+k)!/\gamma_j!). This follows from (1+s)k≤2k−1(1+sk)(1+s)^k\leq2^{k-1}(1+s^k), 0<ε≤10<\varepsilon\leq1, and a coordinate with γj≥|γ|/n\gamma_j\geq|\gamma|/n, exactly as in CT7–CT8. Summing gives a bound by ∥f∥22+ε2k∑j∥(cj*)kf∥22\|f\|_2^2+\varepsilon^{2k}\sum_j\|(c_j^*)^kf\|_2^2. Expand each creator power with every lower commutator term retained, as in CT9. A resulting monomial has total degree d≤kd\leq k, and its actual coefficient εk\varepsilon^k obeys εk≤εd\varepsilon^k\leq\varepsilon^d. The triangle inequality proves (IP24) on Schwartz vectors. The diagonal spectral graph norm is closed, so convergence of all the weighted derivative terms passes it to the corresponding graph domain.

We now prove the estimate for the original qεwq_\varepsilon^w. The coordinate pullback q̂ε(w)=ε−2Mqε(w/ε)\widehat q_\varepsilon(w)=\varepsilon^{-2M}q_\varepsilon(w/\varepsilon) is used only to evaluate its symbol norms and integrals. Its exact derivative comparison is ∂wσq̂ε(w)=ε−2M−|σ|(∂zσqε)(w/ε),|∂wσq̂ε(w)|HS≤Cσh(w)M+|σ|/2.(IP25) \partial_w^\sigma\widehat q_\varepsilon(w) =\varepsilon^{-2M-|\sigma|} (\partial_z^\sigma q_\varepsilon)(w/\varepsilon), \qquad |\partial_w^\sigma\widehat q_\varepsilon(w)|_{\mathrm{HS}} \leq C_\sigma h(w)^{M+|\sigma|/2}. \tag{IP25} Thus this family has uniform S(hM,g)S(h^M,g) bounds, with the original factor qε(z)=ε2Mq̂ε(εz)q_\varepsilon(z)=\varepsilon^{2M}\widehat q_\varepsilon(\varepsilon z) explicit. Applying εxj\varepsilon x_j or εDj\varepsilon D_j on the output uses the full left-action identities CT10, including their coefficients iε/2i\varepsilon/2 and −iε/2-i\varepsilon/2. After d≤kd\leq k such ordered actions, every symbol term is Cεdzρ∂zσqε(z)=Cε2M+d−|ρ|+|σ|wρ∂wσq̂ε(w)|w=εz,|ρ|+|σ|≤d≤k.(IP26) \begin{split} C\varepsilon^d z^\rho\partial_z^\sigma q_\varepsilon(z) &=C\varepsilon^{2M+d-|\rho|+|\sigma|} w^\rho\partial_w^\sigma\widehat q_\varepsilon(w) \big|_{w=\varepsilon z},\\ &\hspace{10mm}|\rho|+|\sigma|\leq d\leq k. \end{split} \tag{IP26} Derivatives hitting a polynomial weight produce lower-degree monomials with the same actual εd\varepsilon^d; none is omitted. The extra exponent d−|ρ|+|σ|d-|\rho|+|\sigma| is nonnegative. In the original 2n2n phase dimensions, changing w=εzw=\varepsilon z contributes ε−2n\varepsilon^{-2n} to the squared integral. By (IP25), each squared ww-norm is bounded by the integral of a constant times ⟨w⟩2|ρ|−4M−2|σ|\langle w\rangle^{2|\rho|-4M-2|\sigma|}. Its worst exponent is 2k−4M≤−2(n+1)<−2n2k-4M\leq-2(n+1)<-2n, so every term is integrable uniformly.

The exact Hilbert–Schmidt formula CI5 and (IP24) now give ∥(Hεk/2⊗Iν)qεw∥2≤C(2π)−n/2ε2M−n.(IP27) \|(H_\varepsilon^{k/2}\otimes I_\nu)q_\varepsilon^w\|_2 \leq C(2\pi)^{-n/2}\varepsilon^{2M-n}. \tag{IP27} One can first use a compact smooth phase cutoff, where every output is Schwartz by CT11a. For fixed ε\varepsilon, the weighted derivative estimates above give convergence of the cutoff symbols and every output-weighted operator in Hilbert–Schmidt norm, by CI5 and the finite CT10 action. The closed spectral graph of Hεk/2H_\varepsilon^{k/2} then identifies the limiting second factor in (IP27) and shows that every output of qεwq_\varepsilon^w lies in its domain. No multiplication of an unverified unbounded operator is used.

The actual factorization on the original Hilbert space is consequently qεw=(Hε−k/2⊗Iν)((Hεk/2⊗Iν)qεw).(IP28) q_\varepsilon^w=(H_\varepsilon^{-k/2}\otimes I_\nu) ((H_\varepsilon^{k/2}\otimes I_\nu)q_\varepsilon^w). \tag{IP28} Use the exact inverse-domain identity in this product, then multiply the two Hilbert–Schmidt bounds (IP23),(IP27) with the trace-ideal inequality. Their exponents are −n-n and 2M−n2M-n, so the result is (IP22), with all original Fourier and fiber factors in the stated constant. W31 puts every finite ordered product containing the genuine remainder uniformly in the same S(hεM,gε)S(h_\varepsilon^M,g_\varepsilon) class; applying this proved estimate gives its trace-norm bound in the same order. In particular the telescoping differences used in (IP21) and the full analytic remainder (DE7) tend to zero in trace norm, as well as in trace.

The complete scaled-oscillator factorization and its original phase-space factors

Equations (IP23)–(IP28) prove both typed factors and every exponent in the diagram. The original symbol, phase coordinates, weight and cutoff remain present in their separate formulas. This is a proved strengthening of the genuine-remainder estimate; it does not assert a trace-norm bound for the full errors.

To connect this genuine remainder estimate to the full errors, retain the cutoff and define the finite symbols sj,ε(M)=(1−ψε)Iν−∑k=1M−1Ck(fj,ε,gj,ε),(f1,ε,g1,ε)=(bε,aε),(f2,ε,g2,ε)=(aε,bε).(IP20) s_{j,\varepsilon}^{(M)} =(1-\psi_\varepsilon)I_\nu -\sum_{k=1}^{M-1}C_k(f_{j,\varepsilon},g_{j,\varepsilon}), \quad (f_{1,\varepsilon},g_{1,\varepsilon})=(b_\varepsilon,a_\varepsilon), \quad (f_{2,\varepsilon},g_{2,\varepsilon})=(a_\varepsilon,b_\varepsilon). \tag{IP20} Equations (IP8),(IP18) give rj,ε=sj,ε(M)−qj,ε(M)r_{j,\varepsilon} =s_{j,\varepsilon}^{(M)}-q_{j,\varepsilon}^{(M)}. For any ordered associative product, the noncommutative telescoping identity u#N−v#N=∑ℓ=0N−1u#ℓ#(u−v)#v#(N−1−ℓ)u^{\# N}-v^{\# N} =\sum_{\ell=0}^{N-1}u^{\#\ell}\#(u-v)\#v^{\#(N-1-\ell)} follows by expanding adjacent differences; it does not interchange factors. Apply it with u=rj,εu=r_{j,\varepsilon} and v=sj,ε(M)v=s_{j,\varepsilon}^{(M)}. Both are uniformly bounded in S(1,gε)S(1,g_\varepsilon), whereas u−v=−qj,ε(M)u-v=-q_{j,\varepsilon}^{(M)} is uniformly in S(hεM,gε)S(h_\varepsilon^M,g_\varepsilon). W31, the fixed-ε\varepsilon trace criterion and (IP17) therefore prove Tr⁡Rj,εN−Tr⁡((sj,ε(M))#N)w=O(ε2M−2n),ind⁡aw=Tr⁡((s1,ε(M))#N)w−Tr⁡((s2,ε(M))#N)w+O(ε2M−2n),N≥n+1.(IP21) \begin{split} &\operatorname{Tr}R_{j,\varepsilon}^{\,N} -\operatorname{Tr} \bigl((s_{j,\varepsilon}^{(M)})^{\# N}\bigr)^w =O(\varepsilon^{2M-2n}),\\ &\operatorname{ind}a^w =\operatorname{Tr} \bigl((s_{1,\varepsilon}^{(M)})^{\# N}\bigr)^w -\operatorname{Tr} \bigl((s_{2,\varepsilon}^{(M)})^{\# N}\bigr)^w +O(\varepsilon^{2M-2n}), \qquad N\ge n+1 . \end{split} \tag{IP21} Each finite-symbol power in (IP21) is trace class for fixed ε\varepsilon: every sj,ε(M)s_{j,\varepsilon}^{(M)} belongs to S(hε,gε)S(h_\varepsilon,g_\varepsilon) with fixed-ε\varepsilon constants, and the argument of (IP11)–(IP13) applies. The displayed errors are uniform in ε\varepsilon, because every telescoping summand contains the one genuinely uniform hεMh_\varepsilon^M factor. This proves a bounded remainder step while keeping the unbounded-support cutoff contribution in the finite main expression. It does not yet identify which of that expression’s differential degrees survive.

8. Differential coefficients in the two original orders

Keep n,ν≥1n,\nu\ge1, h(z)=(1+|z|2)−1h(z)=(1+|z|^2)^{-1}, g=h(|dx|2+|dξ|2)g=h(|dx|^2+|d\xi|^2), and the original square matrix symbol aa, cutoff ψ\psi, and b=ψa−1b=\psi a^{-1} from IP1–IP3. Denote the actual ordered Weyl differential coefficient of W29 by Cℓ(f,g)C_\ell(f,g), so that C0(f,g)=fgC_0(f,g)=fg, the first argument’s matrix derivatives stand to the left of the second argument’s, and CℓC_\ell has exactly ℓ\ell phase derivatives on each factor. Put (f1,g1)=(b,a),(f2,g2)=(a,b),cj,0=(1−ψ)Iν,cj,k=−Ck(fj,gj)(k≥1).(DE1) (f_1,g_1)=(b,a),\qquad(f_2,g_2)=(a,b),\qquad c_{j,0}=(1-\psi)I_\nu,\qquad c_{j,k}=-C_k(f_j,g_j)\quad(k\ge1). \tag{DE1} The scalar ψ\psi is one outside a compact set, so cj,0c_{j,0} is compactly supported. Since fj,gj∈S(1,g)f_j,g_j\in S(1,g), the product rule in the explicit differential formula gives cj,k∈S(hk,g)c_{j,k}\in S(h^k,g) for every k≥1k\ge1. No equality between Ck(b,a)C_k(b,a) and Ck(a,b)C_k(a,b) is assumed for matrix symbols.

For 0<ε≤10<\varepsilon\le1, write 𝒮εf(z)=f(εz)\mathscr S_\varepsilon f(z)=f(\varepsilon z). If f∈S(hi,g)f\in S(h^i,g), the chain rule proves ε2i𝒮εf∈S(hεi,gε)with seminorms uniform in ε,hε(z)=ε2h(εz).(DE2) \varepsilon^{2i}\mathscr S_\varepsilon f \in S(h_\varepsilon^i,g_\varepsilon) \quad\text{with seminorms uniform in }\varepsilon, \qquad h_\varepsilon(z)=\varepsilon^2h(\varepsilon z). \tag{DE2} Indeed an order-qq derivative on the left is bounded by Cε2i+qh(εz)i+q/2=Chε(z)i+q/2C\varepsilon^{2i+q}h(\varepsilon z)^{i+q/2} =C h_\varepsilon(z)^{i+q/2}. Since CℓC_\ell differentiates both of its inputs exactly ℓ\ell times in total, the scaled coefficients obey the exact identity Cℓ(ε2i𝒮εf,ε2j𝒮εg)=ε2(i+j+ℓ)𝒮εCℓ(f,g).(DE3) C_\ell(\varepsilon^{2i}\mathscr S_\varepsilon f, \varepsilon^{2j}\mathscr S_\varepsilon g) =\varepsilon^{2(i+j+\ell)} \mathscr S_\varepsilon C_\ell(f,g). \tag{DE3} Thus powers of ε2\varepsilon^2 count pairs of actual phase derivatives; no factor is absorbed into a completed or normalized symbol.

9. A finite noncommutative star polynomial

Use a formal variable tt only to organize the finite calculation. For coefficient sequences (ui),(vj)(u_i),(v_j), define their degree-kk ordered Weyl polynomial coefficient by (u⋆tv)k=∑i,j,ℓ≥0i+j+ℓ=kCℓ(ui,vj).(DE4) (u\star_t v)_k =\sum_{\substack{i,j,\ell\ge0\\i+j+\ell=k}} C_\ell(u_i,v_j). \tag{DE4} This is a finite sum for every kk, retains the order of matrix multiplication, and is associative through each finite degree because the actual Weyl product is associative and W32 supplies a remainder in the next metric weight. Define Fj,k(N)F_{j,k}^{(N)} as the degree-kk coefficient obtained by applying (DE4) to NN copies of the sequence cj,0,cj,1,…c_{j,0},c_{j,1},\ldots, with left-associative parentheses. Associativity shows that the chosen parentheses do not change the coefficient. Every summand is an ordered differential word in exactly NN of the cj,ic_{j,i}, with total scale degree kk.

We now prove the analytic statement behind this finite algebra. For any fixed integers N≥1N\ge1 and M≥1M\ge1, rj,ε#N=∑k=0M−1ε2k𝒮εFj,k(N)+Ej,ε(N,M),Ej,ε(N,M)∈S(hεM,gε)uniformly in ε.(DE5) r_{j,\varepsilon}^{\# N} =\sum_{k=0}^{M-1} \varepsilon^{2k} \mathscr S_\varepsilon F_{j,k}^{(N)} +E_{j,\varepsilon}^{(N,M)},\qquad E_{j,\varepsilon}^{(N,M)} \in S(h_\varepsilon^M,g_\varepsilon) \ \text{uniformly in }\varepsilon . \tag{DE5} For N=1N=1, IP18 and (DE1) give this identity, including the sign of the two product remainders. Suppose it holds for a given number of factors. Multiply its left and right sides by rj,εr_{j,\varepsilon} in that original order. IP8–IP9 make rj,εr_{j,\varepsilon} uniformly bounded in S(1,gε)S(1,g_\varepsilon), so W31 puts the previous E(N,M)#rE^{(N,M)}\#r uniformly in S(hεM,gε)S(h_\varepsilon^M,g_\varepsilon). Expand the new factor to degree M−1M-1 by the N=1N=1 case; its remainder, multiplied by the finite previous sum, lies in the same weight. For a pair of retained terms of intrinsic degrees i,ji,j with i+j<Mi+j<M, W32 with positive truncation M−i−jM-i-j gives all CℓC_\ell with i+j+ℓ<Mi+j+\ell<M, and puts the rest in S(hεM,gε)S(h_\varepsilon^M,g_\varepsilon). If i+j≥Mi+j\geq M, apply W31 to the entire ordered product, placing it uniformly in S(hεi+j,gε)S(h_\varepsilon^{i+j},g_\varepsilon). Its exact inclusion bound in the required remainder class is pl(v;hεM,gε)≤ε2(i+j−M)pl(v;hεi+j,gε),v∈S(hεi+j,gε),i+j≥M,(DE5a) p_l(v;h_\varepsilon^M,g_\varepsilon) \leq\varepsilon^{2(i+j-M)} p_l(v;h_\varepsilon^{i+j},g_\varepsilon), \quad v\in S(h_\varepsilon^{i+j},g_\varepsilon),\quad i+j\geq M, \tag{DE5a} because hεi+j−M≤ε2(i+j−M)≤1h_\varepsilon^{i+j-M}\leq\varepsilon^{2(i+j-M)}\leq1. Here plp_l is the original directional-derivative seminorm; each of its quotients has exactly this weight ratio. This case uses no negative truncation order and retains the entire high-degree product in the proved remainder. The uniform constants follow from IP5–IP6 and W33; (DE3) identifies the retained term exactly. This proves (DE5) by induction, without treating a formal series as convergent.

10. Compact support of every coefficient through degree nn

Fix the trace-power exponent N=n+1N=n+1 and truncate (DE5) at M=n+1M=n+1. For every 0≤k≤n0\le k\le n, each summand of Fj,k(n+1)F_{j,k}^{(n+1)} contains n+1n+1 original cc-factors. If all those factors had positive intrinsic degree, their degrees would sum to at least n+1>kn+1>k, even before adding the nonnegative Weyl differential degrees. Therefore at least one factor is cj,0=(1−ψ)Iνc_{j,0}=(1-\psi)I_\nu. The CℓC_\ell are local differential expressions, so every derivative of that factor remains supported in the same fixed compact set. Consequently Fj,k(n+1)∈Cc∞(ℝ2n;End⁡ℂν),0≤k≤n.(DE6) F_{j,k}^{(n+1)} \in C_c^\infty(\mathbb R^{2n}; \operatorname{End}\mathbb C^\nu), \qquad 0\le k\le n . \tag{DE6} This is the exact support reason the separate coefficient integrals exist. A decay claim for cj,kc_{j,k} alone at k≤nk\le n would not give the same conclusion.

11. Integrable trace remainder and the exact expansion

The remainder in (DE5) with N=M=n+1N=M=n+1 is uniform in S(hεn+1,gε)S(h_\varepsilon^{n+1},g_\varepsilon). For each fixed ε\varepsilon, the full four-multiindex estimate IP12 proves that its Weyl operator is trace class. The exact matrix Weyl trace formula and IP17 then give |Tr(Ej,ε(n+1,n+1))w|≤(2π)−nνC∫ℝ2nhεn+1(z)dz=O(ε2).(DE7) \left|\operatorname{Tr} (E_{j,\varepsilon}^{(n+1,n+1)})^w\right| \le (2\pi)^{-n}\nu C \int_{\mathbb R^{2n}}h_\varepsilon^{n+1}(z)\,dz =O(\varepsilon^2). \tag{DE7} For 0≤k≤n0\le k\le n, (DE6) makes ε2k𝒮εFj,k(n+1)\varepsilon^{2k}\mathscr S_\varepsilon F_{j,k}^{(n+1)} a compactly supported smooth phase-space symbol for fixed ε\varepsilon. It is trace class by CT1 and has the exact CI12 trace Tr⁡(ε2k𝒮εFj,k(n+1))w=(2π)−nε2k−2n∫ℝ2ntr⁡ℂνFj,k(n+1)(w)dw.(DE8) \operatorname{Tr} \bigl(\varepsilon^{2k} \mathscr S_\varepsilon F_{j,k}^{(n+1)}\bigr)^w =(2\pi)^{-n}\varepsilon^{2k-2n} \int_{\mathbb R^{2n}} \operatorname{tr}_{\mathbb C^\nu} F_{j,k}^{(n+1)}(w)\,dw . \tag{DE8} The Jacobian is ε−2n\varepsilon^{-2n}, not ε−n\varepsilon^{-n}, because both xx and ξ\xi are scaled.

Set the finite coefficients Ak=(2π)−n∫ℝ2ntr⁡ℂν(F1,k(n+1)−F2,k(n+1))(w)dw,0≤k≤n.(DE9) A_k=(2\pi)^{-n} \int_{\mathbb R^{2n}} \operatorname{tr}_{\mathbb C^\nu} \bigl(F_{1,k}^{(n+1)} -F_{2,k}^{(n+1)}\bigr)(w)\,dw, \quad 0\le k\le n . \tag{DE9} They are absolutely defined by (DE6). Insert (DE5)–(DE8) into the already proved IP14 index formula and use IP16 to identify the left side with the original operator index. The result is the full finite asymptotic equality ind⁡aw=∑k=0nAkε2k−2n+O(ε2),ε↓0.(DE10) \operatorname{ind}a^w =\sum_{k=0}^{n}A_k\varepsilon^{2k-2n} +O(\varepsilon^2), \qquad \varepsilon\downarrow0. \tag{DE10} No trace of an individual nonintegrable coefficient was written: each coefficient in (DE9) has compact support, and the remainder has the uniform trace bound (DE7).

12. The exact cancellation and surviving degree

The left side of (DE10) is one fixed integer for all positive ε\varepsilon. If some Ak≠0A_k\ne0 with k<nk<n, choose the smallest such kk. After multiplication by ε2n−2k\varepsilon^{2n-2k}, its term tends to AkA_k, all larger-kk terms and the O(ε2)O(\varepsilon^2) remainder tend to zero, while the constant left side tends to zero. This contradiction proves Ak=0A_k=0. Repeating for each lower kk, then letting ε↓0\varepsilon\downarrow0 in (DE10), gives A0=A1=⋯=An−1=0,ind⁡aw=An.(DE11) A_0=A_1=\cdots=A_{n-1}=0,\qquad \operatorname{ind}a^w=A_n . \tag{DE11} A degree kk term in (DE4) has kk Weyl derivative pairs, hence exactly 2k2k total phase differentiations on the original ordered factor expressions, including derivatives of b=ψa−1b=\psi a^{-1} and of the retained cutoff. Thus (DE11) says precisely that only total differential degree 2n2n contributes after taking the difference of the two full matrix traces. It does not assert that the lower-degree trace of either error power vanishes separately. The cancellation is proved by the exact Fredholm index constancy together with the uniform analytic remainder, rather than by an unsupported formal cancellation.

13. Two examples with complete solutions

Example 1: constant invertible matrix. Take a(z)=A0∈GL⁡ν(ℂ)a(z)=A_0\in\operatorname{GL}_{\nu}(\mathbb C) at every zz, and choose ψ=1\psi=1. The actual cutoff inverse is b=A0−1b=A_0^{-1}. Every positive-order derivative of either factor vanishes, so both ordered Weyl products are the identity. Thus, at every scale, r1,ε=r2,ε=0,F1,k(n+1)=F2,k(n+1)=0(0≤k≤n),ind⁡A0w=An=0.(SD1) r_{1,\varepsilon}=r_{2,\varepsilon}=0,\qquad F_{1,k}^{(n+1)}=F_{2,k}^{(n+1)}=0 \quad(0\le k\le n),\qquad \operatorname{ind}A_0^w=A_n=0 . \tag{SD1} This uses the full n+1n+1-power formula, although the errors happen to vanish already before taking powers.

Example 2: why the full error has no uniform hεh_\varepsilon bound. Keep a(z)=Iνa(z)=I_\nu, but choose a smooth scalar ψ\psi that is zero in a neighborhood of z=0z=0 and one outside a compact set. The hypotheses in (IP2) hold because aa is invertible everywhere. Now b=ψIνb=\psi I_\nu. Weyl multiplication by the constant identity is exact, so d1,ε=d2,ε=0d_{1,\varepsilon}=d_{2,\varepsilon}=0 and r1,ε(0)=r2,ε(0)=Iν,hε(0)=ε2.(SD2) r_{1,\varepsilon}(0)=r_{2,\varepsilon}(0)=I_\nu, \qquad h_\varepsilon(0)=\varepsilon^2 . \tag{SD2} If the full rj,εr_{j,\varepsilon} had a uniform S(hε,gε)S(h_\varepsilon,g_\varepsilon) zeroth seminorm, (SD2) would require 1≤Cε21\le C\varepsilon^2 for every positive ε\varepsilon, which is impossible. For each fixed ε\varepsilon, (IP10) remains true. This distinction is why the finite expansion retains (1−ψε)Iν(1-\psi_\varepsilon)I_\nu separately.

14. Current coefficient and next calculation

Here the complete ordered coefficient AnA_n is defined and proved to equal the analytic index. Identifying that coefficient with its noncommutative exterior boundary integral requires the later programme coefficient calculation; no such identification is used in this lesson.

15. Editorial supplement: the full product weight at this receiving map

The original mathematics above, including its integer truncations and (IP22), is retained. This supplement records the exact receiving map for the later product extension (WG1)–(WG4) in Weyl products, and proves a stronger trace estimate for the same original symbols.

For the original pair of identical metrics gε=hεeg_\varepsilon=h_\varepsilon e, the cross parameter of (W6) is exactly H=hεH=h_\varepsilon. Indeed the symplectic dual is hε−1eh_\varepsilon^{-1}e, so the ratio defining its squared cross parameter is hε2h_\varepsilon^2 on every nonzero direction. The actual product-space multiplier parameter remains hε/4h_\varepsilon/4, the diagonal metric remains 2gε2g_\varepsilon, and the product space still has dimension 4n4n. Therefore the full extra weight is wε(z)=(1+hε(z)/4)4n,1≤wε(z)≤(1+ε2/4)4n≤(5/4)4n.(RP1) w_\varepsilon(z)=(1+h_\varepsilon(z)/4)^{4n},\qquad 1\leq w_\varepsilon(z) \leq(1+\varepsilon^2/4)^{4n}\leq(5/4)^{4n}. \tag{RP1} No factor has been removed from the provider’s actual estimate. For any real s,ts,t, every target derivative order ll, and every nonnegative integer truncation NN, the specialization of (WG4) gives pl(RN(u,v);hεs+t+Nwε,gε)≤CN,l4−N2l/22Jp≤J(u;hεs,gε)p≤J(v;hεt,gε),pl(RN(u,v);hεs+t+N,gε)≤(1+ε2/4)4npl(RN(u,v);hεs+t+Nwε,gε).(RP2) \begin{aligned} &p_l\!\left(R_N(u,v); h_\varepsilon^{s+t+N}w_\varepsilon,g_\varepsilon\right) \leq C_{N,l}\,4^{-N}2^{l/2}2^J p_{\leq J}(u;h_\varepsilon^s,g_\varepsilon) p_{\leq J}(v;h_\varepsilon^t,g_\varepsilon),\\ &p_l\!\left(R_N(u,v); h_\varepsilon^{s+t+N},g_\varepsilon\right) \leq(1+\varepsilon^2/4)^{4n} p_l\!\left(R_N(u,v); h_\varepsilon^{s+t+N}w_\varepsilon,g_\varepsilon\right). \end{aligned} \tag{RP2} Here R0(u,v)=u#vR_0(u,v)=u\#v; RNR_N and its finite coefficients retain the provider’s order. The finite input integer JJ, its Gaussian iteration condition (WG3), and CN,lC_{N,l} use the original structural constants. Those constants are uniform here by (IP5)–(IP6) and the power-weight comparisons. For the second inequality, each original directional-derivative quotient is multiplied by precisely wε(z)w_\varepsilon(z); its supremum is bounded by the displayed upper bound. The definition (B4) differentiates the symbol, with the weight in the denominator, so it does not require differentiating or replacing that weight.

The derivative comparison is also available explicitly. Write Lε(z)=ε−2+|z|2L_\varepsilon(z)=\varepsilon^{-2}+|z|^2. For every real ss and every coordinate multiindex σ\sigma, repeated differentiation of the original hεs=Lε−sh_\varepsilon^s=L_\varepsilon^{-s} gives a finite sum ∂zσhε(z)s=∑p,ρAσ,p,ρ(s)zρLε(z)−s−p,2p−|ρ|=|σ|.(RP3) \partial_z^\sigma h_\varepsilon(z)^s =\sum_{p,\rho} A_{\sigma,p,\rho}(s) z^\rho L_\varepsilon(z)^{-s-p},\qquad 2p-|\rho|=|\sigma|. \tag{RP3} For σ=0\sigma=0, the sole coefficient is A0,0,0=1A_{0,0,0}=1. The exact recursion is ∂zj(zρLε−s−p)=ρjzρ−ejLε−s−p−2(s+p)zρ+ejLε−s−p−1.(RP4) \partial_{z_j}\!\left(z^\rho L_\varepsilon^{-s-p}\right) =\rho_j z^{\rho-e_j}L_\varepsilon^{-s-p} -2(s+p)z^{\rho+e_j}L_\varepsilon^{-s-p-1}. \tag{RP4} The first term is zero when ρj=0\rho_j=0, and all remaining terms, coefficients and signs are retained. This proves (RP3) by induction, because both new terms have 2p−|ρ|=|σ|+12p-|\rho|=|\sigma|+1 with their respective new indices. Since |zj|≤Lε1/2|z_j|\leq L_\varepsilon^{1/2}, every summand has absolute value at most |Aσ,p,ρ(s)|hεs+|σ|/2|A_{\sigma,p,\rho}(s)|h_\varepsilon^{s+|\sigma|/2}. Thus |∂zσhεs|≤Cs,σhεs+|σ|/2,∂zσwε=∑r=04n(4nr)4−r∂zσhεr,|∂zσwε|≤Cn,σhε|σ|/2.(RP5) |\partial_z^\sigma h_\varepsilon^s| \leq C_{s,\sigma}h_\varepsilon^{s+|\sigma|/2},\qquad \partial_z^\sigma w_\varepsilon =\sum_{r=0}^{4n}\binom{4n}{r}4^{-r} \partial_z^\sigma h_\varepsilon^r, \qquad |\partial_z^\sigma w_\varepsilon| \leq C_{n,\sigma}h_\varepsilon^{|\sigma|/2}. \tag{RP5} In the last bound, use hε≤1h_\varepsilon\leq1; for a positive derivative order the r=0r=0 term is exactly zero. The finite polynomial expansion keeps every 4−r4^{-r} and binomial coefficient. In particular the full weight and all its derivatives have uniform S(1,gε)S(1,g_\varepsilon) estimates, rather than being declared absent.

Applying (RP2) in a finite displayed product gives a finite chain of these explicit factors. In particular, if one input is uniform in S(hεM,gε)S(h_\varepsilon^M,g_\varepsilon) and the remaining finitely many inputs are uniform in S(1,gε)S(1,g_\varepsilon), every fixed ordered parenthesization lands uniformly in S(hεM,gε)S(h_\varepsilon^M,g_\varepsilon). At each application the bound includes (1+ε2/4)4n2l/22J(1+\varepsilon^2/4)^{4n}2^{l/2}2^J and the actual source seminorms; the list of derivative orders and input integers can change in a finite chain. This is an explicit inclusion of the provider’s larger target into the original receiving class. For a retained product of intrinsic degree i+j≥Mi+j\geq M, the original (DE5a) bound is still exactly hεi+j−M≤ε2(i+j−M)h_\varepsilon^{i+j-M}\leq\varepsilon^{2(i+j-M)}, after the factor in (RP2). The entire high-degree product remains in the remainder.

The operator product used here has its actual domain. Every S(1,gε)S(1,g_\varepsilon) symbol has bounded coordinate derivatives, since hε≤1h_\varepsilon\leq1, and hence belongs to the classical class S0,00S^0_{0,0} without changing its coordinates or values. The same holds for S(hεM,gε)S(h_\varepsilon^M,g_\varepsilon) when M>0M>0. The exact matrix maps (WO6)–(WO8) give (u#v)w=uwvw(u\#v)^w=u^wv^w on 𝒮(ℝn;ℂν)\mathcal S(\mathbb R^n;\mathbb C^\nu), with each intermediate fiber index summed in its original order. B26 and (RP2) bound both sides on Hν=L2(ℝn;ℂν)H_\nu=L^2(\mathbb R^n;\mathbb C^\nu). Density extends that equality uniquely to all of HνH_\nu. Thus no composition of two unspecified maps to distributions is used.

16. Editorial supplement: every real remainder exponent above the dimension

Theorem. Retain n,ν≥1n,\nu\geq1, the original coordinates z=(x,ξ)z=(x,\xi), the metric and weight (IP4), and the original Weyl quantization (CI4). Let MM be any real number with M>nM>n. Suppose qε(z)∈End⁡ℂνq_\varepsilon(z)\in\operatorname{End}\mathbb C^\nu is smooth and uniform in S(hεM,gε)S(h_\varepsilon^M,g_\varepsilon) for 0<ε≤10<\varepsilon\leq1. Then ∥qεw∥1≤Cn,M,νQε2M−2n,Q=max|σ|≤n+1sup0<ε≤1z∈ℝ2n∥∂zσqε(z)∥HShε(z)M+|σ|/2.(RA1) \|q_\varepsilon^w\|_1 \leq C_{n,M,\nu}\,Q\, \varepsilon^{2M-2n},\qquad Q=\max_{|\sigma|\leq n+1} \sup_{\substack{0<\varepsilon\leq1\\z\in\mathbb R^{2n}}} {\|\partial_z^\sigma q_\varepsilon(z)\|_{\mathrm{HS}} \over h_\varepsilon(z)^{M+|\sigma|/2}}. \tag{RA1} The finite QQ is controlled by finitely many original directional symbol seminorms, with the original finite-dimensional norm comparison. Every fixed finite ordered product described after (RP5) obeys the same exponent, with its proved finite product-seminorm constant. This extends (IP22) while retaining the integer truncation statement (IP18).

Proof of the annular symbol estimates. Fix a smooth radial χ\chi on ℝ2n\mathbb R^{2n}, with 0≤χ≤10\leq\chi\leq1, equal to one for |w|≤1|w|\leq1 and zero for |w|≥2|w|\geq2, and nonincreasing in |w||w|. Such a cutoff follows, for example, from η(t)=0\eta(t)=0 for t≤0t\leq0 and η(t)=e−1/t\eta(t)=e^{-1/t} for t>0t>0, using χ(w)=η(4−|w|2)/(η(4−|w|2)+η(|w|2−1))\chi(w)=\eta(4-|w|^2)/(\eta(4-|w|^2)+\eta(|w|^2-1)); the denominator is positive at every ww. This formula retains both endpoint radii. Put ϕ0(w)=χ(w),ϕj(w)=χ(w/2j)−χ(w/2j−1)(j≥1),Rj=2j(j≥0),δj=ε/Rj,qε,j(z)=qε(z)ϕj(εz),∑j=0Jqε,j(z)=qε(z)χ(εz/2J).(RA2) \begin{gathered} \phi_0(w)=\chi(w),\qquad \phi_j(w)=\chi(w/2^j)-\chi(w/2^{j-1})\quad(j\geq1),\\ R_j=2^j\quad(j\geq0),\qquad \delta_j=\varepsilon/R_j,\qquad q_{\varepsilon,j}(z)=q_\varepsilon(z)\phi_j(\varepsilon z),\\ \sum_{j=0}^Jq_{\varepsilon,j}(z) =q_\varepsilon(z)\chi(\varepsilon z/2^J). \end{gathered} \tag{RA2} This is an exact finite telescoping identity for the original symbol, not a formal replacement. For j≥1j\geq1, the support satisfies Rj/2≤|εz|≤2RjR_j/2\leq|\varepsilon z|\leq2R_j. Define the comparison functions, keeping their factors: Aj=ε2MRj−2M,uε,0(v)=ε−2Mqε(v/ε)χ(v),uε,j(v)=ε−2MRj2Mqε(Rjv/ε)[χ(v)−χ(2v)](j≥1),qε,j(z)=Ajuε,j(δjz).(RA3) \begin{aligned} &A_j=\varepsilon^{2M}R_j^{-2M},\\ &u_{\varepsilon,0}(v) =\varepsilon^{-2M}q_\varepsilon(v/\varepsilon)\chi(v),\\ &u_{\varepsilon,j}(v) =\varepsilon^{-2M}R_j^{2M} q_\varepsilon(R_jv/\varepsilon) [\chi(v)-\chi(2v)]\quad(j\geq1),\\ &q_{\varepsilon,j}(z)=A_j u_{\varepsilon,j}(\delta_j z). \end{aligned} \tag{RA3} The comparison functions do not replace the original Weyl operator. They evaluate its derivatives and phase-space integrals. Their supports lie in |v|≤2|v|\leq2, and for j≥1j\geq1 in 1/2≤|v|≤21/2\leq|v|\leq2. For an original derivative τ\tau, the chain rule and the symbol hypothesis give ∥∂vτ[ε−2MRj2Mqε(Rjv/ε)]∥HS≤QRj2M+|τ|(1+Rj2|v|2)M+|τ|/2≤Q22M+|τ|(1/2≤|v|≤2,j≥1).(RA4) \begin{aligned} &\left\|\partial_v^\tau [\varepsilon^{-2M}R_j^{2M} q_\varepsilon(R_jv/\varepsilon)]\right\|_{\mathrm{HS}}\\ &\quad\leq Q\, {R_j^{2M+|\tau|}\over (1+R_j^2|v|^2)^{M+|\tau|/2}} \leq Q\,2^{2M+|\tau|} \quad(1/2\leq|v|\leq2, j\geq1). \end{aligned} \tag{RA4} Every ε\varepsilon exponent cancels exactly in this comparison: the original derivative contributes ε2M+|τ|\varepsilon^{2M+|\tau|}, and the displayed pullback contributes ε−2M−|τ|\varepsilon^{-2M-|\tau|}. The last inequality uses M+|τ|/2>0M+|\tau|/2>0 and 1+Rj2|v|2≥Rj2|v|21+R_j^2|v|^2\geq R_j^2|v|^2. For j=0j=0, the same first comparison with R0=1R_0=1 is bounded by Q(1+|v|2)−M−|τ|/2≤QQ(1+|v|^2)^{-M-|\tau|/2}\leq Q. The full multiindex product rule now bounds every ∂σuε,j\partial^\sigma u_{\varepsilon,j}, |σ|≤k=n+1|\sigma|\leq k=n+1, uniformly by Q∑τ≤σ(στ)22M+|τ|∥∂σ−τ(χ−χ(2⋅))∥∞Q\sum_{\tau\leq\sigma}\binom\sigma\tau 2^{2M+|\tau|}\|\partial^{\sigma-\tau}(\chi-\chi(2\cdot))\|_\infty for j≥1j\geq1, and by its corresponding χ\chi expression for j=0j=0. Derivatives of χ(2v)\chi(2v) retain the factor 2|σ−τ|2^{|\sigma-\tau|}. Smoothness makes all these finite constants valid at the support boundaries as well. Thus every compact-support comparison norm ∥vρ∂σuε,j∥Lv2;HS\|v^\rho\partial^\sigma u_{\varepsilon,j}\|_{L^2_v;\mathrm{HS}} with |ρ|+|σ|≤k|\rho|+|\sigma|\leq k is bounded by Cn,M,χQC_{n,M,\chi}Q, independently of ε,j\varepsilon,j.

Proof of the typed operator estimate on each annulus. Use the full auxiliary oscillator at the actual δj\delta_j: Hδj=I+δj2∑r=1n(xr2+Dr2),Hδjhγ=(1+nδj2+2δj2|γ|)hγ,Pj=Hδjk/2⊗Iν,∥Pj−1∥22=ν∑r=0∞(r+n−1n−1)(1+nδj2+2δj2r)−k≤Cn,kνδj−2n.(RA5) \begin{gathered} H_{\delta_j}=I+\delta_j^2\sum_{r=1}^n(x_r^2+D_r^2),\qquad H_{\delta_j}h_\gamma =(1+n\delta_j^2+2\delta_j^2|\gamma|)h_\gamma,\\ P_j=H_{\delta_j}^{k/2}\otimes I_\nu, \qquad \|P_j^{-1}\|_2^2 =\nu\sum_{r=0}^\infty\binom{r+n-1}{n-1} (1+n\delta_j^2+2\delta_j^2r)^{-k} \leq C_{n,k}\nu\delta_j^{-2n}. \end{gathered} \tag{RA5} This is (IP23) with 0<δj≤10<\delta_j\leq1. Its domain is the original coefficient domain (IP23a), with ε\varepsilon replaced by δj\delta_j; the identity, nδj2n\delta_j^2, shell multiplicity and every fiber are retained. Equations (IP24) and (CT10) likewise apply with that actual scale. In particular each left output action retains the original identities δjxraw=(δjxra+iδj2∂ξra)w,δjDraw=(δjξra−iδj2∂xra)w.(RA6) \delta_j x_r a^w =\left(\delta_j x_r a+ {i\delta_j\over2}\partial_{\xi_r}a\right)^w, \qquad \delta_j D_r a^w =\left(\delta_j\xi_r a- {i\delta_j\over2}\partial_{x_r}a\right)^w. \tag{RA6} After an ordered word of length d≤kd\leq k, induction using these two identities writes every term, with its actual finite coefficient CC, as Cδjdzρ∂zσqε,j(z)=CAjδjd−|ρ|+|σ|vρ∂vσuε,j(v)|v=δjz,|ρ|+|σ|≤d≤k.(RA7) \begin{aligned} C\delta_j^d z^\rho\partial_z^\sigma q_{\varepsilon,j}(z) &=C A_j\delta_j^{d-|\rho|+|\sigma|} v^\rho\partial_v^\sigma u_{\varepsilon,j}(v) \big|_{v=\delta_j z},\\ &|\rho|+|\sigma|\leq d\leq k. \end{aligned} \tag{RA7} Here zρz^\rho can contain both xx and ξ\xi coordinates. Each new multiplication increases |ρ||\rho| by one, each derivative of the symbol increases |σ||\sigma| by one, and differentiation of an earlier monomial retains its full coefficient and decreases |ρ||\rho| by one. All such terms are kept. Therefore d−|ρ|+|σ|≥0d-|\rho|+|\sigma|\geq0, and its actual δj\delta_j power is at most one. The original phase-space change v=δjzv=\delta_j z has Jacobian δj−2n\delta_j^{-2n}. Equations (RA4)–(RA7), the exact Weyl Hilbert–Schmidt identity (CI5), and (IP24) consequently prove ∥Pjqε,jw∥2≤Cn,M,χ,k(2π)−n/2QAjδj−n.(RA8) \|P_j q_{\varepsilon,j}^w\|_2 \leq C_{n,M,\chi,k}(2\pi)^{-n/2} Q A_j\delta_j^{-n}. \tag{RA8} To verify the unbounded graph before this estimate, qε,jq_{\varepsilon,j} is compact smooth for each fixed positive ε\varepsilon and finite jj. Its full Weyl kernel is Schwartz by (CI4), including the original inverse Fourier factor and coordinate change with absolute Jacobian one. CT11a maps every HνH_\nu input into a Schwartz output. Hence that output lies in the coefficient domain of PjP_j. Apply (IP24) to each vector of an orthonormal input basis, use (RA7) and CI5 for each weighted output, and sum the nonnegative squared norms. Tonelli and the finite sum inequality prove (RA8). This also proves that the displayed composition is an actual Hilbert–Schmidt operator on all of HνH_\nu.

The exact inverse-domain map in (IP23a) now gives, for every f∈Hνf\in H_\nu, qε,jwf=Pj−1(Pjqε,jwf),∥qε,jw∥1≤∥Pj−1∥2∥Pjqε,jw∥2≤Cn,M,χ,kν(2π)−n/2Qε2M−2nRj2n−2M.(RA9) q_{\varepsilon,j}^w f =P_j^{-1}(P_j q_{\varepsilon,j}^w f),\qquad \|q_{\varepsilon,j}^w\|_1 \leq\|P_j^{-1}\|_2\|P_j q_{\varepsilon,j}^w\|_2 \leq C_{n,M,\chi,k}\sqrt\nu(2\pi)^{-n/2}Q \varepsilon^{2M-2n}R_j^{2n-2M}. \tag{RA9} The ideal product inequality is (T7). The last equality of scale factors is exactly Ajδj−2n=ε2MRj−2M(ε/Rj)−2n=ε2M−2nRj2n−2MA_j\delta_j^{-2n} =\varepsilon^{2M}R_j^{-2M}(\varepsilon/R_j)^{-2n} =\varepsilon^{2M-2n}R_j^{2n-2M}. Neither the configuration dimension nn nor the phase dimension 2n2n is exchanged. The ν\sqrt\nu shown in (RA9) comes from the inverse factor; the full symbol fiber norm remains in QQ.

Identification of the sum with the original operator. Since M>nM>n, the exact annular sum is ∑j=0∞Rj2n−2M=11−2−2(M−n),∑j=J+1∞Rj2n−2M=2−2(J+1)(M−n)1−2−2(M−n).(RA10) \sum_{j=0}^\infty R_j^{2n-2M} ={1\over1-2^{-2(M-n)}},\qquad \sum_{j=J+1}^\infty R_j^{2n-2M} ={2^{-2(J+1)(M-n)}\over1-2^{-2(M-n)}}. \tag{RA10} By (RA9) and completeness of the trace ideal, the operator sum converges in trace norm. For each fixed ε\varepsilon, the original symbol lies in L2L^2: its full zeroth estimate bounds its squared norm by Q2hε2MQ^2h_\varepsilon^{2M}, whose integral is finite by (IP17) with 2M>n2M>n. Equation (RA2) converges to qεq_\varepsilon in that original matrix Lz2L^2_z norm by dominated convergence, because χ(εz/2J)→1\chi(\varepsilon z/2^J)\to1 and 0≤χ≤10\leq\chi\leq1. CI5 makes the finite operator sums converge in Hilbert–Schmidt norm to the original qεwq_\varepsilon^w. Trace-norm convergence implies operator-norm convergence, as does Hilbert–Schmidt convergence. Uniqueness of the operator-norm limit identifies both limits exactly. Summing (RA9) with (RA10) proves (RA1), including a bound with the displayed denominator. This argument verifies the original operator rather than assigning a Weyl symbol to a different trace-ideal limit. ▫\square

A two-coordinate phase-space section showing the exact annular supports and trace-ideal factors

The left drawing is the (x1,ξ1)(x_1,\xi_1) section with the other original phase coordinates zero, expressed by the comparison coordinates w=εzw=\varepsilon z. It shows the exact supports Rj/2≤|w|≤2RjR_j/2\leq|w|\leq2R_j for j=1,2j=1,2; their overlap is intentional. The right side records (RA3), (RA5), (RA8), and (RA9)–(RA10), not a numerical approximation to an operator. The proof applies in all 2n2n phase dimensions. The reproducible drawing keeps the original symbol, its factors, both Hilbert–Schmidt maps, and the geometric convergence condition.

17. Sharpness, exact trace, and propagation to the ordered coefficients

The condition M>nM>n is sharp for the whole symbol class. For any real MM, keep the actual family qε(z)=hε(z)MIν=(ε−2+|x|2+|ξ|2)−MIν.(RA11) q_\varepsilon(z)=h_\varepsilon(z)^M I_\nu =(\varepsilon^{-2}+|x|^2+|\xi|^2)^{-M}I_\nu. \tag{RA11} The complete recursion (RP3)–(RP4) with s=Ms=M proves ∥∂σqε(z)∥HS≤νCM,σhε(z)M+|σ|/2\|\partial^\sigma q_\varepsilon(z)\|_{\mathrm{HS}} \leq\sqrt\nu C_{M,\sigma}h_\varepsilon(z)^{M+|\sigma|/2} for every derivative order, uniformly in ε\varepsilon. To pass to the exact directional seminorm, expand each directional derivative into coordinate derivatives; for gε,z(Ti)≤1g_{\varepsilon,z}(T_i)\leq1, each |Ti|≤hε(z)−1/2|T_i|\leq h_\varepsilon(z)^{-1/2}, and the finite coordinate sum cancels precisely the derivative factor hε|σ|/2h_\varepsilon^{|\sigma|/2}. Thus (RA11) satisfies the full stated symbol hypothesis even for M≤nM\leq n. At ε=1\varepsilon=1, it is the original rational symbol of (WT6) in the trace criterion. That result proves failure of trace class for every real M≤nM\leq n, so a universal conclusion in (RA1) cannot hold there.

For completeness the necessity used here is visible in the original Weyl map. Its positive Gaussian Wigner function is exactly 2ne−|u−q|2−|ξ−p|22^n e^{-|u-q|^2-|\xi-p|^2}, with full integral (2π)n(2\pi)^n in (q,p)(q,p). If (RA11) at ε=1\varepsilon=1 were trace class, the absolute rank-one integral bound (CI9)–(CI10), summed over the ν\nu fibers, and nonnegative Tonelli would give ∥q1w∥1≥(2π)−nν∫ℝ2n(1+|x|2+|ξ|2)−Mdxdξ.(RA12) \|q_1^w\|_1 \geq (2\pi)^{-n}\nu \int_{\mathbb R^{2n}}(1+|x|^2+|\xi|^2)^{-M} \,dx\,d\xi. \tag{RA12} The symbol has polynomial growth for every real MM, so its Gaussian pairings are absolutely defined before Tonelli. The radial tail behaves as r2n−1−2Mr^{2n-1-2M}: at M=nM=n the integral of r−1r^{-1} diverges, and for M<nM<n the power also diverges. This proves the needed obstruction with the original Fourier and fiber factors.

For M>nM>n, (RA1) supplies trace class and the original zeroth symbol bound supplies L1L^1, so CI12 gives the full trace and bound Tr⁡qεw=(2π)−n∫ℝ2ntr⁡ℂνqε(x,ξ)dxdξ,|Tr⁡qεw|≤(2π)−nνQ0ε2M−2nπnΓ(M−n)Γ(M),Q0=supε,z∥qε(z)∥ophε(z)M.(RA13) \begin{aligned} \operatorname{Tr}q_\varepsilon^w &=(2\pi)^{-n}\int_{\mathbb R^{2n}} \operatorname{tr}_{\mathbb C^\nu} q_\varepsilon(x,\xi)\,dx\,d\xi,\\ |\operatorname{Tr}q_\varepsilon^w| &\leq (2\pi)^{-n}\nu Q_0 \varepsilon^{2M-2n}\pi^n {\Gamma(M-n)\over\Gamma(M)}, \quad Q_0=\sup_{\varepsilon,z} {\|q_\varepsilon(z)\|_{\mathrm{op}} \over h_\varepsilon(z)^M}. \end{aligned} \tag{RA13} The identity uses both proved hypotheses, not an inference of trace class from L1L^1 alone. For the family (RA11), its trace is exactly (2π)−nνε2M−2nπnΓ(M−n)/Γ(M)>0(2\pi)^{-n}\nu\varepsilon^{2M-2n}\pi^n\Gamma(M-n)/\Gamma(M)>0. Since |Tr⁡T|≤∥T∥1|\operatorname{Tr}T|\leq\|T\|_1, no bound with a strictly larger positive power of ε\varepsilon can hold uniformly for this whole class as ε↓0\varepsilon\downarrow0. Thus the scale exponent in (RA1), as well as the threshold, is sharp.

This extension applies to every actual ordered remainder product through (RP2). It leaves the truncation integer in (IP18), the finite polynomial (DE4), and all the original CkC_k unchanged. For (IP21), keep its exact telescoping order u#ℓ#(u−v)#v#(N−1−ℓ)u^{\#\ell}\#(u-v)\#v^{\#(N-1-\ell)}. The preceding symbol inclusion and (RA1) give trace norm O(ε2M−2n)O(\varepsilon^{2M-2n}) for each of its finitely many terms. In particular the original M=n+1M=n+1 case and (DE5) with N=M=n+1N=M=n+1 give ∥(Ej,ε(n+1,n+1))w∥1=O(ε2).(RA14) \|(E_{j,\varepsilon}^{(n+1,n+1)})^w\|_1 =O(\varepsilon^2). \tag{RA14} Every low-degree coefficient remains compactly supported by the original factor count (DE6), so its individual trace keeps the full (2π)−nε2k−2n(2\pi)^{-n}\varepsilon^{2k-2n} factor (DE8). The constant index and the finite power comparison still give (DE10)–(DE11) without any formal infinite-series assertion. The real exponent concerns a genuine uniform remainder; it cannot be applied to (1−ψε)Iν(1-\psi_\varepsilon)I_\nu or to the full error without its required hypothesis.

18. Editorial supplement: the complete finite ordered coefficient calculation

This provides the direct finite proof of the associativity stated in Section9, and clarifies the complete differential count in Section12. In a binary degree-kk coefficient, the outer contraction degree is only ℓ\ell; the other intrinsic degrees are already inside its inputs. After expanding those inputs the full count is exactly 2k2k, as proved below. The original statements and formulas remain identifiable above. The proof also extends the support assertion from N=n+1N=n+1 to every N≥1N\geq1 and k<Nk<N.

This derivation checks DE4, the differential count in DE11, and the support assertion DE6. It uses the actual coefficient in W26, on the original phase space W=ℝxn⊕ℝξnW=\mathbb R_x^n\oplus\mathbb R_\xi^n, with the original convention σ(P,Q)=q⋅r−p⋅s,P=(p,q),Q=(r,s),D=−i∂,𝒜(P,Q)=12σ(P,Q).(OC1) \sigma(P,Q)=q\cdot r-p\cdot s, \quad P=(p,q),\quad Q=(r,s),\qquad D=-i\partial, \qquad \mathcal A(P,Q)=\frac12\sigma(P,Q). \tag{OC1} All matrix products below have their displayed order. No coefficient sequence is interpreted as a convergent series.

18.1. The exact binary coefficient

For independent phase variables Xr=(xr,ξr)X_r=(x_r,\xi_r) and Xs=(xs,ξs)X_s=(x_s,\xi_s), put Lrs=i𝒜(DXr,DXs)=i2∑a=1n(Dξr,aDxs,a−Dxr,aDξs,a)=i2∑a=1n(∂xr,a∂ξs,a−∂ξr,a∂xs,a).(OC2) \begin{aligned} L_{rs} &=i\mathcal A(D_{X_r},D_{X_s})\\ &=\frac i2\sum_{a=1}^n \left(D_{\xi_{r,a}}D_{x_{s,a}} -D_{x_{r,a}}D_{\xi_{s,a}}\right)\\ &=\frac i2\sum_{a=1}^n \left(\partial_{x_{r,a}}\partial_{\xi_{s,a}} -\partial_{\xi_{r,a}}\partial_{x_{s,a}}\right). \end{aligned} \tag{OC2} The second equality retains W26’s phase σ/2\sigma/2; the third follows from the two factors −i-i in each summand. In particular the sign of the xr,ξsx_r,\xi_s term is positive. If δ2\delta_2 denotes restriction to Xr=Xs=XX_r=X_s=X, W26 is exactly Cℓ(f,g)=1ℓ!δ2Lrsℓ(f(Xr)g(Xs)).(OC3) C_\ell(f,g)=\frac1{\ell!}\delta_2 L_{rs}^{\ell}\bigl(f(X_r)g(X_s)\bigr). \tag{OC3} The coordinate derivatives commute, and their coefficients are scalars. The finite multinomial rule therefore gives the complete formula Cℓ(f,g)(X)=(i2)ℓ∑α,β∈ℕn|α|+|β|=ℓ(−1)|β|α!β!(∂xα∂ξβf)(X)(∂ξα∂xβg)(X).(OC4) C_\ell(f,g)(X) =\left(\frac i2\right)^\ell \sum_{\substack{\alpha,\beta\in\mathbb N^n\\ |\alpha|+|\beta|=\ell}} \frac{(-1)^{|\beta|}}{\alpha!\,\beta!} (\partial_x^\alpha\partial_\xi^\beta f)(X) (\partial_\xi^\alpha\partial_x^\beta g)(X). \tag{OC4} Indeed the multinomial coefficient in LrsℓL_{rs}^\ell is ℓ!/(α!β!)\ell!/(\alpha!\beta!), and OC3 retains its divisor ℓ!\ell!. Every summand differentiates ff exactly ℓ\ell times and gg exactly ℓ\ell times. At ℓ=0\ell=0, OC4 gives fgfg. At ℓ=1\ell=1, it gives i2∑a=1n((∂xaf)(∂ξag)−(∂ξaf)(∂xag)),(OC5) \frac i2\sum_{a=1}^n \bigl((\partial_{x_a}f)(\partial_{\xi_a}g) -(\partial_{\xi_a}f)(\partial_{x_a}g)\bigr), \tag{OC5} which is W29 with its stated symplectic convention. In particular C1(x2,ξ2)=2ixξC_1(x^2,\xi^2)=2ix\xi, preserving the sign in the source’s polynomial example.

18.2. A finite triple coefficient identity

Let f,g,hf,g,h be arbitrary smooth square matrix functions on WW. Write T(X1,X2,X3)=f(X1)g(X2)h(X3),δ3T(X)=T(X,X,X).(OC6) T(X_1,X_2,X_3)=f(X_1)g(X_2)h(X_3), \qquad \delta_3T(X)=T(X,X,X). \tag{OC6} No factor is moved past another. The ordinary chain rule for diagonal restriction says, for each original coordinate zaz_a, ∂za[U(X1,X2)]X1=X2=X=[(∂z1,a+∂z2,a)U(X1,X2)]X1=X2=X.(OC7) \partial_{z_a}\bigl[U(X_1,X_2)\bigr]_{X_1=X_2=X} =\bigl[(\partial_{z_{1,a}}+\partial_{z_{2,a}}) U(X_1,X_2)\bigr]_{X_1=X_2=X}. \tag{OC7} Repeated derivatives obey the repeated version of OC7 because all its derivative operators commute. Apply it to the inner coefficient Cp(f,g)C_p(f,g) in Cq(Cp(f,g),h)C_q(C_p(f,g),h). The outer left derivatives become the sum of derivatives in slots 1 and 2, so their contraction against slot 3 becomes L13+L23L_{13}+L_{23}. Thus, for all p,q≥0p,q\ge0, Cq(Cp(f,g),h)=1p!q!δ3(L13+L23)qL12pT.(OC8) C_q(C_p(f,g),h) =\frac1{p!\,q!}\delta_3 (L_{13}+L_{23})^qL_{12}^pT. \tag{OC8} Apply the same chain rule to the inner coefficient Cq(g,h)C_q(g,h) in the other parenthesization. The outer right derivatives become the sum of derivatives in slots 2 and 3, and hence Cp(f,Cq(g,h))=1p!q!δ3(L12+L13)pL23qT.(OC9) C_p(f,C_q(g,h)) =\frac1{p!\,q!}\delta_3 (L_{12}+L_{13})^pL_{23}^qT. \tag{OC9} Each of L12,L13,L23L_{12},L_{13},L_{23} is a constant coefficient differential operator with scalar coefficients. Every coordinate derivative in one commutes with every coordinate derivative in the others, including derivatives in a shared slot. Consequently these three operators commute on every matrix entry. This commutation reorders differential operators only; it never reorders matrix factors.

Fix m≥0m\ge0. Summing OC8 over p+q=mp+q=m and applying the finite binomial rule gives ∑p+q=mCq(Cp(f,g),h)=δ3∑p+q=mL12p(L13+L23)qp!q!T=δ3∑a+b+c=mL12aL13bL23ca!b!c!T.(OC10) \begin{aligned} \sum_{p+q=m}C_q(C_p(f,g),h) &=\delta_3 \sum_{p+q=m}\frac{L_{12}^p(L_{13}+L_{23})^q}{p!\,q!}T\\ &=\delta_3 \sum_{a+b+c=m}\frac{L_{12}^aL_{13}^bL_{23}^c}{a!\,b!\,c!}T. \end{aligned} \tag{OC10} Summing OC9 over the same total degree gives ∑p+q=mCp(f,Cq(g,h))=δ3∑p+q=m(L12+L13)pL23qp!q!T=δ3∑a+b+c=mL12aL13bL23ca!b!c!T.(OC11) \begin{aligned} \sum_{p+q=m}C_p(f,C_q(g,h)) &=\delta_3 \sum_{p+q=m}\frac{(L_{12}+L_{13})^pL_{23}^q}{p!\,q!}T\\ &=\delta_3 \sum_{a+b+c=m}\frac{L_{12}^aL_{13}^bL_{23}^c}{a!\,b!\,c!}T. \end{aligned} \tag{OC11} Every sum in OC8–OC11 is finite. The two right sides are identical, proving ∑p+q=mCq(Cp(f,g),h)=∑p+q=mCp(f,Cq(g,h)).(OC12) \boxed{\sum_{p+q=m}C_q(C_p(f,g),h) =\sum_{p+q=m}C_p(f,C_q(g,h)).} \tag{OC12} This is an exact differential identity for arbitrary smooth matrix functions. It needs neither analytic Weyl-product associativity nor a remainder estimate nor a uniqueness argument for asymptotic coefficients.

18.3. Associativity of DE4 through every finite degree

For three coefficient sequences (ui),(vj),(wr)(u_i),(v_j),(w_r), DE4 gives ((u⋆tv)⋆tw)k=∑i+j+r+p+q=kCq(Cp(ui,vj),wr),(u⋆t(v⋆tw))k=∑i+j+r+p+q=kCp(ui,Cq(vj,wr)).(OC13) \begin{aligned} ((u\star_t v)\star_t w)_k &=\sum_{i+j+r+p+q=k}C_q(C_p(u_i,v_j),w_r),\\ (u\star_t(v\star_t w))_k &=\sum_{i+j+r+p+q=k}C_p(u_i,C_q(v_j,w_r)). \end{aligned} \tag{OC13} All indices are nonnegative integers. Fix i,j,ri,j,r with i+j+r≤ki+j+r\le k and apply OC12 with m=k−i−j−rm=k-i-j-r. Summing those equalities over the finitely many i,j,ri,j,r proves equality of the two coefficients in OC13. This proves DE4’s associativity for every finite kk.

For any chosen KK, coefficients of degree greater than KK can be discarded throughout the calculation: every intrinsic degree and every contraction degree is nonnegative, so a discarded degree can never enter a coefficient of degree at most KK in a later product. Thus the proof establishes an associative operation on the finite sequences of degrees 0,…,K0,\ldots,K directly.

18.4. The complete ordered N-factor coefficient

For N≥1N\ge1, use independent original phase variables X1,…,XNX_1,\ldots,X_N. For each 1≤r<s≤N1\le r<s\le N keep the operator LrsL_{rs} of OC2, and let δN\delta_N restrict all variables to XX. The coefficient defined recursively from DE4 is exactly Fj,k(N)=∑i1,…,iN≥0,drs≥0∑s=1Nis+∑r<sdrs=k1∏r<sdrs!δN(∏r<sLrsdrs)(cj,i1(X1)⋯cj,iN(XN)).(OC14) F_{j,k}^{(N)} =\sum_{\substack{i_1,\ldots,i_N\ge0,\ d_{rs}\ge0\\ \sum_{s=1}^Ni_s+\sum_{r<s}d_{rs}=k}} \frac1{\prod_{r<s}d_{rs}!} \delta_N\left(\prod_{r<s}L_{rs}^{d_{rs}}\right) \bigl(c_{j,i_1}(X_1)\cdots c_{j,i_N}(X_N)\bigr). \tag{OC14} For N=1N=1 all edge products and sums are empty and OC14 is Fj,k(1)=cj,kF_{j,k}^{(1)}=c_{j,k}. To prove the induction from N−1N-1 to NN, multiply the previous coefficient on the right by the new final coefficient using DE4 in the original order F(N−1)F^{(N-1)} followed by cj,iNc_{j,i_N}. The chain rule replaces each derivative on the diagonal of the first N−1N-1 slots by their sum. An outer contraction of degree qq is therefore 1q!(∑r=1N−1LrN)q=∑d1N,…,dN−1,N≥0∑r=1N−1drN=q∏r=1N−1LrNdrNdrN!.(OC15) \frac1{q!}\left(\sum_{r=1}^{N-1}L_{rN}\right)^q =\sum_{\substack{d_{1N},\ldots,d_{N-1,N}\ge0\\ \sum_{r=1}^{N-1}d_{rN}=q}} \prod_{r=1}^{N-1}\frac{L_{rN}^{d_{rN}}}{d_{rN}!}. \tag{OC15} It commutes with all old edge operators and produces precisely OC14. Associativity already proved in OC13 gives the same formula for every other parenthesization preserving the leaf order.

The full coordinate version records every edge coefficient. Choose αrs,βrs∈ℕn\alpha_{rs},\beta_{rs}\in\mathbb N^n for each r<sr<s, put m=∑r<s(|αrs|+|βrs|),As=∑t>sαst+∑r<sβrs,Bs=∑t>sβst+∑r<sαrs.(OC16) \begin{aligned} m&=\sum_{r<s}(|\alpha_{rs}|+|\beta_{rs}|),\\ A_s&=\sum_{t>s}\alpha_{st}+\sum_{r<s}\beta_{rs},\\ B_s&=\sum_{t>s}\beta_{st}+\sum_{r<s}\alpha_{rs}. \end{aligned} \tag{OC16} Then OC14 expands to Fj,k(N)(X)=∑i1,…,iN≥0,αrs,βrs∈ℕn∑sis+m=k(i2)m(−1)∑r<s|βrs|∏r<sαrs!βrs!∏s=1N(∂xAs∂ξBscj,is)(X).(OC17) F_{j,k}^{(N)}(X) =\sum_{\substack{i_1,\ldots,i_N\ge0,\ \alpha_{rs},\beta_{rs}\in\mathbb N^n\\ \sum_s i_s+m=k}} \left(\frac i2\right)^m \frac{(-1)^{\sum_{r<s}|\beta_{rs}|}} {\prod_{r<s}\alpha_{rs}!\,\beta_{rs}!} \prod_{s=1}^{N} \bigl(\partial_x^{A_s}\partial_\xi^{B_s}c_{j,i_s}\bigr)(X). \tag{OC17} The last product is explicitly ordered by increasing ss. Each edge contributes one derivative to each of its two endpoints, with the two original signs in OC2. Therefore ∑s=1N(|As|+|Bs|)=2m.(OC18) \sum_{s=1}^N(|A_s|+|B_s|)=2m. \tag{OC18} No factorial, phase factor, derivative or matrix order is omitted in OC14 or OC17.

18.5. Exactly 2k derivatives after expanding DE1

Keep the two original pairs (f1,g1)=(b,a)(f_1,g_1)=(b,a), (f2,g2)=(a,b)(f_2,g_2)=(a,b), with b=ψa−1b=\psi a^{-1}. DE1 is cj,0=(1−ψ)Iν,cj,q=−Cq(fj,gj)(q≥1).(OC19) c_{j,0}=(1-\psi)I_\nu, \qquad c_{j,q}=-C_q(f_j,g_j)\quad(q\ge1). \tag{OC19} Let γ∈ℕ2n\gamma\in\mathbb N^{2n} be a full phase multiindex, write (α,β)(\alpha,\beta) for its spatial and frequency components when needed, and use the componentwise relation δ≤γ\delta\le\gamma. From OC4 and the complete matrix-valued product rule, for q≥1q\ge1, ∂γcj,q=−(i2)q∑|α|+|β|=q∑δ≤γ(−1)|β|α!β!(γδ)⋅(∂δ+(α,β)fj)(∂γ−δ+(β,α)gj).(OC20) \begin{aligned} \partial^\gamma c_{j,q} ={}&-\left(\frac i2\right)^q \sum_{|\alpha|+|\beta|=q} \sum_{\delta\le\gamma} \frac{(-1)^{|\beta|}}{\alpha!\,\beta!} \binom\gamma\delta\\ &\quad\cdot \bigl(\partial^{\delta+(\alpha,\beta)}f_j\bigr) \bigl(\partial^{\gamma-\delta+(\beta,\alpha)}g_j\bigr). \end{aligned} \tag{OC20} Both factors remain in the original order fj,gjf_j,g_j. Every summand has exactly |δ+(α,β)|+|γ−δ+(β,α)|=2q+|γ|(OC21) |\delta+(\alpha,\beta)| +|\gamma-\delta+(\beta,\alpha)| =2q+|\gamma| \tag{OC21} phase differentiations. For q=0q=0, ∂γcj,0=(∂γ(1−ψ))Iν,(OC22) \partial^\gamma c_{j,0} =\bigl(\partial^\gamma(1-\psi)\bigr)I_\nu, \tag{OC22} which has |γ||\gamma| differentiations on the original cutoff expression.

In an OC17 summand take γs=(As,Bs)\gamma_s=(A_s,B_s), expand every positive intrinsic coefficient with OC20, and retain OC22 for every zero intrinsic coefficient. The complete resulting word has total differential degree ∑s=1N(2is+|γs|)=2∑s=1Nis+2m=2k.(OC23) \sum_{s=1}^{N}(2i_s+|\gamma_s|) =2\sum_{s=1}^{N}i_s+2m =\boxed{2k}. \tag{OC23} Its full phase factor is the outer (i/2)m(i/2)^m times the intrinsic factors (i/2)is(i/2)^{i_s}, hence (i/2)k(i/2)^k. Its coefficient also keeps the minus sign from each positive intrinsic coefficient, every sign (−1)|β|(-1)^{|\beta|} of each contraction, all edge factorials, all intrinsic factorials and every Leibniz binomial in OC20. No matrix factors are exchanged.

The count in OC23 concerns the complete differential words in the retained expansion. Particular words can evaluate to zero, and different words can cancel. These evaluations do not change the count of the derivatives in their defining expressions.

If the original b=ψa−1b=\psi a^{-1} is expanded further at a point where the inverse is defined, the count still holds. For every full phase multiindex η\eta, the exact product rule is ∂η(ψa−1)=∑ρ≤η(ηρ)(∂ρψ)(∂η−ρa−1).(OC24) \partial^\eta(\psi a^{-1}) =\sum_{\rho\le\eta}\binom\eta\rho (\partial^\rho\psi)(\partial^{\eta-\rho}a^{-1}). \tag{OC24} For λ≠0\lambda\ne0, differentiating aa−1=Iνaa^{-1}=I_\nu yields the full ordered inverse formula ∂λa−1=∑r=1|λ|(−1)r∑λ1,…,λr∈ℕ2n\{0}λ1+⋯+λr=λλ!λ1!⋯λr!a−1(∂λ1a)a−1⋯(∂λra)a−1.(OC25) \partial^\lambda a^{-1} =\sum_{r=1}^{|\lambda|}(-1)^r \sum_{\substack{\lambda_1,\ldots,\lambda_r\in\mathbb N^{2n}\setminus\{0\}\\ \lambda_1+\cdots+\lambda_r=\lambda}} \frac{\lambda!}{\lambda_1!\cdots\lambda_r!} a^{-1}(\partial^{\lambda_1}a)a^{-1} \cdots(\partial^{\lambda_r}a)a^{-1}. \tag{OC25} For λ=0\lambda=0 retain a−1a^{-1} itself. To verify OC25 without assuming it, the product rule gives the recurrence ∂λa−1=−a−1∑0<τ≤λ(λτ)(∂τa)(∂λ−τa−1).(OC26) \partial^\lambda a^{-1} =-a^{-1}\sum_{0<\tau\le\lambda} \binom\lambda\tau(\partial^\tau a) (\partial^{\lambda-\tau}a^{-1}). \tag{OC26} Induction on |λ||\lambda| in OC26 attaches the first nonzero part τ\tau to every ordered composition of λ−τ\lambda-\tau, with coefficient (λτ)(λ−τ)!/(λ2!⋯λr!)=λ!/(τ!λ2!⋯λr!)\binom\lambda\tau(\lambda-\tau)!/(\lambda_2!\cdots\lambda_r!) =\lambda!/(\tau!\lambda_2!\cdots\lambda_r!) and the additional minus sign. The zero remainder gives the r=1r=1 term. This proves every term and coefficient in OC25. Each inverse term has total derivative order ∑l|λl|=|λ|\sum_l|\lambda_l|=|\lambda| on the original aa, and each term of OC24 therefore has total order |ρ|+|η−ρ|=|η||\rho|+|\eta-\rho|=|\eta|, preserving OC23. Where the cutoff construction uses a smooth extension of bb, OC20–OC23 apply to that original smooth bb directly; OC25 is asserted only on the inverse’s actual domain.

18.6. Compact support for every k<N

Let K=supp⁡(1−ψ).(OC27) K=\operatorname{supp}(1-\psi). \tag{OC27} It is a fixed compact set because ψ=1\psi=1 outside a compact set. In every OC14 or OC17 summand, ∑s=1Nis≤k.(OC28) \sum_{s=1}^{N}i_s\le k. \tag{OC28} The number of positive isi_s is at most ∑sis\sum_s i_s, hence at most kk. For k<Nk<N, there are at least N−kN-k slots whose intrinsic degree is zero. In particular at least one ordered factor in that summand is a derivative of cj,0c_{j,0}.

For every γ\gamma, supp⁡(∂γcj,0)⊂K\operatorname{supp}(\partial^\gamma c_{j,0})\subset K: at a point outside KK, cj,0c_{j,0} is identically zero on a neighborhood, and all of its derivatives vanish there. The entire ordered matrix product is therefore zero outside KK, regardless of the position of this factor and regardless of the other matrix factors. The sum is finite, and all factors are smooth. Consequently Fj,k(N)∈Cc∞(W;End⁡ℂν),supp⁡Fj,k(N)⊂K,0≤k<N.(OC29) \boxed{F_{j,k}^{(N)}\in C_c^\infty (W;\operatorname{End}\mathbb C^\nu),\qquad \operatorname{supp}F_{j,k}^{(N)}\subset K, \qquad 0\le k<N.} \tag{OC29} Taking N=n+1N=n+1 recovers DE6 for 0≤k≤n0\le k\le n. This support argument is valid in both original orders independently and uses no equality between Cq(b,a)C_q(b,a) and Cq(a,b)C_q(a,b).

The same full derivative count verifies the symbol weight used in DE5. On the original metric g=h(|dx|2+|dξ|2)g=h(|dx|^2+|d\xi|^2), DE1 gives |∂γcj,i|≤Ci,γhi+|γ|/2|\partial^\gamma c_{j,i}|\le C_{i,\gamma}h^{i+|\gamma|/2}, including i=0i=0 by compact support and positivity of hh on KK. If qq extra derivatives are applied to an OC17 word, the full product rule distributes all qq among its factors. Each resulting term is bounded by Ch∑sis+(∑s|γs|+q)/2=Chk+q/2.(OC30) C h^{\sum_s i_s+(\sum_s|\gamma_s|+q)/2} =C h^{k+q/2}. \tag{OC30} All sums are finite for fixed N,k,qN,k,q, which proves Fj,k(N)∈S(hk,g)F_{j,k}^{(N)}\in S(h^k,g) with the original weight and coordinates.

18.7. The exact zeroth coefficient in both original orders

Degree zero in the full formula (OC14) forces every intrinsic degree and every contraction degree to be zero. Both sequences have the same original zeroth coefficient cj,0=(1−ψ)Iνc_{j,0}=(1-\psi)I_\nu. Therefore, for every integer N≥1N\geq1, F1,0(N)=F2,0(N)=(1−ψ)NIν,F1,0(N)−F2,0(N)=0pointwise on ℝ2n.(OC31) F_{1,0}^{(N)}=F_{2,0}^{(N)}=(1-\psi)^N I_\nu, \qquad F_{1,0}^{(N)}-F_{2,0}^{(N)}=0 \quad\hbox{pointwise on }\mathbb R^{2n}. \tag{OC31} For N=n+1N=n+1, this proves the original A0=0A_0=0 directly, retaining each of its trace and Fourier factors in (DE9). The index-constancy argument (DE10)–(DE11) still proves the cancellation of every other lower coefficient and the exact surviving coefficient AnA_n. No pointwise equality of positive coefficients in the two matrix orders is assumed.