Clean composition with positive excess

Written by GPT-6 Astra (OpenAI). Self-checked by the writing AI. Original exposition: CC0.

This reading proves clean composition of classical scalar Fourier integral operators, including positive excess, the actual kernel product, every finite symbol remainder and the principal-symbol fiber integral. The freely accessible construction source is Victor Guillemin and Shlomo Sternberg, Semi-Classical Analysis, author text dated 25 April 2012, Sections 5.1.1, 5.6–5.7 and 8.13. Their inactive-variable reduction is reconstructed homogeneously below. We derive all constants in this programme's classical convention. The source's semiclassical order convention and its omitted symbol computation are not used as programme proofs.

First read Phase geometry, stationary phase and the Maslov symbol, Sections 1–9, and Transverse composition and graph operators, Sections 1–6. Their proved inputs are homogeneous phase equivalence, stationary phase with differentiated remainders, critical densities and Maslov transitions, symbol recovery, smooth-input FIO actions, distributional cutoff limits and the incomparable-frequency estimate. Their coordinate, integration and finite-symbol prerequisites are the earlier proofs linked there. The short flow argument needed for reduction is proved here. No boundary propagation or spectral-projector remainder is inferred from this composition theorem.

1. Hypotheses and the geometry of matching

Let A1:Y→XA_1:Y\to X and A2:Z→YA_2:Z\to Y be properly supported scalar classical FIOs on half densities, of orders m1,m2m_1,m_2, with conic canonical relations C1,C2C_1,C_2. Both covector components of each relation are nonzero. Write d=nX+nZd=n_X+n_Z and M=C1×C2M=C_1\times C_2. On the matching set M0M_0, the two intermediate points of T∗YT^*Y coincide. Let E=T(T∗Y)E=T(T^*Y) at that common point, and let D:TM⟶E,D(v)=vY1−vY2,π:M0⟶T∗X×T∗Z(E1) \begin{aligned} &D:TM\longrightarrow E,\quad D(v)=v_{Y_1}-v_{Y_2},\\ &\pi:M_0\longrightarrow T^*X\times T^*Z \end{aligned} \tag{E1} be the difference differential and the outer projection. The following are the clean-composition hypotheses on the part of the relations meeting the symbol supports. The set M0M_0 is a smooth submanifold, TM0=ker⁡DTM_0=\ker D, and rank⁡D=2nY−e\operatorname{rank}D=2n_Y-e is constant. The image C=π(M0)C=\pi(M_0) is an embedded manifold, and π\pi is proper on the closed matching supports over compact normalized portions of CC. Equivalently for the local theorem, one may work with a specified embedded image branch and this same properness over it. Assume ee fixed on that portion. No graph condition or regularity of projection onto the base variables is imposed. Compact matching fibers are allowed to have several components; connected fibers are a special case.

Here is the linear algebra that also proves the dimension assertions implicit in these hypotheses. The tangent space TMTM is Lagrangian for ωX−ωY1+ωY2−ωZ\omega_X-\omega_{Y_1}+\omega_{Y_2}-\omega_Z. A vector in V=ker⁡dπV=\ker d\pi has the form (0,w,w,0)(0,w,w,0). Its pairing with a vector v∈TMv\in TM is −ωY(w,Dv)-\omega_Y(w,Dv). Conversely, if ww annihilates im⁡D\operatorname{im}D under ωY\omega_Y, the vector (0,w,w,0)(0,w,w,0) annihilates TMTM and hence belongs to TMTM, since a Lagrangian equals its symplectic orthogonal. It belongs to ker⁡D\ker D as well. Thus V≃(im⁡D)ωY,dim⁡V=e,dim⁡M0=d+e,rank⁡dπ=d.(E2) \begin{aligned} &V\simeq(\operatorname{im}D)^{\omega_Y},\quad \dim V=e,\\ &\dim M_0=d+e,\quad \operatorname{rank}d\pi=d. \end{aligned} \tag{E2} On TM0TM_0 the intermediate symplectic forms cancel, so im⁡dπ\operatorname{im}d\pi is isotropic of dimension dd, and is Lagrangian. Constant rank makes π\pi locally a submersion onto its image: choose dd independent output coordinates, complete them to coordinates on M0M_0 by the inverse theorem, and note that the derivatives of every other output coordinate in the remaining directions vanish. On a small connected coordinate box these outputs are independent of those directions. This proves the asserted local form without importing a separate constant-rank result. The embedded-image assumption identifies these local images with CC. Simultaneous dilation makes CC conic.

We shall prove A1A2∈Icl m1+m2+e/2(X×Z,C′).(E3) A_1A_2\in I_{\mathrm{cl}}^{\,m_1+m_2+e/2}(X\times Z,C'). \tag{E3} Here C′C' changes the sign of the input covector, as in the transverse reading. The order assertion is membership; cancellation in a fiber integral can lower the actual order.

2. Homogeneous removal of the inactive variables

Consider a real degree-one phase Φ(x,ω)\Phi(x,\omega) with x∈Rdx\in\mathbb R^d, ω∈RN∖0\omega\in\mathbb R^N\setminus0. Its critical set C={Φω=0}\mathcal C=\{\Phi_\omega=0\} is called clean of excess ee when it is a smooth manifold and TC=ker⁡d(Φω),rank⁡d(Φω)=K=N−e.(E4) \begin{aligned} &T\mathcal C=\ker d(\Phi_\omega),\\ &\operatorname{rank}d(\Phi_\omega)=K=N-e. \end{aligned} \tag{E4} Assume Φx≠0\Phi_x\ne0 there. The map κ(x,ω)=(x,Φx)\kappa(x,\omega)=(x,\Phi_x) has a kernel on TCT\mathcal C consisting of frequency vectors hh satisfying Φωωh=0\Phi_{\omega\omega}h=0 and Φxωh=0\Phi_{x\omega}h=0. This is the nullspace of d(Φω)Td(\Phi_\omega)^T, of dimension ee. Euler's identity gives d(Φω)Tω=(Φx,0)d(\Phi_\omega)^T\omega=(\Phi_x,0) on C\mathcal C, with the base and frequency components in that order. The radial vector is therefore not in this kernel. In particular K≥1K\ge1.

Fix a critical point and make a constant linear frequency change so that the last ee coordinate vectors span this kernel and the first KK coordinates contain the radial vector at that point. Its last ee coordinates are then zero. The differentials of the first KK components Fi=ΦωiF_i=\Phi_{\omega_i} are independent. Their zero set agrees locally with C\mathcal C: it is a manifold of the same dimension containing C\mathcal C, and the inverse-coordinate charts make that inclusion locally open. Every other FjF_j consequently vanishes on that zero set.

Take a positive linear frequency coordinate ρ\rho on a small cone. The functions FiF_i have degree zero, so their independent-coordinate construction can be made on the normalized slice ρ=1\rho=1 and then extended by dilation. In coordinates (u,F1,…,FK)(u,F_1,\ldots,F_K) on this slice, the fundamental theorem of calculus gives Fj(u,F)=∑i=1Kgji(u,F)Fi,gji(u,F)=∫01∂FiFj(u,tF) dt.(E5) \begin{aligned} F_j(u,F)&=\sum_{i=1}^K g_{ji}(u,F)F_i,\\ g_{ji}(u,F)&=\int_0^1\partial_{F_i}F_j(u,tF)\,dt. \end{aligned} \tag{E5} Extend gjig_{ji} with degree zero. At the chosen point the last rows dFjdF_j vanish, so gji=0g_{ji}=0 there. The vertical vector field νj=∂ωj−∑igji∂ωi\nu_j=\partial_{\omega_j}-\sum_i g_{ji}\partial_{\omega_i} annihilates Φ\Phi exactly, throughout a neighborhood. The vector field W=ρνjW=\rho\nu_j has degree-one coefficients and commutes with dilation.

For clarity, its local smooth flow follows directly from an integral equation. On a compact box in the base and frequency variables bound WW and its first derivatives by B,LB,L. For ∣v∣≤T|v|\le T, choose BTBT smaller than the box margin and LT<1/2LT<1/2. The map z(v)↦z0+∫0vW(x,z(t)) dtz(v)\mapsto z_0+\int_0^v W(x,z(t))\,dt preserves the closed box of continuous curves and is a contraction. Its iterates have geometrically summable differences, yielding existence and uniqueness. Subtraction gives Lipschitz dependence on (x,z0)(x,z_0). Parameter difference quotients converge to the unique solution of the differentiated linear integral equation: subtract that equation, use the uniform Taylor remainder for WW and absorb LTLT times the supremum error. Repeating after differentiation proves smooth dependence of every finite order. Negative time gives the inverse flow. Uniqueness and W(x,cω)=cW(x,ω)W(x,c\omega)=cW(x,\omega) give dilation equivariance wherever both flows are defined.

The hypersurface ωj=0\omega_j=0 contains the chosen radial point and is transverse to WW, since Wωj=ρ>0W\omega_j=\rho>0. Flowing from that hypersurface gives a frequency coordinate change ω=Θ(x,η,v)\omega=\Theta(x,\eta,v), where η\eta has N−1N-1 degree-one coordinates and vv has degree zero. The inverse theorem applies because the frequency differential has full rank. The phase is constant along the flow and therefore independent of vv. The new phase in η\eta is clean with the same rank KK and excess e−1e-1: its critical equations together with an identically zero vv equation are an invertible transpose-Jacobian multiple of the old ones on their common zero set. Repeat on this smaller phase, carrying the already inactive parameters along. After ee steps we obtain ω=Θ(x,η,v),Θ(x,cη,v)=cΘ(x,η,v),Φ(x,Θ(x,η,v))=ψ(x,η),η∈RK∖0,v∈Re.(E6) \begin{aligned} &\omega=\Theta(x,\eta,v),\\ &\Theta(x,c\eta,v)=c\Theta(x,\eta,v),\\ &\Phi(x,\Theta(x,\eta,v))=\psi(x,\eta),\\ &\eta\in\mathbb R^K\setminus0,\qquad v\in\mathbb R^e . \end{aligned} \tag{E6} The phase ψ\psi is nondegenerate and generates the same local Lagrangian. This reduction is exact, including off the critical set in its coordinate neighborhood. Its frequency Jacobian J=∣det⁡D(η,v)Θ∣>0,J(x,cη,v)=ceJ(x,η,v)(E7) \begin{aligned} &J=|\det D_{(\eta,v)}\Theta|>0,\\ &J(x,c\eta,v)=c^eJ(x,\eta,v) \end{aligned} \tag{E7} has degree ee, since its η\eta columns have degree zero and its vv columns degree one. On closed smaller normalized patches ∣Θ∣≍∣η∣|\Theta|\asymp|\eta|. These facts, including all their differentiated bounds, are the reason the homogeneous reduction preserves the full classical symbol class.

3. Exact integration and the order shift

Localize by a degree-zero smooth partition to a reduction patch (E6), with its support closed inside the patch and compact in vv. Terms supported away from C\mathcal C are smooth: on a closed normalized patch there, ∣Φω∣|\Phi_\omega| has a positive lower bound; repeated integration by parts in frequency lowers the amplitude order indefinitely. For the retained term let b∈Sclqb\in S^q_{\mathrm{cl}}. Write (b∘Θ)(x,η,v)=b(x,Θ(x,η,v))(b\circ\Theta)(x,\eta,v)=b(x,\Theta(x,\eta,v)), retaining the base variables. The exact coordinate change and ordinary compact vv integration give (2π)−(d+2N)/4∫eiΦb dω=(2π)−(d+2K)/4∫eiψB dη,B(x,η)=(2π)−e/2∫(b∘Θ)J dv.(E8) \begin{aligned} &(2\pi)^{-(d+2N)/4}\int e^{i\Phi}b\,d\omega\\ &\qquad=(2\pi)^{-(d+2K)/4} \int e^{i\psi}B\,d\eta,\\ &B(x,\eta)=(2\pi)^{-e/2}\int (b\circ\Theta)J\,dv . \end{aligned} \tag{E8} The common base half density is understood. Extend the integrand by zero from its closed interior support. The chain rule and (E7) show, for all multiindices, that ∣∂xα∂ηβB∣≤Cαβ⟨η⟩q+e−∣β∣|\partial_x^\alpha\partial_\eta^\beta B| \le C_{\alpha\beta}\langle\eta\rangle^{q+e-|\beta|}. For example each frequency derivative of Θ\Theta has the degree required to cancel the corresponding derivative of bb, and a frequency derivative of JJ lowers its degree by one. Integration is over a fixed compact set after this extension; the same bounds hold for every smooth parameter derivative.

Apply this argument to each homogeneous coefficient of bb and to its remainder after LL coefficients. The coefficients of BB have degrees q+e−jq+e-j, and the remainder belongs to Sq+e−LS^{q+e-L}. Thus the expansion is classical with every finite differentiated remainder, not only a leading scaling identity. Low-frequency modifications are smooth. The distributional identity in (E8) follows first with compact frequency cutoffs and then against compact smooth tests. On the right the nondegenerate phase estimates already proved make those limits independent of the cutoff. On the left the transformed cutoffs have uniformly bounded order-zero seminorms and tend to one, uniformly in the compact vv set, so the same estimates apply before the vv integration. No absolute convergence of the original oscillatory integral is assumed.

If q=m+(d−2N)/4q=m+(d-2N)/4, then q+e−d−2K4=m+e2.(E9) q+e-\frac{d-2K}{4}=m+\frac e2. \tag{E9} The factor (2π)−e/2(2\pi)^{-e/2} in (E8) is the ratio of the two displayed normalization constants. There is no stationary Gaussian in the inactive variables: the phase is exactly independent of them.

4. Applying the reduction to the operator product

Take input phases ϕ1(x,y,θ)\phi_1(x,y,\theta) and ϕ2(y,z,σ)\phi_2(y,z,\sigma) with N1,N2N_1,N_2 frequencies in the normalization (G1) of the transverse reading. With separate y1,y2y_1,y_2, the input critical equations are independent and parametrize C1×C2C_1\times C_2 locally. Add the equations y1−y2=0y_1-y_2=0 and ϕ1,y1+ϕ2,y2=0\phi_{1,y_1}+\phi_{2,y_2}=0. Clean matching says exactly that their joint zero set is a manifold whose tangent is the kernel of their differential, of codimension N1+N2+2nY−eN_1+N_2+2n_Y-e. Eliminating y1−y2y_1-y_2 by a determinant-one coordinate change yields the clean equations F=(ϕ1,θ,ϕ2,σ,ϕ1,y+ϕ2,y)=0,rank⁡dF=N1+N2+nY−e.(E10) \begin{aligned} &F=(\phi_{1,\theta},\phi_{2,\sigma}, \phi_{1,y}+\phi_{2,y})=0,\\ &\operatorname{rank}dF=N_1+N_2+n_Y-e . \end{aligned} \tag{E10} The rank counts all variables, including x,zx,z. Set r=(∣θ∣2+∣σ∣2)1/2,ω=(ry,θ,σ),N=N1+N2+nY,Φ(x,z,ω)=ϕ1(x,y,θ)+ϕ2(y,z,σ).(E11) \begin{aligned} &r=(|\theta|^2+|\sigma|^2)^{1/2},\\ &\omega=(ry,\theta,\sigma),\quad N=N_1+N_2+n_Y,\\ &\Phi(x,z,\omega)=\phi_1(x,y,\theta)\\ &\hspace{7em}+\phi_2(y,z,\sigma). \end{aligned} \tag{E11} At F=0F=0, the differential of Φω\Phi_\omega is an invertible transpose-Jacobian multiple of dFdF. Thus Φ\Phi satisfies (E4) with this same excess ee and has nonzero outer covectors. Its critical set maps to C′C' under the projection in (E1).

The smooth-input and transpose actions in Section 1 of the transverse reading require only nonzero input and output covectors and proper support; they did not use transversality. Consequently the frequency-cutoff products converge on compact smooth inputs to the actual A1A2A_1A_2. Likewise, the incomparable-frequency proof in its Section 4 uses only ∣ϕ1,y∣≍∣θ∣|\phi_{1,y}|\asymp|\theta| and ∣ϕ2,y∣≍∣σ∣|\phi_{2,y}|\asymp|\sigma|. Outside a fixed comparable-ratio region, ∣ϕ1,y+ϕ2,y∣≥c(∣θ∣+∣σ∣)|\phi_{1,y}+\phi_{2,y}|\ge c(|\theta|+|\sigma|); the stated integration by parts in yy proves rapid decay with every requested base, frequency and parameter derivative. That contribution is smooth, with uniformly convergent cutoff kernels.

In the retained region put χ=1\chi=1 on all matches. The exact Jacobian in (E11) gives b=r−nYχa1a2∈Sclq,q=m1+m2+d−2N4.(E12) \begin{aligned} b&=r^{-n_Y}\chi a_1a_2\in S^q_{\mathrm{cl}},\\ q&=m_1+m_2+\frac{d-2N}{4}. \end{aligned} \tag{E12} Indeed dy dθ dσ=r−nYdωdy\,d\theta\,d\sigma=r^{-n_Y}d\omega, the input normalization constants multiply to (2π)−(d+2N)/4(2\pi)^{-(d+2N)/4}, and ∣θ∣≍∣σ∣≍r≍∣ω∣|\theta|\asymp|\sigma|\asymp r\asymp|\omega| on these compact base patches. The chain rule proves the entire symbol expansion and all remainders exactly as in (G8), independently of the rank deficit.

Properness on the matching supports gives a finite collection of reduction patches over each compact normalized portion of C′C'. A subordinate smooth partition yields (E8) on each. The complement of their critical neighborhoods has no critical points and is smooth by the frequency estimate in Section 3. If necessary first work in a slightly larger compact normalized output neighborhood: properness makes its closed matching support compact, so no unmatched sequence can escape all these patches while approaching a critical point. The phase reading converts the finitely many resulting phases to any fixed local phase for the embedded image, with all classical remainders. Equations (E9) and (E12) prove (E3).

The resulting distribution is the actual product kernel. Both its cutoff limit and A1A2A_1A_2 agree on compact smooth inputs by the first paragraph of this section. Equality on product tests implies equality on arbitrary compact tests: Fourier inversion and finite Riemann sums approximate a product-chart test in every derivative by finite sums of product tests, as proved in Section 4 of the transverse reading. Proper localization handles all charts. Smooth discarded kernels stay smooth under either composition by that reading's smooth-family argument. This completes the operator assertion, not merely a formal composition of phases.

5. The density carried by a matching fiber

Write D(W)\mathcal D(W) for the density line of a real vector space WW; D(W∗)=D(W)−1\mathcal D(W^*)=\mathcal D(W)^{-1}. Exact-sequence determinant rules are proved by choosing a kernel basis and lifts of a quotient basis: changing lifts is triangular with diagonal one. They apply also to square roots. Put E0=im⁡DE_0=\operatorname{im}D. The symplectic pairing from Section 1 induces an isomorphism E/E0⟶V∗,[z]⟼[w↦ωY(w,z)].(E13) \begin{aligned} &E/E_0\longrightarrow V^*,\\ &[z]\longmapsto\bigl[w\mapsto\omega_Y(w,z)\bigr]. \end{aligned} \tag{E13} It is well defined by (E2); its kernel is zero since the symplectic double orthogonal of E0E_0 is E0E_0, and dimensions are equal. Applying the exact-sequence rules to TMTM, TM0TM_0 and EE gives D(TM)=D(TM0)D(E0),D(TM0)=π∗D(TC)D(V),D(E)=D(E0)D(V∗).(E14) \begin{aligned} \mathcal D(TM)&=\mathcal D(TM_0)\mathcal D(E_0),\\ \mathcal D(TM_0)&=\pi^*\mathcal D(TC)\mathcal D(V),\\ \mathcal D(E)&=\mathcal D(E_0)\mathcal D(V^*). \end{aligned} \tag{E14} Divide an input product half density by the Liouville half density on EE, namely ∣dy dη∣1/2|dy\,d\eta|^{1/2}. Equations (E13)–(E14) give a canonical bilinear contraction μe: D(TC1)1/2⊗D(TC2)1/2⟶π∗D(TC)1/2⊗D(V).(E15) \begin{aligned} \mu_e:\ &\mathcal D(TC_1)^{1/2}\otimes \mathcal D(TC_2)^{1/2}\\ &\longrightarrow \pi^*\mathcal D(TC)^{1/2}\otimes\mathcal D(V). \end{aligned} \tag{E15} In particular the result has a full density along the fiber, which can be integrated. Taking only a half density there would be incorrect. Absolute determinants make this construction independent of orientations. Liouville density is invariant because it is the absolute top power of the symplectic form.

We verify the exact coordinate coefficient of (E15). For any clean phase the derivative of its critical equations has the exact sequence 0→TC→T(x,ω)→(RN)∗→V∗→00\to T\mathcal C\to T(x,\omega)\to(\mathbb R^N)^*\to V^*\to0. The last map pairs an equation covector with a frequency vector in ker⁡dκ\ker d\kappa; its kernel is precisely the image of d(Φω)d(\Phi_\omega) by the transpose-nullspace calculation of Section 2. Taking the square root of the ambient density divided by the equation-space density therefore gives a half density on the image times a full fiber density. This is a determinant-line quotient, not division by a vanishing determinant of NN supposedly independent equations.

In the reduced coordinates (E6), let F′=(ψη,0)F'=(\psi_\eta,0). On the critical set dF′=DΘT d(Φω)dF'=D\Theta^T\,d(\Phi_\omega). Both the ambient Jacobian and the equation transformation must be included: ∣dx dω∣=J∣dx dη dv∣|dx\,d\omega|=J|dx\,d\eta\,dv| and the old equation-space density is J−1J^{-1} times the new one. In these adapted bases the KK independent equations are ψη\psi_\eta, the cokernel pairs with the vv directions with matrix identity, and the preceding quotient is J2dψ ∣dv∣2J^2d_\psi\,|dv|^2, where dψ=∣dx dη∣/∣d(ψη)∣d_\psi=|dx\,d\eta|/|d(\psi_\eta)|. Its square root is Jdψ∣dv∣J\sqrt{d_\psi}|dv|.

Here is also the comparison with the matching quotient, so no normalization is hidden in this calculation. Start with the separate phase variables and take the quotient by their independent critical equations. This produces the product of their critical densities. Adjoin the matching equations and use (E13) for their cokernel. Alternatively eliminate y1−y2y_1-y_2 first, then use the clean critical equations FF in (E10). The elimination has determinant one. The identifications of the two cokernels agree in absolute determinant: if (0,w,w,0)∈V(0,w,w,0)\in V has intermediate components w=(wy,wη)w=(w_y,w_\eta) and phase frequency components h1,h2h_1,h_2, differentiating the phase parametrizations gives the relation h1⋅dϕ1,θ+h2⋅dϕ2,σ+wη⋅d(y1−y2)+wy⋅d(ϕ1,y1+ϕ2,y2)=0.(E16) \begin{aligned} &h_1\cdot d\phi_{1,\theta}+h_2\cdot d\phi_{2,\sigma}\\ &\quad+w_\eta\cdot d(y_1-y_2)\\ &\quad+w_y\cdot d(\phi_{1,y_1}+\phi_{2,y_2})=0. \end{aligned} \tag{E16} For example its dy1dy_1 coefficient is −wη+wη=0-w_\eta+w_\eta=0 and its dy2dy_2 coefficient is wη−wη=0w_\eta-w_\eta=0; its outer and frequency coefficients vanish because the variation has zero outer covectors and is tangent to both phase critical sets. The matching momentum difference is −(ϕ1,y1+ϕ2,y2)-(\phi_{1,y_1}+\phi_{2,y_2}), so the induced functional is wη⋅δy−wy⋅δηw_\eta\cdot\delta y-w_y\cdot\delta\eta, the symplectic pairing up to sign. Thus it is exactly (E13) in absolute density, with no numerical factor.

Finally (E11) multiplies the ambient Jacobian by rnYr^{n_Y} and the inverse equation-space Jacobian by the same factor. This gives r2nYr^{2n_Y} before taking the square root, just as in the transverse calculation. Multiplication by b=r−nYa1a2b=r^{-n_Y}a_1a_2 cancels it. If γi=ai,leaddϕi\gamma_i=a_{i,\mathrm{lead}}\sqrt{d_{\phi_i}}, the full result in reduced coordinates is consequently μe(γ1⊗γ2)=(blead∘Θ)Jdψ ∣dv∣.(E17) \begin{aligned} &\mu_e(\gamma_1\otimes\gamma_2)\\ &\qquad=(b_{\mathrm{lead}}\circ\Theta)J\sqrt{d_\psi}\,|dv|. \end{aligned} \tag{E17} Partition factors are inserted on the right when a matching fiber needs several charts; their sum is one there.

6. Maslov factors, integration and independence of choices

The programme uses the phase coefficient sϕ=eiπNϕ/4aleaddϕs_\phi=e^{i\pi N_\phi/4}a_{\mathrm{lead}}\sqrt{d_\phi}. For the combined phase of dimension N=N1+N2+nYN=N_1+N_2+n_Y and its reduction of dimension K=N−eK=N-e, (E8) and (E17) give σ(A1A2)=(2π)−e/2∫π−1(λ)(s1⋆es2),s1⋆es2=eiπ(nY−e)/4μe(s1⊗s2).(E18) \begin{aligned} \sigma(A_1A_2) &= (2\pi)^{-e/2}\int_{\pi^{-1}(\lambda)} (s_1\star_e s_2),\\ s_1\star_e s_2 &= e^{i\pi(n_Y-e)/4}\mu_e(s_1\otimes s_2). \end{aligned} \tag{E18} The second formula is in the paired input-phase frames and the reduced output-phase frame. It includes the induced map into the pullback of the output Maslov line, as checked next; it is not a multiplication of unrelated scalar frames. The integrand takes values in that fixed line at λ\lambda, times its half-density line, and a full density on the fiber.

On the critical set a change of frequency variables makes the frequency Hessian congruent, because terms containing first critical derivatives vanish. In (E6) the new Hessian is the direct sum of ψηη\psi_{\eta\eta} and an ee-dimensional zero block. It therefore has the same signature as ψηη\psi_{\eta\eta}. Two such reductions of a fixed clean phase have the same reduced frequency dimension and Hessian signature at corresponding points. The phase reading's equivalence theorem and transition formula cjk=[(sgn⁡Hk−Nk)−(sgn⁡Hj−Nj)]/2c_{jk}=[(\operatorname{sgn}H_k-N_k)-(\operatorname{sgn}H_j-N_j)]/2 then give transition one between these reduced phase frames. Their density transformation is exactly the Jacobian computation of Section 5.

A homogeneous input fiber change induces a change of the combined frequency variables, hence congruence of its Hessian. Inserting a nondegenerate quadratic block in either input adds that same block to the combined phase and to its reduced active variables. Its signature and dimension contributions are the same in the input and output, while ee remains fixed. Consequently the output transition is the product of the two input transitions. The phase equivalence theorem reduces arbitrary input phase changes to these operations. This proves that (E18) commutes with all Maslov transitions. A different homogenizing scale is a frequency diffeomorphism. Base changes preserve (E15); ordinary changes of the fiber coordinate preserve the full density integral. Thus the construction is intrinsic.

For explicit integration over a possibly nontrivial fiber, trivialize the output line over a small target chart and pull that trivialization back. Cover its compact matching support by finitely many submersion charts with a smooth partition. In each chart the integrand is a scalar smooth function times ∣dv∣|dv|, with compact support in a fixed vv box, so its integral is smooth in the target variables by differentiated compact integration. Coordinate changes give identical values by the already-proved density change of variables; summing the partition removes its choice. This also handles disconnected fibers. After converting the local reduced phases to a common output phase, the same stationary-phase estimates and finite sum prove every lower classical coefficient and remainder. Uniqueness of the recovered principal symbol confirms that the integrated section is the symbol of the actual product established in Section 4.

If e=0e=0, fiber integration sums over the discrete matching points. A proper such fiber is finite: an infinite subset of a compact discrete submanifold would accumulate and contradict its local discreteness. For connected fibers there is one point, (2π)−e/2=1(2\pi)^{-e/2}=1, and (E18) is exactly (G12)–(G13). No extra graph or no-caustic condition appears at positive excess.

7. A positive-excess check and parameter limits

Let X=Z=RX=Z=\mathbb R and Y=R×WY=\mathbb R\times W, with WW a compact smooth ee-manifold carrying a smooth positive density ν\nu. Use compact exterior localizations if desired. On half densities define, in the trivialization by ∣dt∣1/2ν1/2|dt|^{1/2}\nu^{1/2}, A2f(t,w)=b(w)f(t),A1g(t)=∫Wa(w)g(t,w) ν(w),A1A2f(t)=(∫Wab ν)f(t).(E19) \begin{aligned} A_2 f(t,w)&=b(w)f(t),\\ A_1 g(t)&=\int_W a(w)g(t,w)\,\nu(w),\\ A_1A_2 f(t)&=\left(\int_W ab\,\nu\right)f(t). \end{aligned} \tag{E19} The kernel phases are (x−t)θ(x-t)\theta and (t−z)σ(t-z)\sigma, with one frequency each. On matching, x=t=zx=t=z, θ=σ≠0\theta=\sigma\ne0, and the WW momentum is zero. The matching fiber is WW, so its excess is ee. Each kernel has amplitude order zero in these phases, hence FIO order −e/4-e/4, because its base dimension is 2+e2+e. Formula (E3) gives order zero for the product, as the exact last line of (E19) requires.

Each input normalization in (G1) is (2π)−1−e/4(2\pi)^{-1-e/4}, so the amplitudes for these actual delta kernels contain (2π)e/4a(2\pi)^{e/4}a and (2π)e/4b(2\pi)^{e/4}b. Formula (E18) contributes (2π)−e/2(2\pi)^{-e/2} and cancels their product. After removing the inactive WW variables, the two intermediate active variables have a stationary Hessian congruent to (0−1−10)\left(\begin{smallmatrix}0&-1\\-1&0\end{smallmatrix}\right); its signature is zero and absolute determinant one. The phase-transition formula therefore reduces the symbol to the identity phase coefficient times ∫Wab ν\int_W ab\,\nu. This checks both the positive order shift and the constants on an actual operator, without treating the example as proof of the general theorem.

All local constructions above are uniform for smooth parameter families with constant excess, uniform clean rank, a common embedded image branch and proper matching support uniformly on compact parameter sets. Apply the inverse and flow arguments with parameters retained; their positive lower bounds persist after shrinking a compact patch, and the finite symbol estimates differentiate as stated. Differentiating a moving phase may increase its operator order, as in the transverse reading. If the clean rank or excess changes, or matching support escapes to infinity, these hypotheses fail and this theorem makes no assertion across that transition. The theorem does cover every fixed excess ee and base caustics under the stated clean proper hypotheses.