Tangent zoom and quadratic models

The local estimates below apply to every intrinsic ImI^m distribution through the completed earlier graph theorem. The Gaussian-symbol companion proves the exact intrinsic normalization and global Gaussian/Maslov identification.

At a selected nonzero covector, modulation removes a rapidly oscillating linear phase. Spatial magnification then leaves a quadratic phase and a distribution supported on the range of its frequency Hessian. A singular Hessian produces a delta distribution in the complementary directions. An ordinary symbol can retain an oscillating leading coefficient, so the quadratic model need not be a limit with a fixed coefficient.

1. Conventions and the earlier proofs

All distributions and pairings are complex linear. Put Dj=−i∂jD_j=-i\partial_j, f^(ξ)=∫e−ix⋅ξf(x) dx\widehat f(\xi)=\int e^{-ix\cdot\xi}f(x)\,dx, and use inverse factor (2π)−n(2\pi)^{-n}. A test is a smooth compactly supported function. On a fixed compact KK, write

pK,N(ϕ)=max⁡∣α∣≤Nsup⁡∣∂αϕ∣. p_{K,N}(\phi)=\max_{|\alpha|\leq N}\sup|\partial^\alpha\phi|.

A distribution obeys ∣v(ϕ)∣≤CKpK,NK(ϕ)|v(\phi)|\leq C_Kp_{K,N_K}(\phi) for each such KK. Weak convergence means convergence on every fixed test. In this lesson vt=OD′(t−a)v_t=O_{\mathcal D'}(t^{-a}) means that for every fixed compact test support there are C,N,t0C,N,t_0 such that

∣vt(ϕ)∣≤Ct−apK,N(ϕ),t≥t0.(1.1) |v_t(\phi)|\leq Ct^{-a}p_{K,N}(\phi),\qquad t\geq t_0. \tag{1.1}

This includes uniformity over all tests with the indicated seminorm bounded.

The earlier quadratic stationary-phase component proves the integration limit rule Q1, complex Gaussian Q2, Schwartz Fourier estimates Q3, inversion Q4, symmetric diagonalization Q5 and quadratic Fourier identity Q6. Its complete F0 chains supply the elementary calculus, integration, compactness and matrix results used here. The earlier stationary lesson, Appendix A.4, supplies smooth cutoffs. The selected uniform-bound proof U1 includes its full complete-test-space and Baire prerequisites. These are programme proof texts, with their separate licences.

A smooth symbol b∈Sr(Rn)b\in S^r(\mathbb R^n), r∈Rr\in\mathbb R, satisfies

∣∂αb(ξ)∣≤Cα⟨ξ⟩r−∣α∣,⟨ξ⟩=(1+∣ξ∣2)1/2.(1.2) |\partial^\alpha b(\xi)|\leq C_\alpha\langle\xi\rangle^{r-|\alpha|}, \qquad \langle\xi\rangle=(1+|\xi|^2)^{1/2}. \tag{1.2}

Let HH be real, smooth away from zero and positively homogeneous of degree one. Choose a real smooth extension HeH_e agreeing with HH for ∣ξ∣≥R0|\xi|\geq R_0. For example multiply HH by a smooth radial cutoff vanishing near zero. Define the graph distribution by

u(x)=cn∫ei(x⋅ξ−He(ξ))b(ξ) dξ,cn=(2π)−3n/4.(1.3) u(x)=c_n\int e^{i(x\cdot\xi-H_e(\xi))}b(\xi)\,d\xi, \qquad c_n=(2\pi)^{-3n/4}. \tag{1.3}

The integral means u(ϕ)=cn∫e−iHe(ξ)b(ξ)ϕ^(−ξ) dξu(\phi)=c_n\int e^{-iH_e(\xi)}b(\xi)\widehat\phi(-\xi)\,d\xi. It converges absolutely and is bounded by finitely many Schwartz seminorms of ϕ\phi, by (1.2) and Q3. Thus it defines a tempered distribution and, on compact test supports, a distribution. Its Fourier transform is (2π)n/4e−iHeb(2\pi)^{n/4}e^{-iH_e}b, by Q4 and transposition.

The normalization corresponding to the programme's graph order is r=m−n/4r=m-n/4. The earlier prescribed-phase theorem F6 and intrinsic localization K7 supply the localized graph representation of an intrinsic ImI^m distribution modulo an error regular at the selected covector. The forward graph criterion puts its reduced Fourier coefficient in Sm−n/4S^{m-n/4}; multiplying it by (2π)−n/4(2\pi)^{-n/4} gives exactly the amplitude bb in (1.3). T4 below handles the regular remainder. The companion's Z1 records this argument with all normalization factors.

2. The frequency-graph zoom

Fix ξ0≠0\xi_0\ne0, set x0=H′(ξ0)x_0=H'(\xi_0), and choose a real smooth ψ\psi near x0x_0 with ψ(x0)=0\psi(x_0)=0, ψ′(x0)=ξ0\psi'(x_0)=\xi_0. Write

A=H′′(ξ0),B=ψ′′(x0),QA,B(x,η)=x⋅η−12xTBx−12ηTAη.(2.1) A=H''(\xi_0),\quad B=\psi''(x_0),\quad Q_{A,B}(x,\eta)=x\cdot\eta-\tfrac12x^TBx-\tfrac12\eta^TA\eta. \tag{2.1}

Both matrices are real symmetric. Define

UA,B=cn∫eiQA,B(x,η) dη,Wt=t−(2r+n)(ue−it2ψ)(x0+x/t),t≥1.(2.2) U_{A,B}=c_n\int e^{iQ_{A,B}(x,\eta)}\,d\eta, \qquad W_t=t^{-(2r+n)}(u e^{-it^2\psi})(x_0+x/t),\quad t\geq1. \tag{2.2}

For distributions the last expression is the pullback by the affine map x↦x0+x/tx\mapsto x_0+x/t, followed by the scalar factor. Equivalently,

⟨Wt,ϕ⟩=t−2r⟨u,e−it2ψ(y)ϕ(t(y−x0))⟩.(2.3) \langle W_t,\phi\rangle =t^{-2r}\langle u,e^{-it^2\psi(y)}\phi(t(y-x_0))\rangle. \tag{2.3}

For fixed test support this is defined for all sufficiently large tt. In UA,BU_{A,B}, integrate against the test in xx first. The resulting Fourier transform of e−ixTBx/2ϕe^{-ix^TBx/2}\phi decreases faster than every power, so this definition also converges absolutely.

Theorem T1 (the ordinary-symbol comparison). For every symbol (1.2),

Wt−ctUA,B=OD′(t−1),ct=t−2rb(t2ξ0).(2.4) W_t-c_tU_{A,B}=O_{\mathcal D'}(t^{-1}),\qquad c_t=t^{-2r}b(t^2\xi_0). \tag{2.4}

The coefficients ctc_t are bounded; they need not converge. The result holds componentwise for symbols in any fixed finite-dimensional vector space. The weaker conclusion that the difference tends to zero follows without choosing any leading homogeneous coefficient.

Proof. Differentiating H(sξ)=sH(ξ)H(s\xi)=sH(\xi) in ss at one gives Euler's identity H(ξ)=ξ⋅H′(ξ)H(\xi)=\xi\cdot H'(\xi). Introduce

bt(η)=t−2rb(t2ξ0+tη),Rt(η)=He(t2ξ0+tη)−t2H(ξ0)−tx0⋅η,qt(x)=t2(ψ(x0+x/t)−ξ0⋅x/t),Ft(η)=∫ei(x⋅η−qt(x))ϕ(x) dx.(2.5) \begin{split} b_t(\eta)&=t^{-2r}b(t^2\xi_0+t\eta),\\ R_t(\eta)&=H_e(t^2\xi_0+t\eta)-t^2H(\xi_0)-t x_0\cdot\eta,\\ q_t(x)&=t^2\bigl(\psi(x_0+x/t)-\xi_0\cdot x/t\bigr),\\ F_t(\eta)&=\int e^{i(x\cdot\eta-q_t(x))}\phi(x)\,dx. \end{split} \tag{2.5}

The substitution ξ=t2ξ0+tη\xi=t^2\xi_0+t\eta, justified in the absolutely convergent test pairing, and Euler's identity give exactly

⟨Wt,ϕ⟩=cn∫bt(η)e−iRt(η)Ft(η) dη.(2.6) \langle W_t,\phi\rangle =c_n\int b_t(\eta)e^{-iR_t(\eta)}F_t(\eta)\,d\eta. \tag{2.6}

The Jacobian is tnt^n; it cancels the t−nt^{-n} in (2.2).

The integral Taylor formula, obtained by applying the fundamental theorem twice on a segment, gives

qt(x)=∫01(1−s) xTψ′′(x0+sx/t)x ds.(2.7) q_t(x)=\int_0^1(1-s)\,x^T\psi''(x_0+sx/t)x\,ds. \tag{2.7}

On each fixed compact KK, its difference from xTBx/2x^TBx/2 is O(t−1)O(t^{-1}) in every derivative. Indeed subtract ψ′′(x0)\psi''(x_0) and apply the fundamental theorem once more; every derivative of the resulting compact integral is bounded, with its explicit factor 1/t1/t. Product and chain rules then show, by Q3, that for each LL there are N,CN,C with

∣Ft(η)∣≤CpK,N(ϕ)⟨η⟩−L,∣Ft(η)−F∞(η)∣≤Ct−1pK,N(ϕ)⟨η⟩−L,(2.8) |F_t(\eta)|\leq Cp_{K,N}(\phi)\langle\eta\rangle^{-L},\quad |F_t(\eta)-F_\infty(\eta)| \leq Ct^{-1}p_{K,N}(\phi)\langle\eta\rangle^{-L}, \tag{2.8}

where F∞=∫eix⋅η−ixTBx/2ϕ(x) dxF_\infty=\int e^{ix\cdot\eta-ix^TBx/2}\phi(x)\,dx. This argument uses the fully proved Fourier derivative estimates; it does not drop a divergence term from a transposed vector field.

Choose 0<c<∣ξ0∣/20<c<|\xi_0|/2. In ∣η∣≤ct|\eta|\leq ct the frequency segment between t2ξ0t^2\xi_0 and t2ξ0+tηt^2\xi_0+t\eta has size comparable to t2t^2. The symbol derivative bound and one segment integral give

∣bt(η)−ct∣≤Ct−1∣η∣.(2.9) |b_t(\eta)-c_t|\leq Ct^{-1}|\eta|. \tag{2.9}

For large tt, homogeneity and the third-order integral remainder give on that same region

Rt(η)=12ηTAη+Et(η),∣Et(η)∣≤Ct−1∣η∣3.(2.10) R_t(\eta)=\tfrac12\eta^TA\eta+E_t(\eta),\qquad |E_t(\eta)|\leq Ct^{-1}|\eta|^3. \tag{2.10}

All points ξ0+sη/t\xi_0+s\eta/t stay in one compact set away from zero, so the required third derivatives of HH have a common bound.

There is also a global estimate

∣bt(η)∣+∣ct∣≤C⟨η⟩2∣r∣.(2.11) |b_t(\eta)|+|c_t|\leq C\langle\eta\rangle^{2|r|}. \tag{2.11}

For r≥0r\geq0, use ⟨t2ξ0+tη⟩≤Ct2⟨η⟩\langle t^2\xi_0+t\eta\rangle\leq Ct^2\langle\eta\rangle. For r<0r<0, the near region just considered gives a constant bound; outside it use ⟨t2ξ0+tη⟩r≤1\langle t^2\xi_0+t\eta\rangle^r\leq1 and t2∣r∣≤C⟨η⟩2∣r∣t^{2|r|}\leq C\langle\eta\rangle^{2|r|}.

Subtract ct⟨UA,B,ϕ⟩c_t\langle U_{A,B},\phi\rangle from (2.6). Replacing FtF_t by F∞F_\infty costs at most Ct−1pK,N(ϕ)Ct^{-1}p_{K,N}(\phi), by (2.8), (2.11) and an integrable power bound. In ∣η∣≤ct|\eta|\leq ct, (2.9), (2.10) and ∣eis−eiv∣≤∣s−v∣|e^{is}-e^{iv}|\leq|s-v| for real s,vs,v give the same bound after integration against F∞F_\infty. Outside this region, bound both exponential factors by one and use (2.11). Taking L>n+2∣r∣+4L>n+2|r|+4 in (2.8) makes that tail at most Ct−1pK,N(ϕ)Ct^{-1}p_{K,N}(\phi); the rectangular-shell proof is Q3/P18.3. All constants involve finitely many fixed seminorms. This proves (2.4), including its uniform test estimate. Finitely many component estimates give the vector-valued assertion. □\square

Corollary T2 (classical and complex degrees). If b(ξ)=br(ξ)+Sr−1b(\xi)=b_r(\xi)+S^{r-1} at large frequency, with brb_r homogeneous of real degree rr, then

ct=br(ξ0)+O(t−2),Wt⟶br(ξ0)UA,B.(2.12) c_t=b_r(\xi_0)+O(t^{-2}),\qquad W_t\longrightarrow b_r(\xi_0)U_{A,B}. \tag{2.12}

The error of the full distribution is still in general only OD′(t−1)O_{\mathcal D'}(t^{-1}). For a leading term of complex degree r+iνr+i\nu, ν∈R\nu\in\mathbb R, with remainder in Sr−1S^{r-1},

t−2iνWt⟶br+iν(ξ0)UA,B.(2.13) t^{-2i\nu}W_t\longrightarrow b_{r+i\nu}(\xi_0)U_{A,B}. \tag{2.13}

Here t2iν=e2iνlog⁡tt^{2i\nu}=e^{2i\nu\log t}; without this demodulation there need not be a limit. These conclusions follow by evaluating the stated homogeneity at t2ξ0t^2\xi_0 in (2.4). A different smooth low-frequency extension changes (1.3) by an inverse transform of a smooth compact function, hence by a smooth function by differentiation under its compact integral. The next theorem proves that such a change has no effect.

3. Covectors where the distribution is regular

We use the Fourier definition: (x0,ξ0)(x_0,\xi_0), ξ0≠0\xi_0\ne0, is regular for vv if some compact smooth χ\chi, equal to one near x0x_0, has χv^\widehat{\chi v} decreasing faster than every power in an open cone about ξ0\xi_0. No coordinate-invariance theorem about wavefront sets is assumed in the proof that follows.

Lemma T3 (the compact-distribution Fourier bound). A compactly supported distribution has a smooth Fourier transform of polynomial growth, and for every test ff

v(f)=(2π)−n∫v^(ξ)f^(−ξ) dξ.(3.1) v(f)=(2\pi)^{-n}\int\widehat v(\xi)\widehat f(-\xi)\,d\xi. \tag{3.1}

Proof. Choose a smooth compact cutoff χ=1\chi=1 near the support and define v^(ξ)=v(χ(x)e−ix⋅ξ)\widehat v(\xi)=v(\chi(x)e^{-ix\cdot\xi}). Independence of χ\chi follows from the definition of support. The fixed-support finite-order bound gives ∣v^(ξ)∣≤C⟨ξ⟩M|\widehat v(\xi)|\leq C\langle\xi\rangle^M. Difference quotients of the exponential converge with all derivatives on supp⁡χ\operatorname{supp}\chi, by the integral Taylor formula; applying that bound proves smoothness and all differentiated formulas. Invert ff by Q4, multiply by χ\chi, and apply vv. To justify moving vv inside the integral, truncate frequency to a cube and approximate its integral by Riemann sums in the finitely many derivative supremum norms used by vv. Uniform continuity on the two compact sets gives this convergence. The derivative tails are bounded by C∫∣ξ∣>R⟨ξ⟩M∣f^(−ξ)∣ dξ→0C\int_{|\xi|>R}\langle\xi\rangle^M|\widehat f(-\xi)|\,d\xi\to0, by Q3. Linearity and the finite-order estimate now give (3.1). □\square

Theorem T4 (rapid regular zoom). If (x0,ξ0)(x_0,\xi_0) is regular for vv, then for every real KK and real smooth ψ\psi with the same value and first derivative as in Section 2,

tK(ve−it2ψ)(x0+x/t)=OD′(t−L)for every L>0.(3.2) t^K(v e^{-it^2\psi})(x_0+x/t)=O_{\mathcal D'}(t^{-L}) \quad\hbox{for every }L>0. \tag{3.2}

Proof. Fixed compact tests use only a neighborhood on which χ=1\chi=1 for large tt, so replace vv by χv\chi v. For linear ψ(y)=ξ0⋅(y−x0)\psi(y)=\xi_0\cdot(y-x_0), (3.1) and the affine substitution give the pairing

(2π)−ntK+n∫eix0⋅(t2ξ0+tη)χv^(t2ξ0+tη)ϕ^(−η) dη.(3.3) (2\pi)^{-n}t^{K+n}\int e^{ix_0\cdot(t^2\xi_0+t\eta)} \widehat{\chi v}(t^2\xi_0+t\eta)\widehat\phi(-\eta)\,d\eta. \tag{3.3}

For some c>0c>0, ∣η∣<ct|\eta|<ct puts the Fourier argument in the regular cone with norm at least ct2ct^2. Arbitrarily high cone decay makes this part smaller than any specified power of tt, uniformly in finitely many test seminorms. On ∣η∣≥ct|\eta|\geq ct, Lemma T3 bounds the Fourier factor by Ct2M⟨η⟩MCt^{2M}\langle\eta\rangle^M. Q3 bounds ϕ^\widehat\phi by any inverse power with a finite test seminorm. Choose that power greater than K+n+2M+L+M+n+1K+n+2M+L+M+n+1, increasing it to a positive integer if necessary. The tail in (3.3) is O(t−L)O(t^{-L}). For general ψ\psi, replace the test by e−iqtϕe^{-iq_t}\phi. Formula (2.7) bounds every required seminorm of these tests uniformly, so the same estimates apply. This proves the asserted uniform result. □\square

A smooth error is regular at each nonzero covector: its localized smooth compact representative has rapid Fourier decay by Q3. More generally, any remainder satisfying the specified microlocal regularity is harmless in T1 for every real rr, by taking K=−(2r+n)K=-(2r+n) in T4.

4. Nonlinear coordinates, frames and phases

Lemma T5 (tangent pullback). Let κ\kappa be a smooth local diffeomorphism, κ(0)=0\kappa(0)=0, T=κ′(0)T=\kappa'(0). Put κt(x)=tκ(x/t)\kappa_t(x)=t\kappa(x/t). If vt→vv_t\to v weakly as t→∞t\to\infty, then, locally on every fixed compact test support,

κt∗vt⟶T∗v.(4.1) \kappa_t^*v_t\longrightarrow T^*v. \tag{4.1}

This holds both for scalar distributions and for distributional half-densities. Scalar pullback sends a test ϕ\phi to ϕ∘κt−1∣det⁡Dκt−1∣\phi\circ\kappa_t^{-1}|\det D\kappa_t^{-1}|; half-density pullback uses the power 1/21/2 instead of one.

Proof. The segment identity κt(x)=∫01κ′(sx/t)x ds\kappa_t(x)=\int_0^1\kappa'(sx/t)x\,ds proves convergence to TxTx in every derivative on bounded sets. The inverse map is explicitly κt−1(y)=tκ−1(y/t)\kappa_t^{-1}(y)=t\kappa^{-1}(y/t) wherever needed, so the same proof gives convergence to T−1yT^{-1}y. For a fixed test support KK, the images κt(K)\kappa_t(K) lie in one compact set K′K' and in a common domain of these inverses for large tt. This follows directly from the uniform convergence to TT and the expanding scaled domain of the original local inverse. Choose K′K' slightly larger than those images. Extending the transformed tests by zero is smooth, and their supports stay in the interior of K′K'.

The determinant is a polynomial in the derivative entries. Its absolute value is bounded away from zero near the compact sets in question, because the limiting matrix is invertible. Thus the absolute determinant and its positive square root, with all derivatives, converge to the constant factors for T−1T^{-1}. Chain and product rules show convergence of the transformed tests in every derivative on K′K'.

To apply a weakly convergent distribution to these moving tests, take any sequence tj→∞t_j\to\infty. The selected programme theorem U1 gives one finite-order bound for vtjv_{t_j} on K′K', and its moving-test conclusion gives (4.1) along that sequence. If a scalar test pairing failed to converge as t→∞t\to\infty, its failure would supply a sequence tj≥jt_j\geq j bounded away from the proposed limit, a contradiction. Thus (4.1) holds for the parameter itself. This argument does not assume uniform boundedness over unexamined finite intervals of tt. □\square

The same proof gives an O(t−1)O(t^{-1}) difference between κt∗vt\kappa_t^*v_t and T∗vtT^*v_t whenever vtv_t already has a uniform local finite-order bound: the transformed tests differ by O(t−1)O(t^{-1}) in every derivative, using one more segment Taylor formula. T1 supplies such a bound even when ctc_t does not converge. Thus the ordinary-symbol comparison itself, as well as a convergent classical limit, transforms by the tangent linear map.

For a smooth finite-dimensional frame change GG, G(x0+x/t)−G(x0)=O(t−1)G(x_0+x/t)-G(x_0)=O(t^{-1}) in every compact derivative norm. The same finite-order estimate therefore replaces the frame by its value at x0x_0. For half-densities the pullback exponent above follows also from the ordinary-function rule f(κ(x))∣det⁡Dκ(x)∣1/2f(\kappa(x))|\det D\kappa(x)|^{1/2}: substituting variables in the pairing leaves ∣det⁡Dκ−1∣1/2|\det D\kappa^{-1}|^{1/2}. This establishes the convention without assuming a transformation formula for distributions.

In particular, if the original object is u(y)∣dy∣1/2u(y)|dy|^{1/2}, its normalized half-density zoom is t−2mt^{-2m} times the pullback of e−it2ψu(y)∣dy∣1/2e^{-it^2\psi}u(y)|dy|^{1/2} by y=x0+x/ty=x_0+x/t. The pullback contributes t−n/2t^{-n/2}, so its coefficient in ∣dx∣1/2|dx|^{1/2} is exactly WtW_t, since 2m=2r+n/22m=2r+n/2. Under a nonlinear chart change, the two scaled charts differ by κt\kappa_t; T5 therefore gives the tangent half-density rule with the same normalization.

If p(x0)=0p(x_0)=0, p′(x0)=0p'(x_0)=0, and C=p′′(x0)C=p''(x_0), then (2.7) gives t2p(x0+x/t)→xTCx/2t^2p(x_0+x/t)\to x^TCx/2 in every compact derivative norm. Hence changing ψ\psi to ψ+p\psi+p gives

UA,B+C=e−ixTCx/2UA,B.(4.2) U_{A,B+C}=e^{-ix^TCx/2}U_{A,B}. \tag{4.2}

The exact equality also follows directly from the defining test integral. These changes compose by addition of symmetric matrices. Every such matrix is realized by the quadratic function p(y)=(y−x0)TC(y−x0)/2p(y)=(y-x_0)^TC(y-x_0)/2, so the action on phase second derivatives is transitive.

5. The singular quadratic distribution

Let k=rank⁡Ak=\operatorname{rank}A, R=ran⁡AR=\operatorname{ran}A, Z=ker⁡AZ=\ker A, and use the orthogonal split x=(y,z)∈R⊕Zx=(y,z)\in R\oplus Z. Symmetry gives R=Z⊥R=Z^\perp: AηA\eta is orthogonal to each kernel vector, and the dimensions agree by the rank-nullity proof in the earlier linear-algebra chain. The restriction AR:R→RA_R:R\to R is invertible; write BRB_R for the restriction of the quadratic form BB to RR.

Theorem T6. In the orthonormal coordinates just specified,

UA,B(y,z)=(2π)n/4−k/2∣det⁡AR∣−1/2e−iπsgn⁡AR/4ei2yT(AR−1−BR)y δ0(z).(5.1) U_{A,B}(y,z)= (2\pi)^{n/4-k/2}|\det A_R|^{-1/2} e^{-i\pi\operatorname{sgn}A_R/4} e^{\frac i2 y^T(A_R^{-1}-B_R)y}\,\delta_0(z). \tag{5.1}

For k=0k=0, the empty determinant is one and the signature is zero; the formula is (2π)n/4δ0(x)(2\pi)^{n/4}\delta_0(x).

Proof. Diagonalize ARA_R by the proved Q5. Insert e−ε∣η∣2/2e^{-\varepsilon|\eta|^2/2} in the defining integral. On testing, the limit as ε↓0\varepsilon\downarrow0 follows from the integrable Schwartz bound already used in (2.2), by Q1. Each nonzero eigenvalue aa contributes the elementary Gaussian

(2π)1/2(ε+ia)−1/2exp⁡ ⁣(−y22(ε+ia)),(5.2) (2\pi)^{1/2}(\varepsilon+ia)^{-1/2} \exp\!\left(-\frac{y^2}{2(\varepsilon+ia)}\right), \tag{5.2}

where the square root has positive real part. Q2 gives its limit (2π)1/2∣a∣−1/2e−iπsgn⁡(a)/4eiy2/(2a)(2\pi)^{1/2}|a|^{-1/2}e^{-i\pi\operatorname{sgn}(a)/4}e^{iy^2/(2a)}. The factors are uniformly bounded in absolute value for small positive ε\varepsilon, including their prefactors. Each kernel direction contributes (2π/ε)1/2e−z2/(2ε)(2\pi/\varepsilon)^{1/2}e^{-z^2/(2\varepsilon)}. Its integral is 2π2\pi, and its mass outside any fixed neighborhood of zero tends to zero after the substitution z=εwz=\sqrt\varepsilon w, by the Gaussian tail estimate from Q2. Thus in d=n−kd=n-k such directions the product tends to (2π)dδ0(z)(2\pi)^d\delta_0(z).

For completeness, these limits can be combined against a Schwartz test: the kernel Gaussian has uniformly bounded total mass; on a large compact yy-set the other factors converge uniformly and the test is uniformly continuous in zz; outside that yy-set, a sufficiently high power of ⟨y⟩−1\langle y\rangle^{-1} bounds the test uniformly in zz. The uniform bound on the oscillatory factors then makes the remaining tail arbitrarily small. This is the compact-plus-tail argument Q1 and does not require an interchange of two unevaluated oscillatory integrals. Finally multiply by e−ixTBx/2e^{-ix^TBx/2}. At z=0z=0 this is precisely e−iyTBRy/2e^{-iy^TB_Ry/2}; cross terms vanish. Combining powers gives −3n/4+k/2+(n−k)=n/4−k/2-3n/4+k/2+(n-k)=n/4-k/2, proving (5.1). □\square

Multiplication by this chirp preserves the Schwartz test space: each derivative of the chirp is the chirp times a polynomial, by induction with the product rule. Every weighted derivative of its product with a Schwartz function is therefore bounded by finitely many Schwartz seminorms. This also proves continuity of the tempered distribution defined in (5.1).

Differentiating Euler's identity in ξ\xi gives H′′(ξ)ξ=0H''(\xi)\xi=0. Consequently the radial direction ξ0\xi_0 is always in the kernel in this homogeneous graph problem. Assuming the full matrix AA invertible would discard the actual conic case.

6. The two tangent planes

In the cotangent coordinates (x,ξ)(x,\xi), use ω=∑dξj∧dxj\omega=\sum d\xi_j\wedge dx_j, so ω((v,w),(v′,w′))=w⋅v′−w′⋅v\omega((v,w),(v',w'))=w\cdot v'-w'\cdot v. The graph ξ↦(H′(ξ),ξ)\xi\mapsto(H'(\xi),\xi) has tangent plane {(Aη,η)}\{(A\eta,\eta)\}. The graph of dψd\psi has tangent plane {(v,Bv)}\{(v,Bv)\}. Both have dimension nn; their symplectic pairings vanish by symmetry of AA or BB. This verifies the Lagrangian property directly, where a Lagrangian plane means an nn-dimensional subspace with zero restricted symplectic form.

Subtracting the phase gradient applies the linear shear (v,w)↦(v,w−Bv)(v,w)\mapsto(v,w-Bv). It preserves ω\omega, since its additional terms cancel by symmetry of BB, and sends the second plane to {(v,0)}\{(v,0)\}. The first becomes

λA,B={(Aη,η−BAη):η∈Rn}.(6.1) \lambda_{A,B}=\{(A\eta,\eta-BA\eta):\eta\in\mathbb R^n\}. \tag{6.1}

Its parametrization is injective: if both components vanish, the second is η\eta. It therefore still has dimension nn. The generating quadratic QA,BQ_{A,B} yields exactly this plane, because ∂ηQ=x−Aη=0\partial_\eta Q=x-A\eta=0 and ∂xQ=η−Bx\partial_x Q=\eta-Bx. These calculations establish all assertions about the planes without importing a general normal-form theorem.

6.1. The intrinsic symbol and Gaussian line

Put B=ΩΛ1/2⊗L⊗π∗E\mathscr B=\Omega_\Lambda^{1/2}\otimes\mathscr L\otimes\pi^*E, where L\mathscr L is the relative Maslov line, and let s∈Sm+n/4(Λ;B)s\in S^{m+n/4}(\Lambda;\mathscr B) represent the principal-symbol class. The companion Z2–Z4 constructs the canonical map Jρ\mathcal J_\rho to families of distributional half-densities on Tx0XT_{x_0}X, indexed by phase second jets. In a frequency chart it sends β∣dξ∣1/2\beta|d\xi|^{1/2}, evaluated in the geometric Maslov line at the horizontal transversal, to βUA,B∣dv∣1/2\beta U_{A,B}|dv|^{1/2}.

Theorem T8 (intrinsic tangent comparison). Write δR(x,ξ)=(x,Rξ)\delta_R(x,\xi)=(x,R\xi) and at(v)=x0+v/ta_t(v)=x_0+v/t in a local chart. Then

t−2mat∗(e−it2ψu)=Jρ ⁣(t−2m−n/2(δt2∗s)ρ)ψ+OD′(t−1).(6.2) t^{-2m}a_t^*(e^{-it^2\psi}u) =\mathcal J_\rho\!\left( t^{-2m-n/2}(\delta_{t^2}^*s)_\rho \right)_\psi+O_{\mathcal D'}(t^{-1}). \tag{6.2}

The complete proof is Z5: the intrinsic graph representative gives T1, T4 removes its microlocally regular error, and the half-density radial factor is tnt^n, giving precisely ct=t−2rb(t2ξ0)c_t=t^{-2r}b(t^2\xi_0). A lower-order symbol changes the right side by only OD′(t−2)O_{\mathcal D'}(t^{-2}). Z4 and T5 prove coordinate and frame independence, including the second-derivative term of a nonlinear chart. For a leading homogeneous symbol, (6.2) converges to its fibre value under Jρ\mathcal J_\rho. Ordinary symbols need not have that limit. Z6 checks dilation and all signature factors; Z7 gives an exact nonlinear chart example with a nontrivial fourth-root transition.

7. Exercises and complete solutions

The unit sphere flattens under the stated zoom; its model has delta support in the radial direction and a unit-modulus oscillating coefficient in the transverse direction.

Figure 1. The exact two-dimensional case c=ρ=1c=\rho=1, B=0B=0, in Exercise 2. The first panel is the base projection y=(cos⁡θ,sin⁡θ)y=(\cos\theta, \sin\theta); it does not depict the full cotangent graph. The second plots x=t(cos⁡(s/t)−1,sin⁡(s/t))x=t(\cos(s/t)-1,\sin(s/t)), ∣s∣≤2.5|s|\leq2.5, for t=2,8t=2,8, with limiting support x1=0x_1=0. The last panel shows the real and imaginary parts of the coefficient of δ0(x1)\delta_0(x_1), not a pointwise graph of a delta distribution. Its magnitude is one. T1 and T6 prove the zoom and its exact normalization; Section 6 gives the corresponding cotangent plane. The free generating-function comparison is Guillemin–Sternberg §5.15.4. Reproducible figure source.

Exercise 1 — the three scales. Suppose the carrier frequency is of size tat^a, the spatial zoom is t−bt^{-b}, and the frequency window is tct^c, with b,c>0b,c>0. Find the scales for which the mixed term and both generic quadratic terms survive at order one.

Solution. The mixed term scales as tc−bt^{c-b}. A homogeneous degree-one frequency Hessian at frequency taξ0t^a\xi_0 has size t−at^{-a}, so the frequency quadratic scales as t2c−at^{2c-a}. The modulating phase quadratic scales as ta−2bt^{a-2b}. Setting all three exponents to zero gives a=2b=2ca=2b=2c. Conversely those equalities make all three factors one. This concerns the generic quadratic model; a particular vanishing Hessian does not force its absent term to survive.

Exercise 2 — a sphere. For H(ξ)=c∣ξ∣H(\xi)=c|\xi|, ξ0=ρω\xi_0=\rho\omega, ρ>0\rho>0, ∣ω∣=1|\omega|=1, compute (5.1), including n=1n=1.

Solution. Differentiating the positive square root gives x0=cωx_0=c\omega and A=(c/ρ)(I−ωωT)A=(c/\rho)(I-\omega\omega^T). For c≠0c\ne0, R=ω⊥R=\omega^\perp, Z=RωZ=\mathbb R\omega, and k=n−1k=n-1. Writing x=y+zωx=y+z\omega, the answer is

(2π)(2−n)/4(ρ/∣c∣)(n−1)/2e−iπ(n−1)sgn⁡(c)/4ei2yT((ρ/c)I−BR)yδ0(z).(7.1) (2\pi)^{(2-n)/4}(\rho/|c|)^{(n-1)/2} e^{-i\pi(n-1)\operatorname{sgn}(c)/4} e^{\frac i2y^T((\rho/c)I-B_R)y}\delta_0(z). \tag{7.1}

For n=1n=1 all transverse matrices are empty, giving (2π)1/4δ0(x)(2\pi)^{1/4}\delta_0(x). For c=0c=0, the rank is zero in every dimension and (5.1) gives (2π)n/4δ0(x)(2\pi)^{n/4}\delta_0(x).

Exercise 3 — no fixed limit for an ordinary symbol. In one dimension take H=0H=0, x0=0x_0=0, ξ0=1\xi_0=1, ψ(x)=x\psi(x)=x, and a smooth symbol equal to ξr(2+sin⁡log⁡ξ)\xi^r(2+\sin\log\xi) for ξ≥2\xi\geq2 and zero for ξ≤1\xi\leq1. Prove that the zoom does not converge. Preserve the example after cutting off the distribution near zero.

Solution. Repeated differentiation of the expression on ξ≥2\xi\geq2 gives ξr−j\xi^{r-j} times a fixed linear combination of one, sine and cosine. On the compact transition region all derivatives are bounded. Thus b∈Srb\in S^r, while ct=2+sin⁡(2log⁡t)c_t=2+\sin(2\log t) for large tt. Here U0,0=(2π)1/4δ0U_{0,0}=(2\pi)^{1/4}\delta_0. Along tj=eπj+π/4t_j=e^{\pi j+\pi/4} the coefficient is three, and along sj=eπj+3π/4s_j=e^{\pi j+3\pi/4} it is one. T1 and a test with value one at zero give distinct limits.

Let χ\chi be compact smooth and equal to one near zero. The reduced symbol of χu\chi u, obtained by testing (1.3) with χ(x)e−ixξ\chi(x)e^{-ix\xi}, is the absolutely convergent convolution

bχ(ξ)=(2π)−1∫χ^(ζ)b(ξ−ζ) dζ.(7.2) b_\chi(\xi)=(2\pi)^{-1}\int\widehat\chi(\zeta)b(\xi-\zeta)\,d\zeta. \tag{7.2}

Fourier inversion gives (2π)−1∫χ^=χ(0)=1(2\pi)^{-1}\int\widehat\chi=\chi(0)=1. For every jj, differentiate (7.2) under its integrable Schwartz majorant and subtract b(j)(ξ)b^{(j)}(\xi). On ∣ζ∣≤∣ξ∣/2|\zeta|\leq|\xi|/2, ∣ξ∣≥2|\xi|\geq2, the segment derivative estimate bounds the difference by C∣ζ∣⟨ξ⟩r−j−1C|\zeta|\langle\xi\rangle^{r-j-1}. On its complement, ∣b(j)(ξ−ζ)∣+∣b(j)(ξ)∣|b^{(j)}(\xi-\zeta)|+|b^{(j)}(\xi)| is bounded by a fixed polynomial in ⟨ξ⟩\langle\xi\rangle and ⟨ζ⟩\langle\zeta\rangle; arbitrary Schwartz decay of χ^\widehat\chi, together with ∣ζ∣>∣ξ∣/2|\zeta|>|\xi|/2, gives the same desired bound, even when r−j−1<0r-j-1<0. Bounded ξ\xi is handled by the common integrable majorant. Hence bχ−b∈Sr−1b_\chi-b\in S^{r-1}, including every derivative. Its contribution to ctc_t is O(t−2)O(t^{-2}), so the two limits persist for this compactly supported example. No general unproved symbol composition rule is used.

Exercise 4 — a nonlinear half-density pullback. Near zero let κ(x)=ax+dx2\kappa(x)=ax+dx^2, a≠0a\ne0. Compute the pullback of δ0∣dy∣1/2\delta_0|dy|^{1/2}, and compare its tangent zoom.

Solution. Restrict to a neighborhood on which κ\kappa is a diffeomorphism. Its inverse derivative at zero is 1/a1/a, so the test formula in T5 gives κ∗(δ0∣dy∣1/2)=∣a∣−1/2δ0∣dx∣1/2\kappa^*(\delta_0|dy|^{1/2})=|a|^{-1/2}\delta_0|dx|^{1/2}. The scalar factor would instead be ∣a∣−1|a|^{-1}. Now κt(x)=ax+dx2/t\kappa_t(x)=ax+dx^2/t has the same derivative at zero for every tt; its half-density pullback of the delta is exactly the same distribution. In particular its tangent limit agrees, for either sign of aa, without an orientation sign.

Exercise 5 — conjugation. Determine the model of u‾\overline u at the opposite covector and compare (5.1).

Solution. Complex conjugation of a distribution means u‾(ϕ)=u(ϕ‾)‾\overline u(\phi)=\overline{u(\overline\phi)}. Conjugate the absolutely convergent test pairing (1.3) and change ξ\xi to −ξ-\xi. Then H~(ξ)=−H(−ξ)\widetilde H(\xi)=-H(-\xi), b~(ξ)=b(−ξ)‾\widetilde b(\xi)=\overline{b(-\xi)}, and the chosen covector is −ξ0-\xi_0, with the same base point. Its Hessian is −A-A. Taking phase −ψ-\psi gives −B-B, and

UA,B‾=U−A,−B.(7.3) \overline{U_{A,B}}=U_{-A,-B}. \tag{7.3}

This follows either by η↦−η\eta\mapsto-\eta in the test integral or directly from (5.1): the absolute determinant is unchanged and the signature and quadratic exponent change sign. Half-density Jacobian factors are real positive, so they commute with conjugation. A complex bundle frame changes to its conjugate frame. Thus the local Gaussian transition factors also conjugate. The global identification is supplied separately by Z4 of the Gaussian-symbol companion; the local conjugation calculation agrees with it.

Exercise 6 — the first error can be t−1t^{-1}. In one dimension take u=Dqδ0u=D^q\delta_0, q≥0q\geq0, and ψ(x)=ξ0x\psi(x)=\xi_0x, ξ0≠0\xi_0\ne0. Compute the normalized zoom exactly.

Solution. Its Fourier transform is ξq\xi^q, so b=(2π)−1/4ξqb=(2\pi)^{-1/4}\xi^q, r=qr=q, m=q+1/4m=q+1/4. Leibniz's rule inside the delta derivative pairing gives

e−it2ξ0xDqδ0=∑k=0q(qk)(t2ξ0)q−kDkδ0. e^{-it^2\xi_0x}D^q\delta_0 =\sum_{k=0}^q\binom qk(t^2\xi_0)^{q-k}D^k\delta_0.

Indeed the q−kq-k derivatives falling on the exponential cancel the corresponding factors of ii in the distributional derivative pairing. The affine test formula gives (Dkδ0)(x/t)=tk+1Dkδ0(x)(D^k\delta_0)(x/t)=t^{k+1}D^k\delta_0(x). Consequently

Wt=∑k=0q(qk)ξ0q−kt−kDkδ0.(7.4) W_t=\sum_{k=0}^q\binom qk\xi_0^{q-k}t^{-k}D^k\delta_0. \tag{7.4}

The leading term is ξ0qδ0\xi_0^q\delta_0, exactly T2 and (5.1). For q≥1q\geq1, the first error is qξ0q−1t−1Dδ0q\xi_0^{q-1}t^{-1}D\delta_0, which is nonzero: test with a compact smooth function whose derivative at zero is one. The higher terms are O(t−2)O(t^{-2}) on that test. Thus a universal O(t−2)O(t^{-2}) error for the complete zoom is false even for a classical symbol. For q=0q=0, the formula is exact with zero error.

8. Free readings and programme closure

The completed intrinsic graph, localization, phase and principal-symbol proofs are linked in the Gaussian-symbol companion. Its Z1–Z6 prove the remaining graph-to-zoom and global bundle connection. No paid work supplies any argument in this receiving draft.

Original receiving exposition and exercises: GPT-6 Astra (OpenAI), Ultra, 4 October 2026, CC0. The separate earlier programme components retain their own CC BY-SA 4.0, CC0 or GFDL 1.2-only notices.