Escaping Lagrangians on regular energy surfaces

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

Working question: Can an escaping energy family have a single-valued action on a shell with loops? Local closed one-forms need not be exact on a space with loops. The normalized Hamilton construction provides a global action rather than appealing to local primitives. The family must also be an actual embedded Lagrangian, because a phase description is not useful if different rays project to the same undesired data.

A fixed energy selects the frequencies of a free particle. Its outgoing rays form a normal bundle of that energy surface. We construct the corresponding escaping family for a long-range force, prove that it is embedded, and find a globally normalized action. The action remains single-valued even when the energy surface has loops.

Read Hamilton trajectories under a long-range force, including its spatial and mixed time estimates and its local contraction and continuation proofs. The written coordinate inverse proof, (CI1)–(CI3) supplies the local inverses and smooth level-set charts; the finite bump construction supplies compact cutoffs. Section 3 below proves the cotangent form convention, Hamilton preservation of the symplectic form and the isotropic dimension criterion for this Lagrangian flowout directly.

The free texts of Teschl [T] and Oh [O] provide background on characteristic flows.

1. The free geometry and the initial sheet

Let P0P_0 be a real polynomial with ∣P0(ξ)∣→∞|P_0(\xi)|\to\infty as ∣ξ∣→∞|\xi|\to\infty, and let λ\lambda be a regular value. Write

Mλ={ξ:P0(ξ)=λ},v(ξ)=∇P0(ξ),H=P0+VL.(1) \begin{gathered} M_\lambda=\{\xi:P_0(\xi)=\lambda\},\\ v(\xi)=\nabla P_0(\xi),\\ H=P_0+V_L. \end{gathered} \tag{1}

The free shell MλM_\lambda is compact. The long-range term VLV_L is real and smooth, and is a polynomial of any finite order in ξ\xi. Fix an integer κ≥2\kappa\ge2, 0<δ<1/(κ+1)0<\delta<1/(\kappa+1), and assume, on each compact frequency set KK,

∣∂ξα∂xβVL(x,ξ)∣≤Cαβ,K(1+∣x∣2)−m(∣β∣)/2,m(j)=j+δ(0≤j≤κ),m(j)=1+(κ−1+δ)j/κ(j≥κ).(2) \begin{gathered} |\partial_\xi^\alpha\partial_x^\beta V_L(x,\xi)| \le C_{\alpha\beta,K}(1+|x|^2)^{-m(|\beta|)/2},\\ m(j)=j+\delta\quad(0\le j\le\kappa),\\ m(j)=1+(\kappa-1+\delta)j/\kappa\quad(j\ge\kappa). \end{gathered} \tag{2}

The formulas agree at j=κj=\kappa. We use x′=Hξx'=H_\xi, ξ′=−Hx\xi'=-H_x, and the frequency-base cotangent form β=∑jxj dξj\beta=\sum_jx_j\,d\xi_j. All derivative estimates use ordinary real derivatives.

For VL=0V_L=0, the outgoing family is

N+(Mλ)={(t v(ξ),ξ):t>0, ξ∈Mλ}.(3) \mathcal N_+(M_\lambda) =\{(t\,v(\xi),\xi):t>0,\ \xi\in M_\lambda\}. \tag{3}

Its restriction of β\beta is zero: on tangent directions to the shell it is t dP0=0t\,dP_0=0, and the frequency is constant along a ray. For P0(ξ)=∣ξ∣2P_0(\xi)=|\xi|^2, λ>0\lambda>0, these are the outward normal rays of the sphere of radius λ\sqrt\lambda.

For the perturbed energy equation, start far along the free graph. It matters which frequency sheet is selected. Choose a small regular energy collar U={∣P0−λ∣<ε}U=\{|P_0-\lambda|<\varepsilon\}, and define

MT={η∈U:H(Tv(η),η)=λ},ΣT={(Tv(η),η):η∈MT}.(4) \begin{gathered} M_T=\{\eta\in U:H(Tv(\eta),\eta)=\lambda\},\\ \Sigma_T=\{(Tv(\eta),\eta):\eta\in M_T\}. \end{gathered} \tag{4}

Only this near-shell sheet is used below. Here is why that restriction is needed. Take n≥2n\ge2, P0=∣ξ∣2P_0=|\xi|^2, λ=1\lambda=1, and a real smooth compactly supported χ\chi that equals one on ∣x∣≤1|x|\le1, vanishes on ∣x∣≥2|x|\ge2, and satisfies 0≤χ≤10\le\chi\le1. Set VL=χ(x)V_L=\chi(x). In ∣ξ∣<1/(2T)|\xi|<1/(2T), the unrestricted energy equation becomes

∣ξ∣2+χ(2Tξ)−1=∣ξ∣2=0.(5) |\xi|^2+\chi(2T\xi)-1=|\xi|^2=0. \tag{5}

Thus its full frequency set has an isolated point at zero. Its sheet near the unit sphere is exactly the sphere for large TT. The isolated point starts the stationary Hamilton orbit (0,0)(0,0). An energy-one orbit in ∣x∣≤1|x|\le1 has ξ=0\xi=0 and is stationary, so ODE uniqueness prevents other initial orbits from entering that region. The unrestricted flowout has an isolated point and cannot be an nn-dimensional Lagrangian. All decay assumptions hold for this smooth compactly supported potential.

The counterexample shows why an unrestricted initial set is insufficient. The regular-collar construction below retains the required restriction.

Theorem 1.1 (an escaping energy family and its action). For all sufficiently large TT, MTM_T is a compact smooth hypersurface diffeomorphic to MλM_\lambda. Its Hamilton flow exists for every t≥Tt\ge T, has energy λ\lambda, and satisfies

ξ(t)−η=O(T−δ),x(t)=t(v(ξ(t))+z(t)),z(t)=O(t−δ),c t≤∣x(t)∣≤C t.(6) \begin{gathered} \xi(t)-\eta=O(T^{-\delta}),\\ x(t)=t\bigl(v(\xi(t))+z(t)\bigr),\\ z(t)=O(t^{-\delta}),\\ c\,t\le|x(t)|\le C\,t. \end{gathered} \tag{6}

The constants are uniform over the initial sheet. All future frequencies lie in any prescribed neighborhood of MλM_\lambda if TT is sufficiently large. The open future union ΛT+\Lambda_T^+, with t>Tt>T, is an embedded Lagrangian. It has a unique real smooth action ψ\psi, smooth up to its initial boundary, with

dψ=β∣ΛT+,ψ(T,η)=−TVL(Tv(η),η).(7) d\psi=\beta|_{\Lambda_T^+},\qquad \psi(T,\eta)=-T V_L(Tv(\eta),\eta). \tag{7}

In a fixed chart where MλM_\lambda is ξ1=E(ξ′)\xi_1=E(\xi'), the perturbed sheet is ξ1=ET(ξ′)\xi_1=E_T(\xi'), and every multiindex satisfies

∣∂ξ′α(ET−E)∣≤CαT∣α∣−m(∣α∣).(8) |\partial_{\xi'}^\alpha(E_T-E)| \le C_\alpha T^{|\alpha|-m(|\alpha|)}. \tag{8}

There is an analogous negative-normal construction. If MλM_\lambda is empty, the near-shell sheet and its flowout are empty. The proof is given in Sections 2–5.

2. A collar built by the normal flow

If MλM_\lambda is empty, the sheet in a sufficiently small compact regular energy band is empty and its flowout is empty. All local chart assertions are vacuous. Assume now Mλ≠∅M_\lambda\ne\varnothing.

Polynomial properness makes every bounded energy band compact. Since λ\lambda is regular, choose ε>0\varepsilon>0 so small that {∣P0−λ∣≤2ε}\{|P_0-\lambda|\le2\varepsilon\} has no critical point. To justify one such choice, a sequence of critical points with values tending to λ\lambda would be confined to a compact band; its convergent subsequence would give a critical point of value λ\lambda, a contradiction. There is therefore a uniform positive lower bound for ∣∇P0∣|\nabla P_0| on this band.

Use the smooth normal field ν(ξ)=∇P0(ξ)∣∇P0(ξ)∣2.(9) \nu(\xi)=\frac{\nabla P_0(\xi)}{|\nabla P_0(\xi)|^2}. \tag{9} Its flow N(u,s)N(u,s), starting at u∈Mλu\in M_\lambda, obeys P0(N(u,s))=λ+sP_0(N(u,s))=\lambda+s. Uniform local existence on the compact band and continuation while the level stays in that band give it for ∣s∣<ε|s|<\varepsilon. The map N:Mλ×(−ε,ε)⟶{∣P0−λ∣<ε}(10) N:M_\lambda\times(-\varepsilon,\varepsilon) \longrightarrow \{|P_0-\lambda|<\varepsilon\} \tag{10} is a diffeomorphism. For injectivity, equal image points have equal ss from their energy values; reversing the normal flow recovers the same uu. For surjectivity, start at a point of the energy band and run the field backward for time P0(ξ)−λP_0(\xi)-\lambda; the energy reaches λ\lambda and the entire trajectory stays in the compact noncritical band. This also gives the smooth inverse explicitly.

On the compact smaller collar, write the perturbed equation as s+RT(u,s)=0,RT(u,s)=VL(T∇P0(N(u,s)),N(u,s)).(11) \begin{gathered} s+R_T(u,s)=0,\\ R_T(u,s)=V_L(T\nabla P_0(N(u,s)),N(u,s)). \end{gathered} \tag{11} Each derivative of RTR_T of order zero or one is a finite sum with at most one factor TT and one position derivative of VLV_L. Since ∣Tv∣|Tv| is comparable with TT, the stated decay bounds give ∥RT∥C1=O(T−δ)\|R_T\|_{C^1}=O(T^{-\delta}). For large TT, the map s↦−RT(u,s)s\mapsto-R_T(u,s) is a uniform contraction of [−ε/2,ε/2][-\varepsilon/2,\varepsilon/2] into itself. It has a unique smooth root sT(u)s_T(u), with ∣sT∣≤CT−δ|s_T|\le CT^{-\delta}. Every zero in the open collar has that small ss, since its equation itself gives ∣s∣≤CT−δ|s|\le CT^{-\delta}; there are no additional zeros within the collar. Thus MT={ξ:H(Tv(ξ),ξ)=λ}∩{∣P0−λ∣<ε}={N(u,sT(u)):u∈Mλ}(12) \begin{aligned} M_T&=\{\xi:H(Tv(\xi),\xi)=\lambda\}\\ &\cap\{|P_0-\lambda|<\varepsilon\}\\ &=\{N(u,s_T(u)):u\in M_\lambda\} \end{aligned} \tag{12} is one smooth compact hypersurface, diffeomorphic to MλM_\lambda. Section 5 proves every graph derivative estimate on the finite compact chart covering. This constructs the complete sheet in the regular collar.

3. Escape makes the flowout embedded

Start only from ΣT={(T∇P0(η),η):η∈MT}.(13) \Sigma_T =\{(T\nabla P_0(\eta),\eta): \eta\in M_T\}. \tag{13} The uniform spatial flow theorem in Hamilton trajectories under a long-range force gives a smooth future solution for all t≥Tt\ge T. The frequency changes by O(T−δ)O(T^{-\delta}) and z=O(t−δ)z=O(t^{-\delta}), since the initial displacement is zero. The initial sheet is O(T−δ)O(T^{-\delta}) from MλM_\lambda. Hence all future frequencies stay within any given fixed neighborhood of MλM_\lambda, once TT is sufficiently large. They also stay in one compact regular-frequency set on which c0t≤∣x(t)∣≤C0t.(14) c_0t\le |x(t)|\le C_0t. \tag{14} These follow from x=t(∇P0(ξ)+z)x=t(\nabla P_0(\xi)+z), the positive free-velocity lower bound and the bounded free velocity on that compact set. The identity dH/dt=Hx⋅Hξ−Hξ⋅Hx=0dH/dt=H_x\cdot H_\xi-H_\xi\cdot H_x=0 gives conservation of the Hamilton energy H=P0+VLH=P_0+V_L, so every orbit has energy λ\lambda.

In these coordinates the two-form is ω=∑dξj∧dxj\omega=\sum d\xi_j\wedge dx_j, so ω(XH,w)=−dH(w)\omega(X_H,w)=-dH(w). We prove the preservation needed here directly. Order a variation as a=(ax,aξ)a=(a_x,a_\xi), and set

Q=(0−II0),ω(a,b)=aTQb. Q=\begin{pmatrix}0&-I\\I&0\end{pmatrix},\qquad \omega(a,b)=a^{\mathsf T}Qb.

The Hamilton vector is XH=−Q∇HX_H=-Q\nabla H. Two transported variations satisfy a′=−QH′′aa'=-QH''a and b′=−QH′′bb'=-QH''b, with the Hessian evaluated along the orbit. Its symmetry and Q2=−IQ^2=-I give

ddt(aTQb)=aT((−QH′′)TQ+Q(−QH′′))b=aT(−H′′+H′′)b=0. \frac d{dt}(a^{\mathsf T}Qb) =a^{\mathsf T}\bigl((-QH'')^{\mathsf T}Q+Q(-QH'')\bigr)b =a^{\mathsf T}(-H''+H'')b=0.

Thus the smooth Hamilton flow preserves this two-form on every finite interval of its existence. This calculation also fixes its sign without importing a global generating-function theorem. Teschl [T] develops the same Hamiltonian and canonical-coordinate background; the flow and all estimates required at infinity are supplied by the preceding lesson.

The initial graph x=T∇P0(η)x=T\nabla P_0(\eta) is isotropic: its Hessian is symmetric, so the cotangent two-form vanishes on the full graph, hence on its (n−1)(n-1)-dimensional energy sheet. The Hamilton field is transverse to the initial sheet. A tangent vector to the full graph obeys δx=TP0′′(η)δη\delta x=T P_0''(\eta)\delta\eta. The Hamilton vector would obey that identity only if ∇P0+VL,ξ+TP0′′VL,x=0.(15) \nabla P_0+V_{L,\xi}+T P_0''V_{L,x}=0. \tag{15} Its left side is ∇P0+O(T−δ)\nabla P_0+O(T^{-\delta}), because VL,ξ=O(T−δ)V_{L,\xi}=O(T^{-\delta}) and VL,x=O(T−1−δ)V_{L,x}=O(T^{-1-\delta}) on the initial sheet. The free velocity is bounded away from zero, so the identity is impossible for large TT. The finite-time flow is locally invertible by following the same autonomous equation backward along the compact orbit segment. The local existence proof applies after reversing time, and uniqueness makes forward and backward maps inverse. Uniqueness also gives ΦsΦr=Φs+r\Phi_s\Phi_r=\Phi_{s+r} wherever both sides are defined; differentiating in rr at zero shows DΦsXH=XH∘ΦsD\Phi_sX_H=X_H\circ\Phi_s. Thus its invertible derivative transports the initial transversality and all tangent directions. Thus its flowout is an immersed nn-dimensional manifold. The two-form vanishes on the transported initial tangent directions by the Hamilton preservation theorem. Its pairing with the Hamilton direction vanishes because every transported direction is tangent to the energy level, and ω(XH,w)=−dH(w)\omega(X_H,w)=-dH(w). It is therefore Lagrangian at every future point. Indeed the displayed matrix QQ is invertible, so for a subspace LL the symplectic orthogonal has dimension 2n−dim⁡L2n-\dim L: the map w↦ω(w,⋅)∣Lw\mapsto\omega(w,\cdot)|_L onto L∗L^* is surjective by nondegeneracy and extension of a linear functional from a basis of LL. An isotropic nn-plane is contained in its symplectic orthogonal of the same dimension and hence equals it. This is the Lagrangian dimension criterion used here.

To prove global injectivity, rather than only immersion, use the escape coordinate ℓ(x,ξ)=x⋅∇P0(ξ)∣∇P0(ξ)∣2.(16) \ell(x,\xi)=\frac{x\cdot\nabla P_0(\xi)}{|\nabla P_0(\xi)|^2}. \tag{16} It equals TT on the entire initial sheet. Its xx-gradient is v/∣v∣2v/|v|^2, and its frequency derivatives are O(∣x∣)O(|x|) on the compact regular-frequency region. Consequently XHℓ=1+ℓx⋅VL,ξ−ℓξ⋅VL,x=1+O((1+∣x∣2)−δ/2).(17) \begin{gathered} X_H\ell \\ =1+\ell_x\cdot V_{L,\xi} -\ell_\xi\cdot V_{L,x} \\ =1+O((1+|x|^2)^{-\delta/2}). \end{gathered} \tag{17} The trajectory bound ∣x∣≥c0T|x|\ge c_0T makes this at least 1/21/2 everywhere on every future orbit when TT is large. Each such orbit therefore meets the initial section ℓ=T\ell=T only once. If two future parametrizations meet, autonomous ODE uniqueness says that one initial point is a future translate of the other, unless the times agree. Strict increase of ℓ\ell forbids a positive translate between two initial points, and equal time then forces equal initial data. Thus the flowout map is injective.

It is an embedding, not merely an injective immersion. For a convergent sequence of image points, ∣x∣≍t|x|\asymp t bounds the times, and compactness of the initial sheet supplies convergent subsequences of initial data. If the limit point lies in the open future image, strict ℓ>T\ell>T excludes a limit time TT. Every convergent subsequence has the same parameter limit by injectivity, so the inverse is continuous. To see local smoothness of the inverse without an additional rank theorem, choose nn ambient coordinate components for which the derivative of the parametrization has a nonsingular minor. The written inverse theorem applied to these components supplies local parameter coordinates; the other components are smooth functions of them. The already proved continuity of the global inverse excludes other parameter branches from a sufficiently small relative neighborhood. Hence the interior flowout is an embedded Lagrangian, with compact initial boundary if that boundary is included.

The same reasoning works for the negative normal direction by applying the positive construction to −P0,−VL,−λ-P_0,-V_L,-\lambda and reversing the actual time parameter. The free normal becomes −∇P0-\nabla P_0, the same regular collar and properness remain valid, and the same estimates hold with these reversed signs.

4. An action that does not require simple connectivity

The embedded flowout has the global parameters t≥Tt\ge T and initial frequency η∈MT\eta\in M_T. Define ψ(t,η)=−TVL(T∇P0(η),η)+∫Ttx(s,η)⋅ξ′(s,η) ds.(18) \begin{gathered} \psi(t,\eta)\\ = -T V_L(T\nabla P_0(\eta),\eta) \\ +\int_T^t x(s,\eta)\cdot\xi'(s,\eta)\,ds. \end{gathered} \tag{18} This is a real smooth finite-time integral on the whole flowout. Its initial value is the required −TVL-TV_L. On the initial energy sheet, T∇P0(η)⋅dη=T d(P0(η)−λ)=−T d[VL(T∇P0(η),η)].(19) \begin{gathered} T\nabla P_0(\eta)\cdot d\eta \\ =T\,d(P_0(\eta)-\lambda) \\ =-T\,d[V_L(T\nabla P_0(\eta),\eta)]. \end{gathered} \tag{19} The rightmost differential is the total differential along the initial graph, including its position dependence. Thus its initial differential is the restriction of β=∑jxjdξj\beta=\sum_jx_jd\xi_j.

Use local coordinates qq on the initial sheet. Since every orbit has the same conserved energy λ\lambda, Hx⋅xq+Hξ⋅ξq=0.(20) H_x\cdot x_q+H_\xi\cdot\xi_q=0. \tag{20} The exact finite product calculation gives ∂q(x⋅ξ′)−∂t(x⋅ξq)=xq⋅ξ′−x′⋅ξq=−Hx⋅xq−Hξ⋅ξq=0.(21) \begin{gathered} \partial_q(x\cdot\xi')-\partial_t(x\cdot\xi_q) \\ =x_q\cdot\xi'-x'\cdot\xi_q\\ \\ =-H_x\cdot x_q-H_\xi\cdot\xi_q\\ =0. \end{gathered} \tag{21} Integration with the initial differential therefore gives ∂qψ=x⋅ξq\partial_q\psi=x\cdot\xi_q. In the time direction, ∂tψ=x⋅ξ′\partial_t\psi=x\cdot\xi' by definition. These coordinates span the tangent of the embedded flowout, so dψ=β∣ΛT+.(22) d\psi=\beta|_{\Lambda_T^+}. \tag{22} No global simple connectivity of the energy shell is needed. The exact initial differential and the orbit parametrization supply global exactness directly, rather than inferring it merely from a closed form. Uniqueness follows because the difference of two such functions has derivative zero along every orbit and initial value zero. Every point lies on an orbit from the initial section.

5. Every derivative of the initial sheet

Work near ξ0∈Mλ\xi^0\in M_\lambda, choose a coordinate with ∂1P0(ξ0)≠0\partial_1P_0(\xi^0)\ne0, and let ξ1=E(ξ′)\xi_1=E(\xi') be the free sheet. After shrinking to a fixed compact chart, write P0(ξ)−λ=(ξ1−E(ξ′))A(ξ),(23) P_0(\xi)-\lambda=(\xi_1-E(\xi'))A(\xi), \tag{23} where A(ξ)=∫01∂1P0(E(ξ′)+s(ξ1−E(ξ′)),ξ′) ds(24) \begin{gathered} A(\xi)\\ =\int_0^1\partial_1P_0 (E(\xi')+s(\xi_1-E(\xi')),\xi')\,ds \end{gathered} \tag{24} is smooth and bounded away from zero. A fixed regular-velocity tube has ∣∇P0∣≥c0>0|\nabla P_0|\ge c_0>0. The energy equation becomes ξ1−E(ξ′)=FT(ξ),FT(ξ)=−A(ξ)−1VL(T∇P0(ξ),ξ).(25) \begin{gathered} \xi_1-E(\xi')=F_T(\xi),\\ F_T(\xi)=-A(\xi)^{-1}V_L(T\nabla P_0(\xi),\xi). \end{gathered} \tag{25} Put b(k)=k−m(k)b(k)=k-m(k), c(k)=max⁡(b(k),0)c(k)=\max(b(k),0), and let DD denote ordinary coordinate derivatives. Finite chain expansions and monotonicity of j−m(j)j-m(j) give ∣DαFT∣≤CαTb(∣α∣).(26) |D^\alpha F_T|\le C_\alpha T^{b(|\alpha|)}. \tag{26} In particular its value and first derivative are O(T−δ)O(T^{-\delta}). Choose a fixed interval around E(ξ′)E(\xi') within the chart. The map r↦E(ξ′)+FT(r,ξ′)r\mapsto E(\xi')+F_T(r,\xi') preserves the interval and has contraction constant below 1/21/2 for all sufficiently large TT, uniformly in the inner chart. Its geometrically convergent iterates give a unique sheet r=ET(ξ′)r=E_T(\xi'). Difference quotients of the fixed-point equation solve an invertible scalar linear equation; successive differentiation gives a smooth real sheet. The zero-order difference is O(T−δ)O(T^{-\delta}).

The needed exponent partition is b(q)+∑i=1qc(ki)≤b(k),ki≥1,∑iki=k.(27) \begin{gathered} b(q)+\sum_{i=1}^q c(k_i)\le b(k), \\ k_i\ge1,\\ \sum_i k_i=k. \end{gathered} \tag{27} To prove it, put θ=(1−δ)/κ\theta=(1-\delta)/\kappa. Here b(j)=−δb(j)=-\delta for j≤κj\le\kappa, and b(j)=θj−1>0b(j)=\theta j-1>0 for j≥κ+1j\ge\kappa+1. If all input orders are small, monotonicity b(q)≤b(k)b(q)\le b(k) suffices. Otherwise let ℓ≥1\ell\ge1 count the input orders at least κ+1\kappa+1. Their sum is at most k−q+ℓk-q+\ell, hence ∑ic(ki)≤θ(k−q+ℓ)−ℓ.(28) \sum_i c(k_i)\le\theta(k-q+\ell)-\ell. \tag{28} When q≥κ+1q\ge\kappa+1, adding b(q)=θq−1b(q)=\theta q-1 leaves b(k)+ℓ(θ−1)≤b(k)b(k)+\ell(\theta-1)\le b(k). When q≤κq\le\kappa, the difference from b(k)b(k) is at most 1−δ−θq+ℓ(θ−1)≤−δ−θ(q−1)≤0.(29) 1-\delta-\theta q+\ell(\theta-1) \le-\delta-\theta(q-1)\le0. \tag{29} This proves all partitions, including the low total-order case.

For the first derivative, differentiation of ET−E=FT(ET,ξ′)E_T-E=F_T(E_T,\xi') gives, with F1=∂1FT(ET,ξ′)F_1=\partial_1F_T(E_T,\xi'),

(1−F1)D(ET−E)=F1DE+(Dξ′FT)(ET,ξ′). (1-F_1)D(E_T-E) =F_1DE+(D_{\xi'}F_T)(E_T,\xi').

Both terms on the right are O(T−δ)O(T^{-\delta}), and ∣1−F1∣≥1/2|1-F_1|\ge1/2. Thus D(ET−E)=O(T−δ)D(E_T-E)=O(T^{-\delta}), including the derivative of the affine transverse input. This supplies the base case for the higher-order induction.

Inductively, the graph map (ET(ξ′),ξ′)(E_T(\xi'),\xi') has order-jj derivatives O(Tc(j))O(T^{c(j)}): the fixed free graph is bounded, and its difference has the Tb(j)T^{b(j)} estimate. At order k≥2k\ge2, isolate the term containing the highest derivative: (1−∂1FT(ET,ξ′))Dk(ET−E).(30) (1-\partial_1F_T(E_T,\xi'))D^k(E_T-E). \tag{30} All other finite chain terms have outer derivative order q≥2q\ge2 and lower graph derivatives with orders summing to kk, or the first outer derivative times the fixed DkED^kE. The partition inequality bounds the former by CTb(k)CT^{b(k)}; the latter is O(Tb(1))=O(T−δ)O(T^{b(1)})=O(T^{-\delta}), also no larger than CTb(k)CT^{b(k)}. The coefficient is bounded away from zero. Thus every order satisfies ∣Dξ′α(ET−E)∣≤CαT∣α∣−m(∣α∣).(31) |D_{\xi'}^\alpha(E_T-E)|\le C_\alpha T^{|\alpha|-m(|\alpha|)}. \tag{31} Finitely many such charts cover the compact shell. Their unique roots agree with the global normal-graph construction of Section 2, so the constants can be chosen uniformly on each compact inner chart. This finishes Theorem 1.1.

The action will next supply local generating functions. Its global definition fixes their transition constants, while the mixed flow estimates control all their radial and transverse derivatives.

Use the conclusion

Check the embedding and the normalization at the free end. Retain the regular energy collar when restricting the construction; a statement near the shell must not be promoted to arbitrary frequencies.

6. Exercises with complete solutions

Exercise 1 — Basic — free rays, empty shells and critical energies.

Let P0(ξ)=∣ξ∣2P_0(\xi)=|\xi|^2 on Rn\mathbb R^n. Describe the near-shell construction when VL=0V_L=0 for λ>0\lambda>0, λ<0\lambda<0 and λ=0\lambda=0. For positive energy compute the escape coordinate, the restricted cotangent form, and the normalized action. Explain which case is excluded by the regular-value hypothesis.

Solution 1. For λ>0\lambda>0, MλM_\lambda is the sphere of radius λ\sqrt\lambda, and the initial frequency sheet is the same sphere. The initial graph is x(T)=2Tηx(T)=2T\eta. Hamilton's equations give

ξ(t)=η,x(t)=2tη,ℓ(x,ξ)=x⋅2ξ4∣ξ∣2=t.(32) \begin{gathered} \xi(t)=\eta,\\ x(t)=2t\eta,\\ \ell(x,\xi)=\frac{x\cdot2\xi}{4|\xi|^2}=t. \end{gathered} \tag{32}

On tangents to the sphere, β=x⋅dξ=2tη⋅dη=t d∣η∣2=0\beta=x\cdot d\xi=2t\eta\cdot d\eta=t\,d|\eta|^2=0. Its time component is also zero because ξ\xi is constant. The action has initial value zero and remains zero. The outward family is embedded: its frequency is the initial η\eta, and its position determines t=∣x∣/(2λ)t=|x|/(2\sqrt\lambda).

For λ<0\lambda<0, the shell is empty. An empty inverse image satisfies the definition of a regular value, and a sufficiently small energy band is also empty. Thus the near-shell family is empty. For λ=0\lambda=0, the shell consists of the critical point ξ=0\xi=0, since ∇P0(0)=0\nabla P_0(0)=0. This energy is excluded. Its stationary free orbit does not provide an escaping regular-energy family.

Exercise 2 — Intermediate — two long-range initial frequencies on the line.

Take P0(ξ)=ξ2P_0(\xi)=\xi^2, λ=a2\lambda=a^2 with a>0a>0, and VL(x,ξ)=b(1+x2)−δ/2V_L(x,\xi)=b(1+x^2)^{-\delta/2}, where bb is real and 0<δ<1/(κ+1)0<\delta<1/(\kappa+1). Let ηT+\eta_T^+ and ηT−\eta_T^- be the near-shell roots converging to aa and −a-a. Find their first asymptotic corrections and the leading initial action at both roots. Do not assume that a component of the frequency shell has positive dimension.

Solution 2. This potential has all derivatives bounded by Cj(1+∣x∣2)−(j+δ)/2C_j(1+|x|^2)^{-(j+\delta)/2}, which is at least as strong as the required full class. Both free roots are regular. The local sheet theorem on the line gives one root near each of them, with ηT±=±a+O(T−δ)\eta_T^\pm=\pm a+O(T^{-\delta}). Their equation is

(ηT±)2−a2=−b(1+4T2(ηT±)2)−δ/2.(33) (\eta_T^\pm)^2-a^2 =-b(1+4T^2(\eta_T^\pm)^2)^{-\delta/2}. \tag{33}

Because ηT±→±a\eta_T^\pm\to\pm a, Tδ(1+4T2(ηT±)2)−δ/2⟶(2a)−δ.(34) T^\delta(1+4T^2(\eta_T^\pm)^2)^{-\delta/2} \longrightarrow(2a)^{-\delta}. \tag{34}

Factor the left side at the appropriate free root and divide by the other factor. This yields

ηT+=a−b(2a)−δ2aT−δ+o(T−δ),ηT−=−a+b(2a)−δ2aT−δ+o(T−δ).(35) \begin{aligned} \eta_T^+&=a-\frac{b(2a)^{-\delta}}{2a}T^{-\delta} +o(T^{-\delta}),\\ \eta_T^-&=-a+\frac{b(2a)^{-\delta}}{2a}T^{-\delta} +o(T^{-\delta}). \end{aligned} \tag{35}

The prescribed initial action is −Tb(1+4T2(ηT±)2)−δ/2-Tb(1+4T^2(\eta_T^\pm)^2)^{-\delta/2}. At either root it equals −b(2a)−δT1−δ+o(T1−δ).(36) -b(2a)^{-\delta}T^{1-\delta}+o(T^{1-\delta}). \tag{36}

For b=0b=0 all corrections vanish exactly. In dimension one the frequency sheet is a finite set of points; there are no transverse frequency derivatives. Each nonempty component still generates a one-dimensional escaping orbit and has its own prescribed initial action.

Exercise 3 — Intermediate — derivatives at the joining order.

Set κ=3\kappa=3, δ=1/8\delta=1/8, b(j)=j−m(j)b(j)=j-m(j), and c(j)=max⁡(b(j),0)c(j)=\max(b(j),0). Compute b(j)b(j) for 0≤j≤50\le j\le5. State the bounds for derivatives of ET−EE_T-E of orders three, four and five. In the fifth derivative of the fixed-point expression, estimate a term with outer order two and inner orders four and one. Compare this with the permitted fifth-order bound.

Solution 3. For j≤3j\le3, m(j)=j+1/8m(j)=j+1/8, so b(j)=−1/8b(j)=-1/8. For j≥3j\ge3, m(j)=1+1724j,b(j)=724j−1.(37) m(j)=1+\frac{17}{24}j,\qquad b(j)=\frac7{24}j-1. \tag{37}

The two definitions agree at three. Consequently b(4)=16,b(5)=1124.(38) b(4)=\frac16,\qquad b(5)=\frac{11}{24}. \tag{38}

The third difference derivative is O(T−1/8)O(T^{-1/8}), the fourth is O(T1/6)O(T^{1/6}), and the fifth is O(T11/24)O(T^{11/24}). These are estimates for the difference from the fixed free graph; the graph itself has derivatives O(Tc(j))O(T^{c(j)}).

The indicated chain-rule term has an outer second derivative of FTF_T, of size Tb(2)T^{b(2)}, a fourth graph derivative of size Tc(4)T^{c(4)}, and a first graph derivative of bounded size. Its exponent is b(2)+c(4)+c(1)=−18+16=124≤1124=b(5).(39) \begin{gathered} b(2)+c(4)+c(1)=-\frac18+\frac16=\frac1{24} \\ \le\frac{11}{24}=b(5). \end{gathered} \tag{39}

Thus it is smaller than the allowed fifth-order bound. In contrast, a fourth derivative with outer order four and four first graph derivatives may have exponent b(4)=1/6b(4)=1/6. A growing high derivative therefore does not contradict the small value and first derivatives of the perturbed sheet.

Exercise 4 — Advanced — periods on a shell with loops.

Assume n=2n=2, P0=∣ξ∣2P_0=|\xi|^2, λ=1\lambda=1, and the general long-range hypotheses. For a smooth closed curve γ\gamma in MTM_T, let γt\gamma_t be its transported curve in the escaping family at time tt. Prove directly that ∫γtβ=0\int_{\gamma_t}\beta=0 for every t≥Tt\ge T. Explain why closedness of β\beta alone would not prove the existence of a single-valued action on this family.

Solution 4. The initial sheet is diffeomorphic to the circle and can carry loops. Write the curve with a periodic parameter qq, and let x(t,q),ξ(t,q)x(t,q),\xi(t,q) be its Hamilton transport. At the initial time, x(T,q)⋅∂qξ(T,q)=−T ∂q[VL(Tv(ξ(T,q)),ξ(T,q))].(40) \begin{gathered} x(T,q)\cdot\partial_q\xi(T,q) \\ =-T\,\partial_q[V_L(Tv(\xi(T,q)),\xi(T,q))]. \end{gathered} \tag{40}

Its integral over the period is zero. Conservation of the common energy gives ∂t(x⋅ξq)=∂q(x⋅ξt).(41) \partial_t(x\cdot\xi_q) =\partial_q(x\cdot\xi_t). \tag{41}

Integrating this identity over the periodic qq shows that ∂t∫x⋅ξq dq=0\partial_t\int x\cdot\xi_q\,dq=0. The initial zero period therefore remains zero for every tt.

More explicitly the finite-time action ψ(t,η)=−TVL(Tv(η),η)+∫Ttx(s,η)⋅ξs(s,η) ds(42) \begin{gathered} \psi(t,\eta)\\ =-TV_L(Tv(\eta),\eta) \\ +\int_T^t x(s,\eta)\cdot\xi_s(s,\eta)\,ds \end{gathered} \tag{42} is already a globally defined function on the entire parameter cylinder, and its differential is β\beta. A closed one-form on a cylinder need not be exact: on the circle the smooth one-form represented by dθd\theta in angular charts has period 2π2\pi. The energy identity and exact initial differential, rather than a simple-connectivity assumption, eliminate these periods in the escaping family.

Exercise 5 — Advanced — a degenerate component away from the regular shell.

Let n≥2n\ge2, P0(ξ)=∣ξ∣4P_0(\xi)=|\xi|^4, λ=1\lambda=1, and let VL(x,ξ)=χ(x)V_L(x,\xi)=\chi(x), with the cutoff from Section 1. Show that the unrestricted initial frequency set has an isolated point at zero for every TT, while its near-shell sheet is the unit sphere for large TT. Prove that the unrestricted forward flowout contains an isolated stationary point. Check properness, regularity and the decay assumptions.

Solution 5. The polynomial is real and proper, and ∇P0(ξ)=4∣ξ∣2ξ\nabla P_0(\xi)=4|\xi|^2\xi. On the unit sphere this gradient has norm four, so energy one is regular. In ∣ξ∣<(4T)−1/3|\xi|<(4T)^{-1/3}, the initial position 4T∣ξ∣2ξ4T|\xi|^2\xi has norm less than one. The unrestricted initial equation reduces there to ∣ξ∣4+χ(4T∣ξ∣2ξ)−1=∣ξ∣4=0.(43) |\xi|^4+\chi(4T|\xi|^2\xi)-1=|\xi|^4=0. \tag{43}

Thus zero is an isolated point of that frequency set. Near the unit sphere, the initial position lies outside ∣x∣=2|x|=2 for all sufficiently large TT, so the perturbed equation is exactly ∣ξ∣4=1|\xi|^4=1.

At (x,ξ)=(0,0)(x,\xi)=(0,0), both Hamilton components vanish: x′=4∣ξ∣2ξ=0x'=4|\xi|^2\xi=0 and ξ′=−∇χ(0)=0\xi'=-\nabla\chi(0)=0. If any energy-one orbit enters ∣x∣≤1|x|\le1, its energy identity gives ∣ξ∣4+1=1|\xi|^4+1=1, hence ξ=0\xi=0. Every such phase point is stationary. Uniqueness of the smooth ODE prevents a different orbit from reaching it at a finite time. Other initial points have ∣x∣>1|x|>1, so their forward images stay outside that closed position region. In a position neighborhood of zero the unrestricted union therefore contains only its stationary initial point.

A positive-dimensional hypersurface cannot have an isolated point, and an nn-dimensional Lagrangian cannot have this isolated component. The smooth compactly supported potential has all required position derivative bounds for every allowed κ,δ\kappa,\delta, and is independent of frequency. Restricting the initial frequencies to a regular collar of the unit sphere removes the degenerate component and preserves the escaping construction.

References

[T] Gerald Teschl, Ordinary Differential Equations and Dynamical Systems, free author's preliminary edition, April 2012, §§2.2, 2.4 and 2.6. The programme's preceding flow lesson proves the local existence, smooth data dependence and compact continuation used here. Section 3 above supplies the exact symplectic preservation and embedding arguments.

[O] Sung-Jin Oh, Lecture Notes for Math 222A, free evolving lecture notes, University of California, Berkeley, Fall 2023, §2.4.1, pp. 23–24. The Hamilton characteristic equations (2.20) provide the comparison; the global sheet, escaping embedding and normalized action are proved in Sections 2–5 above.