AN-04 proof edition · CC0; exact source credit and separately licensed prerequisites

When conormal regularity and duality force smoothness

This is the complete smoothness argument of AN03-U034, Global boundary operators, compressed wave fronts, and normal extension, Section 7.1 and Sections 8.1–8.6 through (SP15). Original author credit: Codex, September 2026, CC0. Current prerequisite connections and the final identification of the supported representative: AN-04 course-writing task and OpenAI Codex, 5 October 2026, CC0.

The conormal amplitude and Besov proofs supply (C2) and (C16). The supported test spaces, approximation, full dual pairing and intrinsic trace supply (GA1)–(GA2), (GD1)–(GD12) and their actual topologies. The full Fourier Sobolev embedding, Hölder inequality supply the elementary analytic entries. The exact interpolation identity (SP11) is also proved directly below.

The mathematical antecedent is the approved Hörmander III, 2007 eBook, Section 18.3. The same distribution is used throughout the argument.

1. The weighted norm with its complete order shift

Localize in a compact boundary product chart. An element u∈Amu\in\mathcal A^m has, modulo a smooth term, the exact normal conormal oscillatory representation

u(x′,t)=(2π)−1∫Reitτa(x′,τ) dτ,∣Dx′β∂τja(x′,τ)∣≤Cβj⟨τ⟩μ−j,μ=m+n−24.(HS1) u(x',t)=(2\pi)^{-1}\int_{\mathbb R} e^{it\tau}a(x',\tau)\,d\tau, \qquad |D_{x'}^\beta\partial_\tau^j a(x',\tau)| \le C_{\beta j}\langle\tau\rangle^{\mu-j}, \quad \mu=m+\frac{n-2}{4}. \tag{HS1}

This is the original codimension-one shift, including the ambient dimension. The inverse coefficient and D=−i∂D=-i\partial convention are retained. For a normal derivative of order a≥0a\ge0, the amplitude is τaDx′βa\tau^aD_{x'}^\beta a, of order μ+a\mu+a. Split its integral at ∣τ∣=t−1|\tau|=t^{-1}, 0<t<10<t<1. The low-frequency absolute integral is bounded by Ct−μ−a−1Ct^{-\mu-a-1} when μ+a>−1\mu+a>-1, by C(1+∣log⁡t∣)C(1+|\log t|) at equality, and by CC below it. On the high-frequency part integrate by parts in τ\tau N>μ+a+1N>\mu+a+1 times, including the derivatives of a smooth cutoff at ∣τ∣=t−1|\tau|=t^{-1}. Each resulting term is bounded by Ct−μ−a−1Ct^{-\mu-a-1}. The same estimate holds for every tangential derivative, uniformly on a smaller compact chart. Thus the complete safe weighted L2L^2 implication is

∫K0∫0ct2ν∣D′βDnau(x′,t)∣2 dt dx′<∞whenever ν≥0,ν>m+n4+a.(HS2) \int_{K_0}\int_0^c t^{2\nu}|D'^\beta D_n^a u(x',t)|^2 \,dt\,dx'<\infty \quad\text{whenever } \nu\ge0,\quad \nu>m+\frac n4+a. \tag{HS2}

The strict inequality covers the logarithmic case. Equation (HS2) states the weight inside the norm as tνt^\nu; the power inside the squared integral is 2ν2\nu.

2. Smoothness from the dual conormal condition

2.1. Every boundary delta jet and its precise conormal order

For a compact smooth tangential density h(x′)h(x'), the normal delta derivative has the original Fourier representation

h(x′)⊗δ(j)(t)=(2π)−1∫eitτ(iτ)jh(x′) dτ,mj=2−n4+j=12−n4+j.(SP1) h(x')\otimes\delta^{(j)}(t) =(2\pi)^{-1}\int e^{it\tau}(i\tau)^j h(x')\,d\tau, \qquad m_j=\frac{2-n}{4}+j=\frac12-\frac n4+j. \tag{SP1}

The exponent mjm_j follows from the unmodified codimension-one order relation in (C16): the amplitude has order j=mj+(n−2)/4j=m_j+(n-2)/4. The factor (2π)−1(2\pi)^{-1} and iji^j have not been absorbed. Formula (GA1) gives the same order by dyadic estimation, including every tangent derivative. For u∈A′u\in\mathcal A' define gj=(∇nint,ju)∣t=0g_j=(\nabla_n^{\mathrm{int},j}u)|_{t=0} by (GD8)–(GD9). These are tangential distributions at this stage. For a smooth boundary uu, direct differentiation and the distributional delta sign give

u(h⊗δ(j))=(−1)j⟨∂tju∣t=0,h⟩=(−i)j⟨gj,h⟩.(SP2) u(h\otimes\delta^{(j)})=(-1)^j \langle\partial_t^ju|_{t=0},h\rangle =(-i)^j\langle g_j,h\rangle. \tag{SP2}

Both sides are weakly continuous on A′\mathcal A': the left is one of its defining conormal-test pairings by (SP1), and the right is a composition of the weakly continuous corrected derivative and trace maps (GD8)–(GD9). Smooth functions are weakly dense by (GD5), so (SP2) holds for every u∈A′u\in\mathcal A', with no assertion that the ambient uncorrected derivative has the same trace.

2.2. Quantitative scaled test expansion in the actual topology

Fix J=[1,2]J=[1,2] and a compact tangential chart set. For ψ∈Cc∞(X0×(1,2))\psi\in C_c^\infty(X_0\times(1,2)), put

vε(x′,t)=ε−1ψ(x′,t/ε),Mj(x′)=∫12sjψ(x′,s) ds.(SP3) v_\varepsilon(x',t)=\varepsilon^{-1} \psi(x',t/\varepsilon),\qquad M_j(x')=\int_1^2s^j\psi(x',s)\,ds. \tag{SP3}

The test is smooth and supported in the interior for each ε>0\varepsilon>0. Taylor's formula at the actual boundary gives the distributional expansion

vε=∑j=0N−1(−1)jεjj!Mj(x′)⊗δ(j)(t)+RN,ε.(SP4) v_\varepsilon =\sum_{j=0}^{N-1} \frac{(-1)^j\varepsilon^j}{j!} M_j(x')\otimes\delta^{(j)}(t) +R_{N,\varepsilon}. \tag{SP4}

We need its quantitative conormal topology, not only distributional convergence. For every integer N≥1N\ge1, real M>mδ+NM>m_\delta+N with mδ=(2−n)/4m_\delta=(2-n)/4, compact output set KK, and finite tangent seminorm index LL, there are finite L′L', CC such that

pM,K,L(RN,ε)≤CεN∑∣γ∣≤L′∥Dx′,sγψ∥L∞,0<ε≤1.(SP5) p_{M,K,L}(R_{N,\varepsilon}) \le C\varepsilon^N \sum_{|\gamma|\le L'} \|D_{x',s}^\gamma\psi\|_{L^\infty}, \qquad0<\varepsilon\le1. \tag{SP5}

Here the original Besov index is κ=−M−n/4\kappa=-M-n/4; the strict condition is exactly κ+N+1/2<0\kappa+N+1/2<0.

For completeness, Fourier transform (SP4) in tt. Its normal factor is

∫12e−iεsτψ(x′,s) ds−∑j=0N−1(−iετ)jj!Mj(x′).(SP6) \int_1^2e^{-i\varepsilon s\tau}\psi(x',s)\,ds -\sum_{j=0}^{N-1} \frac{(-i\varepsilon\tau)^j}{j!}M_j(x'). \tag{SP6}

For ε∣τ∣≤1\varepsilon|\tau|\le1, Taylor's integral remainder bounds (SP6), with every tangential derivative, by CεN∣τ∣NC\varepsilon^N|\tau|^N; for ε∣τ∣≥1\varepsilon|\tau|\ge1, the Schwartz integral and all retained polynomial terms bound it by C(ε∣τ∣)N−1C(\varepsilon|\tau|)^{N-1}. The tangential Fourier transform decays faster than any power of ∣ξ′∣|\xi'|, with constants controlled by finitely many displayed ψ\psi derivatives. On a full dyadic block ∣(ξ′,τ)∣≍2l|(\xi',\tau)|\asymp2^l, the L2L^2 size in the normal frequency contributes 2l/22^{l/2}. For 2l≤ε−12^l\le\varepsilon^{-1}, multiplying by 2lκ2^{l\kappa} gives CεN2l(κ+N+1/2)C\varepsilon^N2^{l(\kappa+N+1/2)}, bounded by CεNC\varepsilon^N. For 2l≥ε−12^l\ge\varepsilon^{-1}, it gives CεN−12l(κ+N−1/2)C\varepsilon^{N-1}2^{l(\kappa+N-1/2)}, whose maximum is Cε−κ−1/2≤CεNC\varepsilon^{-\kappa-1/2}\le C\varepsilon^N under the same strict condition. The low block is bounded directly by the Taylor remainder. Tangentially dominant blocks gain arbitrary decay from the x′x'-Fourier transform and obey the same inequality. Applying D′D' differentiates ψ\psi. Applying tDttD_t to (SP4) rescales the same test and acts on each δ(j)\delta^{(j)} by its exact eigenvalue tDtδ(j)=i(j+1)δ(j)tD_t\delta^{(j)}=i(j+1)\delta^{(j)}; the Taylor remainder still has the two bounds above. The triangular identity (GA2) therefore supplies every weighted derivative seminorm, proving (SP5) with all lower terms retained.

The dual estimate (GD3), used at this freely chosen high order MM, turns (SP5) into

∣u(RN,ε)∣≤Cu,N,KεN∑∣γ∣≤L′∥Dx′,sγψ∥∞.(SP7) |u(R_{N,\varepsilon})| \le C_{u,N,K}\varepsilon^N \sum_{|\gamma|\le L'} \|D_{x',s}^\gamma\psi\|_\infty. \tag{SP7}

This is an estimate of distributions in (x′,s)(x',s) of a finite negative Sobolev order, since a sufficiently high Sobolev norm controls the finite smooth-test seminorm.

2.3. The scaled Taylor series in distributions

An element u∈A′u\in\mathcal A' is smooth in the open collar only when it also lies in A\mathcal A. Assume that intersection from now on, and set Fε(x′,s)=u(x′,εs)F_\varepsilon(x',s)=u(x',\varepsilon s) for s∈Js\in J. For every ψ\psi in (SP3), the change of variables gives ⟨Fε,ψ⟩=u(vε)\langle F_\varepsilon,\psi\rangle=u(v_\varepsilon). Insert (SP4), then use the exact sign (SP2). The two (−1)j(-1)^j factors cancel, leaving

Fε(x′,s)=∑j=0N−1ijεjsjj!gj(x′)+OH−LN(K0×J)(εN)for every N≥1.(SP8) F_\varepsilon(x',s) =\sum_{j=0}^{N-1} \frac{i^j\varepsilon^j s^j}{j!}g_j(x') +O_{H^{-L_N}(K_0\times J)}(\varepsilon^N) \quad\text{for every }N\ge1. \tag{SP8}

The remainder means (SP7) uniformly on compact tangential sets; LNL_N is finite but may depend on NN. Formula (SP8) is an all-order distribution-valued boundary Taylor series, with the original iji^j from D=−i∂D=-i\partial. No smoothness of gjg_j has been assumed.

2.4. One fixed conormal growth exponent for every derivative

Let u∈Amu\in\mathcal A^m, with the original normal amplitude order μ=m+(n−2)/4\mu=m+(n-2)/4. The frequency split in (HS1)–(HS2) gives, after enlarging a fixed exponent slightly to absorb a possible logarithm, a number A>max⁡(μ+1,0)A>\max(\mu+1,0) such that for every normal derivative order aa, tangential multiindex β\beta, and compact subchart there is a constant CaβC_{a\beta} with

sup⁡(x′,s)∈K0×J∣Dx′β∂saFε(x′,s)∣≤Caβε−A.(SP9) \sup_{(x',s)\in K_0\times J} |D_{x'}^\beta\partial_s^a F_\varepsilon(x',s)| \le C_{a\beta}\varepsilon^{-A}. \tag{SP9}

Indeed ∂saFε=εa∂tau(x′,εs)\partial_s^a F_\varepsilon =\varepsilon^a\partial_t^au(x',\varepsilon s); the conormal amplitude order rises by exactly aa, so its high-frequency bound contains εaε−μ−a−1=ε−μ−1\varepsilon^a\varepsilon^{-\mu-a-1} =\varepsilon^{-\mu-1}. The low-frequency and logarithmic cases satisfy the same enlarged AA. Crucially the exponent AA is independent of aa and β\beta, although the constants depend on them. Thus for every integer KK, ∥Fε∥HK(K0×J)≤CKε−A\|F_\varepsilon\|_{H^K(K_0\times J)} \le C_K\varepsilon^{-A}.

We record the precise smoothing inference used below. Suppose smooth hεh_\varepsilon on a fixed compact coordinate box obey ∥hε∥HK≤CKε−A\|h_\varepsilon\|_{H^K}\le C_K\varepsilon^{-A} for all KK, and for every NN obey ∥hε−h∥H−LN≤CNεN\|h_\varepsilon-h\|_{H^{-L_N}} \le C_N\varepsilon^N, where hh is a distribution. Then hh is smooth. Fix any one positive integer NN, and retain its finite negative order LNL_N. On a compact subbox insert a fixed smooth cutoff before using full Fourier blocks; multiplication preserves the displayed Sobolev bounds. For any a>0a>0, choose ε=2−al\varepsilon=2^{-al}. The low- and high-Sobolev block bounds give

∥Δlh∥2≤CN2l(LN−aN)+CK2l(aA−K).(SP10) \|\Delta_l h\|_2 \le C_N 2^{l(L_N-aN)} +C_K2^{l(aA-K)}. \tag{SP10}

For any requested r≥0r\ge0, choose a>(LN+r+2)/Na>(L_N+r+2)/N, then an integer K>aA+r+2K>aA+r+2. Both exponents in (SP10) are strictly less than −r−2-r-2, so it gives ∥Δlh∥2≤Cr2−l(r+2)\|\Delta_l h\|_2\le C_r2^{-l(r+2)} for every rr. Summing their squared HrH^r weights proves membership in every local Sobolev space, and the Fourier Sobolev estimate proves smoothness. This argument uses the actual LNL_N; it does not assert one negative order valid for all Taylor degrees. The single positive degree N=1N=1 would suffice for this smoothing inference.

2.5. Smoothness of every boundary coefficient

The expansion (SP8) includes lower powers of ε\varepsilon, so apply an exact finite scale cancellation. For any integer N≥1N\ge1, take distinct positive scales λ1,…,λN\lambda_1,\ldots,\lambda_N, for example 1,…,N1,\ldots,N, and the Lagrange interpolation coefficients at zero

cl=∏r≠l−λrλl−λr,∑l=1Nclλlj={1,j=0,0,1≤j<N.(SP11) c_l=\prod_{r\ne l} \frac{-\lambda_r}{\lambda_l-\lambda_r}, \qquad \sum_{l=1}^{N}c_l\lambda_l^j =\begin{cases}1,&j=0,\\0,&1\le j<N. \end{cases} \tag{SP11}

To verify the identity for every displayed degree, the polynomials Ll(z)=∏r≠l(z−λr)/(λl−λr)L_l(z)=\prod_{r\ne l}(z-\lambda_r)/(\lambda_l-\lambda_r) equal one at their own node and zero at the others. For a polynomial of degree less than NN, subtract its interpolation sum: the difference has those NN distinct roots and degree less than NN, so it is zero, by successive division by its linear factors. Evaluation at zero gives (SP11), with every sign shown. This also covers N=1N=1.

For g0g_0, set GN,ε=∑lclFλlεG_{N,\varepsilon}=\sum_lc_lF_{\lambda_l\varepsilon}. Equations (SP8) and (SP11) make GN,ε=g0+OH−LN(εN)G_{N,\varepsilon}=g_0+O_{H^{-L_N}}(\varepsilon^N). Equation (SP9) gives every high Sobolev norm bounded by CK,Nε−AC_{K,N}\varepsilon^{-A}. The dyadic argument (SP10), with the adjustable scale specified there, proves g0∈C∞g_0\in C^\infty.

Proceed by induction. If g0,…,gj−1g_0,\ldots,g_{j-1} are smooth, define on s∈Js\in J

Hj,ε(x′,s)=j!i−js−jε−j[Fε(x′,s)−∑k=0j−1ikεkskk!gk(x′)].(SP12) H_{j,\varepsilon}(x',s) =j!i^{-j}s^{-j}\varepsilon^{-j} \left[F_\varepsilon(x',s) -\sum_{k=0}^{j-1} \frac{i^k\varepsilon^ks^k}{k!}g_k(x')\right]. \tag{SP12}

The factor s−js^{-j} is smooth on JJ, and every term is a smooth function there. Formula (SP8), taken to order N+jN+j, gives Hj,ε=gj+∑l=1N−1εlhl(x′,s)+OH−L(εN)H_{j,\varepsilon}=g_j+\sum_{l=1}^{N-1} \varepsilon^l h_l(x',s)+O_{H^{-L}}(\varepsilon^N), where the hlh_l may still be distributions. The same Lagrange weights (SP11) cancel all intermediate powers. The high derivative bound (SP9), the smooth already-known coefficients and (SP12) give ∥Hj,ε∥HK≤CjKε−(A+j)\|H_{j,\varepsilon}\|_{H^K} \le C_{jK}\varepsilon^{-(A+j)}. The dyadic argument proves gj∈C∞g_j\in C^\infty. Induction proves

gj∈C∞(∂X)for every j≥0.(SP13) g_j\in C^\infty(\partial X) \quad\text{for every }j\ge0. \tag{SP13}

The argument uses the same original uu at every step; it neither assumes nor constructs an unrelated boundary extension.

2.6. Upgrading the whole series and removing the weight

Now every coefficient in (SP8) is smooth. Fix a desired integer NN and a smooth seminorm order rr. Expand (SP8) through some N′>NN'>N, so the remainder is OH−LN′(εN′)O_{H^{-L_{N'}}}(\varepsilon^{N'}). Its high HKH^K norm is at most CKε−AC_K\varepsilon^{-A} by (SP9) and the smooth polynomial terms. Interpolation between H−LN′H^{-L_{N'}} and HKH^K gives the required bound as follows. Choose an integer d>r+n/2d>r+n/2, then N′>N+AN'>N+A, retaining its actual finite L=LN′L=L_{N'}. With K>dK>d, Fourier Hölder gives the HdH^d exponent θN′−(1−θ)A\theta N'-(1-\theta)A, where θ=(K−d)/(K+L)\theta=(K-d)/(K+L). Choose KK so large that θ>(N+A)/(N′+A)\theta>(N+A)/(N'+A); the exponent is then greater than NN. The Fourier Sobolev bound Hd↪CrH^d\hookrightarrow C^r proves OCr(εN)O_{C^r}(\varepsilon^N). The discarded coefficients of degrees N,…,N′−1N,\ldots,N'-1 contribute the same OCr(εN)O_{C^r}(\varepsilon^N). Hence the full exact Taylor statement is

u(x′,εs)=∑j=0N−1ijεjsjj!gj(x′)+OCr(K0×J)(εN)for every N,r.(SP14) u(x',\varepsilon s) =\sum_{j=0}^{N-1} \frac{i^j\varepsilon^js^j}{j!}g_j(x') +O_{C^r(K_0\times J)}(\varepsilon^N) \quad\text{for every }N,r. \tag{SP14}

Taking any fixed interior value, for example s=3/2s=3/2 and ε=2t/3\varepsilon=2t/3, then differentiating in ss and x′x' before that restriction, shows that every ordinary mixed derivative of uu has a continuous boundary limit with the precise Taylor coefficient in (SP14). This proves u∈C∞(X)u\in C^\infty(X). Conversely a smooth boundary function belongs to A\mathcal A by (GD1)–(GD2), and to A′\mathcal A' by the smooth pairing and (GD3). Therefore, as spaces of their actual supported representatives,

A′(X)∩A(X)=C∞(X).(SP15) \mathcal A'(X)\cap\mathcal A(X)=C^\infty(X). \tag{SP15}

3. Identification of the actual supported representative

The preceding bounds give an interior smooth function whose derivatives of every order extend continuously to the boundary. To check the derivative interpretation there, apply the fundamental theorem of calculus on t>0t>0 and pass to zero using the uniform limits. Repeating it for each tangential and normal derivative proves smoothness up to the boundary with precisely the jets in (SP14). In passing from a scaled ss-derivative of order aa to an ordinary tt-derivative, the remainder acquires ϵ−a\epsilon^{-a}; choose its Taylor degree N>aN>a before taking the limit. Thus that step retains every ordinary derivative, not just the scaled ones.

Let vv be this smooth boundary function, represented by zero extension. It belongs to A′\mathcal A' by the smooth-pairing proof and to A\mathcal A by (GD1)–(GD2). The original distribution u−vu-v belongs to A′\mathcal A' and has zero interior restriction. Injectivity (GD4) makes it zero. This proves (SP15) for the actual supported distributions as stated, including the absence of an unnoticed boundary-delta summand.