Boundary reflection and compressed wavefronts

A singular wave meeting a noncharacteristic boundary has two normal covectors above the same tangential covector. At a transverse hyperbolic point, Dirichlet data couple these two branches. If the forcing and the boundary value are smooth there, regularity on either branch forces regularity on the other and at the boundary itself. This is the microlocal reflection law.

We prove that law for an arbitrary scalar second-order differential operator with smooth coefficients and real principal symbol. The lower-order coefficients may be complex. No global time function, positive energy, formal self-adjointness or globally Lorentzian signature is assumed. We then prove propagation along broken bicharacteristics with discrete transverse reflections. The construction of a distribution singular on a prescribed entire broken ray is a subsequent topic; glancing points require a different analysis.

The preceding proofs used here are Sections 7–9 of Higher-order Cauchy problems: roots, jets and propagation: complete simple-root factorization, commuting tangential tests, forced propagation, intrinsic jets and distributional Cauchy uniqueness. Sections 3–5 of Singularities along a real characteristic direction supply propagation on general interior characteristic pieces, including radial and stationary curves. The exact compressed coordinate law is in Stretched kernels and compressed geometry. The normal extension and tangential action proofs retain (GE26), all intrinsic actions and their precise domains; the smooth tester theorem supplies the two-sided criterion. Elliptic boundary regularity from rough data, Sections 2–7, proves the intrinsic normal inclusion, one common collar for all tangentially smoothing normal-jet errors, the full low-regularity bootstrap and the local differential theorem with redundant boundary data. Its equation (EW32) is exactly the elliptic identity used below. These are complete programme proofs at the exact versions in the proof map. References below to the global-boundary lesson mean the current compressed-geometry and normal-extension components just linked; the cited (GE26) is in Section 1 of the latter. The smooth parameter-dependent flow argument in Section 1 of Energy and existence with a timelike Dirichlet boundary also applies to the Hamilton field used below.

1. Two normal roots above one boundary covector

Let XX be a smooth manifold with boundary, of dimension n≥2n\ge2, and let PP be a scalar second-order differential operator with smooth coefficients. Its homogeneous principal symbol p(x,ξ)p(x,\xi) is real. Write qx(ξ,θ)q_x(\xi,\theta) for its symmetric polarization. The boundary is noncharacteristic:

p(x,ν)≠0(x∈∂X,0≠ν∈Nx∗∂X).(RF1) p(x,\nu)\ne0\qquad (x\in\partial X,\quad0\ne\nu\in N_x^*\partial X). \tag{RF1}
All arguments are local. On a connected small boundary patch, multiply PP by a constant sign so that its principal form is positive on boundary conormals. This does not change any wavefront or data-regularity assertion.

We first construct normal coordinates without assuming that qxq_x is nondegenerate on the entire cotangent space. Start with a boundary defining function r0≥0r_0\ge0, and prescribe on the boundary the covector

ξ0(y)=dr0p(y,dr0).(RF2) \xi_0(y)=\frac{dr_0}{\sqrt{p(y,dr_0)}}. \tag{RF2}
It has principal square one and annihilates boundary tangent vectors. Extend the coefficients across the boundary and flow this initial cotangent submanifold by V=Hp/2V=H_p/2. The local existence, uniqueness and smooth dependence proof cited above applies to this smooth vector field on a bounded cotangent coordinate box. It gives a smooth map (s,y)↦(x(s,y),ξ(s,y))(s,y)\mapsto(x(s,y),\xi(s,y)) for a common small interval. At s=0s=0,
Vr0=qy(ξ0,dr0)=p(y,dr0)>0.(RF3) Vr_0=q_y(\xi_0,dr_0)=\sqrt{p(y,dr_0)}>0. \tag{RF3}
The base map (s,y)↦x(s,y)(s,y)\mapsto x(s,y) therefore has an invertible differential: the yy-directions span the boundary and the ss-direction is transverse. The smooth inverse function theorem supplies product coordinates (r,y)(r,y), with r=sr=s, on a smaller neighborhood, and r≥0r\ge0 on XX.

Here is why the transported covector is drdr. Let θ=ξ dx\theta=\xi\,dx be the canonical one-form, using Hp=pξ∂x−px∂ξH_p=p_\xi\partial_x-p_x\partial_\xi. Then

ιVdθ=−12dp,θ(V)=p,LVθ=12dp.(RF4) \iota_Vd\theta=-\tfrac12dp,\qquad \theta(V)=p,\qquad \mathcal L_V\theta=\tfrac12dp. \tag{RF4}
The flow preserves p=1p=1. On the initial submanifold θ\theta vanishes in all yy-directions. The last identity in (RF4) shows that it remains zero in those directions, because the pullback of dpdp there is zero. Its value in the ss-direction is one. Consequently the pulled-back one-form is dsds, and the covector at x(s,y)x(s,y) is drdr. The coordinates yy are constant on these base flow curves, so
p(x,dr)=1,qx(dr,dyj)=0,p(r,y;ρ,η)=ρ2−R(r,y,η).(RF5) p(x,dr)=1,\qquad q_x(dr,dy_j)=0,\qquad p(r,y;\rho,\eta)=\rho^2-R(r,y,\eta). \tag{RF5}

The one-form calculation can also be verified directly, without invoking a formula for its Lie derivative. In a cotangent coordinate chart write Xj=∂yjxX_j=\partial_{y_j}x and Ξj=∂yjξ\Xi_j=\partial_{y_j}\xi along the flow. Differentiating the actual scalar pairing and then differentiating the quadratic Euler identity gives

dds(ξ⋅Xj)=−12px⋅Xj+12∑i,kξi(pξixk(Xj)k+pξiξk(Ξj)k)=12(px⋅Xj+pξ⋅Ξj)=12∂yjp=0,ξ⋅∂sx=12ξ⋅pξ=p=1.(RFA1) \begin{aligned} \frac{d}{ds}(\xi\cdot X_j) &=-\tfrac12p_x\cdot X_j +\tfrac12\sum_{i,k}\xi_i \bigl(p_{\xi_i x_k}(X_j)_k +p_{\xi_i\xi_k}(\Xi_j)_k\bigr)\\ &=\tfrac12\bigl(p_x\cdot X_j+p_\xi\cdot\Xi_j\bigr) =\tfrac12\partial_{y_j}p=0,\\ \xi\cdot\partial_s x&=\tfrac12\xi\cdot p_\xi=p=1. \end{aligned} \tag{RFA1}
Indeed ∑iξipξixk=2pxk\sum_i\xi_i p_{\xi_i x_k}=2p_{x_k} and ∑iξipξiξk=pξk\sum_i\xi_i p_{\xi_i\xi_k}=p_{\xi_k}, obtained by differentiating ∑iξipξi=2p\sum_i\xi_i p_{\xi_i}=2p in xkx_k and ξk\xi_k. The initial tangential pairings are zero. The displayed identities therefore prove that the pulled-back covector takes value zero on every ∂yj\partial_{y_j} and one on ∂s\partial_s, which characterizes dsds. Since the base differential is invertible, this proves ξ=dr\xi=dr in the base chart. Finally Vyj=qx(dr,dyj)=0V y_j=q_x(dr,dy_j)=0, because the transported base coordinates yjy_j are constant on the flow. This establishes both equalities before (RF5) without an invertibility hypothesis on the full quadratic form.

Here RR is a smooth real quadratic form in the tangential covector η\eta. No definiteness of RR is asserted. The leading normal coefficient of the differential operator is now exactly one. Its full expression is Dr2+C(r)Dr+B(r)D_r^2+C(r)D_r+B(r), where D=−i∂D=-i\partial, CC is multiplication by a smooth function and BB is a tangential differential family of order at most two. Its principal tangential symbol is −R-R.

At the boundary set R0(y,η)=R(0,y,η)R_0(y,\eta)=R(0,y,\eta), for η≠0\eta\ne0. The elliptic, hyperbolic and glancing regions are respectively

E={R0<0},H={R0>0},G={R0=0}.(RF6) \mathcal E=\{R_0<0\},\qquad \mathcal H=\{R_0>0\},\qquad \mathcal G=\{R_0=0\}. \tag{RF6}
At q=(y,η)∈Hq=(y,\eta)\in\mathcal H there are exactly two characteristic lifts:
ℓ±(q)=(0,y;±R0(y,η),η),Hpr=2ρ.(RF7) \ell_\pm(q)=(0,y;\pm\sqrt{R_0(y,\eta)},\eta),\qquad H_pr=2\rho. \tag{RF7}
The plus lift enters XX as the Hamilton parameter increases; the minus lift leaves it. These names refer to the chosen Hamilton parameter, which need not be physical time. On either branch, putting λ±=±R\lambda_\pm=\pm\sqrt R gives
p=(ρ−λ±)(ρ+λ±),Hp=2λ±Hρ−λ±on ρ=λ±.(RF8) p=(\rho-\lambda_\pm)(\rho+\lambda_\pm),\qquad H_p=2\lambda_\pm H_{\rho-\lambda_\pm} \quad\text{on }\rho=\lambda_\pm. \tag{RF8}
Thus rr parametrizes both short characteristic germs, with tangential flow generated by ∂r−Hλ±\partial_r-H_{\lambda_\pm}. The sign in (RF8) reverses one Hamilton orientation but leaves its unparameterized germ unchanged.

The root interchange has an invariant formula. If ξ\xi is characteristic over a boundary point and ν≠0\nu\ne0 is a boundary conormal, define

Rξ=ξ−2qx(ξ,ν)p(x,ν)ν.(RF9) \mathscr R\xi =\xi-2\frac{q_x(\xi,\nu)}{p(x,\nu)}\nu. \tag{RF9}
Restriction to Tx∂XT_x\partial X is unchanged. Expanding the quadratic form shows p(x,Rξ)=p(x,ξ)p(x,\mathscr R\xi)=p(x,\xi), while qx(Rξ,ν)=−qx(ξ,ν)q_x(\mathscr R\xi,\nu)=-q_x(\xi,\nu). Applying the formula twice gives R2=I\mathscr R^2=I. Rescaling ν\nu, or multiplying pp by a nonzero scalar, leaves the formula unchanged. In (RF5) it is precisely (ρ,η)↦(−ρ,η)(\rho,\eta)\mapsto(-\rho,\eta). At a glancing lift the two covectors coincide; transversality, rather than merely the existence of this algebraic involution, is what the proof needs.

2. What boundary regularity measures

In product coordinates, an interior covector ρ dr+η dy\rho\,dr+\eta\,dy has compressed components

π(r,y;ρ,η)=(r,y;rρ,η).(RF10) \pi(r,y;\rho,\eta)=(r,y;r\rho,\eta). \tag{RF10}
This is the dual coordinate system of the vector fields r∂r,∂yr\partial_r,\partial_y. At r=0r=0, both finite normal lifts of q=(y,η)q=(y,\eta) project to (0,y;0,η)(0,y;0,\eta). We identify this embedded nonzero tangential bundle with T∗∂X∖0T^*\partial X\setminus0. The intrinsic coordinate law, the conormal test space A\mathcal A, its dual A′\mathcal A', and the compressed wavefront WF⁡b\operatorname{WF}_b are the exact ones in the global-boundary prerequisite. In the interior the compression is invertible and WF⁡b\operatorname{WF}_b is the image of the ordinary wavefront.

Use the noncharacteristic extension class

N(X)={u∈A′(X):WF⁡b(u)∣∂X⊂T∗∂X∖0}.(RF11) \mathcal N(X)=\{u\in\mathcal A'(X): \operatorname{WF}_b(u)|_{\partial X} \subset T^*\partial X\setminus0\}. \tag{RF11}
For u∈Nu\in\mathcal N, intrinsic derivatives and properly supported tangential operators preserve this class; intrinsic normal jets exist. The precise differential and trace inclusions proved in the prerequisites are
WF⁡b(Dru)⊂WF⁡b(u),WF⁡(γku)⊂WF⁡b(u)∣∂X,γku=(Drku)∣r=0.(RF12) \operatorname{WF}_b(D_ru)\subset\operatorname{WF}_b(u), \qquad \operatorname{WF}(\gamma_k u) \subset\operatorname{WF}_b(u)|_{\partial X}, \quad \gamma_k u=(D_r^ku)|_{r=0}. \tag{RF12}
The derivative here is intrinsic. A supported representative satisfies
Dr(H(r)u)=H(r)Dru−iδ(r)γ0u.(RF13) D_r(H(r)u)=H(r)D_ru-i\delta(r)\gamma_0u. \tag{RF13}
Removing the displayed boundary delta is necessary when identifying the intrinsic derivative. We never replace it by the uncorrected derivative of a zero extension.

Two consequences will be used repeatedly. First, for w∈Nw\in\mathcal N, a properly localized tangentially smoothing family maps ww to a function smooth in a sufficiently small compact collar. Indeed tangential smoothing removes every nonzero tangential boundary covector, and (RF11) excludes the other boundary covectors. If interior wavefront points persisted arbitrarily close to the compact boundary support, their normalized compressed covectors would have a convergent subsequence, contradicting closedness of WF⁡b\operatorname{WF}_b. The output is therefore in A\mathcal A as well as A′\mathcal A'; the proved equality A∩A′=C∞\mathcal A\cap\mathcal A'=C^\infty supplies smoothness up to the boundary. Proper support and separated base supports contribute the same type of smoothing terms. This argument concerns a small collar; tangential smoothing alone does not smooth an arbitrary distribution with pure normal singularities.

Second, if a compact boundary cone misses WF⁡b(w)\operatorname{WF}_b(w), every tangential tester supported in a smaller such cone maps ww to a smooth function in a common smaller collar. The exact tangential tester theorem removes the nonzero tangential boundary points, and the preceding closedness argument finishes the collar. These statements apply to every intrinsic jet in (RF12). The common collar can be fixed before choosing the finite collection of errors: Sections 2–3 of the elliptic boundary regularity companion prove that the original input and all its intrinsic normal derivatives have no pure-normal interior covectors on one compact collar. Every tangentially smoothing normal-jet remainder is smooth on that same collar, by the full Fourier kernel estimate (EW7). For a finite family of source tests, use one fixed elliptic source tester and its full parameter-dependent tangential parametrix as in (EW9); the residuals obey that common-collar argument. Thus a separate order-dependent shrinking has not replaced smoothness on one open neighborhood. They give control over all normal covectors seen by the tangential test.

Suppose now that

u,f∈N(X),Pu=f,γ0u=b∈D′(∂X).(RF14) u,f\in\mathcal N(X),\qquad Pu=f,\qquad \gamma_0u=b\in\mathcal D'(\partial X). \tag{RF14}
At an elliptic boundary covector, the exact result is
E∩WF⁡b(u)=E∩(WF⁡b(f)∪WF⁡(b)).(RF15) \mathcal E\cap\operatorname{WF}_b(u) =\mathcal E\cap \bigl(\operatorname{WF}_b(f)\cup\operatorname{WF}(b)\bigr). \tag{RF15}
Here and below the boundary restriction of WF⁡b(f)\operatorname{WF}_b(f) is understood. To apply the proved elliptic-boundary theorem, freeze (RF5) at an elliptic qq. With a=−R0>0a=\sqrt{-R_0}>0, the bounded normal solutions of
(Dr2+a2)v=0(r≥0) (D_r^2+a^2)v=0\quad(r\ge0)
form the one-dimensional space spanned by e−are^{-ar}. The Dirichlet evaluation sends this solution to one, hence is an isomorphism on that space. There is no real characteristic lift. Thus qq is outside the boundary characteristic set in Section 7 of the elliptic boundary regularity companion. Its complete coefficient completion and regularity proof give the inclusion from left to right in (RF15), for this locally elliptic cone even when the original operator is not elliptic elsewhere. The ordinary differential and trace inclusions (RF12) give the reverse inclusion. The input retains the original operator's full lower-order terms and the class N\mathcal N; no elliptic global boundary hypothesis is added to (RF14).

3. Two factorized equations with the same boundary test

Fix q0∈Hq_0\in\mathcal H. Choose a small compact base-cone neighborhood of its positive tangential ray such that

R(r,y,η)≥c∣η∣2>0,0≤r≤ε,(RF16) R(r,y,\eta)\ge c|\eta|^2>0,\qquad 0\le r\le\varepsilon, \tag{RF16}
on the slightly larger cone that will contain both transported germs. This follows by continuity and homogeneity after shrinking. The positive root R\sqrt R is a smooth homogeneous symbol of order one there, with a uniform separation 2R≥2c∣η∣2\sqrt R\ge2\sqrt c|\eta|.

The bounded simple-root calculus applies after a local extension. Choose an angular and base cutoff supported where (RF16) holds, equal to one near the chosen cone, and replace RR outside it by a fixed positive multiple of ∣η∣2|\eta|^2. The convex combination remains positive with a uniform lower bound. Complete it smoothly at bounded frequencies and extend the full lower-order coefficients with proper base cutoffs. This makes a bounded tangential symbol family with two globally separated real principal roots. It agrees with the original complete differential family on the smaller base-cone region. After a tangential localization, the discrepancies have only tangentially smoothing coefficients multiplying the displayed finite normal jets. They act smoothly by Section 2. Hence the local extension is a device for applying the exact calculus, rather than a change to the equation under study.

Section 7 of the higher-order Cauchy lesson, applied to the two simple roots with m=2m=2, gives two factorizations:

P=(Dr−A+)(Dr−A−)+Ω,P=(Dr−A~−)(Dr−A~+)+Ω~,(RF17) \begin{aligned} P&=(D_r-A_+)(D_r-A_-)+\Omega,\\ P&=(D_r-\widetilde A_-)(D_r-\widetilde A_+) +\widetilde\Omega, \end{aligned} \tag{RF17}
microlocally in the chosen tangential cone. The remainders have tangential order minus infinity and no normal derivative. The principal symbols of A±A_\pm and A~±\widetilde A_\pm are ±R\pm\sqrt R. The two full operators of a given sign can differ at order zero. Interchanging factors without recomputing those lower terms would not give (RF17).

For clarity, the receiving use of the factorization algorithm retains its full remainder. If the present defect is a tangential RN∈S2−NR_N\in S^{2-N}, correcting the left root by a symbol of order 1−N1-N cancels its leading value at that root, since the derivative of the principal polynomial there is ±2R\pm2\sqrt R. Exact left polynomial division supplies the order-1−N1-N correction to the right factor, including derivatives of its coefficients in rr. The remaining defect is tangential of order 1−N1-N. Repeating and summing with bounds on every rr-derivative gives actual symbols; comparison with every finite stage makes the final defect belong to all negative tangential orders. This is the m=2m=2 instance of the complete ordered algorithm (HC16)–(HC17), not a formal root factorization.

For this second-order operator the absence of a normal derivative in the remainder has an exact algebraic check. Set A+=−C−A−A_+=-C-A_- after choosing the right factor. Normal differentiation acts on the coefficient family as well as on its input, so the complete operator identity is

(Dr−A+)(Dr−A−)=Dr2−(A++A−)Dr+i∂rA−+A+A−,P−(Dr−A+)(Dr−A−)=B−i∂rA−+CA−+A−2.(RFA2) \begin{aligned} (D_r-A_+)(D_r-A_-) &=D_r^2-(A_++A_-)D_r+i\partial_rA_-+A_+A_-,\\ P-(D_r-A_+)(D_r-A_-) &=B-i\partial_rA_-+CA_-+A_-^2. \end{aligned} \tag{RFA2}
Every product here is an actual ordered tangential operator product. Its normal degree is zero, before any symbolic approximation. The leading tangential equation is a−2−R=0a_-^2-R=0. If the remaining symbol is eNe_N of order 2−N2-N, a correction of order 1−N1-N with leading value −eN/(2a−)-e_N/(2a_-) cancels that defect. The gap estimate (RF16) justifies this division and its differentiated symbol bounds. Derivatives of the correction and lower composition terms have order at most 1−N1-N; its square has that order as well for N≥1N\geq1. Iteration with all normal-parameter derivatives and the proved symbol summation gives a full A−A_-. Reimpose A+=−C−A−A_+=-C-A_- on the summed operators. Formula (RFA2) then still cancels the normal coefficient exactly, and comparison with every finite stage makes the remaining tangential symbol smoothing. For the second ordering start with the right factor A~+\widetilde A_+, whose principal root is +R+\sqrt R, and set A~−=−C−A~+\widetilde A_-=-C-\widetilde A_+. This repeats the same ordered calculation with the other root; it does not interchange two already constructed full factors.

Put

L+=Dr−A+,L~−=Dr−A~−.(RF18) L_+=D_r-A_+,\qquad \widetilde L_-=D_r-\widetilde A_-. \tag{RF18}
There are proper tangential order-zero families Q+(r)Q_+(r) and Q~−(r)\widetilde Q_-(r), with microsupport in arbitrarily small tubes transported by their respective principal flows, such that
[L+,Q+]∈Ψtan−∞,[L~−,Q~−]∈Ψtan−∞,Q+(0)=Q~−(0)=Q0.(RF19) [L_+,Q_+]\in\Psi^{-\infty}_{\mathrm{tan}},\qquad [\widetilde L_-,\widetilde Q_-] \in\Psi^{-\infty}_{\mathrm{tan}},\qquad Q_+(0)=\widetilde Q_-(0)=Q_0. \tag{RF19}
Choose Q0Q_0 elliptic at q0q_0, with arbitrarily small microsupport. The principal commuting equation is
(∂r−Hλ+)q+,0=0,(∂r−Hλ−)q~−,0=0.(RF20) (\partial_r-H_{\lambda_+})q_{+,0}=0,\qquad (\partial_r-H_{\lambda_-})\widetilde q_{-,0}=0. \tag{RF20}
It transports the same initial symbol along the two different flows. To retain the full prescribed initial operator, begin with a symbol expansion of Q0Q_0. At each successive order, solve the inhomogeneous scalar transport equation with the corresponding prescribed initial coefficient; the leading commutator is −i(∂r−Hλ)-i(\partial_r-H_\lambda) on that correction. All other contributions drop one order, because the symbol is scalar and the difference between the full generator and its real order-one principal part has order zero. The parameter-aware asymptotic sum gives the full commuting test, with its initial value equal to Q0Q_0 modulo tangential smoothing. Subtract that initial smoothing difference as a smooth rr-family constant near zero. Its commutator is still tangentially smoothing. Thus the last equality in (RF19) is an exact operator equality. Proper cutoff errors retain the same smoothing qualification on a smaller tube. This is precisely the commuting-test construction in Section 8 of the Cauchy lesson, with its initial condition specified here.

Assume q0q_0 misses the forcing and boundary-data wavefronts. Shrink the tubes and collar so that their boundary closures miss those sets. Define two intrinsic distributions

v=Q+(Dr−A−)u,v~=Q~−(Dr−A~+)u.(RF21) v=Q_+(D_r-A_-)u,\qquad \widetilde v=\widetilde Q_-(D_r-\widetilde A_+)u. \tag{RF21}
All derivatives and tangential actions preserve N\mathcal N. Equations (RF17)–(RF19) give the actual equations
L+v=Q+f−Q+Ωu+[L+,Q+](Dr−A−)u,L~−v~=Q~−f−Q~−Ω~u+[L~−,Q~−](Dr−A~+)u.(RF22) \begin{aligned} L_+v&=Q_+f-Q_+\Omega u +[L_+,Q_+](D_r-A_-)u,\\ \widetilde L_-\widetilde v &=\widetilde Q_-f-\widetilde Q_-\widetilde\Omega u +[\widetilde L_-,\widetilde Q_-] (D_r-\widetilde A_+)u. \end{aligned} \tag{RF22}
Every term on the right is smooth in a smaller collar. For the forcing terms this is the compact-cone tester assertion of Section 2. For the other terms it is tangential smoothing on N\mathcal N, with the indicated intrinsic jets retained.

Here is the support convention for the global first-order Cauchy theorem used below. Extend the two root families to bounded symbols on the whole tangential chart space, as above, and realize them with a uniform proper kernel cutoff. Choose the initial test with compact output base support; its complete transported symbol remains in a fixed compact base set on the short normal interval. Its proper realization and the initial smoothing correction can have output support in a slightly larger fixed compact set. The equations (RF22) then have compact tangential output support in a further fixed set: the first-order operator applied to such a test enlarges it only by the fixed proper kernel range. Include all these input and output sets in the common collar just proved. The full-symbol comparison with the original operator adds only the already controlled tangentially smoothing normal-jet terms. The right sides are therefore smooth functions in every global tangential Sobolev space, with all normal derivatives, on each closed shorter interval. The smooth Cauchy solutions are legitimate inputs to the all-order theorem, while the original tested distributions have compact tangential support. This supplies the compact-time Schwartz test topology in the distributional uniqueness argument; no spatially unbounded smooth datum is silently assigned a finite global Sobolev norm.

At the boundary, (RF19) and the intrinsic trace action give

γ0v~−γ0v=Q0(A−(0)−A~+(0))b.(RF23) \gamma_0\widetilde v-\gamma_0v =Q_0\bigl(A_-(0)-\widetilde A_+(0)\bigr)b. \tag{RF23}
The right side is smooth because bb is regular on the microsupport of Q0Q_0. Both tests must have the same full boundary operator for the normal derivative of uu to cancel exactly. Principal-symbol equality alone would leave an uncontrolled lower-order boundary term.

4. Regularity of either germ determines the reflection

Transverse reflection theorem. Under (RF1) and (RF14), let

q0∈H∖(WF⁡b(f)∪WF⁡(b)).(RF24) q_0\in\mathcal H\setminus \bigl(\operatorname{WF}_b(f)\cup\operatorname{WF}(b)\bigr). \tag{RF24}
There is a sufficiently short source-free collar in which the following three assertions are equivalent: q0∉WF⁡b(u)q_0\notin\operatorname{WF}_b(u); the inward characteristic germ from ℓ+(q0)\ell_+(q_0) is regular; the outward characteristic germ from ℓ−(q0)\ell_-(q_0) is regular. Equivalently, a singular boundary point under (RF24) has singularities on both germs. Each germ remains entirely singular or entirely regular along any connected interior continuation that avoids the forcing wavefront before its next boundary encounter. General interior propagation is Theorem 5.1 of the real-characteristic lesson; the simple-root argument alone is needed only in the small normal collar.

Proof. If q0q_0 is regular in the compressed sense, its regular set contains a compressed conic neighborhood. Both lifted germs converge to this point under (RF10), so their short interior pieces are regular. The general interior propagation theorem extends their regularity along source-free continuations.

Conversely suppose the plus germ is regular. Use the first distribution vv in (RF21). Its equation has smooth right side by (RF22). Away from pure normal covectors, its first-order principal symbol is ρ−R\rho-\sqrt R, so ordinary elliptic regularity confines WF⁡(v)\operatorname{WF}(v) in the small transported tube to that graph. The no-pure-normal and ordinary/tangential compatibility arguments in Section 8 of the Cauchy lesson justify this use of ordinary ellipticity on a compact subcollar with r>0r>0. First shrink the compact input and output collars: (RF11), (RF12) and closedness exclude pure normal interior wavefront covectors of uu and its relevant intrinsic jets there, since such covectors would otherwise accumulate at a forbidden pure compressed normal boundary point. On a compact subcollar away from zero, proper localization consequently has bounded normal slopes and smooth normal slices in a common negative tangential Sobolev scale. Tangential spatially smoothing tails are fully smooth on these localized inputs. Moreover WF⁡(v)⊂WF⁡(u)\operatorname{WF}(v)\subset\operatorname{WF}(u) on this smaller tube, modulo the already smooth terms. A point on the regular plus germ is therefore regular for vv at its graph lift. This interior-slice argument does not assume a uniform ordinary normal-frequency bound all the way to r=0r=0.

Choose a small positive slice r=r1r=r_1. Shrink the initial microsupport of Q0Q_0 so that the entire graph over the transported support of Q+(r1)Q_+(r_1) lies in a regular neighborhood of that point. The other normal lifts over this support are regular by first-order ellipticity. High normal slopes are regular by the preceding localization. The noncharacteristic pullback theorem now gives a smooth boundary-cone slice of vv at r1r_1. Since Q+Q_+ is smoothing outside the tube, the properly localized output slice is fully smooth. This step uses regularity over every normal lift seen by the spatial test, rather than taking a slice from regularity at only one arbitrary covector.

Solve L+w=L+vL_+w=L_+v backward from this smooth slice. The scalar first-order Cauchy theorem at every Sobolev order makes ww smooth down to r=0r=0, and the distributional uniqueness argument in Section 9 of the Cauchy lesson identifies it with vv on the smaller collar. That argument applies after reversing the normal coordinate and retains its compact-time Schwartz test topology. Thus vv is smooth up to the boundary and γ0v\gamma_0v is smooth. Equation (RF23) makes γ0v~\gamma_0\widetilde v smooth as well. The smooth forward Cauchy solution of L~−w~=L~−v~\widetilde L_-\widetilde w=\widetilde L_-\widetilde v agrees with v~\widetilde v by the same distributional uniqueness proof. We have proved

v,v~∈C∞in a smaller collar.(RF25) v,\widetilde v\in C^\infty\quad\text{in a smaller collar}. \tag{RF25}

Recovering uu is a separate step. Both Q+Q_+ and Q~−\widetilde Q_- are elliptic over a common sufficiently small base-cone neighborhood of q0q_0 for rr near zero. Choose a tangential localizer a(r,y,Dy)a(r,y,D_y), elliptic there, with microsupport inside that common region. Proper tangential parametrices for the two tests, applied to (RF21) and (RF25), give

a(Dr−A−)u∈C∞,a(Dr−A~+)u∈C∞.(RF26) a(D_r-A_-)u\in C^\infty,\qquad a(D_r-\widetilde A_+)u\in C^\infty. \tag{RF26}
The parametrix remainders act on the actual N\mathcal N jets, so Section 2 turns them into smooth functions in a smaller collar. Subtract the two displayed outputs. The normal derivative cancels and gives
a(A~+−A−)u∈C∞,σ1(A~+−A−)=2R≠0.(RF27) a(\widetilde A_+-A_-)u\in C^\infty,\qquad \sigma_1(\widetilde A_+-A_-)=2\sqrt R\ne0. \tag{RF27}
An order-minus-one tangential parametrix, with a further smaller conic cutoff, gives a smooth tangential order-zero elliptic test of uu. Its smoothing errors again act smoothly on N\mathcal N. The exact tester criterion removes q0q_0 from WF⁡b(u)\operatorname{WF}_b(u). This proves the converse from the plus germ.

If the minus germ is regular, start with v~\widetilde v, solve its Cauchy equation backward from a small positive slice, and use (RF23) to obtain the initial value of vv. Solve that equation forward and repeat (RF26)–(RF27). This proves the converse from the minus germ, without reversing the order of the full factors. Interior continuation uses the full real-symbol propagation theorem stated in the prerequisites, with m=2m=2; it does not require a normal simple-root chart along the entire continuation. The theorem follows.

The proof also explains why a statement only about γ1u\gamma_1u would be insufficient. The boundary comparison (RF23) transmits initial regularity between the two scalar equations; the equations transmit that regularity throughout a collar; and the elliptic difference of their full outputs recovers uu there. All three steps are needed to reach compressed boundary regularity.

5. Propagation along a broken bicharacteristic

A transverse broken bicharacteristic is a continuous compressed curve γ:I→T~∗X∖0\gamma:I\to\widetilde T^*X\setminus0, with discrete reflection parameters, having the following properties. Between reflection parameters it is the compression of a characteristic HpH_p-trajectory, with its Hamilton orientation and positive scaling consistently retained. At a reflection parameter its two limits have full lifts ξ−\xi_- and ξ+\xi_+ over one boundary base point, satisfying

ξ+=Rξ−,ξ−∣T∂X=ξ+∣T∂X=η≠0,(y,η)∈H.(RF28) \xi_+=\mathscr R\xi_-,\qquad \xi_-|_{T\partial X}=\xi_+|_{T\partial X}=\eta\ne0,\qquad (y,\eta)\in\mathcal H. \tag{RF28}
One piece leaves XX and the next enters it. The two limits of (RF10) agree, so the compressed curve is continuous. A change of smooth positive parameter on an interior piece does not alter the curve. Discreteness here means local finiteness inside II: every compact subinterval contains only finitely many reflections. Accumulation at an endpoint outside II is not an interior reflection in this definition. Glancing contact is excluded.

Define the source set in the compressed bundle by

Σ=WF⁡b(f)∪WF⁡(b),(RF29) \Sigma=\operatorname{WF}_b(f) \cup\operatorname{WF}(b), \tag{RF29}
where the second set lies in the embedded tangential boundary bundle. Let J⊂IJ\subset I be a connected interval with γ(J)∩Σ=∅\gamma(J)\cap\Sigma=\varnothing. Then
γ(J)⊂WF⁡b(u)orγ(J)∩WF⁡b(u)=∅.(RF30) \gamma(J)\subset\operatorname{WF}_b(u) \quad\text{or}\quad \gamma(J)\cap\operatorname{WF}_b(u)=\varnothing. \tag{RF30}

Proof. Pull back the regular set along γ\gamma. It is open because the wavefront is closed. At an interior parameter, Theorem 5.1 of the real-characteristic lesson makes regularity constant on a small interval: PP is properly supported because it is differential, has a real homogeneous principal symbol of order two, and its ordinary forcing wavefront agrees with the compressed one in the interior. That theorem retains arbitrary complex lower terms, and its radial or stationary cases use conicity and flow uniqueness rather than an impossible normal straightening. At a reflection parameter, Section 4 makes regularity on either adjacent germ equivalent to regularity at the boundary and on the other germ. Thus regularity is constant on a small parameter neighborhood there too. The singular set is also open in JJ: a singular interior point propagates in both directions locally, and a singular reflection has both neighboring germs singular. The two disjoint sets are open, cover JJ, and connectedness forces one to be empty. Equivalently, any two parameters in JJ can be joined through finitely many interior pieces and reflections on the compact interval between them. This proves (RF30), including noncompact connected JJ.

The conclusion concerns the whole compressed curve. At a reflection the ordinary interior wavefront has two different normal covectors. It becomes a single continuous curve only after the normal covector is compressed. The elliptic equality (RF15) and the transverse law (RF30) account for the two open regions E\mathcal E and H\mathcal H. Neither supplies propagation at G\mathcal G.

6. Flat waves and a trace that misses singularities

For the flat wave operator on r≥0r\ge0, use coordinates (r,t,z)(r,t,z) and

P=Dr2+∣Dz∣2−Dt2,R=τ2−∣η∣2.(RF31) P=D_r^2+|D_z|^2-D_t^2,\qquad R=\tau^2-|\eta|^2. \tag{RF31}
Its hyperbolic boundary region is ∣τ∣>∣η∣|\tau|>|\eta|. Reflection keeps (τ,η)(\tau,\eta) and changes the sign of ρ\rho. Its glancing set is ∣τ∣=∣η∣|\tau|=|\eta|. The elliptic region is ∣τ∣<∣η∣|\tau|<|\eta|. For a characteristic lift with τ>0\tau>0, physical time has derivative Hpt=−2τH_pt=-2\tau. Consequently the physical-time slope is
drdt=−ρτ.(RF32) \frac{dr}{dt}=-\frac{\rho}{\tau}. \tag{RF32}
The plus normal root is incoming in increasing physical time, although it is inward in increasing HpH_p-parameter. This distinction prevents an orientation error.

In two dimensions the exact Dirichlet wave

u(r,t)=δ(t−r)−δ(t+r)(RF33) u(r,t)=\delta(t-r)-\delta(t+r) \tag{RF33}
satisfies Pu=0Pu=0 and γ0u=0\gamma_0u=0. Its two hypersurfaces are disjoint for r>0r>0, so there is no interior cancellation. Their conormal covectors are
λ(dt−dr)(t=r),λ(dt+dr)(t=−r),λ≠0.(RF34) \lambda(dt-dr)\quad(t=r),\qquad \lambda(dt+dr)\quad(t=-r),\qquad \lambda\ne0. \tag{RF34}
The asserted wavefronts have a direct Fourier verification. A localization of δ(t)\delta(t) at zero has constant Fourier transform equal to the nonzero cutoff value, so both nonzero one-dimensional covector directions occur. Away from zero it vanishes. Tensoring with the smooth constant in the other coordinate gives exactly zero covector in that coordinate, by the complete tensor/projection theorem in the real-characteristic lesson. The invertible linear maps (r,t)↦(r,t−r)(r,t)\mapsto(r,t-r) and (r,t)↦(r,t+r)(r,t)\mapsto(r,t+r), with their exact cotangent transforms, give (RF34). Their distinct interior supports prevent cancellation. Differentiation and these submersion pullbacks are continuous on distributions, so the displayed wave equation follows from the same derivative calculation as for an arbitrary smooth one-variable function.

They compress to the same boundary covectors (t=0;τ=λ)(t=0;\tau=\lambda). The distribution is smooth as an rr-dependent tangential distribution and is extendible across zero by its displayed full-space formula. Since Pu=0Pu=0 and the boundary is noncharacteristic, (GE26) of the global-boundary lesson supplies its intrinsic representative in N\mathcal N. Its intrinsic traces agree with the displayed tangential-distribution limits, by interior injectivity and the trace theorem. Direct differentiation yields

γ1u=2iδ′(t).(RF35) \gamma_1u=2i\delta'(t). \tag{RF35}
The singular normal trace records the reflected pair even though the Dirichlet trace vanishes. Closedness gives the boundary points of its compressed wavefront, while (RF34) gives both neighboring germs. This verifies the theorem with zero forcing and zero boundary value.

The flat reflected wave and the exact compressed normal-coordinate section

The upper panel is the base projection of (RF33)–(RF34), in the positive covector section τ=1\tau=1. Arrows follow increasing physical time, rather than increasing HpH_p-parameter. The incoming lift has ρ=1\rho=1, the outgoing lift ρ=−1\rho=-1, and both have tangential component one. The lower panel plots the exact compressed components ζ=r\zeta=r and ζ=−r\zeta=-r from (RF10); it omits the base time component, which is −r-r and rr respectively. Both full branches meet the single boundary point (r,t;ζ,τ)=(0,0;0,1)(r,t;\zeta,\tau)=(0,0;0,1). These are coordinate sections, with no assertion that the uncompressed covectors agree. Positive dilation recovers the whole positive cone; (RF34) also has its negative cone. The general reflection theorem is Section 4, following the mathematical antecedent of Hörmander III, Theorem 24.2.1, printed 424–425.

There is also a useful warning in the opposite direction. Set a(r)=e−1/r2a(r)=e^{-1/r^2} for r>0r>0, smoothly continued by zero at r=0r=0, and

w(r,t)=a(r)δ(t).(RF36) w(r,t)=a(r)\delta(t). \tag{RF36}
Every intrinsic normal trace is zero, because each derivative of aa is flat at zero. The class assertion can be checked directly. A conormal test of any fixed order, and each of its finitely many required tangential derivatives, has at most a fixed power growth in 1/r1/r on a compact collar. Pairing its tangential restriction with δ(t)\delta(t) and integrating against a(r)a(r) is continuous in that test topology, since aa decays faster than every power. This proves the dual conormal requirement. For an actual compressed tester, the needed uniform estimate can be given directly. Multiply aa by a compact normal cutoff equal to one near zero and extend by zero to negative rr; denote the resulting smooth compact function by AA. Every derivative is zero at zero: successive derivatives of e−1/r2e^{-1/r^2} are that exponential times a polynomial in 1/r1/r, and the series inequality ex≥xj/j!e^x\ge x^j/j! with arbitrarily large jj proves the assertion. Repeated integration by parts makes A^(ρ)\widehat A(\rho) rapidly decreasing with every prescribed power.

Fix a compressed boundary covector whose normal component ζ\zeta is nonzero. Choose a smooth order-zero full compressed symbol c(r,t,η,ζ)c(r,t,\eta,\zeta), elliptic there, with its high-frequency support in ∣η∣≤C∣ζ∣|\eta|\le C|\zeta|. The local realization and full conic cutoff proofs in the global-boundary calculus provide such a tester, with a full symbol supported in a slightly larger such cone modulo a residual symbol. In the left coordinate quantization its action on the compact input A(r)δ(t)A(r)\delta(t), before harmless outer cutoffs, is

(Tcw)(r,t)=(2π)−2∫ei(tη+rρ)c(r,t,η,rρ)A^(ρ) dη dρ.(RFA3) (T_cw)(r,t)=(2\pi)^{-2} \int e^{i(t\eta+r\rho)}c(r,t,\eta,r\rho) \widehat A(\rho)\,d\eta\,d\rho. \tag{RFA3}
The formula is first justified with compact frequency cutoffs. For 0≤r≤r00\le r\le r_0, the main cone has ∣η∣≤C(1+r∣ρ∣)|\eta|\le C(1+r|\rho|). Its η\eta-length is at most C′(1+∣ρ∣)C'(1+|\rho|). Any fixed number of ordinary r,tr,t derivatives of the integrand adds at most a fixed polynomial in ∣η∣+∣ρ∣|\eta|+|\rho|: differentiating rρr\rho contributes a factor ρ\rho, while every base or frequency derivative of the order-zero symbol has its uniform symbol bound. On that cone the absolute differentiated integral is thus bounded by an integral of CN(1+∣ρ∣)N∣A^(ρ)∣C_N(1+|\rho|)^N|\widehat A(\rho)|, finite for each fixed NN. For the residual symbol off the cone, take its η\eta-decay order greater than all requested derivative orders plus two; its integral has the same type of polynomial ρ\rho-bound. These estimates are uniform down to zero. Dominated convergence removes the frequency cutoffs and gives every ordinary derivative continuously up to r=0r=0. The output is smooth.

One may make the operator properly supported by compact input and output coordinate multipliers; take the input multiplier one on the compact support already fixed. The local lacunary realization, its residual discrepancy and the coordinate identification are precisely the full local-to-global operator statements GB:C3–C6. All residual terms act smoothly on this dual input by GB:D1 and the dual/conormal intersection theorem; their complete symbol estimates were also allowed above. A smooth output is an admissible conormal tester output. Hence the selected covector is absent from the compressed wavefront. Since every point outside the equator has ζ≠0\zeta\ne0, this proves w∈Nw\in\mathcal N using the actual compressed quantization, rather than assuming that a logarithmic-coordinate Fourier transform alone defines the boundary wavefront.

Nevertheless ww is singular at every (r,0;0,τ)(r,0;0,\tau), r>0r>0, τ≠0\tau\ne0. These wavefronts approach (0,0;0,τ)(0,0;0,\tau) in the compressed bundle, making that boundary point singular by closedness. Its forcing is

Pw=−a′′(r)δ(t)+a(r)δ′′(t).(RF37) Pw=-a''(r)\delta(t)+a(r)\delta''(t). \tag{RF37}
For every r>0r>0, this is a nonzero combination of δ\delta and δ′′\delta''; its Fourier polynomial is −a′′−aτ2-a''-a\tau^2, whose leading coefficient is nonzero. It is singular in both positive and negative tangential directions. Hence its compressed boundary wavefront contains the same candidate. The source-regularity hypothesis (RF24) fails. Vanishing traces therefore do not imply compressed regularity, and this example does not contradict the theorem.

7. Exercises with complete solutions

Exercise 1. Reflection before normalizing coordinates

At a fixed boundary point let

p(ρ,η)=aρ2+2ρℓ(η)+c(η),a≠0,(RF38) p(\rho,\eta)=a\rho^2+2\rho\ell(\eta)+c(\eta), \qquad a\ne0, \tag{RF38}
where ℓ\ell is linear and cc is quadratic. Compute the two normal roots, the transverse region and the reflection map. Verify all assertions without setting the mixed term to zero.

Solution. The discriminant quantity is Δ(η)=ℓ(η)2−ac(η)\Delta(\eta)=\ell(\eta)^2-a c(\eta). For Δ>0\Delta>0, the roots are

ρ±=−ℓ(η)±Δa,∂ρp(ρ±,η)=±2Δ.(RF39) \rho_\pm=\frac{-\ell(\eta)\pm\sqrt\Delta}{a}, \qquad \partial_\rho p(\rho_\pm,\eta)=\pm2\sqrt\Delta. \tag{RF39}
They are distinct and transverse. Reflection is ρ↦−ρ−2ℓ/a\rho\mapsto-\rho-2\ell/a, fixing η\eta. Its square is the identity and it interchanges the roots. Substitution in (RF38) preserves pp, and it changes the sign of aρ+ℓa\rho+\ell, exactly as (RF9) requires. For Δ<0\Delta<0 there is no real lift. For Δ=0\Delta=0 the root has zero normal velocity and the transverse theorem cannot be applied. Adding a multiple of the defining conormal changes the chosen normal coordinate ρ\rho, but q(ξ,ν)2−p(ν)p(ξ)=Δq(\xi,\nu)^2-p(\nu)p(\xi)=\Delta is unchanged under ξ↦ξ+sν\xi\mapsto\xi+s\nu. Thus these three regions are intrinsic; scaling ν\nu scales the discriminant by its positive square.

Exercise 2. The reflection law for an anisotropic wave

For a positive definite symmetric matrix

G=(ghThK),p(ρ,τ,η)=gρ2+2ρhTη+ηTKη−τ2,(RF40) G=\begin{pmatrix}g&h^T\\h&K\end{pmatrix}, \qquad p(\rho,\tau,\eta)=g\rho^2+2\rho h^T\eta +\eta^TK\eta-\tau^2, \tag{RF40}
find the reflected normal covector and the normal physical-time velocities. Describe the hyperbolic region.

Solution. Positive definiteness gives g>0g>0 and the Schur complement K−hhT/g>0K-hh^T/g>0: for z≠0z\ne0, evaluate the positive quadratic form on (−hTz/g,z)(-h^Tz/g,z). Put

A(τ,η)=τ2−ηT(K−hhT/g)η.(RF41) A(\tau,\eta)=\tau^2-\eta^T(K-hh^T/g)\eta. \tag{RF41}
Completing the square gives p=g(ρ+hTη/g)2−Ap=g(\rho+h^T\eta/g)^2-A. The hyperbolic region is A>0A>0, with roots −hTη/g±A/g-h^T\eta/g\pm\sqrt{A/g}. Formula (RF9) gives ρ′=−ρ−2hTη/g\rho'=-\rho-2h^T\eta/g, rather than generally −ρ-\rho. Since ∂ρp=2(gρ+hTη)\partial_\rho p=2(g\rho+h^T\eta) and ∂τp=−2τ\partial_\tau p=-2\tau, the two normal physical-time velocities are ∓gA/τ\mp\sqrt{gA}/\tau. Here τ≠0\tau\ne0 in the hyperbolic region, by positivity of the Schur complement. The magnitudes agree and the signs reverse. This is reflection in the principal metric; using a Euclidean covector reflection would retain an incorrect mixed term.

Exercise 3. An elliptic Dirichlet singularity is confined to the boundary

Let P=Dr2+∣Dy∣2P=D_r^2+|D_y|^2 on r≥0r\ge0, and let bb be a compactly supported tangential distribution. Define uu by the tangential Fourier formula

u^(r,η)=e−r∣η∣b^(η).(RF42) \widehat u(r,\eta)=e^{-r|\eta|}\widehat b(\eta). \tag{RF42}
Prove the equation, the trace and interior smoothness, then determine its compressed boundary wavefront.

Solution. The Fourier transform of bb has polynomial growth. For r>0r>0, differentiating (RF42) any number of times gives polynomial factors times an exponentially decreasing function at infinity. Fourier inversion is therefore smooth on each strip r≥δ>0r\ge\delta>0, including all tangential and normal derivatives. At r=0r=0, dominated integration against a Schwartz tangential test gives u(r)→bu(r)\to b in distributions, and every normal derivative has a distributional limit with Fourier transform (−∣η∣)kb^(-|\eta|)^k\widehat b. The multiplier at zero frequency is locally integrable with at most polynomial growth, so these limits exist for each kk. Distributional differentiation gives Dr2u^=−∣η∣2u^D_r^2\widehat u=-|\eta|^2\widehat u, proving Pu=0Pu=0 for r>0r>0.

The family is extendible across zero as a distribution: integrate its continuous tangential-distribution values against a compact smooth normal test on r≥0r\ge0, and set the extension zero for r<0r<0. This is continuous on compact smooth tests because the polynomial Fourier bound is uniform for r≥0r\ge0. The equation has zero source, the boundary is noncharacteristic, and this interior distribution is extendible. The exact noncharacteristic extension theorem (GE26) in Section 5 of the global-boundary lesson therefore supplies its unique intrinsic extension in N\mathcal N. No raw derivative of the zero extension is used to assert the equation across the boundary. Here R0=−∣η∣2R_0=-|\eta|^2, so every nonzero tangential covector is elliptic. Apply (RF15) with f=0f=0 and γ0u=b\gamma_0u=b. It gives exactly WF⁡b(u)∣∂X=WF⁡(b)\operatorname{WF}_b(u)|_{\partial X}=\operatorname{WF}(b). Interior smoothness shows that the entire compressed wavefront is this boundary set. The growing normal exponential is excluded by the bounded-solution Dirichlet map used in Section 2; it was not used to manufacture interior regularity.

Exercise 4. A reflected ray in a strip

For P=Dr2−Dt2P=D_r^2-D_t^2 on 0≤r≤L0\le r\le L, L>0L>0, parametrize a characteristic that starts at (t,r)=(0,r∗)(t,r)=(0,r_*), 0<r∗<L0<r_*<L, with tangential covector τ=1\tau=1 and initial ρ=−1\rho=-1, using increasing physical time. Reflect it at each endpoint of the strip. Give its first three reflections and prove the wavefront alternative when both Dirichlet values and the forcing are regular along the resulting finite broken segment.

Solution. The physical velocity is dr/dt=−ρdr/dt=-\rho, so the first piece has r=r∗+tr=r_*+t, ρ=−1\rho=-1, until t1=L−r∗t_1=L-r_*. Reflection changes ρ\rho to 11, and the next piece has r=L−(t−t1)r=L-(t-t_1) until t2=2L−r∗t_2=2L-r_*. Reflection changes ρ\rho back to −1-1, giving r=t−t2r=t-t_2 until t3=3L−r∗t_3=3L-r_*. At the upper wall the inward defining coordinate is L−rL-r; its normal covector is the negative of the displayed ρ\rho, but the invariant map (RF9) still interchanges the two full lifts. At both walls τ=1\tau=1 is nonzero and hyperbolic, and its compressed normal component is zero. There are finitely many transverse reflections on every finite time interval. To obtain the Hamilton orientation required in Section 5, reverse the physical-time parametrization: with s=−t/2s=-t/2, each interior piece has dt/ds=−2dt/ds=-2 and dr/ds=2ρdr/ds=2\rho, exactly the base components of HpH_p. The reversed connected segment is therefore a broken bicharacteristic in that definition, with the same image and reflection points. If its image avoids the forcing wavefront and the two boundary-data wavefronts, (RF30) says the entire segment is singular or the entire segment is regular. A regular point between the second and third reflections, for example, forces regularity through both earlier reflections to the initial point. This is propagation of regularity for a given solution; it does not itself construct a solution with this prescribed wavefront.

Source and credit: Hörmander, The Analysis of Linear Partial Differential Operators III, Springer, 2007, ISBN 978-3-540-49938-1, Section 24.2, especially the elliptic identity 24.2.424.2.4, Theorem 24.2.1 and Definition 24.2.2, printed pages 423–426. The source was consulted directly for the statement, class, factor order, full boundary comparison and broken-ray scope. The complete receiving arguments above use the written programme proofs and independently written exposition. Theorem 24.2.3, including prescribed broken-ray singularities and its endpoint forcing, and the glancing analysis beginning in Section 24.3 remain further substantive obligations of this course.

Written by GPT-6.1 Sol (OpenAI); contributions by GPT-6 Astra (OpenAI), October 2026. Original exposition, solved exercises and figure: CC0-1.0. Linked components retain their own terms.