A normal gauge that keeps the boundary condition
The weak-source lift removes the natural-dual forcing without discarding its boundary action. The next operation removes the first normal derivative from the interior equation. With matrix coefficients these two operations leave a complex Robin condition. We construct the matrix gauge, regularize that condition exactly, and prove the full normal-commutator identity, including its boundary form and its matrix lower-order terms.
The identity first holds on actual twice weakly differentiable approximants. We then prove its weak graph extension using complementary Sobolev dual pairings, including the normal trace. The weak graph identity alone does not supply higher microlocal estimates, incoming-error bounds or strict-diffractive propagation.
1. The problem and the exact earlier proofs
N0. Conventions and inputs. Work in a flat collar X=[0,∞)x×Rzd, with values in CN, finite N, and Lebesgue density. We use Dj=−i∂j; all inner products are linear in their first entry. Tangential adjoints are actual adjoints in Lz2, with x as a parameter. Coefficients are smooth up to the boundary with all derivatives bounded on the localized collar. After extension they may be taken bounded with all derivatives on the half-space. Every function to which the identity is applied has bounded normal support. Tangential pseudodifferential symbols are classical S1,0m with uniform bounds for the required normal derivatives.
The normalized form retained from WSL:L0 is
q(v,ϕ)=(Dxv,Dxϕ)+a,b∑(habDbv,Daϕ)+j∑(ℓjDjv,ϕ)+j∑(mjv,Djϕ)+(cv,ϕ).(NC1)
Here hab=hba are real scalars; ℓj,mj,c are arbitrary complex matrices. Normal indices are denoted by x, tangential indices by a,b, and j ranges over both. The energy domain is either H01(X;CN) or H1(X;CN). After WSL:L4–L8, the right-hand side is an interior distribution g∈Lx2Hz−1 locally. Its microlocal improvement, when assumed, is exactly the one proved there.
For zero-trace tests the interior equation is
Pv=g,P=Dx2+KDx+T,K=ℓx+mx,T=a,b∑DahabDb+a∑ℓaDa+a∑Da∘Mma+c+MDxmx.(NC2)
Multiplication is written M when confusion with composition is possible. WSL:L6–L7 proves the graph trace of the actual flux Dxv+mxv. Thus the natural domain has zero flux in Hz−1/2, whereas the Dirichlet domain has v∣0=0.
Exact dependencies are collected in the proof map. WSL:L6 supplies graph density and the signed Green formula. SY:H1–H2 supplies Fourier inversion, Plancherel and all real Sobolev spaces. SY:H7–H8 supplies ordered composition and actual adjoints, with finite-seminorm remainders; SY:H10 supplies their Sobolev bounds. BF:A0 supplies the two energy domains and their smooth density. These are complete earlier programme proofs, not external theorem citations.
2. An ordered matrix gauge
N1. Construction and inverse. On the chosen collar solve
2DxS+KS=0,S(0,z)=I,equivalentlySx=−2iKS.(NC3)
Put F=−iK/2. A direct construction is
S(x,z)=I+n≥1∑∫0<tn<⋯<t1<xF(t1,z)⋯F(tn,z)dtn⋯dt1.(NC4)
The order of the matrices is part of the formula. If ∥F∥≤M on 0≤x≤a, the n-th term has norm at most (Ma)n/n!. A fixed number k of tangential derivatives gives at most Cknk ordered terms, each bounded by Cknan/n!. These summable bounds justify termwise differentiation and integration. Normal differentiation removes the outermost integral and gives NC3; repeated normal derivatives then follow from that equation. The same bounds are uniform in every fixed finite collection of parameters.
The right equation Hx=−HF, H(0)=I, has the same construction, with right multiplication and its corresponding order. Differentiation gives (HS)x=0, hence HS=I. A square finite-dimensional matrix with a left inverse is injective and therefore surjective, by elementary elimination; thus H=S−1. This proves invertibility and ∥S∥,∥S−1∥≤eMa. Uniqueness also follows directly: the difference of two solutions is bounded by M∫0x∥W(t)∥dt; iterating this inequality n times gives sup∥W∥(Ma)n/n!, which tends to zero. Extend K to vanish past the collar if a half-space operator is wanted; S is then constant in x there. No self-adjointness or simultaneous diagonalization is used.
N2. Transform both slots of the weak form. Set v=Su, ϕ=S−∗ψ, and q′(u,ψ)=q(Su,S−∗ψ). Multiplication by either smooth bounded invertible matrix is an isomorphism on H1 and on H01, by the product rule, the inverse bounds and smooth density. For the transformed antidual use exactly g′(ψ)=g(S−∗ψ). For an interior distribution this is multiplication g′=S−1g.
Write gxx=1, gxa=gax=0, gab=hab. Expanding both slots, without commuting matrices, gives the same scalar principal form and
ℓj′mi′c′=S−1ℓjS−i∑(DiS−1)gijS,=S−1miS+j∑S−1gij(DjS),=S−1cS+j∑S−1ℓj(DjS)−i∑(DiS−1)miS−i,j∑(DiS−1)gij(DjS).(NC5)
The minus signs follow from (DiS−∗)∗=−DiS−1. Conversely, inserting these three coefficients in the form reproduces every term in the expanded NC1, proving equality on smooth inputs and then on both energy domains by continuity. In particular,
mx′=κ=21S−1(mx−ℓx)S,ℓx′=−κ,P′=S−1PS=Dx2−R,(NC6)
where, as an ordered tangential operator,
R=−S−1(MDx2S+KDxS+T∘MS).(NC7)
Indeed Dx2(Su)=SDx2u+2(DxS)Dxu+(Dx2S)u; NC3 cancels the coefficient of Dxu. The principal symbol of R remains r=−∑habηaηb, times I. These formulas include all normal derivatives of the original lower coefficients.
3. Boundary traces through the gauge
N3. The retained Robin condition. We first record the positive trace as well as the earlier negative trace. Put Λ=(1+∣Dz∣2)1/2. For smooth compactly supported u, integrate the derivative of ∥Λ1/2u(x)∥2, using tangential Plancherel:
∥u(0)∥H1/22≤2∥Λu∥L2(X)∥∂xu∥L2(X).(NC8)
Inward translation, convolution and compact cutoffs give smooth density in H1, as in BF:A0. The inequality defines a continuous trace by completion, agrees with ordinary restriction for smooth functions, and shows that the H01 trace is zero. The graph trace in WSL:L6 is continuous for N∈Lx2Lz2, ∂xN∈Lx2Hz−1, and lies in H−1/2.
Smooth multiplication commutes with both traces. For the graph trace, the product rule bounds ∂x(MN)=(∂xM)N+M∂xN in Lx2Hz−1; SY:H10 bounds multiplication on H−1 and H−1/2. Apply the smooth identity to graph-dense approximants. For the positive trace use H1 density. Therefore the exact flux factorization
Dxv+mxv=S(Dxu+κu)(NC9)
implies, since S(0)=I,
{u∣0=0,(Dxu+κu)∣0=0in Hz−1/2,Dirichlet domain,natural domain.(NC10)
Here Dxu itself has the graph trace: the transformed equation gives Dx2u=Ru+g′∈Lx2Hz−1. The second line is a complex Robin condition. An interior normal gauge has not made it the condition Dxu∣0=0.
For later use, the same density argument in H2 applies to each first derivative. Its traces lie in H1/2 by NC8. Tangential differentiation commutes with trace on smooth approximants; Fourier weights then show u∣0∈H3/2, because the sum of its H1/2 norm and those of its tangential first derivatives is equivalent to its H3/2 norm. Also Dxu∣0∈H1/2. Density in H2 follows by inward translation first, then ordinary convolution at a smaller radius and cutoffs, with convergence of every derivative of order at most two. No boundary condition is required for this density statement.
4. Smoothing with an exact boundary correction
N4. Localization and the Gaussian commutators. Choose ω(x,z)=χx(x)χz(z) of compact support, equal to one on a smaller collar, with ∂xω(0,z)=0. Let w=ωu and H=P′u. Its equation is exactly
H0=P′w=ωH+2(Dxω)Dxu+(Dx2ω)u−[R,ω]u.(NC11)
All terms belong to Lx2Hz−1. The last commutator has tangential order one and bounded coefficients. The boundary condition NC10 holds for w with the same κ0=κ(0,z). Extend the localized coefficients with bounded derivatives. No term in NC11 is dropped.
For 0<ϵ≤1 let Jϵ=e−ϵ2∣Dz∣2, a scalar multiplier independent of x, and jϵ(η)=e−ϵ2∣η∣2. Plancherel and dominated convergence give strong convergence to identity in every fixed Hzt, and also in the corresponding Lx2Hzt spaces. Every derivative of the Gaussian has the form ϵ∣α∣pα(ϵη)e−ϵ2∣η∣2. Splitting at ϵ⟨η⟩=1 proves
jϵ is bounded in S0,jϵ−1=O(ϵ1/2) in S1/2.(NC12)
For derivative order zero, use ∣e−t2−1∣≤min(t2,1). For a positive derivative order, the displayed polynomial-Gaussian formula gives the same bound after multiplication by ⟨η⟩∣α∣−1/2: in the first region its positive powers are bounded by a power of ϵ⟨η⟩, and in the second region Gaussian decay bounds every remaining power. These estimates cover every symbol seminorm.
In the ordered composition formula SY:H7, the zero-order product cancels in a commutator with a scalar multiplier. The remainder with one derivative, including its full finite-seminorm estimate and actual operator interpretation SY:H8, therefore gives
[κ0,Jϵ]=O(ϵ1/2) in Ψ−1/2,[R,Jϵ]=O(ϵ1/2) in Ψ3/2.(NC13)
Use jϵ−1, whose commutator is unchanged, to obtain the rates. The ordinary global bounded-coefficient calculus applies to the multiplier itself; proper support is not silently imposed on a Gaussian. In particular SY:H10 proves
∥[κ0,Jϵ]b∥H1≤Cϵ1/2∥b∥H1/2,∥[R,Jϵ]w∥Lx2Hz−1≤Cϵ1/2∥w∥Lx2Hz1/2.(NC14)
N5. Actual H2 approximants and their full equations. In the Robin case put βϵ=[κ0,Jϵ]γw, where γw=w∣0. Choose a fixed smooth compactly supported normal function θ, equal to one near zero, and define
Cϵ(x,z)=−ixθ(x)βϵ(z),wϵ=Jϵw+Cϵ.(NC15)
For the Dirichlet case set Cϵ=0 and use the same Jϵw. For each positive ϵ, tangential smoothing gives Jϵw∈Lx2Hzt and DxJϵw∈Lx2Hzt for every finite t. Moreover
Dx2Jϵw=Jϵ(Rw+H0)∈L2.
Thus all distributional derivatives of total order at most two lie in L2, which is the definition of H2. The same is true of Cϵ: both Jϵγw and Jϵ(κ0γw) have every finite tangential Sobolev order, and multiplication by κ0 preserves them.
The graph trace commutes with Jϵ by its bounded action and graph-density approximation. Hence the Robin defect of Jϵw is precisely βϵ. Since Cϵ∣0=0 and DxCϵ∣0=−βϵ,
(Dx+κ0)wϵ∣0=0,wϵ⟶w in H1(X).(NC16)
The first assertion is exact; the second follows from NC14, NC8 and ∥Cϵ∥H1≤Cθ∥βϵ∥H1. In the Dirichlet case the trace stays zero without a correction. The complete equation is
P′wϵ=JϵH0−[R,Jϵ]w+P′Cϵ⟶H0 in Lx2Hz−1.(NC17)
Indeed NC14 treats the commutator. The normal derivatives of xθ are fixed bounded compactly supported functions, and an order-two tangential operator maps Hz1 to Hz−1; these facts bound P′Cϵ by C∥βϵ∥H1 in the stated graph space. Each individual right-hand side in NC17 is in L2, since wϵ∈H2. No uniform L2 bound as ϵ↓0 has been asserted.
5. The full normal commutator
N6. An exact identity including the matrix lower terms. In this section write P=Dx2−R for the transformed operator. Decompose R=Rh+iRa, with Rh=(R+R∗)/2 and Ra=(R−R∗)/(2i). Both are tangentially self-adjoint; Rh has order two and scalar real principal symbol rI, and Ra has order at most one.
Choose scalar, componentwise operators A=A∗∈Ψ−1, C=C∗∈Ψ0, smoothly depending on x. A dot means ∂x. Set
B=C−2iA˙,Q=ADx+B=21(ADx+DxA)+C,S0=21(ARh+RhA),J0=21(ARh−RhA).(NC18)
Thus Q is formally self-adjoint in the interior, S0∗=S0, and J0∗=−J0. We do not assume that Q preserves a boundary domain. Define the Hermitian operator array
E11E10E00=2A˙,=B˙+iJ0−ARa,E01=E10∗,=21(AR˙h+R˙hA)+i[C,Rh]−(B∗Ra+RaB).(NC19)
Its entries have orders at most 1−j−k. For u∈H2 with bounded normal support, write uj=Dxju in the interior and for its trace at the boundary. Then
2Im(Pu,Qu)XB(u)=B(u)+j,k=0∑1(Ejkuk,uj)X,=(Au1,u1)∂+2Re(Bu0,u1)∂+(S0u0,u0)∂.(NC20)
All boundary operators are evaluated at x=0.
Here is the full calculation. Prove it first for smooth inputs compactly supported in the closed half-space. The signed Green identity is
(f,Dxh)X=(Dxf,h)X−i(f(0),h(0))∂.
Consequently 2Im(h,Qh)X=−(Ah(0),h(0))∂. Apply Green once to (Dx2u,Qu), use [Dx,Q]=−i(A˙Dx+B˙), and then apply this last formula to h=Dxu. The resulting boundary terms are (Au1,u1)+2Re(Bu0,u1); the volume terms are
2(A˙u1,u1)+2Re(B˙u,u1).
For the real part of R, move A tangentially and split ARh=S0+J0. Direct differentiation gives
−2Im(S0u,u1)=−∂x(S0u,u)+(S˙0u,u),−2Im(J0u,u1)=2Re(iJ0u,u1).(NC21)
The first identity uses u1=−iux and the self-adjointness of S0. Integration gives +(S0u0,u0)∂. Since B=C−iA˙/2, the remaining term from Rh is
(i[C,Rh]u,u)−((A˙Rh+RhA˙)u/2,u).
Its second summand cancels precisely the two derivatives of A in S˙0, leaving the first term of E00.
Finally the part iRa of R contributes
−2Re(Rau,Au1+Bu)=−2Re(ARau,u1)−((B∗Ra+RaB)u,u).(NC22)
This proves NC19–NC20 with every ordering shown. For general H2 inputs, the H2 density and traces proved in N3 and the Sobolev bounds SY:H10 make all boundary and volume expressions continuous. In particular u0∈H3/2, u1∣0∈H1/2, and Q:H2→H1. Cutoffs at tangential infinity followed by the same density prove the identity without tangential compact support. No normal trace of an arbitrary H1 derivative has been paired in this extension.
N7. Principal polynomial and the surviving boundary form. Let a=σ−1(A), c0=σ0(C), k=σ1(Ra), so k=k∗ is a matrix. For the tangential Poisson bracket use {f,g}z=fηgz−fzgη, with contraction of indices. SY:H7 applied to NC19 gives
e11e10=e01e00=2axI,=(∂xc0+21{a,r}z)I−ak,=(a∂xr+{c0,r}z)I−2c0k.(NC23)
Indeed i[A,Rh] has principal symbol {a,r}zI; the scalar leading symbols commute with every matrix lower coefficient. The term −iA¨/2 in B˙ is one order below the cross coefficient. All other terms follow by their ordered leading products. With p=ξ2−r and q=aξ+c0, summing the three coefficients yields
j,k′=0∑1ejk′ξj+k′=HpqI−2qk.(NC24)
To check the sign directly, Hpq=2axξ2+(2∂xc0+{a,r}z)ξ+a∂xr+{c0,r}z. The Hermitian matrix k cannot in general be replaced by a scalar or omitted.
In the Dirichlet case the boundary form is exactly (Au1,u1)∂. Choosing A=−T∗T with T∈Ψ−1/2 gives −∥Tu1∥2, including all lower terms. In the Robin case u1=−κ0u0, so instead
B(u)=((κ0∗Aκ0−κ0∗B−B∗κ0+S0)u0,u0)∂.(NC25)
This operator is self-adjoint of order at most one. There is no unconditional sign assertion for NC25. Each wϵ from N5 satisfies the appropriate version exactly. N8 below supplies the weak graph extension of this identity. Uniform estimates at higher microlocal orders are still needed for strict propagation.
N8. The identity on the actual weak graph. Suppose u∈H1(X) has bounded normal support and Pu∈Lx2Hz−1. Coefficients retain all the bounds above. Put b0=γu∈Hz1/2 and b1=γDxu∈Hz−1/2; the latter is defined because Dx2u=Ru+Pu∈Lx2Hz−1. For an operator of order m, SY:H10 gives every required tangential mapping. In particular,
Q:H1(X)⟶Lx2Hz1,∥Qu∥Lx2Hz1≤C∥u∥H1(X),A:Hz−1/2→Hz1/2,B:Hz1/2→Hz1/2,S0:Hz1/2→Hz−1/2.(NC26)
The first bound follows separately from ADxu and Bu: the normal derivative is in L2, and A gains one tangential order. Thus the source pairing and every term of the boundary form have a continuous definition:
2Im⟨Pu,Qu⟩Lx2(H−1,H1)Bweak(u)=Bweak(u)+j,k=0∑1(EjkDxku,Dxju)X,=⟨Ab1,b1⟩+2Re⟨Bb0,b1⟩+⟨S0b0,b0⟩.(NC27)
Every bracket is the linear-first extension of the L2 pairing between the indicated complementary Sobolev spaces. For example its frequency integral is bounded by the product of the H1/2 and H−1/2 norms. SY:H1–H2 proves that this agrees with the distribution pairing. No product of two arbitrary negative-order distributions is being defined.
To prove NC27, take the approximants of N5 with w=u. Its construction and estimates use precisely these global graph norms; the localization NC11 is needed only when starting from a local solution. If no boundary condition is specified, use Jϵu alone. If one is specified, retain the exact Dirichlet or Robin correction. In all cases,
∥uϵ−u∥H1+∥P(uϵ−u)∥Lx2Hz−1⟶0,∥γDx(uϵ−u)∥Hz−1/2⟶0.(NC28)
For the second conclusion, the two graph norms of Dx(uϵ−u) tend to zero: its normal derivative is i(R(uϵ−u)+P(uϵ−u)), and R:Hz1→Hz−1. WSL:L6 gives the trace convergence. NC8 gives convergence of the solution traces in H1/2. All three boundary pairings now converge by NC26. For the volume terms, E00 of order one maps Hz1 to Lz2; E10,E01 of order zero and E11 of order minus one give continuity on H1(X) by Cauchy–Schwarz. The source pairing converges by the first bound in NC26 and the first limit in NC28. Passing to the limit in the already proved NC20 establishes NC27.
For a Dirichlet solution and A=−T∗T, the remaining boundary term is exactly −∥Tb1∥Lz22: T:H−1/2→L2, and adjoint duality follows by Sobolev density from SY:H2 and SY:H8. For a Robin solution, b1=−κ0b0 in H−1/2, and NC25 is valid as the H−1/2,H1/2 pairing of its order-one operator. This proves the identity in both weak domains with the actual complex boundary coefficient. It does not bound that Robin form by a favorable sign, nor supply the higher microlocal estimate.
6. Three solved checks of the mechanism
Exercise 1. A nonunitary matrix gauge
Take E=(0010), K=2iE, mx=M=(i01−2) and ℓx=K−M. Compute the gauge and the Robin matrix.
Solution. Since E2=0, the series NC4 stops at its linear term: S=I+xE, S−1=I−xE. Substitution gives 2DxS+KS=−2iE+2iE=0. Direct ordered multiplication gives
κ=S−1(M−iE)S=(i01−i+x(i+2)−2),ℓx′=−κ,mx′=κ.(NC29)
The off-diagonal coefficient changes with x. Replacing these products by commuting scalars would lose that change. At the boundary κ0=M−iE, and the second line of NC10 keeps its full off-diagonal entry.
Exercise 2. Smoothing a Robin mode
On the circle with normalized measure take κ0(z)=eiz and a boundary value b(z)=einz, n∈Z. Find the defect of Gaussian smoothing and cancel it. The periodic Fourier proof is the same finite-mode calculation as SY:H1; no periodic pseudodifferential theorem is needed in this example.
Solution. Multiplication raises the frequency by one. Thus
[κ0,Jϵ]b=(e−ϵ2n2−e−ϵ2(n+1)2)ei(n+1)z.(NC30)
It is generally nonzero. Multiplying the right-hand side by −ixθ(x) produces zero value and normal derivative equal to the negative of this defect at x=0, exactly as in NC15. For example w=θ1(x)(1−ixκ0(z))einz, with θ1=1 near zero, has trace b, derivative Dxw∣0=−κ0b, and hence the required original Robin condition. This also verifies that the computed defect occurs for an actual smooth input, not just a formal boundary value.
Exercise 3. Check the normal orientation and the imaginary coefficient
Suppress z, take P=Dx2−(r+ib), A=−1, C=c0, with real constants, and let u(x)=eitxf(x), where t∈R and f is any real smooth rapidly decreasing profile. Find both sides of NC20 explicitly.
Solution. Here B=c0, S0=−r, J0=0, E10=E01=b, E00=−2c0b, and E11=0. Since Dxu=eitx(tf−if′),
B(u)=(−t2+2c0t−r)f(0)2−f′(0)2,∑(Ejkuk,uj)=2b(t−c0)∫0∞f2dx.(NC31)
For an independent direct evaluation, remove the common factor eitx. Then Pu=(t2−r)f−f′′−i(2tf′+bf), and Qu=(c0−t)f+if′. The imaginary part of their scalar inner product is
(t2−2c0t+r)ff′+f′′f′+b(t−c0)f2.
Integrating 2ff′=(f2)′ and 2f′′f′=(f′2)′ gives NC31 and NC20. In particular the sign at x=0 is fixed, not a choice of energy orientation. At t=0, the profiles e−x, xe−x, (1+x)e−x give boundary forms −1−r, −1, −r, respectively. Their squared L2 norms are 1/2,1/4,5/4, by integrating by parts: the zeroth moment is 1/2, and each successive moment is j/2 times the preceding one. The middle profile is Dirichlet; the last has zero normal derivative.
Left: the exact real and imaginary parts of the off-diagonal coefficient in NC29. Middle: the exact defect coefficient in NC30 for n=2, with its sign fixed by the order of the commutator. Right: the three boundary forms from Exercise 3 at r=2,t=0; the volume contribution is zero when b=0. These are exact coefficient and profile calculations, not numerical solutions of a diffraction problem. The reproducible figure is tied to N1–N8 and the three solutions.
7. Credit, validation and the remaining estimate
Lars Hörmander, The Analysis of Linear Partial Differential Operators III: Pseudo-Differential Operators, corrected 2007 printing, ISBN 978-3-540-49938-1, Section 24.4, pp. 443–444 and 449, is the source for the scalar normal-gauge reduction, normal energy identity and the use of regularization. The present proof writes out the ordered matrix gauge, the transformed weak form, its complex Robin condition, the boundary-preserving Gaussian correction, the exact Hermitian operator array and its weak graph extension. It does not attribute the full matrix Robin argument to a scalar Dirichlet statement in the book. The earlier source-lift proof and Hilbert coefficient calculus retain their own source and component notices.
The prerequisite proof index gives the earlier results used here. Independent new exposition, code and illustration are CC0-1.0. Linked components retain their own licences.
The uniform commutant estimate for the full weak problem requires bounds for the source and every regularization defect, the Robin boundary contribution, characteristic-square control of energy regularity, incoming terms and a justified higher-order limit on a common microlocal neighborhood. The characteristic-symbol, weak-quadratic and source-lift supplements provide the earlier pieces.