2025-11-23T13:07:17.414241

Robust nonlocal trace and extension theorems

Grube, Kassmann
We prove trace and extension results for Sobolev-type function spaces that are well suited for nonlocal Dirichlet and Neumann problems including those for the fractional $p$-Laplacian. Our results are robust with respect to the order of differentiability. In this sense they are in align with the classical trace and extension theorems.
academic

Robust nonlocal trace and extension theorems

Basic Information

  • Paper ID: 2305.05735
  • Title: Robust nonlocal trace and extension theorems
  • Authors: Florian Grube, Moritz Kassmann (Universität Bielefeld)
  • Classification: math.AP (Analysis of PDEs)
  • Publication Date: May 2023 (v1), June 2023 (v2)
  • Paper Link: https://arxiv.org/abs/2305.05735

Abstract

This paper establishes trace theorems and extension theorems for Sobolev-type function spaces applicable to nonlocal Dirichlet and Neumann problems, including fractional p-Laplacians. These results are robust with respect to the order of differentiability, maintaining consistency with classical trace and extension theorems in this sense.

Research Background and Motivation

Core Problem

The paper studies nonlinear nonlocal equations on bounded domains, such as: (Δ)psu=f in Ω,u=g in RdΩ(-\Delta)_p^s u = f \text{ in } \Omega, \quad u = g \text{ in } \mathbb{R}^d \setminus \Omega

where the fractional p-Laplacian is defined as: (Δ)psu(x)=(1s)p.v.Rdu(x)u(y)p2(u(x)u(y))xyd+spdy(-\Delta)_p^s u(x) = (1-s) \text{p.v.} \int_{\mathbb{R}^d} \frac{|u(x)-u(y)|^{p-2}(u(x)-u(y))}{|x-y|^{d+sp}} dy

Research Significance

  1. Essential characteristics of nonlocal operators: Due to the nonlocal nature of the operator, problem (1.1) requires specifying values u(x)u(x) on RdΩ\mathbb{R}^d \setminus \Omega for well-posedness
  2. Limitations of existing methods: Traditional assumptions such as gWs,p(Rd)g \in W^{s,p}(\mathbb{R}^d) impose unnatural restrictions, since problem (1.1) does not involve any regularity of gg in RdΩ\mathbb{R}^d \setminus \Omega, only requiring certain weighted integrability
  3. Missing continuity: Existing theory lacks continuity as s1s \to 1^-, failing to recover classical local theory

Limitations of Existing Approaches

  • Homogeneous problems (g=0g=0) are relatively simple and have been studied
  • Treatment of nonzero data gg typically assumes gg is regular on the entire space Rd\mathbb{R}^d, such as gWs,p(Rd)g \in W^{s,p}(\mathbb{R}^d)
  • Compromise approaches like gWs,p(Ωϵ)Lp(Rd;(1+x)dspdx)g \in W^{s,p}(\Omega_\epsilon) \cap L^p(\mathbb{R}^d; (1+|x|)^{-d-sp}dx) remain unnatural
  • As s1s \to 1^-, the classical trace space W11/p,p(Ω)W^{1-1/p,p}(\partial\Omega) cannot be recovered

Research Motivation

Construct trace spaces Ts,p(Ωc)T^{s,p}(\Omega^c) suitable for nonlocal variational methods such that:

  1. Naturally characterize the regularity of external data gg
  2. Be robust with respect to parameter ss
  3. Recover classical trace and extension theorems as s1s \to 1^-

Core Contributions

  1. Introduction of new trace spaces: Define function space Ts,p(Ωc)T^{s,p}(\Omega^c) (defined on Ωc\Omega^c) equipped with measure μs(dx):=1Ωc(x)(1s)dxs(1+dx)ds(p1)dx\mu_s(dx) := 1_{\Omega^c}(x)(1-s)d_x^{-s}(1+d_x)^{-d-s(p-1)}dx where dx:=dist(x,Ω)d_x := \text{dist}(x, \partial\Omega)
  2. Proof of robust trace theorem (Theorem 1.2, p>1p>1): The trace operator Trs:Vs,p(ΩRd)Ts,p(Ωc)\text{Tr}_s: V^{s,p}(\Omega|\mathbb{R}^d) \to T^{s,p}(\Omega^c) is continuous linear, and there exists a continuous linear right inverse (extension operator) Exts:Ts,p(Ωc)Vs,p(ΩRd)\text{Ext}_s: T^{s,p}(\Omega^c) \to V^{s,p}(\Omega|\mathbb{R}^d), with continuity constants depending only on Ω\Omega, lower bounds of ss, and upper and lower bounds of pp
  3. Extension to p=1p=1 case (Theorem 1.3): Prove that the trace operator Trs:Vs,1(ΩRd)L1(Ωc;μs)\text{Tr}_s: V^{s,1}(\Omega|\mathbb{R}^d) \to L^1(\Omega^c; \mu_s) is continuous, and there exists a continuous linear right inverse Exts:Ts,1(Ωc)Vs,1(ΩRd)\text{Ext}_s: T^{s,1}(\Omega^c) \to V^{s,1}(\Omega|\mathbb{R}^d)
  4. Limiting behavior theorem (Theorem 1.4): Prove that as s1s \to 1^-:
    • TrsuLp(Ωc;μs)γuLp(Ω)\|\text{Tr}_s u\|_{L^p(\Omega^c;\mu_s)} \to \|\gamma u\|_{L^p(\partial\Omega)} (uW1,p(Rd)u \in W^{1,p}(\mathbb{R}^d))
    • [Trsu]Ts,p(Ωc)[γu]W11/p,p(Ω)[\text{Tr}_s u]_{T^{s,p}(\Omega^c)} \to [\gamma u]_{W^{1-1/p,p}(\partial\Omega)} (uW1,p(Rd)u \in W^{1,p}(\mathbb{R}^d))
    • Recover classical trace spaces and Besov spaces
  5. Application to Dirichlet problems (Corollary 1.7): Using the new trace and extension theory, establish well-posedness and energy estimates for problem (1.1)

Detailed Methodology

Task Definition

Study the mapping from energy space Vs,p(ΩRd)V^{s,p}(\Omega|\mathbb{R}^d) to trace spaces on exterior domain Ωc\Omega^c: Vs,p(ΩRd):={u:RdR measurable[u]Vs,p(ΩRd)<}V^{s,p}(\Omega|\mathbb{R}^d) := \{u: \mathbb{R}^d \to \mathbb{R} \text{ measurable} \mid [u]_{V^{s,p}(\Omega|\mathbb{R}^d)} < \infty\} where the seminorm is defined as: [u]Vs,p(AB)p:=(1s)A×Bu(x)u(y)pxyd+spdxdy[u]_{V^{s,p}(A|B)}^p := (1-s) \iint_{A \times B} \frac{|u(x)-u(y)|^p}{|x-y|^{d+sp}} dx dy

Construction of Trace Spaces

Key Innovation: Design of measure μs\mu_sμs(dx)=1Ωc(x)(1s)dxs(1+dx)ds(p1)dx\mu_s(dx) = 1_{\Omega^c}(x)(1-s)d_x^{-s}(1+d_x)^{-d-s(p-1)}dx

This measure design has the following characteristics:

  1. Behavior near boundary: The (1s)dxs(1-s)d_x^{-s} term ensures that as s1s \to 1^-, μs\mu_s weakly converges to the Hausdorff measure on boundary Ω\partial\Omega
  2. Decay at infinity: The (1+dx)ds(p1)(1+d_x)^{-d-s(p-1)} term ensures appropriate decay at infinity
  3. Parameter robustness: The choice of coefficient (1s)(1-s) makes the measure scaling natural as s1s \to 1^-

Definition of trace space Ts,p(Ωc)T^{s,p}(\Omega^c) (equation 1.6): Ts,p(Ωc):={g:ΩcR measurablegTs,p(Ωc)<}T^{s,p}(\Omega^c) := \{g: \Omega^c \to \mathbb{R} \text{ measurable} \mid \|g\|_{T^{s,p}(\Omega^c)} < \infty\} with norm: gTs,p(Ωc)p:=gLp(Ωc;μs)p+[g]Ts,p(Ωc)p\|g\|_{T^{s,p}(\Omega^c)}^p := \|g\|_{L^p(\Omega^c;\mu_s)}^p + [g]_{T^{s,p}(\Omega^c)}^p

The seminorm is defined as (for p>1p>1): [f,g]Ts,p(Ωc)p:=Ωc×Ωcf(x)f(y)p2(f(x)f(y))(g(x)g(y))((xy+dx+dy)1)d+s(p2)μs(dx)μs(dy)[f,g]_{T^{s,p}(\Omega^c)}^p := \iint_{\Omega^c \times \Omega^c} \frac{|f(x)-f(y)|^{p-2}(f(x)-f(y))(g(x)-g(y))}{((|x-y|+d_x+d_y) \wedge 1)^{d+s(p-2)}} \mu_s(dx)\mu_s(dy)

Key design principles:

  • The denominator term (xy+dx+dy)1(|x-y|+d_x+d_y) \wedge 1 becomes (xy)1(|x-y|) \wedge 1 near the boundary, recovering the classical Slobodeckij seminorm structure
  • The weight in measure μs\mu_s compensates for the dimensional change (from dd dimensions to (d1)(d-1) dimensional boundary)

Proof Strategy for Trace Theorem

Section 3: Proof of trace embedding

The proof is divided into two propositions:

  1. Proposition 3.9: Embedding of LpL^p part TrsuLp(Ωc;μs)CuVs,p(ΩRd)\|\text{Tr}_s u\|_{L^p(\Omega^c;\mu_s)} \leq C \|u\|_{V^{s,p}(\Omega|\mathbb{R}^d)}
  2. Proposition 3.10: Embedding of seminorm part [Trsu]Ts,p(Ωc)CuVs,p(ΩRd)[\text{Tr}_s u]_{T^{s,p}(\Omega^c)} \leq C \|u\|_{V^{s,p}(\Omega|\mathbb{R}^d)}

Key technical tools:

(a) Approximate trace inequality (Theorem 3.5): For uWs,p(Ω)u \in W^{s,p}(\Omega), prove Ωu(x)pτs(dx)+ΩΩu(x)u(y)p((xy+dx+dy)1)d+s(p2)τs(dy)τs(dx)CuWs,p(Ω)p\int_{\Omega} |u(x)|^p \tau_s(dx) + \int_{\Omega} \int_{\Omega} \frac{|u(x)-u(y)|^p}{((|x-y|+d_x+d_y) \wedge 1)^{d+s(p-2)}} \tau_s(dy)\tau_s(dx) \leq C \|u\|_{W^{s,p}(\Omega)}^p

where τs(dx)=(1s)dxs1Ω(x)dx\tau_s(dx) = (1-s)d_x^{-s}1_{\Omega}(x)dx

Proof method:

  • Use interpolation theory of Bessel potential spaces
  • Define α0:=s(1+p)/(2p)\alpha_0 := s(1+p)/(2p), α1:=1+s/(2p)\alpha_1 := 1 + s/(2p)
  • Prove continuity for Hαi,p(Rd)H^{\alpha_i,p}(\mathbb{R}^d) (Lemma 3.4)
  • Obtain result via real interpolation [Hα0,p,Hα1,p]θp=Ws,p(Rd)[H^{\alpha_0,p}, H^{\alpha_1,p}]_\theta^p = W^{s,p}(\mathbb{R}^d)

(b) Hardy-type inequality (Theorem 3.6, p=1p=1 case): (1s)Ωu(x)dxsdxC(uL1(Ω)+s(1s)Ω×Ωu(x)u(y)xyd+sd(x,y))(1-s) \int_{\Omega} \frac{|u(x)|}{d_x^s} dx \leq C\left(\|u\|_{L^1(\Omega)} + s(1-s) \int_{\Omega \times \Omega} \frac{|u(x)-u(y)|}{|x-y|^{d+s}} d(x,y)\right)

Proof strategy:

  • Use local Lipschitz property of boundary
  • Apply optimal Hardy inequality on half-spaces (Theorem B.1)
  • Handle global case via partition of unity

(c) Technical lemmas on distance functions:

  • Lemma 3.7: Convert dxsd_x^{-s} to integrals over Ω\Omega
  • Lemma 3.8: Establish lower bounds using uniform interior cone condition

Proof of Proposition 3.9 (LpL^p part):

  1. Decompose Ωc=Ω1extΩ1ext\Omega^c = \Omega_1^{\text{ext}} \cup \Omega^{\text{ext}}_1
  2. On Ω1ext\Omega_1^{\text{ext}}: Use Lemma 3.8 and Theorem 3.5/3.6
  3. On Ω1ext\Omega^{\text{ext}}_1: Direct estimation using distance function lower bounds

Proof of Proposition 3.10 (seminorm part, p>1p>1):

  1. Decompose integral domain: Ωρext×Ωρext\Omega^{\text{ext}}_\rho \times \Omega^{\text{ext}}_\rho, Ωc×Ωextρ\Omega^c \times \Omega^\rho_{\text{ext}}, Ωextρ×Ωc\Omega^\rho_{\text{ext}} \times \Omega^c
  2. Far-field part (Ωextρ\Omega^\rho_{\text{ext}}): Use Proposition 3.9 and distance lower bounds
  3. Near-field part (Ωρext×Ωρext\Omega^{\text{ext}}_\rho \times \Omega^{\text{ext}}_\rho):
    • Apply Lemma 3.8 twice, introducing intermediate integration variables z,wΩz,w \in \Omega
    • Apply Theorem 3.5 to estimate integral over Ω×Ω\Omega \times \Omega
    • Key: Refined estimation of kernel function a(z,w)a(z,w) (equation 3.19)

Proof Strategy for Extension Theorem

Section 4: Construction of extension operator

Core idea: Weighted averaging extension based on Whitney decomposition

(a) Whitney decomposition:

  • Perform dyadic Whitney decomposition W(Ω)\mathcal{W}(\Omega) of Ω\Omega
  • Each cube QQ satisfies: diam(Q)d(Q,Ω)4diam(Q)\text{diam}(Q) \leq d(Q,\partial\Omega) \leq 4\text{diam}(Q)
  • Define Wκ(Ω)\mathcal{W}_{\leq\kappa}(\Omega) as cubes with side length 2κ\leq 2^\kappa, where κ:=log2(ρ/d)\kappa := \lfloor\log_2(\rho/\sqrt{d})\rfloor, ρ:=inr(Ω)/21/2\rho := \text{inr}(\Omega)/2 \wedge 1/2

(b) Partition of unity: Construct smooth partition of unity {ϕQ}QW(Ω)\{\phi_Q\}_{Q \in \mathcal{W}(\Omega)} satisfying QϕQ=1Ω\sum_Q \phi_Q = 1_\Omega

(c) Measure μ~s\tilde{\mu}_s: Define simplified measure (without decay at infinity): μ~s(dz)=1Ωc(z)1sdzsdz\tilde{\mu}_s(dz) = 1_{\Omega^c}(z)\frac{1-s}{d_z^s} dz and normalization constant: aQ,s:=(μ~s(B6lQ(qQ)))1a_{Q,s} := \left(\tilde{\mu}_s(B_{6l_Q}(q_Q))\right)^{-1}

Lemma 4.1 proves: C1lQsdaQ,sClQsdC^{-1}l_Q^{s-d} \leq a_{Q,s} \leq C l_Q^{s-d}

(d) Definition of extension operator (equation 4.11):

\sum_{Q \in \mathcal{W}_{\leq\kappa}(\Omega)} \phi_Q(x) a_{Q,s} \int_{\Omega^c \cap B_{6l_Q}(q_Q)} g(z) \tilde{\mu}_s(dz) & x \in \Omega \\ g(x) & x \in \Omega^c \end{cases}$$ **Design principles**: - Within cube $Q$, the extension value is the weighted average of $g$ over neighborhood $B_{6l_Q}(q_Q) \cap \Omega^c$ with respect to $\tilde{\mu}_s$ - Weight $a_{Q,s}$ ensures normalization - Partition of unity ensures smooth transition **Key technical lemma** (Lemma 4.2): For $Q_1, Q_2 \in \mathcal{W}_{\leq\kappa-2}(\Omega)$, $x \in Q_1$, $y \in Q_2$: (a) $|\text{Ext}_s(g)(x) - \text{Ext}_s(g)(y)| \leq C J_p(q_{Q_1}, q_{Q_2})$ (b) $|\nabla\text{Ext}_s(g)(x)| \leq C l_{Q_1}^{-1} J_p(q_{Q_1}, q_{Q_2})$ where $$J_p(q_{Q_1}, q_{Q_2}) := \left(a_{Q_1,s} a_{Q_2,s} \int_{B_{30l_{Q_1}}(q_{Q_1})} \int_{B_{30l_{Q_2}}(q_{Q_2})} |g(z_1)-g(z_2)|^p \tilde{\mu}_s(dz_2)\tilde{\mu}_s(dz_1)\right)^{1/p}$$ **Proposition 4.5** ($L^p$ part): $$\|\text{Ext}_s(g)\|_{L^p(\Omega)} \leq \frac{C}{s^{1/p}} \|g\|_{L^p(\Omega_3^{\text{ext}};\mu_s)}$$ Key points of proof: - Use finite overlap property of Whitney cubes - Estimate number of cubes at same scale (equation 4.14) - Change order of summation **Proposition 4.6** (seminorm part): $$[\text{Ext}_s(g)]_{V^{s,p}(\Omega|\mathbb{R}^d)} \leq \frac{C}{(d+s(p-2))^{1/p}s^{2/p}} \|g\|_{T^{s,p}(\Omega^c)}$$ Proof strategy: 1. Decompose $[u]_{V^{s,p}(\Omega|\mathbb{R}^d)}^p$ into four parts (based on $|h| \geq c_1$, $|h| < c_1 2^{-j}$, etc.) 2. Apply different techniques to each part: - Far-field: Direct estimation - Near-field: Use Lemma 4.2 and fundamental theorem - Mid-field: Whitney decomposition technique 3. Key: Lemma 4.3 (measure scaling lemma) and Lemma 4.4 (series summation lemma) ### Technical Innovations 1. **Sophistication of measure design**: - Construction of $\mu_s$ simultaneously handles boundary behavior and decay at infinity - Factor $(1-s)$ ensures correct scaling as $s \to 1^-$ - Weak convergence to surface measure (Lemma 5.1) 2. **Application of interpolation theory**: - Use real interpolation of Bessel potential spaces to avoid direct estimation - Constants independent of $s$ choice (Theorem 3.1) 3. **Improvement of Whitney extension**: - Use measure $\tilde{\mu}_s$ rather than Lebesgue measure - Ensure recovery of classical extension as $s \to 1^-$ 4. **Unified treatment of $p \geq 1$**: - Differences between $p>1$ and $p=1$ handled via different Hardy inequalities - For $p=1$, recover BV space and Besov space $B_1^{0,1}(\partial\Omega)$ ## Experimental Setup **Note**: This is a pure theoretical mathematics paper with no numerical experiments. Main results are rigorous proofs of theorems. ### Theoretical Verification Framework **Verification content**: 1. Continuity of trace and extension operators 2. Dependence of constants on parameters 3. Limiting behavior as $s \to 1^-$ **Methods**: - Constructive proofs (explicit operator definitions) - Precise constant tracking - Weak convergence analysis ## Experimental Results ### Main Theoretical Results **Theorem 1.2** ($1<p<\infty$): - Trace operator $\text{Tr}_s: V^{s,p}(\Omega|\mathbb{R}^d) \to T^{s,p}(\Omega^c)$ is continuous linear - Extension operator $\text{Ext}_s: T^{s,p}(\Omega^c) \to V^{s,p}(\Omega|\mathbb{R}^d)$ is continuous linear - Continuity constants depend only on $\Omega$, lower bounds of $s$, and upper/lower bounds of $p$ **Theorem 1.3** ($p=1$): - Trace operator $\text{Tr}_s: V^{s,1}(\Omega|\mathbb{R}^d) \to L^1(\Omega^c;\mu_s)$ is continuous linear - Extension operator $\text{Ext}_s: T^{s,1}(\Omega^c) \to V^{s,1}(\Omega|\mathbb{R}^d)$ is continuous linear - For $d=1$, constants additionally depend on lower bounds of $1-s$ **Theorem 1.4** (Limiting behavior): For $u \in W^{1,p}(\mathbb{R}^d)$ ($1<p<\infty$): $$\|\text{Tr}_s u\|_{L^p(\Omega^c;\mu_s)} \to \|\gamma u\|_{L^p(\partial\Omega)}$$ $$[\text{Tr}_s u]_{T^{s,p}(\Omega^c)} \to [\gamma u]_{W^{1-1/p,p}(\partial\Omega)}$$ as $s \to 1^-$. For $u \in BV(\mathbb{R}^d)$ ($p=1$): $$\|\text{Tr}_s u\|_{L^1(\Omega^c;\mu_s)} \to \|\gamma u\|_{L^1(\partial\Omega)}$$ For $u \in C_c^{0,1}(\mathbb{R}^d)$: $$[\text{Tr}_s u]_{T^{s,1}(\Omega^c)} \to [\gamma u]_{B_1^{0,1}(\partial\Omega)}$$ **Corollary 1.7** (Well-posedness of Dirichlet problem): For $g \in T^{s,p}(\Omega^c)$, $f \in V^{s,p}(\Omega|\mathbb{R}^d)'$, there exists a unique weak solution $u \in V^{s,p}(\Omega|\mathbb{R}^d)$ to problem (1.1), with $$\|u\|_{V^{s,p}(\Omega|\mathbb{R}^d)} \leq c(\|g\|_{T^{s,p}(\Omega^c)} + \|f\|_{V^{s,p}(\Omega|\mathbb{R}^d)'})$$ ### Key Estimates **Proposition 3.9**: There exists $C = C(\Omega, p_\star, s_\star) > 0$ such that $$\|\text{Tr}_s u\|_{L^p(\Omega^c;\mu_s)} \leq C \|u\|_{V^{s,p}(\Omega|\mathbb{R}^d)}$$ for all $s \in (s_\star, 1)$, $1 \leq p \leq p_\star$ **Proposition 3.10**: There exists $C = C(\Omega, p_\star, p^\star, s_\star) > 0$ such that $$[\text{Tr}_s u]_{T^{s,p}(\Omega^c)} \leq C \|u\|_{V^{s,p}(\Omega|\mathbb{R}^d)}$$ for all $s \in (s_\star, 1)$, $p_\star \leq p \leq p^\star$ **Proposition 4.5**: $$\|\text{Ext}_s(g)\|_{L^p(\Omega)} \leq \frac{C}{s^{1/p}} \|g\|_{L^p(\Omega_3^{\text{ext}};\mu_s)}$$ **Proposition 4.6**: $$[\text{Ext}_s(g)]_{V^{s,p}(\Omega|\mathbb{R}^d)} \leq \frac{C}{(d+s(p-2))^{1/p}s^{2/p}} \|g\|_{T^{s,p}(\Omega^c)}$$ ### Auxiliary Results **Lemma 5.1** (Weak convergence): The family of measures $\{\mu_s\}_{s \in (0,1)}$ weakly converges to surface measure $\sigma$, i.e., for all $f \in C_c(\mathbb{R}^d)$: $$\int_{\mathbb{R}^d} f d\mu_s \to \int_{\partial\Omega} f d\sigma \quad \text{as } s \to 1^-$$ **Corollary 2.3**: Space $T^{s,p}(\Omega^c)$ continuously embeds into $W^{s-1/p,p}(\partial\Omega)$ (when $s \in (1/p, 1)$, $p \in (1,\infty)$), and the embedding is surjective ## Related Work ### Nonlocal Function Space Theory 1. **Early work**: - Servadei-Valdinoci [SV12, SV13, SV14]: Energy spaces for $p=2$ case - Felsinger-Kassmann-Voigt [FKV15]: Dirichlet problems for nonlocal operators 2. **Trace and extension theory**: - Dyda-Kassmann [DK19]: Trace and extension for $V^{s,p}(\Omega|\mathbb{R}^d)$ (using Lebesgue measure, not robust) - Bogdan-Grzywny-Pietruska-Pałuba-Rutkowski [BGPR20]: $p=2$ case based on Poisson kernel - Jonsson-Wallin [JW78, JW84]: Function spaces on doubling measures 3. **$p=2$ and $C^{1,1}$ domains**: - Grube-Hensiek [GH22]: Earlier work of this paper, treating $p=2$ and $C^{1,1}$ domains 4. **General $p$ and Lévy measures**: - Foghem-Kassmann [FK22]: Systematic study of $L^2$-Lévy integrable kernels - Foghem [Fog23b]: Extension to $p>1$ ### Classical Sobolev Space Theory 1. **Classical trace theorems**: - Gagliardo [Gag57]: $W^{1,p}(\Omega) \to W^{1-1/p,p}(\partial\Omega)$ - Slobodeckij [Slo58]: Fractional Sobolev spaces - Nečas [Neč67]: Systematic treatment 2. **Lipschitz domains**: - Grisvard [Gri11]: Elliptic problems on Lipschitz domains - Marschall [Mar87]: Trace and extension under minimal regularity assumptions 3. **BV and Besov spaces**: - Dávila [Dáv02]: Traces of BV functions - Malý-Shanmugalingam-Snipes [MSS18]: Traces and extensions of BV functions ### Improvements over Prior Work Compared to [DK19]: - Introduce measure $\mu_s$ achieving parameter robustness - Recover classical trace theorem (as $s \to 1^-$) Compared to [BGPR20]: - Handle general $p \geq 1$ - Explicitly give norm of trace space - Prove limiting behavior Compared to [GH22]: - Extend to general Lipschitz domains (not just $C^{1,1}$) - Handle full range $p \geq 1$ (not just $p=2$) - Different extension construction (Whitney decomposition vs. Poisson extension) ## Conclusions and Discussion ### Main Conclusions 1. **Complete trace and extension theory**: Establish complete trace and extension theory for energy space $V^{s,p}(\Omega|\mathbb{R}^d)$ applicable to all $p \geq 1$ and bounded Lipschitz domains 2. **Parameter robustness**: All results are robust with respect to $s \in (s_\star, 1)$, with continuity constants depending only on lower bounds of $s$ 3. **Recovery of classical theory**: As $s \to 1^-$, nonlocal trace spaces converge to classical trace spaces: - $p>1$: $T^{s,p}(\Omega^c) \to W^{1-1/p,p}(\partial\Omega)$ - $p=1$: $T^{s,1}(\Omega^c) \to B_1^{0,1}(\partial\Omega)$ (Besov space) 4. **Application value**: Provides natural variational framework for fractional p-Laplacian and other nonlocal problems ### Limitations 1. **Domain regularity**: Requires bounded Lipschitz domains; not applicable to rougher domains (e.g., fractal boundaries) 2. **Peculiarities of $p=1$ case**: - Extension operator norm depends on lower bounds of $1-s$ when $d=1$ - Cannot directly embed into $T^{s,1}(\Omega^c)$ (requires Besov regularity), see Remark 3.11 3. **Measure choice**: Specific form of measure $\mu_s$ depends on problem structure; may require adjustment for other types of nonlocal operators 4. **Optimality of constants**: While continuity is proven, optimality of constants is not discussed 5. **Higher-order case**: Only treats $s \in (0,1)$; does not include $s \geq 1$ ### Future Directions 1. **Nonlocal Neumann problems**: Authors mention that Neumann problems can be studied within this framework with limiting behavior analysis 2. **More general operators**: Extend to more general Lévy measures and nonlinear operators (Remark 1.9) 3. **Numerical methods**: Develop numerical methods based on these theoretical results 4. **Rough domains**: Study corresponding theory on fractal or rougher boundaries 5. **Higher-order spaces**: Study case $s > 1$ ## In-Depth Evaluation ### Strengths 1. **Theoretical completeness**: - Provides complete trace and extension theory, filling important gap in nonlocal analysis - Unified treatment of all cases $p \geq 1$ - Explicit dependence of constants on parameters 2. **Technical innovations**: - Clever design of measure $\mu_s$ simultaneously achieving boundary behavior and decay at infinity - Creative combination of Bessel potential interpolation theory and Whitney extension technique - Careful treatment of $p=1$ case (Hardy inequalities, BV spaces) 3. **Mathematical rigor**: - All theorems have complete proofs - Dependence of constants on parameters is explicit - Multiple technical details handled carefully (e.g., Lemmas 3.7, 3.8) 4. **Continuity and robustness**: - Parameter robustness with respect to $s$ is core contribution - Recovery of classical theory (as $s \to 1^-$) provides theoretical consistency check - Important for applications (e.g., numerical methods) 5. **Writing quality**: - Clear structure, progressing from motivation to applications - Detailed literature review - Well-organized technical lemmas ### Weaknesses 1. **Readability challenges**: - Proof techniques are highly technical, requiring deep background in function space theory - Some key steps (e.g., proof of Proposition 3.10) are very complex - Lack of intuitive explanations or geometric illustrations 2. **Application limitations**: - No numerical experiments validating theory - Few examples of application to concrete problems (only Corollary 1.7) - Explicit constant values unknown (may affect numerical implementation) 3. **Technical limitations**: - Additional dependence on $1-s$ for $p=1$, $d=1$ (Theorem 1.3) indicates theory is not fully unified - Counterexample 3.11 shows certain natural conjectures fail, but lacks deep discussion 4. **Generalizability**: - Limited to scalar problems; vector-valued case not addressed - Extension to time-dependent problems unclear ### Impact 1. **Theoretical contribution**: - Provides fundamental tools for nonlocal analysis - Solves long-standing problem in field (robust trace theory) - Expected to become standard reference in field 2. **Practical value**: - Provides theoretical foundation for numerical methods for fractional PDEs - Valuable for peridynamics and related applications (Remark 1.1) - Well-posedness results (Corollary 1.7) directly applicable 3. **Reproducibility**: - All theorems are constructive - Extension operator has explicit definition (equation 4.11) - Though technically complex, in principle implementable 4. **Follow-up research**: - Already has follow-up work (e.g., Neumann problems) - Provides template for studying other nonlocal operators - May inspire new numerical method designs ### Applicable Scenarios 1. **Theoretical research**: - Variational theory of fractional PDEs - Functional analysis of nonlocal operators - Interpolation theory of function spaces 2. **Numerical analysis**: - Theoretical foundation for finite element methods - Treatment of boundary conditions - Error estimation 3. **Application domains**: - Anomalous diffusion - Image processing (nonlocal filtering) - Peridynamics (nonlocal solid mechanics) - Mathematical finance (Lévy processes) 4. **Not applicable to**: - Numerical calculations requiring explicit constants (constants not given) - Very rough domains (beyond Lipschitz class) - High-dimensional problems (constants may degenerate as $d \to \infty$) ## References **Key references**: 1. [BBM01] Bourgain-Brezis-Mironescu: Convergence $V^{s,p}(\Omega|\mathbb{R}^d) \to W^{1,p}(\Omega)$ 2. [DK19] Dyda-Kassmann: Function spaces and extensions for nonlocal Dirichlet problems 3. [JW84] Jonsson-Wallin: Function spaces on subsets of $\mathbb{R}^n$ 4. [Gag57] Gagliardo: Classical trace theorem 5. [GH22] Grube-Hensiek: Robust trace theory for $p=2$ and $C^{1,1}$ domains 6. [BGPR20] Bogdan et al.: Extensions and traces for nonlocal operators 7. [FK22] Foghem-Kassmann: General framework for nonlocal Neumann problems on domains 8. [Fog23b] Foghem: Stability of complement problem for p-Lévy operators --- **Overall Assessment**: This is a high-quality theoretical mathematics paper solving important problems in nonlocal analysis. The design of measure $\mu_s$ and robustness results are core innovations. While highly technical, it provides fundamental tools for the field and is expected to have broad impact. Main limitations are lack of numerical validation and limited discussion of certain technical peculiarities (e.g., $p=1$, $d=1$ case).