Hamilton fields and subprincipal transport
When the principal symbol of a scalar differential or pseudodifferential operator vanishes along a canonical relation, its leading product with an FIO vanishes. The next symbol contains a derivative along that relation. Because the FIO symbol is a half density, that derivative includes half the divergence of the Hamilton field. The remaining scalar coefficient is the subprincipal symbol.
The exact earlier programme proofs are:
- The restored phase-space lesson fixes the Hamilton sign and proves the conic Lagrangian identities and symplectic reductions. Its finite-coordinate flow companion, Sections 17.1–17.7, proves existence, uniqueness, smooth parameters, the local flow law and coordinate compatibility.
- The restored intrinsic-regularity lesson, Section 3, gives the full frequency normal form; the Gaussian-symbol lesson, Sections 4–7, gives flat Maslov frames, actual density degrees and the exact principal-symbol quotient.
- The restored kernel and composition lesson supplies the smooth mapping, scalar kernel correspondence, proper localization and ordinary composition proofs. Its kernel companion specifies the distribution topologies.
- The included subprincipal coordinate companion retains AN03-U012's complete scalar calculation G11–G14, with its licence and notices. Section S1 proves the first coordinate correction, including every differentiated remainder, using the existing amplitude and coordinate calculus. Thus the invariant coefficient below is an exact programme result.
- Homogeneous symbol estimates, H2–H3 and H6–H8, prove all derivative counts, support-preserving asymptotic sums and cutoff convergence. Density integration, product measure and exponentials and circle periods supply the elementary inputs.
The proof map binds every use to the exact retained source and proof. The primary source is the reprint of Hörmander IV, corrected second printing (1994), Section 25.2 equation 25.2.11, Theorem 25.2.4 and Lemma 25.2.5. The independently written lesson, its calculations and all six original solved exercises are retained. All symbol estimates in this lesson are ordinary S1,0 estimates. The scalar operator is classical with step-one homogeneous expansion. The FIO need not have a homogeneous leading amplitude. All constructions are local over compact base sets, with microlocal cutoffs inside the stated open cones.
1. A vector field differentiates a density as well as its coefficient
For a complex number κ, the density line Ωκ(M) has transition factors Jκ=eκlogJ, with J>0. Their multiplicativity defines the line even for complex κ. In a coordinate chart a section is f∣du∣κ.
Let V=∑jVj∂uj be real, with local flow Ft. Define
LVa=dtdFt∗at=0.(1.1)
In coordinates,
LV(f∣du∣κ)=(Vf+κj∑∂ujVjf)∣du∣κ.(1.2)
Proof. Pullback gives
Ft∗(f∣du∣κ)=f(Ft(u))∣detdFt(u)∣κ∣du∣κ.
At t=0, dF0=I and ∂tdFt∣0=dV. The determinant derivative is trdV; its sign is positive near t=0. Differentiating the positive power gives κtrdV, and differentiating f(Ft(u)) gives Vf. This proves (1.2). Pullback is coordinate independent, so the expression patches on overlaps. ∎
For a complex vector field V=V1+iV2, with real V1,V2, define LV=LV1+iLV2. Formula (1.2) is complex linear in the coefficients and their derivatives, so it remains valid and coordinate independent. No real flow of a complex vector field is asserted.
A Maslov line has locally constant transition factors. In such local frames, differentiate its density coefficient by (1.2). Constant transition factors commute with the derivative, so this defines LV on Maslov-valued densities. For half densities the coefficient of divergence is 1/2.
For a real field and compactly supported smooth half densities, integration by parts also gives
∫⟨LVa,b⟩+∫⟨a,LVb⟩=0.(1.3)
Indeed the integrand coefficient is V(fg)+(divV)fg, the divergence of Vfg. Its integral is zero in each chart; a partition of unity proves the global identity. More explicitly, a partition ∑ℓχℓ=1 is finite near the compact supports. Apply the chart identity to χℓa,b, then sum. The additional coefficients ∑ℓVχℓ vanish, proving the claimed identity with no residual chart-boundary term. Thus the half-density correction is exactly what makes the derivative skew under the density pairing.
2. Vanishing on a Lagrangian makes the Hamilton field tangent
Let C be a homogeneous canonical relation from T∗Y∖0 to T∗X∖0, with the full punctured-cotangent closure of the preceding lesson, and put Λ=C′. Let p be a smooth scalar function on T∗X∖0, homogeneous of degree m, whose restriction to the X projection of C is zero. On T∗(X×Y), lift p by making it independent of the Y variables. Its Hamilton field is
Hp=j=1∑nX(pξj∂xj−pxj∂ξj).(2.1)
Lemma 2.1. The lifted field is tangent to Λ. If p is complex, it takes values in the complexification of TΛ.
Proof. The restriction p∣Λ is zero, so dp(v)=0 for every v∈TΛ. With ω=∑dξ∧dx+∑dη∧dy, our convention is ω(Hp,v)=−dp(v). Therefore Hp∈(TΛ)ω=TΛ, since Λ is Lagrangian. Apply the same argument to the real and imaginary parts when needed. The input-covector reflection transporting Λ to C preserves this lifted field. ∎
The Lie derivative of the symbol is therefore defined along the relation itself. Its dependence on a local coefficient includes the density term from Section 1.
3. The transport theorem and its exact hypotheses
Let P be a properly supported classical scalar pseudodifferential operator of order m on half densities on X. In a coordinate half-density frame, write its left symbol as
P(x,ξ)=p(x,ξ)+r(x,ξ)(modSm−2),(3.1)
where p and r are homogeneous of degrees m and m−1 away from zero. Its intrinsic scalar subprincipal symbol is the imported quantity
c=r+2ij∑pxjξj.(3.2)
The half-density hypothesis is part of the coordinate-invariance statement.
Suppose p vanishes on the X projection of C, and
A∈Im′(X×Y,C′;Ω1/2).
Write a for its principal symbol, of intrinsic symbol order
ν=m′+4nX+nY(3.3)
in MΛ⊗ΩΛ1/2, modulo one lower order.
Theorem 3.1 (scalar subprincipal transport). The kernel of PA belongs to
Im+m′−1(X×Y,C′;Ω1/2),(3.4)
with principal symbol
σ(PA)=i1LHpa+ca,(3.5)
of intrinsic order ν+m−1, modulo one lower order.
The proper support of P makes this operator product meaningful even when A is not properly supported: A maps compact smooth inputs to smooth outputs, and P acts on all smooth functions. The kernel correspondence then defines PA. The theorem asserts its Lagrangian class; proper support of PA additionally follows if A is proper. The proof is microlocal and applies inside any fixed working cone. It concerns scalar half-density operators; bundle systems require their corresponding first-order frame or connection data.
We prove the theorem after establishing a coordinate choice that keeps the two manifolds separate.
4. Almost every diagonal graph complements a Lagrangian plane
Write T∗RN=RxN×RξN, and set
Lb={(x,ξ):ξj=bjxj, 1≤j≤N},b∈RN.(4.1)
These planes are Lagrangian, since their graph matrices are symmetric.
Lemma 4.1 (diagonal complements). For any Lagrangian plane L, the planes L and Lb are transverse for almost every b, with respect to Lebesgue measure.
Proof of existence. Use induction on N, beginning with the zero-dimensional space. In the first symplectic coordinate pair choose b1 so that
v=e1+b1f1∈/L.
Such a choice exists: two distinct slopes in L would put both e1 and f1 in L, contradicting isotropy. Put V=Rv. Its symplectic orthogonal is Vω={ξ1=b1x1}. The quotient Vω/V is the symplectic space of the remaining coordinate pairs. Because v∈/L=Lω, the equation ω(v,w)=0 is a nonzero linear condition on L. Thus L∩Vω has dimension N−1, intersects V trivially, and projects to a Lagrangian plane L in this quotient.
Choose b2,…,bN by induction so that L is transverse to the corresponding diagonal graph. A vector in L∩Lb then has zero image in the quotient, so lies in V∩L=0. Hence L∩Lb=0, proving existence.
Proof of almost-everywhere transversality. Write a basis of L as the columns of (WU). A vector in its intersection with Lb has coefficient vector in the kernel of
W−diag(b)U.
Consequently nontransversality is the zero set of the polynomial
q(b)=det(W−diag(b)U).(4.2)
Existence proves that this polynomial is not identically zero. A one-variable nonzero polynomial of degree d has at most d roots: for a root c, each difference tk−ck is (t−c)∑j=0k−1tk−1−jcj; dividing off this factor reduces the degree, and induction proves the bound. A finite set has measure zero by covering its points with intervals of arbitrarily small total length. A nonzero real polynomial in several variables has a null zero set: by induction on the number of variables, outside the null common zero set of a nonzero coefficient polynomial its one-variable slices have finitely many roots. Fubini on each bounded box, followed by their countable union, proves the assertion. This proves the lemma. ∎
The measure assertion is stronger than merely finding one complement. It follows from the determinant polynomial, rather than from an unspecified genericity assumption.
5. Separate base coordinates give a full-frequency generating function
Fix (x0,ξ0;y0,η0)∈Λ; here η0 is the kernel covector, so the input relation covector is −η0. Both ξ0 and η0 are nonzero. Translate the base coordinates to make x0=y0=0. Lemma 4.1 supplies a diagonal graph complementary to TΛ in the full cotangent tangent space. Its diagonal matrix splits into an X block BX and a Y block BY.
We can realize that graph as the old-coordinate image of the new horizontal tangent plane using separate base changes. Choose a component ξ0,ℓ=0, and define an old base coordinate map x=F(u) with
Fj(u)=uj(j=ℓ),Fℓ(u)=uℓ−2ξ0,ℓuTBXu.(5.1)
It fixes zero and has derivative identity there, so is a local diffeomorphism. Under its cotangent lift the new covector is (dF)Tξ. At the fixed point,
δξnew=δξold−BXδu.(5.2)
Thus a new horizontal variation maps to δξold=BXδu. Use a nonzero component of η0 to construct the same change on Y. These independent maps realize the full diagonal graph. Transversality means that, in the new coordinates, projection of Λ to all its frequency variables is locally invertible.
Write those variables as θ=(ξ,η). The relation becomes a graph
(x,y)=h(θ),(5.3)
where h is homogeneous of degree zero. The conic Lagrangian identity α∣Λ=0, proved in the phase-space lesson, says θ⋅dh=0. Therefore
H(θ)=θ⋅h(θ),dH=h⋅dθ.(5.4)
This is a degree-one generating function, and
x=Hξ,y=Hη.(5.5)
The phase
Φ=x⋅ξ+y⋅η−H(ξ,η)(5.6)
is nondegenerate: the critical equations (x,y)−Hθ=0 have independent differentials because their derivative in (x,y) is identity. On the selected cone its full differential is nonzero, since θ=0.
Let n=nX+nY. The Fourier normal form of the intrinsic lesson gives, microlocally modulo a smooth kernel,
KA(x,y)=(2π)−3n/4∬eiΦb(ξ,η)dξdη,b∈Sm′−n/4.(5.7)
Its symbol is b∣dξdη∣1/2 in the local Maslov unit supplied by this frequency graph. This unit is a flat local frame. The normalization in (5.7) is the general (2π)−(n+2N)/4 with N=n.
One can also remove base dependence directly. For a phase amplitude b0(x,y,θ), Taylor's integral formula about (x,y)=Hθ expresses its difference from the restriction to the critical graph as ∑j((x,y)j−Hθj)bj. In the oscillatory integral, each such factor becomes −Dθjbj, one lower ordinary symbol order. Derivatives of Hθ cost one inverse frequency radius. Repeat the reduction and use the support-preserving asymptotic sum proved in H6. The compact base and angular cutoffs can be chosen inside the original working patch; H6 retains that support. Every finite remainder has the corresponding lower order, and the resulting frequency-only amplitude leaves a smooth remainder. This explains the same normal form without treating a base cutoff as a frequency-independent amplitude. All these normal forms concern an interior cone; compact base cutoffs equal to one near the point are retained for operator estimates.
6. Integration by parts produces the lower product order
We now prove Theorem 3.1 in the coordinates of Section 5. For fixed compact external X- and Y-sets, properness of P confines the intermediate X points to a compact set. Cut A's two base variables off outside neighborhoods of these compact sets. That localized kernel is proper, so the ordinary composition theorem applies. The canonical relation of P is the identity graph: its local kernel phase is (x−x′)⋅ξ, with exactly the ordinary pseudodifferential normalization. Thus composition keeps the relation C, and the cutoffs give the same product on the selected external sets. Smooth kernel errors remain smooth after composition with this localized A, by the smooth-input mapping theorem and its adjoint version. Properness and the exact conic localization proof therefore allow the calculation: inserting base cutoffs equal to one near the point changes the calculation only by a microlocally smooth kernel. Acting on the x plane wave gives its left symbol P(x,ξ), so the product kernel is represented by (5.7) with amplitude P(x,ξ)b(ξ,η).
On the critical graph, p(Hξ,ξ)=0. Extend the local symbols over short base segments if necessary. Taylor's integral formula writes
p(x,ξ)=j∑(xj−Hξj)pj(x,ξ,η),(6.1)
pj=∫01pxj(Hξ+t(x−Hξ),ξ)dt.(6.2)
Each pj is homogeneous of degree m in the joint frequencies on the working cone. Here ∣ξ∣≍∣θ∣, since the X covector is nonzero on its compact angular support. All derivatives have the ordinary symbol bounds: positive frequency derivatives of Hξ have the corresponding inverse-radius loss.
Since
(xj−Hξj)eiΦ=i1∂ξjeiΦ,
integration by parts replaces the principal amplitude by
−j∑Dξj(pjb)=ij∑∂ξj(pjb).(6.3)
These derivatives hold x,y fixed. The resulting amplitude has order m+m′−n/4−1. The rb term has that same order, and the Sm−2 remainder times b has one lower order. Frequency cutoffs and oscillatory-integral continuity justify (6.3); at infinity their derivative errors vanish in the distribution limit by integration by parts. The phase theorem now gives (3.4).
To find its leading symbol, restrict this amplitude to the critical graph. Principal symbols depend on that restriction; every additional base reduction contributes another inverse frequency order. Put
Fj(ξ,η)=pxj(Hξ,ξ)=pj(Hξ,ξ,η).(6.4)
The coefficient of the new symbol is
ij∑Fj∂ξjb+(r+ij∑(∂ξjpj)x=Hξ)b,(6.5)
where r is evaluated at (Hξ,ξ). We next identify its scalar correction invariantly.
7. Half the divergence leaves exactly the subprincipal coefficient
On Λ, the lifted Hamilton field in the frequency coordinates is
W=−j∑Fj∂ξj,Wη=0.(7.1)
Indeed (2.1) gives its frequency components, and Lemma 2.1 gives tangency. More explicitly, differentiating p(Hξ,ξ)=0 in ξ and η shows that the induced derivatives of Hξ,Hη along (7.1) are respectively pξ and zero. These are exactly its base components.
Consequently
i1LW(b∣dξdη∣1/2)=(ij∑Fj∂ξjb+2ij∑∂ξjFjb)∣dξdη∣1/2.(7.2)
The derivative of Fj in (7.2) is a total derivative along the critical graph. This differs from the fixed-base derivative in (6.5).
To compare them, abbreviate, with all terms evaluated on the graph,
S=j∑∂ξjpj,T=j,k∑(∂xkpj)Hξkξj.(7.3)
The chain rule gives
j∑∂ξjFj=S+T.(7.4)
On the other hand, differentiate (6.1) first in xk, then in ξk at fixed x, and restrict to the graph. The factor xj−Hξj then vanishes, leaving
pxkξk=∂ξkpk−j∑(∂xkpj)Hξjξk.
Summing and using the symmetry of the Hessian of H yields
k∑pxkξk=S−T.(7.5)
Subtract (7.2) from (6.5). The remaining multiplication coefficient is
r+iS−2i(S+T)=r+2i(S−T)=r+2ik∑pxkξk=c.(7.6)
This proves (3.5), including its sign and density correction.
The calculation patches: c is an intrinsic scalar by the exact AN-03 half-density result, the Hamilton field is intrinsic, and the Lie derivative uses the flat Maslov transitions and intrinsic density pullback. Thus it gives the same symbol in every working patch. Altering the input symbol by one lower order alters (6.5) by one lower output order, so the formula is well defined on principal-symbol quotient classes.
Finally cotangent dilation pulls the tangent Hamilton field back by the factor tm−1. One can see this directly in (7.1): its coefficients have degree m, while a frequency derivative costs one degree. Lie differentiation therefore raises intrinsic symbol order by m−1. The coefficient c has the same degree m−1. Applied to (3.3), this gives precisely m+m′−1+n/4, the symbol order of (3.4). For the derivative estimates, the coefficients Fj have degree m in frequency coordinates, and their α-derivatives have degree m−∣α∣. The product rule applied to Fj∂ξjb and (∂ξjFj)b therefore gives order m−1+ord(b), with every derivative and finitely many input seminorms. The frequency half-density frame has degree n/2, so this gives the asserted intrinsic order as well. H2–H3 transport these estimates between homogeneous charts; compact base/direction localization and the proved locally finite partitions patch them. The Taylor integrals in Section 6 have the same finite-seminorm bounds. Applying the same argument to a one-order-lower input proves the asserted quotient independence. This completes the theorem. ∎
8. A lifting model shows the density term directly
Take X=Rs×St1 and Y=Rz, and use the frames ∣dsdt∣1/2, ∣dz∣1/2. For a smooth periodic coefficient b(t), define
(Abg)(s,t)=b(t)g(s).(8.1)
Its kernel is b(t)δ(s−z), of order −1/4. The relation has s=z, the same nonzero s,z covector χ, and zero t covector. Coordinates on its kernel Lagrangian are (s,t,χ). The invariant symbol is a fixed normalized Maslov factor times
b(t)∣dsdtdχ∣1/2.(8.2)
The fixed factor is (2π)1/4, from the codimension-one delta normalization, and does not affect the following differential identity.
For real smooth periodic v(t) and a smooth scalar q(t), let
P=21(vDt+Dtv)+q=vDt−2iv′+q.(8.3)
This properly supported differential operator acts on half densities. Its principal symbol is p=v(t)ξt, which vanishes on the relation, and
c=−2iv′+q+2i∂t∂ξtp=q.(8.4)
On the kernel Lagrangian, Hp=v∂t; its divergence in (s,t,χ) is v′. Formula (3.5) gives the coefficient
−i(vb′+21v′b)+qb.(8.5)
Applying (8.3) to (8.1) gives exactly the same coefficient, with no discarded terms. Its kernel order is again −1/4=1−1/4−1. This checks the transport law in an actual properly supported operator model satisfying the nonzero-covector hypotheses.
On an interval where v>0, vanishing of this coefficient is the first-order equation
b′+2vv′b+viqb=0.
Its solutions are
b(t)=Kv(t)−1/2exp(−i∫t0tv(u)q(u)du),(8.6)
with arbitrary complex constant K. Differentiate to verify the equation; conversely multiplying it by v1/2exp(i∫q/v) gives derivative zero. This is a scalar transport solution, not a claim of a global parametrix construction.
9. Exercises with complete solutions
Exercise 9.1 (flow and density weight; introductory). On R2, let V=u∂u−2w∂w. Compute its flow and its Lie derivative on a half density f∣dudw∣1/2. Check the result by direct pullback.
Solution. The flow is Ft(u,w)=(etu,e−2tw), whose positive determinant is e−t. Its divergence is 1−2=−1. Formula (1.2) gives
LV(f∣dudw∣1/2)=(ufu−2wfw−21f)∣dudw∣1/2.
Direct pullback has coefficient e−t/2f(etu,e−2tw). Its derivative at zero is the same displayed expression. Differentiating only f would omit the determinant contribution.
Exercise 9.2 (the exceptional slopes; intermediate). In T∗R2, let L be the graph of B=(0110). Find the slopes (b1,b2) for which L fails to be transverse to the diagonal graph Lb. Find one complement and the dimension of the intersection at an exceptional slope.
Solution. With U=I,W=B, the determinant is
det(B−diag(b1,b2))=b1b2−1.
The exceptional set is the hyperbola b1b2=1, a null subset of the plane. The choice (0,0) is transverse. At an exceptional slope, neither bj is zero, the displayed matrix has determinant zero and rank one, and its kernel has dimension one. Thus the Lagrangian intersection is one-dimensional there. The failure condition concerns the two slopes together.
Exercise 9.3 (realize two separate coordinate jets; intermediate). At a kernel covector with ξ0=2, η0=−3 on R×R, construct separate base maps, fixing zero with derivative one, which carry their new horizontal tangent plane to the old diagonal graph of slopes (5,7). Explain the role of the two nonzero covectors.
Solution. Take
x=F(u)=u−45u2,y=G(w)=w+67w2.
Both maps are locally invertible at zero. Contracting their second derivatives with the fixed covectors gives 2F′′=−5, −3G′′=−7. Therefore
δξnew=δξold−5δu,δηnew=δηold−7δw.
Zero new frequency variations are precisely the old slopes (5,7). A zero fixed covector would contract every second derivative to zero; a base change with derivative identity could then not create an arbitrary slope in that block by this construction. This is why separate changes here use both nonzero components.
Exercise 9.4 (a sign check with a circle mode; intermediate). In the lifting model take b(t)=eiℓt, ℓ∈Z, and P=Dt. Compute PAb, its order and its transport symbol. Include ℓ=0.
Solution. The principal symbol ξt vanishes on the relation, the subprincipal symbol is zero, and the tangent Hamilton field is ∂t, with zero divergence. Hence the transport coefficient is (1/i)b′=ℓb. Direct differentiation gives Dteiℓt=ℓeiℓt, so PAb=ℓAb. For ℓ=0 its exact order is −1/4, as predicted by 1−1/4−1; its leading symbol is nonzero. For ℓ=0, the product is zero and belongs to every lower order. Declared upper order need not be the exact nonzero order.
Exercise 9.5 (a left coefficient is not the subprincipal coefficient; advanced). In (8.3) set v(t)=2+sint, q=0. Compute the left order-zero coefficient, the subprincipal symbol and the action on the lift. Find a nonzero periodic b killed by this product.
Solution. The left order-zero coefficient is −icost/2, while the mixed principal derivative adds +icost/2. Thus the subprincipal symbol is zero. The action coefficient is
−i((2+sint)b′+21costb).
Take b=(2+sint)−1/2. It is smooth, positive and periodic, and
b′=−21cost(2+sint)−3/2, so the two terms cancel exactly. Keeping only vb′, or treating the left order-zero coefficient as the invariant subprincipal symbol, would give an incorrect transport formula.
Exercise 9.6 (periodic transport has an obstruction; advanced). Suppose v>0 and q are real smooth 2π-periodic functions in (8.3). Determine when PAb=0 has a nonzero smooth periodic solution. Give examples with and without such a solution.
Solution. Formula (8.6) is the complete local solution, and extends smoothly along the line because v>0. Over one period the positive factor v−1/2 returns to its original value, while the phase is multiplied by
exp(−i∫02πv(t)q(t)dt).
A nonzero solution is periodic exactly when this factor is one, equivalently when the real integral lies in 2πZ. If it does, all its derivatives are periodic as well, by the smooth periodic differential equation, or by differentiating the explicit formula. If it does not, the only periodic solution is zero. For v=2+sint and q=v, the solution b=v−1/2e−it is periodic. For the same v and q=v/2, its period multiplier is −1, so no nonzero periodic solution exists. The Maslov unit in this lifting model is globally fixed; this calculation asserts this model's scalar obstruction and does not omit an additional arbitrary Maslov holonomy.
References
- [Hörmander IV, §25.2] Lars Hörmander, The Analysis of Linear Partial Differential Operators IV: Fourier Integral Operators, Springer, reprint of the corrected second printing (1994), equation 25.2.11, Theorem 25.2.4 and Lemma 25.2.5.
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).