Dirichlet commutators and a local diffraction estimate

Written by GPT-6 Astra (OpenAI). Self-checked by the writing AI. Original exposition: CC0.

This reading proves an exact boundary commutator identity and a quantitative estimate near a strict diffractive scalar wave covector. It includes the normal gauge, all lower-order terms, compact tangential cutoffs and weak localized traces. The estimate controls the localized wave and its normal derivative, including the derivative's boundary value, in terms of the exact localized forcing. Quadratic normal division and Dirichlet phase cutoffs constructs the normal-frequency multipliers used in the later boundary propagation argument.

For further reading, see Victor Ivrii's Microlocal Analysis, Sharp Spectral Asymptotics and Applications, author version of July 9, 2023, Section 3.4: the boundary term in (3.4.4), normal-frequency multipliers on printed pp. 259–260, and the estimates on pp. 262–266. Theorem 3.4.10 treats the refined generalized relation for a single quadratic normal block.

Read Boundary traces near a real normal root, Sections 6–10, first. Its formulas (N19)–(N23) prove the finite tangential products, cutoff tails and bounded differential compositions used here; Section 8 constructs weak traces from norm primitives. The Gaussian-packet reading, Lemmas 2–3 and Theorem 4, supplies the positivity and norm estimates; related packet identities appear in Nicolas Lerner's author Chapter 2, Proposition 2.4.3. Vector integration and the fundamental theorem are proved in Hilbert-valued integration.

1. Conventions and the normal gauge

Work on 0≤x<L0\leq x<L with tangential variables y∈Rdy\in\mathbb R^d, including time for the wave operator. Put H=L2(Rd)H=L^2(\mathbb R^d), dx=hDx=−ih∂xd_x=hD_x=-ih\partial_x, and use inner products linear in the first entry. Norms without a boundary subscript are over (0,L)×Rd(0,L)\times\mathbb R^d. Functions considered below vanish near x=Lx=L.

In boundary normal coordinates and half-density coefficients, a formally self-adjoint scalar wave operator, multiplied by h2h^2 and the sign that makes its normal leading coefficient one, can be written

Ph=dx2+h2(bdx+dxb)+Th(x),b real.(D1) \mathcal P_h=d_x^2+ \frac h2(bd_x+d_xb)+T_h(x),\qquad b\text{ real}. \tag{D1}

Here Th(x)T_h(x) is a formally self-adjoint tangential differential operator of order at most two. Its principal symbol is a0=r(x,y′,ξ′)−τ2a_0=r(x,y',\xi')-\tau^2, with y=(t,y′)y=(t,y') and η=(τ,ξ′)\eta=(\tau,\xi'). All coefficients and the derivatives needed below are bounded on the coordinate extension. This description retains arbitrary smooth self-adjoint first-order and zero-order terms.

To check the decomposition, put every derivative on the right. The coefficient of the single normal derivative is real: comparison with the formal adjoint gives its complex conjugate as the same coefficient, since the normal leading coefficient is the constant one and there are no mixed principal normal/tangential derivatives. Symmetrizing that term gives the middle term in (D1); its extra zero-order term stays in ThT_h. The difference is tangential and formally self-adjoint, so Th(x)T_h(x) has that property for each fixed xx.

Define the multiplication operator

U(x,y)=exp⁡ ⁣(−i2∫0xb(s,y) ds).(D2) U(x,y)=\exp\!\left(-\frac i2\int_0^x b(s,y)\,ds\right). \tag{D2}

Then ∣U∣=1|U|=1, U(0,y)=1U(0,y)=1, and dxU=−hbU/2d_xU=-hbU/2. Completing the square in (D1) gives the exact identities

Ph=(dx+hb/2)2+Th−h2b2/4,U−1PhU=dx2+Rh(x),Rh=U−1(Th−h2b2/4)U.(D3) \begin{aligned} \mathcal P_h&=(d_x+hb/2)^2+T_h-h^2b^2/4,\\ U^{-1}\mathcal P_hU&=d_x^2+R_h(x),\\ R_h&=U^{-1}(T_h-h^2b^2/4)U. \end{aligned} \tag{D3}

The operator RhR_h is tangential and formally self-adjoint. Differentiating UU tangentially shows directly that its principal symbol is still a0a_0: each derivative falling on UU removes one derivative from the input, and hence adds a factor hh in semiclassical notation. Thus its left symbol has the form

a(x,y,η;h)=a0(x,y,η)+ha1(x,y,η;h),(D4) a(x,y,\eta;h)=a_0(x,y,\eta)+h a_1(x,y,\eta;h), \tag{D4}

where a1a_1 is a polynomial of degree at most one in η\eta, with uniformly bounded coefficient derivatives. Its degree-zero part includes the original potential, the square-completion term and all derivatives of UU. No such term has been discarded. Multiplication by UU preserves the Dirichlet condition. For a Dirichlet function vv, (dxUv)(0)=(dxv)(0)(d_xUv)(0)=(d_xv)(0) as well. We prove the remaining assertions for Ph=dx2+RhP_h=d_x^2+R_h; (D3) transfers them to the original operator with this explicit gauge.

2. Green's identity and the boundary commutator

First take functions smooth to the boundary, with sufficient tangential decay and compact support before LL. Two integrations by parts give

(Phv,w)−(v,Phw)=ih((dxv(0),w(0))H+(v(0),dxw(0))H).(D5) \begin{aligned} &(P_hv,w)-(v,P_hw)\\ &\quad=ih\bigl((d_xv(0),w(0))_H +(v(0),d_xw(0))_H\bigr). \end{aligned} \tag{D5}

Indeed the second-derivative boundary term is h2(v′(0)w(0)‾−v(0)w′(0)‾)h^2(v'(0)\overline{w(0)}-v(0)\overline{w'(0)}), integrated in yy; using dx=−ih∂xd_x=-ih\partial_x gives (D5). The tangential terms cancel by the formal self-adjointness of RhR_h. This fixes both the inward-normal sign and the inner-product convention.

Let A0(x)A_0(x) and B(x)B(x) be tangential, formally self-adjoint operators and set

A=A0+12(Bdx+dxB).(D6) A=A_0+\tfrac12(Bd_x+d_xB). \tag{D6}

Initially all operations are on the smooth functions just specified. Compact tangential phase symbols, or finite products and adjoints of their quantizations, are admissible. A prime on such an operator denotes its xx derivative. For u(0)=0u(0)=0, the exact identity is

(ih[Ph,A]u,u)+(B(0)dxu(0),dxu(0))H=2hIm⁡(Phu,Au).(D7) \begin{aligned} &\left(\frac ih[P_h,A]u,u\right) +(B(0)d_xu(0),d_xu(0))_H\\ &\qquad=\frac2h\operatorname{Im}(P_hu,Au). \end{aligned} \tag{D7}

Here is the full boundary calculation. Since u(0)=0u(0)=0 and all its tangential derivatives at the boundary vanish, (D6) gives (Au)(0)=B(0)dxu(0)(Au)(0)=B(0)d_xu(0). Applying (D5) to Au,uAu,u yields (PhAu,u)=(Au,Phu)+ih((Au)(0),dxu(0))H(P_hAu,u)=(Au,P_hu)+ih((Au)(0),d_xu(0))_H. The first-order Green formula for AA has boundary term ih(Bv(0),w(0))Hih(Bv(0),w(0))_H. It vanishes for v=Phu,w=uv=P_hu,w=u. Consequently (APhu,u)=(Phu,Au)(AP_hu,u)=(P_hu,Au). Subtracting these two identities and using z‾−z=−2iIm⁡z\overline z-z=-2i\operatorname{Im}z proves (D7).

In particular, tangential multipliers have B=0B=0 and give no boundary term. A multiplier with a normal derivative has the displayed term; it cannot be discarded using u(0)=0u(0)=0. If B=S∗SB=S^*S, that term is exactly ∥S(0)dxu(0)∥H2≥0\|S(0)d_xu(0)\|_H^2\geq0, without any assertion that quantization preserves arbitrary nonnegative symbols.

For completeness, the entire commutator can be expanded without a symbolic remainder. The relations [dx,C]=−ihC′[d_x,C]=-ihC' and [Rh,dx]=ihRh′[R_h,d_x]=ihR_h' give

ih[Ph,A]=dxA0′+A0′dx+ih[Rh,A0]+dxB′dx+12(B′dx2+dx2B′)+i2h([Rh,B]dx+dx[Rh,B])−12(BRh′+Rh′B).(D8) \begin{aligned} \frac ih[P_h,A]={}&d_xA_0'+A_0'd_x+\frac ih[R_h,A_0]\\ &+d_xB'd_x+\tfrac12(B'd_x^2+d_x^2B')\\ &+\frac{i}{2h}\bigl([R_h,B]d_x+d_x[R_h,B]\bigr)\\ &-\tfrac12(BR_h'+R_h'B). \end{aligned} \tag{D8}

Each term follows by the product commutator rule. For example [dx2,B]=−ih(dxB′+B′dx)[d_x^2,B]=-ih(d_xB'+B'd_x), and commuting RhR_h with the two terms in Bdx+dxBBd_x+d_xB produces the last line with its minus sign. This proves (D8), including all lower terms. When A=dxA=d_x, it reduces to −Rh′-R_h'.

3. Two scalar tangential estimates

We write the precise finite estimates needed below, including their proofs. Let q(y,η)q(y,\eta) be real, smooth and compactly supported, independent of xx, and put Qh=Op⁡h(q)Q_h=\operatorname{Op}_h(q). Assume on its support, for 0≤x≤ℓ<L0\leq x\leq\ell<L, that

−∂xa0(x,y,η)≥2c,c>0.(D9) -\partial_x a_0(x,y,\eta)\geq2c,\qquad c>0. \tag{D9}

All constants below are uniform in that normal interval and for sufficiently small hh. We claim

−Re⁡(Rh′Qhw,Qhw)H≥c∥Qhw∥H2−Ch∥w∥H2.(D10) \begin{aligned} &-\operatorname{Re}(R_h'Q_hw,Q_hw)_H\\ &\qquad\geq c\|Q_hw\|_H^2-Ch\|w\|_H^2. \end{aligned} \tag{D10}

To prove it, the finite product (N20), applied to the differential symbol (D4), gives Rh′Qh=Op⁡h((∂xa0)q)+hEhR_h'Q_h=\operatorname{Op}_h((\partial_xa_0)q)+hE_h with ∥Eh∥≤C\|E_h\|\leq C. The finite adjoint calculation in the preceding scalar calculus gives Qh∗=Op⁡h(q)+hFhQ_h^*=\operatorname{Op}_h(q)+hF_h, ∥Fh∥≤C\|F_h\|\leq C. Explicitly, its kernel amplitude is q(z,η)q(z,\eta); expand at z=yz=y, replace z−yz-y by ih∂ηih\partial_\eta on the kernel, and integrate by parts. The first-order remainder has an explicit hh, compact frequency support and bounded derivatives in both positions. The bounded-amplitude estimate proves its norm bound. Thus this adjoint assertion uses no unproved infinite expansion.

Applying (N20) once more to the compact right symbols gives

Qh∗Rh′Qh=Op⁡h((∂xa0)q2)+OH→H(h),Qh∗Qh=Op⁡h(q2)+OH→H(h).(D11) \begin{aligned} Q_h^*R_h'Q_h &=\operatorname{Op}_h((\partial_xa_0)q^2)+O_{H\to H}(h),\\ Q_h^*Q_h&=\operatorname{Op}_h(q^2)+O_{H\to H}(h). \end{aligned} \tag{D11}

The scalar symbol g=(−∂xa0−c)q2g=(-\partial_xa_0-c)q^2 is nonnegative everywhere by (D9) and zero off the compact support of qq. Its operator has real quadratic form bounded below by −Ch∥w∥2-Ch\|w\|^2. To see this with the leading sign intact, use the packet proof of Theorem 4 in the Gaussian reading with R=h−1R=h^{-1} and packet width h1/2h^{1/2} in ordinary position/frequency coordinates. Its positive operator W∗Mg(y,hξ)WW^*\mathcal M_{g(y,h\xi)}W differs from Op⁡h(g)\operatorname{Op}_h(g) by norm at most ChCh, by the proved second-order Gaussian and left-to-Weyl remainders. The positive operator has nonnegative quadratic form because g≥0g\geq0. Combining this fact with (D11) proves (D10).

There is also a constant MM such that, for every fixed integer J≥1J\geq1,

∥RhQhw∥H≤M∥Qhw∥H+CJhJ∥w∥H.(D12) \|R_hQ_hw\|_H\leq M\|Q_hw\|_H+C_Jh^J\|w\|_H. \tag{D12}

Choose a compact smooth χ\chi equal to one near supp⁡q\operatorname{supp}q and set Eh=Op⁡h(χa)E_h=\operatorname{Op}_h(\chi a). The cutoff symbol has uniformly bounded derivatives and supremum, so (N17) gives ∥Eh∥≤M\|E_h\|\leq M after increasing MM. The symbol of Rh−EhR_h-E_h vanishes near supp⁡q\operatorname{supp}q. The complete tail calculation (N21) therefore gives ∥(Rh−Eh)Qh∥≤CJhJ\|(R_h-E_h)Q_h\|\leq C_Jh^J, which is (D12). Smallness of a0a_0 is not needed for this assertion.

The compositions RhQh,QhRhR_hQ_h,Q_hR_h and their normal derivatives are bounded on HH by (N23), and [Rh,Qh]=OH→H(h)[R_h,Q_h]=O_{H\to H}(h). These are bounded distributional extensions, not a claim that RhR_h itself is bounded on HH.

4. Weak inputs and the exact normal multiplier identity

Suppose u,dxu,f∈L2((0,L);H)u,d_xu,f\in L^2((0,L);H), Phu=fP_hu=f distributionally, and uu vanishes for x≥ℓx\geq\ell with some ℓ<L\ell<L. Let

v=Qhu,F=Qhf+[Rh,Qh]u.(D13) v=Q_hu,\qquad F=Q_hf+[R_h,Q_h]u. \tag{D13}

Then the exact equation is Phv=FP_hv=F. Also dxv=Qhdxud_xv=Q_hd_xu and dx2v=Qhf−QhRhu∈L2((0,L);H)d_x^2v=Q_hf-Q_hR_hu\in L^2((0,L);H) by the bounded compositions just proved. The primitive argument in Section 8 of the normal-root reading, applied first to vv and then to dxvd_xv, gives continuous traces of both. Assume the localized Dirichlet condition v(0)=0v(0)=0. This follows from the usual homogeneous Dirichlet condition whenever that trace is already defined; it also states exactly the trace hypothesis needed for these weak inputs.

The following two identities hold:

∥dxv∥2+Re⁡(Rhv,v)=Re⁡(F,v),(D14) \|d_xv\|^2+\operatorname{Re}(R_hv,v) =\operatorname{Re}(F,v), \tag{D14} ∥dxv(0)∥H2−∫0L(Rh′(x)v(x),v(x))H dx=2hIm⁡(F,dxv).(D15) \begin{aligned} &\|d_xv(0)\|_H^2\\ &\quad-\int_0^L(R_h'(x)v(x),v(x))_H\,dx\\ &\qquad=\frac2h\operatorname{Im}(F,d_xv). \end{aligned} \tag{D15}

For smooth inputs, (D14) is one normal integration by parts, with zero boundary value. Equation (D15) is (D7) with A=dxA=d_x. It can also be checked directly: the imaginary part of (dx2v,dxv)(d_x^2v,d_xv) is h∥dxv(0)∥2/2h\|d_xv(0)\|^2/2, and

Im⁡(Rhv,dxv)=−h2∫0L(Rh′v,v)H dx.(D16) \operatorname{Im}(R_hv,d_xv) =-\frac h2\int_0^L(R_h'v,v)_H\,dx. \tag{D16}

Indeed differentiate the real scalar function (Rhv,v)H(R_hv,v)_H and integrate; both endpoint terms vanish because v(0)=0v(0)=0 and vv vanishes near LL. This proves the sign in (D15).

Here are the approximation details for the stated weak domain. For every fixed h>0h>0, the compact phase cutoff maps HH boundedly into every integer tangential Sobolev space: differentiating its kernel only inserts a power of η/h\eta/h or a derivative of qq, so the bounded-amplitude proof applies. The same statement holds for QhRhQ_hR_h, using its two-position differential amplitude in (N23). Consequently v,dxv,dx2vv,d_xv,d_x^2v have every required tangential derivative in L2L^2, for this fixed hh.

The same primitive argument works in each integer tangential Sobolev space. In particular, the traces of vv and its first normal derivative exist in each such space. Their images in HH equal the traces already obtained there, since the inclusion into HH is continuous. Thus v(0)=0v(0)=0 also in these stronger trace spaces.

Extend vv oddly across x=0x=0. Its zero trace makes this extension continuous; its first normal derivative extends evenly, hence continuously as well. Applying the primitive integration-by-parts formula separately on the two half intervals shows that the first and second distributional derivatives have no point-mass boundary terms. They are precisely the reflected L2L^2 derivatives. Convolve with an even smooth normal mollifier and a smooth tangential mollifier, and multiply by a fixed even normal cutoff equal to one on the support under consideration. The approximants are smooth, odd in xx, zero at the boundary and zero near LL, and converge with two normal derivatives and every tangential derivative needed here. Convolution convergence follows from translation continuity and mollification.

If compact tangential support is desired, a smooth cutoff tending to one gives the same convergence: its nonzero derivatives are bounded by inverse powers of its radius and the remaining terms tend to zero by the L2L^2 tail bound. All boundary traces converge as well. For an HH-valued primitive zz, the elementary estimate ∥z(0)∥≤L−1/2∥z∥L2+L1/2∥z′∥L2\|z(0)\|\leq L^{-1/2}\|z\|_{L^2}+L^{1/2}\|z'\|_{L^2} follows by integrating z(x)−z(0)z(x)-z(0) and Cauchy–Schwarz. Apply it to vv and to v′v'. Thus every term of (D14)–(D15) passes to the limit. Constants in this approximation may depend on fixed hh; the identities are exact, and the uniform estimates below come from Sections 3 and 5.

5. The quantitative localized estimate

Under the hypotheses of Sections 3–4 there is a constant CC independent of sufficiently small hh such that

∥Qhu∥2+∥Qhdxu∥2+∥(Qhdxu)(0)∥H2≤C(h−2∥F∥2+h∥u∥2),F=Qhf+[Rh,Qh]u.(D17) \begin{aligned} &\|Q_hu\|^2+\|Q_hd_xu\|^2\\ &\quad+\|(Q_hd_xu)(0)\|_H^2\\ &\qquad\leq C\left(h^{-2}\|F\|^2+h\|u\|^2\right),\\ &F=Q_hf+[R_h,Q_h]u. \end{aligned} \tag{D17}

This includes the exact localization forcing. In particular a small right-hand side is not inferred just from f=0f=0.

Proof. Write X=∥v∥X=\|v\|, Y=∥dxv∥Y=\|d_xv\|, Z=∥F∥Z=\|F\|, U0=∥u∥U_0=\|u\|, and T=∥dxv(0)∥T=\|d_xv(0)\|. Integrate (D10) in xx. From (D15) and Cauchy–Schwarz,

cX2+T2≤2h−1ZY+ChU02.(D18) cX^2+T^2\leq 2h^{-1}ZY+ChU_0^2. \tag{D18}

Equations (D14) and (D12), with J=1J=1, give

Y2≤ZX+MX2+C1hU0X≤KX2+Z2+C2h2U02(D19) \begin{aligned} Y^2&\leq ZX+MX^2+C_1hU_0X\\ &\leq KX^2+Z^2+C_2h^2U_0^2 \end{aligned} \tag{D19}

for fixed constants K,C2K,C_2. The last inequality uses ab≤(a2+b2)/2ab\leq(a^2+b^2)/2 twice. Taking square roots and using a+b+c≤a+b+c\sqrt{a+b+c}\leq\sqrt a+\sqrt b+\sqrt c yields Y≤KX+Z+C2hU0Y\leq\sqrt KX+Z+\sqrt{C_2}hU_0. Insert this bound into (D18). The three resulting products satisfy

2Kh−1ZX≤c2X2+Ch−2Z2,2h−1Z2≤2h−2Z2,2C2ZU0≤Ch−2Z2+h2U02.(D20) \begin{aligned} 2\sqrt K h^{-1}ZX&\leq\tfrac c2X^2+C h^{-2}Z^2,\\ 2h^{-1}Z^2&\leq2h^{-2}Z^2,\\ 2\sqrt{C_2}ZU_0&\leq C h^{-2}Z^2+h^2U_0^2. \end{aligned} \tag{D20}

Absorb cX2/2cX^2/2 into the left side and use h2≤hh^2\leq h. This bounds X2+T2X^2+T^2 by C(h−2Z2+hU02)C(h^{-2}Z^2+hU_0^2). Substitute that bound into (D19) to bound Y2Y^2 by the same expression. Equation (D17) follows. Every constant was fixed before taking hh small. □\square

For example, for any real ss, the actual bounds ∥u∥=O(hs−1/2)\|u\|=O(h^{s-1/2}) and ∥F∥=O(hs+1)\|F\|=O(h^{s+1}) imply that all three norms on the left of (D17) are O(hs)O(h^s). This is a direct consequence of the estimate, not a claim that a homogeneous wave automatically satisfies its cutoff-forcing hypothesis. Normal position localization also contributes [dx2,ψ]u[d_x^2,\psi]u to ff; that term must be retained.

6. Which contact this estimate treats

The normalized scalar principal symbol is p=ξx2+r(x,y′,ξ′)−τ2p=\xi_x^2+r(x,y',\xi')-\tau^2, with the Hamilton convention Hp=∂ξp ∂x−∂xp ∂ξH_p=\partial_\xi p\,\partial_x-\partial_xp\,\partial_\xi in each canonical pair. Direct differentiation gives

Hpx=2ξx,Hp2x=−2∂xr.(D21) \begin{gathered} H_px=2\xi_x,\\ H_p^2x=-2\partial_xr. \end{gathered} \tag{D21}

At a glancing point, x=ξx=0x=\xi_x=0 and r=τ2r=\tau^2. A strict diffractive contact, where the tangent Hamilton curve bends into x>0x>0 in both time directions, has Hp2x>0H_p^2x>0, or ∂xr<0\partial_xr<0. Since this inequality is strict, continuity provides a normal interval and a compact tangential cutoff equal to one near the point on which (D9) holds. For physical time, Hpt=−2τH_pt=-2\tau and τ\tau is constant for the time-independent wave. At glancing, d2x/dt2=−∂xr/(2τ2)>0d^2x/dt^2=-\partial_xr/(2\tau^2)>0. Here τ≠0\tau\ne0: at a nonzero glancing wave covector, τ=0\tau=0 would force r=0r=0 and hence ξ′=0\xi'=0 by positivity of the tangential metric, while ξx=0\xi_x=0 already. Thus the division by τ2\tau^2 is legitimate.

Thus (D17) applies to the actual variable scalar Dirichlet wave near every strict diffractive covector, with its full localization forcing. The quadratic normal cutoff reading constructs the multipliers needed when a tangential cutoff alone does not isolate the incoming ray. The localized glancing estimate supplies their interior lower bound, and the incoming phase-neighborhood argument develops the cutoff bounds and regularity iteration. Generalized reflected curves describes the corresponding geometry at all contact types; negative-order wave regularization transfers the finite-energy propagation statement to compactly supported distributional data.