Reducing Gaussian symbols and composing linear relations

Clean composition retains a density along the middle fiber. Its linear model comes from restricting a quadratic phase to the symplectic constraint and integrating the variables eliminated by the quotient. A phase radical records the directions in which that integration has no oscillation. Keeping its density makes the construction finite and invariant, even when the corresponding unweighted integral would diverge.

The exact earlier programme proofs are:

The proof map binds the results to these exact programme proofs. The primary source for this restoration is the reprint of Hörmander III, corrected second printing (1994), Section 21.6 including Theorems 21.6.6–21.6.7. The lesson retains its own explanations, density-first organization, calculations and six solved exercises. It proves the linear symbol maps; analytic composition of variable phases and amplitudes requires the additional estimates and support hypotheses proved in the analytic-composition component.

Let Ω(V)\Omega(V) denote the complex line of translation-invariant densities on a real vector space VV. Its positive real ray determines its real powers Ωs(V)\Omega^s(V), including s=±1/2s=\pm1/2. We use Ω(0)=C\Omega(0)=\mathbb C with positive unit one. The symplectic form is ω=∑dξj∧dxj\omega=\sum d\xi_j\wedge dx_j, with vertical distinguished plane λ0={x=0}\lambda_0=\{x=0\}. Write

G(λ)=M(λ)⊗Ω1/2(λ)(0.1) \mathscr G(\lambda)=M(\lambda)\otimes\Omega^{1/2}(\lambda) \tag{0.1}

for the intrinsic Gaussian symbol line constructed in the preceding lesson.

1. A quadratic phase with a radical gives a density-valued Gaussian

Let Q(x,θ)Q(x,\theta) be any real quadratic form, with x∈Rnx\in\mathbb R^n and θ∈F=RN\theta\in F=\mathbb R^N. Set

R={r∈F:Qxθr=0, Qθθr=0},e=dim⁡R.(1.1) R=\{r\in F:Q_{x\theta}r=0,\ Q_{\theta\theta}r=0\}, \qquad e=\dim R. \tag{1.1}

Thus {0}×R\{0\}\times R is the pure phase-variable part of the full Hessian radical. For every r∈Rr\in R,

Q(x,θ+r)=Q(x,θ).(1.2) Q(x,\theta+r)=Q(x,\theta). \tag{1.2}

The critical space and its image are

CQ={Qθ=0},λQ={(x,Qx):Qθ=0}.(1.3) C_Q=\{Q_\theta=0\},\qquad \lambda_Q=\{(x,Q_x):Q_\theta=0\}. \tag{1.3}

Proposition 1.1. The space λQ\lambda_Q is Lagrangian. The critical map CQ→λQC_Q\to\lambda_Q is onto with kernel {0}×R\{0\}\times R; hence CQ/R→λQC_Q/R\to\lambda_Q is an isomorphism. On any complement F=W⊕RF=W\oplus R, the restriction Q∣Rn×WQ|_{\mathbb R^n\times W} is a nondegenerate phase parametrizing this same plane.

Proof. The transpose of the linear map (x,θ)↦Qθ(x,\theta)\mapsto Q_\theta has kernel RR: its components are exactly the two matrices in (1.1). Its rank is therefore N−eN-e, so dim⁡CQ=n+e\dim C_Q=n+e. A vector in the kernel of the critical map has x=0x=0, Qxθθ=0Q_{x\theta}\theta=0 and Qθθθ=0Q_{\theta\theta}\theta=0, hence lies in RR. The image has dimension nn.

For two vectors (v,w),(v~,w~)∈CQ(v,w),(\widetilde v,\widetilde w)\in C_Q, the critical equations are

Qθxv+Qθθw=0,Qθxv~+Qθθw~=0. Q_{\theta x}v+Q_{\theta\theta}w=0, \qquad Q_{\theta x}\widetilde v+Q_{\theta\theta}\widetilde w=0.

Symmetry of the Hessian gives

ω((v,Qxxv+Qxθw),(v~,Qxxv~+Qxθw~))=0: \omega\big((v,Q_{xx}v+Q_{x\theta}w), (\widetilde v,Q_{xx}\widetilde v+Q_{x\theta}\widetilde w)\big)=0:

the QxxQ_{xx} terms cancel, and the remaining difference is zero by the two critical equations and symmetry of QθθQ_{\theta\theta}. Thus the nn-dimensional image is isotropic and is Lagrangian.

On a complement WW, a pure phase radical vector would also be a vector of RR, so is zero. This proves nondegeneracy of the restricted phase. Equation (1.2) shows that it has the same critical image. ∎

The integral of eiQe^{iQ} over all of FF has an infinite volume factor when e>0e>0. Its appropriate replacement is an element of

I(λQ,H)⊗Ω(R),(1.4) \mathcal I(\lambda_Q,H)\otimes\Omega(R), \tag{1.4}

where HH is the horizontal reference. Fix Lebesgue measure dθd\theta, a scalar aa, and any linear projection T:F→RT:F\to R that is the identity on RR. For χ∈Cc∞(R)\chi\in C_c^\infty(R), consider

UQ[χ](x)=a(2π)−(n+2N)/4∫FeiQ(x,θ)χ(Tθ) dθ ∣dx∣1/2.(1.5) \mathcal U_Q[\chi](x) =a(2\pi)^{-(n+2N)/4} \int_F e^{iQ(x,\theta)}\chi(T\theta)\,d\theta\ |dx|^{1/2}. \tag{1.5}

The integration in the quotient directions is the Fresnel/delta evaluation of the preceding lesson; the radical direction has compact support.

Lemma 1.2. Formula (1.5) is independent of TT and is the pairing of a well-defined translation-invariant I(λQ,H)\mathcal I(\lambda_Q,H)-valued density on RR with χ\chi.

Proof. Use a splitting θ=w+r\theta=w+r. The density exact sequence identifies dθd\theta with a product μWμR\mu_W\mu_R; any reciprocal rescaling of those two factors leaves their product unchanged. By (1.2), the phase is independent of rr. Since T(w+r)=Tw+rT(w+r)=Tw+r, translation invariance gives

∫Rχ(Tw+r) μR(r)=∫Rχ(r) μR(r). \int_R\chi(Tw+r)\,\mu_R(r)=\int_R\chi(r)\,\mu_R(r).

Thus (1.5) is

[a(2π)−(n+2N)/4∫WeiQ(x,w) μW(w) ∣dx∣1/2]∫Rχ(r) μR(r).(1.6) \left[a(2\pi)^{-(n+2N)/4} \int_W e^{iQ(x,w)}\,\mu_W(w)\ |dx|^{1/2}\right] \int_R\chi(r)\,\mu_R(r). \tag{1.6}

The first bracket belongs to the nonzero Gaussian line unless a=0a=0, by Proposition 1.1 and the previous lesson. Changing the complement replaces each ww by w+Lww+Lw, with Lw∈RLw\in R, which leaves the phase unchanged. The density exact sequence accounts for its measure change. Therefore the bracket tensored with μR\mu_R is independent of the splitting as well as of TT. This proves the claim. One may justify the displayed integrations by inserting a damping factor in WW, performing the compact radical integration first, then taking the proved Fresnel boundary value. ∎

With the same quotient measure, the normalized nondegenerate Gaussian for Q∣WQ|_W would use N−eN-e phase variables. Consequently (1.6) is that Gaussian tensored with μR\mu_R, multiplied by

(2π)−e/2.(1.7) (2\pi)^{-e/2}. \tag{1.7}

Here we retain the normalization with the original count NN. Later it is the actual integral over eliminated variables that determines that count. An adjusted clean-phase normalization can absorb (1.7), but that adjustment must be explicit.

2. The density map supplied by symplectic reduction

Let SS have dimension 2n2n, let Δ⊂S\Delta\subset S be isotropic, and set

d=dim⁡Δ,K=λ∩Δ,e=dim⁡K. d=\dim\Delta,\qquad K=\lambda\cap\Delta,\qquad e=\dim K.

Define

SΔ=Δω/Δ,λΔ=(λ∩Δω)/K,λ0Δ=(λ0∩Δω)/(λ0∩Δ).(2.1) S_\Delta=\Delta^\omega/\Delta, \qquad \lambda_\Delta=(\lambda\cap\Delta^\omega)/K, \qquad \lambda_{0\Delta}=(\lambda_0\cap\Delta^\omega)/(\lambda_0\cap\Delta). \tag{2.1}

The images implicit in these quotient formulas are Lagrangian in SΔS_\Delta, by the proved linear reduction proposition. Their dimensions are n−dn-d.

There is a canonical positive isomorphism

dΔ:Ω1/2(λ)⟶Ω1/2(λΔ)⊗Ω(K)⊗Ω−1/2(Δ).(2.2) \mathfrak d_\Delta: \Omega^{1/2}(\lambda) \longrightarrow \Omega^{1/2}(\lambda_\Delta) \otimes\Omega(K)\otimes\Omega^{-1/2}(\Delta). \tag{2.2}

Proof. Put B=λ∩ΔωB=\lambda\cap\Delta^\omega. The two exact sequences give

Ω1/2(λ)≃Ω1/2(B)⊗Ω1/2(λ/B),Ω1/2(B)≃Ω1/2(K)⊗Ω1/2(λΔ).(2.3) \Omega^{1/2}(\lambda) \simeq\Omega^{1/2}(B)\otimes\Omega^{1/2}(\lambda/B), \qquad \Omega^{1/2}(B) \simeq\Omega^{1/2}(K)\otimes\Omega^{1/2}(\lambda_\Delta). \tag{2.3}

For completeness, the density isomorphism of an exact sequence is obtained by lifting a basis of the quotient and adjoining a basis of the kernel. Changing a lift adds kernel columns and has determinant one. The absolute determinant of a block change of basis is the product of the two diagonal absolute determinants. This proves canonicity and positivity.

The pairing

(λ/B)×(Δ/K)⟶R,([v],[w])⟼ω(v,w)(2.4) (\lambda/B)\times(\Delta/K)\longrightarrow\mathbb R, \qquad([v],[w])\longmapsto\omega(v,w) \tag{2.4}

is well-defined and perfect. Its left kernel is BB; its right kernel before quotienting is Δ∩λω=K\Delta\cap\lambda^\omega=K. Both quotients have dimension d−ed-e. Duality therefore identifies

Ω1/2(λ/B)≃Ω−1/2(Δ/K)≃Ω1/2(K)⊗Ω−1/2(Δ).(2.5) \Omega^{1/2}(\lambda/B) \simeq\Omega^{-1/2}(\Delta/K) \simeq\Omega^{1/2}(K)\otimes\Omega^{-1/2}(\Delta). \tag{2.5}

Combine (2.3) and (2.5). The two factors Ω1/2(K)\Omega^{1/2}(K) become the full density Ω(K)\Omega(K), proving (2.2). ∎

This density map uses only exact sequences and the symplectic pairing. It has no Fresnel phase and no factor of 2π2\pi. The Gaussian construction below has the same positive Jacobians, with an additional normalization factor.

3. Adapted coordinates and the quadratic reduction operation

Set

k=dim⁡(Δ∩λ0),j=d−k,l=n−d.(3.1) k=\dim(\Delta\cap\lambda_0),\quad j=d-k,\quad l=n-d. \tag{3.1}

Choose symplectic coordinates split as

x=(x′,y,z),ξ=(ξ′,η,ζ),dim⁡x′=k,dim⁡y=l,dim⁡z=j,(3.2) x=(x',y,z),\qquad\xi=(\xi',\eta,\zeta), \quad\dim x'=k,\quad\dim y=l,\quad\dim z=j, \tag{3.2}

such that

λ0={x=0},Δ={x′=y=η=ζ=0}.(3.3) \lambda_0=\{x=0\},\qquad \Delta=\{x'=y=\eta=\zeta=0\}. \tag{3.3}

Thus (ξ′,z)(\xi',z) are coordinates on Δ\Delta,

Δω={x′=ζ=0},(3.4) \Delta^\omega=\{x'=\zeta=0\}, \tag{3.4}

and (y,η)(y,\eta) are symplectic coordinates on the quotient.

Such coordinates exist with the specified distinguished vertical plane. Choose a basis of Δ∩λ0\Delta\cap\lambda_0, extend it to Δω∩λ0\Delta^\omega\cap\lambda_0, then to λ0\lambda_0. The dimensions of these successive spaces are k,k+l,nk,k+l,n: restriction of the symplectic pairing of Δ\Delta with λ0\lambda_0 has rank jj. Choose jj vectors of Δ\Delta dual to the last vertical basis vectors; they complete Δ∩λ0\Delta\cap\lambda_0 to Δ\Delta and are mutually isotropic. Here is a completion that keeps those prescribed vectors fixed. Write the vertical basis as eie_i, and call the prescribed horizontal vectors fjf_j, with indices in the last block JJ, so ω(ei,fj)=δij\omega(e_i,f_j)=\delta_{ij}. For each remaining index choose a vector gig_i dual to all eae_a; this is possible since the pairing identifies S/λ0S/\lambda_0 with λ0∗\lambda_0^*. Replace gig_i by gi−∑j∈Jω(gi,fj)ejg_i-\sum_{j\in J}\omega(g_i,f_j)e_j. It stays dual to the vertical basis and is now orthogonal to every prescribed fjf_j. Finally replace each such gig_i by gi+12∑a∉Jω(ga,gi)eag_i+\tfrac12\sum_{a\notin J}\omega(g_a,g_i)e_a, using the vectors before this last replacement in the coefficients. The new pairing between indices i,b∉Ji,b\notin J is ω(gi,gb)−12ω(gi,gb)+12ω(gb,gi)=0\omega(g_i,g_b)-\tfrac12\omega(g_i,g_b)+\tfrac12\omega(g_b,g_i)=0. This last correction preserves the pairings with prescribed vectors. Together with the eie_i and fixed fjf_j, the resulting vectors form the required symplectic basis. This gives (3.2)–(3.3), including empty blocks.

Represent an element of G(λ)\mathscr G(\lambda) in the horizontal reference by a nondegenerate quadratic phase:

u(x)=a(2π)−(n+2N)/4∫eiQ(x,θ) dθ ∣dx∣1/2.(3.5) u(x)=a(2\pi)^{-(n+2N)/4} \int e^{iQ(x,\theta)}\,d\theta\ |dx|^{1/2}. \tag{3.5}

Such a phase exists: complete the nonzero vertical constraints of λ\lambda, and write it as q(x)+x⋅Tθq(x)+x\cdot T\theta with TT injective and qq quadratic, as in the previous lesson's signature proof.

Form the new quadratic phase

qΔ(y,z,θ)=Q(0,y,z,θ),(3.6) q_\Delta(y,z,\theta)=Q(0,y,z,\theta), \tag{3.6}

with base variable yy and phase variables (z,θ)(z,\theta). Interpret

g(y)=a(2π)−(n+2N)/4∬eiqΔ(y,z,θ) dz dθ ∣dy∣1/2(3.7) g(y)=a(2\pi)^{-(n+2N)/4} \iint e^{iq_\Delta(y,z,\theta)}\,dz\,d\theta\ |dy|^{1/2} \tag{3.7}

by Lemma 1.2, as a density-valued Gaussian. Equation (3.7) is a quadratic phase operation. It does not assert that an arbitrary distribution admits ordinary restriction to x′=0x'=0 followed by an unweighted pushforward.

Lemma 3.1. The pure phase radical of (3.6) is naturally isomorphic to KK, and its critical image is λΔ\lambda_\Delta. Hence (3.7) belongs to

I(λΔ,HΔ)⊗Ω(K).(3.8) \mathcal I(\lambda_\Delta,H_\Delta)\otimes\Omega(K). \tag{3.8}

Proof. Its critical equations are Qz=Qθ=0Q_z=Q_\theta=0, evaluated at x′=0x'=0. In the original critical map, these are precisely x′=ζ=0x'=\zeta=0, so the original covector lies in λ∩Δω\lambda\cap\Delta^\omega. Projection to (y,η)(y,\eta) gives the reduced plane in (2.1).

For a pure radical vector (z,θ)(z,\theta), set x′=y=0x'=y=0. All derivatives of (3.6), including its yy derivative, vanish. Thus its image in the original critical map has x′=y=η=ζ=0x'=y=\eta=\zeta=0, and lies in λ∩Δ=K\lambda\cap\Delta=K. Conversely every vector of KK has a unique preimage under the original nondegenerate critical map; that preimage is a pure radical vector for (3.6). This is the claimed isomorphism. Lemma 1.2 and Proposition 1.1 now prove (3.8). ∎

Tensor (3.8) with the inverse half density

∣dξ′ dz∣−1/2on Δ,(3.9) |d\xi'\,dz|^{-1/2}\quad\text{on }\Delta, \tag{3.9}

and pass to the intrinsic Gaussian symbol lines. This is the candidate map

RΔ:G(λ)⟶G(λΔ)⊗Ω(K)⊗Ω−1/2(Δ).(3.10) \mathfrak R_\Delta: \mathscr G(\lambda)\longrightarrow \mathscr G(\lambda_\Delta)\otimes\Omega(K) \otimes\Omega^{-1/2}(\Delta). \tag{3.10}

4. Why the reduction operation is independent of its choices

Theorem 4.1 (Gaussian reduction). The map (3.10) is a canonical nonzero linear map. It depends on (S,λ0,Δ)(S,\lambda_0,\Delta) and λ\lambda, and is invariant under symplectic isomorphisms preserving those data. For fixed (S,λ0,Δ)(S,\lambda_0,\Delta), on any smooth family of λ\lambda's with dim⁡(λ∩Δ)\dim(\lambda\cap\Delta) constant, it is a smooth bundle map. Its positive density factor is

(2π)(d−2k−2e)/4dΔ,(4.1) (2\pi)^{(d-2k-2e)/4}\mathfrak d_\Delta, \tag{4.1}

and its unit-modulus factor gives a canonical map

M(λ)⟶MΔ(λΔ).(4.2) M(\lambda)\longrightarrow M_\Delta(\lambda_\Delta). \tag{4.2}

Proof of independence. First change the phase representing the same input Gaussian. Complete squares in the nonzero block of QθθQ_{\theta\theta}. Eliminating hh variables by the exact Fresnel formula multiplies the amplitude by

eiπsgn⁡D/4∣det⁡D∣−1/2 e^{i\pi\operatorname{sgn}D/4}|\det D|^{-1/2}

and changes the normalization from NN to N−hN-h, because the Gaussian factor (2π)h/2(2\pi)^{h/2} supplies exactly that change. The same elimination, with the same factor, can be performed in (3.7); setting x′=0x'=0 does not change the eliminated pure phase block. Thus both operations reduce to a phase linear in its remaining phase variables,

Q(x,θ)=q(x)+x⋅Tθ,(4.3) Q(x,\theta)=q(x)+x\cdot T\theta, \tag{4.3}

where TT is injective. Its constraint subspace V=ker⁡TTV=\ker T^T is the base projection of λ\lambda, and q∣Vq|_V is determined by λ\lambda, since q(x)=ξ⋅x/2q(x)=\xi\cdot x/2 there. Two such injective TT's with the same range differ by an invertible phase-basis change. Two such qq's differ by a quadratic polynomial vanishing on VV. In linear coordinates V={x1=⋯=xs=0}V=\{x_1=\cdots=x_s=0\}, every monomial of that polynomial has a constrained coordinate as a factor. Consequently the difference is x⋅TLxx\cdot T Lx for a linear map LL, and can be absorbed by the translation θ↦θ+Lx\theta\mapsto\theta+Lx. These changes preserve the integral and its measure, both before and after (3.6). The scalar amplitudes agree after the measure changes because the resulting nonzero input Gaussian is fixed. This proves phase independence, including unequal original variable counts.

Next keep the base coordinates and change the horizontal reference by ξ↦ξ−Ax\xi\mapsto\xi-Ax, A=ATA=A^T, preserving (3.3). This preservation requires A(0,0,z)A(0,0,z) to have only a ξ′\xi' component. Hence the yzyz and zzzz blocks of AA vanish. On x′=0x'=0, its quadratic form is independent of zz and is yTAyyyy^TA_{yy}y. The input chirp e−ixTAx/2e^{-ix^TAx/2} in the preceding lesson therefore changes (3.7) by precisely e−iyTAyyy/2e^{-iy^TA_{yy}y/2}, the induced chirp on the reduced space. Its shear on Δ\Delta has determinant one. Thus (3.10) respects the intrinsic Gaussian-line identifications.

Finally keep the horizontal reference, and change dual bases by x=Pvx=Pv, ξ=P−Tν\xi=P^{-T}\nu, preserving the two flags x′=0x'=0 and x′=y=0x'=y=0. The matrix has the block form

P=(P1100P21P220P31P32P33).(4.4) P=\begin{pmatrix}P_{11}&0&0\\P_{21}&P_{22}&0\\P_{31}&P_{32}&P_{33}\end{pmatrix}. \tag{4.4}

At v′=0v'=0, the remaining base is y=P22vyy=P_{22}v_y, and the integrated variable is z=P32vy+P33vzz=P_{32}v_y+P_{33}v_z. Pullback of the input half density and integration over vzv_z give the factor

∣det⁡P∣1/2∣det⁡P33∣−1=∣det⁡P22∣1/2∣det⁡P11∣1/2∣det⁡P33∣−1/2.(4.5) |\det P|^{1/2}|\det P_{33}|^{-1} =|\det P_{22}|^{1/2} |\det P_{11}|^{1/2}|\det P_{33}|^{-1/2}. \tag{4.5}

The first factor is the reduced half-density Jacobian. The other factors are exactly the transformation of Ω−1/2(Δ)\Omega^{-1/2}(\Delta): on Δ\Delta, the old coordinates (ξ′,z)(\xi',z) are (P11−Tν′,P33vz)(P_{11}^{-T}\nu',P_{33}v_z). The radical density is transported to KK by the critical-map isomorphism in Lemma 3.1. Equations (4.5) and (3.9) therefore cancel all remaining basis dependence. A general adapted coordinate change is the composition of this dual-basis change and the preceding shear: a symplectic map preserving the vertical plane has precisely a dual linear base change followed by a symmetric frequency shear. This proves intrinsic independence.

For smoothness, work near one parameter value. Choose a fixed horizontal Lagrangian transverse to both the distinguished vertical and the input plane there. Transversality persists nearby, so in the associated coordinates the input planes have the smooth graph form x=B(p)θx=B(p)\theta, with B(p)=B(p)TB(p)=B(p)^T. Returning to the fixed adapted coordinates only adds a fixed quadratic base form and an invertible base change. Thus they admit a smooth family of nondegenerate quadratic phases with the same number of variables, N=nN=n. The matrix of the reduced critical equations has constant rank because its transpose kernel has dimension ee, by Lemma 3.1. An invertible matrix minor expresses its kernel as a smooth graph: the remaining rows impose no additional condition since that minor already has the full rank. Applying the same argument to the transpose gives a smooth frame and complement for the radical. Quotienting by it gives a smooth nondegenerate phase qˉp\bar q_p. Its critical map is injective, and a fixed nonzero minor selects a smooth basis for its image λΔ(p)\lambda_\Delta(p).

Choose in the reduced space a fixed horizontal plane transverse to that image and to its distinguished vertical at the central parameter. It remains transverse nearby. In this horizontal reference the full Hessian of qˉp\bar q_p, including its base variables, is invertible: a vector in its kernel lies in the phase critical space and maps to a vector in the intersection of the reduced plane with that horizontal; transversality makes the image zero and nondegeneracy then makes the vector zero. The determinant is therefore nonzero nearby, its signature is locally constant, and its inverse and critical-density Jacobians are smooth. The exact Gaussian formula in the preceding lesson now gives smooth nonzero symbol coefficients in this frame. The density identification with KK is a smooth linear isomorphism by Lemma 3.1. The transition laws already proved transport these coefficients to any other frame.

This argument does not require the pure phase Hessian to have constant rank. The distributional smoothness of the corresponding Gaussian frames, even when that rank changes, is proved by Fourier-test differentiation in the linked Gaussian companion Z2. It follows that (3.10) is a smooth bundle map on the stated constant-ee family. At a jump of ee, the radical spaces no longer form a vector bundle of one rank; the theorem makes no smoothness assertion across that jump.

The map is nonzero because removing its radical leaves a nondegenerate quadratic phase, whose Gaussian is nonzero. Its phase and positive factor separate uniquely: the positive-density rays specify the positive factor, and the Maslov transitions have modulus one. It remains to compute that positive factor; Sections 5–6 do so and prove (4.1). ∎

5. Compute the factor when the constraint lies in the plane

Assume first that Δ⊂λ\Delta\subset\lambda, so e=d=k+je=d=k+j. In adapted coordinates, choose reduced horizontal coordinates making λΔ\lambda_\Delta the graph y=Bηy=B\eta, B=BTB=B^T. This is possible by a common complement to λ0Δ\lambda_{0\Delta} and λΔ\lambda_\Delta, lifted to a horizontal complement preserving (3.3).

Then λ\lambda is represented by

Q(x′,y,z;θ′,η)=x′⋅θ′+y⋅η−12ηTBη.(5.1) Q(x',y,z;\theta',\eta) =x'\cdot\theta'+y\cdot\eta-\tfrac12\eta^TB\eta. \tag{5.1}

Indeed inclusion of Δ\Delta forces x′=ζ=0x'=\zeta=0 on λ\lambda; reduction gives y=Bηy=B\eta, and the remaining ξ′\xi' and zz variables are free. Here the phase-variable count is N=k+lN=k+l. Its critical density is

∣dθ′ dη dz∣, |d\theta'\,d\eta\,dz|,

so the positive input symbol is ∣a∣∣dθ′ dη dz∣1/2|a||d\theta'\,d\eta\,dz|^{1/2}.

After setting x′=0x'=0, the variables (θ′,z)(\theta',z) are the radical and identify with K=ΔK=\Delta. The remaining reduced Gaussian has coefficient aa times

(2π)−(n+2k+2l)/4+3l/4=(2π)−(3k+j)/4.(5.2) (2\pi)^{-(n+2k+2l)/4+3l/4} =(2\pi)^{-(3k+j)/4}. \tag{5.2}

Thus, apart from its unit phase, the output of (3.10) is

∣a∣(2π)−(3k+j)/4∣dη∣1/2⊗∣dθ′ dz∣⊗∣dξ′ dz∣−1/2.(5.3) |a|(2\pi)^{-(3k+j)/4}|d\eta|^{1/2} \otimes|d\theta'\,dz|\otimes|d\xi'\,dz|^{-1/2}. \tag{5.3}

The exact-sequence map (2.2) has exactly these positive density frames without the scalar factor. Since

−3k−j=d−2k−2e, -3k-j=d-2k-2e,

this proves (4.1) in the contained case.

6. Compute the factor when the intersection is zero

Now assume λ∩Δ=0\lambda\cap\Delta=0. We need a horizontal complement HH transverse to both λ0\lambda_0 and λ\lambda, preserving the adapted form of Δ\Delta, whose induced reduced horizontal plane is transverse to λΔ\lambda_\Delta. Here is a construction that keeps all those requirements.

Choose V⊂ΔV\subset\Delta complementary to Δ∩λ0\Delta\cap\lambda_0. It is transverse to both λ0\lambda_0 and λ\lambda. Reduce first by VV. In that quotient, the remaining constraint D=Δ/VD=\Delta/V lies in the distinguished vertical plane and is transverse to the reduced λ\lambda. In the further quotient by DD, choose a common horizontal complement to the two reduced Lagrangians. Lift it into a symplectic decomposition of the quotient by VV with coordinates (x′,y;ξ′,η)(x',y;\xi',\eta), where DD is the ξ′\xi' plane and the selected horizontal plane in the further quotient is η=0\eta=0.

For the reduced λ\lambda, its subspace η=0\eta=0 projects bijectively to the x′x' coordinates. To check this, a vector there with x′=0x'=0 projects to the intersection with the selected reduced horizontal plane, so its projection is zero; it then belongs to λ∩D=0\lambda\cap D=0. Also the η\eta map on this Lagrangian has rank ll: its transpose kernel is its intersection with the yy-horizontal plane, which has just been seen to be zero. Therefore its kernel has dimension kk, equal to the number of x′x' coordinates. On this kernel, ξ′=Cx′\xi'=Cx' for a symmetric CC, by isotropy. The plane ξ′=tx′\xi'=t x', η=0\eta=0 is transverse to it if C−tIC-tI is invertible; choose real tt outside the finite eigenvalue set. This plane is a Lagrangian complement to the vertical plane, is transverse to the reduced λ\lambda, and induces the previously selected horizontal plane after reduction by DD. Its preimage under reduction by VV contains VV and is a horizontal complement with all the stated properties in SS.

With that choice, write

λ={x=Bξ},B=BT,ξ=(ξ′,η,ζ).(6.1) \lambda=\{x=B\xi\},\qquad B=B^T, \quad \xi=(\xi',\eta,\zeta). \tag{6.1}

The block B11B_{11} is invertible. Indeed, the constrained equations are ζ=0\zeta=0, (Bξ)′=0(B\xi)'=0; a solution with η=0\eta=0 gives a vector of λΔ\lambda_\Delta in its horizontal complement, hence is zero there and then lies in K=0K=0. Thus B11ξ′=0B_{11}\xi'=0 implies ξ′=0\xi'=0.

The normalized input is

u(x)=c(2π)−3n/4∫ei(x⋅ξ−ξTBξ/2) dξ ∣dx∣1/2.(6.2) u(x)=c(2\pi)^{-3n/4} \int e^{i(x\cdot\xi-\xi^TB\xi/2)}\,d\xi\ |dx|^{1/2}. \tag{6.2}

In (3.7), integration over zz gives (2π)jδ0(ζ)(2\pi)^j\delta_0(\zeta). Completing squares in ξ′\xi' then gives the reduced graph matrix

BΔ=B22−B21B11−1B12(6.3) B_\Delta=B_{22}-B_{21}B_{11}^{-1}B_{12} \tag{6.3}

and the reduced Gaussian coefficient

cΔ=c(2π)(j−k)/4e−iπsgn⁡B11/4∣det⁡B11∣−1/2.(6.4) c_\Delta=c(2\pi)^{(j-k)/4} e^{-i\pi\operatorname{sgn}B_{11}/4}|\det B_{11}|^{-1/2}. \tag{6.4}

All constants follow directly: the zz integral contributes (2π)j(2\pi)^j, the kk-variable Fresnel integral contributes (2π)k/2(2\pi)^{k/2}, and converting from the input coefficient (2π)−3n/4(2\pi)^{-3n/4} to the reduced coefficient (2π)−3l/4(2\pi)^{-3l/4} gives (2π)(j−k)/4(2\pi)^{(j-k)/4}. The negative signature is due to the Hessian −B11-B_{11}.

To compare with (2.2), change coordinates on λ\lambda from (ξ′,η,ζ)(\xi',\eta,\zeta) to (x′,η,ζ)(x',\eta,\zeta). The determinant is ∣det⁡B11∣|\det B_{11}|, so

c∣dξ′ dη dζ∣1/2=c∣det⁡B11∣−1/2∣dx′ dη dζ∣1/2.(6.5) c|d\xi'\,d\eta\,d\zeta|^{1/2} =c|\det B_{11}|^{-1/2}|dx'\,d\eta\,d\zeta|^{1/2}. \tag{6.5}

The coordinates on λ/B\lambda/B are (x′,ζ)(x',\zeta). In the symplectic pairing (2.4) they are dual to (ξ′,z)(\xi',z) on Δ\Delta, up to a sign in the first block, irrelevant to a positive density. The geometric map therefore sends (6.5) to

c∣det⁡B11∣−1/2∣dη∣1/2⊗∣dξ′ dz∣−1/2.(6.6) c|\det B_{11}|^{-1/2}|d\eta|^{1/2} \otimes|d\xi'\,dz|^{-1/2}. \tag{6.6}

Comparison with (6.4) proves (4.1), since j−k=d−2kj-k=d-2k and e=0e=0. It also gives the explicit Maslov map in these graph frames: multiplication by e−iπsgn⁡B11/4e^{-i\pi\operatorname{sgn}B_{11}/4}. A map between symbol lines of different dimensions can have an eighth-root coefficient in these frames; the geometric Maslov bundles themselves retain their fourth-root transition groups.

For a general intersection KK, reduce first by K⊂λK\subset\lambda, then by Δ/K\Delta/K. The second reduced intersection is zero. If k1=dim⁡(K∩λ0)k_1=\dim(K\cap\lambda_0), the second vertical intersection has dimension k−k1k-k_1: a vector of Δ/K\Delta/K in the reduced vertical plane lifts to a vector of (Δ∩λ0)+K(\Delta\cap\lambda_0)+K. The two exponents already computed add to

e−2k1−2e4+(d−e)−2(k−k1)4=d−2k−2e4.(6.7) \frac{e-2k_1-2e}{4} +\frac{(d-e)-2(k-k_1)}4 =\frac{d-2k-2e}{4}. \tag{6.7}

The agreement with direct reduction can be checked without interchanging divergent integrals. We prove the slightly more general assertion for any isotropic flag Δ1⊂Δ\Delta_1\subset\Delta. Choose a vertical basis starting with Δ1∩λ0\Delta_1\cap\lambda_0, extend it through Δ∩λ0\Delta\cap\lambda_0 and Δω∩λ0\Delta^\omega\cap\lambda_0, and then to λ0\lambda_0. In Δ/(Δ∩λ0)\Delta/(\Delta\cap\lambda_0), start with the image of Δ1\Delta_1 and extend a basis; choose representatives of the first vectors in Δ1\Delta_1. Choose the remaining vertical basis dual to these representatives. The completion in Section 3 fixes these isotropic representatives. Thus the adapted coordinates can be split as x′=(a,b)x'=(a,b) and z=(s,t)z=(s,t), so that Δ1\Delta_1 uses precisely the ξa\xi_a and ss coordinates, and Δ\Delta uses all of (ξ′,z)(\xi',z).

The first reduction sets a=0a=0 and treats (s,θ)(s,\theta) as phase variables, with remaining base (b,y,t)(b,y,t). Its radical R1R_1 identifies with λ∩Δ1\lambda\cap\Delta_1. Every vector of R1R_1 annihilates the full Hessian after a=0a=0, so it also annihilates the final restricted phase after b=0b=0, where (t,s,θ)(t,s,\theta) are phase variables. Hence R1⊂RR_1\subset R, the final radical. Choose a complement to R1R_1 in the first phase-variable space and retain its density by Lemma 1.2. Before any Fresnel evaluation, keep the first quotient phase itself as the nondegenerate representation of the first output. Applying the second reduction sets b=0b=0 and adds tt to the phase variables. Its radical is exactly R/R1R/R_1: the restricted phase is already independent of R1R_1, and a vector annihilates the descended full Hessian precisely when any lift annihilates the original restricted full Hessian. The final quotient is consequently the same quadratic form on the same quotient space as direct reduction.

The measures agree by the determinant-one lift rule for density exact sequences. In particular Ω(R1)⊗Ω(R/R1)≃Ω(R)\Omega(R_1)\otimes\Omega(R/R_1)\simeq\Omega(R), and these identifications coincide under the critical maps with the corresponding intersections in λ\lambda. The original prefactor is unchanged while a radical density is split off. If one instead rewrites the intermediate Gaussian with the normalized prefactor for its new base and phase-variable counts, its scalar coefficient acquires the ratio of the old prefactor to the new one; at the next step these ratios telescope to the direct ratio. The inverse half densities of the two constraint spaces likewise combine to Ω−1/2(Δ)\Omega^{-1/2}(\Delta). Thus direct and staged operations are equal as density-valued Gaussian symbols. This proves the phase assertion as well as the absolute factor; neither equality assumes an unweighted integral over a radical.

The same basis-lift construction combines the geometric maps (2.2): the successive exact sequences are refinements of the direct sequences, and their dual quotient bases pair by the same symplectic form. Taking Δ1=K\Delta_1=K, the first intersection is all of KK and the second is zero. The scalar computation (6.7) therefore completes Theorem 4.1 for every KK.

7. Apply reduction to composition of linear canonical relations

Let SiS_i be symplectic of dimension 2ni2n_i, with distinguished Lagrangian λi0\lambda_{i0}, i=1,2,3i=1,2,3. Let

G1⊂S1×S2,G2⊂S2×S3 G_1\subset S_1\times S_2,\qquad G_2\subset S_2\times S_3

be Lagrangian for ω1−ω2\omega_1-\omega_2 and ω2−ω3\omega_2-\omega_3, respectively. Their composition and excess space are

G=G1∘G2={(s1,s3):(s1,s2)∈G1,(s2,s3)∈G2 for some s2}, G=G_1\circ G_2 =\{(s_1,s_3): (s_1,s_2)\in G_1, (s_2,s_3)\in G_2\text{ for some }s_2\},
N={s2:(0,s2)∈G1, (s2,0)∈G2},e=dim⁡N.(7.1) N=\{s_2:(0,s_2)\in G_1,\ (s_2,0)\in G_2\}, \qquad e=\dim N. \tag{7.1}

All intersections here are linear and hence clean. The linear reduction proof shows that GG is Lagrangian, and that the matching space

F={(s1,s2,s2,s3)∈G1×G2} F=\{(s_1,s_2,s_2,s_3)\in G_1\times G_2\}

fits into 0→N→F→G→00\to N\to F\to G\to0.

Theorem 7.1 (linear symbol composition). There are canonical bilinear maps

M(G1)×M(G2)⟶M(G),(7.2) M(G_1)\times M(G_2)\longrightarrow M(G), \tag{7.2}
Ω1/2(G1)×Ω1/2(G2)⟶Ω1/2(G)⊗Ω(N).(7.3) \Omega^{1/2}(G_1)\times\Omega^{1/2}(G_2) \longrightarrow\Omega^{1/2}(G)\otimes\Omega(N). \tag{7.3}

The second map is (2π)−e/2(2\pi)^{-e/2} times the positive exact-sequence/duality map. These maps are smooth on constant-excess families and are compatible with the tensor product Gaussian model.

Proof. In

S=S1×S2×S2×S3,ω=ω1−ω2+ω2−ω3, S=S_1\times S_2\times S_2\times S_3, \qquad \omega=\omega_1-\omega_2+\omega_2-\omega_3,

take λ=G1×G2\lambda=G_1\times G_2 and the middle diagonal constraint

Δ={(0,s2,s2,0):s2∈S2}.(7.4) \Delta=\{(0,s_2,s_2,0):s_2\in S_2\}. \tag{7.4}

It is isotropic. Its orthogonal imposes matching middle vectors, so λ∩Δω=F\lambda\cap\Delta^\omega=F; the quotient is S1×S3S_1\times S_3, the reduced distinguished plane is λ10×λ30\lambda_{10}\times\lambda_{30}, and K=λ∩ΔK=\lambda\cap\Delta identifies with NN. Also

d=2n2,k=n2.(7.5) d=2n_2,\qquad k=n_2. \tag{7.5}

The density factor (4.1) is consequently (2π)−e/2(2\pi)^{-e/2}.

Trivialize the residual factor Ω−1/2(Δ)\Omega^{-1/2}(\Delta) by the positive symplectic volume

ρ2=∣ω2n2/n2!∣(7.6) \rho_2=|\omega_2^{n_2}/n_2!| \tag{7.6}

on S2≃ΔS_2\simeq\Delta. This contraction is intrinsic; no choice of a Euclidean middle measure enters. Density lines on the product identify with tensor products, and the Gaussian line of a direct sum is the tensor product of its Gaussian lines: add the two quadratic phases, multiply their amplitudes, and observe that their normalized powers of 2π2\pi add. Apply Theorem 4.1 to this product. Its unit phase gives (7.2), and its positive factor gives (7.3). Smoothness follows from the constant-intersection assertion of Theorem 4.1. ∎

The output is a half density on GG with a full density on the excess space NN. This is the linear source of the fiber density integrated in analytic clean composition. A translation-invariant density on N≠0N\ne0 does not have a finite total integral; the subsequent analytic theorem must supply properness or compact fiber support and amplitudes.

8. A graph factor removes the excess density

Suppose G1G_1 is the graph of a symplectic isomorphism κ:S2→S1\kappa:S_2\to S_1, so n1=n2n_1=n_2. Then N=0N=0, and

G2⟶G,(s2,s3)⟼(κs2,s3)(8.1) G_2\longrightarrow G,\qquad(s_2,s_3)\longmapsto(\kappa s_2,s_3) \tag{8.1}

is an isomorphism. Lift ρ2\rho_2 to a density ρG1\rho_{G_1} on the graph. If d1=cρG11/2d_1=c\rho_{G_1}^{1/2} and d2∈Ω1/2(G2)d_2\in\Omega^{1/2}(G_2), the positive map (7.3) is

(d1,d2)⟼c (8.1)∗d2.(8.2) (d_1,d_2)\longmapsto c\,(8.1)_*d_2. \tag{8.2}

Proof. A vector of G1×G2G_1\times G_2 has the unique decomposition

(κs2,s2,s2′,s3)=(κ(s2−s2′),s2−s2′,0,0)+(κs2′,s2′,s2′,s3).(8.3) (\kappa s_2,s_2,s_2',s_3) =(\kappa(s_2-s_2'),s_2-s_2',0,0) +(\kappa s_2',s_2',s_2',s_3). \tag{8.3}

The second term is in FF; the first identifies the quotient by FF with the graph G1G_1. Pairing that quotient with the middle diagonal is the symplectic self-duality of S2S_2, up to a sign, which does not affect positive densities. The contraction in (7.6) therefore contracts d1d_1 with ρG1−1/2\rho_{G_1}^{-1/2}, leaving precisely cc. There is no excess factor because e=0e=0. This proves (8.2). The Maslov map (7.2) still transports its line coefficient; (8.2) describes the positive density factor. ∎

For the identity relation, its natural Gaussian phase (x−y)⋅θ(x-y)\cdot\theta has normalization (2π)−n2(2\pi)^{-n_2}. It represents δ(x−y)∣dx dy∣1/2\delta(x-y)|dx\,dy|^{1/2}, with unit critical amplitude. Composing two such identity models gives the same model, so the identity Maslov unit is also preserved. This checks both the no-excess constant and the reflected middle-frequency convention.

9. Exercises with complete solutions

Exercise 9.1 (find the radical; introductory). On Rx×Rθ3\mathbb R_x\times\mathbb R^3_\theta, take

Q=xθ1+12(θ2−aθ3)2,a∈R. Q=x\theta_1+\tfrac12(\theta_2-a\theta_3)^2, \qquad a\in\mathbb R.

Determine RR, the Lagrangian image, and the density-valued Gaussian (1.5) when the scalar amplitude is one.

Solution. The radical is R={(0,ar,r):r∈R}R=\{(0,ar,r):r\in\mathbb R\}, so e=1e=1. The critical equations give x=0x=0, θ2=aθ3\theta_2=a\theta_3, and the output covector is θ1\theta_1; hence λQ={x=0}\lambda_Q=\{x=0\}. Put w=θ2−aθ3w=\theta_2-a\theta_3 and r=θ3r=\theta_3. This change has determinant one. The θ1\theta_1 integral is 2πδ0(x)2\pi\delta_0(x), and the ww integral is (2π)1/2eiπ/4(2\pi)^{1/2}e^{i\pi/4}. With the original normalization (2π)−7/4(2\pi)^{-7/4}, the density-valued Gaussian is

(2π)−1/4eiπ/4δ0(x)∣dx∣1/2⊗∣dr∣. (2\pi)^{-1/4}e^{i\pi/4}\delta_0(x)|dx|^{1/2}\otimes|dr|.

Thus (1.5) is this Gaussian multiplied by ∫Rχ ∣dr∣\int_R\chi\,|dr|. A projection with radical coordinate r+bθ1+cwr+b\theta_1+cw gives the same integral: the rr integration is a translation. The quotient phase has two variables and its normalized Gaussian is (2π)1/4eiπ/4δ0(2\pi)^{1/4}e^{i\pi/4}\delta_0; their ratio is the factor (2π)−1/2(2\pi)^{-1/2} in (1.7).

Exercise 9.2 (vertical and horizontal reduction; intermediate). In a two-dimensional symplectic space, compare the positive factors in (4.1) for each of the following: Δ=λ=λ0\Delta=\lambda=\lambda_0; Δ=λ=H\Delta=\lambda=H, where HH is horizontal; and Δ=λ0\Delta=\lambda_0, λ={x=bξ}\lambda=\{x=b\xi\}, b≠0b\ne0.

Solution. The reduced symplectic space is zero-dimensional in all three cases. In the first, d=k=e=1d=k=e=1, giving (2π)−3/4(2\pi)^{-3/4}. The normalized input point mass is (2π)−3/4∫eixθdθ=(2π)1/4δ0(2\pi)^{-3/4}\int e^{ix\theta}d\theta=(2\pi)^{1/4}\delta_0; its quadratic reduction sets x=0x=0 and retains ∣dθ∣|d\theta|, with coefficient (2π)−3/4(2\pi)^{-3/4}. This is not ordinary evaluation of δ0\delta_0 at zero.

In the second, d=e=1d=e=1, k=0k=0, giving (2π)−1/4(2\pi)^{-1/4}. The normalized input is the constant (2π)−1/4∣dx∣1/2(2\pi)^{-1/4}|dx|^{1/2}; eliminating xx retains ∣dx∣|dx|, and the inverse half density of Δ\Delta gives exactly the stated factor. In the third, d=k=1d=k=1, e=0e=0, giving (2π)−1/4(2\pi)^{-1/4}. The reduced coefficient from (6.4) is

(2π)−1/4e−iπsgn⁡b/4∣b∣−1/2. (2\pi)^{-1/4}e^{-i\pi\operatorname{sgn}b/4}|b|^{-1/2}.

The determinant factor belongs to the geometric density map, and the remaining Fresnel factor is the Maslov map. These examples distinguish the two sources of normalization.

Exercise 9.3 (a coupled Gaussian; intermediate). Let n=2n=2, Δ\Delta be the ξ1\xi_1 axis, and

λ={x=Bξ},B=(2335). \lambda=\{x=B\xi\},\qquad B=\begin{pmatrix}2&3\\3&5\end{pmatrix}.

Reduce the normalized Gaussian with coefficient cc. Give its reduced graph matrix, coefficient and inverse density factor.

Solution. Here k=1,j=0,l=1k=1,j=0,l=1, B11=2B_{11}=2, and K=0K=0. The Schur complement is 5−9/2=1/25-9/2=1/2. Formula (6.4) gives

cΔ=c(2π)−1/4e−iπ/4/2. c_\Delta=c(2\pi)^{-1/4}e^{-i\pi/4}/\sqrt2.

The reduced distribution is

cΔ(2π)−3/4∫ei(yη−η2/4) dη ∣dy∣1/2, c_\Delta(2\pi)^{-3/4} \int e^{i(y\eta-\eta^2/4)}\,d\eta\ |dy|^{1/2},

and its symbol is tensored with ∣dξ1∣−1/2|d\xi_1|^{-1/2}. The factor 1/21/\sqrt2 is the exact-sequence density Jacobian for x1=2ξ1+3ξ2x_1=2\xi_1+3\xi_2; the phase e−iπ/4e^{-i\pi/4} and factor (2π)−1/4(2\pi)^{-1/4} come from Gaussian reduction.

Exercise 9.4 (stage the constraints; advanced). Suppose Δ1⊂Δ\Delta_1\subset\Delta are isotropic. Let di,ki,eid_i,k_i,e_i denote the constraint dimension, vertical-intersection dimension and plane-intersection dimension for the first reduction. Express the corresponding dimensions for the second reduction by Δ/Δ1\Delta/\Delta_1, and check that the normalization exponents add to that for direct reduction by Δ\Delta.

Solution. Write d=dim⁡Δd=\dim\Delta, k=dim⁡(Δ∩λ0)k=\dim(\Delta\cap\lambda_0), e=dim⁡(Δ∩λ)e=\dim(\Delta\cap\lambda). The second dimensions are d−d1,k−k1,e−e1d-d_1,k-k_1,e-e_1. For either plane, a vector of the second intersection lifts to a vector of Δ\Delta differing by an element of Δ1\Delta_1 from a vector in the original plane; its quotient is therefore the quotient of the original intersection by the first intersection. This proves the dimension statements. Then

d1−2k1−2e14+d−d1−2(k−k1)−2(e−e1)4=d−2k−2e4. \frac{d_1-2k_1-2e_1}{4} +\frac{d-d_1-2(k-k_1)-2(e-e_1)}4 =\frac{d-2k-2e}{4}.

The two radical-density factors combine to Ω(λ∩Δ)\Omega(\lambda\cap\Delta) by the density exact sequence, and the two inverse constraint half densities combine to Ω−1/2(Δ)\Omega^{-1/2}(\Delta). Compatible quadratic coordinates give the same Fresnel/delta integration, as in Section 6, so the phase factors compose as well. The arithmetic alone checks the normalization, rather than replacing that invariant phase argument.

Exercise 9.5 (an excess-one constant model; intermediate). Take Si=T∗RS_i=T^*\mathbb R, let LiL_i be the horizontal plane, and compose G1=L1×L2G_1=L_1\times L_2 with G2=L2×L3G_2=L_2\times L_3. Compute (7.3) on the coordinate half densities and verify it with constant Gaussian kernels.

Solution. The composition is L1×L3L_1\times L_3, and N=L2N=L_2, with coordinate yy and excess one. The geometric density map sends

∣dx dy∣1/2⊗∣dy dz∣1/2to∣dx dz∣1/2⊗∣dy∣. |dx\,dy|^{1/2}\otimes|dy\,dz|^{1/2} \quad\text{to}\quad |dx\,dz|^{1/2}\otimes|dy|.

The Gaussian map multiplies this by (2π)−1/2(2\pi)^{-1/2}. Each normalized no-phase kernel has constant coefficient (2π)−1/2(2\pi)^{-1/2}, since its base dimension is two. Their product integrated formally over the middle variable has coefficient (2π)−1∣dy∣(2\pi)^{-1}|dy|. The reduced normalized no-phase kernel has coefficient (2π)−1/2(2\pi)^{-1/2}; dividing yields exactly the additional factor (2π)−1/2∣dy∣(2\pi)^{-1/2}|dy|. Pairing the retained density with a compact χ(y)\chi(y) gives a finite value. Integrating it over all of R\mathbb R without a cutoff would diverge.

Exercise 9.6 (identity graph; advanced). For Si=T∗RqS_i=T^*\mathbb R^q, compose two identity relations. Verify the excess, positive density map and Maslov unit using their delta kernels.

Solution. Matching identity relations forces all three vectors equal, so N=0N=0. The graph half densities are the square roots of the symplectic volume ∣dx dξ∣|dx\,d\xi|, and the graph-case rule (8.2) sends their units to the same unit. The normalized phase (x−y)⋅θ(x-y)\cdot\theta gives δ(x−y)∣dx dy∣1/2\delta(x-y)|dx\,dy|^{1/2}, with unit critical amplitude. Its composition satisfies

∫δ(x−y)δ(y−z) dy=δ(x−z). \int\delta(x-y)\delta(y-z)\,dy=\delta(x-z).

This equality is interpreted by the proper identity operators, or by testing the two kernels successively; it does not multiply unrelated delta distributions on the same variable. No Fresnel factor is introduced by these linear phases, so the identity Maslov unit is preserved. This also checks that the canonical input covector is reflected when relating kernel Lagrangians to operator relations.

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