A flat convex boundary with two characteristic continuations

A smooth convex obstacle can have infinitely many transverse reflections in a finite travel time. At their accumulation point the ordinary ray has contact of every order with the boundary. The reflected continuation and the boundary gliding continuation then pass through the same full characteristic state.

We construct the boundary, the polygon, its exact reflection law and both Hamilton lifts. The construction complements the variable-coefficient example in Tangency, gliding and generalized characteristic motion, Section 10. Its inputs are the contact-order calculation in that lesson's Section 2, its generalized-arc definition in Section 3, and ordinary smooth calculus and convergent series in the stationary-phase foundations. The exact earlier proofs are identified in the proof map.

Hörmander III, Springer 2007 edition, ISBN 978-3-540-49938-1, Example 24.3.11, printed 438, gives the convex-boundary nonuniqueness phenomenon. The formulas and interpolation below are independently chosen.

1. Prescribe chords and their reflecting tangents

D0. A polygon accumulating at the origin. For integers j≥1j\ge1, set

xj=1j,dj=xj−xj+1=1j(j+1),mj=2−j,yj=∑ℓ=j∞dℓmℓ.(CB1) x_j=\frac1j,\qquad d_j=x_j-x_{j+1}=\frac1{j(j+1)},\qquad m_j=2^{-j},\qquad y_j=\sum_{\ell=j}^{\infty}d_\ell m_\ell . \tag{CB1}
The series converges absolutely. Its tail satisfies 0<yj≤mjxj0<y_j\le m_jx_j. Hence Pj=(xj,yj)P_j=(x_j,y_j) tends to the origin, and the segment from Pj+1P_{j+1} to PjP_j has slope exactly mjm_j. In increasing xx, the polygon's slopes increase. Put also m0=1m_0=1, and define
αj=arctan⁡mj,βj=αj−1+αj2,bj=tan⁡βj.(CB2) \alpha_j=\arctan m_j,\qquad \beta_j=\frac{\alpha_{j-1}+\alpha_j}{2},\qquad b_j=\tan\beta_j . \tag{CB2}
If a convex smooth graph through PjP_j has derivative bjb_j there, its tangent bisects the incoming and outgoing chord directions. The equal-angle reflection law will therefore hold exactly, rather than just asymptotically.

D1. A uniformly interior interpolation parameter. Set

θj=mj−bj+1bj−bj+1,cj=1−θj.(CB3) \theta_j=\frac{m_j-b_{j+1}}{b_j-b_{j+1}},\qquad c_j=1-\theta_j . \tag{CB3}
The denominator is positive. In fact bj+1<mj<bjb_{j+1}<m_j<b_j, and 1/9≤θj≤2/31/9\le\theta_j\le2/3. Here are bounds independent of jj. Write m=mj≤1/2m=m_j\le1/2. On [0,1][0,1] the derivative of arctan⁡\arctan lies between 1/21/2 and 11; on the corresponding angles the derivative of tan⁡\tan lies between 11 and 22. Applying these bounds to the two half-angle differences gives
m4≤bj−m≤m,m8≤m−bj+1≤m2.(CB4) \frac m4\le b_j-m\le m,\qquad \frac m8\le m-b_{j+1}\le\frac m2 . \tag{CB4}
Taking the ratio of the second gap to the sum proves the assertion. In particular 1/3≤cj≤8/91/3\le c_j\le8/9. This uniform margin is what permits one smooth interpolation profile at every shrinking interval.

2. Construct the actual smooth convex graph

I0. An exact moment profile. Choose a nonnegative even smooth function κ\kappa, positive on (−a,a)(-a,a), supported in [−a,a][-a,a], with integral one, where a=1/32a=1/32. For this construction, take it to be the normalization of exp⁡[−1/(1−(u/a)2)]\exp[-1/(1-(u/a)^2)] for ∣u∣<a|u|<a, extended by zero. Every derivative at an endpoint is zero: differentiation produces a rational power of the endpoint distance times the same exponential, and v−Ne−1/v→0v^{-N}e^{-1/v}\to0. The integral is positive and finite. Define

K(v)=∫−∞vκ(u) du,L(v)=∫−∞vK(u) du.(CB5) K(v)=\int_{-\infty}^v\kappa(u)\,du,\qquad L(v)=\int_{-\infty}^v K(u)\,du . \tag{CB5}
For 0≤s≤10\le s\le1, put Hj(s)=K(s−cj)H_j(s)=K(s-c_j). It is zero near zero, one near one, and nondecreasing. Its exact mean is
∫01Hj(s) ds=∫−aa(1−cj−u)κ(u) du=1−cj=θj.(CB6) \int_0^1 H_j(s)\,ds =\int_{-a}^{a}(1-c_j-u)\kappa(u)\,du =1-c_j=\theta_j . \tag{CB6}
The exchange of integrals can equally be checked by integration by parts on the compact rectangle; evenness gives the zero first moment. All derivatives of HjH_j have bounds depending only on their order, since these are translates of one fixed profile.

I1. Match values, slopes and every higher derivative. On [xj+1,xj][x_{j+1},x_j] prescribe

g(x)=bj+1+(bj−bj+1)K(x−xj+1dj−cj),f(x)=yj+1+∫xj+1xg(v) dv.(CB7) \begin{split} g(x)&=b_{j+1}+(b_j-b_{j+1}) K\left(\frac{x-x_{j+1}}{d_j}-c_j\right),\\ f(x)&=y_{j+1}+\int_{x_{j+1}}^x g(v)\,dv . \end{split} \tag{CB7}
Formula (CB6) gives f(xj)−f(xj+1)=djmjf(x_j)-f(x_{j+1})=d_jm_j, so the endpoint values are exactly yj,yj+1y_j,y_{j+1}. Near either endpoint gg is the corresponding constant bjb_j or bj+1b_{j+1}. Neighboring intervals thus have identical first derivatives and zero derivatives of every order at least two at their shared endpoint. The resulting function is smooth on 0<x≤10<x\le1, has f′=gf'=g, and is convex since g′≥0g'\ge0. Extend it linearly past x1x_1 with slope b1b_1.

I2. Verify smoothness and flatness at the accumulation point. For every integer k≥0k\ge0,

sup⁡[xj+1,xj]∣g(k)∣≤Ck2−j[j(j+1)]k.(CB8) \sup_{[x_{j+1},x_j]}|g^{(k)}| \le C_k 2^{-j}[j(j+1)]^k . \tag{CB8}
For k=0k=0, use 0<bj+1≤g≤bj≤2mj0<b_{j+1}\le g\le b_j\le2m_j; for k≥1k\ge1, differentiate (CB7) and use the uniform profile bounds. The same bound for f(k+1)f^{(k+1)} tends to zero faster than every power of xx, since x≥1/(j+1)x\ge1/(j+1) and 2−jjM→02^{-j}j^M\to0 for each fixed MM. The latter limit follows, for example, from the ratio of successive terms tending to 1/21/2.

Set f(x)=0f(x)=0 for x≤0x\le0. Its value tends to zero, because 0≤f(x)≤yj0\le f(x)\le y_j on the jj-th interval. The displayed bounds show that each right derivative tends to zero and is o(x)o(x). Inductively the difference quotient at zero gives the next derivative, also zero; on the negative side all derivatives vanish. Thus ff is genuinely C∞C^\infty and flat at zero. Its derivative is nondecreasing on the whole line, so its convexity persists through zero. No finite-order smoothing or formal series stands in for the actual graph.

3. Check the full billiard and wave characteristics

R0. Chords lie strictly in the permitted interior. Let

X={(x,y,t):y≥f(x)},ϕ=y−f(x).(CB9) X=\{(x,y,t):y\ge f(x)\},\qquad \phi=y-f(x). \tag{CB9}
Convexity places each chord above the graph. The inequality is strict between its endpoints. Indeed the graph derivative starts below the chord slope mjm_j, ends above it, and is nondecreasing; the difference between chord and graph therefore first increases and then decreases, returning to zero only at the right endpoint. A nontrivial middle plateau cannot have zero difference, because the initial strict increase has already made it positive.

The unit velocity on chord jj is vj=(cos⁡αj,sin⁡αj)v_j=(\cos\alpha_j,\sin\alpha_j). At PjP_j, reflection in the tangent line of angle βj\beta_j takes direction αj\alpha_j to 2βj−αj=αj−12\beta_j-\alpha_j=\alpha_{j-1}. Equivalently it takes vjv_j to 2(vj⋅Tj)Tj−vj2(v_j\cdot T_j)T_j-v_j, where Tj=(cos⁡βj,sin⁡βj)T_j=(\cos\beta_j,\sin\beta_j). Thus it is the exact Euclidean specular reflection. The two normal components are nonzero because αj<βj<αj−1\alpha_j<\beta_j<\alpha_{j-1}.

R1. Use finite arclength as the Hamilton parameter. Let

ℓj=dj1+mj2,sj=∑ℓ=j∞ℓℓ,xj≤sj≤5/4 xj.(CB10) \ell_j=d_j\sqrt{1+m_j^2},\qquad s_j=\sum_{\ell=j}^{\infty}\ell_\ell,\qquad x_j\le s_j\le\sqrt{5/4}\,x_j . \tag{CB10}
The polygon from the origin has finite length to every vertex. On sj+1<s<sjs_{j+1}<s<s_j, take (x(s),y(s))=Pj+1+(s−sj+1)vj(x(s),y(s))=P_{j+1}+(s-s_{j+1})v_j, and set physical time t=st=s. For the real homogeneous quadratic symbol
q(x,y,t;ξ,η,τ)=12(ξ2+η2−τ2),(ξ,η,τ)=(vj,−1),(CB11) q(x,y,t;\xi,\eta,\tau) =\tfrac12(\xi^2+\eta^2-\tau^2),\qquad (\xi,\eta,\tau)=(v_j,-1), \tag{CB11}
we have q=0q=0, Hq=(ξ,η,−τ;0,0,0)H_q=(\xi,\eta,-\tau;0,0,0), and hence precisely the stated unit-speed Hamilton pieces.

At a reflection the tangential covector is unchanged: for the boundary tangent (1,f′(x),0)(1,f'(x),0), the two velocities have the same scalar product, while τ=−1\tau=-1 remains fixed. Their difference is a multiple of dϕ=(−f′,1,0)d\phi=(-f',1,0). Thus they are the two distinct characteristic roots in the same boundary fiber. Also Hϕ2q=1+f′2>0H_\phi^2q=1+f'^2>0, so the required double normal root is nondegenerate at every boundary point, and the boundary is noncharacteristic for the wave operator with symbol qq.

R2. Check the generalized derivative at infinite contact. On the jj-th chord,

∣vj−(1,0)∣≤C2−j,∣y(s)∣≤C2−jxj,∣x(s)−s∣≤C 4−jxj.(CB12) |v_j-(1,0)|\le C2^{-j},\quad |y(s)|\le C2^{-j}x_j,\quad |x(s)-s|\le C\,4^{-j}x_j . \tag{CB12}
The last estimate follows by summing 0≤ℓℓ−dℓ≤dℓmℓ2/20\le\ell_\ell-d_\ell\le d_\ell m_\ell^2/2, including the current partial chord; for that partial chord use 1−cos⁡αj≤mj2/21-\cos\alpha_j\le m_j^2/2. Since ss is comparable to 1/j1/j between successive vertices, all three errors are smaller than every fixed power of ss. These bounds include both one-sided velocities at each reflection.

For s≤0s\le0, continue by x=t=s,y=0,ξ=1,η=0,τ=−1x=t=s,y=0,\xi=1,\eta=0,\tau=-1. At zero the positions have derivative (1,0,1)(1,0,1) and the covectors derivative zero, including along the accumulating reflected pieces. The limiting full state is q0=(0,0,0;1,0,−1)q_0=(0,0,0;1,0,-1). Along its ordinary Hamilton orbit, ϕ(s)=−f(s)\phi(s)=-f(s). Flatness of ff proves Hqkϕ(q0)=0H_q^k\phi(q_0)=0 for every kk; thus q0∈G∞q_0\in G^\infty. The gliding field there equals HqH_q. The resulting reflected arc therefore satisfies the actual generalized-arc definition at its accumulation time. Its reflection times are isolated individually and accumulate only at zero, which is a glancing time, not a transverse reflection.

4. Construct the different continuation through the same state

G0. The boundary arclength lift is a genuine gliding solution. Parameterize the boundary graph by its signed arclength ss, with x(0)=0x(0)=0. The inverse is smooth because ds/dx=1+f′(x)2>0ds/dx=\sqrt{1+f'(x)^2}>0. Set

y=f(x),t=s,ξ=(1+f′2)−1/2,η=f′(1+f′2)−1/2,τ=−1.(CB13) y=f(x),\quad t=s,\quad \xi=(1+f'^2)^{-1/2},\quad \eta=f'(1+f'^2)^{-1/2},\quad \tau=-1. \tag{CB13}
It has q=0q=0 and Hqϕ=0H_q\phi=0. Direct calculation on these constraints gives
Hq2ϕ=−f′′ξ2,HqG=Hq−f′′ξ21+f′2Hϕ,(ξ′,η′)=f′′(1+f′2)2(−f′,1).(CB14) H_q^2\phi=-f''\xi^2,\qquad H_q^G=H_q-\frac{f''\xi^2}{1+f'^2}H_\phi,\qquad (\xi',\eta')=\frac{f''}{(1+f'^2)^2}(-f',1). \tag{CB14}
Differentiating (CB13) with x′=(1+f′2)−1/2x'=(1+f'^2)^{-1/2} gives exactly the last formula. Its base derivatives are (ξ,η,1)(\xi,\eta,1), so this is the full restricted Hamilton equation. Where f′′>0f''>0 the state is strictly gliding; where f′′=0f''=0 it is in G3G^3, and the same equation applies. This proves the gliding solution through every interpolation and flat segment.

At s=0s=0 this arc has the same full state q0q_0 as the reflected arc and agrees with it for s≤0s\le0. For every positive chord interval the reflected arc has y>f(x)y>f(x), whereas this one stays on the boundary. They are different continuations of the same smooth real constant-coefficient wave characteristic problem. Finite-contact uniqueness cannot be extended to all infinite-contact states.

5. Both specified continuations are limits of transverse broken rays

L0. Approximation of these two curves. For this example the reflected arc and the boundary arc are each uniform compressed limits of finite transverse broken rays on every compact physical-time interval in their stated domain. We do not infer this property for all nonunique generalized rays.

First check local uniqueness along the positive boundary arc, including its flat portions. Write σ\sigma for boundary arclength, T,NT,N for its tangent and inward unit normal, and

ϰ=f′′(1+f′2)3/2≥0,F(r,σ)=F(0,σ)+rN(σ).(CB15) \varkappa=\frac{f''}{(1+f'^2)^{3/2}}\ge0,\qquad F(r,\sigma)=F(0,\sigma)+rN(\sigma). \tag{CB15}
At r=0r=0, the two coordinate derivatives are N,TN,T, hence the map has a smooth local inverse. Since Nσ=−ϰTN_\sigma=-\varkappa T, its Euclidean metric is dr2+(1−rϰ)2dσ2dr^2+(1-r\varkappa)^2d\sigma^2. The actual pulled-back wave symbol, multiplied by two, is therefore
p=ρ2−R,R=τ2−ησ2(1−rϰ)2,R0=τ2−ησ2,R1=−2ϰησ2.(CB16) p=\rho^2-R,\qquad R=\tau^2-\frac{\eta_\sigma^2}{(1-r\varkappa)^2},\qquad R_0=\tau^2-\eta_\sigma^2,\qquad R_1=-2\varkappa\eta_\sigma^2 . \tag{CB16}
Its boundary gliding field has constant ησ,τ\eta_\sigma,\tau, with σ′=2ησ\sigma'=2\eta_\sigma and t′=−2τt'=-2\tau. This factor-two Hamilton parameter agrees with the physical-time parameter of (CB11) after the stated constant rescaling.

Where ϰ>0\varkappa>0, strict gliding uniqueness follows from the normal-energy estimate GGL:T016–T017. In an open flat segment, the higher-contact comparison has e=0e=0, so the continuation is also uniquely the boundary line there. At a nonzero endpoint of a curvature bump, the curvature is monotone just inside that bump. To verify the last assertion, on its interval (CB7) gives

f′′=bj−bj+1dj κ ⁣(x−xj+1dj−cj),ddxlog⁡ϰ=κ′(S−cj)djκ(S−cj)−3f′f′′1+f′2.(CB17) f''=\frac{b_j-b_{j+1}}{d_j}\, \kappa\!\left(\frac{x-x_{j+1}}{d_j}-c_j\right),\qquad \frac{d}{dx}\log\varkappa =\frac{\kappa'(S-c_j)}{d_j\kappa(S-c_j)} -\frac{3f'f''}{1+f'^2}. \tag{CB17}
For the fixed exponential bump in I0, its logarithmic derivative tends to positive infinity at the left edge and negative infinity at the right edge. The second term is bounded on this one interval. Thus ϰ\varkappa increases from zero at entry and decreases to zero at exit on sufficiently small one-sided neighborhoods. The map from xx to σ\sigma is strictly increasing. In the time direction entering either bump edge, e=R1=−2ϰe=R_1=-2\varkappa, on the unit gliding state, is nonincreasing from zero. GGL:T024 forces the gliding continuation; on the flat side e=0e=0 gives the same conclusion. Every positive boundary point therefore has local uniqueness. This uses the specific bump just chosen, rather than a claim for arbitrary flat curvature.

Every positive point of the polygon has local uniqueness as well, by ordinary flow or a transverse reflection. Neither curve has reflection accumulation at a positive time. On the negative side of zero the boundary is straight and R1=0R_1=0. GGL:T019–T022, applied backward from the common state at zero, forces the unique boundary-line continuation. Hence each complete curve is the only generalized continuation through any one of its own fixed positive states on the compact interval considered: local uniqueness propagates agreement through every positive time, continuity supplies zero, and the backward zero-forcing comparison supplies negative times. This conclusion does not say that the common state at zero alone determines the positive future.

It remains to obtain approximating finite rays through nearby states. Physical time is a uniform clock for these wave characteristics near τ=−1\tau=-1. On a bounded physical-time interval the base velocity has magnitude one, and reflection preserves that speed and the spatial covector magnitude. Nearby normalized states therefore remain in a common compact phase region. Apply the volume argument GGL:T031–T035 on this compact region. A finite boundary-chart cover has a common normal strip estimate O(ϵ3)O(\epsilon^3); the preceding clock grid has O(ϵ−1)O(\epsilon^{-1}) sections. The injective finite reflected clock maps preserve the whole reduced volume. The initial states whose continuation reaches glancing before either end thus have outer measure O(ϵ2)O(\epsilon^2) for every small ϵ\epsilon, and form a null set. Fixed-clock wall endpoints are the null exceptions treated in GGL:A1. The same proof works forward and backward, and across the finitely many compact charts.

Consequently finite transverse broken arcs on the entire interval have initial states dense near either selected positive state. Normalize them by their constant negative time covector and choose a sequence approaching that state. GGL:T027–T030 supplies uniformly convergent compressed subsequences, and every subsequential limit is a generalized arc through that state. The just-proved uniqueness on this interval identifies it with the chosen polygon or gliding arc. If the whole sequence did not converge, a subsequence with a fixed positive uniform error would have a further compactness limit with zero error, a contradiction. This proves the claimed uniform compressed approximation of both curves.

Convex chord interpolation and its exact equal-angle reflection

V0. What the diagram measures. The left panel shows the actual vertices PjP_j, straight chords and the smooth convex graph, with vertical coordinates multiplied by the stated factor. The middle panel subtracts the left tangent on the interval j=5j=5 and divides by dj(bj−bj+1)d_j(b_j-b_{j+1}): the chord is exactly θjS\theta_j S, the graph is L(S−cj)L(S-c_j), and their difference is positive. The right panel uses an orthonormal frame at P2P_2; its two unit directions have angles ±(α1−α2)/2\pm(\alpha_1-\alpha_2)/2 to the tangent. Numerical integration is used only to draw the fixed smooth profile; exact mean, convexity, flatness and reflection are proved above. The picture is a spatial section; physical time and covectors are supplied in (CB10)–(CB14).

6. Exercises with complete solutions

Exercise 1. Track the normal reflection gap

For a vertex PjP_j, compute the incoming and outgoing values of HqϕH_q\phi, including their signs. Explain why their magnitudes agree.

Solution. Write δj=(αj−1−αj)/2>0\delta_j=(\alpha_{j-1}-\alpha_j)/2>0. Since f′(xj)=tan⁡βjf'(x_j)=\tan\beta_j, Hqϕ=sin⁡α−tan⁡βjcos⁡α=sin⁡(α−βj)/cos⁡βjH_q\phi=\sin\alpha-\tan\beta_j\cos\alpha =\sin(\alpha-\beta_j)/\cos\beta_j. The incoming value is −sin⁡δj/cos⁡βj-\sin\delta_j/\cos\beta_j, the outgoing value is its positive opposite. Both are nonzero and tend to zero. This is the normal-root exchange for the defining function actually used, rather than a reflection of only the base velocity without its covector.

Exercise 2. Remove a false finite-order conclusion

Show that the reflected arc approaches the tangent state faster than every algebraic power of time, although it is interior on infinitely many intervals in every positive neighborhood of zero.

Solution. On the jj-th interval, s≥1/(j+1)s\ge1/(j+1). Multiplying the bounds (CB12) by s−Ns^{-N} gives at most CN(j+1)N2−jC_N(j+1)^N2^{-j}, or the still smaller analogous 4−j4^{-j} bound. Each tends to zero by its successive-term ratio. The same applies to the limiting one-sided momenta at the vertices. Yet sj↓0s_j\downarrow0, and strict chord separation in R0 gives an open interior interval (sj+1,sj)(s_{j+1},s_j) for every jj. Infinite-order contact therefore does not imply that the reflected curve equals the boundary curve on a positive neighborhood.

Exercise 3. Why a linear derivative interpolation fails

Suppose one replaces Hj(S)H_j(S) by SS. Determine the obstruction to keeping both endpoints of the same chord, and explain why the shifted smooth profile fixes it.

Solution. Linear interpolation has mean 1/21/2; its total height is dj(bj+bj+1)/2d_j(b_j+b_{j+1})/2. The required height is djmjd_jm_j. These agree only if θj=1/2\theta_j=1/2. Here, as j→∞j\to\infty, bj/mj→3/2b_j/m_j\to3/2 and bj+1/mj→3/4b_{j+1}/m_j\to3/4, by differentiability of tan⁡\tan and arctan⁡\arctan at zero. Thus θj→1/3\theta_j\to1/3, and the endpoint error is of the same scale as the height itself. Moreover linear interpolation does not match all higher derivatives at the vertices. The translated fixed profile has exact mean θj\theta_j and is constant near both endpoints, simultaneously preserving the chord and every derivative match.

Independent exposition, proofs, exercises, figure and scripts by GPT-6.1 Sol (OpenAI), October 2026, dedicated under CC0 1.0 Universal. Linked programme materials retain their recorded terms. The source credit concerns Hörmander III, Example 24.3.11; this construction does not assert approximation of every nonunique generalized arc by transverse broken rays.