Nonlinear equations with a generalized fractional Laplacian

We study a generalization of the power of Laplace operator with null Dirichlet conditions by means of the spectral theory and prove several existence results for nonlinear equations with such operators, especially when the problem is resonant. Some regularity results are also obtained.


Introduction
Fractional differential operators are modern tools in the modeling of many phenomena of mechanics, electric circuits and material sciences. However, this notion has several meanings starting from fractional derivatives of Riemann-Liouville, Caputo or other kinds, through pseudo-differential and singular integral operators up to definitions given by the Spectral Theory. Some definitions are equivalent though they seem to be completely different, see [16]. They have rather global nature in a sense that functions being their arguments have to be defined on the whole R k . If such operators are applied in equations on an open bounded set ⊂ R k with the null Dirichlet's boundary condition, then the authors set u = 0 outside . Operators of this kind are generalizations of the Laplace operator and they show some probabilistic motivation (Lévy processes) of models considered [2,16,17,24].
In this paper, we use operator defined by the Dirichlet Laplacian − : L 2 ( ) ⊃ dom(− ) → L 2 ( ) by using the Spectral Theorem and functional calculus given in [8]. When g is a real function defined on the spectrum of − , then g(− ) is a certain self-adjoint operator on L 2 ( ). This approach includes positive powers (− ) s , but operators defined in this way are different from the ones considered above. The arguments in [23] show that eigenfunctions of global (− ) s differ from the eigenfunctions of the Laplacian which are the same for g(− ) for any function g; earlier the eigenvalues of these both operators were compared in [6]. Some existence and regularity results for equations involving spectral fractional Laplacians are given in [1,4,5,14,21,25]. It would be interesting to compare results of B Igor Kossowski igor.kossowski@p.lodz.pl 1 Institute of Mathematics, Lodz University of Technology, ul. Wólczańska 215, 93-005 Łódź, Poland this paper with theorems from the above mentioned articles even for fractional Laplacians. However, [1,4,5,25] consider superlinear nonlinearities and our approach admits only ones with linear or sublinear growth (comp. (3.2)). In [21], properties of the fractional Laplacians with Navier and Dirichlet boundary conditions are studied as linear operators on different Sobolev spaces. The closest to our setting is [14], where the fractional Laplacian is replaced by The operator acts in L 2 ( ) as in our approach and assumptions on parameters give the compactness of the inverse operator. However, the nonlinear term f is potential, since the method used in [14] is variational (here, we use the fixed point and topological degree methods). As in the case of usual Laplace operator, the results obtained by using the variational method need slightly weaker assumptions as obtained by the topological one. We also refer to some applications of fractional Laplacian [3,18] and fractional p-Laplacian [7,11,26,27].
There is no reason to restrict to power functions g(z) = z s , since engineers using fractional operators rely on experimental arguments: numerical results obtained with fractional powers of differential operators are closer to experimental data. It is possible that some deviations from power functions will give even better results. The paper is organised as follows: in the next section we remind some facts concerning the Dirichlet Laplacian and the Spectral Theory. In Sect. 3, the existence of solutions for g(− )u = f (x, u) with invertible g(− ) will be obtained. The proof of these theorems are simple. More interesting is the resonant case, in which this operator has a nontrivial kernel. It happens when g vanishes on some eigenvalue of the Laplacian. The existence theorem is then given with conditions similar to these due to Landesman and Lazer [15] and the proof is based on the coincidence degree technique of Mawhin. The last section is devoted to regularity of solutions to such equations. The results obtained there are not satisfactory in general case when any bounds of the sequence of L 2 and supremum norms of derivatives of every order are not known. In fact, we can use these results for cuboids.
The theorems of this paper are completely analogous to those available for case of usual Dirichlet Laplacians g(s) = s. The existence and uniqueness of a solution obtained by the Schauder Fixed Point Theorem or the Banach Contraction Principle (Theorems 1 and 2) are essentially weaker than the known ones in the case of the fractional Laplacians g(s) = s β , however the novelty of this paper is about the replacement of the power function by an arbitrary one. The main interesting resonant case (Theorem 4) cannot be applied to fractional Laplacians since s β = 0 only for s = 0 which is not an eigenvalue of the Dirichlet Laplacian.
The authors are very grateful to the Reviewers for their helpful comments, corrections of misprints and improvements of the text. Special acknowledgements are due to the Reviewer 1 who noticed essential error in the proof of Theorem 5 and proposed a solution.

Dirichlet-Laplace operator and its generalizations
Let be an open bounded subset of R k with Lipschitzian boundary and let − : L 2 ( ) ⊃ dom(− ) → L 2 ( ) be the weak Dirichlet-Laplacian (compare [14]), i.e. it is Friedrich's extension of the classical Dirichlet-Laplace operator to a self-adjoint operator.
Moreover, for being of class C 2 or being cuboid we have − = (− ) s .
The operator − is self-adjoint and positive, its spectrum σ (− ) consists of positive eigenvalues λ n , n ∈ N, with finite multiplicities being a nondecreasing sequence tending to infinity. If we put each eigenvalue as many times in the sequence as its multiplicity and by e n denote a corresponding normalized in L 2 eigenfunction than the operator can be written as − u = ∞ n=1 λ n u, e n e n . (2.1) The spectral theory [8] provides us, so called, the operational calculus. If g : σ (− ) → R is any function then a new operator g(− ) can be defined by the formula In the special case g(s) = s β this operator is the fractional Laplacian. Generally, the correspondence g → g(− ) has the following algebraic properties: where g 1 , g 2 are real functions on the spectrum, c is a real constant and C ⊃ B means that operator C is an extension of operator B. g(− ) is a self-adjoint operator for any real function g, it is positive, if g takes positive values, it is invertible, if g(λ n ) = 0 for any n. If b := inf n |g(λ n )| > 0, then g(− ) −1 is defined on the whole L 2 ( ) : u, e n e n (2.3) and it is bounded. Moreover, if lim then this inverse operator is compact. However, if (2.4) holds but g vanishes for a finite set G := {n ∈ N : g(λ n ) = 0}, then the resolvent is meromorphic with the pole 0 an it has the form where lim λ→0 |c n (λ)| = ∞ and R 0 is holomorphic and compact in a neighborhood of 0.
We also need some properties of eigenfunctions e n . They are of C ∞ -class in (even analytic, see [12], Cor. 8.11). It is well-known that the first eigenvalue is simple, hence λ 1 < λ 2 , and the corresponding eigenfunction e 1 is positive in .
In the one dimensional case eigenvalues can be found explicitely. Assume for simplicity = (0, π), − u = −u , then λ n = n 2 , e n (x) = √ 2/π sin(nx). In multidimensional case, one can find them for cuboids, some other polytopes or for balls. We write down the eigenvalues and the eigenfunctions in the case = (0, π) k : (one can find eigenvalues and eigenfunctions for general cuboids, since they can be obtained from (0, π) k by using an appropriate homothety, translation and orthogonal mappings) and in the case where λ n,m stands for m-th positive zero of the Bessel function J n , for n = 1, 2, . . . they are double eigenvalues (comp. [13]). In a general case ⊂ R k [13], we only know that if 1 ⊂ 2 , then Moreover, there are two positive constants depending on such that

Elliptic equations with generalized fractional Laplacian
In this section we are interested in the following equations: The last condition is necessary and suffient for the Nemytskii operator defined by f , i.e. f(u)(x) = f (x, u(x)) to act in L 2 ( ) as a continuous mapping and map bounded sets into bounded ones.
where b is a constant from (3.2), then equation (3.1) has a solution.
Proof It is easily seen that (3.1) is equivalent to the existence of a fixed point of operator where · stands for the L 2 -norm. Since g(− ) −1 = sup n |g(λ n )| −1 < b −1 by our assumption, hence Uniqueness of solutions can be obtained under another assumptions.
Theorem 2 Let f be a Carathéodory function, f (·, 0) ∈ L 2 ( ) and there exists L > 0 such that Proof Notice that assumptions on f from this theorem imply (3.2) with a(·) = f (·, 0) and b = L. Moreover and our assumption gives g(− ) −1 f is a contraction on L 2 ( ). The assertion follows from the Banach Contraction Principle.
The above theorems consider equations with invertible linear part g(− ). We can also study resonant equations, when function g vanishes on some eigenvalues. The first result for the Laplace operator was obtained by Landesman and Lazer in 1970 and we generalize it for our operators. The most satisfactory result can be obtained in the case g(λ 1 ) = 0, g satisfies condition (2.4) and inf n≥2 |g(λ n )| > 0.
Lemma 1 Assume f : × R → R is a Carathéodory function satisfying (3.2). If there exist limits lim

4)
uniform with respect to x and f ± ∈ L 2 ( ), then there is a function h ∈ L 2 ( ) such that for any (x, u).
Proof Due to (3.4) there exists c > 0 such that Then the assertion holds with In 1972 Jean Mawhin [19] introduced a topological degree for pairs of maps (L, N ), where L was a linear Fredholm map of index 0, L : X ⊃ domL → Z , and N : X → Z which generalized the Leray-Schauder theory for pairs (I , N ) with N completely continuous maps. We follow two later monographs of this theory [10,20]. Let X and Z be Banach spaces, L : X ⊃ Y → Z be a linear operator such that dim ker L = codim imL > 0. (Such operators are called Fredholm with index 0. More generally, Fredholm operators have finite dimensional kernels, closed images with finite dimensional Z /imL. The index of the operator is the difference of these dimensions.) Let N : X → Z be continuous (and nonlinear). Let P be a linear projector on ker L and Q a linear projector in Z along im L. Denote by K P the inverse of The operator N is called L-compact if Q N and K P (I − Q)N are completely continuous (usually K P is completely continuous and N maps bounded sets into bounded ones).

then (3.1) has a solution.
Proof It is obvious that g(− ) has index 0 with kernel and cokernel being one-dimensional spaces spanned by e 1 . Using the notation of Mawhin's Continuation Theorem where f(u) ≤ h for any u ∈ L 2 ( ) by the lemma. This operator is then compact due to We shall show that equations g(− )u = μf(u) have no solution for μ ∈ (0, 1] with the norm sufficiently large. Suppose the converse: there are μ j ∈ (0, 1] and u j ∈ L 2 , such that Denoting byũ the orthogonal projection of u on ker L ⊥ , we can write u j =ũ j + d j e 1 . Then we have hence the norms ofũ j are uniformly bounded so the sequence |d j | tends to the infinity. We shall show that both possibilities: there is a subsequence d j → +∞ (resp. d j → −∞) lead to a contradiction in the same way. Take the first option into account and let (d j ) stands for the corresponding subsequence. Denote the sequence of integrals of e 2 1 over A j tends to 0. Take any ε > 0. Hence where μ k denotes the k-dimensional Lebesgue measure, thus μ k (A j ) → 0, since ε > 0 is arbitrary. Passing to a subsequence, we can assume Due to the Borel-Cantelli Lemma On the other hand e 1 is positive on , thus 1 2 d j e 1 (x) → +∞ as j → ∞ and, thereforẽ u j (x) + d j e 1 (x) → +∞ a.e. and Let us multiply both sides of equation (3.6) by e 1 and integrate it over and divide by μ j . The left-hand integral vanishes since the image of g(− ) is orthogonal to e 1 . The right-hand side tends to f + e 1 by Lebesgue's Dominated Convergence Theorem, so it is nonzero for sufficiently large j. We have a contradiction. The second option gives the same with f − instead of f + .
Thus, we can take R > 0 so large that all solutions of g(− )u = μf(u) sit in the open ball B(0, R) in L 2 ( ). We should calculate the Brouwer degree of the mapping on the interval (−R, R). But these mappings take values with the same signs as (3.5) on both ends of this interval. Therefore this degree is +1 or −1 in the cases respectively. By Mawhin's Continuation Theorem we get the assertion.
Similar theorem can be obtained replacing λ 1 by any simple eigenvalue λ m . Since the corresponding eigenfunction e m change signs, the main assumption (3.5) should be changed by: the integrals One can also replace the limits in definitions f ± by limes inferior and superior as in [9]. Namely, if u)e 1 and the limit superior and limit inferior exist uniformly with respect to x, then a solution exists.

Some regularity results
There could be obtained two kinds of regularity theorems. Both of them rely on some boundedness properties of the sequence of eigenfunctions e n which can be checked in some cases. However, in general, they are difficult to verify. Assume that u is a solution of (3.1) where function g can vanish for a finite number of λ n and f satisfies (3.2). We consider the question when a solution u belongs to the Sobolev space H p ( ). In this section, we assume that all eigenfunctions are sufficiently smooth up to the boundary.
is given by corresponding scalar products. We need an estimate: for a positive constant C depending on only. In particular, it means that e n ∈ H p ( ) for n ∈ N. The second assumption is the orthogonality of vectors ∇ p e n and ∇ p e m for n = m and any p : ∇ p e n , ∇ p e m = 0, n = m, p ∈ N. The first summands equals 0 for m = n or equals λ n ∇ p−1 e n for m = n by induction and the second ones are given by the assumptions of this lemma.
It is known that if domain is C p,α , then e n ∈ C p,α ( ) (compare Theorem 6.19 in [12]), hence ∇ p e n is finite. In general, ∇ p e n need not to be finite if p ≥ 2. We will not focus on the question of the finiteness of this norm, since (4.1), (4.2) require much more than this problem. We have written eigenvalues and eigenfunctions for normalized cuboids (2.5). The eigenfunctions are from C ∞ ( ). One can easily check that (4.1) holds for cuboids with C = 1 (both sides are then equal). Simple calculations show that (4.2) is also true for these normalized eigenfunctions. General cuboids are images of the normalized one by an appropriate affine mapping being a composition of an orthogonal map, a homothety and a translation. These maps change eigenvalues multiplying them by c k , where c is the ratio of the homothety. The eigenfunctions are the composition of the normalized one with the above affine map. Thus, e n ∈ C ∞ ( ), n ∈ N, and relations (4.1), (4.2) hold for a general cuboid.
We know also the explicite form of eigenfunctions for a disk. However, in this case the problem is with the orthogonality of higher derivatives of different eigenfunctions and we are not able to check the needed conditions. Finally, we want to emphasize that, to the best of our knowledge, conditions (4.1) and (4.2) do not follow from known properties of eigenfunctions for general domains. Only the explicit form of eigenfunctions can be used to obtain estimates (4.1) and (4.2).
Proof Letũ denotes the orthogonal projection of u on ker g(− ) ⊥ and u =ũ + n∈G d n e n .
Since eigenfunctions e n , n / ∈ G, belong to H p ( ), we should only check ifũ ∈ H p ( ). Butũ where all entries of the series sit in H p ( ). By (4.2), for m ≤ p, we have from (4.1) and (4.3) that, for m ≤ p, thus we are done.
Remark If function g is nondecreasing then one can use estimates c 1 n 2/k ≤ λ n ≤ c 2 n 2/k since contains a cube and it is a subset of another cube. Thus the assumptions guaranteeing regularity of both kinds do not require an explicit knowledge of eigenvalues but estimates on the growth of norms ∇ p e n are still needed. The proof is obvious since all partial derivatives ofũ up to the order p are sums of uniformly convergent series of continuous functions -multiples of D α e n . The crucial assumption (4.4) holds for cubes in any dimension.

Examples
Let := (0, π) ⊂ R, − u = −u , then λ n = n 2 , n ∈ N. Consider a fractional power: g(s) := s β , where β > 0. Then inf |g(λ n )| = 1 and we can use Theorem 1 for the following nonlinearities: The function f given by the last two formulas satisfy assumptions of Theorem 2, thus we have exactly one solution in these cases. Now, we consider the same problem with g(s) = s and f (x, u) = u, then all assumptions of Theorem 1 will be satisfied except b < inf n g(λ n ) -here we have the equality instead the sharp inequality. The same statement is true for Theorem 2 and the assumption L < inf n g(λ n ). However the problem has an obvious solution u(x) = sin x. If we change the nonlinear term and put f (x, u) = u − 1, then the assumptions of both Theorem 1 and 2 are satisfied except the sharp inequalities but the problem has no solution.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.