Evolution Equations with Sectorial Operator on Fractional Power Scales

Originating with the famous monograph by Dan Henry, the semigroup approach to evolution problems having a positive sectorial operator in the main part allows us to settle them at various levels of the fractional power scale associated with the main linear operator. This translates into different regularity properties of local solutions to such equations. Specific applications of the abstract results to the 2D surface quasi-geostrophic equation or the fractional chemotaxis system are presented.


Introduction
This paper is devoted to the intensively studied problems of regularity and global in time continuation of the local solutions to semilinear evolution equations. Our approach follows closely the semigroup approach known from the classical monographs by Henry [24] and Pazy [36] and makes use of the important results of Triebel [43] and Amann [1] concerning the construction of the fractional power scales connected with positive operators. Our aim here is to discuss applications of these abstract results for semilinear equations with sectorial positive operator in the main part to a class of new examples. The paper consists of two parts. In the first one we recall known from the literature abstract facts needed in further considerations. A variety of formulas and properties of fractional powers of positive operators and semigroups generated by sectorial operators are first recalled, with detailed presentation to be found e.g. in [1,6,14,24,25,33]. Next we list basic properties of the fractional power spaces, in which the semilinear equations will be studied. In the literature there are known a few papers and monographs dealing with fractional power scales and equations on it. From among of them we select here, [1,5,8,16,23,26,29,31,38,39]. Not only do we discuss in Sect. 2 the one-sided scale associated to non-negative fractional powers of positive operators, but also extend the scale to the negative side by extrapolating the main positive sectorial operator in a reflexive Banach space. This allows us to consider evolution equations in the extrapolated fractional power scale, which in applications may include spaces of continuous functionals on subspaces of classical Sobolev spaces. A brief recapitulation of the local solvability result by Henry [24] (later extended in [6]) in the new context is also mentioned. The abstract technique used in this paper is based on locating the considered problem at an admissible level of the fractional power scale associated with the sectorial or elliptic operator from the main part. Once the level is set, it provides the base space in which we need to place the action of the nonlinearity, specific in a particular example. We need to assure that the nonlinearity acts between suitably chosen levels of that fractional power scale; more precisely, it is Lipschitz continuous on bounded sets as a map between elements of the scale-the phase space and the base space. Then our original studies are reported starting with the global extendibility for small data solutions in Lemma 2. 28.
The second part of the paper is devoted to two examples known from the recent references [3,10,12,44]-the surface quasi-geostrophic equation and a version of the chemotaxis system. Using presented earlier abstract results and constructing suitable estimates, we discuss a possible extent of regularity of the local solutions in these examples. The problems studied here are subject to Dirichlet boundary conditions in bounded regular domains, which is different from examples of Cauchy problems in the whole of R N treated in some of the above-mentioned references, where the complications connected with boundary conditions are absent. Very helpful in these studies is the idea introduced by Giga and Miyakawa [22] in case of the celebrated Navier-Stokes equation. Thus, in some sense, we continue here the direction used in [7], where the results of [22] were used. We also give an application of Lemma 2.28 on small data solutions to the studied examples.
Considering (1.1) in the one-sided fractional power scale generated by A = (− ) ω in X 0 = L 2 ( ), we show in Theorem 3.4 that unique local (X α , X 0 ) solutions of (1.1) exist for α ∈ 1 2ω , 1 and f ∈ L 2 ( ). Knowing that f ∈ L ∞ ( ), they are in fact global in time and bounded as shown in Proposition 3.5. Smoother solutions of 2D quasi-geostrophic equation can be obtained by considering (1.1) in the base space from the positive part of the one-sided scale. In Theorem 3.6 we show that if f ∈ X β = D((− ) ωβ ), 0 < β < 1 4ω , and α ∈ β + 1 2ω , β + 1 , then local (X α , X β ) solutions of (1.1) exist. We use Lemma 2.28 in order to determine in Proposition 3.8 the range of parameters for which these solutions are in fact global in time, provided that the initial data is sufficiently small.
In Sect. 3.3, borrowing the idea from the paper [22] on the Navier-Stokes equation, we consider (1.1) in the spaces of negative part of the extrapolated fractional scale generated by A. We prove in Theorem 3.11 that if f ∈ X β with β ∈ 1 ω − 2, 0 and α ∈ 1 2 β + 1 2ω , β + 1 , then for each θ 0 ∈ X α there exists a unique (X α , X β ) solution of (1.1). Invoking the natural L ∞ ( ) a priori estimate for the 2D quasi-geostrophic equation, in Proposition 3.14 we show that these (X α , X β ) solutions of (1.1) are in fact global in time for f ∈ L ∞ ( ), β ∈ 1 2ω − 1, 0 and 1 2ω < α < β + 1. A similar analysis is made in Sect. 4 devoted to the abstract Cauchy problem u t + Au = F(u), u(0) = u 0 , (1.2) related to the fractional chemotaxis system in a bounded C 2 domain of R N u t + (− ) ω u + ∇ · (u∇v) = 0, x ∈ , t > 0, v + u = 0, x ∈ , t > 0, with a parameter ω ∈ 1 2 , 1 . Local (X α , X 0 ) solutions of (1.2) with α ∈ 1 2ω , 1 are obtained in Theorem 4.1 and Remark 4.2 in X 0 = L p ( ), p > N , and X 0 = L 2 ( ) with N = 3, respectively. The latter ones are global in time for small initial conditions, see Corollary 4.4. The corresponding result on local solvability of fractional chemotaxis equation in smoother spaces is given in Theorem 4.5. Finally, the problem (1.2) considered on the extrapolated part of the fractional powers scale generated by (− ) ω , ω ∈ 1 2 , 1 , in L 2 ( ) with N ≥ 2 is shown in Theorem 4.7 to possess local (X α , X β ) solutions provided that β ∈ (−1, 0] and α > 0 are such that 2α − β ≥ N 4ω and 1 2ω < α − β < 1. There are two notational conventions used throughout the text. For a real number s by s − we understand a number strictly smaller but eventually close to s. For conve-nience we also use the letter c to denote a positive constant that can differ even within one calculation.

Scales of Banach Spaces Associated with Positive Operators
A partial differential equation which can be written in the semilinear form has various interpretations depending on the choice of the base space in which the realization of the linear differential operator A is understood. This choice is dictated mainly by the boundary conditions associated with the equation, if any, by the properties of the realization we want to exploit in the arguments, and by the suitable interplay between the linear part and the qualities of the nonlinear part F with respect to the subspaces of the base space. All these features eventually become determined in a suitable definition of a solution in order to prove its existence for a sufficiently large set of initial conditions and to study its long time behavior as time evolves. One of the successful approaches to treat large classes of partial differential equations emerged forty years ago in the monograph of Henry [24] where a given realization of the linear operator A in (2.1) in a chosen Banach space was a sectorial operator.
where is a path from a + ∞e −iθ to a + ∞e iθ for θ ∈ φ, π 2 and surrounding a in order to avoid [a, ∞). If Re σ (A) > a then for some C > 0 we have

Fractional Powers
Positive sectorial operators are in particular positive operators for which fractional powers can be defined.

Definition 2.4 We say that
In fact, if A is a positive operator, then there exist 0 < ε < π, r > 0 and where is a path from ∞e −iθ to ∞e iθ for θ ∈ [ε, π) and surrounding 0 in order to avoid (−∞, 0], whereas is a path from ∞e −i(π −θ) to ∞e i(π −θ) avoiding [0, ∞). For z = 0, A 0 is defined as an identity operator I on X .
Using Euler's function, these fractional powers can be expressed as for a given n ∈ N and z ∈ C such that 0 < Re z < n, see [14, Theorem 3.0.7]. The operator A −z with Re z > 0 is invertible from X onto im(A −z ) which allows to extend the above definition.

Definition 2.6
For z ∈ C such that Re z > 0 we define the operator The fractional powers of positive operators enjoy many useful properties including the following.

Proposition 2.7 A z with
Re z > 0 is a densely defined closed linear operator and for n ∈ N the operator A n coincides with the nth iterate of A. If z 1 , z 2 ∈ C are such that For α, β ∈ R we also have In the vein of (2.3), the fractional powers of exponents with positive real part can be alternatively expressed as we derive from (2.5) the Balakrishnan formulas (see [33, (3.1)], [25]), that is, for n ∈ N 0 and 0 < Re z < 1 we have Note that by (2.4) in the above formula we can write Lemma 2.8 ([1, (V.1.2.12)]) Let −∞ < γ < β < α < ∞ and A : X ⊇ dom(A) → X be a positive operator in a Banach space X . Then the moment inequality holds, that is, If A is a positive operator in a reflexive Banach space X , then its adjoint A * is also a positive operator, since they have the same resolvent set and resolvent estimate and in (2.2) one can change A into A * and X into X * with the same ε, K 1 and r .

Proposition 2.10 If A is a positive sectorial operator in a Banach space X , then for
Re z > 0 we have Given any γ > 0, e −At : and there exist a > 0 and C γ > 0 such that , if A is a (positive) sectorial operator in a reflexive Banach space X , then A * is also a (positive) sectorial operator in X * and (e −At ) * = e −A * t , t ≥ 0.

One-Sided Fractional Power Scale
Having recalled the main properties of the fractional powers of positive operators, we now observe that they define a scale of Banach spaces and give rise to a scale of positive operators.

Definition 2.11
Let A : X ⊇ dom(A) → X be a positive operator in a Banach space X . For any α ≥ 0 we consider its fractional power A α , which is a densely defined closed linear operator with 0 ∈ ρ(A α ); A 0 = I . We define the fractional power space X α A as dom(A α ) endowed with the equivalent graph norm of A α , that is, The space X α A is a Banach space and a dense subspace of X . If 0 ≤ α < β then X is an isometry between Banach spaces, i.e., By the moment inequality (Lemma 2.8) we obtain a relation for the norms of fractional power spaces: Along with the fractional power spaces, corresponding operators are defined as Thus we get for 0 ≤ α ≤ β is commutative. Observe that by (2.4) we also have which leads to the following result.
Lemma 2.14 For A as in Definition 2.11, A α is a densely defined closed operator in X α and 0 ∈ ρ(A α ). Moreover, ρ(A) ⊆ ρ(A α ) and R(λ : In particular, A α is a positive operator (of the same type as A).

Definition 2.17
Let A : X ⊇ dom(A) → X be a positive operator in a Banach space X . Then the family of triples (X α A , · α,A , A α ) α≥0 constructed above is called a onesided fractional power scale generated by A. Hence (X α A , · α,A ), α ≥ 0, is a family of Banach spaces and A α ∈ Liso(X α+1 If no confusion arises, we drop the subscript A and denote the space by X α and the norm by · α in the one-sided fractional power scale. Lemma 2.18 If A is positive sectorial, then, for a given α ≥ 0, its realization A α in X α A is also positive sectorial. Moreover, e −A α t = e −At | X α A and there exists a > 0 such that for γ ≥ 0 we have with some C γ > 0 A one-sided fractional power scale can be generated by the realization of a given positive operator in a fractional power space.
Then the one-sided fractional power scale generated by A γ coincides with shifted by γ ≥ 0 the one-sided fractional power scale generated by A, that is, A one-sided fractional power scale can be generated by the adjoint operator of a given positive operator in a reflexive Banach space.

Example 2.20
Let A : X ⊇ dom(A) → X be a positive operator in a reflexive Banach space X and let A * : X * ⊇ dom(A * ) → X * denote its adjoint operator. Then α≥0 is a one-sided fractional power scale generated by A * . A one-sided fractional power scale can also be generated by the extrapolated operator of a given positive operator, which we describe below.
Let A : X ⊇ dom(A) → X be a positive operator in a Banach space (X , · X ). Since 0 ∈ ρ(A), we have A −1 ∈ L(X ) and |x| X = A −1 x X , x ∈ X , is another norm on X .
We consider a completion (X −1 , · −1 ) of the normed space (X , |·| X ) via an isometric embedding j : (X , |·| X ) → (X −1 , · −1 ). We call the Banach space X −1 an extrapolated space. Recall that j(X ) is a dense subspace of X −1 and Therefore, using j we consider A with domain in X −1 and values in X −1 , that is, Since dom(A) is dense in X , there exists a unique linear operator which extends A X −1 and such that A −1 • j ∈ Liso(X , X −1 ), i.e., A −1 • j is an isometry from X (endowed with the original norm · X ) onto (X −1 , · −1 ), so We call A −1 the extrapolated operator A in the extrapolated space X −1 . Thus we have so it is understood as a realization of the operator A in X −1 in this sense.
Using the definition of a fractional power of a positive operator and properties of the extrapolated operator, we obtain the following relations.
Note that in the above the choice of j determines X −1 . Without its specification X −1 is only determined up to an isometry. A particular choice of an isometric embed- and (X * ) α A * are reflexive Banach spaces. Moreover, we have dense and continuous is injective and has a dense image in equipped with the norm · (X * ) α A * . Denoting by J : X → X * * the canonical injection, which by the reflexivity of X is an isometry, we define the injective map j α :

Proposition 2.24
For α ≥ 0 the map j α is an isometric embedding between X normed by A −α · X and [(X * ) α A * ] * endowed with the standard norm and Taking α = 1 in the above proposition, we see that the isometric embedding j = j 1 allows to construct one-sided fractional power scale generated by A −1 .

Extrapolated Fractional Power Scale in a Reflexive Banach Space
The equations of type (2.1) can be considered in a Banach space for which the operator A is realized as a positive sectorial operator A : X ⊇ dom(A) → X and hence in any space X α A from the one-sided fractional power scale generated by A, that is, a subspace of the base space X . However, using the extrapolated fractional power scale one can realize the operator even in a bigger space using its adjoint operator, as it has already been observed by many authors, see e.g. [30]. A comprehensive description of the extrapolated fractional power scale can be found e.g. in [38, p. 14], where the presentation follows the abstract approach of [1]. We outline it briefly below.
We consider a positive operator A : X ⊇ dom(A) → X in a reflexive Banach space X . Using the isometric embedding j = j 1 = (i * ,1,0 ) * • J ∈ L(X , [(X * ) 1 A * ] * ), we construct the extrapolated operator A −1 : Thus there exists a one-sided fractional power scale We define a family of triples which we call extrapolated fractional power scale generated by A. The spaces X α , α ≥ −1, are reflexive Banach spaces equipped with the norm All spaces X α , α ≥ −1, are subspaces of X −1 and X β is densely and continuously injected into X α provided that −1 ≤ α ≤ β. Note also that the moment inequality Denoting the injection i β,α : X β d → X α , we see that Using Lemma 2.22, for α ≥ 0 the spaces X α are characterized in terms of the previously defined spaces as In particular, we have The spaces X α for α ∈ (−1, 0) have the following characterization.
Iterating the construction of extrapolated operators, for a chosen N ∈ N, a fractional power scale (X α , , but one has to be careful choosing the appropriate isometric embeddings (see [1, p. 263]). Since in our applications we do not go beyond the space [(X * ) 1 A * ] * , we restrict our presentation to the case N = 1 as described above. Moreover, in the applications we will identify X with j(X ) and use Banach spaces isometric to X α . Finally, to simplify the notation, considering an extrapolated fractional power scale, instead of A α,−1 we just write A α , which should not be confused with an operator from the one-sided fractional power scale.

Solutions of Sectorial Equations
Interpreting the linear differential operator A from (2.1) in a suitable (reflexive) Banach space X as a positive sectorial operator A in X , the construction of fractional power scales described above gives us a natural family of Banach spaces (X α , · α ) α∈J , J being the interval [0, ∞) in the case of one-sided or [−1, ∞) in the case of extrapolated fractional power scale, in which the connected with that operator linear abstract Cauchy problem will be considered in the base space X β , β ∈ J . Such a choice of spaces is specific for the operator A. Very often it will coincide with a family of closed subspaces of the Sobolev spaces, as seen in the examples described further in the text. It is also well-known that even in the simplest case, when A is a realization of a uniformly elliptic operator, the spaces obtained in this way may not in general coincide with the whole standard Sobolev spaces, being their proper closed linear subspaces.
Having established the fractional power scale, we verify that the nonlinearity F in (2.1) gives rise to its Nemitskii operator F : X α → X β with some α, β ∈ J , 0 ≤ α − β < 1, which is Lipschitz continuous on bounded subsets of X α . Thus we consider the semilinear abstract Cauchy problem with the first equation satisfied in the base space X β . Note that the operator A β takes isometrically X β+1 onto X β , cp. (2.10) or (2.14), and we have X β+1 d → X α d → X β . By the theory of D. Henry (see [6,24]) for each u 0 ∈ X α there exists a unique (X α , X β ) solution of (2.15), that is, a function in the class defined on the maximal interval of existence and satisfying the first equation in (2.15) in X β and the initial condition in (2.15) in X α . In particular, the Duhamel formula holds, where e −A β t is the linear analytic semigroup generated by −A β , which coincides with T (t) = e −At or with T (t) = e −A −1 t on X β , depending on the scale we consider, see also [29, (3 .10)-(3.11)]. In particular, we have and for some a > 0 the estimates hold with some C γ,δ > 0, see (2.7). Note that the solutions to the problem (2.15) and, more generally, to (2.16) can be considered in a generalized sense for less regular initial data u 0 than from X α if the linear semigroup satisfies (2.17) in some scale of Banach spaces. For details of this approach we refer the reader to [40] for the linear equations and to [8,38] for the nonlinear problems with applications. In the latter paper there was introduced and thoroughly discussed the notion of a γ -solution to (2.16), (2.17), hence a generalization of the notion of a mild solution to the semilinear problem (2.15), requiring that the nonlinearity fulfills and there exist ρ ≥ 1 and L > 0 such that Under this particular form of local Lipschitz continuity of F more accurate estimates for solutions of the corresponding integral equation (2.16) are possible without assuming analyticity nor continuity of the linear semigroup T (t) at t = 0, that is, for The following definition of the γ -solution is used in that reference.
Using the smoothing action of the integral equation, the authors are able to set initial data in a larger space X γ on the scale than X α . More precisely, they may have u 0 ∈ X γ with the parameter γ satisfying

Domains of Fractional Powers of Negative Dirichlet Laplacian
We next recall the description of the domains of fractional powers  In case of sectorial operators in the Lebesgue space L p ( ), 1 < p < ∞, such description can be found in [24, Theorem 1.6.1] when ∂ ∈ C 2 α , or in more detail in [45, pp. 81, 560]. In particular, for uniformly strongly elliptic operators in the L p ( ) setting under the assumptions of [45,Theorem 16.14] (which are satisfied for the negative Laplacian in part regarding the boundedness of imaginary powers thanks to [37, p. 166]) we have D((− L p ) α ) ⊆ W 2α, p ( ), α ≥ 0, and, in particular, ]. However, those results were formulated under the assumption of C ∞ regularity of the domain and the coefficients of the differential operator involved. A lot of effort was needed to weaken those too high regularity assumptions (see e.g. [13,17,34,37]). An extended discussion of the domains of fractional powers of sectorial operators corresponding to second order uniformly strongly elliptic operators in bounded smooth domains ⊂ R N with Dirichlet or Neumann boundary conditions was also reported in [29] in case of the Hilbert scales.

Remark 2.27
Considering higher fractional powers of the negative Dirichlet Laplacian, we deal with the so called Navier boundary conditions. More precisely, the functions u from the domains of higher powers (− ) k , k = 1, 2, . . ., considered in a bounded domain with ∂ ∈ C 2k , satisfy the Navier boundary conditions, see [35], given by The situation presented in [43,Theorem 5.4.4/2] for a regular elliptic operator A of 2mth order and m boundary conditions B 1 , . . . , B m , and realized as a sectorial positive operator A in the space L p ( ), 1 < p < ∞, is, however, slightly different. Besides the (convenient) C ∞ regularity assumption on ∂ , the result formulated in [43] states that the operator acts as a bijective bounded operator from W 2m+s p,{B j } ( ) onto W s p ( ) (with notation of [43]). However, the functions from the fractional power space D(A α ) for α ∈ N require as many as αm boundary conditions (see [6,Remark 1.3.7]). Additionally, the operator A takes isometrically D(A α+1 ) onto D(A α ) (see (2.10)), making it different from the action of (2.20).

Global Extendibility of Small Data Solutions for Super-Linear Nonlinearity
Let the realization of the linear differential operator A in (2.1) be a selfadjoint and positive definite operator A in a Hilbert space X , so that the (one-sided or extrapolated) fractional power scale X α , α ∈ J , consists of Hilbert spaces normed by · α given in (2.8) or (2.11), cp. also (2.13). For more details of the Hilbert fractional power scale we refer the reader e.g. to [29], [ Assume that the nonlinearity F in (2.1) defines the Nemitskii operator F : X α → X β with some α, β ∈ J , 0 ≤ α − β < 1, which is Lipschitz continuous on bounded subsets of X α . We consider the semilinear problem (2.15) with the first equation satisfied in the base space X β . For each u 0 ∈ X α there exists a unique (X α , X β ) solution of (2.15), defined on the maximal interval of existence and satisfying the first equation in (2.15) in X β and the initial condition in (2.15) in X α . In particular, it satisfies the Duhamel formula (2.16) with T (t) = e −At or T (t) = e −A −1 t being an analytic semigroup for which the estimates (2.17) hold.
We further assume that the following super-linear growth condition holds where a, b ≥ 0, ν ∈ [0, 1) and μ + ν > 1. We next formulate a result on the global in time a priori estimate in X β+ 1 2 valid for small data solutions (see [20,Sect. 6.2] for a prototype of this estimate for the 3D Navier-Stokes equation).

Lemma 2.28
Let condition (2.21) be satisfied with b ≥ 0 such that (2.24) holds and let u be an (X α , X β ) solution of the problem (2.15) with small initial data u 0 ∈ X β+ 1 2 , as specified in (2.25). Then such solutions are bounded a priori in X β+ 1 2 norm.
Proof We carry out the proof for the extrapolated fractional power scale. Under the conditions of the theorem, the extrapolated operator A −1 is selfadjoint and positive definite in the Hilbert space X −1 , see [1, Theorem V.1.5.15]. We apply the operator (A −1 ) β+1 : X β → X −1 to the first equation in (2.15) and take the scalar product ·, · in X −1 with (A −1 ) β+2 u, getting Using the fact that A −1 is selfadjoint and the growth condition (2.21), we obtain Applying to the penultimate term the Young inequality with ε = 1 4 , and to the last term of (2.22) the Cauchy inequality, it extends to with a certain constant C = C(a, μ, ν) > 0.
Next, using the generalized Poincaré inequality (see e.g. [20, p. 115]) where λ 1 > 0 stands for the first positive eigenvalue of the operator A −1 , we find that Analyzing the real function in the right-hand side of the above differential inequality, where y corresponds to u 2 and the minimal value of g b is attained at the point y min = λ 1 we get g b (y min ) < 0, with the graph of g b having two nonnegative zeros; 0 ≤ y 1 < y 2 (note that y 1 = 0 when b = 0). Consequently, for such a small b 2 , it follows from In particular, this lemma guarantees the global in time extendibility of the local X β+ 1 2 , X β solutions of (2.15) corresponding to small initial data u 0 satisfying (2.25) whenever the growth condition (2.21) is satisfied with b ≥ 0 sufficiently small [see (2.24)].
We also recall [15, Theorem 4.2] suitably modified for our purposes.

Theorem 2.29
Assume that the following a priori estimate for the (X α , X β ) solution u(t) of (2.15) satisfying u(0) = u 0 ∈ X α holds in a normed space Y such that Y → X α , that is, there exists a function c : [0, T ) → [0, ∞), 0 < T ≤ ∞, bounded on compact intervals and such that

26)
where τ u 0 denotes the life time of the solution. Furthermore, assume that the following subordination condition holds for the nonlinearity, that is, there exist a nondecreasing function g : [0, ∞) → [0, ∞) and a constant θ ∈ [0, 1) such that Then we have τ u 0 ≥ T .
If T = ∞ in the a priori estimate (2.26), then the solution of (2.15) from Theorem 2.29 exists globally in time. Moreover, if T = ∞ in (2.26) and the function c(t) is bounded on [0, ∞) by some constant, then this solution of (2.15) exists globally in time and is bounded.
As will be seen from the analysis of the examples in Sects. 3 and 4, for many nonlinearities in the physical models (in particular, in the 2D quasi-geostrophic equation) there is a limitation from above of the level on the fractional power scale at which the nonlinearity will be realized. This is connected with the required compatibility conditions to be satisfied, see [28] for considerations concerning this issue. Consequently, for such equations we are not able to construct 'too regular' solutions in the formalism of this paper. But there is also another limitation, from below, for the admissible 'too weak' regularity of solutions, connected with fact that the nonlinearity will not behave well on such less regular solutions. More precisely, the Nemitskii superposition operator (see e.g. [2] or [9, Sect. 2.5]) associated with the nonlinearity will not be well-defined.
In general, it is not easy to relate the nonlinearity with the fractional power scale generated by the linear operator in the equation. However, for example, the semilinear parabolic Dirichlet problem considered in a bounded regular domain ⊂ R N fits well into the scale generated by − because the gradient ∇u is of the same order as (− ) 1 2 u and it is not hard to formulate the admissible growth restriction on the nonlinearity leading to global in time solutions of (2.28). In particular, we can assume that where the admissible growth exponent 2 − ε follows from the moment inequality In low space dimensions N ≤ 3, having the embedding H 2−ε ( ) → L ∞ ( ), it is sufficient to take g : R → R locally Lipschitz continuous. For N ≥ 4, we need to limit eventually the growth of g to stay inside the scale generated by − on L 2 ( ). Alternatively, we can switch into the scale generated by − realized in L p ( ), p > N , to avoid implementing additional growth restrictions on g.

Dirichlet Problem for 2D Quasi-geostrophic Equation
This and the following section are devoted to a discussion of two examples of semilinear equations, which can be solved at various levels of the fractional power scale generated by the positive sectorial operator in their main part. Here we will concentrate on the 2D quasi-geostrophic equation with Dirichlet boundary conditions, whereas in the next section we discuss the fractional chemotaxis equation, with recent results reported in [3,19,20,44]. The approach to these problems is unified and can be used to treat other, even more complicated, problems as well. For instance, in separate papers [7,18] the celebrated 3D Navier-Stokes equation on the extrapolated fractional power scale was discussed in detail. Let us also mention that studying admissible levels of the extrapolated fractional power scale at which the examples can be settled, we use a similar reasoning to the one originated in [22,Lemma 2.2] in case of the Navier-Stokes equation.
In this section we consider the Dirichlet problem in a bounded C 2 domain for the 2D quasi-geostrophic equation (see e.g. [11,19,20]). Frequently studied earlier Cauchy problem in R 2 for that equation, with a positive parameter ω, has the form where the velocity field u = (u 1 , u 2 ) is expressed by the stream function ψ through the relation Here θ represents the potential temperature, κ > 0 is a diffusivity coefficient, f is a free time-independent force, and (− ) ω denotes a nonlocal operator which is an equivalent way to define ω-power of the negative Laplacian in R N , see [27]. Following the studies of e.g. [10,19], we replace in the above − and (− ) ω by the negative Dirichlet Laplacian in L 2 ( ) and its fractional power, respectively, and consider the counterpart of (3.1), (3.2) in the form of the Dirichlet problem for the quasi-geostrophic equation in a bounded domain in R 2 with C 2 boundary. Observe that the Formula (3.2) can be alternatively expressed using the components of the Riesz transform R in the bounded domain in the form

Rescaling the time variable, renaming f and plugging (3.3) into the equation leads to the abstract Cauchy problem, which we consider in this section
with a parameter ω ∈ 1 2 , 1 . Before addressing the question of solvability of (3.4), we recall some basic properties of the Riesz transform. By an analogy to the case of the whole of R N , see e.g. [33,Sect. 12.3,(12.22)], the Riesz transform in a bounded regular domain ⊂ R N is defined through the formula for elements from the domain of σ -power of the negative Dirichlet Laplacian D((− ) σ )), see Sect. 2.5.

Proposition 3.1 The Riesz transforms
where the equivalence of norms in D((− ) s ) and H 2s ( ) was used.
We also note that the Riesz transform is bounded in L p ( ), 1 < p < ∞; see [41,Theorem C] in case of a bounded C 1 domain .

Proposition 3.2 If is a bounded C
Proposition 3.1 shows that the Riesz operators R j switch from the scale D((− ) σ ) into the scale H 2σ ( ), σ ≥ 0. Indeed, they do not preserve the boundary condition (see also Remark 3.3). Thus, if we want to consider (3.4) in the scale D((− ) σ ), there exists a maximal admissible regularity of the considered solutions of (3.4) and we are allowed to choose as the base space the space D((− ) σ ) as far as it coincides with H 2σ ( ), that is, whenever σ < 1 4 , see (2.18). [10] that for regular enough solutions, for example if θ ∈ H 1 0 ( ), the normal component of the velocity vanishes at the boundary R ⊥ θ ·ν | ∂ = 0, because the stream function ψ = −(− ) − 1 2 θ vanishes at the boundary and its gradient is normal to the boundary. In case of a bounded C 2 domain we cannot thus expect that the Riesz transforms R j will vanish at the boundary of .

Remark 3.3 It was observed in
We also remark that if = R N the Riesz transforms R j in the whole of R N commute with the fractional powers of the negative Laplacian, for the proof see [33, (12.24)]. Consequently, we have the commutativity of the fractional powers with partial derivatives in R N . However, there is no similar property for the Riesz operator R in a bounded domain.
Our next goal is to extend the just constructed local (X α , X 0 ) solutions with α ∈ 1 2ω , 1 globally in time. For that purpose, we will use a version of a subordination condition. Indeed, we estimate applying the Hölder inequality with ε > 0 and Proposition 3.2 in which by the Sobolev embedding gives the following subordination condition: The uniform in ω a priori estimate is available for solutions of (3.5), see e.g. [19, (24)]. This a priori estimate in L ∞ ( ) together with the subordination condition (3.8) is sufficient for the global in time extendibility of the local solutions constructed in Theorem 3.4. Proof Applying the moment inequality (2.6) to the subordination condition (3.8), we get where g is a nondecreasing function. Combining this with the a priori estimate (3.9) in L ∞ ( ), we apply Theorem 2.29 and conclude that τ = ∞ and the solution θ is bounded in X α .

Solutions of (3.4) in a Smoother Base Space than L 2 (Ä)
As it was already mentioned, the Riesz operators do not preserve the boundary condition, so if we want to consider the 2D quasi-geostrophic equation in the base space X β with β ≥ 0 from the one-sided fractional power scale generated by A = (− ) ω , our choice is limited to spaces for which X β = D((− ) ωβ ) coincides with H 2ωβ ( ) up to the equivalence of norms, that is, whenever 0 ≤ β < 1 4ω , see (2.18). The operator A in (3.5) is now understood as the realization of (− ) ω in X β . Note that, due to the relation u = R ⊥ θ , the following equality holds (3.10) A procedure, similar to the one that led to Theorem 3.4, allows us to find a local in time solution to (3.5) varying in the phase space X α → D (− ) ωβ+ 1 2 with α such that α ∈ β + 1 2ω , β + 1 and α > 1 2ω . Note that if β > 0 then we can take α = β + 1 2ω . The key point is the Lipschitz continuity of the nonlinearity F on bounded sets of where we used Proposition 3.1, the property that ∇ is a bounded linear operator from H 1+2ωβ ( ) into H 2ωβ ( ) and the fact that H 2ωα ( ) is a Banach algebra, since 2ωα > 1, see [45, (1.89)]. Consequently, given f ∈ X β , we have obtained Theorem 3.6 Let 0 ≤ β < 1 4ω , f ∈ X β and α ∈ β + 1 2ω , β + 1 be such that α > 1 2ω . Then there exists a unique local in time (X α , X β ) solution θ of (3.5) such that where γ < ω(β + 1).
Since for β > 0 the available a priori estimate in L ∞ ( ) is not sufficiently strong to guarantee global solvability of the problem, we will establish below an a priori estimate in H 2ω β+ 1 2 ( ) for small data via Lemma 2.28. This will be carried out only for the values of parameter ω in (3.4) from the interval 2 3 , 1 . An appropriate subordination condition will guarantee that these solutions exist globally in time.

Proof
The (X α , X β ) solution θ of (3.5) defined on the maximal interval of existence was constructed in Theorem 3.6. Note that α > β + 1 2 . Since ω ∈ 2 3 , 1 , we can take r = β + 1 in Lemma 3.7 and get with some ε ∈ (0, 1). Thus by Lemma 2.28 the solution is bounded in the norm of X β+ 1 2 provided that the data is small. Applying Lemma 3.7 with r = α, we obtain the following subordination condition which by Theorem 2.29 implies that the solution θ exists globally in time.

Solutions of (3.4) in a Base Space Larger than L 2 (Ä)
We begin with a simple observation helpful in deriving a range of parameters in the estimates made in this subsection and later in Sect. 4.3.
To examine the global existence in time of these solutions, we observe that a very similar reasoning as above leads to a suitable estimate of the bi-linear form associated with the nonlinearity, which involves the L ∞ ( ) norm.
This allows us to extend the solutions globally in time.

Fractional Chemotaxis System
Our second example is a nonlocal model of chemotaxis, which has recently been popularized by [3] and [44]. The considered problem is a generalization of the known Keller-Segel model of chemotaxis introduced in 1970 to describe the self-organization of a number of biological systems. In [4,21] a motivation and mathematical background of the fractional version of that model was described and in [45, p. 243] a semigroup approach, similar to ours, to the original chemotaxis system was used. Due to the biological interpretation its solution u(t, x) expressed density of the biological individuals (most often bacteria) at time t and location x and v(t, x) denoted density of the chemical attractant produced to incline the colony to aggregate. Required by that interpretation, non-negativity of solutions was discussed e.g. in [45, p. 246] or [3,Sect. 2.4]. Concentrating here on mathematical aspects of the local and global solvability of (4.1), we consider solutions eventually changing their sign; non-negativity is thus an additional property required for solutions appearing in the biological application of that model. We will not study that property here.
Although typically considered as a Cauchy problem in R N , we study it here as a Dirichlet problem in a bounded C 2 domain in the form with a parameter ω ∈ 1 2 , 1 . Though at first sight this looks like a system of a parabolic and an elliptic equation, we will rewrite it as a single semilinear parabolic equation. Indeed, noting that ∇ · (u∇v) = ∇u · ∇v + u v and deriving (formally) v from the second equation of (4.1) as v = (− ) −1 u, we insert it into the first equation of (4.1) to obtain the initial boundary value problem In this section we consider (4.2) as an abstract Cauchy problem with A = (− ) ω , ω ∈ 1 2 , 1 , being a positive sectorial operator in L p ( ) with a chosen p and where − denotes the Dirichlet Laplacian realized in L p ( ). For the description of the domains of its fractional powers (− ) σ , σ > 0, see (2.19). Note also that, by (2.19), if u is in D((− ) σ ), σ > 0, then v = (− ) −1 u satisfies zero Dirichlet boundary condition, which justifies the reformulation of the problem (4.1).

Solutions of (4.3) in L p (Ä) as the Base Space
We choose p > N and consider (4.3) in the base space X = L p ( ). We estimate the terms appearing in the nonlinearity F for any ϕ, ψ ∈ W 1, p 0 ( ) as follows (4.5) These estimates show that both above terms are bi-linear operators from W Hence the nonlinearity F is Lipschitz continuous on bounded sets as an operator between X α = D((− ) ωα ) → W 1, p 0 ( ) with α ≥ 1 2ω and X = L p ( ). Indeed, given u 1 , u 2 ∈ B with B bounded in X α , α ≥ 1 2ω , by (4.4) and (4.5) we have Consequently, we get local solvability of (4.3) in the phase space X α for 1 2ω ≤ α < 1.

Theorem 4.1
Assume that ω ∈ 1 2 , 1 , ∂ ∈ C 2 and u 0 ∈ X α with 1 2ω ≤ α < 1. Then there exists a unique local (X α , X 0 ) solution to the problem (4.3) having the following regularity properties: where γ < ω and τ > 0 stands for the maximal time of existence of that solution. Moreover, the corresponding Duhamel formula is satisfied in X .

Remark 4.2
Considering the problem in dimension N = 3, we can stay within the fractional scale of Hilbert spaces defined by A = (− ) ω , where − is the negative Laplacian realized in the base space X = L 2 ( ). Indeed, this observation is a simple consequence of − being an isometry between D (− ) 3 2 and D (− ) 1 2 , see (2.9), and the following estimates valid for ϕ, ψ ∈ H 1 0 ( ) = D (− ) Hence the nonlinear terms give rise to bi-linear transformations from H 1 0 ( )× H 1 0 ( ) into L 2 ( ) and F satisfies in this case for u 1 , u 2 ∈ H 1 0 ( ) Consequently, if N = 3 we again obtain local solvability of the problem (4.3) in the phase space X α = D((− ) ωα ) with 1 2ω ≤ α < 1, which is a Hilbert space embedded into H 1 0 ( ), and Theorem 4.1 holds true if X = L 2 ( ) and N = 3.
In case of the Cauchy problem in R N for chemotaxis equation a singular steady state or a generalized Chandrasekhar type solution u C (x) exists, see [3, p. 32]. It is further shown there that to a nonnegative, radial and sufficiently regular initial data staying below u C (x) correspond global in time solutions, cp. [3, Theorem 2.3.1]. Even stronger property is known to hold in 1D case for solutions of the fractional Keller-Segel system on the real line with ω ∈ 1 2 , 1 ; all solutions are global in time in that case (see [21] for a more complete analysis).
Here, in case of the Dirichlet boundary condition, we will show that local (X α , X 0 ) solutions with α ∈ 1 2ω , 1 , introduced in Remark 4.2 for N = 3, will be extended globally in time if their initial data come from a certain neighbourhood of zero in the space X 1 2 = D (− ) ω 2 . To this end, observe that for u ∈ D((− ) ω ) = X 1 we have where the following consequences of the Nirenberg-Gagliardo inequality were used Thus we obtain which means that the assumption (2.21) is satisfied with β = 0 and b = 0. Therefore, it follows from Lemma 2.28 that for α ∈ 1 2ω , 1 the (X α , X 0 ) solutions of (4.3) originating in a sufficiently small ball around zero in X see (4.6), and a simplified argument of the proof of Lemma 2.28. Then the (X α , X 0 ) solutions with α ∈ 1 2 , 1 , which correspond to small initial data satisfying where c is the constant at the right hand side of (4.7) and λ 1 denotes the first eigenvalue of − , are bounded a priori in H 1 0 ( ). Consequently, they exist globally in time. In order to prove global extendibility in time of (X α , X 0 ) solutions of (4.3) for α ∈ 1 2ω , 1 and small initial data, it suffices to show a subordination condition involving the X 1 2 norm in which these solutions are bounded a priori. Note that if α ∈ 1 2ω , 1 then for u ∈ X α ∇u · ∇((− ) −1 u) , with any 0 < δ < 2ωα−1 ω(2α−1) ≤ 1, since by the Nirenberg-Gagliardo inequality we have Moreover, we also have , since by the Nirenberg-Gagliardo inequality we get with δ as above. Concluding, we obtain which is a form of a subordination condition as in (2.27). Hence by Theorem 2.29 the local (X α , X 0 ) solutions of the fractional chemotaxis equation (4.3) exist globally in time for initial data with small X 1 2 norm. Note that for α = 1 2ω we have the estimate which combined with the boundedness of the X 1 2 norm of solutions with small initial data implies that F is sublinear in this case. The Duhamel formula and the Volterra type inequality guarantee that these local solutions can be extended globally in time.

More Regular Local Solutions of (4.3)
If − is realized in the space L p ( ) with p > N , we are also able to construct more regular local solutions to the fractional chemotaxis problem (4.3) by considering it in the base space X β = D((− ) ωβ ) = W 2ωβ, p ( ) with β ∈ 0, 1 2ω p , see (2.19). For this purpose, we will estimate the bi-linear form ∇ϕ · ∇((− ) −1 ψ) in the D((− ) σ ) norm for 0 ≤ σ < 1 2 p , that is, in the case when the fractional power space does not involve the boundary condition. To shorten some notations, we set g := ∇(− ) −1 ψ and h := ∇ϕ. We need to estimate g · h D((− ) σ ) when 0 ≤ σ < 1 and estimate the W 1, p ( ) norm, that is, We have whereas the second term is estimated as follows: We also get Consequently, the nonlinearity defines a bi-linear form from W 2, p ( )× W 2, p ( ) into D((− ) σ ), 0 ≤ σ < 1 2 p . Therefore, we obtain the following local solvability result.