Tangent zooms and quadratic models
To inspect a singularity at a particular covector, remove a matching oscillation and magnify a shrinking neighborhood of its base point. For a Lagrangian distribution, the resulting model is quadratic. Its Hessians describe two tangent Lagrangian planes, and its coefficient is the leading amplitude at the selected frequency. Ordinary symbols need not have a limiting leading coefficient, so the general result is an asymptotic comparison rather than an assertion that every zoom converges.
The exact earlier programme proofs are:
- The restored intrinsic-regularity lesson, especially its frequency-graph criterion, and the restored oscillatory-distribution lesson, including its exact localization and Fourier reduction.
- Quadratic stationary phase, Q1–Q6: the compact-plus-tail limit rule, Gaussian factors, Schwartz estimates, Fourier inversion, symmetric diagonalization and the regularized quadratic identity. The full stationary-phase proof supplies the nonstationary integration-by-parts bound with its transposed vector field.
- Wavefront tests and pullbacks, T0–W5, together with the earlier tangent companion's compact-distribution Fourier estimate T3 and rapid regular zoom T4.
- Uniform local distribution bounds and moving tests, U1. Its complete-test-space component, Sections 6, 14.1–14.2 and 19, proves the Baire and completeness inputs at the declared logical base. That separate component retains its GFDL 1.2 licence and notices.
The exact proof map records every dependency. The primary sources for this restoration are the reprints of Hörmander IV, Section 25.1 and Hörmander III, Section 21.6. They supply the mathematical results; this lesson supplies its own exposition, exercises and full programme proofs. The earlier tangent companion remains available, including its uniform OD′(t−1) estimate, Gaussian-symbol identification and reproducible figures.
We keep D=−i∂, the Fourier phase −x⋅ξ, and the inverse coefficient (2π)−n. Half densities are written in the local frame ∣dx∣1/2.
1. Two scales select one tangent problem
Suppose the Lagrangian has the frequency graph
Λ={(H′(ξ),ξ)},
where H is real, smooth and homogeneous of degree one away from zero. Extend it smoothly through low frequencies when writing a Fourier integral. For a compactly supported u∈Im(Rn,Λ), define the normalized symbol
b(ξ)=(2π)−n/4eiH(ξ)u(ξ),b∈Sr,r=m−n/4.(1.1)
Thus, up to a smooth low-frequency term,
u(x)=(2π)−3n/4∫ei(x⋅ξ−H(ξ))b(ξ)dξ.(1.2)
The normalization in (1.1) differs by a constant from the reduced symbol in the preceding lesson. It is chosen so that the coefficient of the quadratic model will later be a half density on Λ.
Fix γ0=(x0,ξ0)∈Λ, ξ0=0. Choose a real smooth function ψ near x0 with
ψ(x0)=0,ψ′(x0)=ξ0.(1.3)
For t→+∞, examine
Wt(x)=t−2m−n/2(ue−it2ψ)(x0+x/t).(1.4)
This defines a distribution on every fixed compact set of the new x variables for sufficiently large t. The factor t−n/2 is the half-density Jacobian of the map x↦x0+x/t; the remaining t−2m measures order at frequency t2ξ0. The spatial scale is t−1, and the frequency deviations from that center have scale t.
Put
A=H′′(ξ0),B=ψ′′(x0),QA,B(x,η)=x⋅η−21xTBx−21ηTAη.(1.5)
Define the tempered quadratic distribution
UA,B(x)=(2π)−3n/4∫eiQA,B(x,η)dη.(1.6)
This is (2π)n/4 times the inverse Fourier transform of e−iηTAη/2, multiplied by e−ixTBx/2. Both operations are defined on tempered distributions: derivatives of each quadratic exponential are polynomially bounded, so multiplication preserves Schwartz space and extends by transposition.
2. The asymptotic comparison and its tail estimate
Theorem 2.1 (tangent zoom). With the preceding hypotheses,
Wt−t−2rb(t2ξ0)UA,B⟶0in D′(Rn).(2.1)
The scalar coefficient is bounded, but need not converge. The result holds for finite-dimensional vector-valued amplitudes componentwise. A smooth localized error in u contributes a term tending to zero faster than every inverse power of t.
Proof. Fix a compactly supported smooth test function f. In (1.2), substitute ξ=t2ξ0+tη. The measure contributes tn, so (1.4) paired with f is
(2π)−3n/4∫bt(η)Ft(η)dη,bt(η)=t−2rb(t2ξ0+tη),(2.2)
where
Ft(η)Et(x,η)=∫eiEt(x,η)f(x)dx,=(x0+x/t)⋅(t2ξ0+tη)−t2ψ(x0+x/t)−H(t2ξ0+tη).(2.3)
These identities can first be taken with frequency cutoffs. The estimates below make their test pairings absolutely convergent and justify removing those cutoffs.
Euler's identity and the graph condition give H(ξ0)=x0⋅ξ0 and H′(ξ0)=x0. For bounded η, homogeneity and Taylor expansion give
H(t2ξ0+tη)=t2H(ξ0)+tx0⋅η+21ηTAη+O(t−1).
The analogous expansion of t2ψ(x0+x/t) has linear term tξ0⋅x and quadratic term xTBx/2. All constant and linear terms cancel in (2.3), leaving
Et(x,η)⟶QA,B(x,η)
uniformly on compact sets, with the needed fixed x derivatives. Consequently Ft→F on compact η sets, where F(η)=∫eiQA,Bfdx.
The tails are uniformly controlled. On the fixed support of f,
∂xEt=tξ0+η−tψ′(x0+x/t)=η+O(1).(2.4)
Every x derivative of order at least two is bounded uniformly: at order k its phase contribution is −t2−kψ(k)(x0+x/t). For large ∣η∣, (2.4) therefore has size at least ∣η∣/2. The full transposed nonstationary vector field from the stationary-phase lesson, including its divergence, gives
∣Ft(η)∣+∣F(η)∣≤CN⟨η⟩−Nfor every N,(2.5)
uniformly in large t. Each bound uses only finitely many derivatives of f.
The symbol variation can now be frozen at the central frequency. In ∣η∣<ct, with c<∣ξ0∣/2, the segment from t2ξ0 to t2ξ0+tη remains in a fixed high-frequency annulus after scaling by t2. The derivative estimate for b gives
∣bt(η)−bt(0)∣≤Ct−1∣η∣.(2.6)
For the complementary region ∣η∣≥ct, we only need a uniform polynomial bound. There t≤C⟨η⟩. If r≥0, the growth bound for b gives ∣bt(η)∣≤C⟨η⟩2r. If r<0, boundedness of b at all frequencies gives ∣bt(η)∣≤C⟨η⟩−2r. Thus
∣bt(η)∣≤C⟨η⟩2∣r∣,∣bt(0)∣≤C.(2.7)
On the small region, (2.5)–(2.6) bound the integral of (bt−bt(0))Ft by C/t. On the complementary region, (2.5)–(2.7), with N>2∣r∣+n+1, give a tail tending to zero. Hence
∫(bt(η)−bt(0))Ft(η)dη⟶0.
Compact convergence and (2.5) also give ∫Ftdη→∫Fdη. Multiplying by the bounded bt(0) proves (2.1).
Finally, a compactly localized smooth distribution has a rapidly decreasing Fourier transform. In the same integral, its transform decreases rapidly near t2ξ0, while (2.4) controls the complementary frequencies. The more general wavefront-regular argument in the next section proves the asserted rapid decay, including the fixed polynomial factor in (1.4). ∎
If b is classical with leading term br homogeneous of the real degree r, then
t−2rb(t2ξ0)⟶br(ξ0),Wt⟶br(ξ0)UA,B.(2.8)
The remainder of symbol order r−1 contributes O(t−2) to this coefficient. A complex homogeneous degree r+iν has the additional factor t2iν; without removing that factor, (2.8) is not a claimed limit when ν=0.
3. A regular covector disappears under every zoom
Proposition 3.1. Let u∈D′ near x0, with (x0,ξ0)∈/WF(u), ξ0=0. For ψ satisfying (1.3), and every real K,
tK(ue−it2ψ)(x0+x/t)⟶0in D′.(3.1)
Proof. Insert a compact cutoff equal to one near x0, small enough that the compact distribution's Fourier transform is rapidly decreasing in a cone about ξ0. On every fixed zoomed test support, this replacement is exact for large t.
First take ψ0(y)=ξ0⋅(y−x0). The Fourier transform in the zoomed x variable is
tK+neix0⋅(t2ξ0+tη)u(t2ξ0+tη).(3.2)
For ∣η∣<ct, with c sufficiently small, the argument stays in the good cone and has size comparable to t2. Its rapid decrease beats the fixed factor tK+n, uniformly on every bounded η set. A compact distribution has a global polynomial Fourier bound. For ∣η∣≥ct, the inequalities t≤C⟨η⟩ and ∣t2ξ0+tη∣≤C⟨η⟩2 therefore bound (3.2) by one fixed polynomial in η, independently of t. The inner region is bounded as well, by choosing enough decay in the good cone. Dominated convergence against Schwartz functions gives convergence to zero in S′.
For the actual ψ, write p=ψ−ψ0. Both its value and first derivative vanish at x0, so
t2p(x0+x/t)⟶21xTψ′′(x0)xin C∞ on compact sets.
The additional exponential multipliers and all their derivatives consequently converge on each compact test support. They preserve the conclusion just obtained. Explicitly, along any sequence tj→∞, the convergent distributions in (3.2) have a common finite-order bound there; applying that bound to the difference between the varying multiplied test and its limit, then using weak convergence on the fixed limit test, proves (3.1). Since this works along every such sequence, it proves the full limit. ∎
This result makes the tangent model microlocal. Changes in a distribution away from the selected wavefront direction cannot change the model, even if they are singular in other directions.
4. Nonlinear coordinates become linear in the limit
Lemma 4.1 (rescaling a change of coordinates). Let Mε(x)=εx, and suppose distributions uε, defined near zero, satisfy
Mε∗uε⟶Vin D′,ε↓0.
If κ is a smooth local diffeomorphism, κ(0)=0, and T=κ′(0), then
Mε∗κ∗uε⟶T∗V.(4.1)
The statement holds for scalar distributions and for half-density distributions, with their respective pullback conventions.
Proof. Put vε=Mε∗uε. The exact conjugated map is
κε=Mε−1∘κ∘Mε,κε(x)=ε−1κ(εx).
Thus the left side of (4.1) is κε∗vε. Taylor's formula and differentiation give κε→T in C∞ on every compact set. The inverse maps are ε−1κ−1(εy), which likewise converge to T−1.
For a fixed compact test support, the transposed pullbacks of that test have their supports in one fixed compact set. For scalars, the transposed test is
f(κε−1(y))∣detDκε−1(y)∣.
For half densities its Jacobian exponent is 1/2. In either case these tests converge in Cc∞ to the corresponding test for T. Along an arbitrary sequence εj↓0, uniform boundedness of the convergent vεj supplies a common finite-order estimate. It bounds their pairing with the difference of the varying test and its limit by a quantity tending to zero. Their pairing with the fixed limit test tends to that of V. This proves (4.1) along every sequence, and hence the result. ∎
After centering charts at x0, the lemma says that an existing tangent limit is a distribution on Tx0X. For half densities, it transforms with the half-density Jacobian of the tangent linear map. Changing a smooth bundle frame contributes only its value at x0, by the same varying-test argument, so vector coefficients belong to the fiber Ex0.
The function ψ is also a choice. If it is replaced by ψ+p, where p(x0)=p′(x0)=0, then the tangent limit is multiplied by
e−ixTp′′(x0)x/2.(4.2)
Equivalently B in the quadratic model becomes B+p′′(x0). These chirp identifications are transitive because their exponents add. They are the elementary change between the linear reference Lagrangian planes η=Bx.
5. Evaluate the model even when the Hessian is singular
Let A be any real symmetric matrix, with rank k. Split orthogonally
Rn=ranA⊕kerA,x=(y,z),AR=A∣ranA.
The restriction AR is invertible, and let BR be the restriction of the quadratic form B to ranA. Exact Fourier inversion in the kernel variables gives a point mass δ0(z). The Fresnel identity in the other k variables gives
UA,B(y,z)=(2π)n/4−k/2∣detAR∣−1/2e−iπsgnAR/4×eiyT(AR−1−BR)y/2δ0(z).(5.1)
When k=0, the empty determinant is one, the empty signature is zero, and (5.1) is (2π)n/4δ0(x). The sign in the Fresnel factor is negative because the frequency Hessian of the phase is −AR. Terms in B containing z vanish on multiplication with δ0(z), so only BR remains.
For the tangent model of a conic graph, Aξ0=0: differentiate Euler's identity for H. Thus its frequency Hessian is always singular at a nonzero covector. Formula (5.1) covers precisely that possibility, and the model is a nonzero Gaussian simple layer on ranA.
The phase in (1.6) has critical equation x=Aη and output covector η−Bx. Its linear Lagrangian is therefore
λA,B={(Aη,η−BAη):η∈Rn}.(5.2)
The tangent plane of the original graph is {(Aη,η)}. Removing the Hessian of ψ applies the canonical shear (x,η)↦(x,η−Bx). This explains why the model is determined by those two tangent planes rather than by higher derivatives of H or ψ.
6. Exercises with complete solutions
Exercise 6.1 (the two scales; introductory). Replace the substitution in the zoom proof by ξ=taξ0+tbη and the spatial scale by x0+t−cx. Determine the relations among a,b,c that keep both the cross term x⋅η and the quadratic frequency term at order one.
Solution. The cross term has factor tb−c, so b=c. By degree-one homogeneity the frequency Hessian at taξ0 is t−aH′′(ξ0), so its quadratic deviation has factor t2b−a; thus a=2b=2c. The phase-removal frequency is ta, and its spatial quadratic term has factor ta−2c=1 as well. Taking c=1 gives the scales t2,t,t−1 used above. If the particular Hessian vanishes, its term imposes no constraint, but these relations give the common nondegenerate quadratic scaling across all graphs.
Exercise 6.2 (a sphere tangent model; intermediate). Take H(ξ)=c∣ξ∣, c=0, ξ0=se1, s>0, and ψ(x)=ξ0⋅(x−ce1). Compute UA,0.
Solution. Here x0=ce1 and
A=sc(I−e1e1T).
It has kernel Re1, rank n−1, and signature (n−1)sgnc. Writing x=(x1,x′), formula (5.1) gives
UA,0=(2π)n/4−(n−1)/2∣c/s∣−(n−1)/2e−iπ(n−1)sgnc/4eis∣x′∣2/(2c)δ0(x1).
The supporting hyperplane is the tangent plane of the sphere ∣x∣=∣c∣ at ce1, in the translated tangent variables. For n=1 the empty-factor convention gives (2π)1/4δ0.
Exercise 6.3 (an ordinary symbol with no zoom limit; intermediate). In dimension one take H=0, ξ0=1, and, for positive high frequencies,
b(ξ)=ξr(2+sin(logξ)),
with a smooth cutoff that removes low and negative frequencies. Show that this is an ordinary symbol and that a compact localization of its inverse phase integral can have no limit in (1.4).
Solution. Each frequency derivative is bounded by Ckξr−k, because differentiating logξ or the power costs one inverse frequency. Thus b∈Sr, with m=r+1/4. The normalized central coefficient is t−2rb(t2)=2+sin(2logt). A compact base cutoff equal to one near zero changes the reduced symbol only by Sr−1, by Fourier reduction for xξ; that change contributes O(t−2) to the coefficient. Take ψ(x)=x. Here A=B=0 and U=(2π)1/4δ0. Along tj=exp(πj+π/4) the coefficient tends to three, and along sj=exp(πj+3π/4) it tends to one. Theorem 2.1 gives these two distinct distributional limits. A test nonzero at zero distinguishes them. The general asymptotic comparison is therefore strictly more general than (2.8).
Exercise 6.4 (a nonlinear chart; intermediate). In dimension one let κ(x)=ax+dx2, a=0, near zero. If a rescaled half-density family has limit Cδ0, compute its limit after changing coordinates by κ.
Solution. The conjugated map is ε−1κ(εx)=ax+dεx2, converging to multiplication by a. The scalar pullback of δ0 is ∣a∣−1δ0. Half-density pullback contributes the Jacobian factor ∣a∣1/2, so the resulting coefficient is C∣a∣−1/2δ0 in the new half-density frame. The quadratic chart coefficient d does not enter the tangent pullback. This concerns the full rescaled family; its phase-removal function must be transformed with the same coordinate change.
Exercise 6.5 (complex conjugation; advanced). Show directly that UA,B=U−A,−B. If the Fourier phase of u is H, find the graph phase and central covector for u.
Solution. Conjugate (1.6) and replace η by −η. The cross term returns to x⋅η, and both quadratic matrices change sign. This proves the model identity, including the opposite Fresnel signature. Since u(ξ)=u(−ξ), the new graph function is H(ξ)=−H(−ξ), and its base gradient is H′(−ξ). The selected point is therefore (x0,−ξ0), and the new normalized amplitude is b(−ξ). The Hessian at −ξ0 is −A; taking phase removal −ψ gives −B. Thus the conjugated zoom has exactly the conjugated quadratic model. Patching this statement globally will also require conjugating the Maslov transition factors.
Exercise 6.6 (an exact distributional zoom; advanced). On R, take u=Dqδ0, with q a nonnegative integer, H=0, and ψ(y)=ξ0y, ξ0=0. Compute (1.4) exactly and identify its limit.
Solution. The Fourier transform of u is ξq, so its intrinsic order is m=q+1/4. Fourier transformation of the rescaled and modulated distribution gives
Wt(η)=t−2q(tη+t2ξ0)q=(ξ0+η/t)q.
Thus the exact formula is
Wt=k=0∑q(kq)ξ0q−kt−kDkδ0.
The limit is ξ0qδ0. In Theorem 2.1, the normalized coefficient is (2π)−1/4ξ0q, and U0,0=(2π)1/4δ0, giving exactly that limit. For q≥1, the first correction is qξ0q−1t−1Dδ0; its order shows why a universal O(t−2) error is not supplied by the classical coefficient's O(t−2) remainder alone.
References
- [Hörmander III, §21.6] Lars Hörmander, The Analysis of Linear Partial Differential Operators III: Pseudo-Differential Operators, corrected second printing, Springer, 1994, §21.6, especially.
- [Hörmander IV, §25.1] Lars Hörmander, The Analysis of Linear Partial Differential Operators IV: Fourier Integral Operators, corrected second printing, Springer, 1994, §25.1, especially Proposition 25.1.7, its wavefront remark and Lemma 25.1.8.
Written by GPT-6.1 Sol (OpenAI), Ultra, September 2026. Restoration and exact prerequisite review: GPT-6 Astra (OpenAI), Ultra, 5 October 2026. Self-checked by the writing AI. Original text: public domain (CC0).