Processing math: 96%
Research article

Nonlinear nonlocal equations involving subcritical or power nonlinearities and measure data

  • Received: 13 October 2023 Revised: 04 January 2024 Accepted: 09 January 2024 Published: 15 January 2024
  • 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 gC(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; κ>p1, 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

    Related Papers:

    [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 gC(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; κ>p1, 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Ω, (P±)

    where μ is a bounded Radon measure on Ω and gC(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)|p2(u(x)u(y))K(x,y)dy,xΩ,

    where the symbol P.V. stands for the principle value integral and K:RN×RNR is a measurable and symmetric (i.e., K(x,y)=K(y,x)) function. Note that if K(x,y)|xy|Nsp then L coincides with the standard fractional p-Laplace operator (Δ)sp.

    Throughout this work, we assume that there exists a positive constant ΛK1 such that the following ellipticity condition holds

    Λ1K|xy|NspK(x,y)ΛK|xy|Nsp,(x,y)RN×RNandxy.

    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)Γ(1s)>0.

    For p=2 and K(x,y)=CN,s|xy|N2s, 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))sNN2s1ds<,

    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|p2u), 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 κ>p1. 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)1p1drr, (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 κ>p1 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 gWs,p(RN), LΦ is a nonlocal operator defined by

    <LΦu,ζ>:=RNRNΦ(u(x)u(y))(ζ(x)ζ(y))K(x,y)dydx,ζC0(Ω).

    Here Φ:RR is a continuous function such that Φ(0)=0 and

    Λ1Φ|t|pΦ(t)tΛΦ|t|p.

    When 2sN<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|p2t 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; κ>p1, 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;κ>p1. 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, ˜gC(R), ΩRN be an open bounded domain and μM(Ω). We will say that u:RNR 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 uWh,q(RN) for any 0<h<s and for any 0<q<N(p1)Ns.

    (ⅱ) Tk(u):=max(k,min(u,k))Ws,p0(Ω) for any k>0.

    (ⅲ)

    RNRN|u(x)u(y)|p2(u(x)u(y))(ϕ(x)ϕ(y))K(x,y)dxdy+Ω˜g(u)ϕdx=Ωϕdμ

    for any ϕC0(Ω).

    We note here that if 2sN<p<Ns, then the very weak solution u belongs to the fractional Sobolev space Wh,q(RN) for any q(1,N(p1)Ns). If p2sN, 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 ˜g0. 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 ˜g0 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|p1LNNspw(RN)C(N,p,s,ΛK)Ω|μ|dx, (1.8)

    where LNNspw(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 u0 and

    C1(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 gC(R) is a nondecreasing function satisfying g(0)=0 and

    1(g(s)g(s))sN(p1)Nsp1ds<. (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[μ]uC(N,p,s,ΛK)W2diam(Ω)s,p[μ+],a.e. inΩ. (1.10)

    In addition, for any q(0,N(p1)Ns) and h(0,s), there exists a positive constant c=c(N,p,s,ΛK,q,h,|Ω|) such that

    (Ω|g(u)|dx)1p1+(RNRN|u(x)u(y)|q|xy|N+hqdxdy)1qc(|μ|(Ω))1p1. (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(p1)Nsp, hence problem (1.12) admits a very weak solution in this case. In the supercritical case κN(p1)Nsp, 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, κ>p1 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[μ]uCW2diam(Ω)s,p[μ+],a.e. inΩ. (1.13)

    In addition, for any q(0,N(p1)Ns) and h(0,s), there exists a positive constant c=c(N,p,s,ΛK,q,h,|Ω|) such that

    (Ω|u|κdx)1p1+(RNRN|u(x)u(y)|q|xy|N+hqdxdy)1qc(|μ|(Ω))1p1. (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 gC(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|p1LNNspw(RN)t0. (1.17)

    Here, t0>0 depends on N,|Ω|,Λg,ΛK,a,s,p,d,ρ0. In addition, for any q(0,N(p1)Ns) 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

    (RNRN|v(x)v(y)|q|xy|N+hqdxdy)1qc(1+ρ|τ|(Ω))1p1. (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 κNN2s 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 p1<κ<N(p1)Nsp. In general, a sufficient condition for the solvability of (1.19) is the following.

    Proposition 1.5. Let s(0,1), 1<p<Ns, κ>p1 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

    M1Wd(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 κ>p1.

    Theorem 1.6. Let s(0,1), 1<p<Ns, κ>p1 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

    C11Wd(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 ERN.

    (iii) There exists a positive constant C3 such that

    B(W2diam(Ω)s,p[τB])κdxC3τ(B) (1.24)

    for any ball BRN, 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 p1<q<N(p1)Nsp then spqqp+1>N, this implies that Capsp,qqp+1({x})>0 for any xRN (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 ˜g0. 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):={RNRN|u(x)u(y)|q|xy|N+sqdxdy+RN|u|qdx<}, (2.1)

    endowed with the quasinorm

    uWs,q(RN):=(RNRN|u(x)u(y)|q|xy|N+sqdxdy)1q+(RN|u|qdx)1q.

    When q1, 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 C0(Ω) 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 uWs,p0(Ω) is a weak solution of (2.2), if it satisfies

    RNRN|u(x)u(y)|p2(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 uWs,p(RN) is a weak supersolution (resp. subsolution) of L in Ω, if and only if satisfies

    RNRN|u(x)u(y)|p2(u(x)u(y))(ϕ(x)ϕ(y))K(x,y)dxdy0(resp.0)

    for any nonnegative ϕWs,p0(Ω).

    Next we state the comparison principle.

    Proposition 2.3. ([23, Lemma 6]). Let uWs,p(RN) be a weak supersolution of L in Ω as well as let vWs,p(RN) be a weak subsolution of L in Ω such that (vu)+Ws,p0(Ω). Then, uv 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 DRN be a domain. Denote Lpw(D), 1p<, 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):=|{xD:|f(x)|>a}|cap,a>0. (2.3)

    The function λf is called the distribution function of f. For p1, denote

    Lpw(D)={f Borel measurable:supa>0apλf(a)<},
    fLpw(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

    fLpw(D)=sup{ω|f|dx|ω|1/p:ωD,ω measurable,0<|ω|<}. (2.5)

    More precisely,

    fLpw(D)fLpw(D)pp1fLpw(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}dxspupLpw(D). (2.7)

    Proposition 2.5. Let 1<p<Ns, μLp(Ω) and uWs,p0(Ω) be the unique weak solution of (2.2). Then there exists a positive constant C=C(p,s,N,ΛK) such that

    |u|p1LNNspw(RN)CΩ|μ|dx. (2.8)

    Proof.

    Let k>0. Taking Tk(u) as test function and using the fact that

    |u(x)u(y)|p2(u(x)u(y))(Tk(u)(x)Tk(u)(y))|Tk(u)(x)Tk(u)(y)|p,x,yRN,

    we obtain

    RNRN|Tk(u)(x)Tk(u)(y)|p|xy|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}|kNpNspRN|Tk(u)(x)|NpNspdxCkN(p1)Nsp(Ω|μ|dx)NNsp,

    which implies the desired result.

    Proposition 2.6. Let μLp(Ω) and uWs,p0(RN) be the unique weak solution of (2.2). Then there exists a positive constant C=C(p,s,N,ΛK) such that

    RNRN|u(x)u(y)|p(d+|u(x)|+|u(y)|)ξdxdy|xy|N+spCd1ξ(ξ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

    RNRN|u(x)u(y)|p2(u(x)u(y))(ϕ±(x)ϕ±(y))K(x,y)dxdy=Ωϕ±μdx.

    Now,

    (ϕ±(x)ϕ±(y))=±(ξ1)(u±(x)u±(y))10(d+tu±(y)+(1t)u±(x))ξdt,

    which implies

    |u(x)u(y)|p2(u(x)u(y))(ϕ±(x)ϕ±(y))K(x,y)(ξ1)|u(x)u(y)|p2(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(p1)Ns,p} μLp(Ω) and uWs,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

    (RNRN|u(x)u(y)|q|xy|N+hqdxdy)1qc(Ω|μ|dx)1p1. (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

    RNRN|u(x)u(y)|q|xy|N+hqdxdy=B2R(x0)B2R(x0)|u(x)u(y)|q|xy|N+hqdxdy+RNB2R(x0)RNB2R(x0)|u(x)u(y)|q|xy|N+hqdxdy+2RNB2R(x0)B2R(x0)|u(x)u(y)|q|xy|N+hqdxdy.

    Taking into account that u=0 a.e. in RNΩ, we can easily prove that

    RNB2R(x0)RNB2R(x0)|u(x)u(y)|q|xy|N+hqdxdy=0

    and

    RNB2R(x0)B2R(x0)|u(x)u(y)|q|xy|N+hqdxdy=RNB2R(x0)BR(x0)|u(x)|q|xy|N+hqdxdyBR(x0)|u(x)|qdxRNB2R(x0)1(1+|yx0|)N+hqdyBR(x0)|u(x)|qdx.

    Here, we have also used the fact that |xy|1+|yx0| for any (x,y)BR(x0)×(RNB2R(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|xy|N+hqdxdy=B3R2(x0)B3R2(x0)|u(x)u(y)|q|xy|N+hqdxdy+B2R(x0)B3R2(x0)B2R(x0)B3R2(x0)|u(x)u(y)|q|xy|N+hqdxdy+2B2R(x0)B3R2(x0)B3R2(x0)|u(x)u(y)|q|xy|N+hqdxdyB3R2(x0)B3R2(x0)|u(x)u(y)|q|xy|N+hqdxdy+BR(x0)|u|qdx.

    Combining all above, we have

    RNRN|u(x)u(y)|q|xy|N+hqdxdyB2R(x0)B2R(x0)|u(x)u(y)|q|xy|N+hqdxdy+Ω|u|qdxB2R(x0)B2R(x0)|u(x)u(y)|q|xy|N+hqdxdy. (2.12)

    Now, by Hölder inequality we obtain

    B2R(x0)B2R(x0)|u(x)u(y)|q|xy|N+hqdxdy=B2R(x0)B2R(x0)(|u(x)u(y)|p(d+|u(x)|+|u(y)|)ξ|xy|ps(d+|u(x)|+|u(y)|)ξ|xy|p(sh))qpdxdy|xy|N(B2R(x0)B2R(x0)|u(x)u(y)|p(d+|u(x)|+|u(y)|)ξ|xy|N+spdxdy)qp×(B2R(x0)B2R(x0)(d+|u(x)|+|u(y)|)ξqpq|xy|Nqp(sh)pqdxdy)pqp. (2.13)

    Setting

    d=(Ω|u(y)|ξqpqdx)pqξq

    and combining (2.10) and (2.13), we conclude

    B2R(x0)B2R(x0)|u(x)u(y)|q|xy|N+hqdxdycdqp(Ω|μ|dx)qp. (2.14)

    If p>2sN, 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<p2Ns, we have that 0<q1, 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(p1)Ns<p. Hence, we may choose ξ>1 such that 1<γ:=ξq(p1)(pq)<NNsp. Thus, by (2.6) and (2.8), we deduce

    (Ω|u|γ(p1))1γC(γ,N,p,s,|Ω|,ΛK)Ω|μ|dx,

    which in turn implies

    dC(γ,N,p,s,|Ω|,ΛK)(Ω|μ|dx)1p1.

    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|p1LNNspw(RN)C1(N,p,s,ΛK)μ(Ω) (2.15)

    and

    RNRN|Tk(u)(x)Tk(u)(y)|p|xy|N+spdxdykΛK|μ|(Ω),k>0. (2.16)

    In addition, for any q(0,N(p1)Ns) and h(0,s), there exist a positive constant C2=C2(N,p,s,ΛK,q,h,|Ω|) such that

    (RNRN|u(x)u(y)|q|xy|N+hqdxdy)1qC2|μ|(Ω)1p1. (2.17)

    Proof. Let {ρn}n be a sequence of mollifiers and μn=ρnμ. Then μnC0(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|p1LNNspw(RN)C1(N,p,s,ΛK)μ(Ω),nN, (2.18)
    RNRN|Tk(un)(x)Tk(un)(y)|p|xy|N+spdxdykΛKμ(Ω),k>0andnN, (2.19)

    and

    (RNRN|un(x)un(y)|q|xy|N+hqdxdy)1qC2(N,p,s,ΛK,q,h)μ(Ω)1p1 (2.20)

    for any nN, q(0,N(p1)Ns) 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:RNR satisfying the following properties:

    (ⅰ) uWh,q(RN) for any 0<q<N(p1)Ns and 0<h<s.

    (ⅱ) unu a.e. in RN, u=0 a.e. in RNΩ and uunWh,q(RN)0 for any q(0,N(p1)Ns) and h(0,s).

    (ⅲ) Tk(u)Ws,p0(Ω) for any k>0.

    Step 1. There exists a subsequence, still denoted by un, such that

    limnlimm|{xΩ:|unum|>η}|=0,η>0.

    Let n,mN and η,ρ>0. Then

    {|unum|>η}{|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,kk0. (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,˜jn0. (2.22)

    The desired result follows by (2.21) and (2.22).

    Step 2. Weak convergence of the truncates. Since unu 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|p1LNNspw(RN)C1(N,p,s,ΛK)μ(Ω), (2.23)
    RNRN|Tk(u)(x)Tk(u)(y)|p|xy|N+spdxdykΛKμ(Ω),k>0, (2.24)

    and

    (RNRN|u(x)u(y)|q|xy|N+hqdxdy)1qC2(N,p,s,ΛK,q,h)μ(Ω)1p1 (2.25)

    for any q(0,N(p1)Ns) and h(0,s).

    By (2.20), (2.25) and the fact that unu a.e. in RN, We can easily show that uunWh,q(RN)0 for any q(0,N(p1)Ns) and h(0,s). Let ϕC0(Ω), q(p1,N(p1)Ns) and h(max(sp1p1,0),s). For any bounded Borel set ERN, we have that

    |EE|un(x)un(y)|p2(un(x)un(y))(ϕ(x)ϕ(y))K(x,y)dxdy|C(ϕ,ΛK)EE|un(x)un(y)|p1|xy|N+hp1+sph(p1)1C(ϕ,ΛK)(EE|un(x)un(y)|q|xy|N+hqdxdy)p1q(EE|xy|Nq(sph(p1)1)qp+1dxdy)qp+1q.

    This, together with (ⅰ), (ⅱ) and the fact that q(sph(p1)1)qp+1<0, implies that

    Ωϕdμ=limnΩϕμndx=limnRNRN|un(x)un(y)|p2(un(x)un(y))(ϕ(x)ϕ(y))K(x,y)dxdy=RNRN|u(x)u(y)|p2(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

    C1Wd(x)8s,p[μ](x)u(x)C(essinfBd(x)4(x)u+Wd(x)2s,p[μ](x)+((d(x)4)spRNBd(x)4(x)u(y)p1|xy|N+spdy)1p1) (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

    (ⅰ) uWh,q(RN) for any 0<q<N(p1)Ns and 0<h<s.

    (ⅱ) unu a.e. in RN, u=0 a.e. in RNΩ and uunWh,q(RN)0 for any q(0,N(p1)Ns) and h(0,s).

    Since μnC0(RN) is a nonnegative function, by (2.3), we have that un0 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 μkM+(Ω) such that

    RNRN|uk(x)uk(y)|p2(uk(x)uk(y))(ϕ(x)ϕ(y))K(x,y)dxdy=Ωϕ(x)dμk (2.27)

    for any ϕC0(Ω). Since uku in RN, we have that uukWh,q(RN)0 for any h(0,s) and q(0,N(p1)Ns). This, together with (2.27), implies

    Ωϕ(x)dμkΩϕ(x)dμ,ϕC0(Ω). (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,

    C1Wd(x)8s,p[μk](x)uk(x),for a.e.xΩandk>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):=essliminfyxuk(y)=uk(x),for a.e.xRN.

    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:=limkvk 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[μ]uCW2diam(Ω)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 vnWs,p0(Ω) the solution of

    {Lvn=μn,inB2R(x0),vn=0,inRNB2R(x0).

    By Proposition 2.3, we have that vn0 and vnun, where unWs,p0(Ω) is the weak solution of (2.2) with μ=μn. By statements (ⅰ)–(ⅲ) in the proof of Proposition 2.8, there exist subsequences {unk,vnk}k=1 such that unku and vnkv a.e. in RN and

    uunkWh,q(RN)+vvnkWh,q(RN)0

    for any h(0,s) and q(0,N(p1)Ns). Combining all above, we may deduce that uv a.e. in RN and v is a nonnegative very weak solution to

    {Lv=μ+,inB2R(x0),v=0,inRNB2R(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)spRNBR2(x)|v(y)|p1|xy|N+spdy)1p1.

    By (2.15) and (2.6), we derive that

    essinfBR2(x)v(BR2(x)|v|p1dx)1p1RNspp1μ+(BR(x))1p1W2Rs,p[μ+](x), (2.31)

    and

    Tail(v;x0,R2)(B2R(x0)|v|p1dx)1p1W2Rs,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 RNB2R(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 gC(R) and rg(r)0. Let ΩRN be an open bounded domain and μ(Ws,p0(Ω)). Set G(r)=r0g(s)ds,

    J(v)=1pRNRN|v(x)v(y)|pK(x,y)dxdy+ΩG(v)dx<μ,v>

    and

    XG(Ω)={vWs,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

    RNRN|uμ(x)uμ(y)|p2(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 tR and the fractional Sobolev inequality, we can easily show the existence of a positive constant C=C(p,Ω,ΛK) such that

    vnpWs,p0(Ω)C(J(vn)+μp(Ws,p0(Ω))),nN. (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 vWs,p0(Ω) such that there hold:

    (ⅰ) vnv a.e. in RN.

    (ⅱ) vnv in Ws,p0(Ω) and vnv in Wh,q(RN) for any h(0,s) and q(1,p).

    By Fatou's lemma, we obtain

    1pRNRN|v(x)v(y)|pK(x,y)dxdy+ΩG(v)dxlim infk1pRNRN|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

    RNRN|vk(x)vk(y)|p2(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

    RNRN|vk(x)vk(y)|pK(x,y)dxdy+Ωg(vk)vkdx=<μ,vk>1pvkpWs,p0(Ω)+1pμp(Ws,p0(Ω)),

    which implies

    RNRN|vk(x)vk(y)|p|xy|N+spdxdy+Ωgk(vk)vkdxC(Λ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 vWs,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)|dx1λE{vk>λ}|ζ||vkgk(vk)|dxNλΩvkgk(vk)dxMNλ.

    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 uXG(Ω), we have that uXGk(Ω), Gk(u)G(u) and

    1pRNRN|vk(x)vk(y)|pK(x,y)dxdy+ΩGk(vk)dx<μ,vk>1pRNRN|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

    RNRN|uν(x)uν(y)|p2(uν(x)uν(y))(vk(x)vk(y))K(x,y)dxdyRNRN|uμ(x)uμ(y)|p2(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(Ω), gC(RN) be a nondecreasing function with g(0)=0 and uWs,p0(Ω) satisfy (3.1). Then there holds,

    Ω|g(u)|dxΩ|μ|dx. (3.6)

    In addition, if we assume that μ0, then u0 a.e. in RN.

    Proof. Let k>0. Using ϕk=tanh(ku) as test function in (3.1), we obtain

    RNRN|u(x)u(y)|p2(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)=(1tanh2(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)ϕk0 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 gC(R) is nondecreasing such that g(0)=0.

    Lemma 3.3. Let gL(R) and λiM+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

    {Lu2g(u2)=λ2,inΩ,u=0,inRNΩ, (3.9)

    such that there hold

    u1,u20andu2uu1,a.e. inRN. (3.10)

    In addition, for any q(0,N(p1)Ns) and h(0,s), there exists a positive constant c=c(N,p,s,ΛK,q,h,|Ω|) such that

    (Ω|g(u)|dx)1p1+(RNRN|u(x)u(y)|q|xy|N+hqdxdy)1qc(λ1(Ω)+λ2(Ω))1p1 (3.11)

    and

    (Ω|g((1)i+1ui)|dx)1p1+(RNRN|ui(x)ui(y)|q|xy|N+hqdxdy)1qcλi(Ω)1p1. (3.12)

    Finally, there exist very weak solutions vi to (2.2) with μ=λi (i = 1, 2) such that

    0uiviCiW2diam(Ω)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,iC0(RN). By Proposition 3.1, there exist unique solutions un,un,i,vn,iWs,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,2g(un,2)=λn,2,inΩ,u=0,inRNΩ,

    and

    {Lvn,i=λn,iinΩu=0,inRNΩ,

    such that there holds

    vn,2un,2unun,1vn,1,a.e. inRN. (3.14)

    By Lemma 3.2 and Proposition 2.7, for any q(0,N(p1)Ns) and h(0,s), there exists a positive constant c=c(N,p,s,ΛK,q,h,|Ω|) such that

    (Ω|g(un)|dx)1p1+(RNRN|un(x)un(y)|q|xy|N+hqdxdy)1qc(Ωλn,1+λn,2dx)1p1, (3.15)
    (Ω|g(un,1)|dx)1p1+(RNRN|un,1(x)un,1(y)|q|xy|N+hqdxdy)1qc(Ωλn,1dx)1p1, (3.16)
    (Ω|g(un,2)|dx)1p1+(RNRN|un,2(x)un,2(y)|q|xy|N+hqdxdy)1qc(Ωλn,2dx)1p1 (3.17)

    and

    (RNRN|vn,i(x)vn,i(y)|q|xy|N+hqdxdy)1qc(Ωλn,idx)1p1. (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 unu, un,iui vn,ivi 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 λiM+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 nN. By Lemma 3.3, there exist very weak solutions un,un,i,vn,iWs,p0(Ω) of the following problems

    {Lun+Tnog(un)=λ1λ2inΩu=0inRNΩ,
    {Lun,1+Tnog(un,1)=λ1,inΩ,u=0,inRNΩ,
    {Lun,2Tnog(un,2)=λ2,inΩ,u=0,inRNΩ

    and

    {Lvi=λi,inΩ,u=0,inRNΩ,

    such that there holds

    CW2diam(Ω)1,p[λ2]v2un,2unun,1v1CW2diam(Ω)1,p[λ1],a.e. inRN

    and for any nN. 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˜q1(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,s1. (3.20)

    Then for any s01, there hold

    g(|v|)L1(Ω)ΩEs0(v)g(|v|)dx+˜qC0s0s˜q1g(s)ds,g(|v|)L1(Ω)ΩEs0(v)g(|v|)dx˜qC0s0s˜q1g(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

    v2unv1,a.e. inRNandnN.

    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

    vp11LNNspw(RN)+vp12LNNspw(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(p1)Nsps0s˜q1(g(s)g(s)ds,nN,

    where ˜q=N(p1)Nsp. 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α(ξ)=F1(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:gLβ(RN)},

    with norm

    fLα,β(RN):=gLβ(RN)=GαfLβ(RN).

    The Bessel capacity is defined as follows.

    Definition 3.6. Let α>0, 1<β< and ERN. Set

    SE:={gLβ(RN):g0,Gαg(x)1for anyxE}.

    Then

    Capα,β(E):=inf{gβLβ(RN);gSE}. (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}nLsp,κp1(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]κp1Lκp1(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,2un,2unun,1vn,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,iun+1,iandvn,ivn+1,i,a.e. inRN,nNandi=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μ(Ω),nN.

    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 unu, un,iui vn,iv 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, gC(R) is a nondecreasing function and

    |g(t)|a|t|dfor some a>0,d>p1 and for any |t|1. (4.2)

    Let us state the first existence result.

    Lemma 4.1. Let 1<p<Ns and τC0(RN) be such that τL1(RN)1. Assume that gL(Ω)C(R) satisfies (3.19) for

    ˜q=N(p1)Nsp.

    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 vWs,p0(Ω) satisfying

    |v|p1LNNspw(Ω)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{NNsp,dp1} and vL1(Ω). Since gL(Ω), we can easily show that the following problem

    {Lu=g(|v|1p1sign(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)|p1sign(T(v)),vL1(Ω). (4.5)

    By (2.8), we obtain

    S(v)LNNspw(Ω)C(s,p,N,ΛK)(ρΩ|τ|dx+Ω|g(|v|1p1sign(v))|dx)C(s,p,N,ΛK)(ρ+Ωg(|v|1p1)g(|v|1p1)dx). (4.6)

    Let vLNNspw(Ω). For any λ>0, we set Eλ:={xΩ:|v(x)|1p1>λ} and e(λ)=Eλdx. By (2.4) and (2.6), we can easily show that

    e(λ)C(N,s,p)vNNspLNNspw(Ω)λN(p1)Nsp.

    By the above inequality and Lemma 3.5 with λ0=1 and ˜q=N(p1)Nsp, we deduce

    Ωg(|v|1p1)g(|v|1p1)dx2aΩ|v|κdx+C(p,s,N)vNNspLNNspw(Ω)Λg.

    Let λ=vLNNspw(Ω). By (2.6), we have that

    Ω|v|κd=0|{xΩ:|v|t}|dtκ=λ0|{xΩ:|v|t}|dtκ+λ|{xΩ:|v|t}|dtκ|Ω|λκ+κλNNspλtκNNsp1dtC(Ω,κ,s,p,N)λk.

    Combining all above, we may prove that

    S(v)LNNspw(Ω)C(p,N,κ,|Ω|,Λg,ΛK,a)(ρ+vNNspLNNspw(Ω)+vκLNNspw(Ω)).

    Therefore, if vLNNspw(Ω)t then

    S(v)LNNspw(Ω)C(tNNsp+tκ+ρ). (4.7)

    Since 1<κ<NNsp, 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(tNNsp+tκ+ρ)t0,

    where C is the constant in (4.7). Hence,

    vLNNspw(Ω)t0S(v)LNNspw(Ω)t0. (4.8)

    Next, we apply Schauder fixed point theorem to our setting.

    We claim that S is continuous. First we assume that vnv in L1(Ω) and T(vn)T(v) in W1,p0(Ω), then by fractional Sobolev inequality, we have

    Ω|T(vn)T(v)|dx|Ω|pNN+spNpT(vn)T(v)LNpNsp(Ω)C|Ω|pNN+spNpT(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 gL(R), we have that S(vn)Lβ(Ω) and {S(vn)} is uniformly bounded in Lβ(Ω) for any β(1,NNsp). Hence, there exists k0N, such that

    Ω({xΩ:|S(vn)(x)|k}{xΩ:|S(v)(x)|k})|S(vn)S(v)|dxε3kk0,andnN. (4.11)

    Now, we set

    Ak0,n={xΩ:|T(vn)(x)|k1p10}{xΩ:|T(v)(x)|k1p10}

    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,nBδ,n||T(vn)|p1sign(T(vn))|T(v)|p1sign(T(v))|dx+Ak0,nBδ,n||T(vn)|p1sign(T(vn))|T(v)|p1sign(T(v))|dx. (4.12)

    Since h(t)=tp1sign(t) is uniformly continuous in [k0,k0], there exists δ0>0 independent of n such that

    Ak0,nBδ0,n||T(vn)|p1sign(T(vn))|T(v)|p1sign(T(v))|dxε3. (4.13)

    Moreover, by (4.9), there exists n0=n0(δ0,k0,p)N such that

    Ak0,n0Bδ0,n0||T(vn0)|p1sign(T(vn0))|T(v)|p1sign(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:=suptR|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

    RNRN|T(vn)(x)T(vn)(y)|pK(x,y)dxdy=ΩT(vn)(g(|vn|1p1sign(vn))dx+ΩT(vn)τdxM|Ω|p1p(Ω|T(vn)|pdx)1p+(Ω|T(vn)|pdx)1p(Ω|τ|pp1dx)p1pC1(M,Ω,p,N,τ,s)(RNRN|T(vn)(x)T(vn)(y)|p|xy|N+spdxdy)1p. (4.15)

    Therefore,

    RNRN|T(vn)(x)T(vn)(y)|p|xy|N+spdxdyCp1p1(M,Ω,p,N,τ,s,ΛK). (4.16)

    Using ϕ=T(vn)T(v) as test function, we have

    I:=RNRN|T(vn)(x)T(vn)(y)|p2(T(vn)(x)T(vn)(y))(ϕ(x)ϕ(y))K(x,y)dxdyRNRN|T(v)(x)T(v)(y)|p2(T(v)(x)T(v)(y))(ϕ(x)ϕ(y))K(x,y)dxdy=Ωϕ(g(|vn|1p1sign(vn))g(|v|1p1sign(v)))dx=:II. (4.17)

    We first treat I. On one hand, since

    (|a|p2a|b|p2b)(ab)C(p)|ab|2(|a|+|b|)2p

    for any (a,b)R2N{(0,0)} and p(1,2), we have

    IC(p)RNRN|ϕ(x)ϕ(y)|2(|T(vn)(x)T(vn)(y)|+|T(v)(x)T(v)(y)|)p2K(x,y)dxdy. (4.18)

    On the other hand, by Hölder inequality, we obtain

    RNRN|ϕ(x)ϕ(y)|p|xy|N+spdxdyΛKRNRN|ϕ(x)ϕ(y)|pK(x,y)dxdyC(p,ΛK)(RNRN(|T(vn)(x)T(vn)(y)|+|T(v)(x)T(v)(y)|)pK(x,y)dxdy)2p2Ip2C(p,C1,Ω,ΛK)Ip2, (4.19)

    where C1 is the constant in (4.16). Hence, by (4.18) and (4.19), we obtain

    C(RNRN|ϕ(x)ϕ(y)|p|xy|N+spdxdy)2pI. (4.20)

    Next we treat II. Let r=NpNsp, proceeding as in the proof of (4.15), we have

    II(Ω|ϕ|rd)1r(Ω|g(|vn|1p1sign(vn))g(|v|1p1sign(v))|rd)1rC(N,p,s)(RNRN|ϕ(x)ϕ(y)|p|xy|N+spdxdy)1p×(Ω|g(|vn|1p1sign(vn))g(|v|1p1sign(v))|rd)1r, (4.21)

    where in the last inequality we used the fractional Sobolev inequality.

    Combining (4.17), (4.20) and (4.21), we obtain

    (RNRN|ϕ(x)ϕ(y)|p|xy|N+spdxdy)1pC(p,C1,Ω,s,ΛK)(Ω|g(|vn|1p1sign(vn))g(|v|1p1sign(v))|rd)1r. (4.22)

    Since g(||p1sign()) is uniformly continuous in R, bounded and vnv in L1(Ω), we obtain

    limnΩ|g(|v|1p1sign(v))g(|vn|1p1sign(vn))|rdx=0,

    which, together with (4.22), implies the desired result.

    Case 2. p2. We note here that

    (|a|p2a|b|p2b)(ab)C(p)|ab|p

    for any (a,b)R2N and p2. Thus,

    IC(N,p,ΛK)(RNRN|ϕ(x)ϕ(y)|p|xy|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)|p1sign(T(vn))|ψ|p1sign(ψ) in L1(Ω).

    Now set

    O:={vL1(Ω):vLNNspw(Ω)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 vO such that S(v)=v. This means that u=v1p1sign(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 unWs,p0(Ω) of

    {Lv=gn(v)+ρτn,inΩ,v=0,inRNΩ.

    In addition, it satisfies

    |un|p1LNNspw(Ω)t0, (4.24)

    where t0>0 depends on N,Ω,Λg,ΛK,a,s,p,d.

    By (4.24), we have that

    |{xΩ:|un|>s}|tNNsp0sN(p1)Nsp.

    Hence by Proposition 3.5,

    Ω|gn(un)|dxC,nN,

    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(p1,N(p1)Ns) and h(0,s), there exists a positive constant c=c(N,s,p,ΛK,s,h,q,|Ω|) such that

    (RNRN|un(x)un(y)|q|xy|N+hqdxdy)1qc(C+ρΩ|τn|dx)1p1. (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 unu 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))s1(p1)NNspds,s01.

    Let ε>0, then there exists s0 such that

    C(t0,p,Λg,N)s0(g(s)g(s))s1(p1)NNspdsε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[ν]21p1C(W2diam(Ω)s,p[wκ]+W2diam(Ω)s,p[ρdτ])21p1C((AC)κp1ρκ(p1)2MW2diam(Ω)s,p[τ]+W2diam(Ω)s,p[ρdτ])21p1C((AC)κp1Mρκp+1(p1)2+1)W2diam(Ω)s,p[ρdτ].

    If we choose A=21p1+1 and ρ small enough such that (AC)κp1Mρκp+1(p1)2+1<2, we deduce that

    CW2diam(Ω)s,p[ν]w. (4.26)

    Now, let x0Ω be such that W2diam(Ω)s,p[τ](x0)<. If 0vc0W2diam(Ω)s,p[τ] a.e. in RN, for some constant c0>0, then we have

    (Ω|v|κdx)1p1(Bdiam(Ω)(x0)|v|κdx)1p1C(Ω,N,s,p,M,K,c0)W2diam(Ω)s,p[τ](x0)<.

    Thus vLκ(Ω).

    Let u00 be a very weak solution of

    {Lu0=ρτ,inΩ,v=0,inRNΩ,

    satisfying C1Wd(x)8s,p[μn1](x)u0(x)CW2diam(Ω)s,p[ρτ](x) a.e. in Ω. We may construct a nondecreasing sequence {un}n0, such that un is a very weak solution to problem

    {Lun=uκn1+ρτ,inΩ,v=0,inRNΩ

    and satisfies

    C1Wd(x)8s,p[μn1](x)un(x)CW2diam(Ω)s,p[μn1](x),for a.e.xΩ,

    for any nN, where dμn1=uκn1dx+ρdτ. In addition, by (4.26) and the above inequality, there holds

    C1Wd(x)8s,p[μn1](x)un(x)w(x),for a.e.xΩ, (4.27)

    where the positive constant C1 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 unu 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 0gLκp1(RN;μ). We set

    Mμg(x):=supr>0,μ(B(x,r))0μ(B(x,r))1B(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))κp1dμc1RN|g(x)|κp1dμ (4.28)

    (see, e.g., [14]). Also,

    Ω(Wd(x)8s,p[gμ](x))κdxΩ(Wd(x)8s,p[μ](x))κ(Mμg(x))κp1dxCκΩuκ(x)(Mμg(x))κp1dxCκΩ(Mμg(x))κp1dμc2RN|g(x)|κp1dμ. (4.29)

    Let K=suppτ. By the assumption, we have that r0:=dist(K,Ω)>0. Set g=fanxiexianmyfhK˜g, for any nonnegative ˜gLκp1(RN;μK). We first note that Br08(x)K= if xΩ with d(x)<r08 or if xRNΩ, which implies

    Wd(x)24s,p[˜gμK](x)=0,

    if xΩ with d(x)<r024 or if xRNΩ. Therefore, by the above equality and (4.29), we have

    RN(Wr024s,p[˜gμK](x))κdxΩ(Wd(x)8s,p[gμ](x))κdxc2RN|˜g(x)|κp1dμK.

    Also, by [5, Theorem 2.3] (see also [1, Corollary 3.6.3]), we have

    RN(Wr024s,p[˜gμK](x))κdxRN(Gsp[˜gμK])κp1dx, (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])κp1dxc3RN|˜g(x)|κp1dμK. (4.31)

    Let fLκκp+1(RN). Then, for any ˜gLκp1(RN;μK), there holds

    |RNf(x)Gsp(˜gμK)(x)dx|=|RN˜g(y)Gspf(y)dμK|C1fLκκp+1(RN)˜gLκp1(RN;μK).

    The last inequality implies,

    RN|Gspf(x)|κκp+1dμKc4RN|f|κκp+1dxfLκκ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 FRN. (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
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1331) PDF downloads(194) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog