Abstract
We give Hardy–Stein and Douglas identities for nonlinear nonlocal Sobolev–Bregman integral forms with unimodal Lévy measures. We prove that the corresponding Poisson integral defines an extension operator for the Sobolev–Bregman spaces. As an application, we obtain the boundedness of the Dirichlet-to-Neumann operator on weighted \(L^p\) spaces. We also show that the Poisson integrals are quasiminimizers of the Sobolev–Bregman forms.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
In 1931 Douglas [25] established a connection of the energy of the harmonic function u on the unit disc B(0, 1) with the “energy” of its boundary trace g, regarded as a function on \([0,2\pi )\):
The formula arose in the study of the so-called Plateau problem—the problem of existence of minimal surfaces posed by J.-L. Lagrange. The identity holds true provided that the left-hand side is finite—for details see, e.g., Chen and Fukushima [15, (2.2.60)]. Thus, under the integrability condition, (1.1) is valid for the solutions of the Dirichlet problem,
In our paper we propose a variant of (1.1), which we call nonlinear nonlocal Douglas identity. The term “nonlocal” means that the Laplace operator \(\Delta \) above is replaced by a nonlocal operator L. Specifically, we adopt the following setting. Let \(d=1,2,\ldots \). Suppose that the function \(\nu :[0,\infty )\rightarrow (0,\infty ]\) is nonincreasing and, with a slight abuse of notation, let \(\nu (z)=\nu (|z|)\) for \(z\in {{{\mathbb {R}}}^d}\). In particular, \(\nu \) is symmetric, i.e., \(\nu (z)=\nu (-z)\), \(z\in {\mathbb {R}}^d\). Assume further that
Thus, \(\nu \) is a strictly positive density function of an infinite isotropic unimodal Lévy measure on \({{{\mathbb {R}}}^d}\) (in short, \(\nu \) is unimodal). For \(u:{{{\mathbb {R}}}^d}\rightarrow {{\mathbb {R}}}\) and \(x\in {{{\mathbb {R}}}^d}\) we let
Here, \(\nu (x,y):=\nu (y-x)\), and the limit exists, e.g., for u in \(C_c^\infty ({\mathbb {R}}^d)\), the smooth functions with compact support. Operators of the form (1.3) are called nonlocal, because the value of \(L \phi (x)\) also depends on the values of \(\phi \) outside of a neighborhood of x. Furthermore, the operators satisfy the maximum principle, meaning that if \(\phi (x_0)=\sup \{\phi (x): x\in {{{\mathbb {R}}}^d}\}\), then \(L\phi (x_0)\le 0\). It is well known that such operators may be used to describe transportation of mass, charge, etc. in elliptic and parabolic equations; especially to pose boundary-value problems.
To our nonlocal setting we bring a judicious way of measuring the smoothness of functions for a given set. Let \(D\subset {{{\mathbb {R}}}^d}\) be open. For the sake of gradual introduction we first consider the quadratic form
Such forms appeared in Servadei and Valdinoci [55, 56], where the set \({{\mathbb {R}}}^d\times {{\mathbb {R}}}^d\setminus D^c\times D^c\) was denoted Q, then in Ros-Oton [52, (3.1)] and Dipierro, Ros-Oton and Valdinoci [24, p. 379]. Similar forms were also used in Felsinger, Kassmann and Voigt [29, Definition 2.1 (ii)]. \( \mathcal {E}_{D} \) is the energy functional of the nonlocal Dirichlet problem
see [55, 56] and Bogdan, Grzywny, Pietruska-Pałuba and Rutkowski [8]. It should be noted that \( \mathcal {E}_{D} \) is better than the vanilla form \( \mathcal {E}_{{\mathbb {R}}^d} \) for solving (1.5), because it allows for more general external conditions g due to the restriction of integration in (1.4) to \(Q={{\mathbb {R}}}^d\times {{\mathbb {R}}}^d{\setminus } D^c\times D^c\), cf. [8, p. 39]. Therefore \( \mathcal {E}_{D} \) constitutes an important step forward in nonlocal variational problems; we refer the reader to [8] for more details and to [55, 56] for applications to nonlinear equations. We note that our results also have consequences for the Dirichlet problem for L on D when \( \mathcal {E}_{{\mathbb {R}}^d} \) is used, see Corollary 4.4 below.
Numerous papers study the nonlocal Dirichlet problem by variational methods for nonlocal operators — in the present setting we should note [29, 52], and Rutkowski [53]. It is known for many Lévy and Lévy-type kernels \(\nu \) and bounded D [29, 53, 8, Section 5] that a unique weak solution of (1.5) exists provided that \(g:D^c\rightarrow {{\mathbb {R}}}\) can be extended to a function \(u\in L^2(D)\) from the Sobolev class
It is therefore important to determine conditions on g that allow for such an extension—in other words—to determine the trace space, say, \( \mathcal {X}_{D} \), of \( \mathcal {V}_{D} \). We note in passing that by [8, Lemma 3.4], the functions from \( \mathcal {V}_{D} \) are automatically square integrable on D. For the fractional Laplacian \(\Delta ^{\alpha /2}:=-(-\Delta )^{\alpha /2}\) (see Sect. 2.1 for a definition) a solution to this problem was proposed by Dyda and Kassmann [26] by using the Whitney decomposition and the method of reflection. In fact, [26, Theorem 3] concerns general p-increments, i.e., \(|u(x)-u(y)|^p\) with \(p\ge 1\).
In [8] we resolved the extension and trace problem for \(p=2\) for a wide class of unimodal Lévy operators by a different approach based on the (quadratic) nonlocal Douglas identity. Namely, [8, Theorem 2.3] asserts that the trace space \( \mathcal {X}_{D} \) consists of functions \(g:D^c\rightarrow D\) for which the following form on \(D^c\) is finite,
Here and afterwards we call
the kernel of interaction via D, or interaction kernel, and \(G_D\) is the Green function of L for D; see Sect. 2.1 for details. We note that \(\gamma _D\) is the nonlocal normal derivative of the Poisson kernel of L, see (6.1) and (2.8) below, similarly as the kernel in the classical Douglas identity, see Bogdan, Fafuła, and Rutkowski [7, Subsection 2.3]. The nonlocal Douglas identity of [8] can be stated as follows,
where \(g:D^c\rightarrow {{\mathbb {R}}}\), \( \mathcal {H}_{D} [g]<\infty \), and \(u= P_D[g]\) is the Poisson integral of g, see Sect. 2.1. Notably, \(P_D[g]\) is a harmonic function of L, so the identity (1.7) explains the energy of a harmonic function by the energy of its external values. In the language of Chen and Fukushima [15, Chapter 5], the right-hand side of (1.7) is the trace form and \(\gamma _D(z,w)\,\textrm{d}z\textrm{d}w\) is the Feller measure for \(( \mathcal {E}_{{\mathbb {R}}^d} , \mathcal {V}_{{\mathbb {R}}^d} )\) on \(D^c\), but the extension and trace problem for \( \mathcal {V}_{D} \) were not investigated in [15]. We also note that Jacob and Schilling [41] studied Douglas identities for nonlocal censored-type Dirichlet forms.
Our present goal is to extend the nonlocal Douglas formula (1.7) to a more general nonlinear case. The possibility of such a setting occurred to us owing to the recent Hardy–Stein identities of Bogdan, Dyda and Luks [6, Theorem 2]. To this end we will use the following notion, the French power:
More precisely, \(x^{\langle \kappa \rangle }=x^\kappa \) if \(x>0\), \(x^{\langle \kappa \rangle }=-|x|^\kappa \) if \(x<0\), and \(0^{\langle \kappa \rangle }=0\). For example, \(x^{\langle 0\rangle }={{\,\textrm{sgn}\,}}(x)\) and \(x^{\langle 2\rangle }\ne x^2\) as functions on \({{\mathbb {R}}}\). In what follows we fix \(1<p<\infty \), the exponent of the “nonlinearity” alluded to in the title of the paper. Our nonlinear nonlocal Douglas identity is as follows:
where \(u=P_D[g]\) and \(g:D^c\rightarrow {{\mathbb {R}}}\). For a precise statement see Theorem 4.1 and Remark 4.2 below, since the result hinges on suitable additional assumptions on \(\nu \), D and g. No analogue of (1.8) seems to exist in the literature for \(p\ne 2\), even for \(\Delta ^{\alpha /2}\). However, related nonlinear forms \(\int u^{\langle p-1\rangle }Lu\) appear often in the literature concerning Markovian semigroups of operators on \(L^p\) spaces, see also (2.27) and (7.3) below. This is because for \(p\in (1,\infty )\) the dual space of \(L^p\) is \(L^{p/(p-1)}\) and for \(u\in L^p\) we have \(u^{\langle p-1\rangle }\in L^{p/(p-1)}\), and \(\int |u|^p=\int |u^{\langle p-1\rangle }|^{p/(p-1)}=\int u^{\langle p-1\rangle }u\). Therefore in view of the Lumer–Phillips theorem, \(u^{\langle p-1\rangle }\) yields a linear functional on \(L^p\) appropriate for testing dissipativity of generators, see, e.g., Pazy [48, Section 1.4]. In this connection we note that Davies [18, Chapter 2 and 3] gives some fundamental calculations with forms and powers. For the semigroups generated by local operators we refer to Langer and Maz’ya [45] and Sobol and Vogt [57, Theorem 1.1]. Liskevich and Semenov [46] use the \(L^p\) setting to analyze perturbations of Markovian semigroups. For nonlocal operators we refer to Farkas, Jacob and Schilling [28, (2.4)], and to the monograph of Jacob [40, (4.294)].
The following variant of (1.8) is also true, see (2.21), (2.22), and (2.24) below,
The integrands in (1.9) come from the second order Taylor remainder of the convex function \(x\mapsto |x|^p\), see (2.13), which leads us to the notion of Bregman divergence; see Sect. 2.2, see also Bregman [12] for the original contribution or Sprung [58]. Bregman divergence is important for statistical learning, see Nielsen and Nock [47] or Frigyik, Gupta and Srivastava [31] and the references therein. The Bregman divergence based on the power function \(|x|^p\) defines the free energy functionals in the studies of Sobolev and Gagliardo–Nirenberg–Sobolev inequalities by Carrillo et al. [14, p. 71] and Bonforte, Dolbeault, Nazaret, and Simonov [11]. It also commonly appears in entropy inequalities, see, e.g., Wang [60].
The present paper indicates further uses of Bregman divergence in PDEs. As we show in Sect. 6, \(\gamma _D\) is the kernel of the Dirichlet-to-Neumann map (6.2) for L. Over the last few years, Dirichlet-to-Neumann map related to nonlocal operators was intensively studied in the context of inverse problems, see, e.g., [3, 17, 33, 34]. The forms in (1.8) are suitable for studying the Dirichlet-to-Neumann map as an operator in \(L^p\). In particular, using our Douglas identity we show that the normalized Dirichlet-to-Neumann operator (6.5) is bounded on a certain weighted \(L^p\) space. Results in this direction were obtained by Vondraček [59], 63 and Foghem and Kassmann [30] for \(p=2\), but even in this case our approach gives new insights.
As another motivation, we mention that the form on the left-hand side of (1.8) with \(D={\mathbb {R}}^d\) is appropriate for studying \(L^p\) properties of Markovian semigroups. For instance, it was used by Bogdan, Jakubowski, Lenczewska, and Pietruska-Pałuba [9] to characterize the contractivity on \(L^p({\mathbb {R}}^d)\) of the semigroups generated by the fractional Laplacian with Hardy-type potentials. The interested reader may find insights into the technique in [9, Lemma 7 and Proof of Theorem 3], or even in (2.27) and (7.3) below.
The paper is organized as follows. Section 2 contains definitions and basic facts. Subsection 2.1 introduces notions from the probabilistic potential theory and Sect. 2.2 introduces our nonlinear setting and novel Sobolev–Bregman spaces \(\mathcal {V}_{D}^{\, p} \) and \(\mathcal {X}_{D}^{\, p} \) defined by the condition of finiteness of the respective sides of (1.8). In (2.19) we collect in one place four (equivalent) approximations for our Bregman divergence, which appear in the literature. In Sect. 3 we generalize the Hardy–Stein identities of [6] and [8] to our present context. This is instrumental for the proof of the Douglas identity in Sect. 4. In Corollary 4.3 we conclude that the Poisson integral \(P_D\) and the restriction to \(D^c\) are the extension and trace operators between the Sobolev–Bregman spaces. In view toward applications in variational problems, in Sect. 5 we prove the Douglas formula with the remainder for the energy of sufficiently regular nonharmonic functions. We also show that harmonic functions are quasi-minimizers of the considered nonlinear nonlocal forms, but in general not minimizers. In Sect. 6 we apply our results for the analysis of the Dirichlet-to-Neumann operator in \(L^p\) for \(p\ge 2\). Finally, in Sect. 7 we give, for \(p\ge 2\), the following result for Poisson integrals \(u=P_D[g]\) and the more usual integral forms based on the p-increments of functions:
It follows that \(g\mapsto P_D[g]\) is an extension operator for nonlocal Sobolev-type spaces \(\mathcal {W}_{D}^{\, p} \), defined by the finiteness of the left-hand side. In the remainder of Sect. 7 we compare \(\mathcal {V}_{D}^{\, p} \) and \(\mathcal {W}_{D}^{\, p} \).
2 Preliminaries
All the considered functions, sets and measures are tacitly assumed to be Borel. When we write \(f\approx g\) (resp. \(f \lesssim g\)), we mean that there is a number \(c>0\), i.e. a constant, such that \((1/c) f(x)\le g(x)\le c f(x)\) (resp. \(f(x) \le c g(x)\)) for all arguments x. Important constants will be capitalized: \(C_1,C_2,\ldots ,\) and their values will not change throughout the paper.
2.1 Processes and potential-theoretic notions
Let L and \(\nu \) be as in the Introduction. Following [8], we additionally assume that:
- (A1):
-
\(\nu \) is twice continuously differentiable on \((0,\infty )\) and there is a constant \(C_1\) such that
$$\begin{aligned} |\nu '(r)|, |\nu ''(r)|\le C_1\nu (r), \qquad r>1. \end{aligned}$$ - (A2):
-
There exist constants \(\beta \in (0,2)\) and \(C_2>0\) such that
$$\begin{aligned} \nu (\lambda r)\le & {} C_2 \lambda ^{-d-\beta }\nu (r),\qquad 0<\lambda , r\le 1, \end{aligned}$$(2.1)$$\begin{aligned} \nu (r)\le & {} C_2 \nu (r+1), \qquad \quad r\ge 1. \end{aligned}$$(2.2)
A prominent representative of unimodal Lévy operators L is the fractional Laplacian \(\Delta ^{\alpha /2}:=-(-\Delta )^{\alpha /2}\). In this case we have \(\nu (x,y) = c_{d,\alpha }|y-x|^{-d-\alpha }\), where \(\alpha \in (0,2)\), \(x,y\in {{{\mathbb {R}}}^d}\), and
We refer the reader to Bogdan and Byczkowski [5], Di Nezza, Palatucci and Valdinoci [22], Garofalo [32], and Kwaśnicki [44] for more information on \(\Delta ^{\alpha /2}\). Clearly, \(\nu (r)=c_{d,\alpha }r^{-d-\alpha }\) satisfies both (A1) and (A2).
Our results depend in part on martingale properties of harmonic functions, so we introduce the Lévy process \((X_t,t\ge 0)\) on \({\mathbb {R}}^d\) whose generator is given by (1.3). Let
the Lévy–Khinchine exponent of \((X_t)\). Since \(\nu ({\mathbb {R}}^d) = \infty \), by Sato [54, Theorem 27.7] and Kulczycki and Ryznar [43, Lemma 2.5], the densities \(p_t(x)\) of \((X_t)\) are continuous on \({\mathbb {R}}^d\setminus \{0\}\) for \(t>0\), and satisfy
For \(t>0\) and \(x,\,y\in {{{\mathbb {R}}}^d}\) denote \(p_t(x,y)=p_t(y-x)\), the transition density of \((X_t)\) considered as Markov process on \({{{\mathbb {R}}}^d}\). Namely, for starting point \(x\in {{{\mathbb {R}}}^d}\), times \(0\le t_1<t_2<\ldots t_n\) and sets \(A_1, A_2,\ldots A_n\subset {{{\mathbb {R}}}^d}\) we let, as usual,
This determines \({{\mathbb {P}}^{x}}\), the distribution of the process \((X_{t})\) starting from x, and \({\mathbb {E}}^{x}\), the corresponding expectation. In the wording of [54, Section 11], \((X_t)\) is the symmetric Lévy process in \({{{\mathbb {R}}}^d}\) with \((0,\nu ,0)\) as the Lévy triplet. Without losing generality we actually assume that each \(X_t\) is the canonical projection \(X_t(\omega )=\omega (t)\) on the space of càdlàg functions \(\omega :[0,\infty )\rightarrow {{{\mathbb {R}}}^d}\). We will also use the standard complete right-continuous filtration \(({\mathcal {F}}_t, t\ge 0)\) to analyze \((X_t)\), see Protter [50, Theorem I.31]. In passing we recall that every Lévy process is a Feller process [50].
Let \(\emptyset \ne D\subset {{\mathbb {R}}}^d\) be an open set. The time of the first exit of X from D is, as usual,
The Dirichlet heat kernel \(p_t^D(x,y)\) is defined by Hunt’s formula, cf. Chung and Zhao [16, Chapter 2.2],
It is the transition density of the process \((X_t)\) killed upon exiting D, i.e.,
for integrable functions f. The Green function of D is the potential of \(p_t^D\):
and by Fubini–Tonelli we have
The Poisson kernel of D for L is defined by
With (A2) for bounded set D we easily see that for all \(x,y\in D\) and \(z\in D^c\) with \({{\,\textrm{dist}\,}}(z, D )\ge \rho >0\),
where comparability constants depend on \(\nu \), D and \(\rho \). Consequently, (2.4) implies
with the same proviso on comparability constants. Note that if D is bounded and \(x\in D\) is fixed, then \({\mathbb {E}}^x\tau _D\) is bounded by a positive constant, see Pruitt [51]. We further note that for \(w,z\in D^c\) the interaction kernel satisfies
Finally, the L-harmonic measure of D for \(x\in {{{\mathbb {R}}}^d}\) is, as usual,
the distribution of the random variable \(X_{\tau _D}\) with respect to \({\mathbb {P}}^x.\)
From the Ikeda–Watanabe formula (see, e.g., Bogdan, Rosiński, Serafin and Wojciechowski [10, Section 4.2]) it follows that \(P_D(x,z)\,\textrm{d}z\) is the part of \(\omega _D^x(\textrm{d}z)\) which results from the discontinuous exit from D (by a jump). Below, by suitable assumptions on D and \(\nu \), we assure that \(P_D\) is the density of the whole harmonic measure, that is
This is true, e.g., if D is bounded, \(\nu \) satisfies (A2), \(|\partial D| = 0\) and \(D^c\) has the property (VDC). The latter means that there is \(c>0\) such that for every \(r>0\) and \(x\in \partial D\),
Here, as usual, \(B(x,r)=\{y\in {{{\mathbb {R}}}^d}: |y-x|<r\}\). For the proof of (2.10) under the above conditions, see [8, Corollary A.2].
Observe that for \(U\subset D\) we have \(p^U\le p^D\) and \(G_U\le G_D\). Therefore, \(P_U(x,z)\le P_D(x,z)\) for \(x\in U\), \(z\in D^c\), and \(\gamma _U(z,w)\le \gamma _D(z,w)\) for \(z,w\in D^c\). These inequalities may be referred to as domain monotonicity. For \(g:D^c\rightarrow {{\mathbb {R}}}\) we define the Poisson extension of g:
and we call \(\int _{D^c} g(z)P_D(x,z)\,\textrm{d}z\) the Poisson integral, as long as it is convergent.
2.2 Function \(F_p\) and related function spaces
We depend on the two humble real functions:
Clearly, \(|x|^\kappa \) is symmetric, \(x^{\langle \kappa \rangle }\) is antisymmetric: \((-x)^{\langle \kappa \rangle }=-x^{\langle \kappa \rangle }\), and their derivatives obey
Recall that \(p>1\). We let
For instance, if \(p=2\), then \(F_2(a,b)=(b-a)^2\), and if \(p=4\), then \(F_4(a,b)=(b-a)^2(b^2+2ab+3a^2)\). As the second-order Taylor remainder of the convex function \(|x|^p\), \(F_p\) is nonnegative. In fact,
see [6, Lemma 6]. In particular, for \(p\ge 2\) we have
Recall that if X is a random variable with the first moment finite and \(a\in {{\mathbb {R}}}\), then
Here we do not exclude the case \({\mathbb {E}} X^2=\infty \), in which case both sides of (2.16) are infinite, hence equal. This variance formula has the following analogue for \(F_p\).
Lemma 2.1
Let \(p> 1.\) Suppose that X is a random variable such that \({\mathbb {E}}|X|<\infty \). Then,
-
(i)
\({\mathbb {E}}F_p({\mathbb {E}}X, X) ={\mathbb {E}}|X|^p-\left| {\mathbb {E}}X\right| ^p\ge 0,\)
-
(ii)
\({\mathbb {E}} F_p(a,X)=F_p(a,{\mathbb {E}}X) +{\mathbb {E}}F_p({\mathbb {E}}X,X) \ge {\mathbb {E}}F_p({\mathbb {E}}X,X),\quad a\in {{\mathbb {R}}},\)
-
(iii)
\({\mathbb {E}}F_p(a,X)= {\mathbb {E}} F_p(b, X) + F_p(a,b) +(pa^{\langle p-1\rangle }-pb^{\langle p-1\rangle })(b-{\mathbb {E}}X), \quad a,b\in {\mathbb {R}}.\)
Proof
The verification is elementary, but we present it to emphasize that the finiteness of the first moment suffices. We have
where \({\mathbb {E}}|X|^p=\infty \) is permitted, too. The expression in (i) is nonnegative by Jensen’s inequality or because \(F_p\) is nonnegative. For all \(a\in {\mathbb {R}}\) we have,
as claimed in (ii). Finally, for all \(a,b\in {\mathbb {R}}\) the right-hand side of (iii) is
which simplifies to the left-hand side of (iii). Needless to say, (ii) is a special case of (iii). \(\square \)
We next propose a simple lemma concerning the p-th moments of random variables, which is another generalization of (2.16).
Lemma 2.2
For every \(p\ge 1\) there exist constants \(0<c_p\le C_p\) such that for every random variable X with \({\mathbb {E}}|X|<\infty \) and every number \(a\in {\mathbb {R}},\)
Proof
If \({\mathbb {E}}|X|^p=\infty \), then all the sides of (2.17) are infinite. Otherwise, by convexity,
For the lower bound we make two observations: \(\left| {\mathbb {E}}X-a\right| ^p\le {\mathbb {E}}|X-a|^p\) (Jensen’s inequality), and
Adding the two, we get that \(|{\mathbb {E}}X-a|^p+{\mathbb {E}} |X-{\mathbb {E}}X|^p\le (1+2^p){\mathbb {E}}|X-a|^p.\) \(\square \)
The function \(F_p(a,b)\) is not symmetric in a, b, but the right-hand side of (2.14) is, so it is natural to consider the symmetrized version of \(F_p\), given by the formula:
We can relate \(H_p\) to a “quadratic” expression as follows.
Lemma 2.3
For every \(p>1\) we have \(F_p(a,b)\approx H_p(a,b)\approx (b^{\langle p/2\rangle } - a^{\langle p/2\rangle })^2\).
Proof
The first comparison follows from (2.14): we have \(F_p(a,b)\approx F_p(b,a)\), hence \(F_p\approx H_p\). As for the second statement, if either a or b are equal to 0, then the expressions coincide up to constants depending on p. If \(a,b\ne 0\), then \(a = tb\) with \(t\ne 0\). Using this representation we see that the second comparison is equivalent to the following:
The latter holds because both sides are continuous and positive except at \(t=1\); at infinity both are power functions with the leading term \(|t|^p\), and at \(t=1\) their ratio converges to a positive constant. \(\square \)
Summarizing, by (2.14) and Lemma 2.3 for each \(p\in (1,\infty )\) we have
It is hard to trace down the first occurrence of such comparisons in the literature. The one-sided inequality \(|b^{p/2}-a^{p/2}|^2\le \frac{p^2}{4(p-1)}(b-a)(b^{p-1}-a^{p-1})\) for \(a,b\ge 0\) can be found in connection with logarithmic Sobolev inequalities, e.g., in Davies [18, (2.2.9)] for \(2<p<\infty \), and Bakry [2, p. 39] for \(p>1\). The opposite inequality \((b-a)(b^{p-1}-a^{p-1})\le (b^{p/2}-a^{p/2})^2\) with \(a,b>0\) and \(p>1\) appears, e.g., in [46, Lemma 2.1].
In fact the following inequalities hold for all \(p\in (1,\infty )\) and \(a,b\in {{\mathbb {R}}}\):
Indeed, if a and b have opposite signs then it is enough to consider \(b=t\ge 1\) and \(a=-1\), and to compare \((t+1)(t^{p-1}+1)=t^p+t^{p-1}+t+1\) with \((t^{p/2}+1)^2=t^p+2t^{p/2}+1\). We have \(t^{p/2}= \sqrt{t^{p-1}t}\le (t^{p-1}+t)/2\), which verifies the left-hand side inequality in (2.20) with constant 1, which is better than \(4(p-1)/p^2\). We further get the right-hand side inequality in (2.20), and the constant 2 suffices, because \(t^{p-1}+t -(t^p+1)=(1-t)(t^{p-1}-1)\le 0\). Note that the constant 2 is not optimal for individual values of p, e.g., for \(p=2\), but the constant 1 does not suffice for \(p\in (1,2)\cup (2,\infty )\) because then \(1\vee (p-1) > p/2\), and so \(t^{p-1} + t > 2t^{p/2}\) for large t.
If a and b have the same sign, then we may assume \(b=ta\), \(a>0\), \(t\ge 1\), and consider the quotient
where \(s=\sqrt{t}\), \(h(s)= s(s^{p-2}-1)/(s^p-1)\). We see that h(s) is strictly positive for \(p>2,\) \(s>1\) and negative for \(p\in (1,2)\). We claim that it decreases in the former case and increases in the latter. The sign of the derivative of h is the same as the sign of the function \(l(s) = -s^{2p-2} + (p-1)s^p - (p-1)s^{p-2} + 1\). Now, since \(l(1) = 0\), the sign of l on \((1,\infty )\) is in turn equal to the sign of \(l'(s) = (p-1)s^{p-3}(-2s^p + ps^2 - (p-2))\), and further equal to the sign of \(-2p(s^{p-1} - s)\). Since the last function is negative on \((1,\infty )\) if \(p>2\) and positive for \(p\in (1,2)\), the claim is proved. Consequently, the function \(s\mapsto h(s)^{2}\) is decreasing on \((1,\infty )\), so we get
and (2.20) follows. The above also shows that the constant \(4(p-1)/p^2\) in (2.20) cannot be improved.
We would like to note that for \(p\ne 2\), \(F_p(a+t,b+t)\) is not comparable with \(F_p(a,b)\). Indeed, for \(a,r>0\) one has \(F_p(a,a+r)\approx r^2(a\vee (a+r))^{p-2} = r^2(a+r)^{p-2}\), which is not comparable with \(F_p(0,r)=r^2\) for large values of a. Here are one-sided comparisons of \(F_p(a,b)\) with the more usual p-increments, see, e.g., Zeidler [61, p. 503].
Lemma 2.4
If \(p\ge 2\) then \(F_p(a,b) \gtrsim |b-a|^p\), and if \(1<p\le 2\), then \(|b-a|^p\gtrsim F_p(a,b)\).
Proof
If \(a=b\), then the inequalities are trivial, so assume that \(a\ne b\) and consider the quotient
Both parts of the statements now follow from the inequality \(|b-a|^r\le 2^r(|a|\vee |b|)^r\), \(r>0.\)
\(\square \)
In analogy to (1.4) for \(u:{{{\mathbb {R}}}^d}\rightarrow {{\mathbb {R}}}\) we define
By the symmetry of \(\nu \) and (2.18),
Of course, \(\mathcal {E}_{D}^{(2)} = \mathcal {E}_{D} \). For \(D={{\mathbb {R}}}^d\) we have
Clearly, for \(p=2\) we retrieve the classical Dirichlet form of the operator L.
Let \(g:D^c\rightarrow {{\mathbb {R}}}\). To quantify the increments of g, we use the form:
The spaces \( \mathcal {V}_{D} \) and \( \mathcal {X}_{D} \) discussed in the Introduction lend themselves to the following generalizations:
and
We call them Sobolev–Bregman spaces, since they involve the Bregman divergence. Our development below indicates that \(\mathcal {V}_{D}^{\, p} \) and \(\mathcal {X}_{D}^{\, p} \) provide a viable framework for nonlocal nonlinear variational problems. In view of (2.22) for all \(u:{{{\mathbb {R}}}^d}\rightarrow {{\mathbb {R}}}\) we have
where
if the integral is well defined, which is the case in (2.27) for \(v=u^{\langle p-1\rangle }\). For clarity we also note that by (2.20), (2.22) and (2.24), we have the comparisons
and
for all \(u:{{{\mathbb {R}}}^d}\rightarrow {{\mathbb {R}}}\) and \(g:D^c\rightarrow {{\mathbb {R}}}\) with the comparability constants depending only on p. Below, however, we focus on genuine equalities.
3 Hardy–Stein identity
We first collect properties of harmonic functions that are needed in the proof of the identity (1.8). We mostly follow [8], so our presentation will be brief. We write \(U\subset \subset D\) if the closure of U is a compact subset of D.
Definition 3.1
We say that the function \(u:{\mathbb {R}}^d\rightarrow {\mathbb {R}}\) is L-harmonic (or harmonic, if L is understood) in D if it has the mean value property inside D, that is: for every open set \(U\subset \subset D\),
If \(u(x) = {\mathbb {E}}^x u(X_{\tau _D})\) for all \(x\in D\), then we say that u is regular harmonic.
In the above we assume that the expectations are absolutely convergent.
The strong Markov property of \((X_t)\) implies that if u is regular L-harmonic in D, then it is L-harmonic in D. By [8, Section 4], if u is L-harmonic in D, then \(u\in L^1_{loc} ({\mathbb {R}}^d) \cap C^2(D),\) Lu(x) can be computed pointwise for \(x\in D\) as in (1.3), and \(Lu(x)=0\) for \(x \in D\). We also note that the Harnack inequality holds for L-harmonic functions (see Grzywny and Kwaśnicki [39, Theorem 1.9]; the assumptions of that theorem follow from (A2)).
We will use the following Dynkin-type lemma, proven in our setting in [8, Lemma 4.11].
Lemma 3.2
Let the set \(U\subset \subset D\) be open and Lipschitz. If \(\int _{{\mathbb {R}}^d} |\phi (y)| (1\wedge \nu (y))\, \textrm{d}y< \infty \) and \(\phi \in C^2({\overline{U}})\), then \(L\phi \) is bounded on \({\overline{U}}\) and for every \(x\in {\mathbb {R}}^d\),
where the integrals converge absolutely.
The following Hardy–Stein formula extends [6, Lemma 8] and [8, Lemma 4.12], where it was proved, for the fractional Laplacian and \(p>1\), and for unimodal operators L and \(p=2\), respectively.
Proposition 3.3
If \(u:{\mathbb {R}}^d\rightarrow {\mathbb {R}}\) is L-harmonic in D, \(p >1 \), and \(U \subset \subset D\) is open Lipschitz, then
Proof
As a guideline, the result follows by taking \(\phi = |u|^p\) in the Dynkin formula (3.1). We combine the methods of [6] and [8]. By [8, Lemma 4.9] if u is harmonic in D, then \(u\in C^2(D)\). Thus, in particular, \(|u|^p\) is bounded in a neighborhood of \({\overline{U}}\). Let \(x\in U\). Consider the complementary cases:
Since \(|u|^p\) is bounded in a neighborhood of \({\overline{U}}\), this dichotomy can be reformulated as
see the end of the proof of [8, Lemma 4.11] and (2.6).
In case (i), we show that the right-hand side of (3.2) is infinite as well. Assume first that \(|u|>0\) on a subset of U of positive measure. Pick \(y\in U\) satisfying \(|u(y)|>0\), and let \(A = \{z\in {U^c}: |u(z)| \ge (2+\sqrt{2})|u(y)|\}\). Now, since \(x,y\in U\) are fixed and \(\nu \) is positive, continuous, and satisfies (2.2), we have \(\nu (x,z) \approx \nu (y,z)\) for \(z\in U^c\). Therefore, by (i),
as well. Furthermore,
and consequently we must have
By the definition of A, for \(z\in A\) we have
By (2.14) and (3.3) we therefore obtain
This is true for all points y in a set of positive Lebesgue measure, which proves that the right-hand side of (3.2) is infinite. If, on the other hand, \(u\equiv 0\) in U, then \(F_p(u(y),u(z)) = c|u(z)|^p\) for all \(z\in {\mathbb {R}}^d\), \(y\in U\), and by (i) the right-hand side of (3.2) is infinite again.
We now consider the case (ii). Thus \({\mathbb {E}}^x|u(X_{\tau _U})|^p<\infty \) and the integrability condition of Lemma 3.2 is satisfied for \(\phi = |u|^p\). We will first prove (3.2) for \(p\ge 2\). Then \(\phi \) is of class \(C^2\) on D, so we are in a position to use Lemma 3.2 and we get
The integral on the right-hand side is absolutely convergent. Furthermore, since u is L-harmonic,
Inserting this to (3.4) gives the statement.
When \(p\in (1,2)\), the function \({{\mathbb {R}}}\ni r\mapsto |r|^p\) is not twice differentiable, and the above argument needs to be modified. We work under the assumption (ii), and we follow the proof of [6, Lemma 3]. Consider \(\varepsilon \in {{\mathbb {R}}}\) and the function \({{\mathbb {R}}}^d\ni x\mapsto (x^2+\varepsilon ^2)^{p/2}\). Let
Since \(1<p<2\), by [6, Lemma 6],
Let \(\varepsilon >0\). We note that \((u^2+\varepsilon ^2)^{p/2}\in C^2(D)\). Also, the integrability condition in Lemma 3.2 is satisfied for \(\phi =(u^2+\varepsilon ^2)^{p/2}\) since it is satisfied for \(\phi =|u|^p\) by (ii), and
see also (1.2). Furthermore, \({\mathbb {E}}_{x}(u(X_{\tau _U})^2+\varepsilon ^2)^{p/2}<\infty \). As in the first part of the proof,
therefore by Lemma 3.2,
From the Dominated Convergence Theorem the left-hand side of (3.9) goes to \({\mathbb {E}}^x|u(X_{\tau _U})|^p<\infty \) as \(\varepsilon \rightarrow 0^+\). Of course, \(F_p^{(\varepsilon )}(a,b)\rightarrow F_p(a,b)\) as \(\varepsilon \rightarrow 0^+\). Furthermore, by Fatou’s lemma and (3.9),
By (3.6) and the Dominated Convergence Theorem, we obtain (3.2) for \(p\in (1,2)\). \(\square \)
As a consequence, we obtain the the Hardy–Stein identity for D, generalizing and strengthening [6, (16)] and [8, Theorem 2.1].
Proposition 3.4
Let \(p>1\) be given. If u is L-harmonic in D and \(x\in D\), then
If u is regular L-harmonic in D, then the left-hand side can be replaced with \({\mathbb {E}}^x |u(X_{\tau _D})|^p\).
Proof
As noted in [8, Remark 4.4], \(\{u(X_{\tau _U}), \, U\subset D\}\) is a martingale ordered by the inclusion of open subsets of D. By domain monotonicity of the Green function and the nonnegativity of \(F_p\), both sides of (3.2) increase if U increases. Since every open set \(U\subset \subset D\) is included in an open Lipschitz set \(U\subset \subset D\), the supremum in (3.10) may be taken over open Lipschitz sets \(U\subset \subset D\). The first part of the statement follows from the Monotone Convergence Theorem.
If additionally u is regular harmonic, then
This is delicate. Indeed, by [8, Remark 4.4], the martingale \(\{u(X_{\tau _U}), \, U\subset \subset D\}\) is closed by the integrable random variable \(u(X_{\tau _D})\). Therefore Lévy’s Martingale Convergence Theorem yields that \(u(X_{\tau _U})\) converges almost surely, and in \(L^1\) to a random variable Z, as \(U\uparrow D\), and we have \(Z = {\mathbb {E}}^x[u(X_{\tau _D})|\sigma (\bigcup _{U\subset \subset D}{\mathcal {F}}_{\tau _U})]\), see, e.g., Dellacherie and Meyer [21, Theorem 31 a,b, p. 26]. We claim that the \(\sigma \)-algebra \(\sigma (\bigcup _{U\subset \subset D}{\mathcal {F}}_{\tau _U})\) is equal to \({\mathcal {F}}_{\tau _D}\). Indeed, by Proposition 25.20 (i),(ii), and Proposition 25.19 (i),(ii) in Kallenberg [42, p. 501], the filtration of \((X_t)\) is quasi-left continuous. Therefore \(\tau _U\) increases to \(\tau _D\) as U increases to D, and our claim follows from Dellacherie and Meyer [20, Theorem 83, p. 136]. Consequently, \(Z = u(X_{\tau _D})\). Now, if \(\sup _{x\in U\subset \subset D} {\mathbb {E}}^x|u(X_{\tau _U})|^p < \infty \), then [21, Theorem 31 c, p. 26] yields (3.11). Else, if the supremum is infinite, then \({\mathbb {E}}^x |u(X_{\tau _D})|^p = \infty \) by Jensen’s inequality, and (3.11) holds, too. \(\square \)
We note in passing that the case \(p=2\) of (3.11) was stated for less general sets D in the first displayed formula following (5.2) in the proof of Theorem 2.3 in [8]. Accordingly, the proof in [8] was easier.
4 The Douglas identity
We now present our main theorem. It is a counterpart of (1.7) with square increments of the function replaced by “increments” measured in terms \(F_p\) or \(H_p.\)
Theorem 4.1
Let \(p>1\). Assume that the Lévy measure \(\nu \) satisfies (A1) and (A2), \(D\subset {\mathbb {R}}^d\) is open, \(D^c\) satisfies (VDC), and \(|\partial D| = 0\).
-
(i)
Let \(g:D^c\rightarrow {{\mathbb {R}}}\) be such that \(\mathcal {H}_{D}^{(p)} [g]<\infty .\) Then \(P_D[g]\) is well-defined and satisfies
$$\begin{aligned} \mathcal {H}_{D}^{(p)} [g]=\mathcal {E}_{D}^{(p)} [P_D[g]]. \end{aligned}$$(4.1) -
(ii)
Furthermore, if \(u:{\mathbb {R}}^d \rightarrow {\mathbb {R}}\) satisfies \(\mathcal {E}_{D}^{(p)} [u] < \infty \), then \(\mathcal {H}_{D}^{(p)} [u|_{D^c}] < \infty \).
Here, as usual, \(u|_{D^c}\) is the restriction of u to \(D^c\), but in what follows we will abbreviate:
and
Remark 4.2
The more explicit expression of the Douglas identity (1.8) stated in the Introduction follows from (4.1), (2.22) and (2.24).
To the best of our knowledge the present Douglas identities are completely new, and our approach is original. The proof of Theorem 4.1 is given below in this section.
Recall the space \(\mathcal {V}_{D}^{\, p} \), defined in (2.25), which is a natural domain of \(\mathcal {E}_{D}^{(p)} \), and the space \(\mathcal {X}_{D}^{\, p} \), defined in (2.26), which is a natural domain of \(\mathcal {H}_{D}^{(p)} \). From Theorem 4.1 we immediately obtain the following trace and extension result in the nonquadratic setting.
Corollary 4.3
Let \(\textrm{Ext }\, g=P_D[g]\), the Poisson extension, and \(\textrm{Tr}\, u=u|_{D^c}\), the restriction to \(D^c\). Then \(\textrm{Ext}:\mathcal {X}_{D}^{\, p} \rightarrow \mathcal {V}_{D}^{\, p} \), \(\textrm{Tr}:\mathcal {V}_{D}^{\, p} \rightarrow \mathcal {X}_{D}^{\, p} \), and \(\textrm{Tr}\,\textrm{Ext}\) is the identity operator on \(\mathcal {X}_{D}^{\, p} \).
It is well justified to call \(\textrm{Ext}\) the extension operator and \(\textrm{Tr}\) the trace operator for \(\mathcal {V}_{D}^{\, p} \).
We next give the Douglas identity for the Poisson extension on D and the form \(\mathcal {E}_{{\mathbb {R}}^d}^{(p)} \) (with the integration over the whole of \({{{\mathbb {R}}}^d}\times {{{\mathbb {R}}}^d}\)).
Corollary 4.4
If \(P_D[|g|]<\infty \) on D, in particular if \(\mathcal {H}_{D}^{(p)} [g] < \infty \), then
We note that the kernel on the right-hand side of the above identity also appears in [15, Theorem 5.6.3] for \(p=2\), but even the form \(\mathcal {E}_{D}^{(2)} \) and the Douglas identity of Theorem 4.1 with \(p=2\) on full domain \(\mathcal {V}_{D}^{\, 2} \) do not appear in [15].
The proof of Theorem 4.1 uses the following lemma, which asserts that the condition \(\mathcal {H}_{D}^{(p)} [g]<\infty \) implies the finiteness of \(P_D[|g|^p]\) and \(P_D[|g|]\) on D.
Lemma 4.5
Suppose that \(g:D^c\rightarrow {{\mathbb {R}}}\) satisfies \(\mathcal {H}_{D}^{(p)} [g]<\infty .\) Then for every \(x\in D\) we have \(\int _{D^c}|g(z)|^p P_D(x,z)\,\textrm{d}z<\infty \). In particular, the Poisson integral of g is well-defined.
Proof
Denote \(I=\int _{D^c}|g(z)|^p P_D(x,z)\,\textrm{d}z\). If \(\mathcal {H}_{D}^{(p)} [g]<\infty \), then
Since \(\nu > 0\), for almost all (hence for some) pairs \((w,x)\in D^c\times D\) we get
For the remainder of the proof, we only consider pairs (w, x) satisfying the above condition.
We will use different approaches for \(p\ge 2\) and \(p\in (1,2).\) Let \(p\ge 2\). From (2.15) we obtain
For \(z\in D^c\), let \(g_n(z)= -n\vee g(z) \wedge n\). Clearly \(|g_n(z)|\le |g(z)|\) and \(|g_n(z)|\nearrow |g(z)|\) when \(n\rightarrow \infty \). Since \(|g_n(z)|\le n, \) the integral \(I_n:=\int _{D^c}|g_n(z)|^p P_D(x,z)\,\textrm{d}z\) is finite. It is also true that the increments of \(g_n\) do not exceed those of g, that is \(|g_n(z)-g_n(w)|\le |g(z)-g(w)|\). Consequently,
The last inequality is obvious for \(p=2\), and follows from Jensen’s inequality if \(p>2\). Thus,
hence the sequence \((I_n)\) is bounded. By the Monotone Convergence Theorem, \(I <\infty \). By Jensen’s inequality we also get \(\int _{D^c}|g(z)| P_D(x,z)\,\textrm{d}z<\infty \). By the Harnack inequality, the finiteness of the Poisson integral of |g| or \(|g|^p\) at any point \(x\in D\) guarantees its finiteness at every point of D, see, e.g., [8, Lemma 4.6], therefore the proof is finished for \(p\ge 2\).
Now let \(p\in (1,2)\). If \(g\equiv 0\) a.e. on \(D^c\), then the statement is trivial. Otherwise, pick \(w\in D^c\) such that \(0<|g(w)|<\infty \). Let \( B=\{z\in D^c: |g(z)| > |g(w)|\}\). We have
Thus, \(P_D[|g|^p](x)<\infty \). The rest of the proof is the same as in the case \(p\ge 2\). \(\square \)
Proof of Theorem 4.1 To prove (i) we let \(\mathcal {H}_{D}^{(p)} [g]<\infty \) and we have (4.2). Let \(u = P_D[g]\). By (2.12), \(u=g\) on \(D^c\). By Lemma 4.5, u is well-defined, and it is regular L-harmonic in D, that is \({\mathbb {E}}^x[u(X_{\tau _D})]=u(x)\) for \(x\in D\), cf. Definition 3.1 and (2.10). In particular, we have \({\mathbb {E}}^x|u(X_{\tau _D})|<\infty \).
For \(x\in D\) consider the integral \(\int _{D^c}F_p(u(w),u(z)) P_D(x,z)\,\textrm{d}z\). By (2.10), \(P_D(x,z)\) is the density of the distribution of \(X_{\tau _ D}\) under \({\mathbb {P}}^x\), hence
By Lemma 2.1 (ii) applied to \(a=u(w)\), \(X=u(X_{\tau _D})\) and \({\mathbb {E}}={\mathbb {E}}^x\), the above expression is equal to
By integrating the first term on the right-hand side of (4.5) against \(\nu (x,w)\,\textrm{d}x\textrm{d}w\) we obtain
For the second term in (4.5) we use Lemma 2.1 (i) and Proposition 3.4:
We integrate the latter expression against \(\nu (x,w)\,\textrm{d}x\textrm{d}w\). By Fubini–Tonelli, (2.4) and (2.10),
Since the sum of (4.6) and (4.7) equals \(p\mathcal {E}_{D}^{(p)} [u]\), we obtain the Douglas identity.
We now prove (ii). It is not obvious how to directly conclude that \(\mathcal {E}_{D}^{(p)} [u] < \infty \) implies \(P_D[|u|]<\infty \) on D, thus we cannot apply Lemma 2.1. Instead we use another approach: by Lemma 2.3, \(\mathcal {E}_{D}^{(p)} [u] < \infty \) is equivalent to \( \mathcal {E}_{D} [u^{\langle p/2\rangle }] < \infty \). By the trace theorem for \(p=2\), see [8, Theorem 2.3], \( \mathcal {H}_{D} [u^{\langle p/2\rangle }] < \infty \). By Lemma 2.3 we get (ii). \(\square \)
5 Douglas and Hardy–Stein identities with remainders
Throughout this section we assume that D is bounded. In the (quadratic) case \(p=2\), under a mild additional assumption on D, the Poisson integral \(P_D[g]\) was shown to be the minimizer of the form \( \mathcal {E}_{D} \) among all Borel functions with a fixed exterior condition \(g\in \mathcal {X}_{D} \) (see [8, Proposition 5.4 and Theorem 5.5]). This needs not be the case when \(p\ne 2\), and in this section we give an example of D and \(g\in \mathcal {X}_{D}^{\, p} \) for which \(P_D[g]\) does not minimize \(\mathcal {E}_{D}^{(p)} \) among functions in \(\mathcal {V}_{D}^{\, p} \) equal to g on \(D^c\). However, \(P_D[g]\) is always a quasiminimizer, if we adopt the following definition:
Definition 5.1
Let \(K\ge 1.\) Function \(u\in \mathcal {V}_{D}^{\, p} \) is a K-quasiminimizer of \(\mathcal {E}_{D}^{(p)} \), if \(\mathcal {E}_{U}^{(p)} [u] \le K\mathcal {E}_{U}^{(p)} [v]\) for every nonempty open set \(U\subset D\) satisfying (VDC) and \(|\partial U| = 0\), and every \(v\in \mathcal {V}_{U}^{\, p} \) equal to u on \(U^c\). We say that u is a quasiminimizer if it is a K-quasiminimizer for some \(K\in [1,\infty )\).
The definition is inspired by the classical one given by Giaquinta and Giusti [35, (5.26)]. To avoid technical complications and to make this digression short we require regular test sets U above. However, to be prudent we note that the choice of admissible sets U may affect the definition of quasiminimizers and should be carefully considered, cf. Giusti [36, Example 6.5]. In the classical PDEs, quasiminimizers display many regularity properties similar to minimizers, see, e.g., Adamowicz and Toivanen [1], DiBenedetto and Trudinger [23], and Ziemer [62]. The main motivation for studying quasiminimizers is the fact that the solution of a complicated variational problem may be a quasiminimizer of a better understood functional see, e.g., [35, Theorem 2.1].
Proposition 5.2
Suppose that the assumptions of Theorem 4.1 are satisfied, D is bounded, and let \(g\in \mathcal {X}_{D}^{\, p} .\) Then \(P_D[g]\) is a K-quasiminimizer of \(\mathcal {E}_{D}^{(p)} \) with K independent of g.
Proof
Fix a subset \(U\subset D\) satisfying (VDC) and \(|\partial U| = 0\), and let \(v\in \mathcal {V}_{U}^{\, p} \) be equal to \(u:=P_D[g]\) on \(U^c\). According to (2.28) we have \(v^{\langle p/2\rangle }\in \mathcal {V}_{U} \) and
with constants independent of U and v. Note that \(v^{\langle p/2\rangle }\) agrees with \(u^{\langle p/2\rangle }\) on \(U^c\). Since \(U^c\) satisfies (VDC), by [8, Proposition 5.4 and Theorem 5.5],
By applying the Douglas identity for the set U, first with exponent 2, and then with exponent p, and by (2.29) we get that the right-hand side of (5.1) is equal to
In the last equality we use the identity \(P_U[u]= u\), see (2.10). The proof is complete. \(\square \)
To prove that Poisson integrals need not be minimizers, we first extend the Hardy–Stein and Douglas identities to functions that are not harmonic. The results are new even for \(p=2\) and \(\Delta ^{\alpha /2}\).
Recall that D is bounded, hence \({\mathbb {E}}^x\tau _D\) is bounded. In what follows by \(\lim \limits _{U\uparrow D}\) we denote the limit over an arbitrary ascending sequence of Lipschitz open sets \(U_n \subset \subset D\) such that \(\bigcup _{n} U_n = D\). Here is an extended version of the Hardy–Stein formula.
Proposition 5.3
Let \(p> 1\) and assume that \(\nu \) satisfies (A1) and (A2). Let \(u:{{{\mathbb {R}}}^d}\rightarrow {{\mathbb {R}}}\). If \(u\in C^2(D)\) and u and Lu are bounded in D, then for every \(x\in D\),
If in addition \(D^c\) satisfies (VDC) and \(|\partial D| = 0\), then \(\lim \limits _{U\uparrow D}{\mathbb {E}}^x |u(X_{\tau _U})|^p={\mathbb {E}}^x |u(X_{\tau _D})|^p\).
Proof
Let \(x\in D\). Since u, Lu, and \({\mathbb {E}}^x \tau _D\) are bounded on D, by (2.3) we get that the integral in (5.3) is finite. Therefore, using the arguments from the proof of Proposition 3.4, in what follows we may and do assume that \(\int _{{\mathbb {R}}^d} |u(x)|^p (1\wedge \nu (x))\, \textrm{d}x< \infty \), because otherwise both sides of (5.2) are infinite. With this in mind we first consider open Lipschitz \(U\subset \subset D\) so large that \(x\in U.\)
Let \(p\ge 2\). Since \(u\in C^2(D),\) we get that \(L|u|^p(x)\) and \({\mathbb {E}}^x|u(X_{\tau _{U}})|^p\) are finite for \(x\in U\), and (3.4) holds. Furthermore, since Lu is finite in D, the following manipulations are justified for \(y\in D\):
Consequently, (3.4) takes on the form
For clarity we note that the left-hand side of (5.5) is finite and the integral in (5.6) is absolutely convergent, so the integral in (5.5) is finite as well.
For \(p\in (1,2)\) we proceed as in the proof of Proposition 3.3, that is, instead of \(|u(x)|^p\) we consider \(\varepsilon > 0\) and the function \(x\mapsto (u(x)^2 + \varepsilon ^2)^{p/2}\). We obtain (cf. (3.8) and (5.4)),
As in the proof of Proposition 3.4, the left-hand side tends to \({\mathbb {E}}^x |u(X_{\tau _U})|^p\) as \(\varepsilon \rightarrow 0^+\). Furthermore, since Lu and u are bounded in D, the integral in (5.8) converges to that in (5.6). Then we apply Fatou’s lemma and the Dominated Convergence Theorem to the integral on the right-hand side of (5.7) and we obtain (5.5) for \(p\in (1,2)\), too.
We let \(U\uparrow D\) in (5.5). By the boundedness of u and Lu in D, the integral in (5.6) tends to the one in (5.3), which is absolutely convergent. The integral on the right-hand side of (5.5) converges to the one on the right-hand side of (5.2) by the domain monotonicity and the Monotone Convergence Theorem. Since the limit on the right-hand side of (5.2) exists, the limit on the left-hand side must exist as well. This proves (5.2).
If \(D^c\) satisfies (VDC) and \(|\partial D| = 0\), then (2.10) holds true. Furthermore, we have
The first term on the right converges to 0 by the boundedness of u on D and the fact that \({\mathbb {P}}^x(\tau _U\ne ~\tau _D)\) decreases to 0 as \(U\uparrow D\) (see the remark preceding (2.10); see also the proof of Lemma 17 in Bogdan [4] and the proof of Lemma A.1 in [8]). The second term converges to \({\mathbb {E}}^x |u(X_{\tau _D})|^p\) by the Monotone Convergence Theorem. Thus the left-hand side of (5.5) tends to \({\mathbb {E}}^x |u(X_{\tau _D})|^p\). \(\square \)
We next provide a Douglas-type identity for a class of nonharmonic functions:
Theorem 5.4
Suppose that the assumptions of Theorem 4.1 hold with the addition that D is bounded. Let \(u:{{{\mathbb {R}}}^d}\rightarrow {{\mathbb {R}}}\) be bounded, \(u\in C^2(D)\), and Lu be bounded in D. Then
where
Proof
Since u is bounded on \({\mathbb {R}}^d\), we have \(\int _{{\mathbb {R}}^d} |u(x)|(1\wedge \nu (x))\, \textrm{d}x< \infty \).
Assume first that \(\mathcal {H}_{D}^{(p)} [u] < \infty \). From Theorem 4.1 we have
By (2.8) and Fubini–Tonelli,
We apply Lemma 2.1 (iii) to \(a=u(w)\), \(b=u(x)\), with \(w\in D^c\) and \(x\in D\), \(X=u(X_{\tau _D})\), and \({\mathbb {E}}={\mathbb {E}}^x\). Note that \({\mathbb {E}}X=P_D[u](x)\). This yields:
After integration, we obtain
Note that every term above is finite. Indeed, by the boundedness of u,
To prove that this is finite, let \(v = u -P_D[u]\). We have \(Lv = Lu = f\in L^\infty (D)\) and \(v= 0\) on \(D^c\). Note that \(v\in C^2(D)\) and \(\int _{{\mathbb {R}}^d} |v(x)|(1\wedge \nu (x))\, \textrm{d}x< \infty \), cf. [8, Lemma 3.6]. Let \(U\subset \subset D\). By Lemma 3.2,
Since u is bounded on \({\mathbb {R}}^d\), we have \({\mathbb {E}}^x u(X_{\tau _U}) \rightarrow {\mathbb {E}}^x u(X_{\tau _D}) = P_D[u](x)\) as \(U\uparrow D\), cf. the last part of the proof of Proposition 5.3. Hence, the boundedness of f, the domain monotonicity, and the Dominated Convergence Theorem yield
This allows us to further estimate \(A_3\):
Since \(A_1(u)\) and \(A_2(u)\) are nonnegative, they must be finite as well, because \(\mathcal {H}_{D}^{(p)} [u] < \infty \). We then have
Thus, by Proposition 5.3 we obtain
where \(A_4(u)\) is the integral in (4.7). Note that \(A_2(u) + A_4(u) = p\mathcal {E}_{D}^{(p)} [u]\). Also, all the expressions in (5.10) are finite, see the discussion of \(A_3(u)\). To finish the proof of (5.9) in the case \(\mathcal {H}_{D}^{(p)} [u] < \infty \), we simply note that \(pA_D(u) = A_1(u) - A_4(u) + A_3(u)\).
The situation \(\mathcal {H}_{D}^{(p)} [u] = \infty \) remains to be considered. Since \(P_D[u]\) is bounded in D, by arguments similar to those in the estimates of \(A_3(u)\) above, we prove that \(A_D(u)\) is finite. Therefore by Theorem 4.1 the identity (5.9) holds with both sides infinite. \(\square \)
Knowing the form of the remainder \(A_D(u)\) in the Douglas identity (5.9), we may provide an example which shows that the Poisson integral need not be a minimizer of \(\mathcal {E}_{D}^{(p)} \) for \(p\ne 2\); it is only a quasiminimizer by Proposition 5.2.
Example 5.5
Let \(p>2\) and consider \(0<R<R_1\) such that \(D\subset \subset B_R\). Define
Since each \(g_n\) is bounded with support separated from D, we have \(g_n\in \mathcal {X}_{D}^{\, p} \cap \mathcal {X}_{D} \); see the discussion following Example 2.4 in [8]. By (2.6) there exists \(c>0\) such that
Furthermore, for every \(U\subset \subset D\) there is \(\epsilon >0\) such that
For \(x\in D\) we let
Obviously \(u_n\) are bounded on \({\mathbb {R}}^d\). We will verify that \(G_D[1] \in C^2(D)\). For this purpose we let f be a smooth, compactly supported, nonnegative function equal to 1 on D. By the Hunt’s formula and Fubini–Tonelli we get
Here G is either the potential kernel or the compensated potential kernel of \((X_t)\); see Grzywny, Kassmann and Leżaj [38, Appendix A] for details. In particular, by [38, Corollary A.3] and [54, Theorem 35.4] G is locally integrable, thus the first term in (5.13) is finite and smooth in D. Since the latter term in (5.13) is a harmonic function, we get that \(G_D[1]\in C^2(D)\). In particular, by [8, Lemma 4.10] and Dynkin [27, Lemma 5.7] we have \(Lu_n = -1\) in D. We are now in a position to apply Theorem 5.4. Fix open \(U\subset \subset D\). We get
We claim that \(A_D(u_n) > 0\) for large n. Indeed, recall that \(G_D[1](x)={\mathbb {E}}^x \tau _D\) is bounded. Since the integrals \(\int _{D^c} g_n(x)\, \textrm{d}x\) are bounded, by (5.11) there is \(M>0\) such that \({\mathbb {E}}^x u_n(X_{\tau _D}) < M\) for every \(x\in D\) and \(n\in {\mathbb {N}}\). Therefore the integral \(\int _{D{\setminus } U}\) in (5.14) is bounded from below, independently of n. Note that \(\int _{D^c}g_n(x)^{p-1}\textrm{d}x\rightarrow \infty \) as \(n\rightarrow \infty \). Thus, by (5.12) we obtain that \(\int _U\rightarrow \infty \) in (5.14) as \(n\rightarrow \infty \). Hence, for sufficiently large n we get that \(A_D(u_n) > 0\), which proves that \(\mathcal {E}_{D}^{(p)} [P_D[u_n]] > \mathcal {E}_{D}^{(p)} [u_n]\) for some n, as needed. The case \(p\in (1,2)\) may be handled similarly, by using \(g_n(z) = ((|z|-R)^{-1} \wedge n){\textbf {1}}_{B_{R_1}{\setminus } B_R}(z)\) and \(u_n = P_D[g_n] - G_D[1]\).
6 Applications to Dirichlet-to-Neumann map
In this section we adopt the assumptions of Theorem 4.1. In addition, we assume that the set D is bounded and \(p\in [2,\infty )\). We define the nonlocal normal derivative as an analogue of the fractional version of Dipierro, Ros-Oton, and Valdinoci [24, (1.2)], see also Vondraček [59]:
Note that the increments of f are integrated on D, but the integral is evaluated for \(z\in D^c\), if convergent. For instance, if \(f\in L^1({\mathbb {R}}^d,1\wedge \nu )\), then \({\mathcal {N}}\!f\in L^1_{loc}(D^c)\).
Assume that \(g\in \mathcal {X}_{D}^{\, p} \). Then \(u=P_D[g]\) solves the Dirichlet problem (1.5). By definition, the Dirichlet-to-Neumann operator \({D\! N}\) maps the exterior condition g to the nonlocal normal derivative \(h:= {\mathcal {N}}u\). So, u solves the Neumann problem
and \({D\! N}:= {\mathcal {N}}\circ P_D\) on \(\mathcal {X}_{D}^{\, p} \). In fact, for almost every \(z\in D^c\),
where we have used the definition of \(\gamma _D\), the fact that \(u=P_D[g]\), and the Fubini–Tonelli theorem (justified by the estimates in the proof of Proposition 6.1). For \(z\in \textrm{Int}(D^c)={{{\mathbb {R}}}^d}\setminus {\overline{D}}\) we let
For example, for \(L = \Delta ^{\alpha /2}\) and (bounded) D of class \(C^{1,1}\), with \(\delta _D(z):= d(z,D)\) we have
Back to general L, we note that sharp estimates of \(\gamma _D\) are known for bounded \(C^{1,1}\) domains and the half-space, see [8, Theorem 2.6 and 6.1]. We next define the normalized Dirichlet-to-Neumann operator, for \(g\in \mathcal {X}_{D}^{\, p} \) and a.e. \(z\in D^c\),
In what follows we give several results for the Dirichlet-to-Neumann operator on \(L^p\). In particular, we show that \({D\! N}\) is well-defined: \(\mathcal {X}_{D}^{\, p} \rightarrow L^p(D^c,m^{1-p})\) and \({\widetilde{{D\! N}}}\) is bounded on \(L^p(D^c,m)\). We also relate the form \(\mathcal {H}_{D}^{(p)} \) to the operator \({D\! N}\) in (6.9).
Proposition 6.1
Assume that \(g\in \mathcal {X}_{D}^{\, p} \). Then \({D\! N}g\in L^p(D^c,m^{1-p})\) and \({\widetilde{{D\! N}}} g\in L^p(D^c,m)\). Furthermore, there exists a constant C, independent of g, such that
Proof
Using (6.3) and Jensen’s inequality we get
Since \(p\ge 2\), we have \(|a-b|^p \le (a-b)^2(|a|^{p-2} + |b|^{p-2})\). So, by (2.19),
which ends the proof. \(\square \)
Proposition 6.2
If \(g\in L^p(D^c,m)\), then \(g\in \mathcal {X}_{D}^{\, p} \) and there is \(C>0\), independent of g, such that
Proof
Following [30, Remark 2.37], we let \({\widetilde{g}}\) be the function g extended to D by 0. Then,
In particular, \({\widetilde{g}}\in \mathcal {V}_{D}^{\, p} \). By Proposition 5.2 we get that there exists a constant C, independent of g, such that \(\mathcal {E}_{D}^{(p)} [P_D[g]] \le C\mathcal {E}_{D}^{(p)} [{\widetilde{g}}]\). Using this, Theorem 4.1, and (6.7), we find that
which proves the result. \(\square \)
Corollary 6.3
The normalized Dirichlet-to-Neumann map \({\widetilde{{D\! N}}}\) is bounded on \(L^p(D^c,m)\).
The following is an analogue of the formula (7.3) below.
Proposition 6.4
Let \(f\in L^p(D^c,m)\) and \(g\in \mathcal {X}_{D}^{\, p} \). Then \(\int _{D^c} |{D\! N}g(z)||f(z)|^{p-1}\, \textrm{d}z< \infty \) and
Furthermore, if \(g\in L^p(D^c,m)\), then
Proof
By Hölder’s inequality with exponents p and \(p' = \frac{p}{p-1}\), and by Proposition 6.1,
It suffices to prove (6.8). By the symmetry of \(\gamma _D\),
The above application of the Fubini–Tonelli theorem is justified by using Hölder’s inequality with exponents p and \(p/(p-1)\), and (6.6); see also (6.4). The first and the last lines above yield (6.8). \(\square \)
Let us discuss related results for \(p=2\). In [59, Proposition 3.2], Vondraček shows that the normalized Dirichlet-to-Neumann operator map is bounded on \(L^2(D^c,m)\); our Corollary 6.3 extends this result to \(L^p\). As observed by Foghem and Kassmann [30, Remark 2.37], the space \(L^2(D^c,m)\) can be smaller than the trace space \( \mathcal {X}_{D} \). In [30, Section 4.4], the authors investigate the Dirichlet-to-Neumann operator for the equation \(Lu = \lambda u + f\), where \(\lambda \in {\mathbb {R}}\) is not a Dirichlet eigenvalue of L in D. They prove the boundedness of the Dirichlet-to-Neumann operator from the trace space into its dual. If we let \({D\! N}_{FK}\) be the Dirichlet-to-Neumann operator defined in [30] for \(\lambda = 0\) and \(f=0\), then using our Douglas identity and \({\widetilde{v}} = u_g = P_D[g]\) in [30, Definition 4.18], for \(g\in \mathcal {X}_{D} \) we get
Here \(\langle \cdot ,\cdot \rangle \) is the pairing between \( \mathcal {X}_{D} \) and its dual, see [30, Section 2.6]. Then, by polarization,
for \(g_1,g_2\in \mathcal {X}_{D} \). Both (6.3) and (6.10) give explicit integral representations for the Dirichlet-to-Neumann operator, which are more direct than (6.2). They were not stated in [30, 59], although similar formulas appear in [59, Section 3 and (4.2)] and the author of [59] was probably aware of the explicit versions.
On an informal level, (6.3) and (6.10) mean that the Dirichlet-to-Neumann map is the negative of the Lévy-type operator on \(D^c\) with jump kernel \(\gamma _D\), and \( \mathcal {H}_{D} \) is the corresponding Dirichlet form. Despite being smaller than \(\mathcal {X}_{D}^{\, 2} \), the space \(L^2(D^c,m)\), used by Vondraček, is suitable for studying the (negative of the) normalized Dirichlet-to-Neumann operator \({\widetilde{{D\! N}}}\) as a generator of a Markov process on \(D^c\). In fact, \(-{\widetilde{{D\! N}}}\) is the generator of the so-called trace process, see [59, (4.2)]. In this connection, the reader may compare (6.9), \({D\! N}\) and \(\mathcal {H}_{D}^{(p)} \) with (7.3), \(-L\) and \(\mathcal {E}_{D}^{(p)} \); see also [9, Lemma 7] for a detailed discussion of \(\mathcal {E}_{{{{\mathbb {R}}}^d}}^{(p)} \) for \(L=\Delta ^{\alpha /2}\).
7 Further discussion
As usual, D is a nonempty open set in \({{{\mathbb {R}}}^d}\). We define
and
Proposition 7.1
If \(p\ge 2\) then (1.10) holds true under the assumptions on D and \(\nu \) from Theorem 4.1, and the Poisson extension acts from \(\mathcal {Y}_{D}^{\, p} \) to \(\mathcal {W}_{D}^{\, p} \).
Proof
Assume that \(g\in \mathcal {Y}_{D}^{\, p} ,\) i.e., the right-hand side of (1.10) is finite. By a simple modification of the proof of [8, Lemma 4.6] we get that \(g\in L^p(D^c, P_D(x,z)\,\textrm{d}z)\) for every \(x\in D\), in particular the Poisson integral \(P_D[g](x)\) converges absolutely. By (2.8), the right-hand side of (1.10) equals
We use Fubini–Tonelli and consider the integral
By Lemma 2.2 we get that for \(x\in D\) and \(w\in D^c\),
We apply Proposition 3.4, to \({\widetilde{u}}(z):=u(z)-u(x).\) It is L-harmonic on D and \({{\widetilde{u}}}(x)=0\), therefore
For \(p\ne 2\) it is not true that \(F_p(a+t,b+t)\) is comparable with \(F_p(a,b)\), but since \(p\ge 2\), by Lemma 2.4 we have \(F_p(a+t,b+t)\ge c |a+t-b-t|^p = c|a-b|^p\). It follows that
and thus
We integrate the inequality on \(D^c\times D\) against \(\nu (w,x)\,\textrm{d}w\textrm{d}x\) as in (4.7), and the right-hand side is
The result follows. \(\square \)
We remark that in general (1.10) fails for \(p\in (1,2)\); see Lemma 7.4 and Example 7.5.
In the remainder of this section we compare \(\mathcal {W}_{D}^{\, p} \) and \(\mathcal {V}_{D}^{\, p} \), see (2.25), by using \(C_c^\infty ({\mathbb {R}}^d)\).
Lemma 7.2
For every \(p>1\) we have \(C_c^\infty ({\mathbb {R}}^d)\subseteq \mathcal {V}_{{\mathbb {R}}^d}^{\, p} \subseteq \mathcal {V}_{D}^{\, p} \).
Proof
The inclusion \(\mathcal {V}_{{\mathbb {R}}^d}^{\, p} \subseteq \mathcal {V}_{D}^{\, p} \) follows from the definition. To prove that \(C_c^\infty ({\mathbb {R}}^d)\subseteq \mathcal {V}_{{\mathbb {R}}^d}^{\, p} \), we let \(\phi \in C_c^\infty ({\mathbb {R}}^d)\). We have
It follows that \(L\phi \) is bounded on \({\mathbb {R}}^d\), cf. (1.3) and (1.2). Thus,
Furthermore, by the Dominated Convergence Theorem and the symmetry of \(\nu \),
By Fubini’s theorem, the substitutions \(z\rightarrow -z\) and \(x\rightarrow x+z\), and the symmetry of \(\nu \),
for every \(\epsilon > 0\). By (2.23), the Monotone Convergence Theorem and the above,
The result follows from (7.2) and (2.25). \(\square \)
The inclusion \(C_c^\infty ({\mathbb {R}}^d)\subseteq \mathcal {V}_{D}^{\, p} \) indicates that the Sobolev–Bregman spaces will be useful in variational problems posed in \(L^p.\)
The situation with the spaces \(\mathcal {W}_{D}^{\, p} \) is more complicated. While for \(p\ge 2\) we have a result similar to that of Lemma 7.2, for \(p\in (1,2)\) it is not so. More precisely, we have the following two lemmas:
Lemma 7.3
For \(p\ge 2\) we have \(C_c^\infty ({\mathbb {R}}^d) \subseteq \mathcal {W}_{{\mathbb {R}}^d}^{\, p} \subseteq \mathcal {W}_{D}^{\, p} \).
Proof
For \(\phi \in C_c^\infty ({\mathbb {R}}^d)\) let \(K = \textrm{supp}\, \phi \). Then we have \(|\phi (x) - \phi (y)| = 0\) on \(K^c\times K^c\) and
It follows that \(\phi \in \mathcal {W}_{{\mathbb {R}}^d}^{\, p} \). The inclusion \(\mathcal {W}_{{\mathbb {R}}^d}^{\, p} \subseteq \mathcal {W}_{D}^{\, p} \) is clear from the definition of the spaces. \(\square \)
Lemma 7.4
Let \(p\in (1,2)\) and assume that for some \(r>0\) we have \(\nu (y) \gtrsim |y|^{-d-p}\) for \(|y|<r\). If \(u\in \mathcal {W}_{D}^{\, p} \) has compact support in \({\mathbb {R}}^d\) and vanishes on \(D^c\), then \(u\equiv 0\).
Results of this type are well-known for the spaces with integration over \(D\times D\), where D is connected. Brezis [13, Proposition 2] shows that any measurable function must be constant in this case; a simpler proof of this fact was given by De Marco, Mariconda and Solimini [19, Theorem 4.1]. Lemma 7.4 follows by taking \(\Omega = {\mathbb {R}}^d\) in the aforementioned results, but we present a different proof. Such facts also hold true in the context of metric spaces, see, e.g., Pietruska-Pałuba [49]. We will see in the proof of Lemma 7.4 that the result reduces to that with \(D={{{\mathbb {R}}}^d}\).
Proof of Lemma 7.4 We may assume that u is bounded, because the p-increments of \((0\vee u)\wedge 1\) do not exceed those of u. Thus, since u is compactly supported, we get that \(u\in L^p({\mathbb {R}}^d)\cap L^2({\mathbb {R}}^d)\). Let
The Hausdorff–Young inequality asserts that for \(u\in L^p({\mathbb {R}}^d)\) we have
where \(p'=\frac{p}{p-1}\), see, e.g., Grafakos [37, Proposition 2.2.16]. We estimate the left-hand side of (1.10) by using (7.4):
By (7.4), \(|{\widehat{u}}(\xi )|^{p'} \, \textrm{d}\xi \) is a finite measure on \({\mathbb {R}}^d\). As we have \(p/p' < 1\), by Jensen and Fubini–Tonelli,
Since \(|1 - e^{2\pi i\xi y}| \ge |\sin 2\pi \xi y|\) and \(\nu (y) \gtrsim |y|^{-d-p}\) for small |y|, the integral is infinite, unless \(u= 0\) a.e. in \({\mathbb {R}}^d\). \(\square \)
As a comment to Lemmas 7.2 and 7.4 we recall that \(\mathcal {V}_{D}^{\, p} \) is defined in terms of \(F_p\). When a is close to b then, regardless of \(p>1\), the Bregman divergence \(F_p(a,b)\) is of order \((b-a)^2\) rather than \(|b-a|^p\). Thus \(\mathcal {V}_{D}^{\, p} \) agrees with the Lévy measure condition (1.2) better than \(\mathcal {W}_{D}^{\, p} \) does.
The following example indicates that the scale of linear spaces \(\mathcal {W}_{D}^{\, p} \) may not be suitable for analysis of harmonic functions when \(p\le 2\):
Example 7.5
Let \(\nu \) and p be as in Lemma 7.4. Let \(B=B(0,1)\) and assume that D is bounded and \({{\,\textrm{dist}\,}}(D,B) > 0\). Then there is \(g\in \mathcal {Y}_{D}^{\, p} \) such that \(u:=P_D[g]\notin \mathcal {W}_{D}^{\, p} \), i.e.,
Let \(g(z)={\textbf {1}}_{B}(z)\) for \(z\in D^c\). Then \(g\in \mathcal {Y}_{D}^{\, p} \), cf. the arguments following [8, Example 2.4]. Clearly, u is bounded in D. By the positivity of \(P_D\) [39, Lemma 2.2], \(u(x) > 0\) for every \(x\in D\). Of course, B, \(D^c\setminus B = B^c\setminus D\) and D form a partition of \({{{\mathbb {R}}}^d}\). Therefore their Cartesian products partition \({{{\mathbb {R}}}^d}\times {{{\mathbb {R}}}^d}\); in fact also \(B^c\times B^c\) and \({\mathbb {R}}^d\times {\mathbb {R}}^d {\setminus } D^c\times D^c\) (see below). Since u vanishes on \(D^c{\setminus } B\), \(u(x)-u(y)\) vanishes on \((D^c{\setminus } B)\times (D^c{\setminus } B)\). It follows that
Define \({\widetilde{u}}=u\) on \(B^c\) and \({\widetilde{u}}=0\) on B. Then, \({{\widetilde{u}}}=u\) on D and \({{\widetilde{u}}} = 0\) on \(D^c\), and
By the boundedness of u, the boundedness of D and the separation of D and B, the last integral is finite. Furthermore, since \({\widetilde{u}}\) is not constant and vanishes on \(D^c\), the left-hand side is infinite by Lemma 7.4. Therefore the left-hand side of (7.6) is infinite, which yields (7.5).
References
Adamowicz, T., Toivanen, O.: Hölder continuity of quasiminimizers with nonstandard growth. Nonlinear Anal. 125, 433–456 (2015)
Bakry D.: L’hypercontractivité et son utilisation en théorie des semigroupes. In: Lectures on Probability Theory (Saint-Flour, 1992), volume 1581 of Lecture Notes in Math., pp. 1–114. Springer, Berlin (1994)
Bhattacharyya, S., Ghosh, T., Uhlmann, G.: Inverse problems for the fractional-Laplacian with lower order non-local perturbations. Trans. Am. Math. Soc. 374(5), 3053–3075 (2021)
Bogdan, K.: The boundary Harnack principle for the fractional Laplacian. Studia Math. 123(1), 43–80 (1997)
Bogdan, K., Byczkowski, T.: Potential theory for the \(\alpha \)-stable Schrödinger operator on bounded Lipschitz domains. Studia Math. 133(1), 53–92 (1999)
Bogdan, K., Dyda, B., Luks, T.: On Hardy spaces of local and nonlocal operators. Hiroshima Math. J. 44(2), 193–215 (2014)
Bogdan, K., Fafuła, D., Rutkowski, A.: The Douglas formula in \(L^p\). arXiv e-prints. arXiv:2207.07431 (2022)
Bogdan, K., Grzywny, T., Pietruska-Pałuba, K., Rutkowski, A.: Extension and trace for nonlocal operators. J. Math. Pures Appl. 9(137), 33–69 (2020)
Bogdan, K., Jakubowski, T., Lenczewska, J., Pietruska-Pałuba, K.: Optimal Hardy inequality for the fractional Laplacian on \(L^p\). J. Funct. Anal. 282(8):Paper No. 109395, 31 (2022)
Bogdan, K., Rosiński, J., Serafin, G., Wojciechowski. Ł.: Lévy systems and moment formulas for mixed Poisson integrals. In: Stochastic analysis and related topics, volume 72 of Progr. Probab., pp. 139–164. Birkhäuser/Springer, Cham (2017)
Bonforte, M., Dolbeault, J., Nazaret, B., Simonov, N.: Stability in Gagliardo–Nirenberg–Sobolev inequalities: flows, regularity and the entropy method. To appear in Memoirs AMS
Bregman, L.M.: A relaxation method of finding a common point of convex sets and its application to the solution of problems in convex programming. Ž. Vyčisl. Mat i Mat. Fiz. 7, 620–631 (1967)
Brezis, H.: How to recognize constant functions. A connection with Sobolev spaces. Uspekhi Mat. Nauk 57(4(346)), 59–74 (2002)
Carrillo, J.A., Jüngel, A., Markowich, P.A., Toscani, G., Unterreiter, A.: Entropy dissipation methods for degenerate parabolic problems and generalized Sobolev inequalities. Monatsh. Math. 133(1), 1–82 (2001)
Chen, Z.-Q., Fukushima, M.: Symmetric Markov Processes, Time Change, and Boundary Theory. London Mathematical Society Monographs Series, vol. 35. Princeton University Press, Princeton, NJ (2012)
Chung, K.L., Zhao, Z.X.: From Brownian motion to Schrödinger’s equation. Grundlehren der Mathematischen Wissenschaften, vol. 312. Springer, Berlin (1995)
Covi, G, Mönkkönen, K., Railo, J., Uhlmann, G.: The higher order fractional Calderón problem for linear local operators: uniqueness. Adv. Math., 399: Paper No. 108246, 29 (2022)
Davies, E.B.: Heat Kernels and Spectral Theory. Cambridge Tracts in Mathematics, vol. 92. Cambridge University Press, Cambridge (1990)
De Marco, G., Mariconda, C., Solimini, S.: An elementary proof of a characterization of constant functions. Adv. Nonlinear Stud. 8(3), 597–602 (2008)
C. Dellacherie and P.-A. Meyer: Probabilities and Potential. North Holland Mathematical Studies, 29. North-Holland Publishing Co., Amsterdam, New York (1978)
Dellacherie, C., Meyer, P.-A., Probabilities and Potential. B, volume 72 of North-Holland Mathematics Studies. North-Holland Publishing Co., Amsterdam. Theory of Martingales. Translated from the French by J. P, Wilson (1982)
Di Nezza, E., Palatucci, G., Valdinoci, E.: Hitchhiker’s guide to the fractional Sobolev spaces. Bull. Sci. Math. 136(5), 521–573 (2012)
DiBenedetto, E., Trudinger, N.S.: Harnack inequalities for quasiminima of variational integrals. Ann. Inst. H. Poincaré Anal. Non Linéaire 1(4), 295–308 (1984)
Dipierro, S., Ros-Oton, X., Valdinoci, E.: Nonlocal problems with Neumann boundary conditions. Rev. Mat. Iberoam. 33(2), 377–416 (2017)
Douglas, J.: Solution of the problem of Plateau. Trans. Am. Math. Soc. 33(1), 263–321 (1931)
Dyda, B., Kassmann, M.: Function spaces and extension results for nonlocal Dirichlet problems. J. Funct. Anal. 277(11), 108134 (2019)
Dynkin, E.B.: Markov processes. Vols. I, II, volume 122 of Die Grundlehren der Mathematischen Wissenschaften, Bände 121. Academic Press Inc., New York, Springer, Berlin (1965)
Farkas, W., Jacob, N., Schilling, R.L.: Feller semigroups, \(L^p\)-sub-Markovian semigroups, and applications to pseudo-differential operators with negative definite symbols. Forum Math. 13(1), 51–90 (2001)
Felsinger, M., Kassmann, M., Voigt, P.: The Dirichlet problem for nonlocal operators. Math. Z. 279(3–4), 779–809 (2015)
Foghem, G., Kassmann, M.: A general framework for nonlocal Neumann problems. arXiv e-prints. arXiv:2204.06793 (2022)
Frigyik, B.A., Srivastava, S., Gupta, M.R.: Functional Bregman divergence and Bayesian estimation of distributions. IEEE Trans. Inform. Theory 54(11), 5130–5139 (2008)
Garofalo, N.: Fractional thoughts. In: New Developments in the Analysis of Nonlocal Operators, volume 723 of Contemp. Math., pp. 1–135. American Mathematical Society, Providence, RI (2019)
Ghosh, T., Rüland, A., Salo, M., Uhlmann, G.: Uniqueness and reconstruction for the fractional Calderón problem with a single measurement. J. Funct. Anal., 279(1):108505, 42 (2020)
Ghosh, T., Salo, M., Uhlmann, G.: The Calderón problem for the fractional Schrödinger equation. Anal. PDE 13(2), 455–475 (2020)
Giaquinta, M., Giusti, E.: On the regularity of the minima of variational integrals. Acta Math. 148, 31–46 (1982)
Giusti, E.: Direct Methods in the Calculus of Variations. World Scientific Publishing Co. Inc, River Edge, NJ (2003)
Grafakos, L.: Classical Fourier Analysis, Volume 249 of Graduate Texts in Mathematics, 2nd edn. Springer, New York (2008)
Grzywny, T., Kassmann, M., Leżaj, Ł.: Remarks on the Nonlocal Dirichlet Problem. Potential Anal. Article in Press (2020)
Grzywny, T., Kwaśnicki, M.: Potential kernels, probabilities of hitting a ball, harmonic functions and the boundary Harnack inequality for unimodal Lévy processes. Stochastic Process. Appl. 128(1), 1–38 (2018)
Jacob, N.: Pseudo Differential Operators and Markov Processes. Fourier Analysis and Semigroups, vol. I. Imperial College Press, London (2001)
Jacob, N., Schilling, R.L.: Some Dirichlet spaces obtained by subordinate reflected diffusions. Rev. Mat. Iberoamericana 15(1), 59–91 (1999)
Kallenberg, O.: Foundations of Modern Probability. Probability and its Applications, 2nd edn. Springer, New York (2002)
Kulczycki, T., Ryznar, M.: Gradient estimates of harmonic functions and transition densities for Lévy processes. Trans. Am. Math. Soc. 368(1), 281–318 (2016)
Kwaśnicki, M.: Ten equivalent definitions of the fractional Laplace operator. Fract. Calc. Appl. Anal. 20(1), 7–51 (2017)
Langer, M., Maz’ya, V.: On \(L^p\)-contractivity of semigroups generated by linear partial differential operators. J. Funct. Anal. 164(1), 73–109 (1999)
Liskevich, V.A., Semenov, Y.A.: Some problems on Markov semigroups. In: Schrödinger operators, Markov Semigroups, Wavelet Analysis, Operator Algebras, vol. 11 of Math. Top., pp. 163–217. Akademie Verlag, Berlin (1996)
Nielsen, F., Nock, R.: Sided and symmetrized Bregman centroids. IEEE Trans. Inform. Theory 55(6), 2882–2904 (2009)
Pazy, A.: Semigroups of linear operators and applications to partial differential equations, volume 44 of Applied Mathematical Sciences, vol. 44. Springer, New York (1983)
Pietruska-Pałuba, K.: Heat kernels on metric spaces and a characterisation of constant functions. Manuscripta Math. 115(3), 389–399 (2004)
Protter, P.E.: Stochastic Integration and Differential Equations, Volume 21 of Stochastic Modelling and Applied Probability. Springer, Berlin (2005). Second edition. Version 2.1, Corrected third printing
Pruitt, W.E.: The growth of random walks and Lévy processes. Ann. Probab. 9(6), 948–956 (1981)
Ros-Oton, X.: Nonlocal elliptic equations in bounded domains: a survey. Publ. Mat. 60(1), 3–26 (2016)
Rutkowski, A.: The Dirichlet problem for nonlocal Lévy-type operators. Publ. Mat. 62(1), 213–251 (2018)
Sato, K.: Lévy Processes and Infinitely Divisible Distributions, Volume��68 of Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge. Translated from the 1990 Japanese original, Revised by the author (1999)
Servadei, R., Valdinoci, E.: Mountain pass solutions for non-local elliptic operators. J. Math. Anal. Appl. 389(2), 887–898 (2012)
Servadei, R., Valdinoci, E.: Variational methods for non-local operators of elliptic type. Discrete Contin. Dyn. Syst. 33(5), 2105–2137 (2013)
Sobol, Z., Vogt, H.: On the \(L_p\)-theory of \(C_0\)-semigroups associated with second-order elliptic operators. I. J. Funct. Anal. 193(1), 24–54 (2002)
Sprung, B.: Upper and lower bounds for the Bregman divergence. J. Inequal. Appl., 12: paper no. 4 (2019)
Vondraček, Z.: A probabilistic approach to a non-local quadratic form and its connection to the Neumann boundary condition problem. Math. Nachr. 294(1), 177–194 (2021)
Wang, F.-Y.: \(\Phi \)-entropy inequality and application for SDEs with jumps. J. Math. Anal. Appl. 418(2), 861–873 (2014)
Zeidler, E.: Nonlinear functional analysis and its applications. II, B. Nonlinear monotone operators. Translated from the German by the author and Leo F, Boron Springer-Verlag, New York (1990)
Ziemer, W.P.: Boundary regularity for quasiminima. Arch. Rational Mech. Anal. 92(4), 371–382 (1986)
Acknowledgements
We thank Tomasz Adamowicz, Włodzimierz Ba̧k, Artur Bogdan, Bartłomiej Dyda, Agnieszka Kałamajska, Moritz Kassmann, Mateusz Kwaśnicki, René Schilling and Enrico Valdinoci for discussions, comments or suggestions.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Laszlo Szekelyhidi.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Krzysztof Bogdan and Katarzyna Pietruska-Pałuba were supported by the NCN Grant 2018/31/B/ST1/03818. Tomasz Grzywny was supported by the NCN Grant 2017/27/B/ST1/01339. Artur Rutkowski was supported by the NCN Grant 2015/18/E/ST1/00239.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Bogdan, K., Grzywny, T., Pietruska-Pałuba, K. et al. Nonlinear nonlocal Douglas identity. Calc. Var. 62, 151 (2023). https://doi.org/10.1007/s00526-023-02458-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s00526-023-02458-x