Поточечные оценки и критерий существования решений сублинейных эллиптических уравнений Текст научной статьи по специальности «Математика»

www.volsu.ru

DOI: https://doi.org/10.15688/mpcm.jvolsu.2017.3.2

UDC 517 LBC 22.161

POINTWISE ESTIMATES OF SOLUTIONS AND EXISTENCE CRITERIA FOR SUBLINEAR ELLIPTIC EQUATIONS

Igor E. Verbitsky

Professor, Department of Mathematics, University of Missouri verbitskyi@missouri.edu Columbia, MO 65211, USA

Abstract. We give a survey of recent results on positive solutions to sublinear elliptic equations of the type —Lu + V uq = f, where L is an elliptic operator in divergence form, 0 < q < 1, f > 0 and V is a function that may change sign, in a domain Q C Rra, or in a weighted Riemannian manifold, with a positive Green's function G. We discuss the existence, as well as global lower and upper pointwise estimates of classical and weak solutions u, and conditions that ensure u e Lr(Q) or u e W 1,P(Q).

Some of these results are applicable to homogeneous sublinear integral equations u = G(uqda) in Q, where 0 < q < 1, and a = — V is a positive locally finite Borel measure in Q. Here G(fda)(x) = fQ G(x,y), f (y) da(y) is an integral operator with positive (quasi) symmetric kernel G on Q x Q which satisfies the weak maximum principle. This includes positive solutions, possibly singular, to sublinear equations involving the fractional Laplacian,

(—A)au = auq, u > 0 in Q,

where 0 < q < 1, 0 < a < n and u = 0 in Qc and at infinity in domains Q C Rra with positive Green's function G.

Key words: sublinear elliptic equations, Green's function, weak maximum principle, fractional Laplacian.

t-

1. Introduction

HJ

Let Q be a connected open set in Rra, or a weighted Riemannian manifold (M,m) (see jn Sec. 2 below). We give a survey of recent results on positive solutions u to the sublinear £ elliptic differential equation

© —Lu + V (x) uq = f, x e Q, (1.1)

where 0 < q < 1, and

Lu = div(^Vw) (1.2)

is an elliptic operator with symmetric coefficients A = (a^-(x))™j=1 in divergence form, V and f are functions in Q, and f > 0. Here V (x) may change sign, and no extra boundary conditions are imposed on V.

One of the main goals is to obtain global pointwise estimates of u in terms of the function h which is defined as the minimal positive solution in Q of the equation —Lh = f, where f = 0. We also consider positive solutions u to the corresponding homogeneous equation

—Lu + V(x) uq = 0, (1.3)

with f = 0. In Sec. 2, we present lower pointwise estimates of solutions to (1.1) and (1.3) obtained in [16; 17].

There are also analogues of these estimates for positive (super) solutions u to sublinear integral equations

u + G(Vuq ) = F, (1.4)

which involve the integral operator Gf (x) = fn G(x,y) f (y) dy with a general positive kernel G: Q x Q ^ (0, If G = Gn is Green's function associated with the elliptic

operator L and F = Gf, then (1.4) is in a sense equivalent to (1.1).

In many cases it is enough to assume that the kernel G is lower semicontinuous, symmetric (or quasi-symmetric) and satisfies the weak maximum principle (see [16]), which includes many non-local operators, for instance, the fractional Laplacian (—A)f. This leads to pointwise estimates of solutions for (1.4) in the case where V = — a > 0, and a is a locally integrable function, or locally finite positive Borel measure in Q. In particular, we obtain sharp pointwise estimates of solutions (possibly singular) and existence results for the fractional Laplacian equation

(—A)au = auq, u > 0 in Q, u = 0 in Qc, (1.5)

in general domains Q C Mra with positive Green's function G for 0 < a < 2, or nice domains (the entire space Mra, or balls, or half-spaces in Rra) for 0 < a < n.

For entire solutions on Q = Rra and V = —a > 0, where a is a locally finite positive Borel measure in Rra, there are more complicated matching upper and lower bounds given in terms of certain nonlinear potentials of Th. Wolff type [8-10]. Earlier pointwise estimates for bounded positive solutions are due to Brezis and Kamin [6]. We consider such results in Sec. 3.

In the case V = —a > 0, there are necessary and sufficient conditions for the existence of positive solutions u e Lq(Q, a) and u e Lqoc(Q, a) to (1.4). The latter class includes all positive solutions 0 < u < da-a.e., and can be characterized either using certain localization techniques developed in [8] in the case L = (—A)f and Q = Rra, or integral inequalities with weights as in [23] in the general case. These results are discussed in Sec. 3 and Sec. 4.

We also present in Sec. 5 criteria for the existence of positive solutions to (1.4) of higher integrability u e Lq+y(Q, a) (y > 0), which are easier to characterize and are useful in applications (see [27]). In particular, the case y =1, i.e., solutions u e Ll+q(Q, a), corresponds to finite energy solutions to sublinear elliptic equations (1.3) in the Dirichlet space u e Wli2(Q) for uniformly elliptic operators L in divergence form.

Similar criteria for y > 0 can be applied to deduce conditions on a which ensure the existence of positive solutions u to (1.3) in Sobolev spaces W 1iP(Q) and Lebesgue spaces Lr(Q,dx) for general domains Q C Rra with positive Green's functions (see [24; 25]). Some earlier existence results of this type for a e Ls(Q,dx) in bounded domains Q are due to Boccardo and Orsina [4; 5]. Such existence theorems under more general assumptions on a which include measures (possibly singular with respect to Lebesgue measure) are discussed in Sec. 6.

2. Lower pointwise estimates of supersolutions

In this section, we assume that Q C Rra is a connected open set, and consider the divergence form operator

n

Lu = div(AVu) (x) = ^ (aij (x) dXju) , (2.1)

i,j=1

where aij are C2 functions in Q, and aij = a^. Assume that L is elliptic, that is, the matrix A = (aij (x)) is positive definite for any x (the uniform ellipticity is not required).

Let us assume that the Dirichlet Green function of L in Q exists. We denote it by GQ (x,y). We set

h (x) = Gf (x)=[ GQ (x,y) f (y) dy,

JQ

where f > 0 (/ = 0), and assume that h (x) < ro for all x e Q (note also that h (x) > 0 in Q), and that the integral

G(hV)(x)= f GQ (x,y) h (y) V (y) dy (2.2)

JQ

is well-defined.

An interesting special case is when f = 1, and h is the solution of the Dirichlet problem:

—Lh = 1 in Q,

!

ь n ЯП (2.3)

h = 0 on о a.

In other words, h(x) = Ex[tq], where tq = inf{t: Xt e Q} is the first exit time from Q of the (rescaled) Brownian motion Xt, and x e Q is a starting point. For bounded C1,1 domains Q, it is well known that h(x) ~ dQ(x), where

dn(x) = dist(x,d Q). (2.4)

One of the main results of [17] in the case 0 < q < 1 is the following global lower estimate for a (super) solution u> 0 to (1.1).

Theorem 2.1. Let Q C Rra be a domain with a positive Green function GQ. Suppose V, f e C(Q), and f > 0, f = 0 in Q. Let u e C2(Q) satisfy the inequality

—Lu + Vuq > f, u > 0 in Q. (2.5)

Set h = Gf and assume that h < ro in Q. Assume also that GQ(hqV)(x) is well-defined for all x e Q. Then the following estimate holds for all x e Q :

u(x) > h(x)

1 - a - 4)

h(x)

i

1-?

. (2.6)

+

This theorem is actually proved in [17] in a more general context of elliptic operators on weighted manifolds (M,m). We recall that if M is a smooth Riemannian manifold with the Riemannian metric tensor g = (g^), then the associated Laplace — Beltrami operator C0 acts on C2 functions u on M and is given in any chart x\, ...,xn by

1 n

C0u = /T-t- ^ (Vd-t gдгзdx.uj .

Here detg is the determinant of the matrix g = (g^), and (gtj) is the inverse matrix of (gij). The Riemannian measure m0 is given in the same chart by

dm0 = \/det g dx\...dxn, so that C0 is symmetric with respect to m0. Using the gradient operator V defined by

(V) = £ 9гзd:

3 = 1

and the divergence div on vector fields F%

1 n

div^=^ £F')

i=\

we have

C0 = div o V.

Now let u be a smooth positive function on M and consider the measure m on M given by

dm = udm0.

The couple (M,m) is called a weighted manifold, and u in this context is called a weight (see [14]). The following operator C

1 1 n

Cu := -div (uVu) = —V dXi (uy/det gg11 dx u) (2.7)

u det g V J )

acting on C2 functions u on M, is called the (weighted) Laplace operator of (M,m). This operator C, which is symmetric with respect to measure m, can be used in place of the operator L in Q C Rn above.

There is also a similar lower estimate for positive (super) solutions to the homogeneous equation (1.3).

Theorem 2.2. Let Q C Rra be a domain with a positive Green function GQ. Suppose V e C(Q), and u e C2(Q) satisfies the inequality

—Lu + Vug > 0, u > 0 in Q, (2.8)

where 0 < q < 1. Assume also that GV(x) is well-defined for all x e Q. Then the following estimate holds for all x e Q :

i

u(x) > [—(1 — q) GV(x)]+-q . (2.9)

Similar global pointwise estimates are obtained in [17] for all exponents q e R (q = 0), but we do not consider the case q e (0,1) here. Also, the preceding estimate remains true for nonnegative solutions u if the integration in the Green's potentials GQ(hqV) and GQV respectively is taken over the set Q+ = {x e Q: u > 0} in place of Q. Note that the differential inequalities (2.5) and (2.8) in Theorem 2.1 and Theorem 2.2 respectively are understood in the classical sense.

In the case V = —a < 0, where a is a locally finite measure in Q, the preceding estimate yields

_ i

u(x) > (1 — q) i GQ(x,y)da(y) , (2.10)

JQ

which holds for all non-trivial weak solutions u (possibly singular).

This estimate was obtained for the Laplacian L = A in the case Q = Rra in [6], without specifying the sharp constant (1 — q)T=9, under some additional assumptions on a > 0 (see also [9]).

There are analogues of this estimate for integral operators with general positive kernels G: Q x Q ^ (0, +<x>] on a locally compact Hausdorff space Q that are lower semicontinuous, symmetric, and satisfy the weak maximum principle with constant b > 1 (see [16]). For a locally finite Radon measure v e ^+(Q), we denote by Gv the potential

Gy(x)= G(x,y) d\(y).

Q

We recall that the kernel G is said to satisfy the weak maximum principle (with constant b > 1) provided that

Gv(x) < M, Vx e supp v Gv(x) < b M, Vx e Q, (2.11)

for any constant M > 0 and v e ^+(Q). When b = 1, we say that G satisfies the strong maximum principle.

These assumptions hold for many non-local operators L of interest in nonlinear analysis, in particular, the fractional Laplacian operators (—A) f in domains Q C Rra. For such kernels G and a e M+(Q), let

G(fda)(x)= G(x,y) f(y) da(y), x e Q,

Q

where f > 0 is a measurable function.

If u> 0 da-a.e. is a supersolution, i.e.,

G(uqda)(x) < u(x) < ro da—a.e., (2.12)

then

i

u(x) > (1 - q)^b-i—\Ga(x)]1-9 da—a.e. in Q. (2.13)

The constant (1 — q)1— in (2.13) in the case b = 1 coincides with that in (2.10). Similar results hold for quasi-symmetric kernels G such that there exists a constant a > 0 for which

G(x, y) < aG(x, y), x,y e Q, (2.14)

with a constant on the right-hand side of (2.13) depending on both a and b. See [2] where quasi-symmetry is established for many elliptic operators in non-divergence form. Properties (2.11) and (2.14) will be used extensively in Sections 4 and 6 below.

3. Matching upper and lower pointwise estimates of solutions

In this section, we discuss recent results established in [8] for entire solutions to quasilinear elliptic equations of the type

—Apu = auq in Rn,

(3.1)

liminf u(x) = 0, u > 0,

where 0 < q < p — 1. Here Ap = div(Vu|Vu|p-2) is the p-Laplace operator, and a > 0 is an arbitrary locally integrable function, or locally finite Borel measure, a e ^+(Rn); if a e LjOc(Rn) we write da = adx.

The main goal is to give necessary and sufficient conditions on a for the existence of weak solutions to (3.1), and obtain matching upper and lower global pointwise bounds for positive solutions. In this approach, key integral inequalities are identified, and new nonlinear potentials of Wolff type are constructed that are intrinsic to a number of similar problems.

This approach is applicable to general quasilinear ^.-Laplace operators of divergence type div ^.(x, Vu) under standard boundedness and structural assumptions on A, as well as to fully nonlinear fc-Hessian operators, and the fractional Laplacian equations (see [8]).

We will restrict ourselves here to the case p = 2 and linear operators div(^.Vu) under the uniform ellipticity assumption on the matrix A e L™(Rn)nxn. In this setup, the fundamental solution of the ^.-Laplacian is known to be pointwise equivalent to that of the Laplacian, which reduces the problem to estimates involving Newtonian potentials. In fact, we will consider below the fractional Laplacian equations,

' (—A)fu = auq in Rn,

(3.2)

lim inf u( x) = 0, u > 0,

for 0 < q < 1 and 0 < a < n; this includes the range a > 2 where the usual maximum principle is not available.

In the classical case p = 2, equation (3.1), or equivalently (3.2) with a =2, provides a model sublinear elliptic problem if 0 < q < 1. It is easy to see that it is equivalent to the integral equation u = N(uqda), where Nu = (—A)-1^ is the Newtonian potential of du = uqda on Rn.

As we emphasize in [8], such problems are governed by the important integral inequality i

(i da) ' < k ||Ap||Li(Rn), (3.3)

\jRn J

for all test functions p e C2(Rra) vanishing at infinity such that —Ap > 0. This inequality turns out to be necessary and sufficient for the existence of a positive solution u e Lq(Rra, a) to (3.1) ifp = 2 and 0 <q< 1.

Inequality (3.3) represents the end-point case of the well-studied (Lp,Lq) trace inequalities for p > 1. A comprehensive treatment of trace inequalities in that case in terms of Sobolev capacities can be found in [22].

To characterize the existence of arbitrary weak solutions u e Lqoc(Rra, a), and obtain sharp global pointwise estimates, a localization procedure is used in [8] which involves the so-called intrinsic Wolff potentials.

More precisely, one needs to consider a localized version of (3.3) where the measure a is replaced with as = Xb a restricted to a ball B = B(x, r), and the corresponding best constant k is denoted by k(B). These constants are used as building blocks in the key tool, a potential of Th. Wolff type (see [1] for the theory and history of Wolff's potentials),

K<x)= r MBi^Z*, , 6 r. (3.4)

J0 rn-2 r

This nonlinear potential, together with the usual Newtonian potential N a, provides sharp bilateral estimates of solutions in the case p = 2 and 0 < q < 1.

This work has been motivated by the results of Brezis and Kamin [6] who proved that (3.1), with p =2 and 0 < q < 1, has a bounded solution u on Rra if and only if Na e L^(Rn); moreover, such a solution is unique, and there exists a constant c > 0 so that

c-1 (na(x)) 1 9 < u(x) < cNa(x), x e Rra (3.5)

The lower bound in (3.5), which is generalized in (2.10) and (2.13), holds for all solutions. On the other hand, the upper bound is valid only for bounded solutions. Some generalizations of the upper bound to possibly singular solutions in the form

u(x) < с

N a(x) + ^N a(x)j

1

1-4

xe

under additional assumptions on a, can be found in [9].

As was pointed out in [6], both the lower and upper estimates in (3.5) are sharp in a sense. However, there is a substantial gap between them. In Theorem 3.1 below this gap is bridged by using Ka + (Na)1—. The corresponding pointwise estimates hold for all, possibly singular solutions, as well as more general nonlinear equations.

In the context of fractional Laplacian equations (3.2), we denote by Ia the Riesz potential of order a e (0, n) on Rra:

= lxla-n -k fx = c(a,'n) (—A)-a fx,

for f e M+(Rra). In what follows, the normalization constant c(a,n) will be dropped for the sake of convenience; in particular, I2f = Nf stands for the Newtonian potential.

We now introduce some elements of nonlinear potential theory which is intrinsic to (3.2). Let B = B (x, r) = {ye R™ : |x - y| < r} be a ball in R™, and let aB = Xb a stand for the measure a restricted to B. By k(B) we denote the least constant in the localized integral inequality

||IaV\\Lq(daB) < k(B) v(R™), Vv e .M + (Rn), (3.6)

where 0 < q < 1 and 0 < a < n.

We observe that the constant k(B) does not change if we restrict ourselves to absolutely continuous v e L+(Rra) on the right-hand side of (3.6).

We define the corresponding nonlinear potential of Wolff type by

W(*)=/" ^ * e R" (3.7)

J 0 r r

The following condition ensures that both K^,qa and Iaa are not identically infinite:

g

l™ [k(B(0, r))j 1-gg dr + l™ a(B(0, r)) dr < ^ (3

j y» a y» j a y»

Clearly, the first part of (3.8) is equivalent to Kœ,qa ^ whereas the second one is equivalent to Iaa ^ (and Iaa < a.e.)

Note that solutions u G Lqoc(Rra,da) to (3.2) are understood in the weak sense. Theorem 3.1. Let 0 < q < 1, 0 < a < n, and a G ^+(Rra).

(i) If (3.8) holds, then there exists a minimal solution u > 0 to (3.2) such that

c~l

i i Ka,ga + (Iaa) 1-* <u <c Ka,qa + (Iaa) ^ , (3.9)

where > 0 is a constant which depends only on , a, and n.

(ii) Conversely, if there exists a nontrivial (super) solution v to (3.2), then (3.8) holds, and v >u, where u is the minimal solution of statement (i).

We observe that neither of the two parts of the condition (3.8) implies the other one. Condition Iaa < <x> a.e. alone is not enough to ensure the existence of a global solution u even if all the local embedding constants k(B) in (3.6) are finite, unless a is radially symmetric (see [9]).

4. Existence of solutions for sublinear integral equations

In this section, G : Q x Q ^ (0, is a lower semicontinuous, quasi-symmetric kernel on Q x Q that satisfies (2.14) along with the weak maximum principle (2.11), on a locally compact Hausdorff space Q. Let us fix a locally finite Radon measure a G ^+(Q).

In [23; 25; 27], necessary and sufficient conditions are obtained for the existence of positive (super) solutions to the sublinear integral equation

u = G(uqda), 0 <u< da - a.e. in Q, (4.1)

The existence of solutions in Lebesgue spaces, u G LY(Q, a), turned out to be related to weighted norm inequalities of the type

||G( f da)\\Lr (q, da) <C \\f\\Ls(n,da), Vf G LS(Q, da), (4.2)

for integral operators G in the case 0 < r < s < s > 1.

The most important and most difficult case is when r = q (0 < r < 1) and s = 1, i.e.,

||G( f da)HLq {n,da) <C ||/||Li {n4u), Vf e L1(Q, da). (4.3)

In the special case of Green's kernels and the corresponding elliptic operators L, (4.3) boils down to a non-standard trace inequality

Ml*(n,da) <C limLi{n4x), V$ e C2(Q), > 0. (4.4)

It will be more convenient, as it often happens in the end-point case = 1, to replace (4.3) with a similar inequality for measures v e M+(Q) in place of L1(Q, da),

llGvllLq{n,da) < C ||v||, Vv e M+(Q), (4.5)

where HvH = v(Q). Clearly, (4.5)=^(4.3) if we set dv = Iflda. The converse is also true in many instances.

The following theorem demonstrates that a positive solution u e Lq(Q, a) to (4.1) exists if and only if (4.5) holds. The same condition characterizes the existence of a positive supersolution u e Lq(Q, a) such that

u >G(uqa), 0 <u< da — a.e. in Q. (4.6)

Theorem 4.1. Let a e M+(Q) and 0 < q < 1. Suppose G is a lower semicontinuous, quasi-symmetric kernel which satisfies the weak maximum principle. Then the following statements are equivalent:

1) There exists a positive constant C = C(a, G) such that (4.5) holds.

2) There exists a supersolution u e Lq(Q,da) such that (4.6) holds.

3) There exists a solution u e Lq(Q, da) to (4.1).

This theorem is proved in [23] for more general non-negative kernels under the additional non-degeneracy assumption on the kernel G, which is in fact necessary for the existence of a positive solution. Non-degeneracy of G is not needed for the existence of a non-trivial positive super-solution.

The next theorem gives some necessary and also close sufficient conditions for (4.5) to hold.

Theorem 4.2. Let a e M+(Q), and 0 < q < 1. Suppose G is a quasi-symmetric, non-degenerate kernel which satisfies the weak maximum principle.

q

1) If (4.5) holds, then Ga e L^(Q, a).

Q 1

2) If Ga e L^'\Q, a), then (4.5) holds.

Here LSi1(Q, a) is the corresponding Lorentz space (see [26]).

It is shown in [23] that, without the assumption that G satisfies the weak maximum principle, the condition

/ (Ga)sda < m (4.7)

Jn

with s = is necessary for the existence of a (super)solution u e Lq(Q, a) to (4.1) only if q e (0, q0], where

^^^ =0.61...

denotes the conjugate golden ratio. For general measures a, the existence of a positive solution u e Lq(Q, a) does not guarantee that (4.7) holds if q e (q0,1), or s = for all q e (0,1), even for symmetric positive kernels G.

Another characterization of the strong type inequality (4.5) can be deduced from Mau-rey's results [21]: it is equivalent to the existence of a nonnegative function F e L1(Q, a) which satisfies

sup fG(X,y)F(,)-i<ia(X) <

yenJn

This is a dual reformulation of (4.5), which does not require G to satisfy the weak maximum principle. In the discrete case where Q consists of a finite number of points, it represents the duality of the two basic concave programming problems (see [1, Sec. 5.7]).

We also give in [23; 25; 27] explicit characterizations of the weak type (1, g)-inequality

HGvh^ndo-) < C||v||, for all v e M+(Q), (4.8)

for any > 0, in terms of energy estimates, as well as capacities.

In [23], we demonstrate how to remove the extra assumption imposed in Theorem 4.1 that a (super)solution u e Lq(Q, a) globally. We prove the following theorem where we only assume that u e Lqoc(a), or equivalently, 0 < u < da-a.e., provided the kernel G satisfies a weak form of the complete maximum principle, or alternatively if G is a quasi-metric kernel (see, e.g., [11; 18;23]).

The weak form (with constant b > 1) of the complete comparison principle, also called the domination principle, says the following:

Gf(x) < v(x) on supp (f), Gf < df — a.e. Gf(x) < b v(x) on Q,

for any positive measure f and positive superharmonic function v.

With a special function m satisfying 0 < m < da-a.e., known as a modifier [11; 18], we can modify the kernel G, so that the modified kernel

K(x, y)= G{X'V] , x,y e Q, (4.9)

m(x) m(y)

satisfies the weak maximum principle (2.11). This makes it possible to apply Theorem 4.1 with K in place of G, and consider u e L\oc(Q, a). A typical modifier that works for general kernels G which satisfy the complete maximum principle is given by

m(x) = min{1,G(x,x0)}, x e Q, (4.10)

where x0 is a fixed pole in Q [19].

We say that G is a quasimetric kernel with quasimetric constant k > 0 provided G is symmetric and d(x, y) := ^ y) satisfies the quasimetric triangle inequality with quasimetric constant k, i.e.,

d(x, y) < K[d(x, z) + d(z, y)], (4.11)

for any x,y,z e Q. Here d(x, y): Q x Q ^ [0,

The kernel G is quasimetrically modifiable with constant k if there exists a modifier m such that

K(x, y)= G(X> (4.12)

m(x)m(y)

is a quasimetric kernel with quasimetric constant k.

Theorem 4.3. Let a e M+(Q) and 0 < q < 1. Suppose G is a quasi-symmetric positive kernel, continuous in the extended sense on Q x Q, which either (a) satisfies the complete maximum principle, or (b) is quasi-metrically modifiable. Then the following statements are equivalent:

1) There exists a positive constant C such that the weighted norm inequality

HGvllL^mda) <c[mdv, for all V e M+(Q), (4.13)

Jn

holds, where the modifier m(x) is given by (4.10) for some x0 e Q.

2) There exists a positive (super) solution u to (4.1) such that u e Lloc(Q, a) (or equivalently 0 < u < da-a.e.)

Theorem 4.3 yields a characterization of the existence of weak solutions u e Lloc(Q, a) to the fractional Laplacian equation (1.5) in general domains Q with positive Green's function G for 0 < a < 2, or the entire space R™, or balls, or half-spaces in R™, for 0 < a < n (quasi-metric kernels G) as discussed in the introduction. In the classical case a =2, such solutions are the so-called very weak solutions for bounded C2-domains Q (see, e.g., [20]).

In the case Q = R™, an alternative criterion for the existence of weak solutions u e e L((oc(Rn, a) to the fractional Laplacian equation (3.2) based on localization of (4.5) is discussed in Sec. 3 above.

5. Higher integrability of solutions

Under the assumptions on the kernel G imposed at the beginning of Sec. 4, let us consider conditions on a which ensure the existence of solutions u G Lr(Q, a) with r > 0 to the sublinear integral equation (4.1). In this case, a simple necessary and sufficient condition in terms of the energy (4.7) was obtained in [27]. In particular, these results were used in [24] to establish conditions for the existence of solutions u G Lr(Q,dx) (now with respect to Lebesgue measure) to the fractional Laplacian equation (1.5) in a general domain Q with positive Green's function G, and also W 1,P(Q) solutions to divergence form sublinear elliptic equations (1.3). These results are discussed in Sec. 6 below.

The following theorem (see [27]) characterizes the existence of solutions u G Lr(Q, a) for > , along with related weighted norm inequalities

\\G(fda)\\Lr-(n,a) <c Wfhnn,*), Vf G Lp(Q, a). (5.1)

Theorem 5.1. Let a G ^+(Q). Suppose G is a positive, lower semicontinuous, quasi-symmetric kernel on Q which satisfies the weak maximum principle.

(i) If 1 < p < and 0 < r < p, then (5.1) holds if and only if

pr

( G a) p-rda < (5.2)

'n

(ii) If 0 < q < 1 and q < r < <x, then there exists a positive (super) solution u e e Lr(Q, da) to (4.1) if and only if (5.1) holds with p = ^, or equivalently by statement (i),

f( Ga)^da < (5.3)

Jn

Remark 5.2. We observe that the "if" parts of statements (i) and (ii) of Theorem 5.1 fail if p = 1, and r = q, respectively. The "only if" parts hold for all 0 < r < p in statement (i), and > 0 in statement (ii).

Remark 5.3. It is known that inequality (5.1) with p = ^ > 1 in the case 0 < q < 1 yields the existence of a positive (super) solution u e Lr(Q, a) for (4.1). This statement follows from a lemma due to Gagliardo [14], and does not require G to be quasi-symmetric or to satisfy the weak maximum principle. However, the converse statement fails without the weak maximum principle (see [23] in the case = ).

Remark 5.4. Without the assumption that G satisfies the weak maximum principle, the "only if" parts of statement (i) (with p = ^ > 1) and statement (ii) (with r > q) of Theorem

5.1 hold only for 0 < r < 1 — q2; see [27]. The case r = q when 0 < q < q0 = is discussed above.

6. Positive solutions in W^(fi) and Lr(Ü)

In this section, we present sharp conditions on the coefficient a which ensure the existence a positive weak solution u G W 1,P(Q) or u G Lr(Q, a) to the sublinear elliptic Dirichlet problem

-div(^Vu) = auq, u> 0 in Q, (6

u = 0 on dQ,

[

where 0 < q < 1, and a is a locally integrable function, or measure, in a domain Q C Rra with nonnegative Green's function G associated with the uniformly elliptic operator —div(^V-), for a symmetric matrix A G L'(Q)nxn. Some of these results obtained in [24] hold for analogous fractional Laplacian equations (1.5) in domains Q with positive Green's function.

More generally, for a lower semicontinuous kernel G: Q x Q ^ (0, where Q C C Rra, we denote by

G(fda)(x)= G(x, y)f(y)da(y), x G Q,

JQ

the corresponding integral operator, where a G M+(Q), the class of positive locally finite Borel measures in Q.

We study positive solutions u G Lr(Q, da), as well as u G Lr(Q,dx), for the integral equation

u = G(uqda), 0 <u< da — a.e. in Q, (6.2)

where 0 < q < 1. If G is Green's kernel associated with the Laplacian in a domain Q C Rra, then (6.2) is equivalent to the sublinear elliptic boundary value problem (6.1) (see [23]).

As above, we assume that the kernel G of the integral operator G is lower semicontinuous, quasi-symmetric, and satisfies a weak maximum principle. In Theorem 6.4 below an additional assumption is imposed on G,

C

G(x, y) ^ \x — y\n-a, x,y G Q, (6.3)

for 0 < a < n; here C is a positive constant which does not depend on x, y.

Such restrictions are satisfied by the Green kernel associated with many elliptic operators, including the fractional Laplacian (—A)a and ^.-Laplacian (with a = 2) in a domain

Q Ç R-

The following theorems and corollaries are proved in [24]. Theorem 6.1. Let a G M+(Q), where Q Ç R- (n > 1) is a domain with a positive Green's function G associated with — div(^V-). Let 0 < q < 1 and y > 0. Then there exists a positive solution u to (6.1) such that

if and only if

\Vu\2uy—1dx< ж (6.4)

'n

(Ga)^da < (6.5)

In

In particular, in the case y =1, the condition

i 1 + 9

( G a) ~9da < (6.6)

Jn

is necessary and sufficient for the existence of a positive solution u G Wg'2(Q) to (6.1).

Theorem 6.1 yields sufficient conditions for u G W01,p (Q) when --1 < p < 2. Here W01'P(Q) stands for the homogeneous Sobolev space defined as the closure of C£°(Q) functions in the norm

r 1

\\u\\w 1.P(n) = ( / iVufdxy. n

Letting p = --++Y with 0 < y < 1 in Theorem 6.1, we obtain the following corollary. Corollary

n > 3. If

n—1+Y

Corollary 6.2. Under the assumptions of Theorem 6.1, suppose that nni < P < 2 and

( G a) (1-9)(«-p) 1 da < (6.7)

n

then there exists a positive solution u G Wq'p(Q) to (6.1).

Corollary 6.3. Under the assumptions of Theorem 6.1, suppose a G LS(Q), where 2) < < s — -+2-q(-~2)■ Then there exists a positive solution u G Wq'p(Q) ilLr(Q) to (6.1) with p = ^ - , and r = ^ o .

1 -~s(1+q) -—2s

In the following theorem we give conditions for the existence of positive solutions u G Lr(Q,dx) for more general kernels G.

Theorem 6.4. Let Q Ç R- (n > 1), and a G M+(Q). Let 0 < q < 1. Suppose G is a positive, lower semicontinuous kernel on Q x Q which is quasi-symmetric and satisfies the weak maximum principle.

(i) If there exists a solution u G Lr(Q,dx), r > 0, to (6.2), then

f ( Ga)dx < (6.8)

n

(ii) Suppose additionally that G satisfies condition (6.3) for some 0 < a < n. If

< r < oo and

n n—a

r(n— a)

( G a) —1 da < (6.9)

n

then there exists a positive solution u e Lr(Q,dx) to (6.2).

Corollary 6.5. Under the assumptions of Theorem 6.4, suppose that a e Ls(Q,dx), where n—q(n—a.) < s < a. Then there exists a positive solution u e Lr(Q,dx) to (6.2), where

r = n(l—g)s n—as '

Corollaries 6.3 and 6.5 for bounded domains Q were obtained earlier by Boccardo and Orsina [4]. We observe that the restrictions on s in these corollaries, and consequently on the exponents y, p and r in Theorem 6.1, Theorem 6.4, and Corollary 6.2, are sharp.

There are also weak-type analogues of these results for Marcinkiewicz spaces Lr'^ as in [5], based on the corresponding weak-type weighted norm inequalities (see [23; 25; 27]). The proof of Theorem 6.1 in [24] makes use of the weighted norm inequality

||G( f da)\\Lr (n>da) <C \\f\\Ls(n,dv), Vf e LS(Q, da), (6.10)

for integral operators G in the case 1 < r < s < considered in Sec. 3.

The proof of Theorem 6.4 uses both (6.10) and a similar weighted norm inequality with Lebesgue measure on the left-hand side,

\\G(f da)\Lr(Q,dx) <C y\Ls(n,da), Vf e Ls(Q, da), (6.11)

for 1 < r < to, 1 < s < +to. Such inequalities in the equivalent dual form

\\Gf Wjs(n,da) <C \\nLr'(Q4X), Vf e Lr'(Q, dx), (6.12)

were characterized in the special case where Q = Rra and G = (—A)—a is the Riesz potential of order a e (0,n), by Maz'ya for r > s > 1, and Maz'ya and Netrusov for 1 < r < s, in terms of Sobolev capacities, and later by Cascante, Ortega and the author in terms of Wolff potentials (see [7;22]).

In [24], we employ an alternative characterization in terms of Green's potentials Ga as in Theorem 5.1 above, which is sufficient for both (6.10) and (6.11) simultaneously.

REFERENCES

1. Adams D.R., Hedberg L.I. Function Spaces and Potential Theory. Berlin; Heidelberg; New York, Springer-Verlag, 1996. 368 p.

2. Ancona A. Some results and examples about the behaviour of harmonic functions and Green's functions with respect to second order elliptic operators. Nagoya Math. J., 2002, vol. 165, pp. 123-158.

3. Berge C., Ghouila-Houri A. Programming, Games and Transportation Networks. London, Methuen, 1965. 248 p.

4. Boccardo L., Orsina L. Sublinear elliptic equations in Ls. Houston Math. J., 1994, vol. 20, pp. 99-114.

5. Boccardo L., Orsina L. Sublinear elliptic equations with singular potentials. Adv. Nonlin. Studies, 2012, vol. 12, pp. 187-198.

6. Brezis H., Kamin S. Sublinear elliptic equation on Rra. Manuscr. Math., 1992, vol. 74, pp. 87-106.

7. Cascante C., Ortega J.M., Verbitsky I.E. On Lp-Lq trace inequalities. J. London Math. Soc., (2006), vol. 74, pp. 497-511.

8. Cao D.T., Verbitsky I.E. Nonlinear elliptic equations and intrinsic potentials of Wolff type. J. Funct. Analysis, 2017, vol. 272, pp. 112-165. DOI: http://dx.doi.org/10.1016/jofa.2016.10.010.

9. Cao D.T., Verbitsky I.E. Pointwise estimates of Brezis-Kamin type for solutions of sublinear elliptic equations. Nonlin. Analysis Ser. A: Theory, Methods & Appl., 2016, vol. 146, pp. 1-19.

10. Dat C.T., Verbitsky I.E. Finite energy solutions of quasilinear elliptic equations with sub-natural growth terms. Calc. Var. PDE, 2015, vol. 52, pp. 529-546.

11. Frazier M., Nazarov F., Verbitsky I. Global estimates for kernels of Neumann series and Green's functions. J. London Math. Soc., 2014, vol. 90, pp. 903-918.

12. Fuglede B. On the theory of potentials in locally compact spaces. Acta Math., 1960, vol. 103, pp. 139-215.

13. Grigor'yan A. Heat Kernel and Analysis on Manifolds. Boston, AMS, 2009. 482 p.

14. Grigor'yan A. Heat kernels on weighted manifolds and applications. Contemp. Math., 2006, vol. 398, pp. 93-191.

15. Grigor'yan A., Hansen W. Lower estimates for a perturbed Green function. J. Anal. Math., 2008, vol. 104, pp. 25-58.

16. Grigor'yan A., Verbitsky I.E. Pointwise estimates of solutions to nonlinear equations for nonlocal operators, preprint, 2017.

17. Grigor'yan A., Verbitsky I.E. Pointwise estimates of solutions to semilinear elliptic equations and inequalities. J. d'Analyse Math. (to appear). URL: https://arxiv.org/abs/1511.03188.

18. Hansen W. Uniform boundary Harnack principle and generalized triangle property. J. Funct. Analysis, 2005, vol. 226, pp. 452-484.

19. Hansen W., Netuka I. On the Picard principle for A + ^ Math. Z., 2012, vol. 270, pp. 783-807.

20. Marcus M., Veron L. Nonlinear Second Order Elliptic Equations Involving Measures. Berlin; Boston, Walter de Gruyter, 2014. 248 p.

21. Maurey B. Theoremes de factorisation pour les opérateurs linéaires à valeurs dans les espaces LP. Paris, Soc. Math. France, 1974. 384 p.

22. Maz'ya V. Sobolev Spaces, with Applications to Elliptic Partial Differential Equations. Berlin, Springer, 2011. 866 p.

23. Quinn S., Verbitsky I.E. A sublinear version of Schur's lemma and elliptic PDE, preprint. URL: https://arxiv.org/abs/1702.02682.

24. Quinn S., Verbitsky I.E. Solutions to sublinear elliptic equations in W 1'P(Q) and Lr(Q) spaces, preprint, 2017.

25. Quinn S., Verbitsky I.E. Weighted norm inequalities of (1,g)-type for integral and fractional maximal operators. Harmonic Analysis, Partial Differential Equations and Applications, in Honor of Richard L. Wheeden, Ser. Applied and Numerical Harmonic Analysis. Basel, Birkhauser, 2017, pp. 217-238.

26. Stein E., Weiss G. Introduction to Fourier Analysis on Euclidean Spaces. Princeton, N.J., Princeton Univ. Press, 1971. 297 p.

27. Verbitsky I.E. Sublinear equations and Schur's test for integral operators, 50 Years with Hardy Spaces, a Tribute to Victor Havin, Ser. Operator Theory: Adv. Appl., vol. 261, Basel, Birkhauser, 2017 (to appear).

ПОТОЧЕЧНЫЕ ОЦЕНКИ И КРИТЕРИЙ СУЩЕСТВОВАНИЯ РЕШЕНИЙ СУБЛИНЕЙНЫХ ЭЛЛИПТИЧЕСКИХ УРАВНЕНИЙ

Игорь Е. Вербицкий

Профессор, Факультет математики, University of Missouri verbitskyi@missouri.edu Columbia, MO 65211, USA

Аннотация. В работе представлен обзор последних результатов о положительных решениях эллиптических уравнений типа —Lu + V и9 = /, где L — эллиптический оператор в дивергентной форме, 0 < q < 1, f > 0 и V — функция, которая может изменять знак, в области Q С Мга или на весовом римановом многообразии с положительной функцией Грина G. Обсуждаются вопросы существования решений, глобальные нижние и верхние поточечные оценки классических и слабых решений u, а также условия, обеспечивающие и Е Lr(Q) или и Е W 1'P(Q).

Некоторые из этих результатов применимы к однородным сублинейным интегральным уравнениям и = G(uqda) in Q, где 0 < q < 1, а а = —V — положительная локально конечная борелевская мера в Q. Здесь G(f da)(x) = = fnG(x, у), f(y) da (у) — интегральный оператор с положительным (квази) симметричным ядром G на Q х Q, который удовлетворяет слабому принципу максимума. Результаты распространяются на положительные, возможно сингулярные, решения сублинейных уравнений, содержащих дробный лапласиан

(—А) аи = а и9, и > 0 в Q,

где 0 < q < 1, 0 < а < п и и = 0 в Qc и на бесконечности в областях Q С Мга с положительной функцией Грина G.

Ключевые слова: сублинейные эллиптические уравнения, функция Грина, слабый принцип максимума, дробный лапласиан.