LOCAL CRITERIA FOR BLOWUP IN TWO-DIMENSIONAL CHEMOTAXIS MODELS PIOTR BILER, TOMASZ CIEŚLAK, GRZEGORZ KARCH, AND JACEK ZIENKIEWICZ arxiv:4.787v3 [math.ap] 5 Apr 26 Abstract. We consider two-dimensional versions of the Keller Segel model for the chemotaxis with either classical Brownian) or fractional anomalous) diffusion. Criteria for blowup of solutions in terms of suitable Morrey spaces norms are derived. Moreover, the impact of the consumption term on the global-in-time existence of solutions is analyzed for the classical Keller Segel system.. Introduction We consider in this paper the following version of the parabolic-elliptic Keller Segel model of chemotaxis in two space dimensions.).2) u t + ) α/2 u+ u v) =, x R 2, t >, supplemented with the initial condition v γv +u =, x R 2, t >,.3) ux,) = u x). Here the unknown variables u = ux,t) and v = vx,t) denote the density of the population and the density of the chemical secreted by the microorganisms, respectively, and the given consumption or degradation) rate of the chemical is denoted by γ. The diffusion operator is described either by the usual Laplacian α = 2) or by a fractional power of the Laplacian ) α/2 with α,2). The initial data are nonnegative functions u L R 2 ) of the total mass.4) M = u x)dx. Our main results include criteria for blowup of nonnegative solutions of problem.).3) expressed in terms of a local concentration of data Theorem 2.), and the existence of global-in-time solutions for the initial condition of an arbitrary mass M and each sufficiently large γ Theorem 2.6). The novelty of these blowup results consists in using local properties of solutions instead of a comparison of the total mass and moments of a solution as was done in e.g. [7], [], [4, 5, 6], [6], [4], [], and [3]. In particular, we complement the result in [] saying that solutions of.).3) with α = 2, fixed γ and sufficiently well concentrated u with M > 8π blow up in a finite time, by showing Date: October 9, 28. 2 Mathematics Subject Classification. 35Q92, 35B44, 35K55. Key words and phrases. chemotaxis, blowup of solutions, global existence of solutions. This work is accepted for publication in DCDS-A.
2 P. BILER, T. CIEŚLAK, G. KARCH, AND J. ZIENKIEWICZ that solutions of that system with u of arbitrary M > and all sufficiently large γ are global-in-time. Many previous works have dealt with the existence of global-in-time solutions with small data in critical Morrey spaces, i.e. those which are scale-invariant under a natural scaling of the chemotaxis model, cf. e.g. [] and [3]. Our criteria for a blowup of solutions with large concentration can be expressed by Morrey space norms see Remark 2.3 below for more details), and we have found that the size of such a norm is critical for the global-in-time existence versus finite time blowup. The analogous question for radially symmetric solutions of the d-dimensional Keller Segel model with d 3 has been recently studied in [5]. 2. Statement of results It is well-known that problem.).3) with α = 2 has a unique mild solution u C[,T);L R 2 )) for every u L R 2 ) and γ. Here, as usual, a mild solution satisfies a suitable integral formulation 5.) of the Cauchy problem.).3) as recalled at the beginning of Section 5. Moreover, given u = ux,t), we define v = + γ) u, see Lemma 3., below. Due to a parabolic regularization effect following e.g. [9, Th. 4.2]), this solution is smooth for t >, hence, it satisfies the Cauchy problem in the classical sense. Moreover, it conserves the total mass.4) ux,t)dx = u x)dx for all t [,T), R 2 R 2 and is nonnegative when u. Proofs of these classical results can be found e.g. in [, 3, 2,, 8], see also Section 5 of this work. Analogous results on local-in-time solutions to the Cauchy problem.).3) with α,2) have been obtained in [3], [3, Th. 2]. To the best of our knowledge, [6, Th..] and a recent [2, Th., Th. 2] are the only results on local-in-time classical solutions of the Cauchy problem.).3) with α,], d 2. Thus, the case iii) of Theorem 2. asserts that such a solution cannot be global-in-time for initial data satisfying 2.2). In our first result, we formulate new sufficient conditions for blowup i.e. nonexistence for all t > ) of such local-in-time solutions of problem.).2). Theorem 2.. Consider u C[,T);L R 2 )) a local-in-time nonnegative classical solution of problem.).3) with a nonnegative u L R 2 ). i) If α = 2, γ = the scaling invariant Keller Segel model), then for each M > 8π the solution u blows up in a finite time. ii) Let α = 2 and γ > the Keller Segel model with the consumption). If M > 8π and if u is well concentrated around a point x R 2, namely, there exists R > such that 2.) e γr u y)dy > 8π and u y)dy < ν { y x <R} { y x R}
2D CHEMOTAXIS MODELS 3 with an explicitly computed small constant ν >, then the solution u blows up in a finite time. iii) Let α,2) and γ the Keller Segel model with fractional diffusion). If there exist x R 2 and R > such that 2.2) R α 2 u y)dy > C and u y)dy < ν, { y x <R} { y x R} for some explicit constants: small ν > and big C >, then the solution u ceases to exists in a finite time. Remark 2.2. The result i) for α = 2 and γ = is, of course, well known, but the proof below slightly differs from the previous ones. The case ii) α = 2 and γ > has been considered in [] but the sufficient conditions for blowup were expressed in terms of globally defined quantities: i.e. mass M > 8π and the moment u x) x 2 dx. Remark 2.3. The case iii) α < 2. Recall that the homogeneous) Morrey space M p R 2 ) is defined as the space of locally integrable functions such that u M p = sup R 2/p ) uy)dy <. R>,x R 2 { y x <R} The first condition in 2.2) is equivalent to a sufficiently large Morrey norm of u in the space M 2/α R 2 ). Indeed, obviously we have u M 2/α R α 2 u y)dy { y x <R} for every x and R >, but also there is x R 2 and R > such that u M 2/α 2R α 2 u y)dy. { y x <R} Thus, our blowup condition in terms of the Morrey norm seems to be new and complementary to that guaranteeing the global-in-time existence of solutions, where smallness of initial conditions in the M 2/α -Morrey norm has to be imposed, cf. prototypes of such results in [, Theorem ] and [3, Remark 2.7]. Remark 2.4. A natural scaling for system.).2) with γ = : u λ x,t) = λ α uλx,λ α t), leads to the equality u λ dx = λ α 2 udx, i.e. mass of rescaled solution u λ can be chosen arbitrarily with a suitable λ >. Thus, the conditions in Theorem 2..iii are insensitive to the actual value of M, so w.l.o.g. we may suppose that M =.
4 P. BILER, T. CIEŚLAK, G. KARCH, AND J. ZIENKIEWICZ Remark 2.5. The second parts of the condition 2.) and 2.2) are not scaling invariant. However, we believe that these assumptions are not necessary for the conclusion in Theorem 2..ii, 2..iii. In fact, one can prove it for α close to 2 by an inspection of methods in [5, 7]. Next, we show that the first condition in the concentration assumptions 2.) is in some sense optimal to obtain a blowup of solutions. In the following theorem we show that for every initial integrable function u, even with its L -norm above 8π, the corresponding mild solution to the model.).3) with α = 2 is global-in-time for all sufficiently large consumption rates γ >. Theorem 2.6. Let α = 2, γ >. For each u L R 2 ), there exists γu ) > such that for all γ γu ) the Cauchy problem.).3) has a global-in-time mild solution satisfying u C[, );L R 2 )). This is a classical solution of system.).2) for t >, and satisfies for each p [, ) the decay estimates 2.3) supt /p ut) p <. t> Thus, for each u not necessarily nonnegative) and γ large enough depending on u L R 2 ), solutions of the Cauchy problem are global-in-time, so there is no critical value of mass which leads to a blowup of solutions. On the other hand, if M > 8π, then for γ the solutions blow up in a finite time, as it is seen from the sufficient conditions for blowup in Theorem 2..ii. 3. Notation and preliminaries In the sequel, q denotes the usual L q R 2 ) norm, and C s are generic constants independent of t, u, z,... which may, however, vary from line to line. Integrals with no integration limits are meant to be calculated over the whole plane. Let us denote by G the Gauss Weierstrass kernel of the heat semigroup e t on L p R 2 ) space 3.) Gx,t) = 4πt) exp ) x 2. 4t As it is well known the convolution with G, denoted by Gt) z = e t z, satisfies the following L q L p estimates 3.2) e t z p Ct /p /q z q and 3.3) e t z p Ct /2+/p /q z q for all q p, t >. Moreover, for each p > and z L R 2 ) the following relation holds 3.4) lim t t /p e t z p =
which is, e.g., noted in [9, Lemma 4.4]. 2D CHEMOTAXIS MODELS 5 Lemma 3.. For every γ >, the operator +γ) solving the Helmholtz equation.2) satisfies 3.5) +γ) z q Cγ /p /q /2 z p, for every p < 2 < q < such that < and some C independent of γ. In the p q 2 critical case = inequality 3.5) also holds provided p >. Moreover, the Bessel p q 2 kernel K γ of +γ) has the following pointwise behavior at and 3.6) 3.7) and satisfies the global one-sided bound K γ x) x as x, 2π x 2 K γ x) C γ x x e as x, 3.8) x K γ x) γ x 2π e. Proof. The proof of inequality 3.5) requires separate arguments in two cases, < p q 2 and =. In the first case, the result is a consequence of inequalities 3.2) and 3.3) p q 2 by representing the operator +γ) as the Laplace transform 3.9) +γ) = e γs e s ds. Indeed, we have the following representation of K γ in the Fourier variables so that K γ z)ξ) = +γ) z q C ξ 2 +γẑξ) = e γs e s ξ 2 ẑξ)ds, e γs s /q /p /2 ds z p Cγ /p /q /2 e s s /2+/q /p ds z p, the latter integral is finite due to the assumption on p and q. When p q = 2, inequality 3.5) follows from the end-point case of the Sobolev inequality ) u q C u p. For properties 3.6), 3.7) and 3.8), see e.g. [, Lemma 3.] and [9, Ch. V, Sec. 6.5].
6 P. BILER, T. CIEŚLAK, G. KARCH, AND J. ZIENKIEWICZ Remark 3.2. Let us notethat thereference [, Theorem 2.9] provides us with precise conditions on radial convolution kernels K leading to a blowup of solutions of general diffusive aggregationequationswiththebrowniandiffusionoftheformu t u+ u K u)) =. They are strongly singular, i.e. they have the singularity at : limsup x x Kx) <, and are of moderate growth at : x Kx) C x 2. Of course, the Bessel kernel K γ is strongly singular in the sense of [], as it is seen from 3.8). 4. Blowup of solutions In this section we prove Theorem 2. using the method of truncated moments which is reminiscent of that in the papers [8], [2]. First, we define the bump function ψ and its rescalings for R > 4.) ψx) = x 2 ) 2 + = { x 2 ) 2 for x <, for x, The function ψ is piecewise C 2 R 2 ), with suppψ = { x }, and satisfies ) x ψ R x) = ψ. R 4.2) 4.3) ψx) = 4x x 2 ) for x <, ψx) = 8+6 x 2 ) 8ψx) 8 for x <. We will use in the sequel the fact that ψ is strictly concave in a neighbourhood of x =. ) Lemma 4.. For each ε, 3, the function ψ defined in 4.) is strictly concave for all x ε. More precisely, ψ satisfies 4.4) Hψ θε)i for all x ε, where Hψ is the Hessian matrix of second derivatives of ψ, θε) = 4 3ε 2 ), and I is the identity matrix. In particular, we have 4.5) θε) ր 4 as ε ց. Proof. For every ξ R 2 the following identity holds ξ Hψξ = 4 ξ 2 x 2) +2x ξ) 2). Thus, by the Schwarz inequality, we have ξ Hψξ 4 ξ 2 3 x 2 ). Next, we recall a well-known property of concave functions. Lemma 4.2. For every function Ψ : R 2 R which is strictly concave on a domain Ω R 2 we have for all x,y R 2 4.6) x y) Ψx) Ψy)) θ x y 2, where θ > is the constant of strict concavity of Ψ on Ω, i.e. satisfying HΨ θi. Proof. By the concavity, we obtain Ψx) Ψy)+ Ψy) x y) θ 2! x y 2.
2D CHEMOTAXIS MODELS 7 Summing this inequality with its symmetrized version with x, y interchanged) leads to the claim. We have the following scaling property of the fractional Laplacian 4.7) ) α/2 ψ R x) = R α ) α/2 ψ ) R, and we notice the following boundedness property of ) α/2 ψ. Lemma 4.3. For every α,2] there exists a constant k α > such that 4.8) ) α/2 ψx) k α. Moreover, ) α/2 ψx) for x. In particular, for α = 2 we have k 2 = 8. Proof. For α = 2, this is an obvious consequence of the explicit form of ψ, hence we assume α,2). To show estimate 4.8) for α,2), it suffices to use the following well-known representation of the fractional Laplacian with α,2) ψx+y) ψx) ) α/2 ψx) = c α P.V. dy y 2+α for certain explicit constant c α >. Now, using the Taylor formula together with the fact that ψ, D 2 ψ L R 2 ), we immediately obtain that the integral on the right-hand side is finite and uniformly bounded in x R 2. Since ψx) and ψx) = for x we have for x. ψx+y) ) α/2 ψx) = c α P.V. dy y 2+α Now, we formulate a crucial inequality in our proof of the blowup result. Lemma 4.4. For the Bessel kernel K γ with γ and a strictly concave function Ψ we have for all x,y on the domain of the strict concavity of Ψ 4.9) K γ x y) Ψx) Ψy)) θ 2π g γ x y ), where θ is the constant of the strict concavity of Ψ introduced in Lemma 4.2, and g γ is a radially symmetric continuous function, such that 4.) K γ x) = x 2π x 2g γ x ). In particular, g γ ) =, the profile of g γ decreases, and g γ x ) Ce γ x. Proof. Combining Lemma 4.2 with equation 4.) and properties 3.6), 3.7) and 3.8) we arrive immediately at the claimed formula. We are in a position to prove our main blowup result.
8 P. BILER, T. CIEŚLAK, G. KARCH, AND J. ZIENKIEWICZ Proof of Theorem 2.. We consider the quantity w R t) = ux,t)ψ R x)dx, a local moment of u.,t), where ψ R x) is defined in 4.) for each R >. Let 4.) M R t) ux,t)dx w R t) { x <R} denote mass of the distribution u contained in the ball { x < R} at the moment t. Now, using equation.) we determine the evolution of w R t) d dt w Rt) = ) α/2 ux,t)ψ R x)dx+ ux,t) vx,t) ψ R x)dx 4.2) = ux,t) ) α/2 ψ R x)dx + ux,t)uy,t) K γ x y) ψ R x) ψ R y) ) dydx, 2 where we applied the formula v = K γ u, and the last expression follows by the symmetrization of the double integral: x y, y x. Since ux,t), by the scaling relation 4.7) and Lemma 4.3, we obtain 4.3) Now, let ε ux,t) ) α/2 ψ R x)dx R α k α { x R} ux,t)dx., 3 ). By Lemma 4., the weight function ψ R in 4.) is concave for x εr with a concavity constant θ = θε). Thus, by Lemma 4.4, we have K γ x y) ψ R x) ψ R y)) R 2θε) 2π g γ x y ) for x, y < εr. Hence, the bilinear term on the right-hand side of 4.2) satisfies 4.4) 2 ux,t)uy,t) K γ x y) ψ R x) ψ R y) ) dydx R 2θε) g γ x y )ux,t)uy,t)dydx+ 4π 2 J, where the letter J denotes the integral J = R 2 R 2 \{ x <εr} { y <εr}) { x <εr}{ y <εr} ux,t)uy,t) K γ x y) ψ R x) ψ R y))dydx.
2D CHEMOTAXIS MODELS 9 We estimate the first integral on the right-hand side of 4.4) in the following way g γ x y )ux,t)uy,t)dydx 4.5) { x <εr}{ y <εr} g γ 2εR) M R t) {εr x R} g γ 2εR)M 2 R t) 2g γ2εr)m R t) ) 2 ux,t)dx {εr x R} g γ 2εR)M R t) 2 2C ε M R t)m w R t)), ψ R x) ux, t) inf { x εr} ψr x) ) dx where C ε = inf { x εr} ψr x) )) = ε 2 ) 2). In the above inequalities we used the fact that g γ is a continuous decreasing function and g γ. Next, since we have the inclusion ) R 2 R 2 \ { x < R} { y εr} { x < εr} { y < εr} ) { x εr} { y < R} ) ) { x R} { y R} andthefactorwith ψ R vanishesontheset { x R} { y R}, weobtainimmediately the estimate J 2CR 2 ψ R y) ux, t)uy, t) ) dx dy inf { y εr} ψr y) { x <R} { y εr} 2R 2 CC ε M R t) uy,t) ψ R y) ) dy 4.6) 2R 2 CC ε M R t)m w R t)), where C = sup z K γ z) D 2 ψ. Finally, estimates 4.3) 4.6) as well as inequality 4.) applied to equation 4.2) lead to the inequalities 4.7) d dt w Rt) R α M R t) k α + θε) ) 4π Rα 2 g γ 2εR)M R t)+cε)r α 2 w R t) M) R α w R t) k α + θε) ) 4π Rα 2 g γ 2εR)w R t)+cε)r α 2 w R t) M), whenever the expression in the parentheses is nonnegative, with Cε) = 3CC ε = 3C ε 2 ) 2). Now, notice that the linear function of w R in the parentheses on the right-hand side of 4.7) is monotone increasing. Thus, if at the initial moment t = the right-hand side of 4.7) is positive, then w R t) will increase indefinitely in time. Consequently, after a moment T = O R α M )) w R thefunctionw ) R t)will become largerthanthetotalmass M. This is a contradiction with the global-in-time existence of a nonnegative solution u since it conserves mass.4).
P. BILER, T. CIEŚLAK, G. KARCH, AND J. ZIENKIEWICZ Now, let us analyze the cases when the right-hand side of inequality 4.7) is strictly positive. Case i): α = 2 and γ =. We recall that by Lemma 4.3 k 2 = 8 holds. For γ =, the Bessel potential K γ should be replaced by the fundamental solution E 2 x) of Laplacian on R 2 which satisfies E 2 x) = x θε) 4π is close to π, so that g 2π x 2 2εR) =. In view of 4.5), the quantity at the expense of taking sufficiently small ε >. Choosing ε > small enough, we get blowup in the optimal range M > 8π. Indeed, if M > 8π, then there exists ε > small, and R Rε) > sufficiently large so that w R ) is sufficiently close to M and we have 8+ π w R)+Cε)w R ) M) >. Case ii): α = 2, γ >. If M > 8π and u is sufficiently well concentrated near θε) the origin, i.e. g 4π γ2εr)w R ) > k 2 = 8 and, at the same time Cε)M w R )) is sufficiently small, then the solution u cannot be global-in-time. Case iii): In the case α < 2, the blowup occurs if for some R > the quantity u x)dx is large enough, and simultaneously u is well concentrated, i.e. R α 2 { x <R} Cε)M w R )) is small. 5. Global existence of large mass solutions In this section we prove Theorem 2.6. Our proof splits naturally into several parts. The first one is a construction of local-in-time mild solutions with initial data in L with an estimate of the existence time uniform in γ. The second step consists in proving the continuation of such a local solution to a global-in-time one satisfying the nonoptimal) decay estimate limsup t t /σ /p ut) p < for each fixed p 4/3,2) and any σ,p). Finally, we will prove a uniform global L bound sup t> ut) < as well as the optimal decay hypercontractive) estimate sup t> t /p ut) p <. First of all, the Cauchy problem.).3) is studied via the integral equation a. k. a. the Duhamel formula) 5.) ut) = e t u +Bu,u)t), whose solutions are called mild solutions of the original Cauchy problem. Here, the bilinear term B is defined as 5.2) Bu, z)t) = e t s) ) us) +γ) zs) ) ds. Then, to solve equation 5.) in a Banach space E,. ) of vector-valued functions, it is sufficient to prove the boundedness of the bilinear form B : E E E 5.3) Bu,z) η u z,
2D CHEMOTAXIS MODELS with a constant η independent of u and z. The first and the second steps toward the proof of Theorem 2.6 are based on a lemma which is convenient to formulate in the following way: Lemma 5.. If Bu,z) η u z and e t u R <, then equation 5.) has a 4η solution which is unique in the ball of radius 2R in the space E. Moreover, these solutions depend continuously on the initial data, i.e. u ũ C e t u ũ ). The detailed proof of Lemma 5. can be found in [3], []. The reasoning involves the Banach contraction theorem, the unique solution being achieved as a limit in E of the sequence of successive approximations 5.4) w t) = e t u, w n+ = w +Bw n,w n ). Step. Local-in-time solutions with the initial data in L. Lemma 5.2. For every u L R 2 ) and p 4 3,2), there exists T > independent of γ, such that equation 5.) has a solution u = ux,t) in the space 5.5) E = {u L loc,t);lp R 2 )) : sup t /p ut) p < }, <t T endowed with the norm 5.6) u sup t /p ut) p <. <t T Proof. Let = 2, so that r,2). Moreover, denote by q a number satisfying r p 2 + =. We estimate the bilinear form B for each t,t) using 3.3) and 3.5) p q r Bu,z)t) p C C t s) /2+/p /r us) +γ) zs) r ds t s) /2+/p 2/p+/2 us) p +γ) zs) q ds Cγ /2 /q+/p ) t s) /p s 2/p ) sup s /p us) p ) sup s /p) zs) p ds, <s t t /p C u z <s t with a constant C > independent of γ > and also of T > ) since 2 q + p =. The last inequality is a consequence of the fact that t s) /p s 2/p ) ds = Ct /p.
2 P. BILER, T. CIEŚLAK, G. KARCH, AND J. ZIENKIEWICZ Finally, given u L R 2 ), by 3.4) we may choose T > so small to have sup t /p e t u p < <t T 4C. Thus, the existence of a solution in E follows by an application of Lemma 5.. Lemma 5.3 The L -bound). The solution constructed in Lemma 5.2 satisfies 5.7) sup ut) <. <t T Proof. Let us take the sequence w n of approximations of u as in 5.4). By Lemma 5.2 we know that there exists a constant C such that 5.8) w n C <. Next, we take any r [ 4 3,2) such that r < p, and interpolate the Lr norm 5.9) w n r w n θ p w n θ, where θ = r p, and thus θ = r. Next, observe that for + r q p p 3.3) we have 5.) w n+ t) u +C u +C u +C u +C u +C t s) /2 wn s) +γ) w n s) ds =, by 5.4) and t s) /2 w n s) r +γ) w n s) q ds t s) /2 w n s) 2 r ds, t s) /2 w n s) 2 θ) p w n s) 2θ ds t s) /2 s 2/r ) w n s) 2θ ds, where we used inequality 3.5) for r 4 3 and 5.8). Now, for t T we take r = 4 3 q = 4 and note that t s) /2 s /2 ds = π = const. Let us define A t) = sup w s),..., A n t) = sup w n s), <s t <s t so that A t) <, and as a consequence of 5.) we arrive at A n+ t) C +C 2 A n t) <, so that where = 2θ <, since r = 4/3 and p < 2. Therefore A n t) is uniformly bounded in L,T) with a bound which depends on C and C 2 but is independent of n) by an easy
2D CHEMOTAXIS MODELS 3 recurrence argument. From the fact that w n u when n we infer that for any t > the family w n,t) converges to u,t) in L loc R2 ). Applying the Fatou lemma we see that ut) liminf n A n t) C for every fixed < t T. Step 2. Global-in-time solutions for γ large. Lemma 5.4. Now we keep p 4 3,2) and take any σ,p). For every u L σ R 2 ), there exists a constant γu ) > such that for all γ γu ) equation 5.) has a unique solution in the new functional space 5.) Ẽ = {u L loc, );L p R 2 )) : supt /σ /p ut) p < }, t> endowed with the norm 5.2) u supt /σ /p ut) p <. t> Proof. Let = 2 + for some suitable σ,p) so that r,2). Moreover, denote r p 2 σ by q a number satisfying + =. Under this choice of parameters we make sure that p q r < and q > 2, so that we can use 3.5) to estimate the bilinear form B p q 2 Bu,z)t) p C C t s) /2+/p /r us) +γ) zs) r ds t s) 2/σ /p us) p +γ) zs) q ds Cγ /2 /q+/p sup <s t t s) /σ /p s 2/p 2/σ ) ) s /σ /p us) p sup s /σ /p zs) p ds t /p /σ Cγ /σ u z. <s t Thus, we obtained inequality 5.3) with the norm defined in 5.2) and with η = Cγ /σ. We may choose γu ) so large to have supt /σ /p e t u p < γ /σ /4C t> for all γ γu ), which is possible due to estimate 3.2). Then, the proof is completed by applying Lemma 5.. Step 3. Proof of Theorem 2.6 and optimal hypercontractive estimates. By Lemma 5.2, we have a local-in-time solution on an interval,t] with T > independent of γ. Moreover, u T 2) L σ R 2 ) L L p. Thus, we may continue this local-in-time solution ut) to the whole half-line, ) choosing γ > sufficiently large. We notice that on the interval T 2,T) solutions obtained in Lemma 5.2 and Lemma 5.4
4 P. BILER, T. CIEŚLAK, G. KARCH, AND J. ZIENKIEWICZ coincide as a consequence of uniqueness assertion of Lemma 5.2. It remains to prove optimal decay estimates. The optimal L -bound. Next, we show that the solution satisfies the uniform global estimate 5.3) sup ut) <. t> For t T, similarly to the proof of Lemma 5.3, we consider a sequence of approximations of u, combine estimates 5.9) and an analog of 5.) this time on T, )) with < σ < 4 < r < p < 2 < q, θ = r p. Moreover, we define 3 p 5.4) ε = σ p, so that t /σ /p = t /p) ε). Again, we arrive at an estimate similar to 5.7) w n+ t) u + t s) /2 w n s) 2 θ) p w n s) 2θ ds, where = 2θ < since we can choose 2 < +. We split the integral on the right-hand r p side into two integrals, over the interval, 2) T and the integral over T,t). The first one 2 is estimated by CT) in view of Lemma 5.2 and 5.7). To estimate the second interval, we notice that T/2 t s) /2 w n s) 2 θ) p 5.2) t CA n t) T/2 w n s) 2θ ds t s) /2 s T 2) 2 θ)/p /σ) ds, but t s) /2 s T 2 θ)/p /σ) ds C t 2) T ) /2 2 /r) ε) 2 T/2 holds with ε as in 5.4). We notice that choosing σ > close enough to and r > 4 3 close enough to 4 we ensure that /2 2 /r) ε) <. We proceed further as in 3 the proof of Lemma 5.3 and arrive at 5.3). Remark 5.5. One can show by a standard method that u C[,T);L R 2 )). Here, it sufficestousethecontinuityofthebilinearformb as, e.g.,intheproofin[2, Theorem.]. The optimal hypercontractive estimate for p >. First, we improve the decay estimates from Lemma 5.4.
2D CHEMOTAXIS MODELS 5 Lemma 5.6. For each p 4 3,2), the solution of.).2) with u L R 2 ) satisfies 5.5) supt /p u,t) p <. t> Proof. By definition 5.) of the space Ẽ, one immediately sees that it is enough to prove 5.5) for t T. By the Duhamel formula 5.), 3.3) and 3.5), we have 5.6) ut) p e t u p +C e t u p +C t s) /2 +/p us) +γ) us) ds t s) /2 +/p us) p +γ) us) q ds Ct /p u +C T/2 t t s) 3/2+/p us) 2 p ds+c t s) 3/2+/p us) 2 p ds. Using 5.) we estimate the second term on the right-hand side of 5.6) as T/2 T/2 t s) 3/2+/p us) 2 p ds C t s) 3/2+/p s 2/p ) ds CT) t T 3/2+/p. 2) Hence for t T and in view of the fact that T/2 5.7) for t T it holds t 2 t T 2 relation 5.6) reads t /p ut) p C u +Ct /2 +Ct /p T/2 t s) 3/2+/p us) 2 p ds. In turn, in view of 5.) and owing to definition of ε in 5.4), for t T we arrive at 5.8) t /p ut) p C u +C+Ct /p Next we use the inequality T/2 T/2 t s) 3/2+/p s T 2) 2 /p) ε) ds C t s) 3/2+/p s T 2) 2 /p) ε) ds. t T 2) /p /2 2 /p) ε), to see that by 5.7) for t T 5.8) yields t /p ut) p C +Ct /2 t T 2) 2 /p) ε) C +Ct /2 2 /p) ε). ) Sincep 4,2),itisenoughtopickupε > smallenoughtoensurethat 2 3 2 p ε) <, and we arrive at 5.5) for t T, Lemma 5.6 is proved.
6 P. BILER, T. CIEŚLAK, G. KARCH, AND J. ZIENKIEWICZ The optimal decay estimate for other p, ). In view of Lemma 5.6, we see that Theorem 2.6 is true for p = and p 4,2). Since 3 for p, 4 ], we have by interpolation 3 ) ϑ, t /p ut) p ut) ϑ t / p ut) p where p < 2, ϑ = p p, and therefore Theorem 2.6 holds also for p, 4]. p 3 Now, we can interpolate estimate 5.5) for p 4,2) and 5.3) to get 5.5) with any 3 p [,2). The last step of the proof of Theorem 2.6 is the extrapolation of the hypercontractive estimates 5.5) for q [2, ). Actually, it is enough to obtain the decay estimate for q 2, ), the remaining case q = 2 will follow by simple interpolation. Taking q 2, ), we know that ut) q Ct /q u +C t s) /2 /r+/q us) +γ) us) r ds 5.9) Ct /q u +C t s) /2 /r+/q us) +γ) σ us) ρ ds Ct /q u +C t s) /2 /r+/q us) 2 σ ds. Here r is chosen in such a way that = 2, so that for σ r σ 2 4,2) we have r,2), 3 r close to 2. At the same time + = and =, the above choice of parameters ρ σ r ρ σ 2 allows us to apply 3.5) to 5.9). Relation 5.5) with σ 4,2) us) 3 σ Cs /σ, applied to 5.9) yields ut) q Ct /q u +C t s) /2 /r+/q s 2/σ ) ds. Since t s) /2 /r+/q s 2/σ ) ds = Ct /2 /r+/q+2/σ, we notice that ut) q Ct /q u +Ct /q holds in view of the relation 2 r + q + 2 σ = q. Thus, the decay estimate for q > 2 is proved.
2D CHEMOTAXIS MODELS 7 Acknowledgments This work was initiated during the visits of T. Cieślak at Uniwersytet Wroc lawski, T.C. would like to express his gratitude for support and hospitality. The first, the third and the fourth authors were supported by the NCN grant 23/9/B/ST/442. The second author was partially supported by the Polish Ministry of Science and Higher Education under the Iuventus Plus grant No. 73/IP3/2/7. The fourth author was also supported by the grant DEC-22/5/B/ST/692. References [] MR33387) P. Biler, The Cauchy problem and self-similar solutions for a nonlinear parabolic equation, Studia Math. 4, 8 25 995). [2] MR34) [.3934/cpaa.25.4.27] P. Biler, I. Guerra, G. Karch, Large global-in-time solutions of the parabolic-parabolic Keller-Segel system on the plane, Comm. Pure Appl. Analysis 4 25), 27 226. [3] MR2643796) [.7/s28-9-48-] P. Biler, G. Karch, Blowup of solutions to generalized Keller Segel model, J. Evol. Equ. 2), 247 262. [4] MR259677) [.88/95-775/22/7/3] P. Biler, G. Karch, P. Laurençot, Blowup of solutions to a diffusive aggregation model, Nonlinearity 22 29), 559 568. [5] [.88/95-775/28/2/4369] P. Biler, G. Karch, J. Zienkiewicz, Optimal criteria for blowup of radial and N-symmetric solutions of chemotaxis systems, Nonlinearity 28 25), 4369 4387. [6] MR66243) P. Biler, W. A. Woyczyński, Global and exploding solutions of nonlocal quadratic evolution problems, SIAM J. Appl. Math. 59 999), 845 869. [7] MR3444)[.464/ba63--6] P. Biler, J. Zienkiewicz, Existence of solutions for the Keller-Segel model of chemotaxis with measures as initial data, Bull. Polish Acad. Sci. Mathematics 63, 25) 4-5. [8] MR222697) A. Blanchet, J. Dolbeault, B. Perthame, Two-dimensional Keller-Segel model: optimal critical mass and qualitative properties of the solutions, Electron. J. Differential Equations 44, 32 pp. 26). [9] MR7289) Y. Giga, T. Miyakawa, H. Osada, Two-dimensional Navier-Stokes flow with measures as initial vorticity, Arch. Rational Mech. Anal. 4 988), 223 25. [] MR28292) [.464/am38-3-] G. Karch, K. Suzuki, Blow-up versus global existence of solutions to aggregation equations, Appl. Math. Warsaw) 38 2), 243 258. [] MR24726) [.7/s28-8-375-6] H. Kozono, Y. Sugiyama, Local existence and finite time blow-up of solutions in the 2-D Keller-Segel system, J. Evol. Equ. 8 28), 353 378. [2] MR972874) M. Kurokiba, T. Ogawa, Finite time blow-up of the solution for a nonlinear parabolic equation of drift-diffusion type, Differ. Integral Equ. 6 23), 427 452. [3] MR32922) P.-G. Lemarié-Rieusset, Small data in an optimal Banach space for the parabolicparabolic and parabolic-elliptic Keller-Segel equations in the whole space, Adv. Diff. Eq. 8 23), 89 28. [4] MR248755) [.7/s22-8-669-] D. Li, J. L. Rodrigo, Finite-time singularities of an aggregation equation in R n with fractional dissipation, Comm. Math. Phys. 287 29), 687 73. [5] MR249379) [.6/j.aim.28..6] D. Li, J. L. Rodrigo, Refined blowup criteria and nonsymmetric blowup of an aggregation equation, Adv. Math. 22 29), 77 738. [6] MR266636) D. Li, J. L. Rodrigo, X. Zhang, Exploding solutions for a nonlocal quadratic evolution problem, Rev. Mat. Iberoam. 26 2), 295 332. [7] MR783582) T. Nagai, Behavior of solutions to a parabolic-elliptic system modelling chemotaxis, J. Korean Math. Soc. 37 2), 72 732. [8] MR887324) [.55/S25583442] T. Nagai, Blowup of nonradial solutions to parabolicelliptic systems modeling chemotaxis in two-dimensional domains, J. Ineq. Appl. 6 2), 37 55. [9] MR2995) E. M. Stein, Singular integrals and differentiability properties of functions, Princeton Mathematical Series, No. 3, Princeton University Press, Princeton, N.J., 97.
8 P. BILER, T. CIEŚLAK, G. KARCH, AND J. ZIENKIEWICZ [2] MR337626)[DOI:.6/j.jde.24.2.33] Y. Sugiyama, M. Yamamoto, K. Kato, Local and global solvability and blow up for the drift-diffusion equation with the fractional dissipation in the critical space, J. Diff. Eq. 258 25), 2983 3. Instytut Matematyczny, Uniwersytet Wroc lawski, pl. Grunwaldzki 2/4, 5-384 Wroc law, Poland E-mail address: Piotr.Biler@math.uni.wroc.pl Instytut Matematyczny, Polska Akademia Nauk, ul. Śniadeckich 8, -956 Warsaw, Poland E-mail address: T.Cieslak@impan.pl Instytut Matematyczny, Uniwersytet Wroc lawski, pl. Grunwaldzki 2/4, 5-384 Wroc law, Poland E-mail address: Grzegorz.Karch@math.uni.wroc.pl Instytut Matematyczny, Uniwersytet Wroc lawski, pl. Grunwaldzki 2/4, 5-384 Wroc law, Poland E-mail address: Jacek.Zienkiewicz@math.uni.wroc.pl