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≥0, with D=−i∂, finite coefficient space CN, 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)
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 v belongs locally to Xt for some finite real t. On one fixed open conic parameter region W, with η=0, assume
χPv∈Lx2Hzm(m∈R),χ(0)γv∈Hzm(m∈R)(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 a, r(a)=0 and rx(a)>0. A short negative-normal interior characteristic germ from a is ordinarily regular for v, and its intervening projection stays inside W.
We prove that one open parameter neighborhood W0 of a satisfies
χv∈Xm(m∈R),γDxv∈H∞ locally on W0∩{x=0}.(RC3)
As throughout the programme, H∞ 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 t; 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∈Xs and Pv∈L2Hs−2 on a localized patch, the equation gives
Dx2v=Rv+Pv∈L2Hs−2,γv∈Hs−1/2,γDxv∈Hs−3/2.(RC4)
The last two statements follow from the proved unequal-order trace inequality LG:L1, applied first with orders s,s−1, and then with s−1,s−2. That proof uses Fourier truncation and the absolutely continuous normal representative. It holds for every real s. 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 w compact in the positive normal direction define
Ew(x)={w(x),3w(−x)−2w(−2x),x≥0,x<0.(RC5)
The coefficients make both the value and first derivative agree at zero: 3−2=1 and −3+4=1. Apply this first to tangentially truncated smooth approximants. For j=0,1,2, a change of variable bounds the negative-half norm of DxjEw, in any fixed tangential Sobolev space, by the sum of the positive-half norms with factors 3 and 2j+1−1/2. In particular the extension is bounded for the three norms
∥w∥L2Hs,∥Dxw∥L2Hs−1,∥Dx2w∥L2Hs−2.(RC6)
This passes to rough inputs: after applying Λs−2, the input is H2 in the normal variable with values in Lz2. 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 H2 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=⟨(ρ,η)⟩ and λ=⟨η⟩, uses the exact identity
Q4λ2s−4=λ2s+2ρ2λ2s−2+ρ4λ2s−4.(RC7)
Hence w∈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), they are
γ(Av)=A(0)γv,γDx(Av)=A(0)γDxv−iAx(0)γv.(RC8)
Indeed Dx(Av)=ADxv−iAxv distributionally. Tangential Fourier truncation, the parameter Sobolev bounds and the two trace limits of (RC4) prove (RC8). Equivalently v,Dxv 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>0. The geometric tubes, full matrix roots A±, and transported matrix test Q− from ISE:I1–I3 depend on the operator and support, not on v's initial Sobolev order. Use exactly those full constructions, including 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)
For v∈Xs, z,h∈L2Hb. The complete right side is ISE's expression
g=Q−Pv+Q−(Pe−P)v−Q−Sv+[Dx−A−,Q−]z.(RC10)
All its terms have every tangential L2Hm order. The first uses (RC2) on a larger support; the other coefficients are full tangentially smoothing families acting on v,Dxv at their known finite orders. Separated input errors have the same property. Since A− has order one, Dxh∈L2Hb−1, and LG:L1 gives the actual trace h(x1)∈Hb−1/2. The uncut factor z has a trace too: its equation has finite order L2Hs−2, including the nonlocalized Pe−P. Formula (RC8) identifies h(x1)=Q−(x1)z(x1).
The exact all-real Sobolev evolution satisfies
h(x)=E−(x,x1)h(x1)+i∫x1xE−(x,y)g(y)dy.(RC11)
To justify this for the actual h, it has an absolutely continuous Hb−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)∗ϕ. Its derivative is iA−∗ϕx. The derivative of the linear-first pairing vanishes, since h′=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 b, with the indices displayed rather than an assumed L2 input.
The source integral is continuous in every tangential Hm. The homogeneous part has the full matrix Cauchy wavefront graph ρ=−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) 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)
on each entire selected incoming support. No strong normal regularity for the source is required.
R3. Retain the complete equation when the original input is rough. Suppose all tests in an open region put v in Xs, for any real s. Choose the two SDC commutants with UMD's sufficiently large damping constant, and a scalar χ equal to one on a larger neighborhood of all their supports. Define
Wϵ=1+ϵ2∣Dz∣2Λs,vϵ=Wϵχv,Pvϵ=gϵ+Zϵvϵ,Zϵ=[P,Wϵ]Wϵ−1,(RC13)
where, exactly,
gϵ=WϵχPv+Wϵ[P,χ]v,[P,χ]v=−2iχxDxv−χxxv−[R,χ]v.(RC14)
For fixed ϵ>0, the weight has order s−2. Thus vϵ∈L2H2, Dxvϵ∈L2H1, and gϵ∈L2H1 locally after outer cutoffs. The exact equation gives Dx2vϵ∈L2. So this actual input is H2, regardless of the sign of s.
Uniformly as ϵ↓0, the same symbol bounds give
∥vϵ∥X0+∥Pvϵ∥L2H−1≤Cs,∥γDxvϵ∥H−3/2≤Cs.(RC15)
The last bound is LG:L1; the value trace is bounded in every Hm by (RC2), (RC8) and the uniform order-s weight. The scalar symbol of Wϵ makes Zϵ uniformly order one with Hermitian part order zero, by the full UMD:D6 remainder proof for every real s.
Every localized source pairing in LG:L3 remains bounded. In fact J∗Wϵ[P,χ] has tangentially smoothing coefficients on both v and Dxv, because the full χ−I symbol is smoothing on the support of J. Its arbitrary-order bound now acts on the finite norms L2Ht,L2Ht−1 instead of H1,L2. The number of integrations can be increased by ∣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ϵ∈H2, not to an unproved energy realization of v. UMD's ordered matrix damping and the two scalar completions then give precisely LG20. The completion Green error pairs the uniformly H−3/2 normal trace with the smooth value trace. Choosing A(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-s coefficients act on the assumed Xs input. We obtain
∥Kvϵ∥X1/2+α=1∑2∥TαγDxvϵ∥2≤Cs.(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χv in X0. 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)
This proves the half-gain for every finite starting order. It does not apply a weak form with an inadmissible test to v.
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<0, assume χv∈Xs locally. The complete cutoff equation (RC14) has forcing L2Hs−1. Use DDB:B2's full elliptic auxiliary cylinder with the stable Dirichlet evaluation isomorphism. Its extension error has no Dx2 term and is tangentially smoothing on the finite-order input. R1 puts Av∈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)
The last map uses its proved gain of one ordinary order and Q≥λ. Thus one use of the actual identity yields Av∈H(2,s−1)⊂Xs+1, for every real s. 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+v has tangentially regular forcing and actual smooth initial value Q+(0)γv. The uniqueness argument of (RC11), now started at its finite negative trace order and with the positive root, gives Q+v∈CxHm for every m. The equation gives its normal derivative in L2Hm−1. Full tangential division proves all Xm 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 s, its order-minus-two parametrix maps the cutoff source L2Hs−1 into L2Hs+1, with the derivative in L2Hs. 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λ, for every real q,
Q≍λ,λq+∣ρ∣λq−1≤CqQq.(RC19)
Consequently ordinary Hq regularity on that finite-normal cone yields both Xq components after the full low-normal cutoff. Terms in which Dx 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≥1 in DDB:B4 by the explicit cone comparison (RC19). Conversely Xq implies ordinary Hq on any smaller finite-normal cone by Qq≍λq, for every real q.
6. Iterate over one open region, starting where the distribution lies
R5. The finite-order Dirichlet theorem. Choose the accessible open region W0 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 W0. Those are geometric facts about the continuous broken cap map, proved there for all three minimum cases.
Let P(s) mean that every scalar tangential test compactly supported in W0 puts v in Xs. The finite-order hypothesis gives P(t). Suppose P(s) holds. At every glancing boundary point, R2 supplies the whole incoming supports after the damping constant for s 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 W0 gives ordinary Hs+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 Dxv, whose finite order is known. Thus
P(s)⟹P(s+1/2),P(t+j/2)(j=0,1,2,…).(RC20)
Order monotonicity proves the first assertion in (RC3). For its normal trace, choose any boundary test with support in W0, enlarge it once within W0, and use the equation together with the established arbitrarily high tangential orders. Formula (RC4) then gives every Hm 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.
R6. A comparison distribution with verified hypotheses. Now let u∈H1, Pu=f∈L2, and suppose every tangential localization of f has all L2Hm 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)
No selfadjointness is imposed on M. Assume a regular negative-normal incoming germ for u. 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/2. Choose nested compact boundary cutoffs so that the retained datum hc equals h microlocally on an open boundary region V 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 L, define
e=E+Dhc=K+Lhc,v=u−e.(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 e on the output patch. In particular, on compact subpatches there is a finite m0 such that e,Dxe,Dx2e∈CxHz−m0. This also follows from ADT:D7's finite-loss bounds for the fixed collection of derivatives j≤2. Hence e∈X−m0, and v 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)
Here S 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, 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+ consists of the positive-normal branch and its glancing limits. On a sufficiently short normalized collar, rx>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 e is ordinarily regular on a small negative-normal cap, for every compact distributional hc. The original regular germ makes v regular there as well.
Shrink the common parameter source region so that its boundary portion lies in V. Equations (RC23) give exactly (RC2), and the finite starting order and incoming cap were just verified. R5 therefore gives all tangential mixed orders for v and smooth actual traces on one fixed smaller boundary region. This comparison uses only the actual hc; 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)
Every term is its actual distributional trace. In particular, multiplication by M acts on the value trace in its displayed order. ABI:B4–B7 proves the localized two-sided conormal inverse Q=JQL, including the full M, with no symmetry or commutativity requirement. Its receiving estimate, for nested tests B0,B1, is
∥B0hc∥Hm≤Cm(∥B1B+hc∥Hm−2/3+∥hc∥H−A).(RC25)
Choose any fixed finite A 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 m. Thus (RC24) makes hc, and hence the actual h, microlocally smooth on one fixed smaller boundary region.
AWT:W5's uniform boundary-collar assertion now applies to K+Lhc. 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 e is therefore tangentially smooth there, with every normal derivative, despite other singularities of the retained boundary datum. Together with R5 for v, this proves
χu∈Xm(m∈R) on one fixed neighborhood of the strict anchor.(RC26)
Since u 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+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=H1, with scalar real wave principal symbol, all complex lower terms including derivatives on tests, and the actual source in V∗. Suppose its smooth source front is absent on an open compressed region. BPL:L7–L8 gives a boundary-preserving actual H1 representative U with PU=F∈L2, 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)
The statement holds on one fixed smaller parameter region. The actual value and normal traces of Z are tangentially smooth there, and Z is regular in every original compressed H1 order on one fixed cone. Thus W∈H1 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 W. Adding Z, undoing the ordered gauge, and undoing BPL and the coordinate transformations gives absence of the original smooth H1-based boundary front at the strict anchor.
The other germ follows by conjugating the entire normalized equation. With D=−i∂, conjugation changes D to −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−Mu=−β.(RC28)
The natural-dual functional conjugates by applying it to the conjugated test and conjugating the value. The H1 domain is preserved. An outgoing regular germ at tangential covector η∗ becomes the negative-normal incoming germ at −η∗. The proved theorem applies with Robin matrix −M. Returning gives the opposite-germ conclusion without a selfadjointness assumption.
Let FV=WFbV,∞(u) on the original source-regular region. The present result for V=H1, together with DDB:B7 for V=H01, proves at a strict diffractive point q, with its unique characteristic lift q,
q∈FV⟹π(exp(tHp)q)∈FV(∣t∣<ε).(RC29)
Choose ε 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 q 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+cx2 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).
Solution. The polynomial is
3w(−x)−2w(−2x)=a+bx−5cx2(x<0).(RC30)
Its value and first derivative match a,b, so the first and second distributional derivatives have no boundary delta terms. The second derivative can jump from −10c to 2c, which is allowed in L2 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; the first derivative has no delta because the values agree, but the second derivative has the term 2bδ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. 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
−37,−611,−34,−65,−31,61,32,67.(RC31)
Seven gains reach 7/6≥1, hence X1=H1 after any compact test supported in W0. Six gains reach only 2/3. Every implication in (RC20) concerns all points of that same W0; 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−Dy1Dy2, λ>0, μ dual to y1, and b=−μλ−1/3. Let F be the proved zero-free outgoing Airy function, φ=F′/F, and let M be a constant two-by-two matrix with M2=0. Compute the conormal row of
e(x,η)=F(b)F(−xλ2/3+b)h(η).(RC32)
Solution. Its value is h, and Dxe∣0=iλ2/3φ(b)h=:nh. Hence the full row is (nI+M)h, and, since n=0,
(nI+M)−1=n−1I−n−2M.(RC33)
Both ordered products are identity by M2=0. Dropping M leaves the nonzero inverse error n−1M. For a solution with smooth β, the comparison gives (nI+M)h=β−Cv; the inverse is applied to that complete row. Conjugation changes it as in (RC28), not merely by conjugating n while leaving the Robin coefficient fixed.
11. The geometry and the proved comparison
F0. Coordinates and meanings. The left panel is the exact (t,x) projection of the affine tangent orbit in Exercise 3 at μ=0,λ=1, with Hamilton parameter s: x=s2, normal covector ρ=s, and physical time t=y1+y2=−s−2s3/3. The omitted spatial coordinate is y1−y2=−s+2s3/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/4, x=1/16, t=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.