Proof. For every finite \(r\ge2\), the smooth orbit foliation is \(r\)-normally hyperbolic: \(d\phi_tX=X\circ\phi_t\) has uniformly bounded norm in both time directions, whereas the normal strong directions contract exponentially in the appropriate direction. Its tangent field is smooth, so the Lipschitz tangent-field hypothesis of Hirsch et al. (1970, sec. 3, pp. 1017–1018) applies. Their stable and unstable manifolds of orbit leaves form \(C^r\) laminations, hence give weak plaque charts continuous in the \(C^r\) topology. Compactness supplies uniform bounds on smaller plaques, including lower bounds for the differentials of their parametrizations.

Here the contact form recovers the strong plaques with the same finite-order control. Invariance and forward contraction give \(\alpha(E^s)=0\) and \(d\alpha(E^s,E^s)=0\); also \(\iota_Xd\alpha=0\). Thus \(\alpha\) restricts to a closed form on each weak stable plaque. In a \(C^r\) plaque chart its pullback \(\beta\) is \(C^{r-1}\). Its local primitive \(\tau\), normalized at the center, satisfies \(d\tau=\beta\) and is therefore \(C^r\) with uniform bounds. On the plaque, \(X\tau=1\). The uniform implicit function theorem makes its central level a \(C^r\) disk tangent to \(E^s\). A path in this disk has its length exponentially contracted by forward flow, so its points lie in the genuine strong stable leaf. Local uniqueness of the \(C^1\) stable manifold identifies the disk with the strong plaque. Reverse time for the unstable case.

The neighborhoods just obtained may initially depend on \(r\). Fix instead a sufficiently small uniform \(C^1\) strong plaque domain, written as a graph over its central tangent plane. For each \(r\), a fixed sufficiently long iterate sends this entire domain into the smaller uniformly controlled \(C^r\) plaque at the image center. Pull back that chart by the smooth inverse iterate. The chart differential and the fixed inverse iterate have uniform lower bounds; projection to the original tangent plane also has a uniformly invertible differential, by the fixed \(C^1\) graph bound. The implicit function theorem therefore gives uniform \(C^r\) bounds on the original domain. The iterate and bounds may depend on \(r\); the domain does not. These operations preserve continuous dependence in \(C^r\). Thickening the strong plaques by a fixed bounded flow interval gives the corresponding weak charts. This proves leafwise regularity and continuous leaf jets, without transverse differentiability. ◻

Proof. Choose a sufficiently long iterate so that the first derivative on the contracted plaques has norm at most \(\rho<1\), and use uniform plaque charts whose finite smooth norms are bounded. The first derivative of a composition of \(j\) plaque maps is bounded by \(\rho^j\). For the second derivative the chain rule gives a sum bounded by a constant times \[\sum_{i=0}^{j-1}\rho^{j-1-i}\rho^{2i} \le C\rho^{j-1}.\] At derivative order \(r\), the chain rule separates the term linear in the order-\(r\) derivative of the preceding composition from terms containing at least two derivatives of lower positive orders. Induction bounds the latter by \(C_r\rho^{2i}\) at time \(i\); the same geometric sum proves the assertion, after changing \(C_r\). Boundedly many omitted iterates only change constants. For weak plaques, use the chart \((y,t)\mapsto\phi_t(y)\) over a strong plaque with \(t\) in a fixed compact interval. The dynamics commutes with the flow and carries \(t\) unchanged. The strong-plaque estimates and smoothness of these bounded flow boxes prove the last assertion. ◻

Proof. The unstable Jacobian rate in the horizontal graph volume is \(J\). Changing graph volume changes that rate by a flow coboundary. Ruelle’s entropy inequality gives \(h_\mu(\phi_1)\le\int J\,d\mu\) for every flow-invariant probability measure, while Pesin’s formula gives equality for smooth Liouville volume (Ruelle 1978; Pesin 1977). The variational principle for pressure proves the first assertion.

Assume entropy equality. Then \(m_L\) is the equilibrium state both for \(-J\) and for the constant potential \(-h\), each with pressure zero. The flow is topologically transitive, since its ergodic Liouville measure has full support. Thermodynamic formalism (Bowen and Ruelle 1975) gives upper and lower Gibbs bounds for ordinary forward-time Bowen balls; we use the precise formulation in Climenhaga (2024, preprint version 2, Theorem 3.1 and Equation (1.5)). For a sufficiently small fixed Bowen radius valid for both potentials, apply both estimates to the same orbit ball centered on a periodic orbit of period \(\ell\), observed for \(k\ell\) units of time. Their comparison bounds \[\left|k\int_0^\ell(J(\phi_t v)-h)\,dt\right|\] by \(\log Q_J+\log Q_h\), where \(Q_J,Q_h\) are the two Gibbs constants, independent of the positive integer \(k\). Dividing by \(k\) proves that the integral vanishes on every periodic orbit. Livšic’s theorem (Livšic 1972) yields (2) with \(F\) Hölder and differentiable along the flow.

If \(x,y\) lie on one local strong unstable plaque, Hölder continuity and backward contraction give \[ F(y)-F(x) =\int_{-\infty}^{0} \bigl(J(\phi_t y)-J(\phi_t x)\bigr)\,dt. \tag{4}\] Every positive leaf derivative of the integrand decays exponentially by Lemma 3 and the uniform finite leaf bounds for \(J\). Differentiation under the integral proves the strong-unstable assertion. The flow equation supplies the flow derivatives and mixed weak-leaf derivatives. Uniform bounds of arbitrarily high order, combined with continuity of the functions and the uniform plaque charts, show by interpolation that every fixed leaf derivative varies continuously also when the plaque varies. This conclusion asserts continuity of the leaf jets, not transverse differentiability.

Let \(\operatorname{vol}_{v^\perp}\) be the oriented normal volume. Define a horizontal form by \[u_v\bigl((q_1,q'_1),\ldots,(q_m,q'_m)\bigr) =e^{F(v)}\operatorname{vol}_{v^\perp} (q'_1-Uq_1,\ldots,q'_m-Uq_m), \qquad \iota_Xu=0.\] It annihilates \(E^u\). Jacobi’s equation together with (1) gives \[\frac{d}{dt}(q'-Uq)=-U(q'-Uq).\] Taking determinants and using (2) proves the first transport identity. Let \(\iota(v)=-v\) be the flip and take \(w=\iota^*u\), up to a fixed nonzero normalization. Since \(\iota\phi_t=\phi_{-t}\iota\), it has the asserted opposite transport and leaf regularity. The two graph planes are transverse, so \(u\wedge w\) is nowhere zero on the contact bundle. Its ratio with \((d\alpha)^m\) is continuous and invariant. Ergodicity makes the ratio constant almost everywhere, and full support of \(m_L\) makes it constant everywhere. ◻

Proof. Write the section map as \((A(x,y),B(x,y))\) and its height shift as \(\rho(x,y)\). It carries the stable axis to the stable axis, so \(B(x,0)=0\) and \(D_xB(x,0)=0\). Symplecticity gives \[ D_yB(a,0)=D_xA(a,0)^{-t}. \tag{8}\] Indeed the derivative on that axis is block upper triangular, and preservation of the pairing between the two coordinate planes gives this identity. Put \[Y(y)=B(P(y),y),\qquad P_*(Y(y))=A(P(y),y).\] Then \(Y(0)=0\) and \(DY(0)=D_yB(a,0)\), independently of \(DP(0)\). Consequently formal inversion of the \(j\)-jet of \(Y\) uses only powers of the base-dependent invertible matrix \(D_yB(a,0)\) as denominators. Successively differentiating the displayed equations proves that the transformed graph jet through order \(j\) is polynomial in the hidden tensors, with smooth coefficients in \(a\).

For the coefficient, let \(M_*=DP_*(Y)\) and \(x=P(y)\). Pullback on the \(x\)-directions yields the determinant \(\det(D_xA-M_*D_xB)\). Our height convention says that the pullback of the target section form is \(e^{-hT+h\rho}\) times the source section form. Hence the transformed log record is given by \[\begin{align*} r_*(Y)={}&r(y)+c(x,y)-c_*(A,Y) \tag{9}\\ &-\log|\det(D_xA-M_*D_xB)|-hT+h\rho(x,y) +\text{constant}. \end{align*}\] At \(y=0\) the determinant is \(\det D_xA(a,0)\), independent of every hidden tensor. The right side has value zero for actual states by actual transport and the definition of the backgrounds. It therefore has value zero for every alternative with that base \(a\).

There is no missing \((j+1)\)st graph derivative in (9). In \(j\) derivatives of \(M_*D_xB(P(y),y)\), the term with all \(j\) derivatives falling on \(M_*\) is multiplied by \(D_xB(a,0)=0\). Every surviving term uses at most \(j-1\) derivatives of \(M_*\), hence at most \(j\) derivatives of \(P_*\). All derivatives of the logarithm have only powers of the nonzero base determinant as denominators. This proves closure at order \(j\) and polynomial dependence for the density jet.

This computation is a universal finite-jet formula: it applies to any smooth symplectic section map preserving the stable axis, with its smooth height shift. A composition of any length is again such a map. The maximal number of occurrences of the hidden tensors in this universal formula depends only on \(j\), not on the section map. This proves the asserted uniform degree bound for arbitrary compositions. Applying the formula to the inverse map proves invertibility and the same finite-level assertions for the inverse. For one time step all required section derivatives and inverse axis determinants are uniformly controlled on smaller compact charts. The universal formulas then prove all the asserted raw finite-order bounds.

It remains to identify the linear diagonal blocks. At the central state write \(S=D_xA(0,0)\) and \(U=D_yB(0,0)=S^{-t}\). The base derivative is \(S\). Modulo lower graph and density orders, the highest graph tensor transforms by \[p_j\longmapsto S\,p_j(U^{-1}\,\cdot,\ldots,U^{-1}\,\cdot),\] and the highest scalar tensor transforms by \[s_j\longmapsto s_j(U^{-1}\,\cdot,\ldots,U^{-1}\,\cdot).\] The latter may also receive a linear contribution from \(p_j\). Thus these quotient diagonal blocks are, respectively, \(S\otimes\operatorname{Sym}^jS\) and \(\operatorname{Sym}^jS\), with a triangular coupling. Their exponents are exactly the sums stated in the lemma. The multiplicative ergodic theorem applied to the invariant order filtration gives no other exponents. Since \(j\geq1\), every new weight is at least \(j\chi_*\). ◻

Proof. Use the measurable direct-sum norm \(\|\sum_i z_i\|_\oplus=\max_i\|z_i\|\) on \(E\). Only measurability of this norm is needed in the argument below; no tempered comparison with the original ambient norm is required. Write \(E=E_<\oplus E_r\) and put \(K=S\cap E_<\) and \(I=\operatorname{pr}_{E_r}S\). In quotient norms, \(S/K\) is the graph of a measurable linear map \(L:I\to E_</K\). Its transport obeys \[\|L_{f^n v}\| \leq C_\epsilon(v)^2 e^{(\lambda_{r-1}-\lambda_r+2\epsilon)n}\|L_v\|.\] Here the operator on the quotient has norm at most that on \(E_<\), and the inverse on \(I\) has norm at most the inverse on \(E_r\). Choose \(2\epsilon<\lambda_r-\lambda_{r-1}\). Thus \(\|L_{f^n v}\|\to0\) almost everywhere. For every \(a>0\), invariance of the probability measure and dominated convergence give \[\mu\{\|L\|>a\} =\int 1_{\{\|L_{f^n v}\|>a\}}\,d\mu(v)\longrightarrow0.\] Consequently \(L=0\) almost everywhere, so \(S=(S\cap E_<)\oplus(S\cap E_r)\). Induction proves the assertion. ◻

Proof. For \(\gamma<0\) the solution is \[q_n=e_n+\sum_{j\geq1} (L_{n-1}\cdots L_{n-j})e_{n-j};\] for \(\gamma>0\) it is \[q_n=-\sum_{j\geq0} (L_n^{-1}\cdots L_{n+j}^{-1})e_{n+j+1}.\] Both expressions solve (11) by shifting the summation index. To verify convergence and temperedness, choose \(0<3\delta<|\gamma|\). For example, in the first expression the \(j\)th term has norm at most \[C_\delta(v)E_\delta(v) e^{2\delta|n|}e^{(\gamma+3\delta)j},\] after slow regularization of the linear estimates and the input. The geometric series converges. Since \(\delta\) is arbitrarily small, its sum is tempered; the second expression has the same proof with \(-\gamma\) in place of \(\gamma\). Partial sums are measurable. The difference of two solutions is a homogeneous orbit. A nonzero such orbit grows exponentially in the future if \(\gamma>0\) and in the past if \(\gamma<0\), contradicting temperedness. ◻

Proof. Fix a finite level and suppress its index. The derivative \(A_v=DG_v(0)\) has a measurable Oseledets splitting by weights \(d\), with exponent \(-d\). Oseledets projections commute with every tower projection: the latter intertwines the derivative cocycles, and a vector with one exponent can project only to that exponent or to zero. Use adapted norms, temperedly equivalent to the raw norms, in which \(\|A_v\|\leq e^{-\chi_*/2}\), decreasing the exponent slack if necessary.

We first construct the formal conjugacy. Suppose it has been constructed through ordinary degree \(b-1\). The discrepancy in degree \(b\) is a homogeneous polynomial. Its component with output weight \(d\) and input weights \(d_1,\ldots,d_b\) lies in the cocycle of homogeneous maps with linear transport \[Q\longmapsto A_v Q(A_v^{-1}\,\cdot,\ldots,A_v^{-1}\,\cdot).\] This cocycle has the single exponent \(d_1+\cdots+d_b-d\). If it is nonzero, Lemma 9 solves the homological equation by a unique tempered correction. If it is zero, choose the correction to be zero and retain the discrepancy in the normal map. Induction gives measurable tempered coefficients at each finite order. There are no resonances after \(b>d_{\max}/d_{\min}\), so the resulting normal map \(N_v\) is a finite weight-preserving polynomial. Its linear diagonal is \(A_v\), and the weight induction just given supplies its polynomial inverse.

The construction respects the tower without choosing complements inside an equal-weight block. Project the degree-\(b\) equation to a lower level. Since raw projection is linear and intertwines the derivatives, this is precisely the lower-level homological equation. On each nonresonant part uniqueness in Lemma 9 identifies the projected correction with the lower correction; on each resonant part both corrections are zero and the projected discrepancy is retained. Induction proves compatibility of all formal changes and normal maps.

We now pass from a formal change to an exact smooth change. This step is included to specify the domains and finite-order estimates; it is the measurable contraction mechanism underlying the usual nonstationary normal-form results (Melnick 2019; Kalinin and Sadovskaya 2017). Let \(\Psi_v^{[b]}\) be the polynomial truncation of the constructed formal change through degree \(b\). Once \(b\) exceeds the largest resonant degree, the defect \[ R_v^{[b]}=\Psi_{fv}^{[b]}G_v-N_v\Psi_v^{[b]} \tag{13}\] vanishes through degree \(b\). On a fixed raw domain its \(C^q\) norms are tempered for every fixed \(b,q\). This follows directly from the raw bounds and the finitely many tempered polynomial coefficients.

Choose a positive radius \(r_v\) on which \(\|DG_v\|\leq e^{-\kappa}\) for a fixed \(\kappa>0\). The raw second-derivative bounds and tempered norm comparison allow \(r_v\) to be chosen tempered. Slowly decrease these radii so that \[r_{fv}\geq e^{-\delta}r_v,\qquad 0<\delta<\kappa.\] Then \(G_v(B(0,r_v))\subset B(0,r_{fv})\). On this forward admissible family of balls, \[ \|G_v^{(n)}(z)\|\leq e^{-\kappa n}\|z\|. \tag{14}\] For each fixed \(q\), the derivatives through order \(q\) of \(G_v^{(n)}\) have tempered bounds with arbitrarily small exponential slack. To see this directly, differentiate \(G_v^{(n+1)}=G_{f^nv}\circ G_v^{(n)}\). At derivative order \(q\) the term containing \(D^qG_v^{(n)}\) is multiplied by \(DG_{f^nv}\); all other terms are products of derivatives of lower order and a tempered fixed-order derivative of \(G_{f^nv}\). Induction on \(q\) and summation of the contracting scalar recurrence give the claim. One can obtain decay as well, but the stated subexponential upper bounds suffice below.

There is a finite exponent \(K\), depending on this state level but independent of \(b\), such that all coefficients of \(N_v^{(-n)}\) are bounded by \(C_\epsilon(v)e^{(K+\epsilon)n}\). Here \(N_v^{(-n)}\) denotes the inverse of \(N_v^{(n)}\) from the space at \(f^nv\) back to the space at \(v\). To check this bound, consider the finite-dimensional space of coordinate monomials of a given weight \(d\). Substitution by \(N\) is triangular by ordinary degree; every diagonal tensor block has the single exponent \(-d\). The triangular-array estimate therefore bounds forward substitution by \(e^{(-d+\epsilon)n}\) and its inverse by \(e^{(d+\epsilon)n}\), with tempered prefactors. Only weights up to \(d_{\max}\) occur in the coordinate functions of \(N^{\pm n}\), so one can take \(K=d_{\max}\), with additional slack. The same exponential bound controls any fixed derivatives of these polynomials on a fixed bounded set, since their degrees are bounded independently of \(n\).

Define \[ Q_{v,n}^{[b]}=N_v^{(-n)}\circ\Psi_{f^nv}^{[b]} \circ G_v^{(n)}. \tag{15}\] Consecutive terms differ by applying \(N_v^{(-n-1)}\) to two arguments whose difference is \(R_{f^nv}^{[b]}\circ G_v^{(n)}\). Both arguments and each of their fixed derivatives are bounded with arbitrarily small exponential slack; the arguments themselves converge to zero. Taylor’s formula for (13), (14), and the chain rule give, for \(b\geq k\), \[ \bigl\|R_{f^nv}^{[b]}\circ G_v^{(n)}\bigr\|_{C^k} \leq C_{b,k,\epsilon}(v) e^{[-\kappa(b+1-k)+\epsilon]n}. \tag{16}\] Indeed, a term with \(q\leq k\) derivatives falling on \(R^{[b]}\) contains its Taylor factor of order \(b+1-q\), evaluated at the exponentially contracting argument, and only a fixed finite number of derivatives of \(G_v^{(n)}\). All remaining costs are tempered. Using the integral difference formula for the polynomial \(N_v^{(-n-1)}\) and differentiating it through order \(k\) now yields \[ \|Q_{v,n+1}^{[b]}-Q_{v,n}^{[b]}\|_{C^k} \leq C_{b,k,\epsilon}(v) e^{[K-\kappa(b+1-k)+\epsilon]n}. \tag{17}\] The slack has been enlarged to absorb the finitely many chain-rule factors; it is still arbitrarily small for fixed \(b,k\).

Choose \(b_0\) so that \(\kappa(b_0+1)>K\). The series of differences converges in \(C^0\) on the forward admissible domain. For any fixed \(k\), choose a larger \(b\) with \(\kappa(b+1-k)>K\); it converges in \(C^k\) there. These limits agree. In fact \(\Psi^{[b']}-\Psi^{[b]}\) vanishes through degree \(b\) when \(b'>b\), and the same estimate with \(k=0\) shows \(Q_{v,n}^{[b']}-Q_{v,n}^{[b]}\to0\) whenever \(b\geq b_0\). Thus the common limit \(\Psi_v\) is smooth. Increasing \(b\) does not introduce a new convergence-domain assumption: each \(\Psi^{[b]}\) is a polynomial, and every fixed derivative of the raw maps is bounded on the original raw domains. On compact subsets of the forward interior, the preceding estimates apply directly. If a fixed-order auxiliary estimate uses a smaller tempered radius at future centers, (14) and slow variation eventually put the arguments inside it.

The identity relating \(Q_{v,n+1}^{[b]}\) and \(Q_{fv,n}^{[b]}\circ G_v\) proves (12) on taking limits. For any fixed Taylor order \(a\), take \(b\geq a\) sufficiently large for \(C^a\) convergence. Every telescoping difference has zero Taylor jet through degree \(b\), so \(j^a\Psi_v=j^a\Psi_v^{[b]}\). In particular \(D\Psi_v(0)=\mathrm{id}\). Summing (17) gives tempered fixed-order bounds. The inverse function theorem, first with the \(C^2\) bound to obtain a tempered injectivity radius and then differentiated to any fixed order, gives the stated inverse bounds. For example, shrink until \(\|D\Psi_v-\mathrm{id}\|\leq1/2\); inverse derivatives are then finite sums of products of derivatives of \(\Psi_v\) and inverse matrices of norm at most two. These estimates use only finitely many tempered quantities at every fixed order. Finally, projection commutes with each expression (15), after choosing a sufficiently large common truncation for the two levels. It therefore commutes with their germ limits. ◻

Proof. Fix a truncation \(b_0\) for which the \(C^0\) limit in (15) converges. By Lemma 6, at fixed raw \(a\) the map \(G_v^{(n)}\) has polynomial hidden degree bounded by a number \(D_j\) independent of \(n\), and its base component is independent of the hidden variables. The polynomial \(\Psi_{f^nv}^{[b_0]}\) has degree at most \(b_0\). The inverse normal iterate is weight preserving, so it has degree at most \(L_j=d_{\max}/d_{\min}\), rounded down. Consequently the hidden degree of every \(Q_{v,n}^{[b_0]}\) is at most \(D_jb_0L_j\). Polynomials of bounded degree form a closed finite-dimensional space under uniform convergence on a ball: evaluation at a fixed unisolvent finite collection of points recovers all coefficients by an invertible matrix. The \(C^0\) limit is therefore polynomial of this bounded hidden degree. Its coefficients are smooth in \(a\), either by the same finite evaluation formula on a smaller product neighborhood or by the already established smoothness of \(\Psi_v\).

Compatibility with level zero gives the triangular expression \[\Psi_v(a,\eta)=(\psi_v(a),F_v(a,\eta)),\] where \(\psi_v\) is a smooth local diffeomorphism and \(F_v\) is polynomial in \(\eta\). At fixed \(a\), its hidden derivative is invertible near the central state. The real analytic inverse function theorem in the hidden variables therefore gives an analytic solution \(\eta=H_v(a,\zeta)\) of \(F_v(a,\eta)=\zeta\). The ordinary smooth parameter inverse function theorem identifies this with a jointly smooth function. Moreover, differentiating the implicit identity in a base variable gives \[\partial_aH_v =-(D_\eta F_v)^{-1}\,\partial_aF_v,\] evaluated at \((a,H_v(a,\zeta))\). The right side is analytic in \(\zeta\) because \(F_v\) and every base derivative of \(F_v\) are polynomial in \(\eta\). Repeated implicit differentiation expresses each mixed derivative as a finite sum of products of such polynomials, lower derivatives of \(H_v\), and \((D_\eta F_v)^{-1}\). Induction proves the same analyticity for all fixed derivatives. Substituting \(a=\psi_v^{-1}(z_0)\) proves the statement in normalized base coordinates. If desired, a common small complex hidden neighborhood for each fixed real base point is obtained by polynomial complexification and implicit inversion; no complexification of the base is involved. ◻

Proof. The assertion about new coordinates is the quotient assertion in Lemma 6, preserved by the compatible grading in Lemma 10. Positivity bounds the ordinary degree of a monomial of weight at most \(D\) by \(D/\chi_*\), and only levels \(j\leq D/\chi_*\) can contribute independent variables. At these finitely many levels the dimensions are finite. Substitution is triangular by ordinary degree on each fixed-weight array; its diagonal consists of tensor products of the derivative blocks with total exponent \(-d\). The triangular-array estimates prove the asserted forward and inverse bounds. Tensor operations add or subtract the indicated exponents. Restriction preserves an upper bound and the inverse upper bound supplies the lower bound; the same two estimates pass to quotient norms. This proves the subquotient assertion. ◻

Proof. Let \(v'\) lie on the stable axis of \(v\). In raw coordinates the change of center is a smooth graph-and-density transformation \(C\); both charts use the same stable plaque. Its conjugate by \(\Psi_{v,j},\Psi_{v',j}\) is a smooth local map \(Q\) near the state representing the plaque through \(v'\). Take \(v'\) sufficiently close that this state is in a forward-admissible normal domain at \(v\). It is the zero state in the chart at \(v'\). Shrinking a neighborhood of it therefore places the neighborhood in forward-admissible domains for both descriptions.

Let \(Q_n\) be the corresponding change at \(f^n v,f^n v'\). Exact functoriality of graph transformation gives, on the chosen neighborhood, \[ Q=(N_{v',j}^{(n)})^{-1}\circ Q_n\circ N_{v,j}^{(n)}. \tag{18}\] At the future actual state, the raw chart changes have bounded smooth norms of each prescribed finite order. The normal changes and their inverses have tempered norms and radii along each of the two regular orbits. The arguments of \(Q_n\) in (18), minus its future actual center, decrease exponentially on a fixed smaller initial neighborhood. Consequently they eventually fit every prescribed finite-order tempered domain.

For a fixed integer \(b\), replace \(Q_n\) by its Taylor polynomial of degree \(b\) at the future actual center. Its error after substitution is \(O(e^{-\kappa(b+1)n+\epsilon n})\), for a fixed \(\kappa>0\) and arbitrarily small \(\epsilon>0\). On the bounded range in question, postcomposition by \((N_{v',j}^{(n)})^{-1}\) costs at most \(e^{Cn+\epsilon n}\), where \(C\) depends on the state level but not on \(b\). Choose \(\kappa(b+1)>C+1\). The resulting polynomial approximations converge uniformly to \(Q\).

Their degrees are bounded independently of \(n\): at a fixed state level both iterated normal maps and their inverses have bounded weighted, hence ordinary, polynomial degree, and \(b\) has been fixed. The vector space of polynomials of this bounded degree is closed under locally uniform convergence on an open set. Thus \(Q\) is polynomial. Applying the same argument to the inverse change gives a polynomial inverse. The two polynomial compositions equal the identity on an open set and therefore identically. Projection compatibility follows first for raw changes and then for the normal changes.

For a fixed flow shift, use its induced raw graph transformation. It commutes with \(f\), fixes the corresponding central states, and has fixed finite-order smooth bounds. The same argument applies; closeness of the two phase-space centers under iteration is unnecessary in this case. ◻

Proof. For each finite array, its evaluation span is measurable: evaluate the continuous actual graph on a countable dense set in successively smaller base balls and take the stabilized finite-dimensional span. Its annihilator gives the measurable space \(I_D\) of polynomial relations of ordinary degree at most \(D\). In an ambient space \(\mathbb C^N\), the dimension \(d\) of the closure is the largest size of a coordinate subset \(S\) for which \(I_D\cap\mathbb C[z_i:i\in S]\) is zero for every \(D\). Indeed, a maximal-dimensional irreducible component has a transcendence-basis subset among its coordinate generators. Conversely, if every component has a nonzero relation in the coordinates \(S\), the product of those finitely many relations lies in the full relation ideal. This dimension test uses only countably many finite-dimensional intersection ranks; the dimension–transcendence-degree identity is The Stacks Project Authors (2026, Tag 00P0).

The spans of the gradients at zero of the spaces \(I_D\) stabilize in \((\mathbb C^N)^*\). Choose the first degree attaining their maximal rank and then a basis by a fixed pivot rule; these are measurable choices. Every irreducible component contains zero: the connected group of positive weighted dilations fixes the finite set of components, and contraction to zero preserves each closed component. Thus the local dimension at zero is \(d\). The tangent-space criterion makes smoothness at zero equivalent to gradient rank \(N-d\); smoothness also forces irreducibility, since all components meet there (The Stacks Project Authors 2026, Tags 00KU, 0B8X, and 00NP). This proves measurability of the central smoothness and irreducibility condition without selecting components measurably. Projection ranks on the smooth germs are finite matrix tests. Nonzero minors and sufficiently small coordinate neighborhoods may then be chosen by countable tests. The components used below are only existential choices at a fixed center.

The dimensions of the closures are invariant under \(f\) and hence constant almost everywhere. Work in one such regular stable-axis chart. Zariski density of the actual graph supplies actual points outside the singular and intersection loci of the maximal-dimensional components. It supplies them in arbitrarily small neighborhoods, since the closure is the closure of a germ. Avoiding these closed algebraic exceptional loci gives a relatively open subset of the continuous actual graph. Its regular points have full conditional leaf measure. At any such point, polynomial re-centering maps its closure into the selected component. The almost-everywhere equality of dimensions makes this closure the whole component. Its center is smooth and its closure is irreducible.

Let \(G\) be the measurable invariant set where the central closure is irreducible and smooth at its center. To pass from the preceding leaf argument to ambient measure, suppose that \(G\) were null. Choose the regular center above so that its stable plaque meets \(G\) in conditional measure zero, as allowed by the conull convention in the preliminaries. The relatively open actual-state patch just constructed instead gives a positive leaf-measure set of regular centers in \(G\), a contradiction. Thus \(G\) has positive ambient measure; invariance and ergodicity make it conull.

For the projection assertion, proceed through the levels. The projection of the actual germ is the lower actual germ, so its image is Zariski dense in the lower irreducible variety. In characteristic zero the dominant projection has full rank on a nonempty open part of the smooth locus. The same actual-point and re-centering argument gives a positive leaf-measure patch of centers with full projection rank. Applying the same conditional-null contradiction to this invariant rank condition proves the central submersion assertion almost everywhere. These are the usual smooth-locus and generic-rank facts for affine varieties; see Hartshorne (1977) and The Stacks Project Authors (2026, Tags 0B8X, 07ND, and 02K4).

The complex formulation introduces no additional component problem. A selected real actual point outside the other components belongs to both its component and the complex conjugate of that component. The two components therefore coincide, so the selected component is defined over the reals. The selected point is nonsingular; its real locus is Zariski dense and locally smooth of real dimension equal to the complex dimension (Mangolte 2019, Theorem 2.2.9(2) and its proof). Density of the actual evaluations themselves comes from the closure argument above. Finally, irreducibility and equal dimension make the local inclusions supplied by polynomial re-centering equalities of germs. Intersect the resulting conull sets over the countably many levels and finite tests. ◻

Proof. The tangent spaces at zero are graded, and their projections are surjective. Choose graded linear retractions onto them, compatibly with the projections. The local graph over the tangent space is invariant under positive weighted dilation. Taylor expansion at zero then makes each output coordinate a weighted homogeneous polynomial: the Taylor remainder, pulled toward zero by dilation and divided by the fixed output weight, tends to zero once the Taylor order is sufficiently high. Choose successive graded complements for the tangent kernels. Their weights are among those of the new raw derivative blocks; the minimum such weight tends to infinity by the raw jet computation. Finally a monomial of ordinary degree \(q\) in parameters of weights at most \(D_0\) has weight at most \(qD_0\), proving the last assertion. ◻

Proof. At each regular source choose an evaluation basis on an arbitrarily small actual-state neighborhood inside a forward-admissible domain. Such a basis exists by the definition of the germ span. Countable dense tests on the actual graph give a measurable choice. Restrict to a positive-measure set where the domain size, basis norms, and inverse conditioning are uniformly controlled. Pull \(b\) back by the monomial substitution cocycle. Its norm is at most \(C_\epsilon e^{(d+\epsilon)n}\) by the single-exponent estimate on its weight block. At negative returns to this set, express it in the chosen basis and transport the resulting identity forward.

Every basis state is actual, and its transport is an actual state in the target chart. Forward-admissibility and the positive minimum weight give exponential convergence to zero. Projection surjectivity follows because evaluation vectors project to evaluation vectors on the same actual germ. Equation (19) gives the assertion weight by weight. Larger cutoffs may alter the conditioning constant and return set, but not the exponent \(d\) of the desired vector. ◻

Proof. Extend the function smoothly off \(V_j\) if necessary. Fix a weight \(d\) in its Taylor series, and take a Taylor polynomial of ordinary degree \(b\) large enough that its remainder on transported states is \(O(e^{-(d+1)n})\). This is possible because all normal variables decrease at some fixed exponential rate. Include every monomial of this polynomial in the evaluation array, and apply Lemma 17 to a pure weight-\(d\) vector. All other weights are killed exactly. The coefficient cost is \(e^{(d+\epsilon)n}\), whereas the smooth remainder tends to zero faster. Thus the weight-\(d\) Taylor polynomial annihilates the whole evaluation span and vanishes on \(V_j\). Repeat for the finitely many weights needed at each ordinary order.

For the second assertion first use continuity to include all actual states on the local graph. Re-center at any nearby typical actual point using Lemma 14; its closure is the same variety germ by Lemma 15. The first assertion makes the function flat there in every tangent direction. Its restriction to the fiber is analytic, hence zero on that connected fiber patch. Typical base labels are dense and have full conditional measure. Continuity therefore gives zero throughout a product neighborhood. The argument can be applied separately to any fixed coefficient. ◻

Proof. Use the paired charts and density backgrounds of Section 3. Transport one actual unstable plaque from a negative return and one actual stable plaque from a positive return. Write their equations as \(x=P(y)\) and \(y=P^-(x)\), with log-density records \(r(y),r^-(x)\). At their unique nearby intersection the actual identity is \[ e^{r(y)+r^-(x)}\det\bigl(I-DP(y)DP^-(x)\bigr)=1. \tag{20}\] The two backgrounds cancel at the same point. Base labels and slopes tend exponentially to zero, so the implicit intersection derivative is uniformly invertible.

Fix a desired pair of weights \((d,d')\). First truncate spatial graphs and density records at a large order \(R\). Actual plaque derivatives have uniform bounds of every fixed order. Spatial truncation changes the intersection and the left side of (20) by at most \(C_R(|x|+|y|)^{R-1}\); the loss of one derivative accounts for the slopes. Since the return lengths \(n,n'\) are comparable, this is smaller than \(e^{-(d+1)n-(d'+1)n'}\) when \(R\) is sufficiently large. The truncated intersection calculation is a smooth function of finitely many state coordinates. Taylor-expand it to sufficiently high ordinary order that the new remainder has the same bound. All low biweight components stabilize as these orders increase, because the formal intersection displacements have zero constant term and every state weight is positive.

Take product evaluation arrays containing every term of this finite polynomial. Pure-weight extraction on the two arrays has coefficient cost at most \(e^{(d+\epsilon)n+(d'+\epsilon)n'}\). The exact identity on actual plaque pairs and the two remainder bounds force its \((d,d')\) component to vanish on the product of the two varieties. This proves the formal identity at every finite order.

Now restrict to the two zero-base fibers. Their intersection equations have the solution \(x=y=0\) identically, so only the two slopes and prescribed zero-order densities remain. Arbitrary low-level fiber germs lift to every higher finite level used in the test. Fiber analyticity upgrades the formal identity to an identity near the actual pair. Finally bilinearity gives the same constant on affine combinations. ◻

Proof. Put \(q=\lfloor W/w_0\rfloor\). The parameter-zero term is the actual central graph, kept exactly. In the graph-parameter equation solve the constant term first. Each successive parameter degree is found by inverting the same central graph Jacobian and substituting already known lower degrees. At degree at most \(q\), Taylor substitution uses at most \(q\) derivatives in the source spatial variable. Taking \(k\) target derivatives adds finitely many derivatives; the density Jacobian adds one. For example, any \(r\ge k+2q+4\) suffices if all graph, inverse, and density operations are differentiated together.

The composed section map, its height function, and the inverse central graph parametrization have derivatives through every fixed order bounded by \(e^{A_r l}\) for some \(A_r\). This follows by induction from bounded one-step derivatives and the chain rule; the number of derivatives is fixed before \(l\) varies. Substitution of \(N^{-l}z\) on the weight-\(W\) array has norm at most \(e^{(W+\epsilon)l}\), by its triangular single-weight estimates. Every coefficient produced by the preceding nilpotent calculation is a fixed finite sum of products of these quantities, input derivatives, and central inverse Jacobians. The number of factors is bounded in terms of \(W,k\). Differences are estimated by replacing one factor at a time. Input bounds are tempered, so, after choosing smaller exponential slack for their finitely many factors, all these costs are bounded by \(e^{Cl+\epsilon l}\). The spatial Taylor degree does not enter the number of parameter operations or the required derivative order \(r\). ◻

Proof. Fix \(W\). Work first in the ambient normal-state tower and restrict to the varieties afterward. By Corollary 12, choose a level \(j_0\) after which every new kernel weight exceeds \(W\). For a fixed \(j\ge j_0\), let \(\pi:\mathcal E_j\to\mathcal E_{j_0}\) be the projection and let \(\Pi_j,\Pi_{j_0}\) be the spectral projections onto weights at most \(W\). If \(\mathcal R_j\) is a bounded raw right inverse of \(\pi\), then \[\pi\Pi_j\mathcal R_j=\Pi_{j_0}\pi\mathcal R_j=\Pi_{j_0}.\] Since \(\ker\pi\) has no weight at most \(W\), the restriction of \(\Pi_j\mathcal R_j\) to \(\operatorname{im}\Pi_{j_0}\) is the inverse of the low-weight projection. Its norm is tempered. Thus the stabilized low-weight inputs supply arbitrarily high spatial jets with tempered bounds at each fixed higher level, without a bound on a measurable choice of parameters for \(V_j\).

At time \(-l\), take the spatial Taylor polynomial of degree \(R\) of the prescribed variation, modulo weights greater than \(W\), and push it to the target, using \(N^{-l}z\) as source parameter. Keep the parameter-zero central graph and its density exactly, rather than Taylor-truncating them. Denote the resulting approximation by \(\Phi^{R}_{v,l}\). The preimage of a fixed small target central patch has spatial size at most \(C e^{-\gamma l}\), for a fixed \(\gamma>0\).

Compare consecutive approximations at source time \(-l\), first applying one step to the earlier source polynomial. Exact raw state dynamics makes their spatial Taylor jets agree through degree \(R\). In particular, the crossing density transformation does not require the next shape derivative. All higher derivatives needed for this comparison are tempered for fixed \(R,W\). Taylor’s formula at the contracted source argument and Lemma 20 give \[ \|\Phi^R_{v,l+1}-\Phi^R_{v,l}\|_{C^k} \le C_{R,W,k,\epsilon}(v) e^{-\{\gamma(R+1-r(W,k))-C(W,k)-\epsilon\}l}. \tag{21}\] Choose \(R\) so that the exponent in braces is positive. The differences are summable in \(C^k\). Comparing two spatial cutoffs whose jets agree through the smaller cutoff gives the same estimate. Thus increasing \(R\) after obtaining \(C^0\) convergence does not change the limit. One may then increase \(R\) again to obtain convergence in each \(C^k\). The limits for different \(W\) are compatible by the same uniqueness estimate.

Every prescribed finite spatial jet is retained by taking \(R\) sufficiently high; convergence in that derivative order proves Taylor matching. Shifting the source sequence by one step proves exact formal equivariance. Finally, all estimates in (21) hold at shifted regular centers with tempered prefactors. Summing the geometric series after reducing the slack shows that each fixed coefficient norm of the limit is tempered. Only a fixed weight cutoff and a fixed spatial norm are being bounded at a time. ◻

Proof. Let the source be at time \(-n\) and consider only times from \(-\lfloor n/2\rfloor\) to zero. On every fixed finite normal array, a monomial of weight \(d\) in the transported state is bounded by \[ C_\epsilon e^{-dn/2+\epsilon n} \tag{22}\] throughout this last half. Here and below the finite arrays are fixed before taking \(n\) large, and slack is divided among their finitely many coefficients. Slow regularization bounds every fixed tempered conversion factor uniformly by \(C_\epsilon(v)e^{\epsilon n}\) on this interval.

We first establish a common geometric patch. The actual crossings remain on the stable axis. At a crossing, the inverse unstable parameter differential over a block is the transpose of the contracting stable-axis differential, because \(D_xB(x,0)=0\) and symplecticity gives \(D_yB(x,0)=D_xA(x,0)^{-t}\). Choose a fixed block length making this inverse norm less than \(1/4\), after the bounded chart conversions. Uniform finite derivatives of actual plaques and of this fixed block map extend the inverse contraction to a fixed small parameter ball with norm less than \(1/2\). During the last half the crossings tend uniformly to zero, so the same ball is available for all sufficiently large \(n\). Shorter intermediate blocks have bounded cost. A fixed small final patch therefore has true inverse arguments of size \[ |y_{-\lfloor n/2\rfloor}|\le C e^{-\gamma_1 n} \tag{23}\] at halfway, for a fixed \(\gamma_1>0\).

Fix the low-order comparison domain before choosing \(W\). It consists of graph jets through \(k+2\) and density jets through \(k+1\) in a compact neighborhood of the bounded central jets, with graph Jacobians bounded away from singularity. The true jets lie there after shrinking the patch. Let \(L_k\ge2\) bound the one-step jet transformation and its first derivatives on a slightly larger such domain. This number uses only the fixed low orders and the uniform raw geometric bounds.

For each fixed \(W\), every positive-parameter term of the truncated profile tends uniformly to zero in these low norms during the last half: its spatial coefficient is tempered, while its state monomial has positive weight and satisfies (22). Thus the substituted profiles lie in the same comparison domain for all \(n\ge n_0(W)\). This is the decisive point: the threshold depends on \(W\), but the amplification constant \(L_k\) does not.

Compare the actual and substituted profiles at halfway. Choose a spatial order \(S\) and a finite tower containing their spatial jets through \(S\). Expand the smooth inverse normal change to ordinary degree \(b\), where \((b+1)w_0>2W\). All terms of weight at most \(W\) agree by Lemma 21. Every remaining polynomial term has weight greater than \(W\) and is bounded using (22). The ordinary Taylor remainder satisfies the same bound by the choice of \(b\). Consequently the differences of spatial Taylor coefficients through order \(S\) are at most \(C e^{-Wn/4}\), after sufficiently small slack. The spatial Taylor remainder at the true inverse argument is bounded by \(C e^{-\gamma_1(S-k)n+\epsilon n}\), using (23), uniform actual plaque derivatives, and tempered profile coefficients. Choose \[ \gamma_1(S-k)>W/2+1. \tag{24}\] The initial error in the graph/density norm at halfway is then at most \(C e^{-Wn/4}\).

At every last-half step, the substituted truncated profiles satisfy the exact graph/density recurrence modulo weights greater than \(W\). For fixed \(W\) and each fixed parameter derivative order, the real graph transformation of the finite polynomial profile is defined on a temperedly small parameter ball and has tempered derivatives there. Indeed its central graph Jacobian is uniformly invertible; the finite profile coefficient bounds give a tempered radius on which it remains invertible, and implicit differentiation uses only finitely many tempered coefficients. Exponentially small transported states eventually lie in these balls, uniformly over the last-half interval. Taylor-expand that smooth recurrence to ordinary degree \(b\) with \((b+1)w_0>2W\). Formal equivariance kills its terms through weight \(W\), and (22) bounds all remaining terms and the smooth remainder by \(C e^{-Wn/4}\). The same estimate holds in the fixed low derivative norms. The constant and the large-\(n\) threshold may depend on \(W,S\), which have now been fixed.

These errors can be compared at the true successive arguments. The approximate transformed argument differs from the true one by a bounded multiple of the source point error. Moving its jets to the true argument costs at most another bounded factor, because one additional spatial derivative lies in the fixed comparison domain. Increasing \(L_k\) once, if necessary, the successive errors \(E_i\) satisfy \[E_{i+1}\le L_k E_i+C e^{-Wn/4}.\] There are at most \(n\) steps. Hence the final error is at most \[ C n L_k^n e^{-Wn/4} \le C' e^{-\{W/4-\log L_k-1\}n}. \tag{25}\] Finally choose \(W>4(K+\log L_k+2)\), then \(S\) by (24), then the requisite finite level and ordinary Taylor degrees. This order of choices proves the result. All estimates were uniform on the fixed final subpatch. ◻

Proof. Fix the arrival chart and a polynomial relation of its finite-level variety. Re-centering is a smooth graph transformation followed by an implicit equation for the crossing. At parameter zero it describes the actual central unstable plaque, and its crossing derivative is invertible. Substituting \(\Phi_v\) therefore gives a well-defined parameter series with smooth spatial coefficients.

Fix one weight \(d\) of this series. Evaluate on actual states transported from negative returns as in Lemma 17. At the arrival point the transported actual plaques satisfy the chosen polynomial relation exactly. By Lemma 22, the substituted profiles approximate all graph and density jets needed for the test to accuracy \(e^{-Kn}\) for arbitrary prescribed \(K\). The arrival chart is fixed, so re-centering and polynomial evaluation have fixed smooth comparison constants. Taylor-expand the resulting composition in positive-weight variables through sufficiently high ordinary order, then retain weights through a large cutoff. On transported states its discarded terms again decrease as fast as required, by the positive minimum weight and (22) at the final time.

Choose \(K>d+1\) and take an evaluation array containing every term of this finite polynomial. Pure-weight extraction has cost \(e^{(d+\epsilon)n}\), so its weight-\(d\) component vanishes on the departure variety. Increasing the level is harmless by the surjectivity in Lemma 17. The same test, applied to the re-centered zero-order log-density minus its prescribed background, proves that equality too; omitting this entry from the state causes no loss. Repeat at every finite weight and arrival level.

On any fixed homogeneous parameter section, an ordinary coefficient uses only bounded weights by Lemma 16. Thus all compositions in the statement are ordinary formal compositions at each finite test. No analyticity in the spatial variable, convergence in parameters, or regular dependence of the arrival normal chart on the center has been used. ◻

Proof. Fix a spatial Taylor order. Each fixed \(y\)-order in (26) uses only finitely many shape and density jets, including one additional shape derivative for the conormal slope. Increasing the number of derivatives in \(t_j\) does not increase that required state level: it differentiates the same functions on that level.

For any \(d+1\) multi-indices \(\nu_0,\ldots,\nu_d\), including zero, form the wedge of the columns \(\partial_{t_j}^{\nu_i}S\) while keeping \((x,y)\) fixed. Substitute \(x=P_{a,t_j}(y)\) afterward. Every spatial Taylor coefficient of this wedge is a function on a finite \(V_j\), smooth in \(a\) and analytic in \(t_j\). We claim that it vanishes.

First center the calculation at a nearby typical actual state. Polynomial stable-axis re-centering identifies its variety germ with the present finite-level germ. In the re-centered chart choose a homogeneous tower section over a level containing all jets needed for the tested spatial coefficient. This section covers the entire finite-level parameter germ after the invertible chart comparison. Higher choices do not affect the tested spatial coefficient at displacement zero.

Apply Lemma 21 in this re-centered chart. At each small finite displacement on its actual unstable plaque, the parameter-zero graph is actual. The derivative of its crossing with respect to the base parameter is invertible near displacement zero, so the fixed-point constraint can be solved there as a formal series in the finite parameters. At typical displaced centers, Lemma 23 places all these constrained variations in the arrival centered fibers. Their augmented-form derivatives belong to the centered \(d\)-dimensional span. Every \((d+1)\)-fold wedge therefore vanishes.

Each such derivative has smooth dependence on spatial displacement. Typical displaced points are dense on this plaque, so the wedge vanishes smoothly along a patch. Its tested Taylor coefficient at zero is precisely the coefficient under consideration at the chosen actual state: finite-level re-centering preserves the needed jets, and the homogeneous lift covers their finite parameter germ. Thus this coefficient vanishes at every nearby requisite typical actual state. Lemma 18 now makes it vanish throughout a product neighborhood in \(V_j\). At a fixed spatial order all multi-indices can use the same finite level and connected fiber patch. The conclusion applies to each of their analytic coefficient functions there.

Choose \(d\) independent derivative columns at the central state, and let \(B\) be their formal continuation in the section-point variables. A corresponding \(d\)-minor is a unit. Substitute the formal solution of (26) in the wedge identities just proved. They show that every derivative column lies in the free rank-\(d\) module generated by \(B\). Differentiating a basis column in a hidden direction gives another derivative column; hence \[\partial_{t_j}B=B C_j\] for a formal coefficient matrix \(C_j\). Normalize the selected minor of \(B\) to the identity. Its hidden derivative is then zero, since the same relation has zero rows on that minor. This proves that the span is independent of hidden variables.

The chain rule under an invertible hidden reparametrization gives the same derivative span in both directions. Higher-level parameters are also hidden parameters, so changing a higher lift cannot alter the span. For each fixed ambient order choose a sufficiently high finite level; the unit-minor characterization makes the resulting jets consistent. Finally, one may retain the base variable \(x\) as an actual smooth variable while expanding only in \(y\) in (26). A smooth reference hidden section then gives smooth coefficients on the stable axis, whose base derivatives are exactly their higher ambient coefficients. ◻

Proof. Stable-axis agreement at nearby typical centers follows from Lemma 14 and the intrinsic constrained-span description. Smooth stable-axis coherence was proved in Lemma 24.

For the unstable axis take a homogeneous tower section over a level containing the base variables and the parameters furnishing the \(d\) independent columns. Use \(\Phi_v\) to carry the same constrained basis calculation to a small finite \(y\). The columns and every one of their formal point coefficients depend smoothly on \(y\). Their central minor remains invertible after shrinking the patch. At a typical arrival center, Lemma 23 places their full graph-state series in every finite arrival variety, including the prescribed zero-order density. For a prescribed point-jet order, take an arrival state level containing the required graph and density derivatives, including the extra graph derivative in the conormal form. Taylor re-expansion of the smooth spatial profile about this fixed arrival center is exactly the graph reconstructed from those state jets. Fix a full formal arrival point \(q=(x',y')\) and a finite jet order. Write \(S_{\mathrm{arr}}(q,t')\) for the arrival constrained augmented family from (26). In an arrival submersion chart, the state inclusion and the point constraint give, at the chosen finite order, \[S_{\mathrm{tr}}(q,t) =S_{\mathrm{arr}}\bigl(q,\Theta(q,t)\bigr),\] where \(t\) denotes the departure hidden parameters, \(\Theta(q,t)\) the induced arrival hidden parameters, and both augmented families are expressed in the arrival chart. Only finitely many state levels enter this identity. Differentiating in \(t\) while keeping \(q\) fixed shows by the chain rule that every transported column belongs to the arrival constrained span. That span also has rank \(d\), so the retained independent columns give the whole bundle. Inclusion and equality of ranks suffice; no surjectivity between infinite towers is required. Equality at each finite order proves full ambient coherence on this unstable-axis patch.

All operations commute with raw graph transformation, normal comparison, and the appropriate conormal height factors. The constant height multiplier at a fixed point is independent of the alternative being spanned. Thus the constructed jets are equivariant for dynamics and changes of section. Regular centers related by a fixed flow shift have compatible jets by Lemma 14.

To cover an entire typical stable leaf, fix a finite collection of orders. Local smooth representatives have a positive radius at almost every center. Choose forward returns on which their radius is bounded below and the required finite construction bounds are bounded above. Every compact portion of the central stable leaf contracts into these patches; pull the representatives back. Two pulled representatives agree on the dense set of typical points and hence everywhere. The same argument backward in time covers the entire unstable leaf. Repeat one finite collection of orders at a time. Agreement of the orders and their derivative identities follows on overlaps by density and smoothness. This does not require a single radius valid at all orders. Measurability of the radii used here is established below. ◻

Proof. Raw leaf charts and actual finite plaque jets have measurable, indeed continuous local transverse dependence. The adapted normal changes have measurable finite jets and domains. All subsequent finite constructions use polynomial spans, kernels, nonzero minors, implicit coordinate charts, or limits of measurable coefficient functions. Countable choices of candidate minors and coordinate balls give measurable versions. The coefficients of the resulting formal fields are therefore measurable.

It remains to justify the patch bounds without presupposing a smooth transverse choice of representative. In a fixed countable plaque box, the smooth representative on a typical plaque agrees almost everywhere with these measurable coefficients. Its norm on a compact subpatch is the essential supremum of the same coefficient norms there; this is measurable by conditional integration. Derivative norms are tested in the same way using the higher coherent coefficients. For a radius bound, work in a fixed Grassmann chart with a chosen minor bounded away from zero and impose bounded jets and Taylor difference estimates through an order strictly higher than the differentiability being requested, between almost every pair of points on a smaller plaque ball. These tests are measurable by Fubini and countably many radii and bounds. Existing smooth coherent representatives satisfy them locally. Conversely, such bounds extend the tested jets continuously to the closure of a dense valid set and give their stipulated finite differentiability and coherence. A dense countable set of mutually valid points may be selected by successive Fubini choices. These tests provide the measurable positive radii required in the globalization proof.

On weak plaques use the whole-strong-leaf representatives and smooth flow transport. At points for which both constructions are regular, their jets agree by the flow comparison already proved. Strong-leaf conditional agreement and flow Fubini then give agreement almost everywhere on the weak plaque; agreement at every flow time is unnecessary. Compact subpatches have finite smooth norms. Conditional essential suprema make these bounds measurable on the family of weak plaques as well.

An almost-everywhere finite measurable bound is bounded on a positive-measure subfamily of any positive-measure plaque-label set. The same applies to the reciprocal of a positive radius and to any finite list of bounds. Absolute continuity of the strong and weak foliations and a countable cover allow all conditional agreement requirements to be imposed simultaneously, also at integer iterates. These statements concern bounds along plaques; they do not assert that the formal fields themselves have tempered norms along every regular orbit. ◻

Proof of Theorem 13. The unstable affine field is supplied by Lemma 24. Its central value contains the actual augmented form and hence the actual form. Lemmas 25 and 26 give the full ambient coherence, globalization, and measurable finite-order bounds. Restore the height dependence \(e^{-h\tau}\).

Apply the same construction with time reversed to obtain \(\mathcal A^s\), restoring height dependence \(e^{h\tau}\). One may choose charts measurably in flip pairs, with the axes interchanged; all graph comparisons preserve the construction, so it is compatible with flip. Proposition 19 gives the constant cross pairing of the centered affine spans. The transformation laws follow from the exact raw form transformation and the fixed-point nature of the span. This proves every assertion of the theorem. ◻

Proof. The symplectic pairing between Lyapunov spaces of exponents \(\lambda,\mu\) vanishes unless \(\lambda+\mu=0\). To see this, transport the pairing along the orbit and use the two-sided Lyapunov estimates; when the sum is nonzero, transport in the direction in which the product of the two norms tends to zero. A measurable symplectic Oseledets frame therefore pairs opposite spaces and has tempered finite-dimensional change-of-frame costs. Centered Darboux coordinates with these first derivatives give the initial constant symplectic form. The field \(D\) preserves \(\omega_0\).

Suppose that the section map has already been made weight-preserving through degree \(j-1\), with \(j\ge2\). Factor out its linear part. The non-weight-preserving part of its degree-\(j\) term is a homogeneous symplectic vector field. More explicitly, the symplectic pullback identity in degree \(j-1\) has a linear term \(\mathcal L_V\omega_0\) and terms determined by lower degrees. Those lower terms have weight zero; its nonzero-weight part is therefore \(\mathcal L_{V_{\ne0}}\omega_0=0\). A homogeneous symplectic vector field is Hamiltonian, by the polynomial Poincare lemma.

On each nonzero-weight component solve the nonresonant homogeneous coordinate-change equation, and take zero correction on the zero-weight component. Explicitly, write the degree-\(j\) error after factoring the linear map as \(V_n\), and the coordinate correction as \(Q_n\). Conjugating by the correction gives the equation \[Q_{n+1}=A_{n*}Q_n-A_{n*}V_n, \qquad (A_{n*}Q)(q)=A_nQ(A_n^{-1}q).\] The operator \(A_{n*}\) preserves Hamiltonian vector fields and has the single exponent \(-\gamma\) on their \(D\)-weight-\(\gamma\) module. Thus Lemma 9 applies for \(\gamma\ne0\). At this degree the nonzero exponent differences form a finite set bounded away from zero. The appropriate forward or backward Green series consequently gives a tempered correction. The formal time-one map of its Hamiltonian vector field is symplectic and changes no lower degree. Induction gives compatible coordinates at all orders. This argument imposes no discreteness condition on the union of the weights over all degrees.

The radial primitive \(\theta_0\) satisfies \(\mathcal L_D\theta_0=0\). A symplectic coordinate change has a strict lift because the difference of its pulled-back primitive and the original primitive is formally exact. Normalize its primitive at zero. For the return map, strictness gives \(dr=\theta_0-F^*\theta_0\). Both terms are \(D\)-invariant, so \(d(Dr)=0\); evaluation at zero gives \(Dr=0\).

Finally, a graded section map preserves each signed linear axis: a monomial involving only negative-weight inputs cannot have positive output weight, and conversely. The actual stable germ is a graph over the negative axis with tempered finite jets. If the first nonzero graph coefficient had degree \(j\), its invariance would give a homogeneous recurrence whose exponent is a positive output weight minus a sum of negative input weights. This exponent is nonzero. A homogeneous recurrence of nonzero exponent has no nonzero two-sided tempered solution. Thus this coefficient vanishes; induction removes every graph coefficient. The unstable argument uses the opposite signs. The strong leaves are horizontal. On either signed linear axis \(\theta_0\) vanishes, so their centered horizontal lifts have height zero as formal germs. ◻

Proof. For clarity, consider the finite-dimensional space of augmented section jets through a fixed order \(N\). Its subspace consisting of jets valued in \(\mathcal B^u\) is measurable and invariant under the return cocycle. In the normalized charts that cocycle acts on the form component by \[ \beta\longmapsto e^{hT-hr}F^*\beta, \tag{34}\] and on the scalar component by composition with \(F\). For \(\mathcal B^s\) the multiplier in (34) is \(e^{-hT+hr}\). These are pullbacks from the next center to the present center; they may equivalently be viewed as a cocycle over the backward base map.

For an \(m\)-form monomial \(q^a\,dq_{i_1}\wedge\cdots\wedge dq_{i_m}\), its ordinary weight is \[\sum_i a_i\lambda_i+\lambda_{i_1}+\cdots+\lambda_{i_m}.\] Add \(hT\) for the \(u\) form component and subtract \(hT\) for the \(s\) form component. Scalar monomials have their unshifted ordinary weight. The return cocycle preserves these augmented weights: \(F\) is graded and \(r\) has weight zero. Each block has the single Lyapunov exponent indicated by its weight. Indeed its diagonal ordinary-degree blocks are tensor representations of the linear return maps with precisely that exponent. The remaining terms are triangular in ordinary degree; at a fixed truncation there are only finitely many such extensions. Their coefficients are tempered, and the finite iterated sums in an upper triangular product add only subexponential costs to the common exponent.

Lemma 8 now splits the measurable invariant jet subspace into these weight blocks. The ambient finite jet cocycle of the smooth return maps on the compact phase space has bounded one-step coefficients, and its formal change of coordinates is tempered at the fixed order. No tempered bound on the hull subspace itself is being assumed. The infinitesimal diagonal action is \((\mathcal L_D+hT, D)\) on the two augmented \(u\) slots and \((\mathcal L_D-hT,D)\) on the \(s\) slots. It therefore preserves both augmented subbundles through order \(N\). Since \(N\) was arbitrary, Equations (28)–(29) show that \(D-T\partial_\tau\) preserves both formal affine fields. The identity \(d(\theta_0(D))=-\iota_D\omega_0\) proves (33).

For any Hamiltonian \(p\), \(\{H,p\}=Dp\). Thus the asserted diagonalization is literal monomial diagonalization. To project onto a prescribed weight \(\gamma\), choose on each finite Taylor quotient the polynomial spectral projector onto \(\gamma\) for \(D\). Applied to a fixed formal series these projectors eventually agree at every fixed coefficient, and converge coefficientwise to its weight-\(\gamma\) part. A closed invariant subspace contains this limit. A continuous finite-dimensional quotient is detected at a finite Taylor order, so the same projectors commute with passage to that quotient and give the claimed lifts. The argument applies also to closed invariant subspaces before taking such a quotient.

All coordinate weights in degree one are nonzero. A weight-zero Hamiltonian therefore has zero differential at zero and belongs to \(\mathfrak k\). Evaluation at zero sees only the constant monomial, which has weight zero. ◻

Proof. Since \([D,Y_p]=0\) and \(Y_p(0)=0\), the derivative \(A_p\) commutes with \(D\). It preserves the signed planes, which are the actual strong planes by Lemma 27. At the center, the actual form \(u\) annihilates \(E^u\) and is a volume form on \(E^s\). Its infinitesimal transform under \(Z_p\) is therefore \[\bigl(\mathop{\mathrm{tr}}(A_p|_{E^s})-h\xi(p)\bigr)u.\] Preservation of the affine field says that this vector belongs to the direction space of \(\mathcal A^u(0)\). Pairing any such direction with the actual opposite form \(w(0)\) gives zero, by the constant cross-pairing property in Theorem 13. Since \(u(0)\wedge w(0)\ne0\), the displayed scalar is zero. Complexification causes no change in this argument. ◻

Proof. Use the background extensions of Lemmas 25 and 26. Write \(K_N(a)\) for the finite space of \((N+2)\)-jets satisfying Equation (30) through degree \(N\) at center \(a\). Isotropy, when required, is an additional linear condition. For \(N+2\ge r\), set \[V_{r,N}(a)=\pi_rK_N(a).\] Both \(K_N\) and the map defining this image are finite-dimensional linear algebra in finitely many background jets. Their matrices are smooth on the axis.

The ranks of the defining matrix \(A_N\) and the augmented matrix obtained by stacking \(A_N\) with \(\pi_r\) are invariant under the return dynamics, because finite symmetry conditions and jet projection are natural under centered coordinate changes. By ergodicity these ranks are constant almost everywhere, simultaneously for all \(r,N\) under consideration. These two ranks determine \(\dim K_N\) and \(\dim V_{r,N}\). At a regular center a nonzero maximal minor remains nonzero on a small axis patch. Dense regular agreement prevents a larger rank anywhere on that patch: otherwise an open set of larger rank would meet the dense regular set. Thus both ranks, and consequently the kernel and image ranks, are constant near the center. Constant-rank linear algebra gives smooth kernel and image bundles there.

This local assertion extends to the entire typical axis. Fix the finite list of matrices just used. Their local patches have measurable positive radii after also recording finite coefficient bounds and a lower bound for the chosen nonzero minors. Such data are provided by Lemma 26. Some positive-measure set has these radii bounded below. Returns of a typical center to that set, in the direction contracting the axis, bring every fixed compact portion of the axis into one of these patches. Pullback preserves the finite equations and both ranks. This proves the constant-rank conclusion on that compact portion. Exhausting the axis and intersecting over the countable list of finite tests gives the assertion for the whole axis. Forward time is used for the stable axis and backward time for the unstable axis.

For fixed \(r\), the spaces \(V_{r,N}(a)\) decrease with \(N\) inside a finite-dimensional vector space. At regular centers all their dimensions are the same deterministic sequence of integers. Choose \(N(r)\) after its eventual stabilization; a first adjacent equality would not suffice. The preceding whole-axis rank argument applies to every \(N\), so for every point of that axis \[ V_r(a):=\bigcap_N V_{r,N}(a)=V_{r,N(r)}(a). \tag{36}\] It remains to identify this intersection with full solutions.

We use the elementary linear compactness principle: a system of linear equations in coefficient variables, each equation involving only finitely many variables, has a solution if every finite subsystem has a solution. One proof regards each equation as a finitely supported row vector. Finite consistency makes the specified right sides a well-defined linear functional on the span of the rows; extend it to the direct sum of all coefficient coordinates. Its values on the coordinate vectors give a solution. If \(v\in V_r(a)\), impose the prescribed \(r\)-jet \(v\) together with all formal symmetry equations. Every finite subsystem is solvable by the defining intersection, so this principle gives a full solution with that jet. Conversely every full solution projects into each \(V_{r,N}\). This proves the claimed equality.

Equation (36) makes \(V_r\) a smooth bundle. A full solution lifting any of its elements also supplies a lift at order \(r+1\), so \(V_{r+1}\to V_r\) is surjective in every fiber. It is a smooth map of constant-rank bundles. All arguments apply unchanged to the isotropy condition or a finite list of linear conditions with the stated smoothness and equivariance. ◻

Proof. Fix a regular center \(a(0)=0\), and a smooth path \(a(t)\) in one of its actual axes. Fix a desired equation order \(r\). By Lemma 30, extend the \((r+3)\)-jet of its particular isotropy \(H\) to a smooth section \(p_t\) of the full isotropy projection bundle on that axis. Each \(p_t\) is represented by a Taylor polynomial in \(z=q-a(t)\). It satisfies the symmetry equations through degree \(r+1\) and \(Y_{p_t}(a(t))=0\).

Differentiate while holding the original coordinate \(q\) fixed. For a moving-center series \(C_t(z)\) this differentiation is \[\partial_t^{\mathrm{fixed}}C_t =\partial_t C_t-\sum_i a_i'(t)\partial_{z_i}C_t.\] Full ambient holonomicity of the background states that its fixed-coordinate derivative vanishes, coefficient by coefficient. The same is true for smooth frames and quotient matrices obtained from it by the Grassmann-chart operations. Differentiating the linear equations through degree \(r+1\) therefore shows that \(\partial_t^{\mathrm{fixed}}p_t\) satisfies them through degree \(r\). The extra degree is what compensates for differentiation of the moving Taylor center; the order-two Hamiltonian operator accounts for the choice of the \((r+3)\)-jet.

Differentiating isotropy by the ordinary chain rule gives \[ Y_{\partial_t^{\mathrm{fixed}}p_t}(0) =-dY_H(0)a'(0). \tag{37}\] At the chosen center, \(dY_H(0)\) is conjugate to the signed dilation \(D\). It is invertible and preserves each actual signed tangent plane. Using paths in the two axes therefore gives, at every prescribed finite equation order, a solution with any prescribed tangent value. Finite sums handle a value with both signs. Apply the linear compactness principle from Lemma 30, with this tangent value prescribed, to obtain a full formal solution. No compatibility of the separately chosen finite families is needed.

Finally, the dimension of the evaluation image is the difference between the dimensions of the full symmetry and isotropy projections at first jet order. Those dimensions have the same constant values on the entire typical axes by Lemma 30. Surjectivity at their regular points consequently implies surjectivity throughout their smooth extensions. ◻

Proof. Transitivity and the Hamiltonian bracket formula give \[\xi(\{p,q\})=\omega_0(Y_p(0),Y_q(0)).\] Its radical is exactly the kernel of tangent evaluation. This is both the asserted symplectic identification and the assertion \(\mathfrak k=\operatorname{stab}_{\mathfrak h}\xi\).

Suppose an ideal \(I\) lies in \(\mathfrak k\) and \(p\in I\). Every iterated bracket of its projected field with fields from \(\mathfrak h\) vanishes at zero. If \(Y_p\) had a first nonzero homogeneous term of degree \(j\ge1\), brackets with symmetries having prescribed constant tangent values would differentiate that term. After \(j\) such brackets a nonzero constant value would result. At each step the other bracket term has higher order, so it cannot cancel this leading derivative. This contradicts \(I\subset \mathfrak k\). Hence \(Y_p=0\) formally and \(p\) is constant. A constant Hamiltonian has weight zero; Lemma 29 then gives \(hp=0\). Thus \(I=0\).

For \(p\in\mathfrak k\) and any \(q\in\mathfrak h\), evaluation of the projected bracket gives \(Y_{\{p,q\}}(0)=-dY_p(0)Y_q(0)\). Apply this first to \(p=H\) to obtain the action \(-D\), and then use Lemma 29 on the stable plane to obtain Equation (38). ◻

Proof. Choose symmetries \(p_1,\ldots,p_{2m}\) whose projected values form a basis. Write \(G_t\) for the composition of their formal flows with parameters \(t_1,\ldots,t_{2m}\), and \(F(t)=G_t(0)\). Its linear term is invertible, so the formal inverse function theorem gives \(t=t(q)\). Transport \(V\) by \(dG_t(0)\) and substitute \(t(q)\). All operations are well-defined in formal displacement parameters of positive order.

We verify independence and invariance. Two such parameterized translations reaching the same point differ by a formal family fixing zero. Its right logarithmic derivative belongs coefficientwise to the projected symmetry algebra: composition and inversion give sums of conjugated generators, and their formal adjoint series are iterated brackets. Since the family fixes zero, these derivatives belong to the projected isotropy. Its derivative at zero consequently preserves \(V\), by the linear differential equation starting at the identity. This proves independence. Including the formal flow of any further symmetry in the comparison proves invariance. Conversely any invariant field must be transported from its value by these spanning flows, which proves uniqueness. The same argument applies to invariant tensors and to the asserted compatibility operations.

For involutivity, the Levi bracket of an invariant distribution at zero, evaluated on the values of \(p,q\in\mathfrak l\), is the negative of the projected bracket \(\{p,q\}\) modulo \(V\). One can check this without assuming that the fundamental fields themselves are tangent everywhere: choose tangent sections with the same initial values and use that their brackets with symmetry fields remain tangent. Subtracting the two identities gives the claimed Levi bracket. It vanishes since \(\mathfrak l\) is a subalgebra. Invariance and the spanning formal translations make the obstruction vanish as a full formal tensor. ◻

Proof. Work first on a smooth axis patch in a fixed section, and let \(a(t)\) be a path there. By Proposition 31 and Lemma 30, the full finite symmetry projections surject smoothly onto prescribed tangent values. Choose a smooth section with value \(a'(t)\) at the first required order. Lift it successively through the surjective smooth bundle maps of Lemma 30, using smooth right inverses chosen from bundle metrics on the axis. This constructs a compatible family of formal symmetries \(p_t\) whose every coefficient is smooth and for which \[Y_{p_t}(a(t))=a'(t).\] The choice can be made smoothly in auxiliary path parameters by the same finite-dimensional right inverses. There is no assertion of a bound uniform in the formal order.

Integrate these symmetries in coordinates relative to the moving center. In the section the generating field is \[Y_{p_t}(a(t)+z)-a'(t),\] which vanishes at \(z=0\). Its ordinary differential equation therefore closes at every finite Taylor order: the order-\(N\) transported jet depends only on the order-\(N\) generator jet and lower orders. Solving these compatible finite-dimensional equations produces a smooth formal coordinate change. For the strict contact lift, choose the center height \(\beta(t)\) by \[\beta'(t)=p_t(a(t))-\theta_{a(t)}(a'(t)).\] The lifted generator relative to \((a(t),\beta(t))\) then also vanishes at its formal origin. The resulting formal changes are strict.

They intertwine the augmented affine fields. Indeed differentiate the pullback of their frame jets, modulo the corresponding bundle. Moving-center differentiation is the formal derivative in the center velocity, by whole-axis holonomicity and the prescribed smooth height dependence. The relative generator subtracts that same velocity. The remaining term is its infinitesimal symmetry equation (30), and is zero in the bundle quotient. An internal change of frame is permitted. This verification is made at each finite order, where it is ordinary differentiation of a finite system. At the endpoint, the center height \(\beta\) introduces the reciprocal constants \(e^{-h\beta}\) and \(e^{h\beta}\) in the two affine fields. They do not change \(\mathfrak h\). Strictness and transport of the center preserve both isotropy and the evaluation functional \(\xi\).

An intrinsic construction at the centers is thus constant under these algebra identifications. Proposition 33 commutes with the identifications, so its entire formal field is transported and has smooth coefficients. To verify holonomicity, let \(T_t\) denote its moving-center formal tensor, projected to the section. For a distribution the following identities are interpreted in a Grassmann chart. The transport identity gives \[\partial_t T_t+\mathcal L_{Y_{p_t}-a'(t)}T_t=0.\] Homogeneous invariance gives \(\mathcal L_{Y_{p_t}}T_t=0\); hence \[\partial_tT_t=\mathcal L_{a'(t)}T_t.\] For the constant displacement vector \(a'(t)\) this is exactly the coefficientwise Taylor coherence identity. The argument for lifted tensors includes the center height velocity in the same way.

For a finite intrinsic list, use the smooth family of paths to trivialize the structures over a whole axis patch. A continuous selection into that finite list is locally constant in this trivialization, and the preceding argument applies. This patch formulation also proves the assertion along paths that need not themselves consist of regular points. The whole-axis extensions of Lemma 30 permit these patches everywhere on the axis; uniqueness and regular density make their jets agree on overlaps. Naturality proves equivariance under the return map and under regular flow shifts. The prescribed height dependence, strong-axis coherence, and flow Fubini then give the weak-axis extension and its almost-everywhere agreement. ◻

Proof. Let \(A\) be a closed abelian ideal. For \(p\in A\cap\mathfrak h_0\) we have \(p\in\mathfrak k\) and \((\mathop{\mathrm{ad}}p)^2=0\) on \(\mathfrak h\). Consequently its induced trace on \(E_+\) vanishes, and (39) gives \(\xi(p)=0\). The functional \(\xi\) kills every nonzero weight because \(H\) stabilizes it. Weight projection therefore gives \(\xi(A)=0\). Since \([\mathfrak h,A]\subset A\), every element of \(A\) stabilizes \(\xi\); hence \(A\subset\mathfrak k\), which forces \(A=0\).

Consider a decreasing chain of nonzero closed ideals contained in a given nonzero closed ideal. Their images in \(E\) are nonzero: an ideal with zero image would be contained in \(\mathfrak k\). These images are nested finite-dimensional spaces. Choose a member whose image has minimal dimension, and a nonzero vector in that image. On the cofinal subchain below this member, the fibers over this vector form a nested family of nonempty closed affine subspaces of a linearly compact space. Their intersection is nonempty. The intersection of the chain is consequently a nonzero closed ideal. The usual minimal-element argument now gives a minimal nonzero closed ideal.

For a minimal ideal \(I\), the closure of \([I,I]\) is a closed ideal of \(\mathfrak h\). It is nonzero by the first paragraph, so it equals \(I\). Also \(\xi(I)\ne0\): otherwise the ideal property would put \(I\) in \(\mathfrak k\). By continuity and topological perfection, there exist \(a,b\in I\) such that \(\xi([a,b])\ne0\). Their values in \(E\) span a symplectic two-plane.

Distinct minimal ideals commute. Indeed their bracket is contained in their intersection, and a nonzero intersection would make them equal. The two-planes just constructed for distinct ideals are thus mutually symplectically orthogonal. There are at most \(\dim(E)/2\) of them. Finally, the common centralizer is a closed ideal. If nonzero, it contains a minimal ideal; this ideal would commute with itself, contradicting the first paragraph. ◻

Proof. The continuous centroid of a simple factor \(S\) is scalar. One may use the linearly compact Schur lemma, as in Fattori and Kac (2002, Lemma 2.2). In the listed realizations it also follows by commuting a centroid map with the grading in the \(W\) and contact cases, or with the degree-zero \(\mathfrak{sl}\) or \(\mathfrak{sp}\) action in the other cases. A grading eigenspace in the former cases, and a finite-multiplicity isotypic component in the latter, is a nonzero finite-dimensional centroid-invariant space. Such components exist because the homogeneous divergence-free polynomial fields and homogeneous Hamiltonians form irreducible modules with distinct highest weights as their degrees vary. An eigenvalue on this finite-dimensional space gives a nonzero kernel for the centroid map minus that scalar. This kernel is an ideal, so simplicity makes the map scalar.

A continuous centroid map of \(I\), restricted to constant \(S\)-series and tested coefficientwise against brackets with constants, is multiplication by a scalar series \(a(z)\). To see that this determines the map everywhere, subtract this multiplication map. For each monomial \(z^\nu\) its value on \(z^\nu S\) is zero: write elements of \(S\) as limits of finite sums of brackets and use \(z^\nu[a,b]=[z^\nu a,b]\), placing the centroid map on the constant entry. Polynomial series are dense, so continuity finishes the argument. Thus the centroid is exactly \(\mathcal O_e\).

A derivation acts on this centroid by commutator and hence induces a continuous derivation of \(\mathcal O_e\), namely a formal parameter vector field. For continuity, evaluate the commutator on a fixed nonzero constant \(s\in S\) and apply a continuous linear functional on \(S\) taking value one on \(s\). Continuity on \(I\) then gives continuity of the induced scalar-series operator. After subtracting the parameter field, the remaining derivation is \(\mathcal O_e\)-linear. Restriction to constant \(S\)-series, coefficient by coefficient, gives derivations of \(S\). Density again gives the decomposition and its uniqueness. This proves (42); it is also the continuous-derivation formula of Bakalov et al. (2001, Proposition 6.12(ii)).

The kernel of the action on all minimal ideals is their common centralizer, which is zero by Lemma 36. The displayed derivation spaces are linearly compact in their formal coefficient topology: factor normalizers are closed spaces of formal vector fields, finite simple factors give finite matrix spaces, and completed coefficient products preserve linear compactness. Continuity in the asserted jet topologies follows from continuity of bracketing. More explicitly, the normalizer realization embeds continuously into the product of its bracket evaluations on factor elements. It is injective and has continuous inverse on its image: a continuous injective map from a linearly compact space to a Hausdorff linearly topologized space has closed image and induces its original topology. Apply this fact first to the normalizer and then to the faithful map from \(\mathfrak h\). ◻

Proof. We first explain measurability without choosing a measurable classification of minimal ideals. Replace the countably many completed-measurable background coefficients by Borel versions on a common conull Borel set, refined under the countable group generated by \(f\) and flip. For this purpose, coefficient sequences are coded in the countable product of the usual standard Borel space \(\mathbb C\). This measurable structure is distinct from the formal linear topology with discrete scalars used in the algebraic arguments. The symmetry algebra is given by countably many linear jet equations. For the intrinsic condition (40), fix an output jet order. It asks for an input jet order such that every power of the adjoint annihilates the corresponding input kernel after output projection. For each power this is a finite-jet linear inclusion test; the full formal solution images are obtainable by the stabilized finite linear tests from Lemma 30. The quantifiers over orders and powers are countable. Thus membership in \(\mathfrak l\) is a Borel relation between the base point and coefficient sequences. Its finite-dimensional projections are analytic, hence universally measurable (Tserunyan 2025, Proposition 12.2 and Corollary 14.10). To recover the tangent plane measurably, for each \(d\) code \(d\) independent projected vectors. Their analytic relation has domain \(A_d\), the set where the dimension is at least \(d\), and admits a universally measurable tuple selection (Moschovakis 2009, sec. 4E.9). Restrict this already chosen tuple to \(A_d\setminus A_{d+1}\) and take its span. This gives the Grassmannian-valued plane for the completed Liouville measure; the exact-dimension strata need not be analytic. Subsequent countable coefficient choices can again be replaced by Borel versions before further coding. No Borel selector, or common Borel refinement invariant under every real flow time, is required.

Once the finite-dimensional plane \(L\) is known, brackets and finite jet lifts determine the representation on \(B\). The socle of its dual is specified by finite-dimensional invariant-subspace and rank conditions. Indeed invariance is a closed Grassmannian incidence condition for a finite basis of the represented matrix algebra; reducibility is a finite union of projections along compact nested Grassmannians. Thus irreducibility is Borel, and the preceding independent-tuple procedure selects the span of the irreducible submodules. These operations give measurable planes and, where needed, measurable choices of finite jets. The homogeneous formal extensions are computed to each requested order from finitely many such jet operations. All resulting constructions are natural for the dynamics. Their dimensions are therefore almost everywhere constant by ergodicity.

Conjugation preserves the intrinsic definitions and the real \(H\)-sign spaces; hence the complex planes descend to real ones. Flip intertwines the two affine-form data, preserves the underlying symmetry structure up to the harmless sign of the contact form, and consequently preserves \(P'\). It exchanges the actual stable and unstable parts. Since it preserves smooth volume and these two ranks are almost everywhere constant, the ranks are equal. Their sum is the positive rank of \(P'\), so both are nonzero.

For smoothness and coherence, use the formal identifications along whole axes supplied by Proposition 34. They carry \(\mathfrak h,\mathfrak k,\xi\), hence all the intrinsic objects just constructed, smoothly along the axis. The transported plane \(P'\) is preserved by the projections to the two actual sign bundles at typical points. These projections are continuous and uniformly transverse. By density they preserve the smooth extension of \(P'\) everywhere on the axis patch. Their restrictions are continuous idempotents, so their ranks are locally constant; their images give continuous sign subbundles.

In this trivialization, Proposition 41 places each sign subspace, at typical points, in a fixed finite list of full isotypic sums. The finite list is closed, so the same holds at all points by density. A continuous choice in that list is locally constant. The same reasoning applies to the finite-case sign spaces in \(E\). Proposition 34, including its version for local smooth families of paths, now gives smooth full formal jets and their holonomicity along the whole axes. The homogeneous extensions are invariant under \(\mathfrak h\) by Proposition 33, and thus under \(H\). Their center values are the prescribed actual sign planes. Finally, flow shifts intertwine all these intrinsic objects; the section fields and their horizontal lifts can be taken time-independent within a flow box. This supplies the weak-axis extensions with the required agreement at regular translates. ◻

Proof. The algebraic assertions are Propositions 35 and 41; the real, measurable and coherent conclusions are Proposition 45. The subalgebra \(\mathfrak l\) gives an involutive homogeneous formal distribution by Proposition 33. Coisotropy gives its characteristic distribution. Invariance and ergodicity make the alternative almost everywhere constant. ◻

Proof. On a transversal section, formal Frobenius supplies independent first integrals \(z_1,\ldots,z_d\) of \(L\), where \(d=\dim P\). Since \(L\) is coisotropic, the Hamiltonian fields \(X_{z_i}\) lie in \(L\) and span \(P\). Furthermore, \[\{z_i,z_j\}=dz_j(X_{z_i})=0, \qquad [X_{z_i},X_{z_j}]=0.\] Declare this frame parallel. A different system of first integrals has the form \(\widehat z=h(z)\), with invertible differential, and \[X_{\widehat z_i} =\sum_j \frac{\partial h_i}{\partial z_j}(z)X_{z_j}.\] The change-of-frame matrix is constant on every characteristic leaf. Thus the connection is independent of the first integrals; its torsion and curvature vanish in the displayed commuting frame. Since \(P\) is isotropic, the bracket of two horizontal lifts has no additional Reeb component: that component is the negative of their \(d\alpha\) pairing. Hence the horizontal lifts commute as well and give the asserted lifted affine structure.

We verify the parallelism of the sign fields. A symmetry preserves the ring of first integrals of \(L\), so it induces a formal field on their base. At the origin, the base tangent space is \(B\). For \(p\in\mathfrak l\), the induced base field vanishes there, and the negative of its linearization is the bracket representation on \(B\). The Hamiltonian translation frame identifies \(P\) with the cotangent space of this base. Its linear transformation under that symmetry is consequently the transferred dual representation on \(B^*\).

Every projected symmetry preserves both \(\omega\) and \(L\), and hence preserves \(P=L^\omega\). If \(p\in\mathfrak u\), its projected value belongs to \(P\). In formal Frobenius coordinates for \(P\), its flow fixes the characteristic leaf through the origin: its transverse component depends only on transverse coordinates and vanishes at their origin. On this leaf its action on the parallel frame is the fixed cotangent derivative at the base origin. This action preserves both \(P^s_a\) and \(P^u_a\), since these are full \(\mathfrak u\)-invariant submodules. The values of \(\mathfrak u/\mathfrak k\) span \(P_a\). Its formal flows are thus transitive along the characteristic leaf, and the homogeneous sign fields are constant in the parallel frame. A strict lift may also translate height; the horizontal structure and its parallel frame are independent of that constant height translation. This proves parallelism for the horizontal lifts.

The connection and its exponential are natural under the maps preserving the data, hence under \(f\). They are computed to any given order by finite formal operations: solve the formal Frobenius equations and then the affine geodesic and horizontal-lift equations in an additional time parameter. Coherence and measurable finite bounds are preserved by these operations on a fixed nonsingular coordinate patch.

Finally use the signed coordinates of Lemma 27 and the dilation of Proposition 28. In these coordinates the actual central unstable and stable plaque Taylor series are the positive and negative axes. The formal affine exponential is dilation-equivariant. A Taylor monomial in only positive-weight initial coordinates cannot have a negative-weight output component, and conversely for negative initial coordinates. Its pure-sign restriction therefore has no component toward the opposite axis. Horizontality fixes its height to that of the actual strong plaque. This is an identity of Taylor series and makes no assertion about their convergence. ◻

Proof. The extended values of \(P^s\) belong to the actual stable tangent space everywhere on the axis, by agreement on a dense set and continuity. Derivatives in \(P^s\) directions are actual derivatives within this smooth axis. Holonomicity identifies them with the corresponding formal derivatives. Thus the formal Frobenius identity becomes the actual involutivity identity there. The same reasoning applies to torsion, curvature, and preservation of the two sign bundles: all their derivatives, when restricted as asserted, are stable-tangential derivatives. Frobenius now gives actual smooth integral slices, and the restricted connection is flat and torsion-free. On small simply connected slices, parallel transport defines parallel frames. In smooth slice coordinates, its ordinary differential equations have coefficients depending smoothly on the transverse label, which gives the stated smooth choice of frames. The identities extend over nonregular points of the axis through the smooth extensions. Equivariance likewise extends by dense agreement. ◻

Proof. In bounded coordinates \(z\) on a stable plaque, the affine geodesic equation on \(P^s\) has the form \[\dot z=v,\qquad \dot v=A(z)(v,v),\qquad v\in P^s_z.\] The acceleration coefficients are finite-order operations on the formal connection coefficients and become smooth actual coefficients on the stable plaque by coherence. The geodesic spray is tangent to this velocity subbundle, since the connection preserves \(P^s\). One can extend the coefficients smoothly to all stable velocities by a smooth projection to \(P^s\) when applying the ordinary existence theorem.

Their suprema on fixed compactly contained stable plaques are measurable and finite at almost every center. Restrict to a positive-measure set on which they have a common bound \(C_A\). The geometric plaque charts and their metric comparisons already have uniform bounds. If \(|v(0)|\) is sufficiently small, the inequality \(|\dot v|\le C_A|v|^2\) gives \[|v(t)|\le 2|v(0)|,\qquad 0\le t\le1.\] Thus the solution remains in the interior of the chosen plaque and its velocity remains in a compact coordinate neighborhood; ordinary continuation gives existence through time \(1\). Integration gives \[z(1)-z(0)=v(0)+O(|v(0)|^2)\] with a uniform constant. Shrink the input radius so that the coordinate displacement is bounded below by \(|v(0)|/2\). Ambient distance in the bounded geometric chart provides the corresponding lower bound for intrinsic stable distance. The length estimate gives the upper bound. Uniform norm comparisons give the constants in (49). ◻

Proof. If the frequency is \(p>0\), then \(i/n_i\to p\). Consequently \(n_{i-1}/n_i\to1\), which is the assertion. The last statement is the ergodic theorem for the indicator of \(K\). ◻

Proof. Choose \(a\) whose forward visits to \(K\) have positive frequency. For any initial vector \(x\in P^s_a\), contraction in (48) makes \(|df^n x|\le\eta_0\) at all sufficiently large good visits. Solve the local affine geodesic there through time \(1\) and pull it back by \(f^{-n}\). Equivariance of the actual restricted connection makes the resulting path an affine geodesic with initial velocity \(x\) at \(a\). Uniqueness of the ordinary geodesic equation makes the construction independent of the chosen sufficiently large visit. Applying it to \(tx\) for arbitrary finite \(t\) proves completeness. On a neighborhood of any fixed input, one sufficiently large good visit works for all inputs in that neighborhood, so the resulting exponential is smooth.

For properness, let \(|x|\to\infty\) and let \(n=n_i\) be the first good visit with \(|df^n x|\le\eta_0\). Then \(n\to\infty\): the inverse norms of the finitely many earlier linear maps are bounded. For large \(x\) there is a previous good visit \(n_{i-1}\), and minimality gives \(|df^{n_{i-1}}x|>\eta_0\). The global one-step inverse derivative bound and Lemma 50 give \[ |df^n x| \ge C^{-1}\eta_0 e^{-K(n_i-n_{i-1})} =e^{-o(n)}. \tag{50}\] The constants in this inequality do not depend on the varying initial vector. By equivariance and (49), \[d_s\bigl(f^na,f^n\operatorname{Exp}^s_a x\bigr) \ge c_1 |df^n x|.\] If the endpoints \(\operatorname{Exp}^s_a x\) remained in a fixed compact subset \(C_0\) of \(W^s(a)\) along such a sequence, their intrinsic distances from \(a\) would be bounded. Stable path-length contraction would instead give \[d_s\bigl(f^na,f^n\operatorname{Exp}^s_a x\bigr) \le C_{C_0}e^{-b_sn},\] contradicting (50). This proves properness. For \(x_0\ne0\), the path \(t\mapsto\operatorname{Exp}^s_a(tx_0)\) is the complete affine geodesic with initial velocity \(x_0\) and is proper as \(t\to+\infty\). Neither this argument nor the preceding construction assumes injectivity of the affine exponential. ◻

Proof. The strong and weak stable plaque estimates follow from Lemma 3.

For the ambient estimate, truncated substitution on polynomials of degree at most \(j\) turns the jet of a composition into a product of fixed-size matrices with uniformly bounded one-step norms. Their products grow at most exponentially. Apply this to \(f\) and \(f^{-1}\). Lastly, contraction of stable path lengths gives \[d_s(f^r a,f^r a')\le C e^{-b_sr}d_s(a,a').\] A ball of radius smaller than \(C^{-1}\) times the intrinsic distance from \(\overline D_0\) to the complement of \(D_1\) in its leaf about a point of \(f^{-r}D_0\) therefore stays in \(f^{-r}D_1\). ◻

Proof. Leibniz differentiation affects at most \(j\) factors. There are at most \(C_j(r+1)^j\) choices of factors and distributions of the derivatives among them. Each term has at most \(j\) differentiated factors, whose norms contribute a constant depending on \(j\), and at most \(r\) undifferentiated factors. Finally \((r+1)^j\le C_j e^r\). Include the terminal vector as one additional factor. ◻

Proof. Taylor interpolation on a fixed unisolvent polynomial stencil, rescaled to size \(t\), gives \[|D^j a(x)|\le C_{j,r} \left(t^{-j}\|a\|_{C^0} +t^{r-j}\|a\|_{C^r}\right), \qquad 0<t<t_0,\] where \(t_0\) is determined by the interior room. Choose \(t=(\|a\|_{C^0}/\|a\|_{C^r})^{1/r}\) when this is at most \(t_0\), and otherwise use \(t_0\). This proves (51); the exponential consequence is immediate. Uniform coordinate changes preserve the estimate. ◻

Proof. Choose finitely many product boxes for the strong stable and weak unstable foliations, with inner boxes covering \(SM\). Choose the outer boxes so that their stable plaque diameters are smaller than \(\rho/2\). All boxes can be taken compactly inside slightly larger ones. The product measure class in each box is the product of its transversal Lebesgue classes: this follows by stable disintegration and absolute continuity of weak unstable holonomy between stable transversals.

For almost every weak unstable plaque, the required tensor jets have smooth coherent representatives on its full compact crossing plaque, agreeing with the global regular data almost everywhere. Their suprema there are finite and measurable by Lemma 26. In each box choose a positive-measure family of such weak unstable plaques with a common finite bound. Product measure equivalence and Fubini show that almost every stable plaque intersects this family, in positive stable measure, at regular agreeing points. There are only finitely many boxes, so one bound \(B_k\) works for all of them.

The properties that a local stable plaque has such intersections and that global regular agreement holds along it discard an ambient null set. Refine typical centers so that their whole stable leaves meet each of these null sets in intrinsic measure zero. This is obtained by local absolute continuity at all integer iterates and the exhaustion \[W^s(a)=\bigcup_{r\ge0} f^{-r}W^s_{\mathrm{loc}}(f^r a).\] Now take any radius-\(\rho\) disk, including one centered at a nonregular point of such a leaf. Its inner half contains a point with all the preceding good-plaque properties. An inner product box containing that point has its entire outer stable crossing plaque inside the disk, by the diameter choice. That plaque contains the required bounded-jet regular point. Countable intersection treats all \(k\) and all integer translates. The same argument in lifted boxes proves the assertion on \(S\widetilde M\). ◻

Proof. At the orbit points \(s,f s,\ldots,f^ns\), center the one-step inverse maps by geometric exponential and logarithm charts. Their finite jets act by substitution on polynomials of degree at most \(k\). This realizes the finite jet of their composition as a matrix product of fixed size. The undifferentiated one-step matrix norms are uniformly bounded by a constant depending only on \(k\). Locally, fix a small stable ball about the point being estimated and hold fixed geometric frames on its successive forward images. Contraction puts these images inside uniformly bounded ambient charts; use the same intermediate frame in the two adjacent factors. These local choices require no global bundle trivialization.

Differentiate these matrices in \(s\) along the moving stable plaque. Each coefficient is a uniformly smooth geometric one-step coefficient evaluated at one of the forward points \(f^i s\). Lemma 52 gives bounds for every fixed derivative order of these forward plaque maps, uniformly in \(i,n\). Thus the differentiated matrix factors satisfy the hypotheses of Lemma 53. The terminal order-\(k\) jets of \(\mathcal F\) at \(f^ns\in D^+\) are smooth and have bounded derivatives on an outer compact enlargement; their composition with \(f^n\) has the same property. That lemma proves (53). This argument differentiates bounded geometric jet matrices; it does not differentiate measurable normal coordinates. ◻

Proof. Fix \(k\). Smooth finite-jet equivariance supplies a number \(D_k\ge0\) with the following property: transporting the jets of \(\mathcal E\) through order \(k\) forward by \(r\) steps costs at most \(C_k e^{D_kr}\) when the starting jets are bounded. Indeed, the output finite jets of \(f^r\) and the inverse linear tangent map on the input have at most exponential growth, by Lemma 52; only finitely many products occur at the fixed order \(k\).

Choose \[0<\epsilon\le\min\left\{1, \frac{b_uk}{8(D_k+1)}\right\}, \qquad \delta=\frac{b_s\epsilon}{4}, \qquad r_n=\lfloor\epsilon n\rfloor.\] Consider an intrinsic stable disk of radius \(e^{-\delta n}\) inside the moving outer patch. Its inverse image under \(f^{r_n}\) contains an intrinsic disk of radius \[C^{-1}e^{b_sr_n-\delta n}.\] This follows by applying the stable path-length contraction to any point in the latter disk. The radius tends to infinity. Hence Lemma 55 supplies a regular point in that inverse image where the required jets through order \(k\) are bounded. Pushing this point forward \(r_n\) times gives a point in the original disk where those jets are bounded by \(C_k e^{D_kr_n}\). Consequently such points form an \(O(e^{-\delta n})\)-dense set on every inner moving patch.

At each of these points, \(\|df^{-n}Y(f^ns)\|\le C_Ye^{-b_un}\) because \(Y\) takes values in the actual unstable bundle. Formula (52), homogeneous of degree \(k\) in this input, therefore gives \[ |G_{n,k}(s)| \le C_k e^{D_k\epsilon n-b_ukn} \le C_k e^{-7b_ukn/8}. \tag{55}\] The bounded-point test used jets only through order \(k\); the interpolation orders chosen next do not enter that test.

To estimate values between the sample points, take a fixed unisolvent stencil for polynomials of degree less than \(p\), scaled to radius \(t_n=e^{-\delta n/2}\) about the point being estimated. Perturb each stencil point to a point satisfying (55). The relative perturbation is \(O(e^{-\delta n/2})\), so the interpolation matrices and their inverses remain uniformly bounded. Taylor’s formula and Lemma 56 yield \[\|G_{n,k}\|_{C^0} \le C_{p,k}\left(e^{-7b_ukn/8} +e^{(A_k-\delta p/2)n}\right).\] Choose \(p\) with \(\delta p/2\ge A_k+3b_uk/4\). We obtain \[ \|G_{n,k}\|_{C^0}\le C_k e^{-3b_ukn/4}. \tag{56}\] Uniform interior room is provided by Lemma 52; nested intermediate patches allow both this step and the next one.

For a prescribed \(j\), apply Lemma 54 with (56) and the naive bound at an order \(q>j\). Choose \(q\) so large that \[\frac{j}{q}\left(A_k+\frac{3b_uk}{4}\right) \le\frac{b_uk}{4}.\] The resulting derivative bound is \(C_{j,k}e^{-b_ukn/2}\). Finitely many small values of \(n\) can be included by increasing the prefactors. Thus, for example, \(c_0=b_u/4\) works for every fixed \(j,k\). The choices of \(\epsilon,\delta,p,q\) may depend on the orders, which is consistent with the assertion. ◻

Proof. Fix a compact set of \(y\) values. For every fixed \(R\), the displacement \(V_{n,R}\) tends uniformly to zero by Proposition 57, so the exponential chart expressions are defined for all sufficiently large \(n\).

First compare consecutive approximants at time \(-n\). With \(s=f^{-n}a\) and \(s_-=f^{-1}s\), the centered one-step map is \[\Psi_{s_-}(z)=\log_s\bigl(f(\exp_{s_-}z)\bigr).\] It is uniformly smooth, vanishes at zero, and has uniformly bounded derivatives of each fixed order. Exact formal equivariance says that the degree-\(R\) Taylor truncation in \(y\) of \(\Psi_{s_-}(V_{n+1,R}(s_-,y))\) is \(V_{n,R}(s,y)\). Taylor-expand \(\Psi_{s_-}\) in \(z\) through spatial order \(R\). Its remainder has size at most \(C_R|V_{n+1,R}|^{R+1}=O(e^{-c_0(R+1)n})\). Among the terms in the finite polynomial part, every omitted parameter monomial has degree at least \(R+1\). Each degree-\(k\) input coefficient costs \(e^{-c_0kn}\), so every such product has the same bound. Consequently \[\left|\Psi_{s_-}(V_{n+1,R})-V_{n,R}\right| \le C_R e^{-c_0(R+1)n}.\] Since \(f^n\) has Lipschitz norm at most \(Ce^{An}\), with \(A\) independent of \(R\), \[ \|F_{n+1,R}-F_{n,R}\|_{C^0} \le C_R e^{-\{c_0(R+1)-A\}n}. \tag{59}\] For a fixed sufficiently large \(R\), this is summable. Completeness of the lifted geometric metric gives a uniform continuous limit. If \(R'>R\), the backward difference of the two truncations consists of degrees \(R+1,\ldots,R'\), and its forward push tends to zero by the same estimate. All sufficiently large fixed cutoffs therefore give one and the same limit \(F\).

We next justify differentiating the limit. For every fixed \(j\), Proposition 57 gives \[\|V_{n,R}\|_{C^j(s,y)}\le C_{R,j}e^{-c_0n}.\] Derivatives in \(y\) introduce only constants depending on the finite cutoff and the compact parameter set. In \((a,y)\) variables, the base substitution \(s=f^{-n}a\) costs at most an exponential with rate depending on \(j\). Applying the geometric exponential and then \(f^n\) has the same kind of fixed-order cost. Thus \[ \|F_{n,R}\|_{C^j} \le C_{R,j}e^{B_jn}, \tag{60}\] where \(B_j\) is independent of \(R\). This independence concerns the exponential rate only; \(C_{R,j}\) need not be bounded as \(R\to\infty\).

The already obtained continuous limit allows us to use common output charts on finitely many compactly contained parameter patches. Apply Lemma 54 to consecutive differences, using (59) and (60) at an order \(q>j\). The decay exponent in the \(C^j\) estimate is at least \[\left(1-\frac jq\right)\{c_0(R+1)-A\} -\frac jq B_q.\] For fixed \(j,q\) it tends to infinity with \(R\). Choosing a fixed \(R\) large enough gives summable \(C^j\) differences and any desired exponential tail bound. This proves smoothness and (57). The derivatives through any fixed parameter order of \(F_{n,R}\) at \(y=0\) agree with the formal jets of \(\mathcal F\) once \(R\) is at least that order, by equivariance. Smooth convergence therefore identifies all Taylor coefficients of \(F\) with those prescribed by \(\mathcal F\).

First the based logarithm in (58) is defined for large \(n\). Indeed, if \(\operatorname{Lip}(f^{-n})\le Ce^{K_1n}\), use the \(C^0\) estimate (57) with a rate \(M>K_1+c_0\). It gives \[d\bigl(f^{-n}F(a,y),\exp_{f^{-n}a}V_{n,R}\bigr) \le C_R e^{-(M-K_1)n}.\] The second point is at distance \(O(e^{-c_0n})\) from \(f^{-n}a\), so both points lie in a fixed based-logarithm ball. Composition with \(f^{-n}\) and this logarithm, including the derivatives of their base \(f^{-n}a\), now costs only a fixed-order exponential by Lemma 52. Choose a larger decay rate in (57) to absorb it; this proves (58). In particular, \[d\bigl(f^{-n}F(a,y),f^{-n}a\bigr)\le C e^{-c_0n}\] on each fixed compact parameter set, for all sufficiently large \(n\). The entire backward tail has this estimate. The local unstable-manifold characterization then places the pair in one local unstable plaque at large negative time; invariance gives \(F(a,y)\in W^u(a)\).

A change of geometric logarithm is a uniformly smooth change of displacement coordinates fixing zero. Its finite Taylor remainder has exactly the degree-weight bound used in (59); its forward push vanishes for a sufficiently large fixed cutoff. Thus the limit is independent of these charts. At a fixed base point, a change of input frame amounts to an invertible linear substitution in the same initial vector, and produces the same limit. The same comparison applies on overlapping base patches. Shifting the backward sequence by one step proves equivariance under \(f\). Taking increasing compact parameter sets now defines the asserted map for all finite \(y\). ◻

Proof. Flatten the actual unstable plaque in uniformly bounded geometric coordinates \((x,z)\), with the plaque given by \(z=0\). Orthogonal projection is a smooth operation on the fixed-rank Grassmannian, so its full formal jets inherit coherence and measurable finite bounds. Write \(Q_\beta(x)\) for its pure normal formal derivatives of multi-order \(\beta\) at \((x,0)\). Define \[\widehat Q(x,z) =\sum_{|\beta|\le N}\frac{z^\beta}{\beta!}Q_\beta(x).\] On the plaque this equals \(Q_0(x)\), the actual projection for its smoothly extended subbundle. Holonomicity identifies every mixed derivative through total order \(N\) at the center with the corresponding formal derivative. Its \(C^r\) norm is controlled by a finite collection of tangential derivatives, through order \(r\), of the \(Q_\beta\), together with bounded geometric chart factors. Those are finite on compactly contained fixed plaques by whole-axis smoothness. Their bounds are measurable; increasing finite bounds exhaust a conull set of centers. One bound therefore holds on a positive-measure set. Use slightly smaller plaques to obtain an ambient neighborhood with the claimed uniform geometric size. No common bound or radius for all \(N,r\) is required. ◻

Proof. In the formal affine coordinates of Lemma 47, varying the base in \(P^s\) with covariantly constant input produces the parallel \(P^s\) variation at the translated point. The horizontal lifts have zero bracket curvature. Thus the formal derivative of \(\mathcal F\) in such a base direction is annihilated by the formal projection \(Q\) onto the perpendicular of \(P^s\). This identity only uses the first derivative of the input in that base direction. It therefore applies to our actual parallel frame by axis coherence.

Fix a regular base point \(a\), an input in a fixed compact parameter set, and \(x\in P^s_a\). We may require \(a\) to be backward recurrent to the positive-measure sets in Lemma 59 for every finite pair of orders. This is a conull condition by ergodicity, and holds for a dense set of points on each of the typical stable axes in use. Set \[s=f^{-n}a,\qquad x_n=df^{-n}x, \qquad W_{n,R}(s,y)=\exp_s V_{n,R}(s,y).\] There is a fixed \(K\) with \(|x_n|\le C e^{Kn}|x|\). Equivariance carries the parallel input frame to a parallel frame on the corresponding pulled-back stable slices.

At a good backward return, choose the extension \(\widehat Q_n\) from Lemma 59, with its Taylor jet and derivative bounds taken to sufficiently high orders for a fixed cutoff \(R\). Keep this chosen ambient matrix field fixed when varying \((s,y)\) in the following test. The expression \[\widehat Q_n\bigl(W_{n,R}(s,y)\bigr) D_sW_{n,R}(s,y)[x_n]\] has zero Taylor coefficients in \(y\) through degree \(R\) by the formal parallelism identity. Its Taylor remainder on a compact parameter set satisfies \[ \left|\widehat Q_n(W_{n,R})D_sW_{n,R}[x_n]\right| \le C_R e^{Kn-c_0(R+1)n}|x|. \tag{62}\] Here is the order accounting, including the base derivative. In bounded coordinates put \(H(s,V)=\exp_s V\), and keep \(\widehat Q_n\) fixed. Writing \(V=V_{n,R}(s,y)\) and \(W=H(s,V)\) gives \[\widehat Q_n(W)D_sW[x_n] =A_s(V)x_n+B_s(V)D_sV[x_n],\] where \(A_s(V)=\widehat Q_n(H(s,V))D_sH(s,V)\) and \(B_s(V)=\widehat Q_n(H(s,V))D_VH(s,V)\). Take these geometric derivatives before Taylor expansion in \(V\). Both degree-\(R\) remainders are \(O(|V|^{R+1})\); the second is also multiplied by \(D_sV[x_n]\). Each positive parameter-degree coefficient of \(V\), and of its first stable base derivative, has cost \(e^{-c_0kn}\) in degree \(k\), by Proposition 57. The polynomial terms through parameter degree \(R\) cancel by the formal identity. Every product of parameter degree at least \(R+1\) has total cost at most \(e^{-c_0(R+1)n}\). Their spatial Taylor remainders have the same cost since \(V=O(e^{-c_0n})\). The base variation contributes \(|x_n|\), and no further order-dependent exponential. All remaining constants are finite because \(R\) is fixed and the extension norms are bounded on the chosen return set. This proves (62).

For a prescribed final error rate, first choose a fixed \(R\) large enough for both (62) and (58), and then use the return set providing the required finite extension norms for that \(R\). Recurrence was imposed simultaneously for all finite orders, so these choices do not require a return set uniform in \(R\). By (58) in \(C^1(a,y)\), the point \(f^{-n}F(a,y)\) and its base derivative can replace \(W_{n,R}\) and \(D_sW_{n,R}[x_n]\) with an error decaying at any prescribed exponential rate, upon taking a sufficiently large fixed cutoff. The uniformly bounded test derivatives and the at-most-exponential size of the vectors preserve that assertion. The true backward point is in the local unstable plaque of \(s\) for large \(n\): the whole further backward tail is uniformly close to the tail of \(s\), as shown in Proposition 58. On that actual plaque \(\widehat Q_n\) is the true projection, regardless of whether the polynomial approximant lies in the plaque.

It follows, by increasing the fixed cutoff, that for any prescribed \(L>0\) along the corresponding sequence of good returns, \[\operatorname{dist}\left( df^{-n}D_aF(a,y)[x],\, P^s_{f^{-n}F(a,y)}\right) \le C_L e^{-Ln}|x|.\] The subbundle appearing here is the smooth extension on this actual unstable plaque; its values are contained in the actual \(E^s\) everywhere by continuity. Push forward by \(df^n\) and choose \(L\) larger than its fixed exponential norm bound. The distance from \(D_aF(a,y)[x]\) to \(E^s_{F(a,y)}\) tends to zero, proving (61) at these regular bases. Smoothness of \(F\) and of \(P^s\) along the base axis, together with continuity of \(E^s\), extends the statement to every base point.

Finally any finite segment of a smooth curve tangent to \(E^s\) has its length exponentially contracted under positive iteration. Its endpoints consequently belong to one strong stable leaf. Apply this to translated base paths, and use \(F(a,y)\in W^u(a)\) for their other sides. ◻

Proof. Choose a typical lifted center \(v\) and a nonzero vector in \(P^s_v\). Proposition 51 gives a proper affine ray \(a(t)\), \(t\ge0\), in \(W^s(v)\), with \(a(0)=v\). Parallel-transport a sufficiently small nonzero vector \(V_0\in P^u_v\) along this ray using the partial connection in Lemma 48. On every finite segment this is an ordinary smooth linear differential equation and therefore has a solution \(V(t)\) throughout that segment. These solutions agree on overlaps, so \(V(t)\) is defined for all finite \(t\ge0\).

Use Proposition 58 to set \[z(t)=F_{a(t)}(V(t)),\] where the subscript denotes the intrinsic initial-vector version of that construction. Local parallel frames on the stable slices show that this is a smooth path; compatibility on overlaps makes it globally well defined. Proposition 60 gives \[z(t)\in W^u(a(t))\cap W^s(z(0)) \qquad(t\ge0).\] Only finite parameter values on compact ray segments have been used, so no boundedness of \(V(t)\) as \(t\to\infty\) is needed.

Let \(\eta\) be the future endpoint of \(v\), and \(\eta'\) the future endpoint of \(z(0)\). Because \(D_VF_v(0)\) is the inclusion of \(P^u_v\), a sufficiently small nonzero \(V_0\) has \(z(0)\ne v\). These two vectors belong to one strong unstable leaf. If their future endpoints also coincided, uniqueness of the connecting geodesic and equality of their past Busemann levels would force the vectors to coincide. Therefore \(\eta'\ne\eta\).

Write \(\zeta_t\) for the past endpoint of \(a(t)\). The past-endpoint map identifies the full strong stable leaf \(W^s(v)\) with \(\partial_\infty\widetilde M\setminus\{\eta\}\): its inverse takes the geodesic from \(\zeta\) to \(\eta\) at the fixed \(b_\eta(v)\) level. Existence, uniqueness, and continuity are exactly the endpoint facts just recalled. Properness of \(a(t)\) therefore implies \(\zeta_t\to\eta\). Indeed otherwise a subsequence of its past endpoints would stay in a compact subset of the boundary minus \(\eta\), whose inverse image in \(W^s(v)\) is compact.

Figure 1 records these actual leaf incidences for a finite segment of the two paths.

We claim \[ b_{\zeta_t}(a(t))\longrightarrow-\infty. \tag{63}\] The sum \(b_\eta+b_{\zeta_t}\) is convex and has zero gradient on the geodesic with these endpoints, so it attains its minimum there. Let \(c(r)\) be the ray from \(o\) toward \(\eta\). For each fixed \(r>0\), \[b_\eta(a(t))+b_{\zeta_t}(a(t)) \le b_\eta(c(r))+b_{\zeta_t}(c(r)).\] As \(t\to\infty\) the right-hand side tends to \(-2r\), by endpoint continuity of normalized Busemann functions. The term \(b_\eta(a(t))=b_\eta(v)\) is fixed. Since \(r\) is arbitrary, (63) follows.

On the other hand, \(z(t)\) has ordered endpoints \((\zeta_t,\eta')\) and fixed future Busemann level \(b_{\eta'}(z(0))\). Their limiting pair \((\eta,\eta')\) is off the diagonal. The corresponding unit tangent vectors therefore converge to the vector with this limiting pair and that level; in particular they remain in a compact subset of \(S\widetilde M\). Endpoint continuity now makes \(b_{\zeta_t}(z(t))\) bounded. But \(z(t)\in W^u(a(t))\) means that the two vectors have the same past endpoint and past Busemann level, so \[b_{\zeta_t}(z(t))=b_{\zeta_t}(a(t)).\] This contradicts (63) and excludes the infinite-dimensional alternative. ◻

Proof. Use the graded section coordinates of Lemma 27 and the dilation of Proposition 28. Write \(x\) for the positive variables and \(y\) for the negative variables. The actual central unstable plaque has formal projection \(y=0\). Let \(D_+\) and \(D_-\) denote the restrictions of \(D\) to the two coordinate spaces. The positive formal field furnished by Proposition 45 has value \(E^u\) at the origin and is invariant under the dilation. Its restriction to \(y=0\) is a graph, with graph matrix \(B(x)\) from the positive to the negative space. Dilation invariance gives \[B(e^{tD_+}x)=e^{tD_-}B(x)e^{-tD_+}.\] Diagonalize \(D\) over the real graded spaces. A monomial \(x^\nu\) in a matrix entry with input weight \(\lambda_i>0\) and output weight \(\mu_j<0\) would have to satisfy \[\sum_a\nu_a\lambda_a=\mu_j-\lambda_i<0.\] The left side is nonnegative, so every coefficient vanishes. Thus \(B=0\). This is an identity of formal series; it needs no convergence. Adding the flow direction gives formal tangency to the weak unstable plaque.

In ordinary smooth coordinates, represent that plaque by a graph \(z=g(v)\) with prescribed crossing \(g(0)\). Tangency to the formal plane field gives a system \[\partial_{v_i}g(v)=A_i(v,g(v)).\] Its constant coefficient determines the first derivatives of \(g\); after \(d\) differentiations, it determines derivatives of order \(d+1\) from the field jet and the already determined graph coefficients. There is no existence assertion hidden in this recursion: the actual central plaque already provides the formally tangent solution. The recursion merely expresses each of its graph coefficients by finitely many operations on the ambient formal field jet.

The full ambient jets of that field are smooth and holonomic on each entire typical stable axis, by Proposition 45. Hence the recursively obtained graph coefficients are smooth there and have the same tangential coherence. Initially the equality with actual plaque coefficients holds at the regular crossing points. It holds everywhere on the axis by continuity: the actual plaques depend continuously in \(C^0\) on their crossings, have uniformly bounded local derivatives of each finite order, and therefore depend continuously in each finite \(C^r\) norm by interior interpolation. Regular crossings are dense on the axis by conditional full measure. This also explains why no transverse smoothness of the actual foliation was used. Flow transport proves the weak-axis assertion, with compatibility on overlapping flow boxes. ◻

Proof. At any fixed Busemann height, the strong stable slice in \(\mathcal B\) is a smooth transversal to the weak unstable plaques. Use coordinates \(a\) on this slice and auxiliary variables \(z\) along the plaques. By Lemma 62, the parameterized actual plaque Taylor graphs form a formal map \[\Psi(a,z),\qquad \Psi(a,0)=a,\] whose coefficients in \(z\) are smooth in \(a\), with all their tangential derivatives. Within a weak unstable plaque, the strong plaque through its crossing is tangent to \(\ker\alpha\). Its graph coefficients are determined recursively from the weak plaque jets and the smooth contact form, as in Lemma 62; the actual strong plaque supplies the solution. These finite operations preserve smooth dependence on \(a\). We may therefore choose \(z=(u,t)\) and write \[\Psi(a,u,t)=\phi_t\Gamma(a,u),\qquad \Gamma(a,0)=a,\] where \(\Gamma(a,u)\) is the actual strong unstable plaque Taylor map through \(a\). Only its finite Taylor coefficients and their smooth dependence on \(a\) are used. The differential of \((a,z)\mapsto\Psi(a,z)\) at \(z=0\) is invertible, by transversality. Formal inversion therefore determines \(a=a(p)\) and the scalar formal function \(q_0|_{\mathrm{slice}}(a(p))\) to every finite order. Concretely, to compute order \(r\), take a sufficiently high finite Taylor polynomial of \(\Psi\) in \(z\) and invert it as an ordinary local smooth map. The resulting order-\(r\) jet is independent of all higher coefficients. Thus this construction uses only finitely many jet operations at each stage.

Smooth dependence of the coefficients on \(a\) gives tangential holonomicity of the inverse jets. The construction is independent of coordinates because it describes the crossing of an actual plaque with the chosen transversal. Moving the height slice by the flow transports both the actual plaque jets and the prescribed labels; formal inversion commutes with this transport. This proves compatibility at different heights. The inverse label is independent of \(u,t\), and \(\partial_t\Psi=X\circ\Psi\) gives the full formal identity \(Xq=0\). Differentiating scalar jet equivariance in time therefore gives the holonomic derivative identity along \(X\); combined with coherence on each height slice, this proves full holonomicity along \(\mathcal B\). In particular, the scalar jet identity is \[j^r_p q=(j^r_{\phi_t p}q)\circ j^r_p\phi_t, \qquad p,\phi_t p\in\mathcal B,\] where both \(q\)-jets in this display mean the formal jets just constructed. ◻

Proof. The crossing \(a(v)\) in (64) depends continuously on \(v\), so \(a(K)\) is compact in \(A\). Put \(\zeta=v_-\) and choose the continuous bounded time shift \[\tau(v)=b_\zeta(\pi v)-b_\zeta(\pi a(v)).\] The vector \(z(v)=\phi_{\tau(v)}a(v)\) has the same past endpoint and the same past Busemann height as \(v\); hence it is on its strong unstable leaf. Moreover their intrinsic strong unstable distance is uniformly bounded on \(K\). To see the uniformity explicitly, join their base points by a spatial geodesic segment of bounded length and project it onto their common past horosphere along the Busemann gradient flow. The required flow times are bounded by the segment length. Uniform bounds for the Busemann Hessian bound the length of this projected path and its lift to the normal graph. This uses the fixed negatively curved metric, not regularity of the endpoint map.

Backward strong unstable contraction now gives \[d^u(f^{-n}v,f^{-n}z(v))\leq Ce^{-bn}.\] Write \(f^{-n}v=(x_n,v_n)\) and \(f^{-n}z(v)=N_\eta(y_n)\). The uniform local bounds of \(N_\eta\) imply \(d(N_\eta(x_n),N_\eta(y_n))\leq Ce^{-bn}\) and the asserted angular bound. If \[\epsilon_n=b_\eta(x_n)-b_\eta(y_n),\] then \(|\epsilon_n|\leq Ce^{-bn}\). The two vectors \(\phi_{\epsilon_n}N_\eta(x_n)\) and \(N_\eta(y_n)\) belong to the same strong stable leaf and have intrinsic distance \(O(e^{-bn})\). Forward strong stable contraction, and the unchanged small time shift, show that \(f^nN_\eta(x_n)\) stays within \(O(e^{-bn})\) in the intrinsic weak stable metric of the compact family \(z(v)\). This proves the compact-axis assertion.

For completeness, the local crossing assertion has a uniform transversality bound. At \(t=0\), the weak unstable plaque intersected with \(b_\eta=r\) has tangent \(E^u\). In Jacobi coordinates the differential of its angular coordinate in a spatial direction \(a\) is \((U(v)+U(-v))a\), which is uniformly positive definite. Uniform smooth bounds on individual plaques and the inverse function theorem therefore give actual graphs \[ s=P_{s_0,r}(t),\qquad P_{s_0,r}(0)=s_0, \tag{65}\] on a fixed radius, with uniformly bounded derivatives in \(t\) of each fixed order. No derivative with respect to \(s_0\) is asserted.

The nearby crossing of the plaque through \(f^{-n}v\) is especially explicit. If \(c_n=b_\eta(x_n)\), then \[c_n=nT-\tau(v)+\epsilon_n, \qquad a_n(v)=\phi_{-c_n}a(v), \qquad f^na_n(v)=\phi_{\tau(v)-\epsilon_n}a(v).\] These identities follow from its past endpoint, future endpoint, and height. They prove compactness of its forward image. On the local graph between angular values \(0\) and \(t(f^{-n}v)\), the base points move by \(O(e^{-bn})\) on this same horosphere. Their axis vectors have that small intrinsic stable distance from \(a_n(v)\); contraction gives the remaining compactness assertion. Replacing \(K\) by a slightly larger compact subset of \(\mathcal U\) proves the neighborhood assertion. ◻

Proof. Transport the full ambient scalar jet of order \(l\) from \(f^np\) to \(p\), using Lemma 63. The one-step transport is a finite-dimensional linear pullback-jet matrix \(M_l(p)\). In uniformly bounded geometric charts these matrices and every fixed-order weak-axis derivative are uniformly bounded: the one-step ambient map is smooth, and the weak-axis charts have the uniform leafwise bounds just explained. The transport over \(n\) steps is the product of these matrices at \(p,fp,\ldots,f^{n-1}p\).

When this product is differentiated \(j\) times along the initial weak axis, at most \(j\) factors are differentiated. Each derivative of a factor composed with \(f^i|_{\mathcal B}\) is bounded by a constant depending on \(l,j\), uniformly in \(i\). The remaining factors are bounded by a fixed \(M_l\geq1\). The number of terms is at most a constant times \((1+n)^j\). The final jet field on \(L\), and its derivatives composed with \(f^n|_{\mathcal B}\), have bounded derivatives. Hence the product rule gives \(C_{L,l,j}(1+n)^jM_l^n\). Absorb the polynomial in \(e^n\), with its \(j\)-dependent prefactor. Taking angular coefficients and changing uniform frames preserve the resulting exponent \(A_l=1+\log M_l\). This is Lemma 53 on the fixed bundle \(J^l\) of scalar \(l\)-jets: differentiation changes the matrix coefficient functions, not the jet order. Using ambient derivative bounds for \(f^n\) in its place would lose the independence from \(j\). ◻

Proof. We induct on \(l\), proving all derivative orders at each stage. For \(l=0\), the identity \(q_0(p)=q_0(f^np)\) and the bounded weak-axis derivatives of forward iterates give bounded derivatives on \(L_n\).

At each fixed height, formal constancy on the actual graph (65) gives the coefficientwise identity \[ \sum_{k\geq0}q_k(P_{s_0,r}(t),r)[t^k]=q_0(s_0,r). \tag{68}\] Every occurrence of a coefficient on the left is Taylor-expanded at \(t=0\). The identity follows directly from the formal crossing inversion of Lemma 63 and its tangential holonomicity. In degree \(l\), its only term containing \(q_l\) is \(q_l(s_0,r)[t^l]\). All other terms involve derivatives in \(s\) of \(q_k\) with \(k<l\) and order at most \(l-k\), multiplied by polynomials in the \(t\)-derivatives of \(P_{s_0,r}\) through order \(l\). These graph derivatives are uniformly bounded. The induction hypothesis and polarization of the homogeneous degree-\(l\) identity therefore give \[\|q_l\|_{C^0(L_n)}\leq C_{L,l,\varepsilon}e^{\varepsilon n} \quad\hbox{for every }\varepsilon>0.\] This argument estimates only values of the next coefficient. In particular it never differentiates the actual graph with respect to its crossing point.

We upgrade values to arbitrary weak-axis derivatives using interpolation on the full weak axis, including the height coordinate. On a fixed interior chart the elementary Taylor difference estimate is \[\|D^ju\|_{C^0} \leq C_{j,r}\bigl(h^{-j}\|u\|_{C^0} +h^{r-j}\|u\|_{C^r}\bigr), \qquad 0<h<h_0,\] where the norms on the right are on a slightly larger chart. The uniform interior room makes its constants independent of \(n\); see also Lemma 54. Fix \(j\geq1\) and target \(\varepsilon>0\). Take \(\delta=\varepsilon/(4j)\) and \(h=e^{-\delta n}\), use the value bound with exponent \(\varepsilon/4\), and choose a fixed integer \(r>j+A_l/\delta\). By (66), the two terms then have exponents at most \(\varepsilon/2\) and \(A_l-\delta(r-j)<0\), respectively. This proves (67). The choice of \(r\) is possible precisely because \(A_l\) does not depend on \(r\). The argument works for every compact \(L\), using a larger compact for interior estimates, so the induction is complete. ◻

Proof. Fix nested compact neighborhoods \(K\Subset K^+\Subset\mathcal U\). For large \(n\), Lemma 64 puts \(f^{-n}K^+\) in the angular coordinate neighborhood of the axis. For a fixed finite cutoff \(R\) define the smooth function, near \(K^+\), \[ p_{n,R}(v)=\sum_{l=0}^R q_l\bigl(\rho(f^{-n}v)\bigr) \bigl[t(f^{-n}v)^l\bigr]. \tag{69}\] These expressions are intrinsic tensor contractions, so no seams from angular frames occur. They use a finite polynomial at each \(n\); no infinite formal series is evaluated.

We first estimate their error without differentiating \(q\). At a pulled back point, use its actual plaque graph \(s=P_{s_0,r}(t)\) and its genuine crossing \((s_0,r,0)\). The function of the angular parameter \[H(t)=\sum_{l=0}^R q_l(P_{s_0,r}(t),r)[t^l]-q_0(s_0,r)\] has zero Taylor coefficients through degree \(R\) by (68). Along the segment from zero to the angular parameter in question, all underlying axis points lie over the compact set supplied by Lemma 64. Lemma 66, the uniform derivatives of the single actual graph, and Taylor’s remainder formula give \[|H(t)|\leq C_{K,R,\varepsilon}e^{\varepsilon n}|t|^{R+1}.\] The actual label equals \(q_0(s_0,r)\) on that plaque, and is flow invariant. Consequently \[ \|p_{n,R}-q\|_{C^0(K^+)} \leq C_{K,R,\varepsilon} e^{-(b(R+1)-\varepsilon)n}. \tag{70}\] Only derivatives of the smooth axis coefficients and of an individual actual plaque entered this estimate.

For each fixed \(j\), the same coefficients have subexponential \(C^j\) bounds. A polynomial of degree at most \(R\) in \(|t|<1\) therefore has a \(C^j\) bound \(C_{K,R,j,\varepsilon}e^{\varepsilon n}\) in the geometric angular coordinates. The ordinary ambient bounds for \(f^{-n}\) and the chain rule give \[ \|p_{n,R}\|_{C^j(K^+)} \leq C_{K,R,j,\varepsilon}e^{(B_j+\varepsilon)n}, \tag{71}\] where \(B_j\) depends on \(j\) and the smooth dynamics, but is independent of \(R\). Indeed differentiating the angular polynomial at most \(j\) times only introduces \(R\)-dependent combinatorial constants; its coefficients already have arbitrary subexponential slack. All exponential losses come from derivatives of \(f^{-n}\) of orders at most \(j\). This independence from \(R\) is the second essential exponent separation in the proof.

Here are the convergence quantifiers. Fix a desired derivative order \(j\geq1\) and a desired positive convergence rate \(\lambda\). Set \(r=2j+2\), take \(\varepsilon=1\), and choose one finite \(R\) large enough that \[\bigl(b(R+1)-1\bigr)(1-j/r) -(B_r+1)j/r>\lambda.\] The differences \(p_{n+1,R}-p_{n,R}\) have the \(C^0\) decay in (70) and the \(C^r\) growth in (71). Interior interpolation on the fixed pair \(K\Subset K^+\) gives \[\|p_{n+1,R}-p_{n,R}\|_{C^j(K)}\leq C e^{-\lambda n}.\] Thus the sequence converges in \(C^j(K)\) to its already identified uniform limit \(q\). Since \(j\) is arbitrary, \(q\) is smooth. There is no requirement that a single cutoff work for all \(j\). If a single sequence converging in the \(C^\infty\) topology is desired, choose cutoffs \(R_k\) which work through order \(k\), and choose \(n_k\) so large that \(\|p_{n_k,R_k}-q\|_{C^k(K)}\leq2^{-k}\). This diagonal sequence has that property and still evaluates only finite polynomials. ◻

Proof. Strong unstable holonomy between smooth transversals is absolutely continuous with locally bounded positive Jacobian; applying the same statement to inverse holonomy supplies the reciprocal bound. This is the bounded-Jacobian theorem for horospheric foliations (Ballmann 1995, Appendix, Theorem 5.1 and Definition 3.3). We explain the passage to the weak foliation.

Thicken each transversal for a short interval by the flow, using a fixed time translate of the target if necessary. The resulting smooth \((m+1)\)-dimensional submanifolds are transversals to the strong unstable foliation. On suitably small flow boxes, strong holonomy has the form \[H(a,t)=(h(a),t+\tau(a)).\] This follows from flow invariance and uniqueness of the crossing; \(\tau\) need only be continuous. Choose a fixed short interval \(I\) so that \(H\) is defined on \(E\times I\) for every \(E\) in a smaller source patch. In product coordinates, its image has the interval \(I+\tau(h^{-1}(y))\) above each \(y\in h(E)\). Fubini’s theorem gives product volume \(|I|\mathop{\mathrm{Vol}}(h(E))\). The smooth volume densities of both flow-box parametrizations are bounded above and below. The strong holonomy comparison applied to \(E\times I\) proves the displayed weak comparison. A finite chain of boxes proves the same assertion along any compact holonomy path. ◻

Proof. Fix a base point \(x\) on the universal cover. On its direction sphere, the restriction \(q|_{S_x\widetilde M}\) is locally the weak unstable holonomy to \(A\), followed by the smooth coordinate map \(q_A\). Both are genuine \(m\)-dimensional transversals: a weak unstable tangent vector with zero spatial projection is zero, and \(TA=E^s\) is complementary to \(E^u\oplus\mathbb RX\). The restriction is a local homeomorphism, by the endpoint parametrizations. By Proposition 67 it is smooth, and by Lemma 68 it satisfies a positive lower volume comparison on every smaller patch.

Its derivative cannot be singular. Use smooth coordinate charts on the source sphere and the target, with volume densities bounded above and below. If the derivative at a point had rank less than \(m\), choose a unit target direction perpendicular to its image. Differentiability then places the image of a coordinate ball of radius \(r\) in a box with widths \(O(r)\) in the other \(m-1\) directions and \(o(r)\) in that direction. Its volume is \(o(r^m)\), contradicting the fixed positive lower volume comparison for these balls. Thus \(dq\) is surjective at every point of \(\mathcal U\). The smooth bundle \(\ker dq\) has dimension \(m+1\) and contains the tangent spaces of the actual weak unstable leaves, since \(q\) is constant on each such smooth leaf. The dimensions agree, so \[\ker dq=E^u\oplus\mathbb RX \quad\hbox{on }\mathcal U.\] Intersecting with the smooth contact distribution gives the smooth strong unstable bundle \[E^u=(\ker dq)\cap\ker\alpha.\]

The coordinate patches of the full typical strong stable axis cover all past endpoints except \(\eta\). A second typical axis with a distinct future endpoint covers the omitted endpoint. Such an axis exists because typical points have full smooth measure, whereas a single fixed-future weak leaf has positive codimension. The preceding argument therefore gives smoothness of \(E^u\) everywhere on the cover. The flip \(v\mapsto-v\) exchanges the strong bundles, and hence gives smoothness of \(E^s\). These are the actual invariant bundles, so they are deck invariant and descend smoothly to \(SM\). ◻