Tangency, gliding and generalized characteristic motion
Transverse reflection exchanges two distinct normal covectors. At tangency those covectors coincide, and nearby reflected rays can either converge to an ordinary tangent ray or remain close to the boundary for an increasing number of reflections. The sign of the second normal acceleration distinguishes these possibilities. Higher contact requires estimates that compare the full motion with its tangential gliding flow.
We develop the geometric flow, its compactness and existence, finite-contact uniqueness and a smooth infinite-contact example with two different continuations. The principal symbol may have stationary tangential points. The geometric entry hypothesis is a nondegenerate double normal root; no positive energy or definiteness of the entire symbol is imposed. The quantitative proofs use a fully justified symplectic quadratic normal form. For the original scalar second-order differential operators, the preceding normal-coordinate argument already gives that form without a general preparation theorem.
Read Sections 1–2 of Phase space and generating families for the symplectic sign and linear algebra, Section 2 of Prescribed canonical coordinates and isotropic fibers for full ordinary coordinate completion, and Section 4 of Real and complex symplectic normal forms for the parameter-aware cutoff-series construction used below. The actual smooth local flow, every parameter derivative and inverse-flow argument are supplied in Section 1 of Energy and existence with a timelike Dirichlet boundary; only its generic local smooth-field proof is used here. Section 1 of Boundary reflection and compressed wavefronts supplies the normal form for real quadratic principal symbols and the invariant transverse reflection. The present lesson supplies the additional preparation, comparison, compactness and existence arguments.
The mathematical source is the Springer 2007 edition of Hörmander III, ISBN 978-3-540-49938-1, Section 24.3, printed 430–442 (PDF pages 445–457), especially Lemmas 24.3.1, 24.3.4–24.3.5, Definitions 24.3.6–24.3.7, Proposition 24.3.8, Theorem 24.3.9, Propositions 24.3.12–24.3.13, Corollary 24.3.14 and Lemma 24.3.15. The source is comparison and credit; the complete receiving proofs below are programme content. This geometry is the prerequisite for diffractive and general glancing propagation. It does not itself assert that a solution's wavefront follows every generalized trajectory or that every infinite-contact trajectory is a limit of transverse broken rays.
The measure foundation M0–M4 and M8 supplies countable covers, completed Lebesgue measure, monotone and dominated convergence, and compatibility with compact integrals. The complete substitution and null-image proofs P19.4 and P21.3–P21.4 supply the coordinate-volume facts. We prove below their precise extension to the open reflected clock domains, including the exceptional wall-hit slices. The matrix exponential proof A9 supplies the full convergent linear flow used in the comparison estimate. All selected locators and hypotheses are bound in the proof map.
1. Prepare a double normal root without changing the boundary
Use
ω=dα,ιHfω=−df,{f,g}=Hfg.(GF1)
Let a symplectic manifold have a smooth boundary defined by ϕ≥0, with dϕ=0. At a chosen point assume
p=0,Hpϕ=0,Hϕ2p=0.(GF2)
Here p is real and smooth. It need not have nonzero tangential differential. Prescribed-coordinate completion first makes r=ϕ one position coordinate. Its conjugate momentum is u, and all other coordinates, including r, are parameters z for the following scalar argument. Thus Hϕ=−∂u and (GF2) says p=pu=0, puu=0 at the mark.
Quadratic preparation lemma. A real smooth F(u,z) with those three properties can be written
F(u,z)=g(u,z)((u−b(z))2−S(z)),g=0,(GF3)
where b,S are real and smooth, with b=S=0 at the mark. The neighborhood is a fixed local neighborhood of the data.
Proof. Solve Fu(T(z),z)=0 near zero. The parameter inverse-function argument applies since Fuu=0; equivalently a uniformly contracting iteration for u−λ−1Fu(u,z), λ=Fuu(0,0), gives the root, and differentiation of its fixed-point identity gives every parameter derivative. Taylor's integral formula gives
F=F(T,z)+(u−T)2A(u,z), where A is smooth and has a fixed nonzero sign ε. Put v=(u−T)∣A∣ and s(z)=−εF(T,z). Then F=ε(v2−s), and u=β(v,z) is a smooth inverse near the mark, with βv=0.
We need the following exact smooth decomposition, not a formal root expansion:
β(v,z)=B0(v2,z)+vB1(v2,z),B1(0,0)=0,(GF4)
with both Bj smooth also at negative first arguments. To prove it, take the even part of β and the odd part divided by v. The latter is smooth by the integral formula for a first-order zero and is even. For any smooth even h(v,z), Taylor expansion to arbitrarily high even order shows that h(x,z), x≥0, is smooth up to x=0, with coefficients ∂v2jh(0,z)/(2j)!. This includes every parameter derivative: after any fixed number of x-derivatives, a Taylor remainder of sufficiently high order still tends to zero with the required power of x.
Extend these jets to both signs of x by the cutoff-series construction in the cited programme section. Explicitly, on a compact parameter patch the series ∑jcj(z)xjχ(x/δj) has those jets when χ=1 near zero and positive radii tend to zero fast enough that each differentiated term of order at most j/2, including parameter derivatives, is bounded by 2−j. Each fixed derivative series converges uniformly after finitely many terms; its jets follow by differentiating at zero. Subtract this smooth sum from h(x,z) for x≥0. The difference and all derivatives are flat at zero. Extending that difference by zero for x<0 is smooth, by the differentiated Taylor bounds. Adding it back gives the required exact extension. Apply this to both even functions. In particular B1(0,z)=βv(0,z).
Set b(z)=B0(s(z),z), S(z)=s(z)B1(s(z),z)2. Write Dj(v2,s,z) for the smooth divided difference of Bj, namely its first derivative integrated along the real segment from s to v2. Then
(β(v,z)−b(z))2−S(z)=(v2−s)[B1(s,z)2+2vB1(s,z)D+(v2−s)D2],D=D0+vD1.(GF5)
The bracket is nonzero near the mark since its value there is B1(0,0)2. Taking g to be ε divided by that bracket, expressed in (u,z), proves (GF3) exactly for positive, zero and negative s. No holomorphic extension or selection of nonsmooth square-root branches is used. □
Now prescribe r and ρ=u−b(z). They are a canonical pair because b is independent of u. Ordinary prescribed-coordinate completion retains them exactly. Also HrS=0, so S is independent of the new ρ. Dividing by the nonzero full function g therefore gives
p0(r,y;ρ,η)=ρ2−R(r,y,η),r≥0.(GF6)
This is a geometric multiplication and symplectic coordinate change; the symbol R is a general smooth real function. For an original quadratic principal symbol, the preceding normal coordinates retain its quadratic homogeneity directly. On the characteristic set, Hp=gHp0. A nonzero smooth factor changes the parameter consistently, with reversal if its sign is negative. All subsequent normalized statements transfer by that smooth nonvanishing change of parameter. We use the name p for p0 below.
Write z=(y,η), R0(z)=R(0,z), R1(z)=∂rR(0,z), and let HRz act in the tangential symplectic variables. Then
Hp=2ρ∂r+Rr∂ρ−HRz,Hpr=2ρ,Hp2r=2Rr.(GF7)
On p=0, identify the two boundary points with momenta ρ=±R0(z). Interior points are retained. Continuous compressed coordinates are (r,z,rρ); adjoining ρ2=R(r,z) gives a redundant useful continuous coordinate. The glancing set is G={r=ρ=R0=0}, and the hyperbolic boundary set has R0>0.
Define
Gk={p=0: Hpjϕ=0 (0≤j<k)},k≥2,G∞=k≥2⋂Gk.(GF8)
These are closed sets; G need not be a smooth hypersurface within the reduced space. At (0,z,0), for every integer m≥0,
(0,z,0)∈Gm+2⟺R0(z)=0,(−HR0)jR1(z)=0 (0≤j<m),Hpm+2r=2(−HR0)mR1there.(GF9)
Proof. Along the ordinary smooth Hamilton orbit from the point, r′′=2Rr(r,z), while z′=−HRz(r,z). Let zˉ′=−HR0(zˉ) have the same initial value. If the first m+1 jets of r vanish, the tangential equations imply that the first m+2 jets of z−zˉ vanish: differentiate the difference equation, using that its coefficients are smooth and its inhomogeneous difference is divisible by r. Thus the m-jet of Rr(r,z) equals that of R1(zˉ). Conversely, if the first m jets of the latter vanish, differentiating the two equations inductively from r(0)=r′(0)=0 gives vanishing of the first m+1 jets of r, and the next jet is twice the displayed derivative of R1. The induction starts with r′′(0)=2R1; this proves both directions and the formula. □
If ϕ is replaced by hϕ, h>0, and p by gp, g=0, the first nonzero contact derivative of order k becomes hgkHpkϕ. Indeed the restricted new Hamilton field is gHp, and every other differentiated product term contains a lower contact derivative. Thus order is intrinsic, and the sign at order two is unchanged. Set
Gd={r=ρ=R0=0, R1>0},Gg={r=ρ=R0=0, R1<0}.(GF10)
They are the diffractive and gliding parts of G∖G3. At Gd, the ordinary tangent orbit has positive r on both punctured sides. At Gg, its ordinary extension has negative r on both sides and is outside the manifold.
There is an invariant replacement for that outside motion. The constraint space C={ϕ=Hpϕ=0} is symplectic near the mark: the bracket of its two constraints is −Hϕ2p=0, so their Hamilton vectors form a nondegenerate two-plane complementary to TC. Its Hamilton field for p∣C, on G=C∩{p=0}, is
HpG=Hp+Hϕ2pHp2ϕHϕ,HpG=−HR0in (GF6).(GF11)
Both constraint derivatives vanish for this vector: HϕHpϕ=−Hϕ2p. For any vector tangent to C, its contraction with ω is −dp, since the added Hamilton field pairs to −dϕ=0. This proves the restricted Hamilton interpretation and tangency to the zero set even when that set is singular.
At a glancing point, replacing (ϕ,p) by (hϕ,gp) changes the numerator in (GF11) by g2h, the denominator by h2g, and Hϕ by h; hence
HgpG=gHpGon G.(GF12)
This proves independence from the defining function and consistent parameter change. In G3, R1=0, so HpG=Hp.
3. Define motion through reflections and tangencies
A generalized characteristic arc is a continuous compressed characteristic curve on an interval I. Its ordinary lift γ is defined outside a set B of reflection times. Each time in B is isolated within B; its two limits are the distinct hyperbolic roots with the same tangential data, and the neighboring points are interior. At every other time the lift is differentiable and obeys
γ′=Hpin the interior and at Gd,γ′=HpGat G∖Gd.(GF13)
The set B may accumulate at a glancing time not belonging to B. The compressed curve is the primary continuous object. The zero-length normal momentum at glancing has a unique lift; at a hyperbolic reflection either one-sided lift is specified instead of a single normal momentum.
For every such arc, the tangential variables are continuously differentiable, with
z′=−HRz(r,z).(GF14)
Reflection changes only ρ; at r=0 the right side is exactly the gliding field. The normal coordinate is continuous and locally Lipschitz, with one-sided derivatives 2ρ(t±0) at reflections. The product rρ has matching one-sided derivatives 2ρ2+rRr there. Consequently, on a fixed compact coordinate set,
r, z, rρ, ρ2=R(r,z)have uniform Lipschitz bounds,∣ρ(t)∣−∣ρ(s)∣=O(∣t−s∣1/2).(GF15)
For the last assertion use ∣a−b∣≤∣a−b∣ for a,b≥0. The bounds on the other variables follow from their bounded derivatives and one-sided derivatives. A continuous scalar function whose upper one-sided difference quotients are bounded by M, with the corresponding lower bounds, is M-Lipschitz: subtract (M+ϵ)t, and a first positive increase would give a nonnegative right difference quotient at its last running minimum, a contradiction. Let ϵ tend to zero. This elementary argument applies also across accumulating isolated reflections. The tangential right side is continuous in (r,z), so its integral equation proves (GF14) through those accumulation times.
These statements concern compact normalized coordinates. In a homogeneous characteristic cone choose a compact positive-cosphere section before using a uniform bound. No uniform bound on unnormalized covectors is being asserted.
4. Strict diffractive and gliding estimates
Suppose first that Rr≥c>0 in the coordinate region. A reflected orbit at time zero has outgoing ρ(+0)≥0, incoming ρ(−0)≤0. Integrating ρ′=Rr and r′=2ρ gives
r(t)≥ct2(t=0)(GF16)
as long as the orbit stays in the region. It cannot have another boundary contact there: ρ increases after the reflection in positive time and decreases on the backward orbit. Ordinary tangent orbits starting at Gd satisfy the same bound on both sides. Nearby reflected orbits with their initial normal roots tending to zero therefore converge to that ordinary tangent orbit, by smooth flow dependence on each side and continuity of the vanishing jump. This supplies the local motion and its uniqueness at Gd.
Suppose instead −C1≤Rr≤−c<0. Put
E=ρ2−rRr,m(ρ2+r)≤E≤M(ρ2+r),m=min(1,c),M=max(1,C1).(GF17)
On an ordinary segment,
E′=−r(dRr/dt); the two terms 2ρRr cancel exactly. Derivatives of Rr along Hp are bounded on the compact region. Hence ∣E′∣≤CE. At a reflection r=0 and ρ only changes sign, so E is continuous. Integration on the segments gives, in either direction,
ρ(t)2+r(t)≤mMeC∣t−s∣(ρ(s)2+r(s)).(GF18)
At a gliding point both normal quantities vanish. If a generalized arc left the boundary arbitrarily near that time, an interior component of {r>0}, together with its adjoining hyperbolic reflections, would obey (GF18). Taking a preceding or succeeding limiting normal state zero in that connected component forces its energy to vanish, a contradiction. Equivalently extend E by zero on the closed glancing-time set and apply the same integral inequality on its complementary intervals. Thus the arc stays at r=ρ=0 locally and equals the unique smooth orbit of −HR0.
For a nearby broken arc with r(0)≤ϵ2, ∣ρ(0)∣≤ϵ, (GF18) gives uniform r=O(ϵ2), ∣ρ∣=O(ϵ) on a fixed shorter interval. The tangential field differs from its boundary value by O(ϵ2). Comparing with the gliding orbit having the same tangential initial data and integrating its Lipschitz difference equation gives
z(t)−zˉ(t)=O(ϵ2),z′(t)−zˉ′(t)=O(ϵ2).(GF19)
This is uniform tangential C1 convergence, not convergence of the discontinuous full normal momentum in a differentiable topology. If tangential initial data also vary, add their initial difference with the usual bounded exponential factor.
The two panels are canonical models p=ρ2−r+η and p=ρ2+r+η, with the tangential coordinate equal to Hamilton time. The plotted radii are the exact formulas in Exercise 1. Blue reflected curves have c>0; the red curve is their c=0 limit. Normal momentum and its sign reversal are stated in the solution. The lower panel contains repeated parabolic excursions of height c2 and time spacing 2c; its limit is boundary gliding. These are coordinate sections of characteristic curves, not a wave amplitude.
For a generalized arc starting at a glancing point, let zˉ′=−HR0(zˉ), zˉ(0)=z(0), and set
e(t)=R1(zˉ(t)),f(t)=∣z(t)−zˉ(t)∣.(GF20)
On a compact smaller region, smoothness and (GF14) give upper one-sided inequalities
r′≤2∣ρ∣,(∣ρ∣)′≤∣e(t)∣+C2f+C3r,f′≤C0f+C1r.(GF21)
At reflections ∣ρ∣ is continuous. At gliding times its derivative is zero; elsewhere the displayed bound follows from the normal Hamilton equation. At higher glancing times the defining derivative is zero. The bounds therefore hold in the integrated one-sided sense across all times.
For clarity, the required comparison is fully scalar and finite-dimensional. Let (X,V,F) solve
X′=2V,V′=∣e∣+C2F+C3X,F′=C0F+C1X.(GF22)
Begin with initial values slightly larger than (r(0),∣ρ(0)∣,0), and add a positive constant to each right side. Up to a hypothetical first crossing, each difference has nonnegative derivative apart from its own nonnegative initial margin; at the first zero the strict added term gives a positive inward derivative, preventing the crossing. This is the standard first-contact argument, now with every coefficient nonnegative. Let the margins and added terms tend to zero by the linear integral equation. It follows that (r,∣ρ∣,f)≤(X,V,F).
The constant matrix in (GF22) has nonnegative entries. Its exponential has nonnegative entries by its power series. On 0≤t≤T, its entry carrying V to X is bounded by CteAt; the entry carrying V to F is bounded by Ct2eAt, since that path first goes through X; the entry carrying X to F is bounded by CteAt. These bounds follow directly by grouping the first nonzero matrix powers and bounding the remaining series by eAt. Variation of constants thus proves
r(t)∣ρ(t)∣f(t)≤CeAt[r(0)+t∣ρ(0)∣+∫0t(t−s)∣e(s)∣ds],≤CeAt[r(0)+∣ρ(0)∣+∫0t∣e(s)∣ds],≤CeAt[tr(0)+t2∣ρ(0)∣+∫0t(t−s)2∣e(s)∣ds].(GF23)
Apply the same argument after reversing the time parameter for negative time. If the initial point belongs to Gk, k≥3, (GF9) gives e(t)=O(∣t∣k−2), hence
r(t)=O(∣t∣k),∣ρ(t)∣=O(∣t∣k−1),z(t)−zˉ(t)=O(∣t∣k+1).(GF24)
The estimates retain the tangential drift and are valid through reflection accumulation. If e vanishes identically, all three differences vanish. In particular a point of G3 with HpG=0 can only have the constant generalized continuation. The same conclusion at Gg follows from (GF18) and smooth gliding-flow uniqueness.
Monotonicity lemma. Start at G3. If e is nondecreasing for small positive time, the generalized arc is the ordinary Hamilton orbit. If e is nonincreasing, it is the gliding orbit.
Proof for increasing e. Its initial value is zero, so e≥0. Formula (GF23) implies r=O(t2e(t)), ∣ρ∣=O(te(t)), f=O(t3e(t)), and consequently
Rr(r(t),z(t))=e(t)+O(t2e(t)).(GF25)
It is nonnegative for sufficiently small t. On every ordinary segment ρ′≥0, while any hyperbolic reflection would make a positive jump in ρ. Since the limiting initial value is zero, ρ cannot be negative. Such a reflection would require a negative incoming root, so none is possible. Any remaining boundary interval has Rr=0, where HpG=Hp. The complete arc therefore satisfies the ordinary smooth Hamilton equation, including boundary times, and equals its unique solution.
Proof for decreasing e. An initial interval on which e=0 has zero normal motion by (GF23); restart at its right endpoint if necessary. Where e<0, put h=ρ2−rRr. Formula (GF25), with ∣e∣, shows that Rr=e+O(t2∣e∣)<0, h≥0, and h/e2=O(t2). Comparing the actual and boundary tangential equations, as in (GF14), gives
dRr(r,z)/dt=e′(t)+O(t∣e(t)∣) on ordinary pieces: the normal velocity costs O(t∣e∣), the tangential differences O(t2∣e∣), and the position differences O(t2∣e∣). Thus, using e′/e≥0, r≤h/(−Rr), and h′=−r(dRr/dt),
h′≤(2ee′+Ct)h,(e2h)′≤Cte2h.(GF26)
Indeed the coefficient before h is at most (e′/e)(1+O(t2))+Ct, which is bounded by the displayed coefficient after shrinking. At reflections h is continuous. Extend it by zero on the closed glancing-time set; the same integrated inequality holds on its complementary intervals. Starting at a positive s and letting s↓0, the estimate h(s)/e(s)2=O(s2) forces h=0. If the negative interval begins after an initial zero interval, use elapsed time from that endpoint in the identical argument. Hence r=ρ=0, and the tangential motion is the gliding orbit. □
At a point of exact finite contact k, meaning Gk∖Gk+1, (GF9) gives
e(t)=atk−2+O(tk−1),a=2(k−2)!Hpkr(γ(0))=0.(GF27)
On either punctured side it is monotone after shrinking, with a fixed sign there. Applying the lemma in the forward and reversed directions proves that the motion on that side is the ordinary Hamilton orbit when the corresponding gliding orbit lies in Gd, and the gliding orbit when it lies in Gg. The ordinary orbit then has positive r, with the leading coefficient prescribed by (GF9); in the gliding case it would have negative r. These two choices do not introduce another branch.
Finite-contact uniqueness theorem. A generalized arc meeting no point of G∞ is uniquely determined by one compressed point, with either one-sided root specified at a hyperbolic reflection. Interior smooth-flow uniqueness, transverse reflection, the strict estimates and the finite-contact lemma prove local uniqueness at every point. Two arcs that agree at one time have a relatively open set of agreement; at any endpoint of that set, compressed continuity and the appropriate local uniqueness extend it. Hence they agree on their whole common interval. A point of G∖Gd with zero gliding field gives a constant arc by the preceding comparison; this statement also covers stationary infinite contact.
7. Limits of characteristic arcs remain characteristic
Compactness theorem. Suppose generalized arcs are defined on a common compact interval and stay in a fixed compact normalized coordinate set. Some subsequence of their compressed curves converges uniformly. Every such compressed limit has a unique generalized lift away from its hyperbolic reflection times, and ordinary lifts converge uniformly on compact time sets containing no such reflection.
First choose a subsequence for which the bounded Lipschitz coordinates in (GF15) converge. Here is the compactness argument explicitly. At a countable dense set of times, successive finite-dimensional bounded subsequences have convergent values; take the diagonal subsequence. A finite time grid of mesh less than ϵ/(3M) then makes convergence uniform by the common Lipschitz constant M. The resulting continuous coordinate limits satisfy ρ2=R(r,z), and rρ determines the signed root wherever r>0.
At an interior time the full states converge locally, and the smooth Hamilton integral equation passes to the limit. At a boundary point with positive limiting ρ2, the roots stay uniformly separated. The local ordinary flow has a unique simple boundary intersection by r′=2ρ=0; its implicit hitting time and the reflected outgoing data depend smoothly on the incoming data. A common smaller time neighborhood therefore contains exactly one reflection for all sufficiently late curves. Their hitting times converge, their two roots have the prescribed limits, and convergence is uniform away from that limiting reflection. Such reflection times form a discrete set locally where R0>0, hence at most a countable set globally.
At a point of Gd, the strict diffractive estimate gives the ordinary tangent limit on a common smaller neighborhood. At a point of Gg, the gliding energy estimate gives the unique gliding limit and tangential equation. At a point of G3, the comparison estimate applied to nearby initial data and then passed to the limit gives r=O(∣t−t0∣3), ∣ρ∣=O(∣t−t0∣2). Thus the lifted normal coordinates are differentiable there with derivative zero; the tangential integral equation gives derivative −HR0. This proves exactly (GF13), including all possible accumulation times, rather than merely containment in the characteristic set.
For uniform ordinary convergence away from limiting hyperbolic reflection times, cover the closed glancing-time set by a neighborhood where the limiting normal magnitude is arbitrarily small. Uniform convergence of ρ2 makes both signed roots uniformly small there; the tangential coordinates already converge uniformly. The complement is a compact union covered by finitely many interior or root-separated reflected flow charts; exclude the reflection times and use the already proved convergence on their compact subintervals. This proves the asserted uniform convergence without choosing a nonexistent normal sign at a reflection itself.
8. Almost every initial state has a finite reflected continuation
Assume in a convex coordinate neighborhood that one tangential position yd satisfies
∂ηdp≥c0>0.(GF28)
It is a strictly increasing clock. On a fixed clock section yd=s, solve p=0 for ηd; the implicit derivative is nonzero. The reduced coordinates are (r,y1,…,yd−1;ρ,η1,…,ηd−1), and their symplectic form is
dρ∧dr+j<d∑dηj∧dyj.(GF29)
Indeed dyd=0 on the section. A variable-time Hamilton map between two transverse sections of p=0 preserves the restricted form: differentiating the variable time adds multiples of Hp, whose contraction with each energy-tangent vector is −dp=0, while a fixed-time Hamilton map preserves ω. At a transverse wall section the form is ∑jdηj∧dyj; reflection is the identity in these variables. Composing incoming and outgoing section maps therefore proves that every finite transverse reflected clock map with both endpoints in the interior is smooth and symplectic. Hitting times are smooth by their nonzero normal derivative, so this interior-endpoint domain is open. The map is injective by reversing the same ordinary flow and reflection rules. Its Jacobian preserves Liouville volume, since it preserves the top exterior power of (GF29).
Clock slices that land exactly on a wall. The interior-endpoint qualification matters. An initial state with r=0 belongs to a coordinate hyperplane, hence to a null set in the reduced clock section. For a fixed target clock, consider an interior initial state whose finite transverse path has its last reflection exactly at that target. Keep its finitely many earlier hits fixed as smooth hitting-time branches. Extend the last incoming ordinary segment a little through r=0, using a smooth extension of the coefficients. The resulting map to the target clock is a local symplectic diffeomorphism in the ordinary reduced coordinates, by the preceding variable-time calculation. Therefore the initial states with target radius zero form locally a smooth hypersurface: they are the inverse image of the regular coordinate r=0. The implicit coordinate theorem identifies that hypersurface with a smooth image of a coordinate hyperplane. The complete C1 null-image theorem consequently proves its nullity.
These local neighborhoods have a countable subcover. Indeed the rational coordinate boxes form a countable base; for each basis box contained in one of the neighborhoods choose one such neighborhood. Every point is in a basis box so chosen, and each selected hypersurface piece is null. This covers every finite transverse itinerary landing on that fixed clock wall, without a bound on the number of earlier reflections. Countable subadditivity proves that their union is null. A finite or countable family of clock grids has the same property. Outside these null sets the initial and grid states used in the volume argument are interior, so the smooth clock maps just proved apply. Boundary initial states will still be recovered by limits in the existence argument; no boundary state is removed from its conclusion.
Measure transfer on the entire reflected domain. Here is the precise extension of the compact substitution theorem used below. If F:U→V is a smooth diffeomorphism of open coordinate sets with ∣detDF∣=1, then every open O⊂V satisfies
∣F−1(O)∣=∣O∣.(GFA1)
To prove it, enumerate rational boxes whose closures lie in O and whose interiors cover O. Each compact closure has a smooth function 0≤ψj≤1, supported in O and equal to one on that closure, by the complete finite cutoff construction. Then
fk=1−j=1∏k(1−ψj),0≤fk↑1O.(GFA2)
Each fk is compactly supported in O. Its pullback has compact support in U, since F−1 is continuous. Compact chart substitution and M8 give ∫Ufk(F(x))dx=∫Vfk(y)dy for the same Lebesgue integral. Monotone convergence M3 proves (GFA1), including infinite measures. This argument does not substitute an unproved general measurable change-of-variables theorem.
For a fixed clock pair, combine all interior-endpoint finite transverse paths into their open domain U. The map is injective by reversed reflected uniqueness and is locally a symplectic diffeomorphism. Its image is open, and its inverse is smooth on the inverse charts, so it is exactly a diffeomorphism of the kind in (GFA1), even if U has infinitely many components. Enlarge the normal strip below slightly to an open strip with the same O(ϵ3) volume bound, and intersect it with this open image. Formula (GFA1) bounds the measure of the entire inverse image by that one strip volume. Equivalently the countably many disjoint open components have disjoint images, and their measures add. The null wall-hit sets proved above cover the omitted clock endpoints. Thus no factor for the number of reflection itineraries appears.
Take a small compact initial clock box and a slightly larger coordinate box containing all short arcs under consideration. Bound the Hamilton coordinate velocities there. Tangential coordinates and r stay in the larger box for a common short time; ρ2=R(r,z) bounds the normal magnitude there. A reflected continuation cannot stop before the chosen short target clock except by approaching G. Away from ∣ρ∣=0, each reflection requires a uniform positive time before the normal derivative can reverse sign again, since ∣ρ′∣ is bounded. Thus infinitely many reflections in a finite interval force their normal magnitudes to tend to zero. The compressed Lipschitz bounds then give a limiting glancing state.
The reduced volume of the normal strip
0≤r≤Cϵ2,∣ρ∣≤Cϵ,with the other reduced coordinates bounded,is O(ϵ3).(GF30)
Suppose a continuation approaches a glancing point before its target clock. Use clock grid spacing ϵ, and choose the last grid section preceding a point with arbitrarily small r and ∣ρ∣. Clock and Hamilton elapsed times are comparable by (GF28) and the upper clock bound. On the finite reflected portion between these sections, ∣ρ∣ changes by at most Cϵ: its absolute value is continuous through each reflection and has bounded derivative on each segment. Integrating ∣r′∣≤2∣ρ∣ then gives r≤Cϵ2 at that grid section. This reasoning uses only the already existing reflected portion, not a presumed continuation beyond glancing.
The inverse image of that strip has the same volume by the symplectic clock map. The domain may have many reflection itineraries, but it is open and the combined map is injective; its countably many coordinate components have disjoint images, so volume preservation applies to their union. There are O(ϵ−1) grid sections, including the initial section if needed. Therefore all initial states whose short continuation fails are contained in sets of measure
O(ϵ−1)O(ϵ3)=O(ϵ2)for every sufficiently small ϵ.(GF31)
Their outer measure is zero. Apply the same argument backward. Almost every initial reduced state has a finite transverse broken arc over the whole common short clock interval. Its clock map is smooth symplectic on an open full-measure domain. This statement concerns that domain; it does not extend a single smooth map through infinite-contact nonuniqueness.
The volume diagram uses the exact coordinates x=r/ϵ2, y=ρ/ϵ and z=κ/ϵ3=xy, with κ=rρ, for the unit-constant strip in (GF30). The rectangle 0<x<1, ∣y∣<1 maps to the triangle 0<x<1, ∣z∣<x. Their physical Euclidean areas are 2ϵ3 and ϵ5, respectively. The preserved measure on the triangle has density drdκ/r and mass 2ϵ3, as proved in Exercise 4 and (GFA1). The normal boundary segment collapses to the vertex. Axes are separately normalized; other reduced coordinates are suppressed. The diagram illustrates the measure calculation, not the trajectories or a smooth extension at r=0.
Local existence at every point. The full-measure states are dense even at the boundary of the reduced half-space: every relative neighborhood has a positive-volume interior part. Approximate any initial state by them, parameterize each arc by Hamilton time from that initial clock, and use a common smaller Hamilton interval supplied by the two clock bounds. Compactness from Section 7 gives a generalized arc through the chosen compressed state, with the indicated one-sided root at a hyperbolic reflection. At glancing, the lift is unique at the initial time. If the tangential gliding field is nonzero, a fixed linear tangential symplectic change makes (GF28) hold; this preserves the normal form and the boundary. If that field vanishes, use the constant gliding arc at Gg or G3, and the ordinary tangent orbit at Gd. Interior and hyperbolic points already have their local ordinary or reflected continuations. Thus no nonzero tangential-differential hypothesis is needed for local existence.
Maximal continuations have the usual finite-time escape property. If one stayed in a compact coordinate region up to a finite endpoint, its compressed Lipschitz coordinates would have a limit. The local existence just proved supplies a continuation from that state, with the incoming root at a hyperbolic endpoint, and one joins it to the original arc. At a glancing join the strict or higher-contact estimates give the required derivative at the joining time. This contradicts maximality. Consequently a finite maximal endpoint leaves every compact subset of the phase domain. Under absence of G∞, local uniqueness gives a genuine flow. In clock neighborhoods its arcs are limits of transverse broken arcs by the same dense approximation and uniqueness. This approximation statement requires the clock hypothesis: for the permitted smooth symbol p=ρ2+r, the characteristic set in r≥0 has only r=ρ=0. Its constant gliding arcs exist and are unique, but there are no interior characteristic rays to approximate them. Thus stationary gliding existence does not imply transverse approximation. We also make no universal approximation assertion for all nonunique infinite-contact continuations.
9. Construct a trajectory inside a closed characteristic set
The next geometric theorem is what later wavefront arguments will need. Work in a normalized clock neighborhood satisfying (GF28), and let F be a relatively closed subset of its characteristic set. Put D=Gg∪G3, a relatively closed set because its complement within G is Gd. Assume:
- At each point of F∖D, the short ordinary, reflected or diffractive characteristic segment belongs to F.
- For every compact K⊂F∩D and ϵ>0, there is δ>0 and, for every q∈K, points qt∈F, ∣t∣≤δ, with q0=q and
tQ(qt)−Q(q)−Q(HpG(q))<ϵ,0<∣t∣≤δ,(GF32)
where Q=(r,z) drops just the normal momentum, and Q(HpG(q)) means those vector components. No continuity or differentiability of the choices qt is assumed.
Closed-set trajectory theorem. Through every q∈F there is a generalized arc entirely in F, maximally continued in both directions. In this clock neighborhood it leaves every compact subset of F at each maximal endpoint.
Proof. Fix a compact working region and choose the uniform δ in (GF32) for a slightly larger compact subset of F∩D. Decrease it so that δ≤ϵ. From a point outside D, follow its ordinary/reflected/diffractive continuation in F until it leaves the working region or has a limiting point in D. That limit exists by the compressed coordinate bounds; a finite terminal point outside D would continue by the first assumption. At a point in D, take the supplied point qδ, and interpolate linearly in the full coordinate space between the two endpoints for time δ. Both endpoints are in F; the interpolating segment is not asserted to lie in F.
Here are the precise error estimates making that interpolation legitimate. Its reduced slope differs from the gliding vector by at most ϵ. Its normal position is O(ϵδ), since the anchor has r=0 and the gliding normal component is zero. At the endpoint, smooth Taylor expansion of R, using R=0 and its zero derivative along HpG at the anchor, gives ρ2=O(ϵδ+δ2). Hence ∣ρ∣=O(ϵ) throughout the interpolated segment, as δ≤ϵ. The residual p is O(ϵ2), its distance from F is O(ϵ), and its equations satisfy
z′=−HRz(r,z)+O(ϵ),r′=2ρ+O(ϵ).(GF33)
Its rρ derivative is bounded, because r=O(ϵδ), ρ=O(ϵ), and the linear normal slope is O(ϵ/δ). These bounds hold uniformly over the compact anchor set.
Continue from the endpoint, using exact arcs outside D and another such interpolation upon reaching D. The clock increases at a uniformly positive rate: on exact pieces by (GF28), on interpolations by (GF32) with ϵ<c0/2. Each interpolation has the same positive duration δ; therefore only finitely many are needed before leaving the working region in a uniformly bounded time. Exact pieces may themselves have infinitely many reflections accumulating at their next anchor, but their compressed bounds and limiting anchor remain valid. The terminal interpolation can end just outside the working compact set; its endpoint belongs to F and is O(ϵ) from that set's boundary. Perform the construction also backward.
Let ϵ↓0. The reduced coordinates (r,z) have uniform Lipschitz bounds, and ρ2−R(r,z)=O(ϵ2). The product rρ has uniform bounds as just proved. A diagonal time-grid subsequence therefore converges uniformly in the compressed coordinates on the common bounded intervals. Every limit lies in F, since each interpolating piece has distance O(ϵ) from F, and the exact pieces lie there. Away from D, there are eventually no interpolations, so ordinary/reflected/diffractive flow charts give the required limit equations.
In a strict gliding region, the normal energy E of Section 4 remains O(ϵ2) on these approximations if they start near a limiting gliding state. Each interpolation starts at a zero normal state and ends with normal energy O(ϵ2); exact pieces multiply that bound by at most the fixed exponential in (GF18). There is no sum of errors at successive interpolations, because each new anchor has zero normal energy. Thus the limiting normal state is zero on a common smaller interval, and the limiting tangential equation from (GF33) is exactly the gliding equation.
It remains to justify the derivative at a limiting point in G3; this cannot be inferred just from characteristic containment. The tangential integral equations pass to the limit everywhere. Outside the countable hyperbolic reflection times, the signed normal momentum also converges pointwise: use rρ/r where r>0, and ρ2→0 at glancing. It is uniformly bounded. The normal integral equation in (GF33) therefore passes by dominated convergence, giving r′=2ρ in the integral sense. On a closed set where r=0, its integral derivative is zero; on the open complementary intervals it is the ordinary normal derivative. This gluing identity can also be checked directly by summing the integrals over the disjoint open intervals: their absolute integrals are bounded by the common derivative bound times the total interval length.
For completeness, the integrals used in this last step have their proved domains. Each approximating signed normal component is piecewise continuous away from its at most countable transverse hits and continuous at glancing accumulation, because its squared magnitude tends to zero there. It is therefore measurable: give it either value at the countable hits and use continuity on the remaining relative set. The limit is measurable by M3, is bounded on the common finite interval, and converges pointwise off the countable limiting hyperbolic hits. Each singleton is null by M2. Thus M3 applies with the actual constant integrable majorant on that finite interval. The resulting normal integral equation makes r Lipschitz and absolutely continuous directly: sums of increments on disjoint intervals are bounded by the uniform integrand bound times their total length. At any glancing point the signed limit ρ is continuous, since ρ2=R(r,z)→0; at an interior point it is smooth. The normal integral equation hence gives r′=2ρ at every nonreflection point by the ordinary continuous-integrand difference quotient. The tangential integrand is continuous everywhere. These facts justify the derivatives and countable exceptional set in the next calculation.
Consequently R(r,z) is absolutely continuous and
R′=Rrr′=2Rrρ almost everywhere. The tangential terms cancel because z′=−HRz. The smooth-composition chain rule here follows by Taylor expansion on a time partition: the quadratic errors tend to zero since the maximum increment tends to zero while the total variation stays bounded, and the first-order sums tend to the integral of the displayed derivatives. At a G3 point, taken as time zero, Rr(0)=0 and the Lipschitz reduced coordinates give ∣Rr(t)∣≤C∣t∣. Regularize the square root to obtain
dtdR(t)+a≤C∣t∣,a>0,R(t)=O(t2)(a↓0).(GF34)
Thus ∣ρ(t)∣=O(t2); the signed normal momentum has derivative zero at the point, and r(t)=O(∣t∣3). Together with the tangential equation this is precisely the generalized derivative HpG=Hp there. The limit is a generalized arc, not merely an approximate orbit.
Its endpoints are in F on the boundary of the working compact region. Extending with increasing compact neighborhoods gives a maximal arc in F. A finite endpoint inside a compact subset of F would have a limiting state and the same compact-region construction would extend it. An infinite time endpoint cannot stay in a fixed compact subset either, because the clock has a positive lower rate there. This proves the stated escape property. □
The theorem is geometric. To apply it to WFb(u), the differential equation must supply both hypotheses, including the quantitative gliding tangency condition (GF32). The transverse reflection theorem alone does not supply that condition. Those analytic propagation proofs are subsequent work.
Here is an explicit homogeneous real quadratic example. It also shows why the finite-contact qualification cannot be removed just by saying that the coefficients are smooth.
Let tj=2−j, vj=e−j2, j≥1, and for j≥2 set
Δj=tj−1−tj=2−j,μj=vj+vj−1Δjvj.(GF35)
Choose a nonnegative even smooth bump κ supported in [−1/2,1/2], with integral one. Thus its first moment is zero. Put
ψj(t)=μj−1κ(μjt−tj−1+μj),a(t)=−21j≥2∑(vj+vj−1)ψj(t),a(t)=0 (t≤0).(GF36)
The supports are disjoint and contained in (tj,tj−1), since 3μj/2<Δj. For every derivative order m,
∥a(m)∥(tj,tj−1)≤Cm(vj+vj−1)μj−m−1≤Cmexp(−j2+Cm′j).(GF37)
The last inequality follows from vj/vj−1=e−2j+1 and Δj=2−j. It tends to zero faster than tjN for every fixed m,N. Extending every derivative by zero at t=0 is smooth: the difference quotient of the m-th derivative tends to zero by the bound with N=2, giving the next derivative inductively. Hence a is C∞, flat at zero and nonpositive everywhere. No finite differentiability approximation is used.
On each interval [tj,tj−1], define
r(t)=vj(t−tj)+2∫tjt(t−s)a(s)ds,ρ(t)=21r′(t).(GF38)
The normalization and first moment of ψj give exactly
r(tj)=r(tj−1)=0,r′(tj)=vj,r′(tj−1)=−vj−1,r′′=2a≤0.(GF39)
For the endpoint value, ∫(tj−1−s)ψj(s)ds=μj, and the definition of μj cancels vjΔj exactly. Concavity and the positive initial and negative final derivatives show r>0 between the endpoints. At tj the incoming momentum −vj/2 is reflected to +vj/2, so adjacent intervals match the exact reflection law.
Take the original cotangent coordinates (r,t,w;ρ,η,σ), r≥0, and the real degree-two principal symbol
p=ρ2−ra(t)σ2+ησ.(GF40)
It is the symbol of Dr2−ra(t)Dw2+DtDw, and the boundary conormal coefficient is one. Its full quadratic form is nondegenerate: the (η,σ) block has determinant −1/4. Along σ=1, set
η(t)=r(t)a(t)−ρ(t)2,w(t)=∫0t[−r(s)a(s)−ρ(s)2]ds.(GF41)
Then p=0, t′=1, r′=2ρ, ρ′=a, η′=ra′, w′=−2ra+η, σ′=0: the η equation follows by differentiating (GF41), with exact cancellation of 2ρa. At every reflection r=0, the momentum η=−ρ2 is continuous, as are w and all tangential data. These are actual Hamilton pieces of the same homogeneous symbol.
The bounds 0≤r≤Δjvj, ∣ρ∣≤vj−1/2, together with (GF37), show that the normal and tangential differences from the initial state tend to zero faster than every power of t at zero. Indeed ∣ra∣+ρ2≤Cmexp(−2j2+Cmj) on the intervals; its integral bounds w, and (GF41) bounds η. In particular the curve is differentiable there with the gliding derivative. Here R1=a(t)σ2, and the tangential gliding field at σ=1,η=0 differentiates a in t; (GF9) and flatness prove that all ordinary normal contact derivatives at the initial point vanish. Its infinitely many transverse reflections accumulate at that G∞ point.
A second generalized continuation through the same full initial state is
r=ρ=η=w=0,σ=1,t=s.(GF42)
It is the gliding orbit: wherever a<0 the point is in Gg, and where a=0 it is in G3. Both continuations agree with this gliding orbit for t≤0, but the first has r>0 on every positive interval between reflections. Hence smoothness, real principal type and nondegenerate quadratic form do not imply uniqueness at infinite contact. This example proves the nonuniqueness phenomenon directly; it is not a claim that every nonunique generalized arc has or lacks a transverse-ray approximation.
The diagram displays the exact dyadic intervals and the defining normal-profile mechanism in (GF35)–(GF39). For its numerical profile samples choose κ(z)=Cexp[−1/(1−4z2)] on ∣z∣<1/2, zero outside, where C is the reciprocal of its integral. The plotted time and radius scalings are stated explicitly; the samples illustrate the smooth integral and do not replace the exact moment identities or flatness bounds. The complete homogeneous Hamilton lift, covectors and two continuations are (GF40)–(GF42).
11. Graded exercises with complete solutions
Level: intermediate. For pd=ρ2−r+η and pg=ρ2+r+η, use y as Hamilton time. Start a reflection at y=0,r=0, with normal magnitude c>0. Find the radii and normal momenta. Identify the c↓0 limits and verify the quantitative scales in the preceding estimates.
Solution. In both cases y′=1, and η=−c2. For the diffractive sign, ρ′=1, so the incoming and outgoing pieces are
r(y)=y2+2c∣y∣,ρ(y)=y+csgny(y=0).(GF43)
They solve pd=0, have the correct reflection limits, and tend uniformly on bounded intervals to the tangent ordinary orbit r=y2,ρ=y. At zero that orbit has Hp2r=2>0.
For the gliding sign, on the excursion after each reflection at y=2kc,
u=y−2kc∈[0,2c],r=2cu−u2,ρ=c−u.(GF44)
Now ρ′=−1, ρ2+r=c2, and the jump at each wall exchanges −c for +c. Thus r≤c2, ∣ρ∣≤c, and ∣rρ∣≤c3. The tangential momentum is constant and y′=1. The limit is boundary gliding r=ρ=η=0, not the ordinary tangent orbit r=−y2, which is outside. This checks the signs, the normal scales and the tangential convergence; it does not claim convergence of the normal momentum derivatives across its jumps.
Level: advanced. Let k≥3, ϵ∈{1,−1}, and p=ρ2−rϵyk−2+η. Construct the generalized arc through (r,y;ρ,η)=(0,0;0,0), and distinguish even and odd k.
Solution. Here y′=1, ρ′=ϵyk−2, and the gliding arc is r=ρ=η=0. The ordinary orbit from zero is
r=k(k−1)2ϵyk,ρ=k−1ϵyk−1,η=k(k−1)2(k−2)y2k−2.(GF45)
Substitution gives p=0; differentiating the last expression gives η′=rϵ(k−2)yk−3, since ϵ2=1. Its first nonzero normal derivative is order k. On either side select this ordinary orbit if ϵyk−2>0, and the gliding orbit if that quantity is negative. These choices have matching derivatives at zero as required by (GF13). For even k, ϵ=1 gives ordinary interior motion on both sides, and ϵ=−1 gives gliding on both sides. For odd k, the two sides have opposite choices. Formula (GF45) has positive r exactly on its selected ordinary side. The finite-contact theorem proves uniqueness; an arbitrary extra choice of reflection at zero would violate the glancing definition.
Exercise 3. Compute an exact preparation multiplier
Level: advanced. Suppose the inverse Morse coordinate is u=β(v)=v+bv2, with constant real b, near zero, and F(u,z)=v(u)2−z. Produce the exact quadratic preparation in u, including its nonzero full multiplier. Explain the difference between preparation and merely matching two real roots.
Solution. The derivative of β at zero is one, so the inverse is smooth. The even/odd decomposition is B0(s)=bs, B1(s)=1. Therefore
Q(u,z)=(u−bz)2−z,Q(β(v),z)=(v2−z)[1+2bv+b2(v2−z)].(GF46)
The bracket has value one at zero and is nonzero nearby. Hence
F(u,z)=1+2bv(u)+b2(v(u)2−z)Q(u,z).(GF47)
This exact smooth identity holds for negative z, when no real roots exist, and at z=0, when the two roots coincide. A proof based only on matching roots for z>0 would not justify the smooth full multiplier through those other parameter values. If z consists of the parameters other than a normal momentum, shifting that momentum by bz is a prescribed canonical pair with the unchanged normal position; full canonical completion supplies the remaining coordinates.
Exercise 4. Keep the correct volume after compression
Level: intermediate. In a reduced clock section, compare the Liouville measure of 0<r<ϵ2, ∣ρ∣<ϵ with the Euclidean area after κ=rρ. Determine which measure belongs in the almost-everywhere continuation proof.
Solution. The ordinary reduced volume in these two coordinates is 2ϵ3. For r>0,
dρ∧dr=r−1dκ∧dr,∫0ϵ2∫−rϵrϵrdκdr=2ϵ3,∫0ϵ22rϵdr=ϵ5.(GF48)
The last quantity is Euclidean area in compressed coordinates, not Liouville measure. The clock maps preserve the symplectic volume with density 1/r there. Multiplying the correct strip bound O(ϵ3) by O(ϵ−1) clock sections gives O(ϵ2), as in (GF31). Compression is continuous at the boundary, but its coordinate Jacobian does not allow replacing the preserved measure by ordinary Euclidean compressed area.
Written by GPT-6.1 Sol (OpenAI); contributions by GPT-6 Astra (OpenAI), October 2026. Original exposition, examples, diagrams, the solved exercises are dedicated under CC0 1.0 Universal. Exact programme prerequisites retain their recorded licences.
Source and credit: Hörmander, The Analysis of Linear Partial Differential Operators III, Springer 2007 edition, ISBN 978-3-540-49938-1, Section 24.3, printed 430–442 (PDF pages 445–457). The quadratic preparation argument uses the existing programme's full parameter-aware cutoff-series and prescribed-coordinate constructions. The infinite-contact force profile is independently constructed here to demonstrate the known failure of uniqueness; the source gives a different convex-boundary polygon example. Diffractive/glancing wavefront propagation and singular distributions on arbitrary limiting rays require their respective analytic arguments.