Let s∈(0,1), 1<p<Ns and Ω⊂RN be an open bounded set. In this work we study the existence of solutions to problems (E±) Lu±g(u)=μ and u=0 a.e. in RN∖Ω, where g∈C(R) is a nondecreasing function, μ is a bounded Radon measure on Ω and L is an integro-differential operator with order of differentiability s∈(0,1) and summability p∈(1,Ns). More precisely, L is a fractional p-Laplace type operator. We establish sufficient conditions for the solvability of problems (E±). In the particular case g(t)=|t|κ−1t; κ>p−1, these conditions are expressed in terms of Bessel capacities.
Citation: Konstantinos T. Gkikas. Nonlinear nonlocal equations involving subcritical or power nonlinearities and measure data[J]. Mathematics in Engineering, 2024, 6(1): 45-80. doi: 10.3934/mine.2024003
[1] | María Ángeles García-Ferrero, Angkana Rüland . Strong unique continuation for the higher order fractional Laplacian. Mathematics in Engineering, 2019, 1(4): 715-774. doi: 10.3934/mine.2019.4.715 |
[2] | Federico Cluni, Vittorio Gusella, Dimitri Mugnai, Edoardo Proietti Lippi, Patrizia Pucci . A mixed operator approach to peridynamics. Mathematics in Engineering, 2023, 5(5): 1-22. doi: 10.3934/mine.2023082 |
[3] | Ko-Shin Chen, Cyrill Muratov, Xiaodong Yan . Layered solutions for a nonlocal Ginzburg-Landau model with periodic modulation. Mathematics in Engineering, 2023, 5(5): 1-52. doi: 10.3934/mine.2023090 |
[4] | Patrizia Pucci, Letizia Temperini . On the concentration–compactness principle for Folland–Stein spaces and for fractional horizontal Sobolev spaces. Mathematics in Engineering, 2023, 5(1): 1-21. doi: 10.3934/mine.2023007 |
[5] | Boumediene Abdellaoui, Ireneo Peral, Ana Primo . A note on the Fujita exponent in fractional heat equation involving the Hardy potential. Mathematics in Engineering, 2020, 2(4): 639-656. doi: 10.3934/mine.2020029 |
[6] | Catharine W. K. Lo, José Francisco Rodrigues . On the obstacle problem in fractional generalised Orlicz spaces. Mathematics in Engineering, 2024, 6(5): 676-704. doi: 10.3934/mine.2024026 |
[7] | José Antonio Vélez-Pérez, Panayotis Panayotaros . Wannier functions and discrete NLS equations for nematicons. Mathematics in Engineering, 2019, 1(2): 309-326. doi: 10.3934/mine.2019.2.309 |
[8] | Francesca Colasuonno, Fausto Ferrari, Paola Gervasio, Alfio Quarteroni . Some evaluations of the fractional p-Laplace operator on radial functions. Mathematics in Engineering, 2023, 5(1): 1-23. doi: 10.3934/mine.2023015 |
[9] | Nguyen Cong Phuc, Igor E. Verbitsky . Uniqueness of entire solutions to quasilinear equations of p-Laplace type. Mathematics in Engineering, 2023, 5(3): 1-33. doi: 10.3934/mine.2023068 |
[10] | Masashi Misawa, Kenta Nakamura, Yoshihiko Yamaura . A volume constraint problem for the nonlocal doubly nonlinear parabolic equation. Mathematics in Engineering, 2023, 5(6): 1-26. doi: 10.3934/mine.2023098 |
Let s∈(0,1), 1<p<Ns and Ω⊂RN be an open bounded set. In this work we study the existence of solutions to problems (E±) Lu±g(u)=μ and u=0 a.e. in RN∖Ω, where g∈C(R) is a nondecreasing function, μ is a bounded Radon measure on Ω and L is an integro-differential operator with order of differentiability s∈(0,1) and summability p∈(1,Ns). More precisely, L is a fractional p-Laplace type operator. We establish sufficient conditions for the solvability of problems (E±). In the particular case g(t)=|t|κ−1t; κ>p−1, these conditions are expressed in terms of Bessel capacities.
Let Ω⊂RN be an open bounded domain, s∈(0,1) and 1<p<Ns. In this article we are concerned with the existence of very weak solutions to the quasilinear nonlocal problems
{Lu±g(u)=μ,inΩ,u=0,inRN∖Ω, | ( |
where μ is a bounded Radon measure on Ω and g∈C(R) is a nondecreasing function such that g(0)=0. Here, the nonlocal operator L is defined by
Lu(x):=P.V.∫RN|u(x)−u(y)|p−2(u(x)−u(y))K(x,y)dy,∀x∈Ω, |
where the symbol P.V. stands for the principle value integral and K:RN×RN→R is a measurable and symmetric (i.e., K(x,y)=K(y,x)) function. Note that if K(x,y)≡|x−y|−N−sp then L coincides with the standard fractional p-Laplace operator (−Δ)sp.
Throughout this work, we assume that there exists a positive constant ΛK≥1 such that the following ellipticity condition holds
Λ−1K|x−y|−N−sp≤K(x,y)≤ΛK|x−y|−N−sp,∀(x,y)∈RN×RNandx≠y. |
In addition, we denote by Mb(Ω) the space of Radon measures on RN such that μ(RN∖Ω)=0, as well as by M+b(Ω) its positive cone.
Let
CN,s:=22sπ−N2sΓ(N+2s2)Γ(1−s)>0. |
For p=2 and K(x,y)=CN,s|x−y|−N−2s, operator L reduces to the well-known fractional Laplace operator (−Δ)s and the problem P+ becomes
{(−Δ)su+g(u)=μ,inΩu=0,inRN∖Ω. | (1.1) |
When g satisfies the subcritical integral condition
∫∞1(g(s)−g(−s))s−NN−2s−1ds<∞, |
Chen and Véron [9] showed that problem (1.1) admits a unique very weak solution for any μ∈Mb(Ω). In addition they showed that problem (1.1) with g(u)=|u|κ−1u(κ>1) possesses a very weak solution if and only if μ is absolutely continuous with respect to Bessel capacity CL2s,κ′, i.e., μ vanishes on compact set E of Ω satisfying Cap2s,κ′(E)=0 (see (3.21) for the definition of the Bessel capacities). Their approach is based on the properties of the Green Kernel associated with fractional Laplace operator (−Δ)s in Ω.
In the local theory and more precisely when Lu=−Δpu=−div(|∇u|p−2∇u), related problems have been studied in [4,5,6,15,30,31,32]. In particular, in the power case, i.e.,
{−Δpu+|u|κ−1u=μ,inΩ,u=0,on∂Ω, | (1.2) |
Bidaut-Véron, Nguyen and Véron [5] established that if μ∈Mb(Ω) is absolutely continuous with respect to the Bessel capacity Capp,κκ−p+1, then there exists a renormalized solution to problem (1.2) with κ>p−1. A main ingredient in the proof of this result is the pointwise estimates for p-superharmonic functions in Ω. These pointwise estimates are expressed in terms of the truncated Wolff potentials WR1,p[μ] (see, e.g., [17,19,20,31]). We recall here that the truncated Wolff potential is given by
WRα,p[μ](x):=∫R0(|μ|(Br(x))rN−αp)1p−1drr, | (1.3) |
for any R>0 and α∈(0,N) such that p∈(1,Nα). Conversely, Bidaut-Véron [4] showed that if problem (1.2) with κ>p−1 admits a renormalized solution, then μ is absolutely continuous with respect to the Bessel capacity Capp,κκ−p+1+ε, for any ε>0.
Phuc and Verbitsky [31] showed that if τ∈M+b(Ω) has compact support in Ω, then the problem
{−Δpu−|u|κ=ρτ,inΩ,u=0,on∂Ω, | (1.4) |
admits a nonnegative renormalized solution for some ρ>0, if and only if, there exists a positive constant C such that
τ(K)≤CCapp,κκ−p+1(K), | (1.5) |
for any compact K⊂Ω. Moreover, they showed that (1.5) is equivalent to
W2diam(Ω)1,p[(W2diam(Ω)1,p[τ])κ]≤CW2diam(Ω)1,p[τ],a.e. inΩ, |
for some positive constant C>0.
Recently, a great attention has been drawn to the study of the fractional p-Laplacian or more general nonlocal operators (see for example [2,11,12,18,21,22,23,24,25,26,27,28,29]). More precisely, Kuusi, Mingione and Sire [26] dealt with the problem
{LΦu=μ,inΩ,u=g,inRN∖Ω, | (1.6) |
where g∈Ws,p(RN), LΦ is a nonlocal operator defined by
<LΦu,ζ>:=∫RN∫RNΦ(u(x)−u(y))(ζ(x)−ζ(y))K(x,y)dydx,∀ζ∈C∞0(Ω). |
Here Φ:R→R is a continuous function such that Φ(0)=0 and
Λ−1Φ|t|p≤Φ(t)t≤ΛΦ|t|p. |
When 2−sN<p, they show the existence of a very weak solution to (1.6), which they called SOLA (Solutions obtained as limits of approximations). They also showed local pointwise estimates for SOLA to (1.6) in terms of the truncated Wolff Potential WRs,p[μ]. In the particular case Φ(t)=|t|p−2t and g=0, the existence of very weak solutions was established in [2] for any 1<p<Ns.
The objective of this work is to determine the subcritical integral conditions on g, which ensure the existence of very weak solutions to problems (P±). In addition, in the power case, i.e., g(u)=|u|κ−1u; κ>p−1, we aim to find sufficient conditions, expressed in terms of Bessel capacities like above, for the solvability of (P±).
Let us mention here that our work is inspired by the article [5] for problem (P+) and by the articles [30,31] for problem (P−) with g(u)=|u|κ−1u;κ>p−1. However, due to the presence of the nonlocal operator, new essential difficulties arise which complicate drastically the study of problems (P±).
In order to state our main results, we need to introduce the notion of the very weak solutions.
Definition 1.1. Let s∈(0,1), 1<p<Ns, ˜g∈C(R), Ω⊂RN be an open bounded domain and μ∈M(Ω). We will say that u:RN→R is a very weak solution to the problem
{Lu+˜g(u)=μ,inΩ,u=0,inRN∖Ω, | (1.7) |
if ˜g(u)∈L1loc(Ω) and if the following conditions are valid:
(ⅰ) u=0 a.e. in RN∖Ω and u∈Wh,q(RN) for any 0<h<s and for any 0<q<N(p−1)N−s.
(ⅱ) Tk(u):=max(−k,min(u,k))∈Ws,p0(Ω) for any k>0.
(ⅲ)
∫RN∫RN|u(x)−u(y)|p−2(u(x)−u(y))(ϕ(x)−ϕ(y))K(x,y)dxdy+∫Ω˜g(u)ϕdx=∫Ωϕdμ |
for any ϕ∈C∞0(Ω).
We note here that if 2−sN<p<Ns, then the very weak solution u belongs to the fractional Sobolev space Wh,q(RN) for any q∈(1,N(p−1)N−s). If p≤2−sN, the space Wh,q(RN) in the above definition is no longer a fractional Sobolev space, however it is defined in the same way (see (2.1)).
In Section 2, we discuss the existence and main properties of the very weak solutions of problem (1.7) with ˜g≡0. Particularly, in the spirit of [26], we show the existence of a SOLA u satisfying statements (ⅰ)–(ⅲ) of the above definition (see Proposition 2.8). The approximation sequence consists of solutions of (1.7) with ˜g≡0 and smooth data. In addition, we prove that these solutions satisfy a priori estimates (2.8) and (2.11). As a result, we establish that the very weak solution satisfies (2.11) and
‖|u|p−1‖∗LNN−spw(RN)≤C(N,p,s,ΛK)∫Ω|μ|dx, | (1.8) |
where ‖⋅‖∗LNN−spw(RN) has been defined in (2.4) and is related to the Marcinkiewicz spaces. Finally, when μ∈M+b(Ω), we construct this solution (see Propositions 2.9 and 2.10) such that u≥0 and
C−1(N,p,s,ΛK)Wd(x)8s,p[μ](x)≤u(x)≤C(N,p,s,ΛK)W2diam(Ω)s,p[μ](x),a.e. inΩ, |
where d(x)=dist(x,∂Ω). The lower estimate in the above display can be obtained as a consequence of [26, estimate (1.25)]. The upper estimate in the above display is an application of [21, Theorem 5.3] and (1.8).
Using the above properties of the very weak solutions and the fact that if u,g satisfies (1.8) and (1.9) respectively then g(u)∈L1(Ω), we obtain the following result.
Theorem 1.2. Let s∈(0,1), 1<p<Ns, μ∈Mb(Ω). We assume that g∈C(R) is a nondecreasing function satisfying g(0)=0 and
∫∞1(g(s)−g(−s))s−N(p−1)N−sp−1ds<∞. | (1.9) |
Then there exist a very weak solution u to problem (P+) satisfying (1.8) and
−C(N,p,s,ΛK)W2diam(Ω)s,p[μ−]≤u≤C(N,p,s,ΛK)W2diam(Ω)s,p[μ+],a.e. inΩ. | (1.10) |
In addition, for any q∈(0,N(p−1)N−s) and h∈(0,s), there exists a positive constant c=c(N,p,s,ΛK,q,h,|Ω|) such that
(∫Ω|g(u)|dx)1p−1+(∫RN∫RN|u(x)−u(y)|q|x−y|N+hqdxdy)1q≤c(|μ|(Ω))1p−1. | (1.11) |
We note here that the integral conditions (1.9) and (1) coincide for p=2. In addition, in the corresponding local case, the integral condition (1.9) with s=1 ensures the existence of the associated renormalized solutions (see [32, Theorem 5.1.2 and (5.1.40)]).
Let us consider problem (P+) with a power absorption, i.e.,
{Lu+|u|κ−1u=μ,inΩ,u=0,inRN∖Ω. | (1.12) |
We first notice that the function g(t)=|t|κ−1t with k>0 satisfies (1.9) if and only if 0<κ<N(p−1)N−sp, hence problem (1.12) admits a very weak solution in this case. In the supercritical case κ≥N(p−1)N−sp, the sufficient condition for the solvability of problem (1.12) is expressed in terms of the Bessel capacity Capsp,κκ−p+1 as follows.
Theorem 1.3. Let s∈(0,1), 1<p<Ns, κ>p−1 and μ∈Mb(Ω). In addition we assume that μ is absolutely continuous with respect to the Bessel capacity Capsp,κκ−p+1. Then there exists a very weak solution u to problem (1.12) such that
−CW2diam(Ω)s,p[μ−]≤u≤CW2diam(Ω)s,p[μ+],a.e. inΩ. | (1.13) |
In addition, for any q∈(0,N(p−1)N−s) and h∈(0,s), there exists a positive constant c=c(N,p,s,ΛK,q,h,|Ω|) such that
(∫Ω|u|κdx)1p−1+(∫RN∫RN|u(x)−u(y)|q|x−y|N+hqdxdy)1q≤c(|μ|(Ω))1p−1. | (1.14) |
In view of the discussion on the existence of solutions to problem (1.4), we expect that the existence phenomenon occurs for (P−) only for measures μ∈Mb(Ω) with small enough total mass. Indeed, using the Schauder fixed point theorem and sharp weak Lebesgue estimates, we prove the following existence result for any μ∈Mb(Ω) with small enough total mass.
Theorem 1.4. Let s∈(0,1), 1<p<Ns and τ∈Mb(Ω) be such that |τ|(Ω)≤1. Assume that g∈C(R) is a nondecreasing function satisfying (1.9) and
|g(s)|≤a|s|dfor somea>0,d>1and for any|s|≤1. | (1.15) |
Then there exists a positive constant ρ0 depending on N,|Ω|,Λg,ΛK,a,s,p,d,|Ω| such that for every ρ∈(0,ρ0) the following problem
{Lv=g(v)+ρτ,inΩ,v=0,inRN∖Ω, | (1.16) |
admits a very weak solution v satisfying
‖|v|p−1‖∗LNN−spw(RN)≤t0. | (1.17) |
Here, t0>0 depends on N,|Ω|,Λg,ΛK,a,s,p,d,ρ0. In addition, for any q∈(0,N(p−1)N−s) and h∈(0,s), there exists a positive constant c depending only on N,p,s,Λg,ΛK,q,h,|Ω|,a,d,ρ0 and t0, such that
(∫RN∫RN|v(x)−v(y)|q|x−y|N+hqdxdy)1q≤c(1+ρ|τ|(Ω))1p−1. | (1.18) |
In the linear case, i.e., p=2, problem (P−) with L=(−Δ)s was thoroughly studied in [7]. More precisely, the authors in [7] showed that the same existence result occurs provided g satisfies (1) and (1.15).
Problem (P−) with g(t)=|t|κ−1t and μ∈M+b(Ω) becomes
{Lv=|v|κ−1v+ρτ,inΩ,v=0,inRN∖Ω. | (1.19) |
When p=2, problem (P−) with L=(−Δ)s and τ=δ0 was studied in [8]. Here δ0 denotes the dirac measure concentrated at a point x0∈Ω. In particular, the authors in [8] established that if κ≥NN−2s and u is a nonnegative solution of (1.19) then ρ=0. Concerning problem (1.19), conditions (1.9) and (1.15) are satisfied if κ belongs to the subcritical range, that is when p−1<κ<N(p−1)N−sp. In general, a sufficient condition for the solvability of (1.19) is the following.
Proposition 1.5. Let s∈(0,1), 1<p<Ns, κ>p−1 and τ∈M+b(Ω) be such that
W2diam(Ω)s,p[(W2diam(Ω)s,p[τ])κ]≤MW2diam(Ω)s,p[τ],a.e.inΩ, | (1.20) |
for some positive constant M. Then problem (1.19) admits a nonnegative very weak solution u for some ρ>0. Furthermore, there holds
M−1Wd(x)8s,p[μ](x)≤u(x)≤MW2diam(Ω)s,p[ρτ](x),for a.e.x∈Ω, | (1.21) |
where dμ=uκdx+ρdτ and the positive constant M depends only on C,N,p,q,ΛK.
Finally, inspired from Phuc and Verbitsky's ideas in [30,31], we establish the following existence result in the whole range κ>p−1.
Theorem 1.6. Let s∈(0,1), 1<p<Ns, κ>p−1 and τ∈M+b(Ω) with compact support in Ω. Then the following statements are equivalent.
(i) Problem (1.19) admits a nonnegative very weak solution uρ for some ρ>0 such that
C−11Wd(x)8s,p[μ](x)≤uρ(x)≤C1W2diam(Ω)s,p[ρτ](x),for a.e.x∈Ω, | (1.22) |
where dμ=uκdx+ρdτ and for some constant C1>0. (ii) There exists a positive constant C2 such that
τ(E)≤C2Capsp,κκ−p+1(E) | (1.23) |
for any Borel set E⊂RN.
(iii) There exists a positive constant C3 such that
∫B(W2diam(Ω)s,p[τ⌊B])κdx≤C3τ(B) | (1.24) |
for any ball B⊂RN, where dτ⌊B=χBdτ.
(iv) There exists a positive constant C4 such that
W2diam(Ω)s,p[(W2diam(Ω)s,p[τ])κ]≤C4W2diam(Ω)s,p[τ],a.e.inΩ. |
We note here that if p−1<q<N(p−1)N−sp then spqq−p+1>N, this implies that Capsp,qq−p+1({x})>0 for any x∈RN (see [1, Section 2.6]). Hence, the statement (ⅱ) in the above theorem is always satisfied in the subcritical range.
Section 2 is devoted to the study of the very weak solutions to problem (1.7) with ˜g≡0. In Section 3, we discuss problem (P+) as well as Theorems 1.2 and 1.3 are proved in Subsections 3.2 and 3.3 respectively. In section 4, we deal with problem (P−). More precisely, we prove Theorem 1.4 in Subsection 4.1 and demonstrate Proposition 1.5 and Theorem 1.6 in Subsection 4.2.
We start with the definition of the fractional spaces, which will be used frequently in this work. For any s∈(0,1) and q>0, we denote by Ws,q(RN) the fractional space
Ws,q(RN):={∫RN∫RN|u(x)−u(y)|q|x−y|N+sqdxdy+∫RN|u|qdx<∞}, | (2.1) |
endowed with the quasinorm
‖u‖Ws,q(RN):=(∫RN∫RN|u(x)−u(y)|q|x−y|N+sqdxdy)1q+(∫RN|u|qdx)1q. |
When q≥1, Ws,q(RN) is a Banach space and is called fractional Sobolev space. Finally, for any p>1, we denote by Ws,p0(Ω) the closure of C∞0(Ω) in the norm ‖⋅‖Ws,p(RN) and by (Ws,p0(Ω))∗ its dual space.
In this subsection, we introduce the notion of the weak solution of the following problem
{Lu=μ,inΩ,u=0,inRN∖Ω, | (2.2) |
where μ∈(Ws,p0(Ω))∗. In addition, when μ∈Lp′(Ω), we establish a priori estimates, which will be used in the construction of the very weak solutions of the above problem with measure data.
Definition 2.1. Let s∈(0,1), p>1, and μ∈(Ws,p0(Ω))∗. We will say that u∈Ws,p0(Ω) is a weak solution of (2.2), if it satisfies
∫RN∫RN|u(x)−u(y)|p−2(u(x)−u(y))(ϕ(x)−ϕ(y))K(x,y)dxdy=<μ,ϕ>,∀ϕ∈Ws,p0(Ω). |
Let us now give the definition of weak supersolutions of L in Ω.
Definition 2.2. Let s∈(0,1) and p>1. We will say that u∈Ws,p(RN) is a weak supersolution (resp. subsolution) of L in Ω, if and only if satisfies
∫RN∫RN|u(x)−u(y)|p−2(u(x)−u(y))(ϕ(x)−ϕ(y))K(x,y)dxdy≥0(resp.≤0) |
for any nonnegative ϕ∈Ws,p0(Ω).
Next we state the comparison principle.
Proposition 2.3. ([23, Lemma 6]). Let u∈Ws,p(RN) be a weak supersolution of L in Ω as well as let v∈Ws,p(RN) be a weak subsolution of L in Ω such that (v−u)+∈Ws,p0(Ω). Then, u≥v a.e. in RN.
In view of the proof [11, Theorem 2.3], we may obtain the following existence result.
Proposition 2.4. For any μ∈(Ws,p0(Ω))∗ there exists a unique weak solution of (2.2).
In order to state the first a priori estimate for the weak solution of (2.2), we need to give the definition and the main properties of Marcinkiewicz spaces. Let D⊂RN be a domain. Denote Lpw(D), 1≤p<∞, the weak Lp space (or Marcinkiewicz space) defined as follows. A measurable function f in D belongs to this space if there exists a constant c such that
λf(a):=|{x∈D:|f(x)|>a}|≤ca−p,∀a>0. | (2.3) |
The function λf is called the distribution function of f. For p≥1, denote
Lpw(D)={f Borel measurable:supa>0apλf(a)<∞}, |
‖f‖∗Lpw(D)=(supa>0apλf(a))1p. | (2.4) |
The ‖.‖∗Lpw(D) is not a norm, but for p>1, it is equivalent to the norm
‖f‖Lpw(D)=sup{∫ω|f|dx|ω|1/p′:ω⊂D,ω measurable,0<|ω|<∞}. | (2.5) |
More precisely,
‖f‖∗Lpw(D)≤‖f‖Lpw(D)≤pp−1‖f‖∗Lpw(D). | (2.6) |
Notice that,
Lpw(D)⊂Lr(D),∀r∈[1,p). |
From (2.4) and (2.6), one can derive the following estimate which is useful in the sequel.
∫{|u|≥s}dx≤s−p‖u‖pLpw(D). | (2.7) |
Proposition 2.5. Let 1<p<Ns, μ∈Lp′(Ω) and u∈Ws,p0(Ω) be the unique weak solution of (2.2). Then there exists a positive constant C=C(p,s,N,ΛK) such that
‖|u|p−1‖∗LNN−spw(RN)≤C∫Ω|μ|dx. | (2.8) |
Proof.
Let k>0. Taking Tk(u) as test function and using the fact that
|u(x)−u(y)|p−2(u(x)−u(y))(Tk(u)(x)−Tk(u)(y))≥|Tk(u)(x)−Tk(u)(y)|p,∀x,y∈RN, |
we obtain
∫RN∫RN|Tk(u)(x)−Tk(u)(y)|p|x−y|N+spdxdy≤ΛKk∫Ω|μ|dx. | (2.9) |
Now, by the above inequality and the fractional Sobolev inequality we have
|{|u(x)|≥k}|=|{|Tk(u)(x)|≥k}|≤k−NpN−sp∫RN|Tk(u)(x)|NpN−spdx≤Ck−N(p−1)N−sp(∫Ω|μ|dx)NN−sp, |
which implies the desired result.
Proposition 2.6. Let μ∈Lp′(Ω) and u∈Ws,p0(RN) be the unique weak solution of (2.2). Then there exists a positive constant C=C(p,s,N,ΛK) such that
∫RN∫RN|u(x)−u(y)|p(d+|u(x)|+|u(y)|)ξdxdy|x−y|N+sp≤Cd1−ξ(ξ−1)∫Ω|μ|dx | (2.10) |
for any ξ>1 and d>0.
Proof. The proof is very similar to that of [26, Lemma 3.1] (see also [25, Lemma 8.4.1]). For the sake of convenience we give it below.
Set ϕ±:=±(d1−ξ−(d+u±)1−ξ). Using ϕ± as test function we obtain
∫RN∫RN|u(x)−u(y)|p−2(u(x)−u(y))(ϕ±(x)−ϕ±(y))K(x,y)dxdy=∫Ωϕ±μdx. |
Now,
(ϕ±(x)−ϕ±(y))=±(ξ−1)(u±(x)−u±(y))∫10(d+tu±(y)+(1−t)u±(x))−ξdt, |
which implies
|u(x)−u(y)|p−2(u(x)−u(y))(ϕ±(x)−ϕ±(y))K(x,y)≥(ξ−1)|u(x)−u(y)|p−2(u±(x)−u±(y))2(d+|u(y)|+|u(x)|)−ξ. |
Combining all above we can easily reach the desired result.
We conclude this subsection by the following a priori estimate for the weak solutions of (2.2) in the whole range p>1.
Proposition 2.7. Let ¯q=min{N(p−1)N−s,p} μ∈Lp′(Ω) and u∈Ws,p0(RN) be the unique weak solution of (2.2). For any q∈(0,¯q) and h∈(0,s), there exists a positive constant c depending only on N,s,p,ΛK,q and |Ω| such that
(∫RN∫RN|u(x)−u(y)|q|x−y|N+hqdxdy)1q≤c(∫Ω|μ|dx)1p−1. | (2.11) |
Proof. The proof is an adaptation of the argument in [26, Lemma 3.2]. Let R=diam(Ω) and x0∈Ω. First, we note that
∫RN∫RN|u(x)−u(y)|q|x−y|N+hqdxdy=∫B2R(x0)∫B2R(x0)|u(x)−u(y)|q|x−y|N+hqdxdy+∫RN∖B2R(x0)∫RN∖B2R(x0)|u(x)−u(y)|q|x−y|N+hqdxdy+2∫RN∖B2R(x0)∫B2R(x0)|u(x)−u(y)|q|x−y|N+hqdxdy. |
Taking into account that u=0 a.e. in RN∖Ω, we can easily prove that
∫RN∖B2R(x0)∫RN∖B2R(x0)|u(x)−u(y)|q|x−y|N+hqdxdy=0 |
and
∫RN∖B2R(x0)∫B2R(x0)|u(x)−u(y)|q|x−y|N+hqdxdy=∫RN∖B2R(x0)∫BR(x0)|u(x)|q|x−y|N+hqdxdy≈∫BR(x0)|u(x)|qdx∫RN∖B2R(x0)1(1+|y−x0|)N+hqdy≈∫BR(x0)|u(x)|qdx. |
Here, we have also used the fact that |x−y|≈1+|y−x0| for any (x,y)∈BR(x0)×(RN∖B2R(x0)), where the implicit constants in the last estimate depend only on R. Similarly, we have
∫B2R(x0)∫B2R(x0)|u(x)−u(y)|q|x−y|N+hqdxdy=∫B3R2(x0)∫B3R2(x0)|u(x)−u(y)|q|x−y|N+hqdxdy+∫B2R(x0)∖B3R2(x0)∫B2R(x0)∖B3R2(x0)|u(x)−u(y)|q|x−y|N+hqdxdy+2∫B2R(x0)∖B3R2(x0)∫B3R2(x0)|u(x)−u(y)|q|x−y|N+hqdxdy≈∫B3R2(x0)∫B3R2(x0)|u(x)−u(y)|q|x−y|N+hqdxdy+∫BR(x0)|u|qdx. |
Combining all above, we have
∫RN∫RN|u(x)−u(y)|q|x−y|N+hqdxdy≈∫B2R(x0)∫B2R(x0)|u(x)−u(y)|q|x−y|N+hqdxdy+∫Ω|u|qdx≈∫B2R(x0)∫B2R(x0)|u(x)−u(y)|q|x−y|N+hqdxdy. | (2.12) |
Now, by Hölder inequality we obtain
∫B2R(x0)∫B2R(x0)|u(x)−u(y)|q|x−y|N+hqdxdy=∫B2R(x0)∫B2R(x0)(|u(x)−u(y)|p(d+|u(x)|+|u(y)|)ξ|x−y|ps(d+|u(x)|+|u(y)|)ξ|x−y|p(s−h))qpdxdy|x−y|N≤(∫B2R(x0)∫B2R(x0)|u(x)−u(y)|p(d+|u(x)|+|u(y)|)ξ|x−y|N+spdxdy)qp×(∫B2R(x0)∫B2R(x0)(d+|u(x)|+|u(y)|)ξqp−q|x−y|N−qp(s−h)p−qdxdy)p−qp. | (2.13) |
Setting
d=(∫Ω|u(y)|ξqp−qdx)p−qξq |
and combining (2.10) and (2.13), we conclude
∫B2R(x0)∫B2R(x0)|u(x)−u(y)|q|x−y|N+hqdxdy≤cdqp(∫Ω|μ|dx)qp. | (2.14) |
If p>2−sN, without loss of generality, we may assume that q>1. Therefore, we may apply the fractional Sobolev inequality to d as in the proof of [26, Lemma 3.2] to obtain the desired result.
If 1<p≤2−Ns, we have that 0<q≤1, therefore, we can not apply the fractional Sobolev inequality to d. To overcome this difficulty we use (2.8) instead of fractional Sobolev inequality. More precisely, let 1<p<Ns, then 0<q<N(p−1)N−s<p. Hence, we may choose ξ>1 such that 1<γ:=ξq(p−1)(p−q)<NN−sp. Thus, by (2.6) and (2.8), we deduce
(∫Ω|u|γ(p−1))1γ≤C(γ,N,p,s,|Ω|,ΛK)∫Ω|μ|dx, |
which in turn implies
d≤C(γ,N,p,s,|Ω|,ΛK)(∫Ω|μ|dx)1p−1. |
The desired result follows by (2.12), (2.14) and the above inequality.
In this subsection, we construct a very weak solution to problem (2.2) which possesses several important properties, such as it satisfies pointwise estimates in terms of Wolff's potential. These estimates play an important role in the study of problems (P±).
We start with the following existence result.
Proposition 2.8. Let 1<p<Ns and μ∈Mb(Ω). Then there exists a very weak solution to (2.2) satisfying
‖|u|p−1‖∗LNN−spw(RN)≤C1(N,p,s,ΛK)μ(Ω) | (2.15) |
and
∫RN∫RN|Tk(u)(x)−Tk(u)(y)|p|x−y|N+spdxdy≤kΛK|μ|(Ω),∀k>0. | (2.16) |
In addition, for any q∈(0,N(p−1)N−s) and h∈(0,s), there exist a positive constant C2=C2(N,p,s,ΛK,q,h,|Ω|) such that
(∫RN∫RN|u(x)−u(y)|q|x−y|N+hqdxdy)1q≤C2|μ|(Ω)1p−1. | (2.17) |
Proof. Let {ρn}n be a sequence of mollifiers and μn=ρn∗μ. Then μn∈C∞0(RN) and μn⇀μ weakly in RN. We denote by un the weak solution of (2.2) with μ=μn.
By (2.8), (2.9) and (2.11), there exist positive constants C1 and C2 such that
‖|un|p−1‖∗LNN−spw(RN)≤C1(N,p,s,ΛK)μ(Ω),∀n∈N, | (2.18) |
∫RN∫RN|Tk(un)(x)−Tk(un)(y)|p|x−y|N+spdxdy≤kΛKμ(Ω),∀k>0andn∈N, | (2.19) |
and
(∫RN∫RN|un(x)−un(y)|q|x−y|N+hqdxdy)1q≤C2(N,p,s,ΛK,q,h)μ(Ω)1p−1 | (2.20) |
for any n∈N, q∈(0,N(p−1)N−s) and h∈(0,s). In the spirit of the proof of [10, Theorem 3.4], we will show that the existence of a subsequence (still denoted by {un}) and a function u:RN→R satisfying the following properties:
(ⅰ) u∈Wh,q(RN) for any 0<q<N(p−1)N−s and 0<h<s.
(ⅱ) un→u a.e. in RN, u=0 a.e. in RN∖Ω and ‖u−un‖Wh,q(RN)→0 for any q∈(0,N(p−1)N−s) and h∈(0,s).
(ⅲ) Tk(u)∈Ws,p0(Ω) for any k>0.
Step 1. There exists a subsequence, still denoted by un, such that
limn→∞limm→∞|{x∈Ω:|un−um|>η}|=0,∀η>0. |
Let n,m∈N and η,ρ>0. Then
{|un−um|>η}⊂{|Tk(un)|>k}∪{|Tk(um)|>k}∪{|Tk(un)−Tk(um)|>η}. |
By (2.18) and (2.7), there exists k0>0 such that
|{|Tk(un)|>k}|+|{|Tk(um)|>k}|≤ρ2,∀k≥k0. | (2.21) |
By (2.19), the fractional Sobolev embedding theorem (see e.g., [13, Corollary 7.2]) and the fact that Ws,p0(Ω) is a reflexive Banach space, we may prove the existence of a subsequence Tk0(unj) of Tk0(un) such that Tk0(unj)→vk0 in Lp(RN) and a.e. in RN as well as Tk0(unj)⇀vk0 in Ws,p0(Ω). Hence,
|{|Tk0(unj)−Tk0(um˜j)|>η}|≤ρ2,∀j,˜j≥n0. | (2.22) |
The desired result follows by (2.21) and (2.22).
Step 2. Weak convergence of the truncates. Since un→u a.e. in RN, we have that Tk(un)→Tk(u) a.e. in RN. Furthermore, by (2.19) and the fractional Sobolev embedding theorem, we can find a subsequence {Tk(unj)}∞j=1 such that Tk(unj)→vk in Lp(RN) and Tk(unj)⇀vk in Ws,p0(Ω). Since vk=Tk(u) a.e. in RN, we have that Tk(u)∈Ws,p0(Ω). This implies that the limit does not depend on the subsequence. Hence, for the same subsequence un of the Step 1, we have that
Tk(un)⇀Tk(u)inWs,p0(Ω),∀k>0. |
Furthermore, by (2.18)–(2.20) and Fatou's lemma, we have that
‖|u|p−1‖∗LNN−spw(RN)≤C1(N,p,s,ΛK)μ(Ω), | (2.23) |
∫RN∫RN|Tk(u)(x)−Tk(u)(y)|p|x−y|N+spdxdy≤kΛKμ(Ω),∀k>0, | (2.24) |
and
(∫RN∫RN|u(x)−u(y)|q|x−y|N+hqdxdy)1q≤C2(N,p,s,ΛK,q,h)μ(Ω)1p−1 | (2.25) |
for any q∈(0,N(p−1)N−s) and h∈(0,s).
By (2.20), (2.25) and the fact that un→u a.e. in RN, We can easily show that ‖u−un‖Wh,q(RN)→0 for any q∈(0,N(p−1)N−s) and h∈(0,s). Let ϕ∈C∞0(Ω), q∈(p−1,N(p−1)N−s) and h∈(max(sp−1p−1,0),s). For any bounded Borel set E⊂RN, we have that
|∫E∫E|un(x)−un(y)|p−2(un(x)−un(y))(ϕ(x)−ϕ(y))K(x,y)dxdy|≤C(ϕ,ΛK)∫E∫E|un(x)−un(y)|p−1|x−y|N+hp−1+sp−h(p−1)−1≤C(ϕ,ΛK)(∫E∫E|un(x)−un(y)|q|x−y|N+hqdxdy)p−1q(∫E∫E|x−y|−N−q(sp−h(p−1)−1)q−p+1dxdy)q−p+1q. |
This, together with (ⅰ), (ⅱ) and the fact that q(sp−h(p−1)−1)q−p+1<0, implies that
∫Ωϕdμ=limn→∞∫Ωϕμndx=limn→∞∫RN∫RN|un(x)−un(y)|p−2(un(x)−un(y))(ϕ(x)−ϕ(y))K(x,y)dxdy=∫RN∫RN|u(x)−u(y)|p−2(un(x)−un(y))(ϕ(x)−ϕ(y))K(x,y)dxdy. |
The proof is complete.
In the next theorem, we establish a priori pointwise estimates for a certain nonnegative very weak solution of problem (2.2) with μ∈M+b(Ω).
Proposition 2.9. Let 1<p<Ns and μ∈M+b(Ω). Then there exist a nonnegative very weak solution u of (2.2) and a positive constant C depending only on N,s,p,ΛK such that
C−1Wd(x)8s,p[μ](x)≤u(x)≤C(essinfBd(x)4(x)u+Wd(x)2s,p[μ](x)+((d(x)4)sp∫RN∖Bd(x)4(x)u(y)p−1|x−y|N+spdy)1p−1) | (2.26) |
for a.e. x∈Ω.
Proof. Let u be the solution constructed in Proposition 2.8 and {un} be the sequence defined in Proposition 2.8 such that
(ⅰ) u∈Wh,q(RN) for any 0<q<N(p−1)N−s and 0<h<s.
(ⅱ) un→u a.e. in RN, u=0 a.e. in RN∖Ω and ‖u−un‖Wh,q(RN)→0 for any q∈(0,N(p−1)N−s) and h∈(0,s).
Since μn∈C∞0(RN) is a nonnegative function, by (2.3), we have that un≥0 a.e. in RN. Hence, by [23, Lemma 7], uk,n=min(un,k) is a nonnegative weak supersolution. By properties (ⅰ) and (ⅱ), we may show that uk=min(u,k) is a nonnegative weak supersolution. Hence, there exists a nonnegative Radon measure μk∈M+(Ω) such that
∫RN∫RN|uk(x)−uk(y)|p−2(uk(x)−uk(y))(ϕ(x)−ϕ(y))K(x,y)dxdy=∫Ωϕ(x)dμk | (2.27) |
for any ϕ∈C∞0(Ω). Since uk→u in RN, we have that ‖u−uk‖Wh,q(RN)→0 for any h∈(0,s) and q∈(0,N(p−1)N−s). This, together with (2.27), implies
∫Ωϕ(x)dμk→∫Ωϕ(x)dμ,∀ϕ∈C∞0(Ω). | (2.28) |
Now, we remark that, in view of the proof of [26, Theorem 1.3], we may apply [26, estimate (1.25)] to uk. Hence,
C−1Wd(x)8s,p[μk](x)≤uk(x),for a.e.x∈Ωand∀k>0. |
Letting k→∞ in the above inequality and using some elementary manipulations, we may obtain the lower estimate in (2.26).
For the upper estimate in (2.26), by [23, Theorem 9], we have that
vk(x):=essliminfy→xuk(y)=uk(x),for a.e.x∈RN. |
Hence, vk is a lower semicontinuous functions in Ω and a nonnegative weak supersolution. By [23, Theorem 12], vk is (s,p)-superharmonic function in Ω (see [23, Definition 1] for the definition of (s,p)-superharmonic function). This, together with [23, Lemma 12], implies that v:=limk→∞vk is (s,p)-superharmonic function in Ω and v=u a.e. in RN. The desired result follows by applying [21, Theorem 5.3] to v and the fact that v=u a.e. in RN.
Proposition 2.10. Let μ∈Mb(Ω). Then there exists a very weak solution u of (2.2) and a positive constant C depending only on N,s,p and ΛK such that
−CW2diam(Ω)s,p[μ−]≤u≤CW2diam(Ω)s,p[μ+],a.e. inΩ. | (2.29) |
Proof. Let u be the solution constructed in Proposition 2.8 and x0∈Ω. Set R=diam(Ω), μn=ρn∗μ and μ⊕n=ρn∗μ+. We denote by v⊕n∈Ws,p0(Ω) the solution of
{Lv⊕n=μ⊕n,inB2R(x0),v⊕n=0,inRN∖B2R(x0). |
By Proposition 2.3, we have that v⊕n≥0 and v⊕n≥un, where un∈Ws,p0(Ω) is the weak solution of (2.2) with μ=μn. By statements (ⅰ)–(ⅲ) in the proof of Proposition 2.8, there exist subsequences {unk,v⊕nk}∞k=1 such that unk→u and v⊕nk→v⊕ a.e. in RN and
‖u−unk‖Wh,q(RN)+‖v⊕−v⊕nk‖Wh,q(RN)→0 |
for any h∈(0,s) and q∈(0,N(p−1)N−s). Combining all above, we may deduce that u≤v⊕ a.e. in RN and v⊕ is a nonnegative very weak solution to
{Lv⊕=μ+,inB2R(x0),v⊕=0,inRN∖B2R(x0). |
In addition, in view of the proof of Proposition 2.9, there exists a positive constant C=C(p,s,ΛK,N) such that
u(x)≤v⊕(x)≤C(WRs,p[μ+](x)+essinfBR2(x)v⊕+Tail(v⊕;x,R2)),for a.e.x∈Ω, | (2.30) |
where
Tail(v⊕;x,R2)=((R2)sp∫RN∖BR2(x)|v⊕(y)|p−1|x−y|N+spdy)1p−1. |
By (2.15) and (2.6), we derive that
essinfBR2(x)v⊕≲(−∫BR2(x)|v⊕|p−1dx)1p−1≲R−N−spp−1μ+(BR(x))1p−1≲W2Rs,p[μ+](x), | (2.31) |
and
Tail(v⊕;x0,R2)≲(−∫B2R(x0)|v⊕|p−1dx)1p−1≲W2Rs,p[μ+](x),∀x∈Ω, | (2.32) |
where the implicit constants in (2.31) and (2.32) depend only on p,s,ΛK,N. The inequalities in (2.32) follow by the fact that v⊕=0 in RN∖B2R(x0) and μ(RN∖Ω)=0.
Combining (2.30)–(2.32), we obtain the upper bound in (2.29).
The proof of the lower bound in (2.29) is similar and we omit it.
We assume that g∈C(R) and rg(r)≥0. Let Ω⊂RN be an open bounded domain and μ∈(Ws,p0(Ω))∗. Set G(r)=∫r0g(s)ds,
J(v)=1p∫RN∫RN|v(x)−v(y)|pK(x,y)dxdy+∫ΩG(v)dx−<μ,v> |
and
XG(Ω)={v∈Ws,p0(Ω):G(v)∈L1(Ω)}. |
Theorem 3.1. Let s∈(0,1), p>1 and μ∈(Ws,p0(Ω))∗. Then, there exists a minimizer uμ of J in XG(Ω). Furthermore, uμ is a weak solution of J, in the sense of
∫RN∫RN|uμ(x)−uμ(y)|p−2(uμ(x)−uμ(y))(ζ(x)−ζ(y))K(x,y)dxdy+∫Ωg(uμ)ζdx=<μ,ζ> | (3.1) |
for any ζ∈Ws,p0(Ω)∩L∞(Ω).
If g is nondecreasing the solution uμ is unique and the mapping μ↦uμ is nondecreasing.
Proof. We adapt the argument used in the proof of [15, Theorem 5.1]. Let {vn} be a minimizing sequence. Taking in to account that G(t)≥0 for any t∈R and the fractional Sobolev inequality, we can easily show the existence of a positive constant C=C(p,Ω,ΛK) such that
‖vn‖pWs,p0(Ω)≤C(J(vn)+‖μ‖p′(Ws,p0(Ω))∗),∀n∈N. | (3.2) |
This implies that vn is uniformly bounded in Ws,p0(Ω). Thus, by the fractional Sobolev embedding theorem (see e.g., [13, Corollary 7.2]) and the fact that Ws,p0(Ω) is a reflexive Banach space, we may prove the existence of a subsequence, still denoted by {vn} and a function v∈Ws,p0(Ω) such that there hold:
(ⅰ) vn→v a.e. in RN.
(ⅱ) vn⇀v in Ws,p0(Ω) and vn→v in Wh,q(RN) for any h∈(0,s) and q∈(1,p).
By Fatou's lemma, we obtain
1p∫RN∫RN|v(x)−v(y)|pK(x,y)dxdy+∫ΩG(v)dx≤lim infk→∞1p∫RN∫RN|vk(x)−vk(y)|pK(x,y)dxdy+∫ΩG(vk)dx. |
Hence v is a minimizer. If g is nondecreasing, the uniqueness of the minimizer follows by the fact that J is strictly convex.
We next show (3.1). Let vk be the minimizer of J associated with gk=max(−k,min(g,k)). Then, in view of the proof of [11, Theorem 2.3], vk satisfies
∫RN∫RN|vk(x)−vk(y)|p−2(vk(x)−vk(y))(ζ(x)−ζ(y))K(x,y)dxdy+∫Ωgk(vk)ζdx=<μ,ζ> | (3.3) |
for any ζ∈Ws,p0(Ω). Taking vk as test function, we have
∫RN∫RN|vk(x)−vk(y)|pK(x,y)dxdy+∫Ωg(vk)vkdx=<μ,vk>≤1p‖vk‖pWs,p0(Ω)+1p′‖μ‖p′(Ws,p0(Ω))∗, |
which implies
∫RN∫RN|vk(x)−vk(y)|p|x−y|N+spdxdy+∫Ωgk(vk)vkdx≤C(ΛK,p)‖μ‖p′(Ws,p0(Ω))∗=:M. | (3.4) |
By the above inequality, we may deduce that there exists a subsequence, still denoted by {vk} and a function v∈Ws,p0(Ω) such that they satisfy statements (ⅰ) and (ⅱ).
Let ζ∈L∞(Ω) with ‖ζ‖L∞(Ω)=N and E⊂Ω be a Borel set. Then, for any λ>0, we have
∫E∩{|vk|>λ}|ζgk(vk)|dx≤1λ∫E∩{vk>λ}|ζ||vkgk(vk)|dx≤Nλ∫Ωvkgk(vk)dx≤MNλ. |
Also,
∫E∩{|uk|≤λ}|ζgk(vk)|dx≤|E|Nsup{|g(t)|:|t|≤λ}. |
Let ε>0, λ=2MNε and δ=ε2Nsup{|g(t)|:|t|≤2MNε}+1. Then for any Borel set E⊂Ω with |E|<δ, we have
∫E|ζgk(vk)|dx<ε. |
Thus, by Vitali's theorem, we conclude
∫Ωgk(vk)ζdx→∫Ωg(v)ζdx. | (3.5) |
Combining all above, we obtain that v satisfies (3.1).
Now for any u∈XG(Ω), we have that u∈XGk(Ω), Gk(u)≤G(u) and
1p∫RN∫RN|vk(x)−vk(y)|pK(x,y)dxdy+∫ΩGk(vk)dx−<μ,vk>≤1p∫RN∫RN|u(x)−u(y)|pK(x,y)dxdy+∫ΩGk(u)dσ−<μ,u>, |
where Gk(r)=∫r0gk(s)ds. By the above inequality and Fatou's Lemma, we deduce that v is a minimizer of J in XG(Ω).
Let g be nondecreasing and uν be the minimizer of J associated with ν∈(Ws,p0(Ω))∗, such that ν≤μ. Then, using vk=min{(uν−uμ)+,k} as test function, we have that
∫RN∫RN|uν(x)−uν(y)|p−2(uν(x)−uν(y))(vk(x)−vk(y))K(x,y)dxdy−∫RN∫RN|uμ(x)−uμ(y)|p−2(uμ(x)−uμ(y))(vk(x)−vk(y))K(x,y)dxdy=−∫Ω(g(uν)−g(uμ))vkdx+<ν−μ,vk>≤0. |
Letting k→∞ in the above inequality and then proceeding as in the proof of [23, Lemma 6], we obtain that uν≤uμ a.e. in RN.
When μ∈Lp′(Ω), we derive the following result which will be useful in the next subsection.
Lemma 3.2. Let μ∈Lp′(Ω), g∈C(RN) be a nondecreasing function with g(0)=0 and u∈Ws,p0(Ω) satisfy (3.1). Then there holds,
∫Ω|g(u)|dx≤∫Ω|μ|dx. | (3.6) |
In addition, if we assume that μ≥0, then u≥0 a.e. in RN.
Proof. Let k>0. Using ϕk=tanh(ku) as test function in (3.1), we obtain
∫RN∫RN|u(x)−u(y)|p−2(u(x)−u(y))(ϕk(x)−ϕk(y))K(x,y)dxdy+∫Ωg(u)ϕkdx=∫Ωμϕkdx. |
If ∞>u(x)>u(y)>−∞, then there exists ξ∈(u(y),u(x)) such that
ϕk(x)−ϕk(y)=(1−tanh2(kξ))(u(x)−u(y))≥c(ξ,k)(u(x)−u(y)). |
Combining the last two displays, we can easily obtain that
∫Ωg(u)ϕkdx≤∫Ω|μ|dx. |
Since g(u)ϕk≥0 a.e. in Ω, by Fatou's lemma and the above inequality, we can easily deduce (3.6).
In this subsection, we always assume that s∈(0,1), 1<p<Ns and g∈C(R) is nondecreasing such that g(0)=0.
Lemma 3.3. Let g∈L∞(R) and λi∈M+b(Ω) (i=1,2). Then there exist very weak solutions u,ui (i=1,2) to problems
{Lu+g(u)=λ1−λ2,inΩ,u=0,inRN∖Ω, | (3.7) |
{Lu1+g(u1)=λ1,inΩ,u=0,inRN∖Ω | (3.8) |
and
{Lu2−g(−u2)=λ2,inΩ,u=0,inRN∖Ω, | (3.9) |
such that there hold
u1,u2≥0and−u2≤u≤u1,a.e. inRN. | (3.10) |
In addition, for any q∈(0,N(p−1)N−s) and h∈(0,s), there exists a positive constant c=c(N,p,s,ΛK,q,h,|Ω|) such that
(∫Ω|g(u)|dx)1p−1+(∫RN∫RN|u(x)−u(y)|q|x−y|N+hqdxdy)1q≤c(λ1(Ω)+λ2(Ω))1p−1 | (3.11) |
and
(∫Ω|g((−1)i+1ui)|dx)1p−1+(∫RN∫RN|ui(x)−ui(y)|q|x−y|N+hqdxdy)1q≤cλi(Ω)1p−1. | (3.12) |
Finally, there exist very weak solutions vi to (2.2) with μ=λi (i = 1, 2) such that
0≤ui≤vi≤CiW2diam(Ω)s,p[λi],a.e. inΩ, | (3.13) |
where Ci is a positive constant depending only on p,s,ΛK and N.
Proof. Let {ρn}∞1 be a sequence of mollifiers and λn,i=ρn∗λi. Then λn,i∈C∞0(RN). By Proposition 3.1, there exist unique solutions un,un,i,vn,i∈Ws,p0(Ω) to the following problems
{Lun+g(un)=λn,1−λn,2,inΩ,u=0,inRN∖Ω, |
{Lun,1+g(un,1)=λn,1,inΩ,u=0,inRN∖Ω, |
{Lun,2−g(−un,2)=λn,2,inΩ,u=0,inRN∖Ω, |
and
{Lvn,i=λn,iinΩu=0,inRN∖Ω, |
such that there holds
−vn,2≤−un,2≤un≤un,1≤vn,1,a.e. inRN. | (3.14) |
By Lemma 3.2 and Proposition 2.7, for any q∈(0,N(p−1)N−s) and h∈(0,s), there exists a positive constant c=c(N,p,s,ΛK,q,h,|Ω|) such that
(∫Ω|g(un)|dx)1p−1+(∫RN∫RN|un(x)−un(y)|q|x−y|N+hqdxdy)1q≤c(∫Ωλn,1+λn,2dx)1p−1, | (3.15) |
(∫Ω|g(un,1)|dx)1p−1+(∫RN∫RN|un,1(x)−un,1(y)|q|x−y|N+hqdxdy)1q≤c(∫Ωλn,1dx)1p−1, | (3.16) |
(∫Ω|g(−un,2)|dx)1p−1+(∫RN∫RN|un,2(x)−un,2(y)|q|x−y|N+hqdxdy)1q≤c(∫Ωλn,2dx)1p−1 | (3.17) |
and
(∫RN∫RN|vn,i(x)−vn,i(y)|q|x−y|N+hqdxdy)1q≤c(∫Ωλn,idx)1p−1. | (3.18) |
Furthermore, in view of the proof of (2.16), we have that Tk(un),Tk(un,i),Tk(vn,i)∈Ws,p0(Ω) and satisfy (2.19) with μ=λ1+λ2.
Since the sequences {λn,i}n are uniformly bounded in Mb(Ω), as in the proof of Proposition 2.8, we may show that there exist subsequences, still denoted by the same index, such that un→u, un,i→ui vn,i→vi in Wh,q(RN) and a.e. in RN. In addition, we may prove that Tk(u),Tk(ui),Tk(vi)∈Ws,p0(Ω) for any k>0. Finally, by dominated convergence theorem, we deduce that g(un)→g(u), g(un,1)→g(u1), g(−un,2)→g(−u2) in L1(Ω). Hence, combining all above, we can easily show that u,ui are very weak solutions of problems (3.7)–(3.9) respectively and vi are very weak solutions of problem (2.2) with μ=λi (i=1,2). By proceeding as in the proof of Proposition 2.10 and using (3.14), we derive (3.13).
Estimates (3.11) and (3.12) follow by (3.15), (3.16) and Fatou's lemma.
Lemma 3.4. Let λi∈M+b(Ω) for i=1,2. We also assume that g((−1)1+iCW2Rs,p[λi])∈L1(Ω), where C is the constant in Proposition 2.10. Then the conclusion of Lemma 3.3 holds true.
Proof. Let Tn(t)=max(−n,min(t,n)) for any n∈N. By Lemma 3.3, there exist very weak solutions un,un,i,vn,i∈Ws,p0(Ω) of the following problems
{Lun+Tnog(un)=λ1−λ2inΩu=0inRN∖Ω, |
{Lun,1+Tnog(un,1)=λ1,inΩ,u=0,inRN∖Ω, |
{Lun,2−Tnog(−un,2)=λ2,inΩ,u=0,inRN∖Ω |
and
{Lvi=λi,inΩ,u=0,inRN∖Ω, |
such that there holds
−CW2diam(Ω)1,p[λ2]≤−v2≤−un,2≤un≤un,1≤v1≤CW2diam(Ω)1,p[λ1],a.e. inRN |
and for any n∈N. The rest of the proof can proceed similarly to the proof of Lemma 3.3 and we omit it.
Proposition 3.5. Assume
Λg:=∫∞1s−˜q−1(g(s)−g(−s))ds<∞ | (3.19) |
for ˜q>0. Let v be a measurable function defined in Ω. For s>0, set
Es(v):={x∈Ω:|v(x)|>s}ande(s):=|Es(v)|. |
Assume that there exists a positive constant C0 such that
e(s)≤C0s−˜q,∀s≥1. | (3.20) |
Then for any s0≥1, there hold
‖g(|v|)‖L1(Ω)≤∫Ω∖Es0(v)g(|v|)dx+˜qC0∫∞s0s−˜q−1g(s)ds,‖g(−|v|)‖L1(Ω)≤−∫Ω∖Es0(v)g(−|v|)dx−˜qC0∫∞s0s−˜q−1g(−s)ds. |
Proof. The proof is very similar to the one of [16, Lemma 5.1] and we omit it.
Proof of Theorem 1.2. Let λ1=μ+ and λ2=μ−. By Lemma 3.3, there exist very weak solutions un,vi of the following problems
{Lun+Tnog(un)=λ1−λ2,inΩ,u=0,inRN∖Ω |
and
{Lvi=λi,inΩ,u=0,inRN∖Ω, |
such that there holds
−v2≤un≤v1,a.e. inRNand∀n∈N. |
Also, taking into consideration that g in nondecreasing with g(0)=0, we may show that Tk(un),Tk(vi) satisfy (2.19) with μ=λ1+λ2. In addition, by (2.15), there holds
‖vp−11‖∗LNN−spw(RN)+‖vp−12‖∗LNN−spw(RN)≤C1(N,p,s,ΛK)(λ1(Ω)+λ2(Ω)). |
By (2.7) and Proposition 3.5, we have that |Tnog(un)|≤g(v1)−g(−v2) and
‖Tnog(un)‖L1(Ω)≤‖g(v1)‖L1(Ω)+‖g(−v2)‖L1(Ω)≤(g(s0)−g(s0))|Ω|+˜qC1(N,p,s,ΛK,Λg)(λ1(Ω)+λ2(Ω))N(p−1)N−sp∫∞s0s−˜q−1(g(s)−g(−s)ds,∀n∈N, |
where ˜q=N(p−1)N−sp. The desired result follows by proceeding as in the proof of Lemma 3.3.
In order to prove Theorem 1.3, we need to introduce some notations concerning the Bessel capacities, we refer the reader to [1] for more detail. For α∈R we define the Bessel kernel of order α by Gα(ξ)=F−1(1+|.|2)−α2(ξ), where F is the Fourier transform of moderate distributions in RN. For any β>1, the Bessel space Lα,β(RN) is given by
Lα,β(RN):={f=Gα∗g:g∈Lβ(RN)}, |
with norm
‖f‖Lα,β(RN):=‖g‖Lβ(RN)=‖G−α∗f‖Lβ(RN). |
The Bessel capacity is defined as follows.
Definition 3.6. Let α>0, 1<β<∞ and E⊂RN. Set
SE:={g∈Lβ(RN):g≥0,Gα∗g(x)≥1for anyx∈E}. |
Then
Capα,β(E):=inf{‖g‖βLβ(RN);g∈SE}. | (3.21) |
If SE=∅, we set Capα,β(E)=∞.
In the sequel, we denote by L−α,β′(RN) the dual of Lα,β(RN) and we set
Gα[μ](x)=∫RNGα(x,y)dμ(y),∀μ∈M(RN). |
Proof of Theorem 1.3. Since μ is absolutely continuous with respect to the capacity Capsp,κκ−p+1, the measures μ+,μ− have the same property. Thus, by [5, Theorem 2.5] (see also [3]), there are nondecreasing sequences {μ±n}n⊂L−sp,κp−1(RN)∩M+b(RN) with compact support in Ω, such that they converge to μ± in the narrow topology. Furthermore, by [5, Theorem 2.3] (see also [1, Corollary 3.6.3]),
‖W2diam(Ω)α,p[μ±n]‖κLκ(RN)≈‖Gsp[μ±n]‖κp−1Lκp−1(RN)<∞. |
By Lemma 3.4, there exist solutions un,un,i,vi to the problems
{Lun+|un|κ−1un=λn,1−λn,2,inΩ,u=0,inRN∖Ω, | (3.22) |
{Lun,1+|un,1|κ−1un,1=λn,1,inΩ,u=0,inRN∖Ω, | (3.23) |
{Lun,2+|un,2|κ−1un,2=λn,2,inΩ,u=0,inRN∖Ω, | (3.24) |
and
{Lvn,i=λn,i,inΩ,u=0,inRN∖Ω, |
such that there holds
−vn,2≤−un,2≤un≤un,1≤vn,1,a.e. inRN. | (3.25) |
Furthermore, in view of the proof of Lemmas 3.3 and 3.4, the sequences {un,i},{vn,i} satisfy (3.15)–(3.18) with g(t)=|t|κsign(t), λn,1=μ+n and λn,2=μ−n, as well as they can be constructed such that
un,i≤un+1,iandvn,i≤vn+1,i,a.e. inRN,∀n∈Nandi=1,2. | (3.26) |
By (3.15) and (3.16) with g(t)=|t|κsign(t), λn,1=μ+n and λn,2=μ−n, we have
∫Ω|un,1|κd≤μ+(Ω)and∫Ω|un,2|κd≤μ−(Ω),∀n∈N. |
By (3.15)–(3.18) with g(t)=|t|κsign(t), λn,1=μ+n and λn,2=μ−n, there are subsequences, still denoted by the same index, such that un→u, un,i→ui vn,i→v in Wh,q(RN) and a.e. in RN. In addition, Tk(u),Tk(ui),Tk(vi)∈Ws,p0(RN) and
∫Ω|u1|kdx≤μ+(Ω),and∫Ω|u2|kdx≤μ−(Ω). |
Therefore, by dominated convergence theorem, we obtain that |un|κ→|u|κ, |un,1|κ→|u1|κ, |un,2|κ→|u2|κ in L1(Ω). This, implies that u,ui are very weak solutions of problems (3.7)–(3.9) respectively and vi are very weak solution of problem (2.2) with μ=λi, where λ1=μ+ and λ2=μ−.
Estimate (1.13) follows by (3.25) and (3.13). Estimate (1.14) follows by (3.15) with g(t)=|t|κsign(t), λn,1=μ+n, λn,2=μ−n and Fatou's lemma.
In this subsection, we investigate the existence of solutions to the following problem
{Lv=g(v)+ρτ,inΩ,v=0,inRN∖Ω, | (4.1) |
where ρ>0, g∈C(R) is a nondecreasing function and
|g(t)|≤a|t|dfor some a>0,d>p−1 and for any |t|≤1. | (4.2) |
Let us state the first existence result.
Lemma 4.1. Let 1<p<Ns and τ∈C∞0(RN) be such that ‖τ‖L1(RN)≤1. Assume that g∈L∞(Ω)∩C(R) satisfies (3.19) for
˜q=N(p−1)N−sp. |
In addition, we assume that g is nondecreasing and satisfies (4.2).
Then there exists a positive constant ρ0 depending on N,Ω,Λg,ΛK,a,d,p,s such that for every ρ∈(0,ρ0), problem (4.1) admits a weak solution v∈Ws,p0(Ω) satisfying
‖|v|p−1‖LNN−spw(Ω)≤t0, | (4.3) |
where t0>0 depends on N,Ω,Λg,ΛK,a,d,p,s.
Proof. We shall use Schauder fixed point theorem to show the existence of a positive weak solution of (4.1).
Let 1<κ<min{NN−sp,dp−1} and v∈L1(Ω). Since g∈L∞(Ω), we can easily show that the following problem
{Lu=g(|v|1p−1sign(v))+ρτ,inΩ,u=0,inRN∖Ω, | (4.4) |
admits a unique weak solution T(v)∈Ws,p0(Ω). We define the operator S by
S(v):=|T(v)|p−1sign(T(v)),∀v∈L1(Ω). | (4.5) |
By (2.8), we obtain
‖S(v)‖LNN−spw(Ω)≤C(s,p,N,ΛK)(ρ∫Ω|τ|dx+∫Ω|g(|v|1p−1sign(v))|dx)≤C(s,p,N,ΛK)(ρ+∫Ωg(|v|1p−1)−g(−|v|1p−1)dx). | (4.6) |
Let v∈LNN−spw(Ω). For any λ>0, we set Eλ:={x∈Ω:|v(x)|1p−1>λ} and e(λ)=∫Eλdx. By (2.4) and (2.6), we can easily show that
e(λ)≤C(N,s,p)‖v‖NN−spLNN−spw(Ω)λ−N(p−1)N−sp. |
By the above inequality and Lemma 3.5 with λ0=1 and ˜q=N(p−1)N−sp, we deduce
∫Ωg(|v|1p−1)−g(−|v|1p−1)dx≤2a∫Ω|v|κdx+C(p,s,N)‖v‖NN−spLNN−spw(Ω)Λg. |
Let λ=‖v‖LNN−spw(Ω). By (2.6), we have that
∫Ω|v|κd=∫∞0|{x∈Ω:|v|≥t}|dtκ=∫λ0|{x∈Ω:|v|≥t}|dtκ+∫∞λ|{x∈Ω:|v|≥t}|dtκ≤|Ω|λκ+κλNN−sp∫∞λtκ−NN−sp−1dt≤C(Ω,κ,s,p,N)λk. |
Combining all above, we may prove that
‖S(v)‖LNN−spw(Ω)≤C(p,N,κ,|Ω|,Λg,ΛK,a)(ρ+‖v‖NN−spLNN−spw(Ω)+‖v‖κLNN−spw(Ω)). |
Therefore, if ‖v‖LNN−spw(Ω)≤t then
‖S(v)‖LNN−spw(Ω)≤C(tNN−sp+tκ+ρ). | (4.7) |
Since 1<κ<NN−sp, there exist t0>0 and ρ0>0 depending on |Ω|,Λg,p,κ,N,a such that for any t∈(0,t0] and ρ∈(0,ρ0), the following inequality holds
C(tNN−sp+tκ+ρ)≤t0, |
where C is the constant in (4.7). Hence,
‖v‖LNN−spw(Ω)≤t0⟹‖S(v)‖LNN−spw(Ω)≤t0. | (4.8) |
Next, we apply Schauder fixed point theorem to our setting.
We claim that S is continuous. First we assume that vn→v in L1(Ω) and T(vn)→T(v) in W1,p0(Ω), then by fractional Sobolev inequality, we have
∫Ω|T(vn)−T(v)|dx≤|Ω|pN−N+spNp‖T(vn)−T(v)‖LNpN−sp(Ω)≤C|Ω|pN−N+spNp‖T(vn)−T(v)‖W1,p0(Ω)→0. | (4.9) |
Let k>0 and ε>0, then
∫Ω|S(vn)−S(v)|dx=∫{x∈Ω:|S(vn)(x)|≤k}∩{x∈Ω:|S(v)(x)|≤k}|S(vn)−S(v)|dx+∫Ω∖({x∈Ω:|S(vn)(x)|≤k}∩{x∈Ω:|S(v)(x)|≤k})|S(vn)(x)−S(v)(x)|dx. | (4.10) |
By (4.6) and the fact that g∈L∞(R), we have that S(vn)∈Lβ(Ω) and {S(vn)} is uniformly bounded in Lβ(Ω) for any β∈(1,NN−sp). Hence, there exists k0∈N, such that
∫Ω∖({x∈Ω:|S(vn)(x)|≤k}∩{x∈Ω:|S(v)(x)|≤k})|S(vn)−S(v)|dx≤ε3∀k≥k0,andn∈N. | (4.11) |
Now, we set
Ak0,n={x∈Ω:|T(vn)(x)|≤k1p−10}∩{x∈Ω:|T(v)(x)|≤k1p−10} |
and Bδ,n={x∈Ω:|T(v)(x)−T(vn)(x)|≤δ}. Then, we have that
∫Ω∩{x∈Ω:|S(vn)|≤k0}∩{x∈Ω:|S(v)|≤k0}|S(vn)−S(v)|dx=∫Ak0,n∩Bδ,n||T(vn)|p−1sign(T(vn))−|T(v)|p−1sign(T(v))|dx+∫Ak0,n∖Bδ,n||T(vn)|p−1sign(T(vn))−|T(v)|p−1sign(T(v))|dx. | (4.12) |
Since h(t)=tp−1sign(t) is uniformly continuous in [−k0,k0], there exists δ0>0 independent of n such that
∫Ak0,n∩Bδ0,n||T(vn)|p−1sign(T(vn))−|T(v)|p−1sign(T(v))|dx≤ε3. | (4.13) |
Moreover, by (4.9), there exists n0=n0(δ0,k0,p)∈N such that
∫Ak0,n0∖Bδ0,n0||T(vn0)|p−1sign(T(vn0))−|T(v)|p−1sign(T(v))|dx≤ε3. | (4.14) |
Hence, combining (4.9)–(4.14), we obtain that S(vn)→S(v) in L1(Ω).
Therefore, it is enough to show that T(vn)→T(v) in Ws,p0(Ω). In order to prove this, we will consider two cases.
Case 1. 1<p<2. Let M:=supt∈R|g(t)|. We will show that T(vn)→T(v) in Ws,p0(Ω). Since T(vn),T(v)∈Ws,p0(Ω) are weak solutions of (4.4) with vn and v respectively, we have
∫RN∫RN|T(vn)(x)−T(vn)(y)|pK(x,y)dxdy=∫ΩT(vn)(g(|vn|1p−1sign(vn))dx+∫ΩT(vn)τdx≤M|Ω|p−1p(∫Ω|T(vn)|pdx)1p+(∫Ω|T(vn)|pdx)1p(∫Ω|τ|pp−1dx)p−1p≤C1(M,Ω,p,N,τ,s)(∫RN∫RN|T(vn)(x)−T(vn)(y)|p|x−y|N+spdxdy)1p. | (4.15) |
Therefore,
∫RN∫RN|T(vn)(x)−T(vn)(y)|p|x−y|N+spdxdy≤Cp−1p1(M,Ω,p,N,τ,s,ΛK). | (4.16) |
Using ϕ=T(vn)−T(v) as test function, we have
I:=∫RN∫RN|T(vn)(x)−T(vn)(y)|p−2(T(vn)(x)−T(vn)(y))(ϕ(x)−ϕ(y))K(x,y)dxdy−∫RN∫RN|T(v)(x)−T(v)(y)|p−2(T(v)(x)−T(v)(y))(ϕ(x)−ϕ(y))K(x,y)dxdy=∫Ωϕ(g(|vn|1p−1sign(vn))−g(|v|1p−1sign(v)))dx=:II. | (4.17) |
We first treat I. On one hand, since
(|a|p−2a−|b|p−2b)(a−b)≥C(p)|a−b|2(|a|+|b|)2−p |
for any (a,b)∈R2N∖{(0,0)} and p∈(1,2), we have
I≥C(p)∫RN∫RN|ϕ(x)−ϕ(y)|2(|T(vn)(x)−T(vn)(y)|+|T(v)(x)−T(v)(y)|)p−2K(x,y)dxdy. | (4.18) |
On the other hand, by Hölder inequality, we obtain
∫RN∫RN|ϕ(x)−ϕ(y)|p|x−y|N+spdxdy≤ΛK∫RN∫RN|ϕ(x)−ϕ(y)|pK(x,y)dxdy≤C(p,ΛK)(∫RN∫RN(|T(vn)(x)−T(vn)(y)|+|T(v)(x)−T(v)(y)|)pK(x,y)dxdy)2−p2Ip2≤C(p,C1,Ω,ΛK)Ip2, | (4.19) |
where C1 is the constant in (4.16). Hence, by (4.18) and (4.19), we obtain
C(∫RN∫RN|ϕ(x)−ϕ(y)|p|x−y|N+spdxdy)2p≤I. | (4.20) |
Next we treat II. Let r=NpN−sp, proceeding as in the proof of (4.15), we have
II≤(∫Ω|ϕ|rd)1r(∫Ω|g(|vn|1p−1sign(vn))−g(|v|1p−1sign(v))|r′d)1r′≤C(N,p,s)(∫RN∫RN|ϕ(x)−ϕ(y)|p|x−y|N+spdxdy)1p×(∫Ω|g(|vn|1p−1sign(vn))−g(|v|1p−1sign(v))|r′d)1r′, | (4.21) |
where in the last inequality we used the fractional Sobolev inequality.
Combining (4.17), (4.20) and (4.21), we obtain
(∫RN∫RN|ϕ(x)−ϕ(y)|p|x−y|N+spdxdy)1p≤C(p,C1,Ω,s,ΛK)(∫Ω|g(|vn|1p−1sign(vn))−g(|v|1p−1sign(v))|r′d)1r′. | (4.22) |
Since g∘(|⋅|p−1sign(⋅)) is uniformly continuous in R, bounded and vn→v in L1(Ω), we obtain
limn→∞∫Ω|g(|v|1p−1sign(v))−g(|vn|1p−1sign(vn))|r′dx=0, |
which, together with (4.22), implies the desired result.
Case 2. p≥2. We note here that
(|a|p−2a−|b|p−2b)(a−b)≥C(p)|a−b|p |
for any (a,b)∈R2N and p≥2. Thus,
I≥C(N,p,ΛK)(∫RN∫RN|ϕ(x)−ϕ(y)|p|x−y|N+spdxdy). |
By using a similar argument to the one in Case 1, we may show that T(vn)→T(v) in Ws,p0(Ω).
Next we claim that S is compact. Indeed, let {vn} be a sequence in L1(Ω) then by (4.16), we obtain that T(vn) is uniformly bounded in Ws,p0(Ω). Hence there exists a subsequence still denoted by {T(vn)} such that T(vn)⇀ψ in Ws,p0(Ω) and T(vn)→ψ a.e. in RN. Furthermore, in view of (4.6), we can easily show that S(vn)=|T(vn)|p−1sign(T(vn))→|ψ|p−1sign(ψ) in L1(Ω).
Now set
O:={v∈L1(Ω):‖v‖LNN−spw(Ω)≤t0}. | (4.23) |
Then O is a closed, convex subset of L1(Ω) and by (4.8), S(O)⊂O. Thus we can apply Schauder fixed point theorem to obtain the existence of a function v∈O such that S(v)=v. This means that u=v1p−1sign(v) is a solution of (4.1) satisfying (4.3).
Proof of Theorem 1.4. Let {ρn}∞n=1 be a sequence of mollifiers. Set τn=ρn∗τ and gn=max(−n,min(g,n)). Then gn satisfies (1.9) with the same constant Λg. Thus, there exists a weak solution un∈Ws,p0(Ω) of
{Lv=gn(v)+ρτn,inΩ,v=0,inRN∖Ω. |
In addition, it satisfies
‖|un|p−1‖LNN−spw(Ω)≤t0, | (4.24) |
where t0>0 depends on N,Ω,Λg,ΛK,a,s,p,d.
By (4.24), we have that
|{x∈Ω:|un|>s}|≤tNN−sp0s−N(p−1)N−sp. |
Hence by Proposition 3.5,
∫Ω|gn(un)|dx≤C,∀n∈N, |
where C depends only on N,Ω,Λg,ΛK,a,s,p,d and t0. This, together with Proposition 2.7, implies that for any q∈(p−1,N(p−1)N−s) and h∈(0,s), there exists a positive constant c=c(N,s,p,ΛK,s,h,q,|Ω|) such that
(∫RN∫RN|un(x)−un(y)|q|x−y|N+hqdxdy)1q≤c(C+ρ∫Ω|τn|dx)1p−1. | (4.25) |
Therefore, in view of the proof of Proposition 2.8, we may show that there exists a subsequence, still denoted by the same notation, such that un→u in Wh,q(RN) and a.e. in RN. Now, we will show that gn(un)→g(u) in L1(Ω). We will prove it by using Vitali's convergence theorem. Let E⊂Ω be a Borel set. Then, by Lemma 3.5 and (4.24), we have
∫E|gn(un)|dx≤∫Ω|g(un)|dx≤(g(s0)−g(−s0))|E|+C(t0,p,Λg,N)∫∞s0(g(s)−g(−s))s−1−(p−1)NN−spds,∀s0≥1. |
Let ε>0, then there exists s0 such that
C(t0,p,Λg,N)∫∞s0(g(s)−g(−s))s−1−(p−1)NN−spds≤ε2. |
Set δ=ε2(1+g(s0)−g(−s0))>0. Then for any Borel set E with |E|≤δ, we have
g(s0)|E|≤ε2. |
Hence, by the last three inequalities, we may invoke Vitali's convergence theorem in order to prove that gn(un)→g(u) in L1(Ω).
In view of the proof of Proposition 2.8, we may deduce that u is a very weak solution of (1.16). Furthermore, by Fatou's lemma, we can easily show that u satisfies (1.17) and (1.18).
Proof of Proposition 1.5. Let w=ACW2diam(Ω)s,p[ρτ], where C is the constant in (2.29) and A>1 is a constant that will be determined later. Set dν=wκdx+ρdτ, then by (1.20), we obtain
CW2(diam(Ω))s,p[ν]≤21p−1C(W2diam(Ω)s,p[wκ]+W2diam(Ω)s,p[ρdτ])≤21p−1C((AC)κp−1ρκ(p−1)2MW2diam(Ω)s,p[τ]+W2diam(Ω)s,p[ρdτ])≤21p−1C((AC)κp−1Mρκ−p+1(p−1)2+1)W2diam(Ω)s,p[ρdτ]. |
If we choose A=21p−1+1 and ρ small enough such that (AC)κp−1Mρκ−p+1(p−1)2+1<2, we deduce that
CW2diam(Ω)s,p[ν]≤w. | (4.26) |
Now, let x0∈Ω be such that W2diam(Ω)s,p[τ](x0)<∞. If 0≤v≤c0W2diam(Ω)s,p[τ] a.e. in RN, for some constant c0>0, then we have
(∫Ω|v|κdx)1p−1≤(∫Bdiam(Ω)(x0)|v|κdx)1p−1≤C(Ω,N,s,p,M,K,c0)W2diam(Ω)s,p[τ](x0)<∞. |
Thus v∈Lκ(Ω).
Let u0≥0 be a very weak solution of
{Lu0=ρτ,inΩ,v=0,inRN∖Ω, |
satisfying C−1Wd(x)8s,p[μn−1](x)≤u0(x)≤CW2diam(Ω)s,p[ρτ](x) a.e. in Ω. We may construct a nondecreasing sequence {un}n≥0, such that un is a very weak solution to problem
{Lun=uκn−1+ρτ,inΩ,v=0,inRN∖Ω |
and satisfies
C−1Wd(x)8s,p[μn−1](x)≤un(x)≤CW2diam(Ω)s,p[μn−1](x),for a.e.x∈Ω, |
for any n∈N, where dμn−1=uκn−1dx+ρdτ. In addition, by (4.26) and the above inequality, there holds
C−1Wd(x)8s,p[μn−1](x)≤un(x)≤w(x),for a.e.x∈Ω, | (4.27) |
where the positive constant C−1 depends only on N,p,s,q. Finally, un satisfies (2.15)–(2.17) with dμ=wκdx+ρdτ.
Proceeding as in the proof of Proposition 2.8, we may show that there exists a subsequence, still denoted by {un}, such that un→u a.e. in RN and u is a very weak solution of problem (1.19). By (4.27) and Fatou's Lemma, we obtain estimate (1.21). The proof is complete.
Proof of Theorem 1.6. We will first prove that (i) implies (ii) by using some ideas from [30]. Without loss of generality we assume that ρ=1. Extend μ to whole RN by setting μ(RN∖Ω)=0.
Let 0≤g∈Lκp−1(RN;μ). We set
Mμg(x):=supr>0,μ(B(x,r))≠0μ(B(x,r))−1∫B(x,r)g(y)dμ. |
It is well known that there exists a positive constant c1 depending only on N,p,κ such that
∫RN(Mμg(x))κp−1dμ≤c1∫RN|g(x)|κp−1dμ | (4.28) |
(see, e.g., [14]). Also,
∫Ω(Wd(x)8s,p[gμ](x))κdx≤∫Ω(Wd(x)8s,p[μ](x))κ(Mμg(x))κp−1dx≤Cκ∫Ωuκ(x)(Mμg(x))κp−1dx≤Cκ∫Ω(Mμg(x))κp−1dμ≤c2∫RN|g(x)|κp−1dμ. | (4.29) |
Let K=suppτ. By the assumption, we have that r0:=dist(K,∂Ω)>0. Set g=fanxiexianmyfhK˜g, for any nonnegative ˜g∈Lκp−1(RN;μ⌊K). We first note that Br08(x)∩K=∅ if x∈Ω with d(x)<r08 or if x∈RN∖Ω, which implies
Wd(x)24s,p[˜gμ⌊K](x)=0, |
if x∈Ω with d(x)<r024 or if x∈RN∖Ω. Therefore, by the above equality and (4.29), we have
∫RN(Wr024s,p[˜gμ⌊K](x))κdx≤∫Ω(Wd(x)8s,p[gμ](x))κdx≤c2∫RN|˜g(x)|κp−1dμ⌊K. |
Also, by [5, Theorem 2.3] (see also [1, Corollary 3.6.3]), we have
∫RN(Wr024s,p[˜gμ⌊K](x))κdx≈∫RN(Gsp[˜gμ⌊K])κp−1dx, | (4.30) |
where the implicit constant depends only on s,p,N,κ and r0.
Hence, combining the last two displays, we may show that there exists a positive constant c3=c3(N,p,s,κ,r0) such that
∫RN(Gsp[˜gμ⌊K])κp−1dx≤c3∫RN|˜g(x)|κp−1dμ⌊K. | (4.31) |
Let f∈Lκκ−p+1(RN). Then, for any ˜g∈Lκp−1(RN;μ⌊K), there holds
|∫RNf(x)Gsp∗(˜gμ⌊K)(x)dx|=|∫RN˜g(y)Gsp∗f(y)dμ⌊K|≤C1‖f‖Lκκ−p+1(RN)‖˜g‖Lκp−1(RN;μ⌊K). |
The last inequality implies,
∫RN|Gsp∗f(x)|κκ−p+1dμ⌊K≤c4∫RN|f|κκ−p+1dx∀f∈Lκκ−p+1(RN). |
By [1, Theorem 7.2.1], the above inequality is equivalent to
μ⌊K(F)≤c5Capsp,κκ−p+1(F), | (4.32) |
for any compact F⊂RN. (1.23) follows by the above inequality and the fact that \tau\leq \mu_{\lfloor K}.
Next, we prove that (ⅱ) implies (ⅲ). We note that proceeding as above, in the opposite direction, we may prove that (1.23) implies
\int_{{\mathbb R}^N}({\mathbb G}_{sp}[ \tilde g\tau])^\frac{ \kappa}{p-1}{{\mathrm{d}}} x \leq c_3\int_{{\mathbb R}^N} |\tilde g(x)|^\frac{ \kappa}{p-1}{{\mathrm{d}}} \tau,\quad\forall \tilde g\in L^\frac{ \kappa}{p-1}({\mathbb R}^N;\tau). |
By (4.30) and taking \tilde g = fanxiexian_myfh_B, we can easily show that there exists a positive constant C depending only on N, s, p, \Omega such that
\int_{{\mathbb R}^N} (W_{s,p}^{2 \mathrm{diam}\,( \Omega)}[\tau_{\lfloor B}])^ \kappa {{\mathrm{d}}} x \leq C\tau(B). |
We will show that (ⅲ) implies (ⅳ). Let R = 2 \mathrm{diam}\, (\Omega) and C_3 be the constant in (1.24). In the spirit of the proof of [31, Theorem 2.10], we need to prove that there exists a positive constant c_0 = c_0(N, p, \kappa, s, C_3, R, \tau(\Omega)) > 0 such that
\begin{align} \tau(B_t(x))\leq c_0t^{\frac{ \kappa(N-sp)-N(p-1)}{ \kappa-p+1}} \end{align} | (4.33) |
for any t\leq R and \forall x\in \Omega.
Concerning the proof of the above inequality, we first note that for any y\in B_t(x) and t\leq \frac{R}{4}, there holds
\begin{align*} W_{s,p}^{R}[\tau_{\lfloor B_t(x)}](y)& = \int_0^{R}\left(\frac{\tau(B_r(y)\cap B_t(x))}{r^{N-sp}}\right)^{\frac{1}{p-1}}\frac{{{\mathrm{d}}} r}{r}\geq \int_{2t}^{4t}\left(\frac{\tau(B_r(y)\cap B_t(x))}{r^{N-sp}}\right)^{\frac{1}{p-1}}\frac{{{\mathrm{d}}} r}{r}\\ &\geq C(N,p,s)\left(\frac{\tau(B_t(x))}{t^{N-sp}}\right)^{\frac{1}{p-1}}. \end{align*} |
By the above inequality, we deduce
\begin{align} \begin{aligned} t^NC^ \kappa(N,p,s)\left(\frac{\tau(B_t(x))}{t^{N-sp}}\right)^{\frac{ \kappa}{p-1}}&\leq \int_{B_t(x)}(W_{s,p}^{R}[\tau_{\lfloor B_t(x)}](y))^ \kappa {{\mathrm{d}}} y\\ &\leq C_3\tau(B_t(x)),\quad\forall t\in(0\frac{R}{4}], \end{aligned} \end{align} | (4.34) |
where in the last inequality we used (1.24). This implies (4.33).
For any x\in \Omega and t < R, we set
\nu_t(x): = \int_{B_t(x)}\bigg(\int_0^t\left(\frac{\tau(B_r(y))}{r^{N-sp}}\right)^{\frac{1}{p-1}}\frac{{{\mathrm{d}}} r}{r}\bigg)^ \kappa{{\mathrm{d}}} y |
and
\mu_t(x): = \int_{B_t(x)}\bigg(\int_t^{R}\left(\frac{\tau(B_r(y))}{r^{N-sp}}\right)^{\frac{1}{p-1}}\frac{{{\mathrm{d}}} r}{r}\bigg)^ \kappa{{\mathrm{d}}} y. |
Then we can easily prove that
\begin{align} W_{s,p}^{R}[(W_{s,p}^{R}[\tau])^ \kappa]\leq C(q,p)\bigg(\int_0^{R}\left(\frac{ \nu_t(x)}{t^{N-sp}}\right)^\frac{1}{p-1}\frac{{{\mathrm{d}}} t}{t} +\int_0^{R}\left(\frac{ \mu_t(x)}{t^{N-sp}}\right)^\frac{1}{p-1}\frac{{{\mathrm{d}}} t}{t}\bigg). \end{align} | (4.35) |
Now, we treat the first term on the right hand in (4.35). By (1.24), we have
\begin{align} \begin{aligned} \nu_t(x)& = \int_{B_t(x)}\bigg(\int_0^t\left(\frac{\tau(B_r(y)\cap B_{2t}(x))}{r^{N-sp}}\right)^{\frac{1}{p-1}}\frac{{{\mathrm{d}}} r}{r}\bigg)^ \kappa{{\mathrm{d}}} y\leq C\tau(B_{2t}(x)), \end{aligned} \end{align} | (4.36) |
which implies
\begin{align} \int_0^{R}\left(\frac{ \nu_t(x)}{t^{N-sp}}\right)^\frac{1}{p-1}\frac{{{\mathrm{d}}} t}{t}\leq CW_{s,p}^{2R}[\tau](x). \end{align} | (4.37) |
Next, we treat the second term on the right hand in (4.35). First we note that
\begin{align*} \mu_t(x)&\leq \int_{B_t(x)}\bigg(\int_t^{R}\left(\frac{\tau(B_{2r}(x))}{r^{N-sp}}\right)^{\frac{1}{p-1}}\frac{{{\mathrm{d}}} r}{r}\bigg)^ \kappa{{\mathrm{d}}} y\\ &\leq C(N) t^N\bigg(\int_t^{2R}\left(\frac{\tau(B_{r}(x))}{r^{N-sp}}\right)^{\frac{1}{p-1}}\frac{{{\mathrm{d}}} r}{r}\bigg)^ \kappa = :C(N)t^N \mu_{1,t}^ \kappa(x), \end{align*} |
which implies
\begin{align*} \int_0^{R}\left(\frac{ \mu_t(x)}{t^{N-sp}}\right)^\frac{1}{p-1}\frac{{{\mathrm{d}}} t}{t}&\leq C(N,p) \int_0^{R} \mu_{1,t}^\frac{ \kappa}{p-1}(x)t^{\frac{sp}{p-1}-1}{{\mathrm{d}}} t\\ & = C(N,p,s,q)\bigg( \mu_{1,R}^{\frac{ \kappa}{p-1}}(x)R^{\frac{sp}{p-1}}+ \int_0^{R}\left( \mu_{1,t}(x)\right)^{\frac{ \kappa}{p-1}-1}t^{\frac{sp}{p-1}} \left(\frac{\tau(B_t(x))}{t^{N-sp}}\right)^{\frac{1}{p-1}}\frac{{{\mathrm{d}}} t}{t}\bigg), \end{align*} |
where we have used integration by parts in the last equality. By (4.33), we have
\mu_{1,t}(x) = \int_t^{2R}\left(\frac{\tau(B_{r}(x))}{r^{N-sp}}\right)^{\frac{1}{p-1}}\frac{{{\mathrm{d}}} r}{r}\leq C t^{-\frac{sp}{ \kappa-p+1}}. |
Combining the last two displays, we obtain
\begin{align} \int_0^{R}\left(\frac{ \mu_t(x)}{t^{N-sp}}\right)^\frac{1}{p-1}\frac{{{\mathrm{d}}} t}{t}\leq C(N,p,s, \kappa,R)\bigg(\tau(B_{2R}(x))^{\frac{ \kappa}{(p-1)^2}}+W_{s,p}^{R}[\tau](x)\bigg). \end{align} | (4.38) |
The desired result follows by (4.35), (4.37), (4.38) and the fact that
\tau(B_{2R}(x))^{\frac{ \kappa}{(p-1)^2}} \leq \tau( \Omega)^{\frac{ \kappa-p+1}{(p-1)^2}}\tau(B_{\frac{R}{2}}(x))^{\frac{1}{p-1}}\leq C(R,N,p,\tau,s, \kappa)W_{s,p}^{R}[\tau](x),\quad\forall x\in \Omega. |
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The author wishes to thank Professor L. Véron for useful discussions. The author would like to thank the anonymous referee for a careful reading of the manuscript and helpful comments. The research project was supported by the Hellenic Foundation for Research and Innovation (H.F.R.I.) under the "2nd Call for H.F.R.I. Research Projects to support Post-Doctoral Researchers" (Project Number: 59).
The author declares no conflict of interest in this paper.
[1] | D. R. Adams, L. I. Hedberg, Function spaces and potential theory, Grundlehren der mathematischen Wissenschaften, Vol. 314, Springer, 1996. https://doi.org/10.1007/978-3-662-03282-4 |
[2] |
B. Abdellaoui, A. Attar, R. Bentifour, On the fractional p-laplacian equations with weights and general datum, Adv. Nonlinear Anal., 8 (2019), 144–174. https://doi.org/10.1515/anona-2016-0072 doi: 10.1515/anona-2016-0072
![]() |
[3] |
P. Baras, M. Pierre, Singularités éliminables pour des équations semi-linéaires, Ann. Inst. Fourier, 34 (1984), 185–206. https://doi.org/10.5802/aif.956 doi: 10.5802/aif.956
![]() |
[4] |
M. F. Bidaut-Véron, Removable singularities and existence for a quasilinear equation with absorption or source term and measure data, Adv. Nonlinear Stud., 3 (2003), 25–63. https://doi.org/10.1515/ans-2003-0102 doi: 10.1515/ans-2003-0102
![]() |
[5] |
M. F. Bidaut-Véron, Q. H. Nguyen, L. Véron, Quasilinear Lane-Emden equations with absorption and measure data, J. Math. Pures Appl., 102 (2014), 315–337. https://doi.org/10.1016/j.matpur.2013.11.011 doi: 10.1016/j.matpur.2013.11.011
![]() |
[6] |
L. Boccardo, T. Gallouët, Nonlinear elliptic and parabolic equations involving measure data, J. Funct. Anal., 87 (1989), 149–169. https://doi.org/10.1016/0022-1236(89)90005-0 doi: 10.1016/0022-1236(89)90005-0
![]() |
[7] | H. Chen, P. Felmer, L. Véron, Elliptic equations involving general subcritical source nonlinearity and measures, arXiv, 2014. https://doi.org/10.48550/arXiv.1409.3067 |
[8] |
H. Chen, A. Quaas, Classification of isolated singularities of nonnegative solutions to fractional semilinear elliptic equations and the existence results, J. Lond. Math. Soc., 97 (2018), 196–221. https://doi.org/10.1112/jlms.12104 doi: 10.1112/jlms.12104
![]() |
[9] |
H. Chen, L. Véron, Semilinear fractional elliptic equations involving measures, J. Difer. Equations, 257 (2014), 1457–1486. https://doi.org/10.1016/j.jde.2014.05.012 doi: 10.1016/j.jde.2014.05.012
![]() |
[10] | G. Dal Maso, F. Murat, L. Orsina, A. Prignet, Renormalized solutions of elliptic equations with general measure data, Ann. Sc. Norm. Super. Pisa, Cl. Sci., 28 (1999), 741–808. |
[11] |
A. Di Castro, T. Kuusi, G. Palatucci, Local behavior of fractional p-minimizers, Ann. Inst. H. Poincaré Anal. Non Linéaire, 33 (2016), 1279–1299. https://doi.org/10.1016/j.anihpc.2015.04.003 doi: 10.1016/j.anihpc.2015.04.003
![]() |
[12] |
A. Di Castro, T. Kuusi, G. Palatucci, Nonlocal Harnack inequalities, J. Funct. Anal., 267 (2014), 1807–1836. https://doi.org/10.1016/j.jfa.2014.05.023 doi: 10.1016/j.jfa.2014.05.023
![]() |
[13] |
E. Di Nezza, G. Palatucci, E. Valdinoci, Hitchhiker's guide to the fractional Sobolev spaces, Bull. Sci. Math., 136 (2012), 521–573. https://doi.org/10.1016/j.bulsci.2011.12.004 doi: 10.1016/j.bulsci.2011.12.004
![]() |
[14] |
R. Fefferman, Strong differentiation with respect to measures, Amer. J. Math., 103 (1981), 33–40. https://doi.org/10.2307/2374188 doi: 10.2307/2374188
![]() |
[15] |
K. T. Gkikas, Quasilinear elliptic equations involving measure valued absorption terms and measure data, JAMA, 2023. https://doi.org/10.1007/s11854-023-0321-0 doi: 10.1007/s11854-023-0321-0
![]() |
[16] |
K. Gkikas, P. T. Nguyen, Semilinear elliptic Schrödinger equations involving singular potentials and source terms, Nonlinear Anal., 238 (2024), 113403. https://doi.org/10.1016/j.na.2023.113403 doi: 10.1016/j.na.2023.113403
![]() |
[17] | J. Heinonen, T. Kilpeläinen, O. Martio, Nonlinear potential theory of degenerate elliptic equations, Oxford Mathematical Monographs, Oxford Science Publications, 1993. |
[18] |
A. Iannizzotto, S. Mosconi, M. Squassina, Global Hölder regularity for the fractional p-Laplacian, Rev. Mat. Iberoam., 32 (2016), 1353–1392. https://doi.org/10.4171/rmi/921 doi: 10.4171/rmi/921
![]() |
[19] |
T. Kilpeläinen, J. Malý, The Wiener test and potential estimates for quasilinear elliptic equations, Acta Math., 172 (1994), 137–161. https://doi.org/10.1007/BF02392793 doi: 10.1007/BF02392793
![]() |
[20] | T. Kilpeläinen, J. Malý, Degenerate elliptic equations with measure data and nonlinear potentials, Ann. Scuola Norm. Sup. Pisa Cl. Sci., 19 (1992), 591–613. |
[21] |
M. Kim, K. A. Lee, S. C. Lee, The Wiener criterion for nonlocal Dirichlet problems, Commun. Math. Phys., 400 (2023), 1961–2003. https://doi.org/10.1007/s00220-023-04632-w doi: 10.1007/s00220-023-04632-w
![]() |
[22] |
J. Korvenpää, T. Kuusi, E. Lindgren, Equivalence of solutions to fractional p-Laplace type equations, J. Math. Pures Appl., 132 (2019), 1–26. https://doi.org/10.1016/j.matpur.2017.10.004 doi: 10.1016/j.matpur.2017.10.004
![]() |
[23] |
J. Korvenpää, T. Kuusi, G. Palatucci, Fractional superharmonic functions and the Perron method for nonlinear integro-differential equations, Math. Ann., 369 (2017), 1443–1489. https://doi.org/10.1007/s00208-016-1495-x doi: 10.1007/s00208-016-1495-x
![]() |
[24] |
J. Korvenpää, T. Kuusi, G. Palatucci, The obstacle problem for nonlinear integro-differential operators, Calc. Var. PDEs, 55 (2016), 63. https://doi.org/10.1007/s00526-016-0999-2 doi: 10.1007/s00526-016-0999-2
![]() |
[25] | T. Kuusi, G. Mingione, Y. Sire, Regularity issues involving the fractional p-Laplacian, In: Recent developments in nonlocal theory, Berlin: De Gruyter, 2017,303–334. https://doi.org/10.1515/9783110571561-010 |
[26] |
T. Kuusi, G. Mingione, Y. Sire, Nonlocal equations with measure data, Commun. Math. Phys., 337 (2015), 1317–1368. https://doi.org/10.1007/s00220-015-2356-2 doi: 10.1007/s00220-015-2356-2
![]() |
[27] |
E. Lindgren, P. Lindqvist, Perron's method and Wiener's theorem for a nonlocal equation, Potential Anal., 46 (2017), 705–737. https://doi.org/10.1007/s11118-016-9603-9 doi: 10.1007/s11118-016-9603-9
![]() |
[28] |
E. Lindgren, P. Lindqvist, Fractional eigenvalues, Calc. Var. PDEs, 49 (2014), 795–826. https://doi.org/10.1007/s00526-013-0600-1 doi: 10.1007/s00526-013-0600-1
![]() |
[29] |
G. Palatucci, The dirichlet problem for the p-fractional laplace equation, Nonlinear Anal., 177 (2018), 699–732. https://doi.org/10.1016/j.na.2018.05.004 doi: 10.1016/j.na.2018.05.004
![]() |
[30] |
N. C. Phuc, I. E. Verbitsky, Singular quasilinear and Hessian equations and inequalities, J. Funct. Anal., 256 (2009), 1875–1906. https://doi.org/10.1016/j.jfa.2009.01.012 doi: 10.1016/j.jfa.2009.01.012
![]() |
[31] | N. C. Phuc, I. E. Verbitsky, Quasilinear and Hessian equations of Lane-Emden type, Ann. Math., 168 (2008), 859–914. |
[32] | L. Véron, Local and global aspects of quasilinear degenerate elliptic equations, Quasilinear elliptic singular problems, World Scientific Publishing Co. Pte. Ltd., 2017. https://doi.org/10.1142/9850 |