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 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,0S_{1,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 κ\kappa, the density line Ωκ(M)\Omega^\kappa(M) has transition factors Jκ=eκlog⁡JJ^\kappa=e^{\kappa\log J}, with J>0J>0. Their multiplicativity defines the line even for complex κ\kappa. In a coordinate chart a section is f∣du∣κf|du|^\kappa.

Let V=∑jVj∂ujV=\sum_jV_j\partial_{u_j} be real, with local flow FtF_t. Define

LVa=ddtFt∗a∣t=0.(1.1) \mathcal L_V a=\left.\frac{d}{dt}F_t^*a\right|_{t=0}. \tag{1.1}

In coordinates,

LV(f∣du∣κ)=(Vf+κ∑j∂ujVj f)∣du∣κ.(1.2) \mathcal L_V(f|du|^\kappa) =\left(Vf+\kappa\sum_j\partial_{u_j}V_j\,f\right)|du|^\kappa. \tag{1.2}

Proof. Pullback gives

Ft∗(f∣du∣κ)=f(Ft(u))∣det⁡dFt(u)∣κ∣du∣κ. F_t^*(f|du|^\kappa) =f(F_t(u))|\det dF_t(u)|^\kappa |du|^\kappa.

At t=0t=0, dF0=IdF_0=I and ∂tdFt∣0=dV\partial_t dF_t|_0=dV. The determinant derivative is tr⁡dV\operatorname{tr}dV; its sign is positive near t=0t=0. Differentiating the positive power gives κtr⁡dV\kappa\operatorname{tr}dV, and differentiating f(Ft(u))f(F_t(u)) gives VfVf. This proves (1.2). Pullback is coordinate independent, so the expression patches on overlaps. ∎

For a complex vector field V=V1+iV2V=V_1+iV_2, with real V1,V2V_1,V_2, define LV=LV1+iLV2\mathcal L_V=\mathcal L_{V_1}+i\mathcal L_{V_2}. 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\mathcal L_V on Maslov-valued densities. For half densities the coefficient of divergence is 1/21/2.

For a real field and compactly supported smooth half densities, integration by parts also gives

∫⟨LVa,b⟩+∫⟨a,LVb⟩=0.(1.3) \int\langle\mathcal L_Va,b\rangle +\int\langle a,\mathcal L_Vb\rangle=0. \tag{1.3}

Indeed the integrand coefficient is V(fg‾)+(div⁡V)fg‾V(f\overline g)+(\operatorname{div}V)f\overline g, the divergence of Vfg‾Vf\overline g. Its integral is zero in each chart; a partition of unity proves the global identity. More explicitly, a partition ∑ℓχℓ=1\sum_\ell\chi_\ell=1 is finite near the compact supports. Apply the chart identity to χℓa,b\chi_\ell a,b, then sum. The additional coefficients ∑ℓVχℓ\sum_\ell V\chi_\ell 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 CC be a homogeneous canonical relation from T∗Y∖0T^*Y\setminus0 to T∗X∖0T^*X\setminus0, with the full punctured-cotangent closure of the preceding lesson, and put Λ=C′\Lambda=C'. Let pp be a smooth scalar function on T∗X∖0T^*X\setminus0, homogeneous of degree mm, whose restriction to the XX projection of CC is zero. On T∗(X×Y)T^*(X\times Y), lift pp by making it independent of the YY variables. Its Hamilton field is

Hp=∑j=1nX(pξj∂xj−pxj∂ξj).(2.1) H_p=\sum_{j=1}^{n_X} \left(p_{\xi_j}\partial_{x_j}-p_{x_j}\partial_{\xi_j}\right). \tag{2.1}

Lemma 2.1. The lifted field is tangent to Λ\Lambda. If pp is complex, it takes values in the complexification of TΛT\Lambda.

Proof. The restriction p∣Λp|_\Lambda is zero, so dp(v)=0dp(v)=0 for every v∈TΛv\in T\Lambda. With ω=∑dξ∧dx+∑dη∧dy\omega=\sum d\xi\wedge dx+\sum d\eta\wedge dy, our convention is ω(Hp,v)=−dp(v)\omega(H_p,v)=-dp(v). Therefore Hp∈(TΛ)ω=TΛH_p\in(T\Lambda)^\omega=T\Lambda, since Λ\Lambda is Lagrangian. Apply the same argument to the real and imaginary parts when needed. The input-covector reflection transporting Λ\Lambda to CC 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 PP be a properly supported classical scalar pseudodifferential operator of order mm on half densities on XX. In a coordinate half-density frame, write its left symbol as

P(x,ξ)=p(x,ξ)+r(x,ξ)(modSm−2),(3.1) P(x,\xi)=p(x,\xi)+r(x,\xi)\pmod{S^{m-2}}, \tag{3.1}

where pp and rr are homogeneous of degrees mm and m−1m-1 away from zero. Its intrinsic scalar subprincipal symbol is the imported quantity

c=r+i2∑jpxjξj.(3.2) c=r+\frac i2\sum_j p_{x_j\xi_j}. \tag{3.2}

The half-density hypothesis is part of the coordinate-invariance statement.

Suppose pp vanishes on the XX projection of CC, and

A∈Im′(X×Y,C′;Ω1/2). A\in I^{m'}(X\times Y,C';\Omega^{1/2}).

Write aa for its principal symbol, of intrinsic symbol order

ν=m′+nX+nY4(3.3) \nu=m'+\frac{n_X+n_Y}{4} \tag{3.3}

in MΛ⊗ΩΛ1/2M_\Lambda\otimes\Omega_\Lambda^{1/2}, modulo one lower order.

Theorem 3.1 (scalar subprincipal transport). The kernel of PAPA belongs to

Im+m′−1(X×Y,C′;Ω1/2),(3.4) I^{m+m'-1}(X\times Y,C';\Omega^{1/2}), \tag{3.4}

with principal symbol

σ(PA)=1iLHpa+c a,(3.5) \sigma(PA)=\frac1i\mathcal L_{H_p}a+c\,a, \tag{3.5}

of intrinsic order ν+m−1\nu+m-1, modulo one lower order.

The proper support of PP makes this operator product meaningful even when AA is not properly supported: AA maps compact smooth inputs to smooth outputs, and PP acts on all smooth functions. The kernel correspondence then defines PAPA. The theorem asserts its Lagrangian class; proper support of PAPA additionally follows if AA 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ξNT^*\mathbb R^N=\mathbb R^N_x\times\mathbb R^N_\xi, and set

Lb={(x,ξ):ξj=bjxj, 1≤j≤N},b∈RN.(4.1) L_b=\{(x,\xi):\xi_j=b_jx_j,\ 1\leq j\leq N\}, \qquad b\in\mathbb R^N. \tag{4.1}

These planes are Lagrangian, since their graph matrices are symmetric.

Lemma 4.1 (diagonal complements). For any Lagrangian plane LL, the planes LL and LbL_b are transverse for almost every bb, with respect to Lebesgue measure.

Proof of existence. Use induction on NN, beginning with the zero-dimensional space. In the first symplectic coordinate pair choose b1b_1 so that

v=e1+b1f1∉L. v=e_1+b_1f_1\notin L.

Such a choice exists: two distinct slopes in LL would put both e1e_1 and f1f_1 in LL, contradicting isotropy. Put V=RvV=\mathbb Rv. Its symplectic orthogonal is Vω={ξ1=b1x1}V^\omega=\{\xi_1=b_1x_1\}. The quotient Vω/VV^\omega/V is the symplectic space of the remaining coordinate pairs. Because v∉L=Lωv\notin L=L^\omega, the equation ω(v,w)=0\omega(v,w)=0 is a nonzero linear condition on LL. Thus L∩VωL\cap V^\omega has dimension N−1N-1, intersects VV trivially, and projects to a Lagrangian plane L‾\overline L in this quotient.

Choose b2,…,bNb_2,\ldots,b_N by induction so that L‾\overline L is transverse to the corresponding diagonal graph. A vector in L∩LbL\cap L_b then has zero image in the quotient, so lies in V∩L=0V\cap L=0. Hence L∩Lb=0L\cap L_b=0, proving existence.

Proof of almost-everywhere transversality. Write a basis of LL as the columns of (UW)\binom{U}{W}. A vector in its intersection with LbL_b has coefficient vector in the kernel of

W−diag⁡(b)U. W-\operatorname{diag}(b)U.

Consequently nontransversality is the zero set of the polynomial

q(b)=det⁡(W−diag⁡(b)U).(4.2) q(b)=\det\bigl(W-\operatorname{diag}(b)U\bigr). \tag{4.2}

Existence proves that this polynomial is not identically zero. A one-variable nonzero polynomial of degree dd has at most dd roots: for a root cc, each difference tk−ckt^k-c^k is (t−c)∑j=0k−1tk−1−jcj(t-c)\sum_{j=0}^{k-1}t^{k-1-j}c^j; 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)∈Λ(x_0,\xi_0;y_0,\eta_0)\in\Lambda; here η0\eta_0 is the kernel covector, so the input relation covector is −η0-\eta_0. Both ξ0\xi_0 and η0\eta_0 are nonzero. Translate the base coordinates to make x0=y0=0x_0=y_0=0. Lemma 4.1 supplies a diagonal graph complementary to TΛT\Lambda in the full cotangent tangent space. Its diagonal matrix splits into an XX block BXB_X and a YY block BYB_Y.

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\xi_{0,\ell}\neq0, and define an old base coordinate map x=F(u)x=F(u) with

Fj(u)=uj(j≠ℓ),Fℓ(u)=uℓ−uTBXu2ξ0,ℓ.(5.1) F_j(u)=u_j\quad(j\neq\ell),\qquad F_\ell(u)=u_\ell-\frac{u^TB_Xu}{2\xi_{0,\ell}}. \tag{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ξ(dF)^T\xi. At the fixed point,

δξnew=δξold−BXδu.(5.2) \delta\xi_{\rm new}=\delta\xi_{\rm old}-B_X\delta u. \tag{5.2}

Thus a new horizontal variation maps to δξold=BXδu\delta\xi_{\rm old}=B_X\delta u. Use a nonzero component of η0\eta_0 to construct the same change on YY. These independent maps realize the full diagonal graph. Transversality means that, in the new coordinates, projection of Λ\Lambda to all its frequency variables is locally invertible.

Write those variables as θ=(ξ,η)\theta=(\xi,\eta). The relation becomes a graph

(x,y)=h(θ),(5.3) (x,y)=h(\theta), \tag{5.3}

where hh is homogeneous of degree zero. The conic Lagrangian identity α∣Λ=0\alpha|_\Lambda=0, proved in the phase-space lesson, says θ⋅dh=0\theta\cdot dh=0. Therefore

H(θ)=θ⋅h(θ),dH=h⋅dθ.(5.4) H(\theta)=\theta\cdot h(\theta),\qquad dH=h\cdot d\theta. \tag{5.4}

This is a degree-one generating function, and

x=Hξ,y=Hη.(5.5) x=H_\xi,\qquad y=H_\eta. \tag{5.5}

The phase

Φ=x⋅ξ+y⋅η−H(ξ,η)(5.6) \Phi=x\cdot\xi+y\cdot\eta-H(\xi,\eta) \tag{5.6}

is nondegenerate: the critical equations (x,y)−Hθ=0(x,y)-H_\theta=0 have independent differentials because their derivative in (x,y)(x,y) is identity. On the selected cone its full differential is nonzero, since θ≠0\theta\neq0.

Let n=nX+nYn=n_X+n_Y. 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) K_A(x,y)=(2\pi)^{-3n/4} \iint e^{i\Phi}b(\xi,\eta)\,d\xi\,d\eta, \qquad b\in S^{m'-n/4}. \tag{5.7}

Its symbol is b∣dξ dη∣1/2b|d\xi\,d\eta|^{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(2\pi)^{-(n+2N)/4} with N=nN=n.

One can also remove base dependence directly. For a phase amplitude b0(x,y,θ)b_0(x,y,\theta), Taylor's integral formula about (x,y)=Hθ(x,y)=H_\theta expresses its difference from the restriction to the critical graph as ∑j((x,y)j−Hθj)bj\sum_j((x,y)_j-H_{\theta_j})b_j. In the oscillatory integral, each such factor becomes −Dθjbj-D_{\theta_j}b_j, one lower ordinary symbol order. Derivatives of HθH_\theta 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 XX- and YY-sets, properness of PP confines the intermediate XX points to a compact set. Cut AA'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 PP is the identity graph: its local kernel phase is (x−x′)⋅ξ(x-x')\cdot\xi, with exactly the ordinary pseudodifferential normalization. Thus composition keeps the relation CC, and the cutoffs give the same product on the selected external sets. Smooth kernel errors remain smooth after composition with this localized AA, 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 xx plane wave gives its left symbol P(x,ξ)P(x,\xi), so the product kernel is represented by (5.7) with amplitude P(x,ξ)b(ξ,η)P(x,\xi)b(\xi,\eta).

On the critical graph, p(Hξ,ξ)=0p(H_\xi,\xi)=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) p(x,\xi)=\sum_j(x_j-H_{\xi_j})p_j(x,\xi,\eta), \tag{6.1}
pj=∫01pxj(Hξ+t(x−Hξ),ξ) dt.(6.2) p_j=\int_0^1 p_{x_j}\bigl(H_\xi+t(x-H_\xi),\xi\bigr)\,dt. \tag{6.2}

Each pjp_j is homogeneous of degree mm in the joint frequencies on the working cone. Here ∣ξ∣≍∣θ∣|\xi|\asymp|\theta|, since the XX covector is nonzero on its compact angular support. All derivatives have the ordinary symbol bounds: positive frequency derivatives of HξH_\xi have the corresponding inverse-radius loss.

Since

(xj−Hξj)eiΦ=1i∂ξjeiΦ, (x_j-H_{\xi_j})e^{i\Phi} =\frac1i\partial_{\xi_j}e^{i\Phi},

integration by parts replaces the principal amplitude by

−∑jDξj(pjb)=i∑j∂ξj(pjb).(6.3) -\sum_jD_{\xi_j}(p_jb) =i\sum_j\partial_{\xi_j}(p_jb). \tag{6.3}

These derivatives hold x,yx,y fixed. The resulting amplitude has order m+m′−n/4−1m+m'-n/4-1. The rbr b term has that same order, and the Sm−2S^{m-2} remainder times bb 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) F_j(\xi,\eta)=p_{x_j}(H_\xi,\xi) =p_j(H_\xi,\xi,\eta). \tag{6.4}

The coefficient of the new symbol is

i∑jFj∂ξjb+(r+i∑j(∂ξjpj)∣x=Hξ)b,(6.5) i\sum_jF_j\partial_{\xi_j}b +\left(r+i\sum_j (\partial_{\xi_j}p_j)\big|_{x=H_\xi}\right)b, \tag{6.5}

where rr is evaluated at (Hξ,ξ)(H_\xi,\xi). We next identify its scalar correction invariantly.

7. Half the divergence leaves exactly the subprincipal coefficient

On Λ\Lambda, the lifted Hamilton field in the frequency coordinates is

W=−∑jFj∂ξj,Wη=0.(7.1) W=-\sum_jF_j\partial_{\xi_j},\qquad W\eta=0. \tag{7.1}

Indeed (2.1) gives its frequency components, and Lemma 2.1 gives tangency. More explicitly, differentiating p(Hξ,ξ)=0p(H_\xi,\xi)=0 in ξ\xi and η\eta shows that the induced derivatives of Hξ,HηH_\xi,H_\eta along (7.1) are respectively pξp_\xi and zero. These are exactly its base components.

Consequently

1iLW(b∣dξ dη∣1/2)=(i∑jFj∂ξjb+i2∑j∂ξjFj b)∣dξ dη∣1/2.(7.2) \frac1i\mathcal L_W\bigl(b|d\xi\,d\eta|^{1/2}\bigr) =\left(i\sum_jF_j\partial_{\xi_j}b +\frac i2\sum_j\partial_{\xi_j}F_j\,b\right) |d\xi\,d\eta|^{1/2}. \tag{7.2}

The derivative of FjF_j 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) S=\sum_j\partial_{\xi_j}p_j, \qquad T=\sum_{j,k}(\partial_{x_k}p_j)H_{\xi_k\xi_j}. \tag{7.3}

The chain rule gives

∑j∂ξjFj=S+T.(7.4) \sum_j\partial_{\xi_j}F_j=S+T. \tag{7.4}

On the other hand, differentiate (6.1) first in xkx_k, then in ξk\xi_k at fixed xx, and restrict to the graph. The factor xj−Hξjx_j-H_{\xi_j} then vanishes, leaving

pxkξk=∂ξkpk−∑j(∂xkpj)Hξjξk. p_{x_k\xi_k} =\partial_{\xi_k}p_k -\sum_j(\partial_{x_k}p_j)H_{\xi_j\xi_k}.

Summing and using the symmetry of the Hessian of HH yields

∑kpxkξk=S−T.(7.5) \sum_kp_{x_k\xi_k}=S-T. \tag{7.5}

Subtract (7.2) from (6.5). The remaining multiplication coefficient is

r+iS−i2(S+T)=r+i2(S−T)=r+i2∑kpxkξk=c.(7.6) r+iS-\frac i2(S+T) =r+\frac i2(S-T) =r+\frac i2\sum_kp_{x_k\xi_k}=c. \tag{7.6}

This proves (3.5), including its sign and density correction.

The calculation patches: cc 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−1t^{m-1}. One can see this directly in (7.1): its coefficients have degree mm, while a frequency derivative costs one degree. Lie differentiation therefore raises intrinsic symbol order by m−1m-1. The coefficient cc has the same degree m−1m-1. Applied to (3.3), this gives precisely m+m′−1+n/4m+m'-1+n/4, the symbol order of (3.4). For the derivative estimates, the coefficients FjF_j have degree mm in frequency coordinates, and their α\alpha-derivatives have degree m−∣α∣m-|\alpha|. The product rule applied to Fj∂ξjbF_j\partial_{\xi_j}b and (∂ξjFj)b(\partial_{\xi_j}F_j)b therefore gives order m−1+ord⁡(b)m-1+\operatorname{ord}(b), with every derivative and finitely many input seminorms. The frequency half-density frame has degree n/2n/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×St1X=\mathbb R_s\times S^1_t and Y=RzY=\mathbb R_z, and use the frames ∣ds dt∣1/2|ds\,dt|^{1/2}, ∣dz∣1/2|dz|^{1/2}. For a smooth periodic coefficient b(t)b(t), define

(Abg)(s,t)=b(t)g(s).(8.1) (A_bg)(s,t)=b(t)g(s). \tag{8.1}

Its kernel is b(t)δ(s−z)b(t)\delta(s-z), of order −1/4-1/4. The relation has s=zs=z, the same nonzero s,zs,z covector χ\chi, and zero tt covector. Coordinates on its kernel Lagrangian are (s,t,χ)(s,t,\chi). The invariant symbol is a fixed normalized Maslov factor times

b(t)∣ds dt dχ∣1/2.(8.2) b(t)|ds\,dt\,d\chi|^{1/2}. \tag{8.2}

The fixed factor is (2π)1/4(2\pi)^{1/4}, from the codimension-one delta normalization, and does not affect the following differential identity.

For real smooth periodic v(t)v(t) and a smooth scalar q(t)q(t), let

P=12(vDt+Dtv)+q=vDt−i2v′+q.(8.3) P=\frac12(vD_t+D_tv)+q =vD_t-\frac i2v'+q. \tag{8.3}

This properly supported differential operator acts on half densities. Its principal symbol is p=v(t)ξtp=v(t)\xi_t, which vanishes on the relation, and

c=−i2v′+q+i2∂t∂ξtp=q.(8.4) c=-\frac i2v'+q+\frac i2\partial_t\partial_{\xi_t}p=q. \tag{8.4}

On the kernel Lagrangian, Hp=v∂tH_p=v\partial_t; its divergence in (s,t,χ)(s,t,\chi) is v′v'. Formula (3.5) gives the coefficient

−i(vb′+12v′b)+qb.(8.5) -i\left(vb'+\frac12v'b\right)+qb. \tag{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-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>0v>0, vanishing of this coefficient is the first-order equation

b′+v′2vb+iqvb=0. b'+\frac{v'}{2v}b+\frac{iq}{v}b=0.

Its solutions are

b(t)=Kv(t)−1/2exp⁡(−i∫t0tq(u)v(u) du),(8.6) b(t)=K v(t)^{-1/2} \exp\left(-i\int_{t_0}^t\frac{q(u)}{v(u)}\,du\right), \tag{8.6}

with arbitrary complex constant KK. Differentiate to verify the equation; conversely multiplying it by v1/2exp⁡(i∫q/v)v^{1/2}\exp(i\int 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\mathbb R^2, let V=u∂u−2w∂wV=u\partial_u-2w\partial_w. Compute its flow and its Lie derivative on a half density f∣du dw∣1/2f|du\,dw|^{1/2}. Check the result by direct pullback.

Solution. The flow is Ft(u,w)=(etu,e−2tw)F_t(u,w)=(e^tu,e^{-2t}w), whose positive determinant is e−te^{-t}. Its divergence is 1−2=−11-2=-1. Formula (1.2) gives

LV(f∣du dw∣1/2)=(ufu−2wfw−12f)∣du dw∣1/2. \mathcal L_V(f|du\,dw|^{1/2}) =\left(u f_u-2w f_w-\frac12f\right)|du\,dw|^{1/2}.

Direct pullback has coefficient e−t/2f(etu,e−2tw)e^{-t/2}f(e^tu,e^{-2t}w). Its derivative at zero is the same displayed expression. Differentiating only ff would omit the determinant contribution.

Exercise 9.2 (the exceptional slopes; intermediate). In T∗R2T^*\mathbb R^2, let LL be the graph of B=(0110)B=\begin{pmatrix}0&1\\1&0\end{pmatrix}. Find the slopes (b1,b2)(b_1,b_2) for which LL fails to be transverse to the diagonal graph LbL_b. Find one complement and the dimension of the intersection at an exceptional slope.

Solution. With U=I,W=BU=I,W=B, the determinant is

det⁡(B−diag⁡(b1,b2))=b1b2−1. \det(B-\operatorname{diag}(b_1,b_2))=b_1b_2-1.

The exceptional set is the hyperbola b1b2=1b_1b_2=1, a null subset of the plane. The choice (0,0)(0,0) is transverse. At an exceptional slope, neither bjb_j 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\xi_0=2, η0=−3\eta_0=-3 on R×R\mathbb R\times\mathbb 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)(5,7). Explain the role of the two nonzero covectors.

Solution. Take

x=F(u)=u−54u2,y=G(w)=w+76w2. x=F(u)=u-\frac54u^2, \qquad y=G(w)=w+\frac76w^2.

Both maps are locally invertible at zero. Contracting their second derivatives with the fixed covectors gives 2F′′=−52F''=-5, −3G′′=−7-3G''=-7. Therefore

δξnew=δξold−5δu,δηnew=δηold−7δw. \delta\xi_{\rm new}=\delta\xi_{\rm old}-5\delta u, \qquad \delta\eta_{\rm new}=\delta\eta_{\rm old}-7\delta w.

Zero new frequency variations are precisely the old slopes (5,7)(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ℓtb(t)=e^{i\ell t}, ℓ∈Z\ell\in\mathbb Z, and P=DtP=D_t. Compute PAbPA_b, its order and its transport symbol. Include ℓ=0\ell=0.

Solution. The principal symbol ξt\xi_t vanishes on the relation, the subprincipal symbol is zero, and the tangent Hamilton field is ∂t\partial_t, with zero divergence. Hence the transport coefficient is (1/i)b′=ℓb(1/i)b'=\ell b. Direct differentiation gives Dteiℓt=ℓeiℓtD_t e^{i\ell t}=\ell e^{i\ell t}, so PAb=ℓAbPA_b=\ell A_b. For ℓ≠0\ell\neq0 its exact order is −1/4-1/4, as predicted by 1−1/4−11-1/4-1; its leading symbol is nonzero. For ℓ=0\ell=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+sin⁡tv(t)=2+\sin t, q=0q=0. Compute the left order-zero coefficient, the subprincipal symbol and the action on the lift. Find a nonzero periodic bb killed by this product.

Solution. The left order-zero coefficient is −icos⁡t/2-i\cos t/2, while the mixed principal derivative adds +icos⁡t/2+i\cos t/2. Thus the subprincipal symbol is zero. The action coefficient is

−i((2+sin⁡t)b′+12cos⁡t b). -i\left((2+\sin t)b'+\frac12\cos t\,b\right).

Take b=(2+sin⁡t)−1/2b=(2+\sin t)^{-1/2}. It is smooth, positive and periodic, and b′=−12cos⁡t(2+sin⁡t)−3/2b'=-\frac12\cos t(2+\sin t)^{-3/2}, so the two terms cancel exactly. Keeping only vb′v b', 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>0v>0 and qq are real smooth 2π2\pi-periodic functions in (8.3). Determine when PAb=0PA_b=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>0v>0. Over one period the positive factor v−1/2v^{-1/2} returns to its original value, while the phase is multiplied by

exp⁡(−i∫02πq(t)v(t) dt). \exp\left(-i\int_0^{2\pi}\frac{q(t)}{v(t)}\,dt\right).

A nonzero solution is periodic exactly when this factor is one, equivalently when the real integral lies in 2πZ2\pi\mathbb 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+sin⁡tv=2+\sin t and q=vq=v, the solution b=v−1/2e−itb=v^{-1/2}e^{-it} is periodic. For the same vv and q=v/2q=v/2, its period multiplier is −1-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

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).