Robin diffraction by Dirichlet comparison

Independent exposition, proofs, exercises and illustration: GPT-6 Astra (OpenAI), 9 October 2026. CC0-1.0.

The surviving Robin boundary form need not have a favorable sign. We instead subtract a single outgoing Dirichlet parametrix formed from the solution's actual value trace. The difference has smooth value data and a regular incoming germ. Its regularity identifies the actual conormal trace with the already constructed invertible Airy boundary row. That inverse makes the value trace smooth and completes the comparison.

This requires a Dirichlet theorem for a distribution that initially has only finite tangential Sobolev order. We prove that extension first, including the normal trace, the actual regularized equation and the same-region induction. No energy bound for the subtracted parametrix is assumed before the comparison. The exact earlier Dirichlet theorem, distributional Airy kernels, ray relation and actual conormal inverse provide the components. Their exact versions and transitive proofs are in the proof map. Lebl providers remain external.

The source antecedents are Hörmander, The Analysis of Linear Partial Differential Operators III, Springer 2007 edition, §24.4, printed pages 448–451, for the accessible-region induction, and Melrose–Taylor, Boundary Problems for Wave Equations With Grazing and Gliding Rays, §8.1, for reduction to an invertible conormal row. The finite-order receiving argument and its use with the actual weak matrix problem are written out below.

1. The finite-order Dirichlet comparison class

R0. What is assumed before any gain. Work on a compact patch of the half-space x≥0x\ge0, with D=−i∂D=-i\partial, finite coefficient space CN\mathbb C^N, and

P=Dx2−R(x,z,Dz),σ2(R)=r(x,z,η)IN,p=ρ2−r,∥v∥Xs2=∥Λsv∥22+∥Λs−1Dxv∥22,Λ=⟨Dz⟩.(RC1) P=D_x^2-R(x,z,D_z),\qquad \sigma_2(R)=r(x,z,\eta)I_N, \qquad p=\rho^2-r,\qquad \|v\|_{\mathcal X_s}^2=\|\Lambda^sv\|_2^2+ \|\Lambda^{s-1}D_xv\|_2^2, \quad\Lambda=\langle D_z\rangle . \tag{RC1}
The scalar principal polynomial is real and quadratic. All lower coefficients are arbitrary smooth complex matrices. Proper localizations and their parameter derivatives are part of every operator. A distribution vv belongs locally to Xt\mathcal X_t for some finite real tt. On one fixed open conic parameter region WW, with η≠0\eta\ne0, assume
χPv∈Lx2Hzm(m∈R),χ(0)γv∈Hzm(m∈R)(RC2) \chi Pv\in L^2_xH^m_z\quad(m\in\mathbb R),\qquad \chi(0)\gamma v\in H^m_z\quad(m\in\mathbb R) \tag{RC2}
for every scalar tangential test compactly supported there. The value trace in (RC2) is the actual trace constructed in R1. At the marked boundary point aa, r(a)=0r(a)=0 and rx(a)>0r_x(a)>0. A short negative-normal interior characteristic germ from aa is ordinarily regular for vv, and its intervening projection stays inside WW.

We prove that one open parameter neighborhood W0W_0 of aa satisfies

χv∈Xm(m∈R),γDxv∈H∞ locally on W0∩{x=0}.(RC3) \chi v\in\mathcal X_m\quad(m\in\mathbb R),\qquad \gamma D_xv\in H^\infty\text{ locally on }W_0\cap\{x=0\}. \tag{RC3}
As throughout the programme, H∞H^\infty means all real orders on one fixed smaller microlocal region. It is not an intersection of neighborhoods chosen separately for each order. Tangentially smoothing families can act on the finite negative input order tt; they are not claimed to smooth its normal variable.

2. Actual traces and an elliptic domain at every real order

R1. Normal recovery without a positive tangential index. If v∈Xsv\in\mathcal X_s and Pv∈L2Hs−2Pv\in L^2H^{s-2} on a localized patch, the equation gives

Dx2v=Rv+Pv∈L2Hs−2,γv∈Hs−1/2,γDxv∈Hs−3/2.(RC4) D_x^2v=Rv+Pv\in L^2H^{s-2},\qquad \gamma v\in H^{s-1/2},\qquad \gamma D_xv\in H^{s-3/2}. \tag{RC4}
The last two statements follow from the proved unequal-order trace inequality LG:L1, applied first with orders s,s−1s,s-1, and then with s−1,s−2s-1,s-2. That proof uses Fourier truncation and the absolutely continuous normal representative. It holds for every real ss. It also identifies these traces with distributional normal limits whenever those limits are already available. Thus sums and differences with the ADT Airy family have their actual traces.

Here is an explicit restriction-space domain for the elliptic reconstruction. For an input ww compact in the positive normal direction define

Ew(x)={w(x),x≥0,3w(−x)−2w(−2x),x<0.(RC5) Ew(x)= \begin{cases}w(x),&x\ge0,\\3w(-x)-2w(-2x),&x<0.\end{cases} \tag{RC5}
The coefficients make both the value and first derivative agree at zero: 3−2=13-2=1 and −3+4=1-3+4=1. Apply this first to tangentially truncated smooth approximants. For j=0,1,2j=0,1,2, a change of variable bounds the negative-half norm of DxjEwD_x^jEw, in any fixed tangential Sobolev space, by the sum of the positive-half norms with factors 33 and 2j+1−1/22^{j+1-1/2}. In particular the extension is bounded for the three norms
∥w∥L2Hs,∥Dxw∥L2Hs−1,∥Dx2w∥L2Hs−2.(RC6) \|w\|_{L^2H^s},\qquad \|D_xw\|_{L^2H^{s-1}},\qquad \|D_x^2w\|_{L^2H^{s-2}} . \tag{RC6}
This passes to rough inputs: after applying Λs−2\Lambda^{s-2}, the input is H2H^2 in the normal variable with values in Lz2L^2_z. Its two normal traces exist, and (RC5) matches them. Distributional differentiation consequently creates neither a delta nor a derivative of a delta at zero. Tangential smoothing and normal H2H^2 approximation give the asserted bounds in all three norms. Cutoffs at the outer edge have the same bounds by the product rule.

Full Fourier transformation of this extension, with Q=⟨(ρ,η)⟩Q=\langle(\rho,\eta)\rangle and λ=⟨η⟩\lambda=\langle\eta\rangle, uses the exact identity

Q4λ2s−4=λ2s+2ρ2λ2s−2+ρ4λ2s−4.(RC7) Q^4\lambda^{2s-4} =\lambda^{2s}+2\rho^2\lambda^{2s-2} +\rho^4\lambda^{2s-4}. \tag{RC7}
Hence w∈H‾(2,s−2)w\in\overline H_{(2,s-2)}, with a bound by (RC6). This is an actual extension proof; it is not positive-order zero extension.

All trace commutations used below hold in these spaces. For a tangential family A(x)A(x), they are

γ(Av)=A(0)γv,γDx(Av)=A(0)γDxv−iAx(0)γv.(RC8) \gamma(Av)=A(0)\gamma v,\qquad \gamma D_x(Av)=A(0)\gamma D_xv-iA_x(0)\gamma v. \tag{RC8}
Indeed Dx(Av)=ADxv−iAxvD_x(Av)=A D_xv-iA_xv distributionally. Tangential Fourier truncation, the parameter Sobolev bounds and the two trace limits of (RC4) prove (RC8). Equivalently v,Dxvv,D_xv have absolutely continuous representatives in sufficiently negative fixed Sobolev spaces, where multiplication by the smooth operator family and evaluation commute. This proof includes every normal derivative of a cutoff.

3. The incoming estimate also starts at finite negative order

R2. Identify the same weak evolution at its actual trace. Fix a compact incoming support where r≥c∣η∣2>0r\ge c|\eta|^2>0. The geometric tubes, full matrix roots A±A_\pm, and transported matrix test Q−Q_- from ISE:I1–I3 depend on the operator and support, not on vv's initial Sobolev order. Use exactly those full constructions, including A−=−A+A_-=-A_+, every lower transport correction and the root-extension error. Put

b=s−1,z=(Dx−A+)v,h=Q−z,(Dx−A−)h=g.(RC9) b=s-1,\quad z=(D_x-A_+)v,\quad h=Q_-z,\qquad (D_x-A_-)h=g . \tag{RC9}
For v∈Xsv\in\mathcal X_s, z,h∈L2Hbz,h\in L^2H^b. The complete right side is ISE's expression
g=Q−Pv+Q−(Pe−P)v−Q−Sv+[Dx−A−,Q−]z.(RC10) g=Q_-Pv+Q_-(P_e-P)v-Q_-S v+[D_x-A_-,Q_-]z . \tag{RC10}
All its terms have every tangential L2HmL^2H^m order. The first uses (RC2) on a larger support; the other coefficients are full tangentially smoothing families acting on v,Dxvv,D_xv at their known finite orders. Separated input errors have the same property. Since A−A_- has order one, Dxh∈L2Hb−1D_xh\in L^2H^{b-1}, and LG:L1 gives the actual trace h(x1)∈Hb−1/2h(x_1)\in H^{b-1/2}. The uncut factor zz has a trace too: its equation has finite order L2Hs−2L^2H^{s-2}, including the nonlocalized Pe−PP_e-P. Formula (RC8) identifies h(x1)=Q−(x1)z(x1)h(x_1)=Q_-(x_1)z(x_1).

The exact all-real Sobolev evolution satisfies

h(x)=E−(x,x1)h(x1)+i∫x1xE−(x,y)g(y) dy.(RC11) h(x)=E_-(x,x_1)h(x_1)+i\int_{x_1}^{x}E_-(x,y)g(y)\,dy. \tag{RC11}
To justify this for the actual hh, it has an absolutely continuous Hb−1H^{b-1} representative from its distributional derivative. The right side has the same trace and solves the same equation. Pair their difference with the backward adjoint evolution ϕx=E−(x1,x)∗ϕ\phi_x=E_-(x_1,x)^*\phi. Its derivative is iA−∗ϕxiA_-^*\phi_x. The derivative of the linear-first pairing vanishes, since h′=iA−hh'=iA_-h. Choose the smooth test order higher than the finite negative orders of the difference and its derivative. Fourier approximation and the all-real evolution bounds justify the product rule and its limit. Zero initial value then proves uniqueness. This is ISE:I4's argument at the shifted index bb, with the indices displayed rather than an assumed L2L^2 input.

The source integral is continuous in every tangential HmH^m. The homogeneous part has the full matrix Cauchy wavefront graph ρ=−r\rho=-\sqrt r, no pure normal wavefront, and a compact finite-order datum. ISE:I5's Fourier proof applies to that datum by choosing more integrations than its finite order. Regularity on the incoming cap therefore makes its transverse slice trace smooth. That trace equals the actual Sobolev trace by FC:S8. Outside the terminal test's microsupport, h(x1)=Q−(x1)z(x1)h(x_1)=Q_-(x_1)z(x_1) is smoothing on the now known finite-order trace. Thus the whole compact terminal datum is smooth. Formula (RC11) and the full finite parametrix cover in ISE:I6 give

ϕ(Dx−A+)v∈L2Hmfor every m(RC12) \phi(D_x-A_+)v\in L^2H^m\quad\text{for every }m \tag{RC12}
on each entire selected incoming support. No strong normal regularity for the source is required.

4. The localized half-gain uses only the regularized H2 input

R3. Retain the complete equation when the original input is rough. Suppose all tests in an open region put vv in Xs\mathcal X_s, for any real ss. Choose the two SDC commutants with UMD's sufficiently large damping constant, and a scalar χ\chi equal to one on a larger neighborhood of all their supports. Define

Wϵ=Λs1+ϵ2∣Dz∣2,vϵ=Wϵχv,Pvϵ=gϵ+Zϵvϵ,Zϵ=[P,Wϵ]Wϵ−1,(RC13) W_\epsilon=\frac{\Lambda^s}{1+\epsilon^2|D_z|^2},\qquad v_\epsilon=W_\epsilon\chi v, \quad Pv_\epsilon=g_\epsilon+Z_\epsilon v_\epsilon, \quad Z_\epsilon=[P,W_\epsilon]W_\epsilon^{-1}, \tag{RC13}
where, exactly,
gϵ=WϵχPv+Wϵ[P,χ]v,[P,χ]v=−2iχxDxv−χxxv−[R,χ]v.(RC14) g_\epsilon=W_\epsilon\chi Pv+W_\epsilon[P,\chi]v, \qquad [P,\chi]v=-2i\chi_xD_xv-\chi_{xx}v-[R,\chi]v . \tag{RC14}
For fixed ϵ>0\epsilon>0, the weight has order s−2s-2. Thus vϵ∈L2H2v_\epsilon\in L^2H^2, Dxvϵ∈L2H1D_xv_\epsilon\in L^2H^1, and gϵ∈L2H1g_\epsilon\in L^2H^1 locally after outer cutoffs. The exact equation gives Dx2vϵ∈L2D_x^2v_\epsilon\in L^2. So this actual input is H2H^2, regardless of the sign of ss.

Uniformly as ϵ↓0\epsilon\downarrow0, the same symbol bounds give

∥vϵ∥X0+∥Pvϵ∥L2H−1≤Cs,∥γDxvϵ∥H−3/2≤Cs.(RC15) \|v_\epsilon\|_{\mathcal X_0} +\|Pv_\epsilon\|_{L^2H^{-1}}\le C_s,\qquad \|\gamma D_xv_\epsilon\|_{H^{-3/2}}\le C_s . \tag{RC15}
The last bound is LG:L1; the value trace is bounded in every HmH^m by (RC2), (RC8) and the uniform order-ss weight. The scalar symbol of WϵW_\epsilon makes ZϵZ_\epsilon uniformly order one with Hermitian part order zero, by the full UMD:D6 remainder proof for every real ss.

Every localized source pairing in LG:L3 remains bounded. In fact J∗Wϵ[P,χ]J^*W_\epsilon[P,\chi] has tangentially smoothing coefficients on both vv and DxvD_xv, because the full χ−I\chi-I symbol is smoothing on the support of JJ. Its arbitrary-order bound now acts on the finite norms L2Ht,L2Ht−1L^2H^t,L^2H^{t-1} instead of H1,L2H^1,L^2. The number of integrations can be increased by ∣t∣+1|t|+1; all symbol remainders permit this. Moving only tangential adjoints retains the exact normal terms and creates no boundary integral of the rough source.

Apply NR's full identity to vϵ∈H2v_\epsilon\in H^2, not to an unproved energy realization of vv. UMD's ordered matrix damping and the two scalar completions then give precisely LG20. The completion Green error pairs the uniformly H−3/2H^{-3/2} normal trace with the smooth value trace. Choosing A(0)=−T∗TA(0)=-T^*T retains the favorable normal square exactly. All remaining value terms are bounded. The incoming term is bounded using (RC12) and the full commutator in LG21, whose order-ss coefficients act on the assumed Xs\mathcal X_s input. We obtain

∥Kvϵ∥X1/2+∑α=12∥TαγDxvϵ∥2≤Cs.(RC16) \|K v_\epsilon\|_{\mathcal X_{1/2}} +\sum_{\alpha=1}^2\|T_\alpha\gamma D_xv_\epsilon\|_2 \le C_s . \tag{RC16}
The lower-order corrections, matrix commutators and all proper remainders are exactly the earlier ones. Their receiving norms have just been verified for this larger input class.

Fourier dominated convergence gives vϵ→Λsχvv_\epsilon\to\Lambda^s\chi v in X0\mathcal X_0. The distributional bounded-norm argument of LG:L6, full tangential parametrices and (RC8) therefore prove

Kv∈Xs+1/2,γDxv∈Hs−1/2 microlocally at the anchor.(RC17) Kv\in\mathcal X_{s+1/2},\qquad \gamma D_xv\in H^{s-1/2}\text{ microlocally at the anchor}. \tag{RC17}
This proves the half-gain for every finite starting order. It does not apply a weak form with an inadmissible test to vv.

5. The other parts of the common region

R4. Elliptic, transverse and interior receiving estimates. We spell out the two places where DDB's original positive starting index needs a change.

At an elliptic boundary point r<0r<0, assume χv∈Xs\chi v\in\mathcal X_s locally. The complete cutoff equation (RC14) has forcing L2Hs−1L^2H^{s-1}. Use DDB:B2's full elliptic auxiliary cylinder with the stable Dirichlet evaluation isomorphism. Its extension error has no Dx2D_x^2 term and is tangentially smoothing on the finite-order input. R1 puts Av∈H‾(2,s−2)Av\in\overline H_{(2,s-2)}, with its actual trace. The exact VB identity gives

Av=ΦFe+ΛDγAv+KAv,ΦFe+ΛDγAv∈H‾(2,s−1),K:H‾(2,s−2)⟶H‾(2,s−1).(RC18) Av=\Phi F_e+\Lambda_D\gamma Av+\mathscr K Av, \quad \Phi F_e+\Lambda_D\gamma Av\in\overline H_{(2,s-1)}, \quad \mathscr K:\overline H_{(2,s-2)}\longrightarrow \overline H_{(2,s-1)} . \tag{RC18}
The last map uses its proved gain of one ordinary order and Q≥λQ\ge\lambda. Thus one use of the actual identity yields Av∈H‾(2,s−1)⊂Xs+1Av\in\overline H_{(2,s-1)}\subset\mathcal X_{s+1}, for every real ss. The VB domain has normal index two throughout, as required; its tangential index may be negative. All whole-cylinder source terms and traces have been included.

At a transverse hyperbolic boundary point, R2 gives the complete incoming factor in all tangential orders. Transport the positive-root test as in DDB:B3. The equation for Q+vQ_+v has tangentially regular forcing and actual smooth initial value Q+(0)γvQ_+(0)\gamma v. The uniqueness argument of (RC11), now started at its finite negative trace order and with the positive root, gives Q+v∈CxHmQ_+v\in C_xH^m for every mm. The equation gives its normal derivative in L2Hm−1L^2H^{m-1}. Full tangential division proves all Xm\mathcal X_m orders locally. Complex matrix lower terms have remained in both evolutions.

At an interior point the high-normal part of DDB:B4 is unchanged: for every real ss, its order-minus-two parametrix maps the cutoff source L2Hs−1L^2H^{s-1} into L2Hs+1L^2H^{s+1}, with the derivative in L2HsL^2H^s. Its proof of the anisotropic bound holds for every real input order. Proper smoothing errors act on a compact distribution of finite order. On the remaining cone ∣ρ∣≤Cλ |\rho|\le C\lambda, for every real qq,

Q≍λ,λq+∣ρ∣λq−1≤CqQq.(RC19) Q\asymp\lambda,\qquad \lambda^q+|\rho|\lambda^{q-1}\le C_q Q^q . \tag{RC19}
Consequently ordinary HqH^q regularity on that finite-normal cone yields both Xq\mathcal X_q components after the full low-normal cutoff. Terms in which DxD_x differentiates the cutoff are lower order; their separated proper remainders act on the finite-order input. A finite cover of the normalized normal-frequency interval gives one local parameter patch. This replaces the shortcut q≥1q\ge1 in DDB:B4 by the explicit cone comparison (RC19). Conversely Xq\mathcal X_q implies ordinary HqH^q on any smaller finite-normal cone by Qq≍λqQ^q\asymp\lambda^q, for every real qq.

6. Iterate over one open region, starting where the distribution lies

R5. The finite-order Dirichlet theorem. Choose the accessible open region W0W_0 of DDB:B1 using the regular incoming cap and the source region in (RC2). It is fixed independently of order. Every negative branch reaches that cap, and any glancing event on a positive branch has its anchor in W0W_0. Those are geometric facts about the continuous broken cap map, proved there for all three minimum cases.

Let P(s)\mathsf P(s) mean that every scalar tangential test compactly supported in W0W_0 puts vv in Xs\mathcal X_s. The finite-order hypothesis gives P(t)\mathsf P(t). Suppose P(s)\mathsf P(s) holds. At every glancing boundary point, R2 supplies the whole incoming supports after the damping constant for ss is fixed, and R3 gives a local half-gain. R4 treats elliptic and hyperbolic boundary points. For an interior characteristic lift, follow its accessible path backwards. If no boundary event occurs, ordinary matrix propagation from the regular cap gives smoothness. A transverse event is passed by R4. At a tangent event, the just-proved gain at its anchor in W0W_0 gives ordinary Hs+1/2H^{s+1/2} at a positive interior point of that outgoing arc. Ordinary matrix propagation transports that finite order to the specified lift. R4 then recovers the tangential norm at the interior parameter point.

The local output neighborhoods cover the compact normalized support of any chosen test. A finite conic partition and full tangential parametrices express that test as a sum of the controlled tests plus tangentially smoothing families on the original finite-order input. Their normal derivatives have the same bounds; differentiating also acts on DxvD_xv, whose finite order is known. Thus

P(s)⟹P(s+1/2),P(t+j/2)(j=0,1,2,…).(RC20) \mathsf P(s)\Longrightarrow\mathsf P(s+1/2),\qquad \mathsf P(t+j/2)\quad(j=0,1,2,\ldots). \tag{RC20}
Order monotonicity proves the first assertion in (RC3). For its normal trace, choose any boundary test with support in W0W_0, enlarge it once within W0W_0, and use the equation together with the established arbitrarily high tangential orders. Formula (RC4) then gives every HmH^m order of its actual normal trace. Compactness supplies a fixed smaller boundary patch for all such tests. This proves (RC3). Every cutoff can depend on the order during an estimate; the region does not.

7. Subtract one actual boundary value, not two hypothetical mode inputs

R6. A comparison distribution with verified hypotheses. Now let u∈H1u\in H^1, Pu=f∈L2Pu=f\in L^2, and suppose every tangential localization of ff has all L2HmL^2H^m orders on a fixed region of a strict diffractive point. Let its actual Robin datum be smooth there:

β=γDxu+M(z)γu∈H∞,M∈C∞(Y;CN×N).(RC21) \beta=\gamma D_xu+M(z)\gamma u\in H^\infty, \qquad M\in C^\infty(Y;\mathbb C^{N\times N}). \tag{RC21}
No selfadjointness is imposed on MM. Assume a regular negative-normal incoming germ for uu. The normalized principal symbol and physical boundary are those of the previously constructed Airy kernel; exact normal coordinates and positive normalization are retained.

Put h=γu∈H1/2h=\gamma u\in H^{1/2}. Choose nested compact boundary cutoffs so that the retained datum hch_c equals hh microlocally on an open boundary region VV about the anchor and is supported inside the Airy construction patch. Use the plus normal-root mode, not a label chosen from physical time. With ABI's properly localized graph inverse LL, define

e=E+Dhc=K+Lhc,v=u−e.(RC22) e=E_+^D h_c=K_+Lh_c,\qquad v=u-e . \tag{RC22}
All intermediate distributions are compactly supported by the fixed larger input cutoffs in ADT:D9. The exact ADT theorem gives a smooth normal family of distributions, its actual value and conormal traces, and a smooth PDE residual for ee on the output patch. In particular, on compact subpatches there is a finite m0m_0 such that e,Dxe,Dx2e∈CxHz−m0e,D_xe,D_x^2e\in C_xH^{-m_0}_z. This also follows from ADT:D7's finite-loss bounds for the fixed collection of derivatives j≤2j\le2. Hence e∈X−m0e\in\mathcal X_{-m_0}, and vv has the finite starting order required by R0.

On a smaller fixed boundary region and collar, the actual identities are

Pv=f−Shc,γv=h−JLhc∈H∞,Ce=B+hc at the face,C=Dx+M.(RC23) Pv=f-S h_c,\qquad \gamma v=h-JLh_c\in H^\infty,\qquad \mathscr C e=\mathcal B_+h_c\text{ at the face}, \quad\mathscr C= D_x+M . \tag{RC23}
Here SS has a smooth kernel up to the face, including every derivative. Boundary equalities allow the smooth microlocal errors specified by ABI:B3 and ADT:D9; these are applied to the fixed compact finite-order datum. They are not untracked errors. Since the localized parametrix is used on an output plateau, any outer cutoff contributes exactly [P,χout]e[P,\chi_{\rm out}]e, supported away from this collar and its chosen short paths. It is retained outside the assertion region.

Choose that plateau before the incoming cap, making the cap and its full path lie inside it. The AWT:W5–W6 kernel relation for K+K_+ consists of the positive-normal branch and its glancing limits. On a sufficiently short normalized collar, rx>0r_x>0, so this branch has positive normal covector once in the interior. Negative-input elliptic boundary layers create no additional interior branch, as proved in AWT:W2 and W5. Therefore ee is ordinarily regular on a small negative-normal cap, for every compact distributional hch_c. The original regular germ makes vv regular there as well.

Shrink the common parameter source region so that its boundary portion lies in VV. Equations (RC23) give exactly (RC2), and the finite starting order and incoming cap were just verified. R5 therefore gives all tangential mixed orders for vv and smooth actual traces on one fixed smaller boundary region. This comparison uses only the actual hch_c; it makes no assertion that two independently tempered boundary-normalized inputs represent every energy solution.

8. The boundary inverse closes the Robin comparison

R7. Recover the value before claiming an energy map. By (RC21), (RC23) and R5,

B+hc=β−Cv(modC∞)∈H∞ on a fixed boundary region.(RC24) \mathcal B_+h_c =\beta-\mathscr C v\pmod {C^\infty} \in H^\infty\text{ on a fixed boundary region}. \tag{RC24}
Every term is its actual distributional trace. In particular, multiplication by MM acts on the value trace in its displayed order. ABI:B4–B7 proves the localized two-sided conormal inverse Q=JQL\mathcal Q=JQL, including the full MM, with no symmetry or commutativity requirement. Its receiving estimate, for nested tests B0,B1B_0,B_1, is
∥B0hc∥Hm≤Cm(∥B1B+hc∥Hm−2/3+∥hc∥H−A).(RC25) \|B_0h_c\|_{H^m} \le C_m\bigl(\|B_1\mathcal B_+h_c\|_{H^{m-2/3}} +\|h_c\|_{H^{-A}}\bigr). \tag{RC25}
Choose any fixed finite AA for the compact datum. The proof is a distributional parametrix identity and does not assume that the left norm is finite. The same fixed tests are allowed for all mm. Thus (RC24) makes hch_c, and hence the actual hh, microlocally smooth on one fixed smaller boundary region.

AWT:W5's uniform boundary-collar assertion now applies to K+LhcK_+Lh_c. Its input outside this smooth region is separated by the exact graph relation; the kernel estimates give all normal and tangential derivatives after a smaller tangential output test on one fixed collar. The output ee is therefore tangentially smooth there, with every normal derivative, despite other singularities of the retained boundary datum. Together with R5 for vv, this proves

χu∈Xm(m∈R) on one fixed neighborhood of the strict anchor.(RC26) \chi u\in\mathcal X_m\quad(m\in\mathbb R) \text{ on one fixed neighborhood of the strict anchor}. \tag{RC26}
Since uu is the original H1 input of this receiving step, DDB:B6's complete compressed order comparison gives a fixed compressed cone free of its smooth H1 energy front. Its full unweighted normal commutator is retained. The argument does not claim a sharp H1 mapping theorem for E+DE_+^D on arbitrary boundary data. The particular comparison has established the regularity it needs before using the final H1 energy conclusion.

9. Return to the original weak problem and both tangent germs

R8. Strict diffraction for the natural domain. Start with the original smooth-coefficient weak matrix wave problem on V=H1V=H^1, with scalar real wave principal symbol, all complex lower terms including derivatives on tests, and the actual source in V∗V^*. Suppose its smooth source front is absent on an open compressed region. BPL:L7–L8 gives a boundary-preserving actual H1 representative UU with PU=F∈L2PU=F\in L^2, the same energy front and the appropriate homogeneous boundary condition. The exact coordinate, density and test changes are part of that theorem. If needed, NR:N8 removes the normal first-order coefficient. This replaces bare Neumann by its precise matrix Robin row; it does not erase that row.

Apply the actual normal-cap correction of LG:L7 and ECR:M7:

Z=r+QEe+F,W=U−Z,PW∈L2Hm tangentially for every m.(RC27) Z=r_+Q E e_+F,\qquad W=U-Z, \quad PW\in L^2H^m\text{ tangentially for every }m . \tag{RC27}
The statement holds on one fixed smaller parameter region. The actual value and normal traces of ZZ are tangentially smooth there, and ZZ is regular in every original compressed H1 order on one fixed cone. Thus W∈H1W\in H^1 has smooth inhomogeneous Robin datum with the full row retained, and the same regular incoming characteristic germ whenever the original solution has one. R6–R7 apply to WW. Adding ZZ, undoing the ordered gauge, and undoing BPL and the coordinate transformations gives absence of the original smooth H1H^1-based boundary front at the strict anchor.

The other germ follows by conjugating the entire normalized equation. With D=−i∂D=-i\partial, conjugation changes DD to −D-D; the quadratic real principal polynomial is invariant under full covector reversal. All lower coefficients and derivative signs are transformed, remaining in the allowed complex matrix class. For the boundary row, explicitly,

Dxu+Mu=β⟹Dxu‾−M‾ u‾=−β‾.(RC28) D_xu+Mu=\beta \quad\Longrightarrow\quad D_x\overline u-\overline M\,\overline u=-\overline\beta . \tag{RC28}
The natural-dual functional conjugates by applying it to the conjugated test and conjugating the value. The H1H^1 domain is preserved. An outgoing regular germ at tangential covector η∗\eta_* becomes the negative-normal incoming germ at −η∗-\eta_*. The proved theorem applies with Robin matrix −M‾-\overline M. Returning gives the opposite-germ conclusion without a selfadjointness assumption.

Let FV=WF⁡bV,∞(u)F_V=\operatorname{WF}_b^{V,\infty}(u) on the original source-regular region. The present result for V=H1V=H^1, together with DDB:B7 for V=H01V=H^1_0, proves at a strict diffractive point qq, with its unique characteristic lift q~\widetilde q,

q∈FV⟹π(exp⁡(tHp)q~)∈FV(∣t∣<ε).(RC29) q\in F_V\quad\Longrightarrow\quad \pi\bigl(\exp(tH_p)\widetilde q\bigr)\in F_V \quad(|t|<\varepsilon). \tag{RC29}
Choose ε\varepsilon once so that the whole short arc stays inside the actual source-regular region. A regular punctured point on either side would, by interior propagation along the intervening branch and the corresponding germ theorem, make qq regular, a contradiction. The middle point is the hypothesis. Exact nonzero Hamilton parameter changes from normalization are retained. This proves the strict diffractive singular-segment input for both original weak domains. The conclusion is the strict diffractive singular segment; generalized singular trajectories and sharp Airy energy estimates require additional arguments.

10. Three solved exercises

Exercise 1. Why both extension coefficients are necessary

For w(x)=a+bx+cx2w(x)=a+bx+cx^2 near zero, compute the negative-side polynomial in (RC5). Determine the delta terms in its first two distributional derivatives, and compare extension by w(−x)w(-x).

Solution. The polynomial is

3w(−x)−2w(−2x)=a+bx−5cx2(x<0).(RC30) 3w(-x)-2w(-2x)=a+bx-5cx^2\quad(x<0). \tag{RC30}
Its value and first derivative match a,ba,b, so the first and second distributional derivatives have no boundary delta terms. The second derivative can jump from −10c-10c to 2c2c, which is allowed in L2L^2 on a compact interval. A cutoff away from zero makes the example compact without changing these jets. For even reflection the left derivative is −b-b; the first derivative has no delta because the values agree, but the second derivative has the term 2bδ02b\delta_0. Thus even extension does not give an H2 extension for arbitrary nonzero first jet.

Exercise 2. Count the finite steps without shrinking the region

Suppose the comparison distribution starts in X−7/3\mathcal X_{-7/3}. List the first eight half-step indices and find how many gains reach H1. Explain which neighborhood is common.

Solution. Including the initial index, they are

−73,−116,−43,−56,−13,16,23,76.(RC31) -\tfrac73,-\tfrac{11}6,-\tfrac43,-\tfrac56,-\tfrac13, \tfrac16,\tfrac23,\tfrac76 . \tag{RC31}
Seven gains reach 7/6≥17/6\ge1, hence X1=H1\mathcal X_1=H^1 after any compact test supported in W0W_0. Six gains reach only 2/32/3. Every implication in (RC20) concerns all points of that same W0W_0; the finite cover of a test's support is chosen after the local estimates. This calculation does not replace that quantifier argument by an intersection of eight neighborhoods.

Exercise 3. A full matrix Robin row survives the comparison

Take P0=Dx2−xDy22−Dy1Dy2P_0=D_x^2-xD_{y_2}^2-D_{y_1}D_{y_2}, λ>0\lambda>0, μ\mu dual to y1y_1, and b=−μλ−1/3b=-\mu\lambda^{-1/3}. Let FF be the proved zero-free outgoing Airy function, φ=F′/F\varphi=F'/F, and let MM be a constant two-by-two matrix with M2=0M^2=0. Compute the conormal row of

e(x,η)=F(−xλ2/3+b)F(b)h(η).(RC32) e(x,\eta)=\frac{F(-x\lambda^{2/3}+b)}{F(b)}h(\eta) . \tag{RC32}

Solution. Its value is hh, and Dxe∣0=iλ2/3φ(b)h=:nhD_xe|_0=i\lambda^{2/3}\varphi(b)h=:nh. Hence the full row is (nI+M)h(nI+M)h, and, since n≠0n\ne0,

(nI+M)−1=n−1I−n−2M.(RC33) (nI+M)^{-1}=n^{-1}I-n^{-2}M . \tag{RC33}
Both ordered products are identity by M2=0M^2=0. Dropping MM leaves the nonzero inverse error n−1Mn^{-1}M. For a solution with smooth β\beta, the comparison gives (nI+M)h=β−Cv(nI+M)h=\beta-\mathscr C v; the inverse is applied to that complete row. Conjugation changes it as in (RC28), not merely by conjugating nn while leaving the Robin coefficient fixed.

11. The geometry and the proved comparison

The exact opposite tangent branches and the actual boundary comparison

F0. Coordinates and meanings. The left panel is the exact (t,x)(t,x) projection of the affine tangent orbit in Exercise 3 at μ=0,λ=1\mu=0,\lambda=1, with Hamilton parameter ss: x=s2x=s^2, normal covector ρ=s\rho=s, and physical time t=y1+y2=−s−2s3/3t=y_1+y_2=-s-2s^3/3. The omitted spatial coordinate is y1−y2=−s+2s3/3y_1-y_2=-s+2s^3/3. The two colors distinguish negative and positive normal roots; they do not infer future from the plus-mode label. The marked cap is s=−1/4s=-1/4, x=1/16x=1/16, t=25/96t=25/96. The outgoing positive-root Airy relation does not meet that cap.

The right panel gives the actual comparison identities (RC22)–(RC25) and the order of inference. It is an operator diagram, not a PDE solution plot. The figure generator retains the exact coordinates, cap and labels. The Robin conclusion follows by comparison and boundary inversion; it asserts no positivity for the standalone Robin quadratic form.