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 restored phase-space lesson, Sections 2, 5 and 7: symplectic bases, common complements, linear reduction, critical maps and canonical relations.
- The restored Gaussian-symbol lesson, Sections 2–5: Gaussian lines, critical densities, the integer Maslov cocycle and exact quadratic normalization. The restored tangent-zoom lesson, Section 5, supplies the singular Gaussian and delta formula. The Gaussian companion, Z2, proves smooth dependence through changes of Hessian rank.
- Quadratic stationary phase, Q1–Q6, proves the limiting integrals, Fourier inversion, diagonalization and Fresnel constants. Finite calculus and inverse maps supplies smooth matrix inversion; compact bumps and density integration, Appendix A.4–A.6, supplies the compact tests and exact-sequence density conventions.
- The existing canonical-composition component, C1–C4, and Maslov-composition component, G1–G5, retain their full alternative proofs of the matching sequence, quadratic presentations, coordinate invariance and associativity. Their geometric density map is unscaled; the Gaussian normalization here multiplies it by (2π)−e/2. Their frames sϕ, rather than the dimension-shifted frames eϕ, give the Gaussian phase convention used here.
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) denote the complex line of translation-invariant densities on a real vector space V. Its positive real ray determines its real powers Ωs(V), including s=±1/2. We use Ω(0)=C with positive unit one. The symplectic form is ω=∑dξj∧dxj, with vertical distinguished plane λ0={x=0}. Write
G(λ)=M(λ)⊗Ω1/2(λ)(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,θ) be any real quadratic form, with x∈Rn and θ∈F=RN. Set
R={r∈F:Qxθr=0, Qθθr=0},e=dimR.(1.1)
Thus {0}×R is the pure phase-variable part of the full Hessian radical. For every r∈R,
Q(x,θ+r)=Q(x,θ).(1.2)
The critical space and its image are
CQ={Qθ=0},λQ={(x,Qx):Qθ=0}.(1.3)
Proposition 1.1. The space λQ is Lagrangian. The critical map CQ→λQ is onto with kernel {0}×R; hence CQ/R→λQ is an isomorphism. On any complement F=W⊕R, the restriction Q∣Rn×W is a nondegenerate phase parametrizing this same plane.
Proof. The transpose of the linear map (x,θ)↦Qθ has kernel R: its components are exactly the two matrices in (1.1). Its rank is therefore N−e, so dimCQ=n+e. A vector in the kernel of the critical map has x=0, Qxθθ=0 and Qθθθ=0, hence lies in R. The image has dimension n.
For two vectors (v,w),(v,w)∈CQ, the critical equations are
Qθxv+Qθθw=0,Qθxv+Qθθw=0.
Symmetry of the Hessian gives
ω((v,Qxxv+Qxθw),(v,Qxxv+Qxθw))=0:
the Qxx terms cancel, and the remaining difference is zero by the two critical equations and symmetry of Qθθ. Thus the n-dimensional image is isotropic and is Lagrangian.
On a complement W, a pure phase radical vector would also be a vector of R, so is zero. This proves nondegeneracy of the restricted phase. Equation (1.2) shows that it has the same critical image. ∎
The integral of eiQ over all of F has an infinite volume factor when e>0. Its appropriate replacement is an element of
I(λQ,H)⊗Ω(R),(1.4)
where H is the horizontal reference. Fix Lebesgue measure dθ, a scalar a, and any linear projection T:F→R that is the identity on R. For χ∈Cc∞(R), consider
UQ[χ](x)=a(2π)−(n+2N)/4∫FeiQ(x,θ)χ(Tθ)dθ ∣dx∣1/2.(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 T and is the pairing of a well-defined translation-invariant I(λQ,H)-valued density on R with χ.
Proof. Use a splitting θ=w+r. The density exact sequence identifies dθ with a product μWμR; any reciprocal rescaling of those two factors leaves their product unchanged. By (1.2), the phase is independent of r. Since T(w+r)=Tw+r, translation invariance gives
∫Rχ(Tw+r)μR(r)=∫Rχ(r)μ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)
The first bracket belongs to the nonzero Gaussian line unless a=0, by Proposition 1.1 and the previous lesson. Changing the complement replaces each w by w+Lw, with Lw∈R, which leaves the phase unchanged. The density exact sequence accounts for its measure change. Therefore the bracket tensored with μR is independent of the splitting as well as of T. This proves the claim. One may justify the displayed integrations by inserting a damping factor in W, performing the compact radical integration first, then taking the proved Fresnel boundary value. ∎
With the same quotient measure, the normalized nondegenerate Gaussian for Q∣W would use N−e phase variables. Consequently (1.6) is that Gaussian tensored with μR, multiplied by
(2π)−e/2.(1.7)
Here we retain the normalization with the original count N. 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 S have dimension 2n, let Δ⊂S be isotropic, and set
d=dimΔ,K=λ∩Δ,e=dimK.
Define
SΔ=Δω/Δ,λΔ=(λ∩Δω)/K,λ0Δ=(λ0∩Δω)/(λ0∩Δ).(2.1)
The images implicit in these quotient formulas are Lagrangian in SΔ, by the proved linear reduction proposition. Their dimensions are n−d.
There is a canonical positive isomorphism
dΔ:Ω1/2(λ)⟶Ω1/2(λΔ)⊗Ω(K)⊗Ω−1/2(Δ).(2.2)
Proof. Put B=λ∩Δω. The two exact sequences give
Ω1/2(λ)≃Ω1/2(B)⊗Ω1/2(λ/B),Ω1/2(B)≃Ω1/2(K)⊗Ω1/2(λΔ).(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)
is well-defined and perfect. Its left kernel is B; its right kernel before quotienting is Δ∩λω=K. Both quotients have dimension d−e. Duality therefore identifies
Ω1/2(λ/B)≃Ω−1/2(Δ/K)≃Ω1/2(K)⊗Ω−1/2(Δ).(2.5)
Combine (2.3) and (2.5). The two factors Ω1/2(K) become the full density Ω(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π. 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)
Choose symplectic coordinates split as
x=(x′,y,z),ξ=(ξ′,η,ζ),dimx′=k,dimy=l,dimz=j,(3.2)
such that
λ0={x=0},Δ={x′=y=η=ζ=0}.(3.3)
Thus (ξ′,z) are coordinates on Δ,
Δω={x′=ζ=0},(3.4)
and (y,η) are symplectic coordinates on the quotient.
Such coordinates exist with the specified distinguished vertical plane. Choose a basis of Δ∩λ0, extend it to Δω∩λ0, then to λ0. The dimensions of these successive spaces are k,k+l,n: restriction of the symplectic pairing of Δ with λ0 has rank j. Choose j vectors of Δ dual to the last vertical basis vectors; they complete Δ∩λ0 to Δ and are mutually isotropic. Here is a completion that keeps those prescribed vectors fixed. Write the vertical basis as ei, and call the prescribed horizontal vectors fj, with indices in the last block J, so ω(ei,fj)=δij. For each remaining index choose a vector gi dual to all ea; this is possible since the pairing identifies S/λ0 with λ0∗. Replace gi by gi−∑j∈Jω(gi,fj)ej. It stays dual to the vertical basis and is now orthogonal to every prescribed fj. Finally replace each such gi by gi+21∑a∈/Jω(ga,gi)ea, using the vectors before this last replacement in the coefficients. The new pairing between indices i,b∈/J is ω(gi,gb)−21ω(gi,gb)+21ω(gb,gi)=0. This last correction preserves the pairings with prescribed vectors. Together with the ei and fixed fj, the resulting vectors form the required symplectic basis. This gives (3.2)–(3.3), including empty blocks.
Represent an element of G(λ) in the horizontal reference by a nondegenerate quadratic phase:
u(x)=a(2π)−(n+2N)/4∫eiQ(x,θ)dθ ∣dx∣1/2.(3.5)
Such a phase exists: complete the nonzero vertical constraints of λ, and write it as q(x)+x⋅Tθ with T injective and q quadratic, as in the previous lesson's signature proof.
Form the new quadratic phase
qΔ(y,z,θ)=Q(0,y,z,θ),(3.6)
with base variable y and phase variables (z,θ). Interpret
g(y)=a(2π)−(n+2N)/4∬eiqΔ(y,z,θ)dzdθ ∣dy∣1/2(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′=0 followed by an unweighted pushforward.
Lemma 3.1. The pure phase radical of (3.6) is naturally isomorphic to K, and its critical image is λΔ. Hence (3.7) belongs to
I(λΔ,HΔ)⊗Ω(K).(3.8)
Proof. Its critical equations are Qz=Qθ=0, evaluated at x′=0. In the original critical map, these are precisely x′=ζ=0, so the original covector lies in λ∩Δω. Projection to (y,η) gives the reduced plane in (2.1).
For a pure radical vector (z,θ), set x′=y=0. All derivatives of (3.6), including its y derivative, vanish. Thus its image in the original critical map has x′=y=η=ζ=0, and lies in λ∩Δ=K. Conversely every vector of K 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)
and pass to the intrinsic Gaussian symbol lines. This is the candidate map
RΔ:G(λ)⟶G(λΔ)⊗Ω(K)⊗Ω−1/2(Δ).(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,Δ) and λ, and is invariant under symplectic isomorphisms preserving those data. For fixed (S,λ0,Δ), on any smooth family of λ's with dim(λ∩Δ) constant, it is a smooth bundle map. Its positive density factor is
(2π)(d−2k−2e)/4dΔ,(4.1)
and its unit-modulus factor gives a canonical map
M(λ)⟶MΔ(λΔ).(4.2)
Proof of independence. First change the phase representing the same input Gaussian. Complete squares in the nonzero block of Qθθ. Eliminating h variables by the exact Fresnel formula multiplies the amplitude by
eiπsgnD/4∣detD∣−1/2
and changes the normalization from N to N−h, because the Gaussian factor (2π)h/2 supplies exactly that change. The same elimination, with the same factor, can be performed in (3.7); setting x′=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)
where T is injective. Its constraint subspace V=kerTT is the base projection of λ, and q∣V is determined by λ, since q(x)=ξ⋅x/2 there. Two such injective T's with the same range differ by an invertible phase-basis change. Two such q's differ by a quadratic polynomial vanishing on V. In linear coordinates V={x1=⋯=xs=0}, every monomial of that polynomial has a constrained coordinate as a factor. Consequently the difference is x⋅TLx for a linear map L, and can be absorbed by the translation θ↦θ+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, A=AT, preserving (3.3). This preservation requires A(0,0,z) to have only a ξ′ component. Hence the yz and zz blocks of A vanish. On x′=0, its quadratic form is independent of z and is yTAyyy. The input chirp e−ixTAx/2 in the preceding lesson therefore changes (3.7) by precisely e−iyTAyyy/2, the induced chirp on the reduced space. Its shear on Δ has determinant one. Thus (3.10) respects the intrinsic Gaussian-line identifications.
Finally keep the horizontal reference, and change dual bases by x=Pv, ξ=P−Tν, preserving the two flags x′=0 and x′=y=0. The matrix has the block form
P=P11P21P310P22P3200P33.(4.4)
At v′=0, the remaining base is y=P22vy, and the integrated variable is z=P32vy+P33vz. Pullback of the input half density and integration over vz give the factor
∣detP∣1/2∣detP33∣−1=∣detP22∣1/2∣detP11∣1/2∣detP33∣−1/2.(4.5)
The first factor is the reduced half-density Jacobian. The other factors are exactly the transformation of Ω−1/2(Δ): on Δ, the old coordinates (ξ′,z) are (P11−Tν′,P33vz). The radical density is transported to K 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)θ, with B(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=n. The matrix of the reduced critical equations has constant rank because its transpose kernel has dimension e, 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. Its critical map is injective, and a fixed nonzero minor selects a smooth basis for its image λΔ(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, 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 K 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-e family. At a jump of e, 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 Δ⊂λ, so e=d=k+j. In adapted coordinates, choose reduced horizontal coordinates making λΔ the graph y=Bη, B=BT. This is possible by a common complement to λ0Δ and λΔ, lifted to a horizontal complement preserving (3.3).
Then λ is represented by
Q(x′,y,z;θ′,η)=x′⋅θ′+y⋅η−21ηTBη.(5.1)
Indeed inclusion of Δ forces x′=ζ=0 on λ; reduction gives y=Bη, and the remaining ξ′ and z variables are free. Here the phase-variable count is N=k+l. Its critical density is
∣dθ′dηdz∣,
so the positive input symbol is ∣a∣∣dθ′dηdz∣1/2.
After setting x′=0, the variables (θ′,z) are the radical and identify with K=Δ. The remaining reduced Gaussian has coefficient a times
(2π)−(n+2k+2l)/4+3l/4=(2π)−(3k+j)/4.(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)
The exact-sequence map (2.2) has exactly these positive density frames without the scalar factor. Since
−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. We need a horizontal complement H transverse to both λ0 and λ, preserving the adapted form of Δ, whose induced reduced horizontal plane is transverse to λΔ. Here is a construction that keeps all those requirements.
Choose V⊂Δ complementary to Δ∩λ0. It is transverse to both λ0 and λ. Reduce first by V. In that quotient, the remaining constraint D=Δ/V lies in the distinguished vertical plane and is transverse to the reduced λ. In the further quotient by D, choose a common horizontal complement to the two reduced Lagrangians. Lift it into a symplectic decomposition of the quotient by V with coordinates (x′,y;ξ′,η), where D is the ξ′ plane and the selected horizontal plane in the further quotient is η=0.
For the reduced λ, its subspace η=0 projects bijectively to the x′ coordinates. To check this, a vector there with x′=0 projects to the intersection with the selected reduced horizontal plane, so its projection is zero; it then belongs to λ∩D=0. Also the η map on this Lagrangian has rank l: its transpose kernel is its intersection with the y-horizontal plane, which has just been seen to be zero. Therefore its kernel has dimension k, equal to the number of x′ coordinates. On this kernel, ξ′=Cx′ for a symmetric C, by isotropy. The plane ξ′=tx′, η=0 is transverse to it if C−tI is invertible; choose real t outside the finite eigenvalue set. This plane is a Lagrangian complement to the vertical plane, is transverse to the reduced λ, and induces the previously selected horizontal plane after reduction by D. Its preimage under reduction by V contains V and is a horizontal complement with all the stated properties in S.
With that choice, write
λ={x=Bξ},B=BT,ξ=(ξ′,η,ζ).(6.1)
The block B11 is invertible. Indeed, the constrained equations are ζ=0, (Bξ)′=0; a solution with η=0 gives a vector of λΔ in its horizontal complement, hence is zero there and then lies in K=0. Thus B11ξ′=0 implies ξ′=0.
The normalized input is
u(x)=c(2π)−3n/4∫ei(x⋅ξ−ξTBξ/2)dξ ∣dx∣1/2.(6.2)
In (3.7), integration over z gives (2π)jδ0(ζ). Completing squares in ξ′ then gives the reduced graph matrix
BΔ=B22−B21B11−1B12(6.3)
and the reduced Gaussian coefficient
cΔ=c(2π)(j−k)/4e−iπsgnB11/4∣detB11∣−1/2.(6.4)
All constants follow directly: the z integral contributes (2π)j, the k-variable Fresnel integral contributes (2π)k/2, and converting from the input coefficient (2π)−3n/4 to the reduced coefficient (2π)−3l/4 gives (2π)(j−k)/4. The negative signature is due to the Hessian −B11.
To compare with (2.2), change coordinates on λ from (ξ′,η,ζ) to (x′,η,ζ). The determinant is ∣detB11∣, so
c∣dξ′dηdζ∣1/2=c∣detB11∣−1/2∣dx′dηdζ∣1/2.(6.5)
The coordinates on λ/B are (x′,ζ). In the symplectic pairing (2.4) they are dual to (ξ′,z) on Δ, up to a sign in the first block, irrelevant to a positive density. The geometric map therefore sends (6.5) to
c∣detB11∣−1/2∣dη∣1/2⊗∣dξ′dz∣−1/2.(6.6)
Comparison with (6.4) proves (4.1), since j−k=d−2k and e=0. It also gives the explicit Maslov map in these graph frames: multiplication by e−iπsgnB11/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 K, reduce first by K⊂λ, then by Δ/K. The second reduced intersection is zero. If k1=dim(K∩λ0), the second vertical intersection has dimension k−k1: a vector of Δ/K in the reduced vertical plane lifts to a vector of (Δ∩λ0)+K. The two exponents already computed add to
4e−2k1−2e+4(d−e)−2(k−k1)=4d−2k−2e.(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⊂Δ. Choose a vertical basis starting with Δ1∩λ0, extend it through Δ∩λ0 and Δω∩λ0, and then to λ0. In Δ/(Δ∩λ0), start with the image of Δ1 and extend a basis; choose representatives of the first vectors in Δ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) and z=(s,t), so that Δ1 uses precisely the ξa and s coordinates, and Δ uses all of (ξ′,z).
The first reduction sets a=0 and treats (s,θ) as phase variables, with remaining base (b,y,t). Its radical R1 identifies with λ∩Δ1. Every vector of R1 annihilates the full Hessian after a=0, so it also annihilates the final restricted phase after b=0, where (t,s,θ) are phase variables. Hence R1⊂R, the final radical. Choose a complement to R1 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=0 and adds t to the phase variables. Its radical is exactly R/R1: the restricted phase is already independent of R1, 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), and these identifications coincide under the critical maps with the corresponding intersections in λ. 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(Δ). 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, the first intersection is all of K and the second is zero. The scalar computation (6.7) therefore completes Theorem 4.1 for every K.
7. Apply reduction to composition of linear canonical relations
Let Si be symplectic of dimension 2ni, with distinguished Lagrangian λi0, i=1,2,3. Let
G1⊂S1×S2,G2⊂S2×S3
be Lagrangian for ω1−ω2 and ω2−ω3, respectively. Their composition and excess space are
G=G1∘G2={(s1,s3):(s1,s2)∈G1,(s2,s3)∈G2 for some s2},
N={s2:(0,s2)∈G1, (s2,0)∈G2},e=dimN.(7.1)
All intersections here are linear and hence clean. The linear reduction proof shows that G is Lagrangian, and that the matching space
F={(s1,s2,s2,s3)∈G1×G2}
fits into 0→N→F→G→0.
Theorem 7.1 (linear symbol composition). There are canonical bilinear maps
M(G1)×M(G2)⟶M(G),(7.2)
Ω1/2(G1)×Ω1/2(G2)⟶Ω1/2(G)⊗Ω(N).(7.3)
The second map is (2π)−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,
take λ=G1×G2 and the middle diagonal constraint
Δ={(0,s2,s2,0):s2∈S2}.(7.4)
It is isotropic. Its orthogonal imposes matching middle vectors, so λ∩Δω=F; the quotient is S1×S3, the reduced distinguished plane is λ10×λ30, and K=λ∩Δ identifies with N. Also
d=2n2,k=n2.(7.5)
The density factor (4.1) is consequently (2π)−e/2.
Trivialize the residual factor Ω−1/2(Δ) by the positive symplectic volume
ρ2=∣ω2n2/n2!∣(7.6)
on S2≃Δ. 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π 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 G with a full density on the excess space N. This is the linear source of the fiber density integrated in analytic clean composition. A translation-invariant density on N=0 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 G1 is the graph of a symplectic isomorphism κ:S2→S1, so n1=n2. Then N=0, and
G2⟶G,(s2,s3)⟼(κs2,s3)(8.1)
is an isomorphism. Lift ρ2 to a density ρG1 on the graph. If d1=cρG11/2 and d2∈Ω1/2(G2), the positive map (7.3) is
(d1,d2)⟼c(8.1)∗d2.(8.2)
Proof. A vector of G1×G2 has the unique decomposition
(κs2,s2,s2′,s3)=(κ(s2−s2′),s2−s2′,0,0)+(κs2′,s2′,s2′,s3).(8.3)
The second term is in F; the first identifies the quotient by F with the graph G1. Pairing that quotient with the middle diagonal is the symplectic self-duality of S2, up to a sign, which does not affect positive densities. The contraction in (7.6) therefore contracts d1 with ρG1−1/2, leaving precisely c. There is no excess factor because e=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)⋅θ has normalization (2π)−n2. It represents δ(x−y)∣dxdy∣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, take
Q=xθ1+21(θ2−aθ3)2,a∈R.
Determine R, 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}, so e=1. The critical equations give x=0, θ2=aθ3, and the output covector is θ1; hence λQ={x=0}. Put w=θ2−aθ3 and r=θ3. This change has determinant one. The θ1 integral is 2πδ0(x), and the w integral is (2π)1/2eiπ/4. With the original normalization (2π)−7/4, the density-valued Gaussian is
(2π)−1/4eiπ/4δ0(x)∣dx∣1/2⊗∣dr∣.
Thus (1.5) is this Gaussian multiplied by ∫Rχ∣dr∣. A projection with radical coordinate r+bθ1+cw gives the same integral: the r integration is a translation. The quotient phase has two variables and its normalized Gaussian is (2π)1/4eiπ/4δ0; their ratio is the factor (2π)−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; Δ=λ=H, where H is horizontal; and Δ=λ0, λ={x=bξ}, b=0.
Solution. The reduced symplectic space is zero-dimensional in all three cases. In the first, d=k=e=1, giving (2π)−3/4. The normalized input point mass is (2π)−3/4∫eixθdθ=(2π)1/4δ0; its quadratic reduction sets x=0 and retains ∣dθ∣, with coefficient (2π)−3/4. This is not ordinary evaluation of δ0 at zero.
In the second, d=e=1, k=0, giving (2π)−1/4. The normalized input is the constant (2π)−1/4∣dx∣1/2; eliminating x retains ∣dx∣, and the inverse half density of Δ gives exactly the stated factor. In the third, d=k=1, e=0, giving (2π)−1/4. The reduced coefficient from (6.4) is
(2π)−1/4e−iπsgnb/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=2, Δ be the ξ1 axis, and
λ={x=Bξ},B=(2335).
Reduce the normalized Gaussian with coefficient c. Give its reduced graph matrix, coefficient and inverse density factor.
Solution. Here k=1,j=0,l=1, B11=2, and K=0. The Schur complement is 5−9/2=1/2. Formula (6.4) gives
cΔ=c(2π)−1/4e−iπ/4/2.
The reduced distribution is
cΔ(2π)−3/4∫ei(yη−η2/4)dη ∣dy∣1/2,
and its symbol is tensored with ∣dξ1∣−1/2. The factor 1/2 is the exact-sequence density Jacobian for x1=2ξ1+3ξ2; the phase e−iπ/4 and factor (2π)−1/4 come from Gaussian reduction.
Exercise 9.4 (stage the constraints; advanced). Suppose Δ1⊂Δ are isotropic. Let di,ki,ei 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, and check that the normalization exponents add to that for direct reduction by Δ.
Solution. Write d=dimΔ, k=dim(Δ∩λ0), e=dim(Δ∩λ). The second dimensions are d−d1,k−k1,e−e1. For either plane, a vector of the second intersection lifts to a vector of Δ differing by an element of Δ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
4d1−2k1−2e1+4d−d1−2(k−k1)−2(e−e1)=4d−2k−2e.
The two radical-density factors combine to Ω(λ∩Δ) by the density exact sequence, and the two inverse constraint half densities combine to Ω−1/2(Δ). 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∗R, let Li be the horizontal plane, and compose G1=L1×L2 with G2=L2×L3. Compute (7.3) on the coordinate half densities and verify it with constant Gaussian kernels.
Solution. The composition is L1×L3, and N=L2, with coordinate y and excess one. The geometric density map sends
∣dxdy∣1/2⊗∣dydz∣1/2to∣dxdz∣1/2⊗∣dy∣.
The Gaussian map multiplies this by (2π)−1/2. Each normalized no-phase kernel has constant coefficient (2π)−1/2, since its base dimension is two. Their product integrated formally over the middle variable has coefficient (2π)−1∣dy∣. The reduced normalized no-phase kernel has coefficient (2π)−1/2; dividing yields exactly the additional factor (2π)−1/2∣dy∣. Pairing the retained density with a compact χ(y) gives a finite value. Integrating it over all of R without a cutoff would diverge.
Exercise 9.6 (identity graph; advanced). For Si=T∗Rq, 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=0. The graph half densities are the square roots of the symplectic volume ∣dxdξ∣, and the graph-case rule (8.2) sends their units to the same unit. The normalized phase (x−y)⋅θ gives δ(x−y)∣dxdy∣1/2, with unit critical amplitude. Its composition satisfies
∫δ(x−y)δ(y−z)dy=δ(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
- [Hörmander III, §21.6] Lars Hörmander, The Analysis of Linear Partial Differential Operators III: Pseudo-Differential Operators, corrected second printing, Springer, 1994, §21.6 clean quadratic discussion and Theorems 21.6.6–21.6.7.
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).