arxiv: v1 [math.ap] 26 Mar 2019 Peter Y. H. Pang

Transcription

1 Asymptotic behavior of solutions to a tumor angiogenesis model with chemotaxis haptotaxis arxiv: v1 [math.ap] 26 Mar 219 Peter Y. H. Pang Department of Mathematics, National University of Singapore 1 Lower Kent Ridge Road, Republic of Singapore Yifu Wang School of Mathematics and Statistics, Beijing Institute of Technology Beijing 181, People s Republic of China Abstract This paper studies the following system of differential equations modeling tumor angiogenesis in a bounded smooth domain R N (N = 1,2): α p t = p p( c+ρ w)+λp(1 p), x,t >, 1+c c t = c c µpc, x,t >, w t = γp(1 w), x,t >, where α,ρ,λ,µ and γ are positive parameters. For any reasonably regular initial data (p,c,w ), we prove the global boundedness (L -norm) of p via an iterative method. Furthermore, we investigate the long-time behavior of solutions to the above system under an additional mild condition, and improve previously known results. In particular, in the one-dimensional case, we show that the solution (p, c, w) converges to (1,,1) with an explicit exponential rate as time tends to infinity. Key words: Angiogenesis, chemotaxis, haptotaxis, boundedness, asymptotic behavior 21 Mathematics Subject Classification: 35A1, 35B35, 35K57, 35Q92, 92C17 Corresponding author. wangyifu@bit.edu.cn 1

2 1 Introduction As a physiological process, angiogenesis involves the formation of new capillary networks sprouting from a pre-existing vascular network and plays an important role in embryo development, wound healing and tumor growth. For example, it has been recognized that capillary growth through angiogenesis leads to vascularization of a tumor, providing it with its own dedicated blood supply and consequently allowing for rapid growth and metastasis. The process of tumor angiogenesis can be divided into three main stages (which may be overlapping): (i) changes within existing blood vessels; (ii) formation of new vessels; and (iii) maturation of new vessels. Over the past decade, a lot of work has been done on the mathematical modeling of tumor growth; see, for example, [2, 3, 5, 6, 17, 28, 29] and the references cited therein. In particular, the role of angiogenesis in tumor growth has also attracted a great deal of attention; see, for example, [1, 7, 14, 24, 27] and the references cited therein. For example, in Levine et al. [14], a system of PDEs using reinforced random walks was deployed to model the first stage of angiogenesis, in which chemotactic substances from the tumor combine with the receptors on the endothelial cell wall to release proteolytic enzymes that can degrade the basal membrane of the blood vessels eventually. In this paper we consider a variation of the model proposed in [2], namely, α p t = p p( c+ρ w)+λp(1 p), x,t >, 1+c c t = c c µpc, x,t >, w t = γp(1 w), x,t >, p ν p( α c 1+c ν +ρ w ν ) = c =, x,t >, ν p(x,) = p (x), c(x,) = c (x), w(x,) = w (x), x, (1.1) in a bounded smooth domain R N (N = 1,2), where, in addition to random motion, the existing blood vessels endothelial cells p migrate in response to the concentration gradient of a chemical signal c (called Tumor Angiogenic Factor, or TAF) secreted by tumor cells as well as the concentration gradient of non-diffusible glycoprotein fibronectin w produced by the endothelial cells [21]. The former directed migration is a chemotatic process, whereas the latter is a haptotatic process. In this model, it is assumed that the endothelial cells 2

3 proliferate according to a logistic law, that the spatio-temporal evolution of TAF occurs through diffusion, natural decay and degradation upon binding to the endothelial cells, and that the fibronectin is produced by the endothelial cells and degrades upon binding to the endothelial cells. For the remainder of this paper, we shall assume that the initial data satisfy the following: (p,c,w ) (C 2+β ()) 3 is nonnegative for some β (,1) with p, p ν p α c ( 1+c ν +ρ w ν ) = c (1.2) ν =. The present paper focuses on the global existence and asymptotic behavior of classical solutions to (1.1). Let us look at two subsystems contained in (1.1). The first is a Keller Segel-type chemotaxis system with signal absorption: p t = p (p c)+λp(1 p), x,t >, c t = c pc, x,t >. (1.3) It is known that, unlike the standard Keller Segel model, (1.3) with λ = possesses global, bounded classical solutions in two-dimensional bounded convex domains for arbitrarily large initial data; while in three spatial dimensions, it admits at least global weak solutions which eventually become smooth and bounded after some waiting time [32]. In the high-dimensional setting, it has been proved that global bounded classical solutions exist for suitably large λ >, while only certain weak solutions are known to exist for arbitrary λ > [13]. Another delicate subsystem of (1.1) is the haptotaxis-only system obtained by letting α = in (1.1): p t = p ρ (p w)+λp(1 p), x,t >, w t = γp(1 w), x,t >. (1.4) Here, since the quantity w satisfies an ODE without any diffusion, the smoothing effect on the spatial regularity of w during evolution cannot be expected. To the best of our knowledge, unlike the study of chemotaxis systems, the mathematical literature on haptotaxis systems is comparatively thin. Indeed, the literature provides only some results on global solvability 3

4 in various special models, and the detailed description of qualitative properties such as long-time behaviors of solutions is available only in very particular cases (see, for example, [8, 18, 19, 2, 31, 35, 36, 39]). More recently, some results on global existence and asymptotic behavior for certain chemotaxis haptotaxis models of cancer invasion have been obtained (see, for example, [16, 22, 23, 3, 33, 34, 37]). Particularly, Hillen et al. [9] have shown the convergence of a cancer invasion model in one-dimensional domains and the result has been subsequently extended to higher dimensions [16, 34, 37]. In [21], in two spatial dimensions, the authors showed the global existence and long-time behavior of classical solutions to (1.1) when the initial data (p,c,w ) satisfy either w > 1 or w 1 L () < δ for some δ > (see Lemma 5.8 of [21]). Generalizing this result, our first main result establishes that, for any choice of reasonably regular initial data (p,c,w ), the L -norm of p is globally bounded. This is done via an iterative method. Theorem 1.1. Letα,ρ,λ,µandγ be positiveparameters. Thenfor anyinitial data (p,c,w ) satisfying (1.2), the problem (1.1) possesses a unique classical solution (p, c, w) comprising nonnegative functions in C( [, )) C 2,1 ( (, ) such that p(,t) L () C for all t >. Next, we investigate the asymptotic behavior of solutions to (1.1). Under an additional mild condition on the initial data w, we will show that the solution (p,c,w) converges to the spatially homogeneous equilibrium (1,, 1) as time tends to infinity. Theorem 1.2. Let α,ρ,λ,µ and γ be positive parameters, and suppose that (1.2) is satisfied and w > 1 1 ρ. Then the solution (p,c,w) C( [, )) C 2,1 ( (, ) of (1.1) satisfies lim p(,t) 1 L t r () + c(,t) W 1,2 () + w(,t) 1 L r () = (1.5) for any r 2. In particular, if N = 1, then for any ǫ (,min{λ 1,1,γ,λ}) there exists C(ǫ) > such that p(,t) 1 L () C(ǫ)e (min{λ 1,1,γ,λ} ǫ)t, (1.6) c(,t) W 1,2 () C(ǫ)e (1 ǫ)t, (1.7) 4

5 w(,t) 1 W 1,2 () C(ǫ)e (γ ǫ)t, (1.8) where λ 1 > is the first nonzero eigenvalue of in with the homogeneous Neumann boundary condition. The main mathematical challenge of the full chemotaxis haptotaxis system is the strong coupling between the migratory cells p and the haptotactic agent w. This strong coupling hasanimportant effect onthespatial regularity ofpandw. Infact, thelack ofregularization effect in the spatial variable in the w-equation and the presence of p therein demand tedious estimates on the solution. The key ideas behind our results are as follows: As pointed out in [34], the variable transformation z := pe ρw plays an important role in the examination of global solvability for the full chemotaxis haptotaxis model in the twoand higher-dimensional setting. However, due to the presence of the additional chemotaxis term in our model, this approach is not directly applicable to our problem. Instead, in the derivation of Theorem 1.1, we introduce the variable transformation q := p(c + 1) α e ρw as in [21], and thereby ensure that q(,t) is bounded in L n () for any finite n (see Lemma 2.4). It is essential to our approach to derive a bound for + t+τ 2 from q2m+1 t q2m the bound of t+τ (m = 1,2,...) by making appropriate use of (2.3) (2.4) in Lemma t q2m 2.3 (see (2.7) below). with The crucial idea of the proof of Theorem 1.2 is to show d dt F(p(t),w(t))+ 1 p p w 2 C p w 1 +C p c 2 2 p 2 F(p,w) = κ w 2 + p(lnp 1)+ p(w 1) γκ p(w 1) 2 for some C > (see (3.1)), on the basis of the global boundedness of the L -norm of p provided by Theorem 1.1 and under the mild assumption that w > 1 1. Furthermore, in ρ the one-dimensional case, with the help of higher regularity estimates of z with z = pe ρw (see (3.29)), we show that p(,t) converges to 1 uniformly in as t. From this, we derive the exponential convergence of solutions as desired. 5

6 2 Proof of Theorem 1.1 In this section, we first recall the following estimates for the heat semigroup (e τ ) τ in R N under the Neumann boundary condition. We shall omit the proof thereof and refer the interested readers to [4, 1, 38]. Lemma 2.1. Let (e τ ) τ be the Neumann heat semigroup in and λ 1 > the first nonzero eigenvalue of with homogeneous Neumann boundary condition. Then there exist positive constants k 1,k 2,k 3 such that: i) If 1 q p, then for all ϕ L q () with ϕ =, e τ ϕ L p () k 1 (1+τ N 2 (1 q 1 p ) )e λ 1τ ϕ L q (); ii) If 1 q p, then for all ϕ L q (), e τ ϕ L p () k 2 (1+τ 1 2 N 2 (1 q 1 p ) )e λ 1τ ϕ L q (); iii) If 1 q p, then for all ϕ C 1 ( ;R N ) with ϕ ν = on, e τ ϕ L p () k 3 (1+τ 1 2 N 2 (1 q 1 p ) )e λ 1τ ϕ L q (). Next, we recall the following result on local existence and uniqueness of classical solutions to (1.1) as well as a convenient extensibility criterion, which follows from Theorem 3.1, Lemma 5.9 and Theorem 5.1 of [21]. Lemma 2.2. ([21]) Let R N be a smooth bounded domain. There exists T max (, ] such that the problem (1.1) possesses a unique classical solution satisfying (p,c,w) (C( [,T max )) C 2,1 ( (,T max )) 3. Moreover, for any s > N +2, if T max < +. lim sup tրt max p(,t) W 1,s () (2.1) From now on, let (p,c,w) be the local classical solution of (1.1) on (,T max ) provided by Lemma 2.1, and τ := min{1, Tmax 6 }. The following basic but important properties of the solution to (1.1) can be directly obtained via standard arguments. 6

7 Lemma 2.3. ([21]) There exists a positive constant C independent of time such that t p(,t) c 1 := max{ p, }, e 2s p 2 ds C for all t (,T max ), (2.2) c(t) c L ()e t, t+τ t p 2 c 1 (1+ 1 λ ) for all t (,T max τ), (2.3) c(t) 2 + t ( c 2 + c 2 ) C for all t (,T max ), (2.4) w(t) max{ w L (),1} for all t (,T max ). (2.5) As the proof of Theorem 1.1 in the one-dimensional case is similar to that for two dimensions, henceforth in this section, we shall focus on the case N = 2. First, we shall show that p remains bounded in L n () for any finite n. We note that the L n ()-bound in Lemma 3.1 of [21] depends on the time variable. Lemma 2.4. For any n (1, ), there exists a positive constant C(n,τ) independent of t, such that p(,t) L n () C(n,τ) for all t (,T max ). Proof. Let q := p(c+1) α e ρw. As in the proof of Lemma 3.1 in [21], we infer that for any m = 1,2,... there exist constants c(m) > depending upon m and c 1 > such that d q 2m (c+1) α e ρw + q 2m 1 2 c(m)( c 2 +1) q 2m +c(m)( q 2m ) 2 +c 1. (2.6) dt Next, we use induction to show q 2m + t+τ Taking m = 1 in (2.6), we get d q 2 (c+1) α e ρw + q 2 c(1)( dt t q 2m 1 2 C(m). (2.7) c q 2 ) q 2 (c+1) α e ρw +c 1, (2.8) which implies that for the functions y(t) = q2 (c+1) α e ρw and a(t) = c(1)( c q2 ), we have dy dt a(t)y +c 1. On the other hand, for any given t > τ, it follows from (2.3) that there exists some t [t τ,t] such that y(t ) c 1 τ (1+ 1 ). Hence by ODE comparison argument we get λ t t t y(t) y(t )e a(s)ds +c 1 e t s a(τ)dτ ds c 2. (2.9) t 7

8 In this inequality, we have taken t = if t τ and noticed that t t τ a(s)ds c 3 for all t < T max by Lemma 2.2. Combining (2.8) with (2.9), one can see that (2.7) is indeed valid for m = 1. Now, suppose that (2.7) is valid for an integer m+1 = k 2, i.e., q 2k 1 + t+τ By the Gagliardo Nirenberg inequality in two dimensions t q 2k 2 2 C(k). (2.1) z 4 L 4 () c 3 z 2 L 2 () z 2 L 2 () +c 4 z 4 L 2 (), and hence q 2k c 3 q 2k 2 2 q 2k 1 +c 4 ( q 2k 1 ) 2. (2.11) Integrating (2.11) between t and t+τ and taking (2.1) into account, we have t+τ t q 2k c 5 (k), (2.12) which implies that for any t τ, there exists some t [t τ,t] such that q2k (t ) c 6. At this point, let y(t) := q2k (c+1) α e ρw and b(t) = c(k)( c q2k ). Then (2.6) can be rewritten as By the argument above, one can obtain dy dt + q 2k 1 2 b(t)y +c 1. q 2k + t+τ t q 2k 1 2 C, (2.13) and thereby conclude that (2.7) is valid for all integers m 1. The proof of Lemma 2.3 is now complete in view of the boundedness of the weight (c+1) α e ρw. To establish a priori estimates of p(,t) L (), we need some fundamental estimates for the solution of the following problem: c t = c c+f, x,t >, c =, x,t >, ν c(x,) = c (x), x. (2.14) 8

9 Lemma 2.5. ([15, Lemma 2.2]) Let T >, r (1, ). Then for each c W 2,r () with c ν = on and f Lr (,T;L r ()), (2.14) has a unique solution c W 1,r (,T;L r ()) L r (,T;W 2,r ()) given by c(t) = e t e t c + t e (t s) e (t s) f(s)ds, t [,T], where e t is the semigroup generated by the Neumann Laplacian, and there is C r > such that t t e rs c(x,s) r dxds C r e rs f(x,s) r dxds+c r v W 2,r (). (2.15) Now applying these estimates to control the cross-diffusive flux appropriately, we can derive the boundedness of p in (,T max ). Lemma 2.6. There exists a constant C > independent of t such that p(,t) L () C for all t (,T max ). Proof. We will only give a sketch of the proof, which is similar to that of Lemma 3.13 of [21]. For k max{2, p L ()}, let q k = max{q k,} and k (t) = {x : q(x,t) > k}. Multiplying the equation of q by q k, we obtain d qk 2 dt (c+1)α e ρw +2 q k 2 +2 qk 2 +8 qk 2 (c+1)α e ρw c 1 qk 3 +c 1k qk 2 +c 1k 2 q k +c 1 (qk 2 +kq k) c (2.16) for some c 1 > independent of k. By the boundedness of q in L n () for any n > 1, the Gagliardo Nirenberg inequality and Young inequality, we obtain c 1 q k 3 L 3 () 1 4 q k 2 H 1 () +c 2 q k L 1 (), c 1 k q k 2 L 2 () 1 4 q k 2 H 1 () +c 2k 2 q k L 1 (), c 1 qk c q k 2 H 1 () +c 2 q k 2 L 2 () c 2 L 2 (), c 1 k q k c 1 4 q k 2 H 1 () +c 2k 2 (1+ c 8 L 8 () ) k

10 Inserting the above estimates into (2.16), we have d qk 2 dt (c+1)α e ρw + q k 2 + qk 2 +8 q 2 k (c+1)α e ρw c 2 q k 2 L 2 () c 2 L 2 () +c 2k 2 (1+ c 8 L 8 () ) k 3 2 +(c 1 +2c 2 )k 2 q k L 1 () c 2 q k 2 L 2 () c 2 L 2 () q k 2 H 1 () +c 3k 4 (1+ c 8 L 8 () ) k 3 2. On the other hand, according to the relation between distribution functions and L p integrals (see e.g. (2.6) of [26]), we can see that (r +1) Hence taking into account Lemma 2.3, we get and thus (k 1) 16 k (t) < d dt k k 1 c 2 c 2 L 2 () s r s (t) ds = q(t) r+1 L r+1 (). s 16 s (t) ds < qk 2 (c+1)α e ρw +8 Therefore if h(t) = 8 c 2 c 2 L 2 (), then q 2 k(c+1) α e ρw c 4 e t s 16 s (t) ds 1 17 q(,t) 17 L 17 () q 2 k (c+1)α e ρw (2.17) qk(c+1) 2 α e ρw +c 4 (1+ c 8 L 8 () ) k 5 4. h(s)ds t (1+ c 8 L 8 () )e s h(σ)dσ k (s) 5 4 ds. Furthermore, since e t h(s)ds = e 8t e c t 2 c 2 L 2 ds () c 5 e 8t by Lemma 2.2 and e s h(σ)dσ e 8s, we get t qk 2 c 6 e 8(t s) (1+ c 8 L 8 () ) k(s) 5 4 ds c 6 t e 8(t s) (1+ c 8 L 8 () )ds sup k (t) 5 4. t To estimate the integral term in the right-side of the above inquality, we apply Lemma 2.5 with r = 8 and Lemma 2.4 to get t e 8(t s) c 8 L 8 () ds c 7 and thus q2 k c 8(sup k (t) ) 5 4. t On the other hand, q2 k (t) j (t) q2 k (t) (j k)2 j (t) for j > k. Consequently (j k) 2 sup t j (t) c 8 (sup k (t) ) 5 4 k (t). t 1

11 According to Lemma B.1 of [11], there exists k < such that k (t) = for all t (,T max ). Therefore q(,t) L () k for any t (,T max ) and thereby the proof is complete. Proof of Theorem 1.1. By the boundedness of p in L ((,T max ),L ()) from Lemma 2.6 and a bootstrap argument as in [21], we can see that the global existence of classical solutions to (1.1) is an immediate consequence of Lemma 2.2, i.e., T max =. Indeed, supposed that T max <, then by Lemma 3.15 and Lemma 3.19 of [21], we can see that for any s > N +2 and t T max c(,t) W 1,s () + w(,t) W 1,s () C. Further by Lemma 3.2 of [21], we have p(,t) W 1,s () C which contradicts (2.1) and thus implies that T max =. Moreover, since τ := min{1, Tmax 6 } = 1, there exists a constant C > independent of time t such that p(,t) L () C for all t by retracing the proofs of Lemma 2.4 and Lemma 2.6. This completes the proof of Theorem Proof of Theorem 1.2 In this section, on the basis of the L -bound of p provided by Theorem 1.1, we shall look at the asymptotic behavior of the solution (p,c,w) of the problem (1.1). 3.1 L r -convergence of solutions in two dimensions When either w > 1 or w 1 L () < δ for some δ >, the authors of [21] removed the time dependence of the L -bound of p (see Lemma 5.8 of [21]) and thereby investigated the asymptotic behavior of solutions to (1.1). In this subsection, on the basis of the L -boundof p being independent of time as provided by Theorem 1.1, we shall derive the same estimates as in Lemma 5.6 and Lemma 5.7 of [21] under the weaker assumption that w > 1 1. ρ We shall show that the solution (p,c,w) to (1.1) converges to the homogeneous steady state (1,,1) as t. 11

12 Before going into the details, let us first collect some useful related estimates. It should be noted that no other assumptions on the initial data (p,c,w ) are made except for reasonable regularity, i.e., (1.2). Lemma 3.1. ([21, Lemmas 3.4, 5.1, 5.2, 3.8]) Let (p,c,w) be the global, classical solution of (1.1). Then c(t) 2 e 2t ( p(1 p) max{ p, }/λ; (3.1) p w 1 w 1 L 1 (); (3.2) p c 2 < (3.3) c 2 +µ 2 c 2 L ()max{ Lemma 3.2. Under the assumptions of Theorem 1.1, we have Proof. We know that c solves the linear equation p, }(t+ 1 λ ) ). (3.4) sup c(t) W 1, () C. (3.5) t c t = c c+f under the Neumann boundary condition with f := µpc. Since p, we know that c(x,t) c L ()e t by the standard sub-super solutions method. On the other hand, by Theorem 1.1, sup p(t) L () c 1, which readily implies that f L ((, );L ()) c 1. t Now upon a standard regularity argument we can deduce the desired result. For the reader s convenience, we only give a brief sketch of the main ideas, and would like refer to the proof of Lemma 1 in [12] or Lemma 4.1 in [1] for more details. Indeed, according to the variationof-constants formula of c, we have for t > 2 c(,t) = e (t 1)( 1) c(,1)+ So by Lemma 2.1(ii), we infer that t 1 e (t s)( 1) f(,s)ds. t c(,t) L () 2k 2 c(,1) L 1 () +k 2 (1+(t s) 1 2 )e (t s)(λ 1 +1) f(,s) L ()ds 2k 2 c(,1) L 1 () +k 2 c 1 (1+σ 1 2 )e σ dσ. 1 12

13 Lemma 3.3. ([21, Lemma 5.4]) Let (p,c,w) be the global, classical solution of (1.1). Then for every t and κ >, where and d F(p(t),w(t)) = G(p(t),w(t),c(t)), (3.6) dt F(p,w) = κ w 2 + p(lnp 1)+ p(w 1) γκ p(w 1) 2, G(p,w,c) = + p 2 + p α 1+c p c+ p (2αγκ(1 w)+αρ) 1+c c w p(1 p)lnp p 2 (1 w) p(1 p)(w 1) 2. (ρ 2 2γκ+2ργκ(1 w))p w 2 +λ +λρ p(1 p)(w 1)+γρ +2γ 2 κ p 2 (w 1) 2 λγκ Lemma 3.4. If w > 1 1, then there exists κ > such that ρ G(p,w,c) 1 p 2 2 p 1 p w 2 +C 2 p w 1 +C p c 2 (3.7) for some C >. and Proof. By the Hölder and Young inequalities we have for some c 1 >. α 1+c p c 1 p 2 + α2 p c 2, 2 p 2 p (2αγκ(1 w)+αρ) 1+c c w 1 2 p w 2 +c 1 p c 2 As w > 1 1 ρ, we can find some ε 1 > such that ρ(1 w ) + 1 ε 1, where (1 w ) + = max{,1 w }. Hence from the w-equation in (1.1), it follows that 1 w = (1 w )e γ t p(s)ds, (3.8) 13

14 and thus (ρ 2 2γκ+2ργκ(1 w))p w 2 (ρ 2 2γκ+2ργκ(1 w ) + )p w 2 (ρ 2 2γκε 1 )p w 2 p w 2 if we pick κ > sufficiently large such that ρ 2 2γκε 1 < 1. Denote the lower-order terms of G(p,w,c) by θ(p,w), i.e., θ(p,w) =: λ p(1 p)lnp+λρ p(1 p)(w 1)+γρ +2γ 2 κ p 2 (w 1) 2 λγκ p(1 p)(w 1) 2. Since s(1 s)lns for s, we get θ(p,w) λρ p(1 p)(w 1)+γρ p 2 (1 w) +2γ 2 κ p 2 (w 1) 2 λγκ p(1 p)(w 1) 2 c 2 ( p L (), w 1 L ()) p w 1, p 2 (1 w) which, alongwith p(,t) L () C fromtheorem1.1and w(,t) 1 L () w 1 L () from (3.8), yields θ(p,w) c 3 The desired result (3.7) then immediately follows. Lemma 3.5. If w > 1 1 ρ, then p w 1. sup w(t) 2 + t p 2 p + p w 2 <. (3.9) Proof. Combining Lemmas 3.3 and 3.4, we have d dt F(p(t),w(t))+ 1 p p w 2 C 2 p 2 p w 1 +C p c 2. (3.1) Hence (3.9) follows upon integration on the time variable, and using (3.2) and (3.3). Lemma 3.6. If w > 1 1, then for any r 2 ρ lim p(,t) p(t) L t r () =, (3.11) 14

15 where p(t) = 1 p(,t), and lim p(t) 1 =, (3.12) t lim w(,t) 1 L t r () =. (3.13) Proof. The proofs of (3.11) and (3.12) are similar to those of Lemma of [21] respectively. However, for the reader s convenience, we only give a brief sketch of (3.12). In fact, from (1.1) and the Poincaré Wirtinger inequality, it follows that p t = λ(p p 2 1 (p p) 2 p 2 λp(1 p c 1 ). p Hence by (3.9), we get p(t) p exp{λt λ c 2 exp{λt λ t t p(s)ds c 1 λ p(s)ds}, p 2 ds} p which means that (3.12) is valid due to either p(t) 1 or p(t) in Lemma 5.1 of [21]. Indeed, supposed that p(t), then there exists t > 1 such that p(t) 1 and thus 2 t p(s)ds t + t p(s)ds for all t t 2. Therefore we arrive at p(t) c 3 e λt 2 for all t t, which contradicts p(t). Now we turn to show (3.13). Invoking the Poincaré inequality in the form ϕ(x) 1 ϕ(y)dy 2 dx C p ϕ 2 dx for all ϕ W 1,2 () for some C p >, one can find that for all j N j+1 j j+1 p(s) p(s) 2 L 2 () ds C p p(s) 2 L 2 () ds j j+1 C p sup p(t) L () t j p(s) 2 ds, p(s) (3.14) which, along with (3.9) and Theorem 1.1, shows that as j. j+1 j p(x,s) p(s) 2 dsdx = 15 j+1 j p(s) p(s) 2 L 2 () ds (3.15)

16 Now defining p j (x) := j+1 j p(x,s) p(s) 2 ds, x,j N, (3.15) tells us that p j in L 1 () as j. There exist a certain null set Q and a subsequence (j k ) k N N such that j k and p jk (x) for every x \Q as k. Restated in the original variable, this becomes for every x \Q as k. jk +1 j k p(x,s) p(s) 2 ds (3.16) Therefore, from (3.8) and p(x,t), it follows that for any x \Q w(x,t) 1 w 1 L ()exp{ γ [t] p(x,s)ds} w 1 L ()exp{ γ m(t) k= w 1 L ()exp{γ m(t) k= jk +1 j k jk +1 j k p(x,s)ds} p(x,s) p(s) ds γ m(t) k= w 1 L ()exp{γ m(t) k= ( j k +1 j k p(x,s) p(s) 2 ds) 1 2 jk +1 j k p(s)ds} γ m(t) k= jk +1 j k p(s)ds}, (3.17) where m(t) := max k N {j k + 1,[t]}. Furthermore, by (3.12), there exists k N such that jk +1 j k p(s)ds 1 2 for all k k. Hence by the fact that m(t) as t and (3.16), we obtain that w(x,t) 1 almost everywhere in as t. On the other hand, as w(x,t) 1 w 1 L (), the dominated convergence theorem ensures that (3.13) holds for any r (2, ). Remark ) It is observed that since W 1,2 () L () is invalid in the twodimensional setting, p(s) p(s) 2 L 2 () in (3.14) cannot be replaced by p(s) p(s) 2 L (), and thus we cannot infer that lim t w(,t) 1 L () =, even though we have established that all the related estimates of (p,c,w) in [21] continue to hold under the milder condition imposed on the initial data w. 2) Similarly to the remark above, we note that, even though w(,t) W 1,n () C(T) for any n 2 and t T, we are not able to infer the global estimate sup t w(t) 2+ε C. Otherwise, we would be able to apply regularity estimates for bounded solutions of semilinear parabolic equations (see [25] for instance) to obtain the Hölder estimates of 16

17 p(x,t) in (1, ), and thereby conclude lim t p(,t) 1 L () =. As things stand at the moment, we are only able to infer convergence in L r. 3.2 L -convergence of solutions with exponential rate in one dimension It is observed that the results, in particular Lemma 3.6, in the previous subsection are still valid in the one-dimensional case. Moreover, in the one-dimensional setting, the weak convergence result in Lemma 3.6 can be improved via a bootstrap argument. In fact, we shall derive some a priori estimates of (p,c,w) and thereby demonstrate that (p,c,w) converges to (1,,1) in L () as t. Furthermore, by a regularity argument involving the variationof-constants formula for p and smoothing L p L q type estimates for the Neumann heat semigroup, we will show that p(,t) 1 decays exponentially in L (). As pointed out in the Introduction, the main technical difficulty in the derivation of Theorem 1.2 stems from the coupling between p and w. Indeed, the lack of regularization effect in the space variable in the w-equation and the presence of p there demand tedious estimates of the solution. The following lemma plays a crucial role in establishing the uniform convergence of p as t (see Lemma 3.11). Thought the proof thereof only involves elementary analysis, we give a full proof here for the sake of the reader s convenience since we could not find a precise reference covering our situation. Lemma 3.7. Let k(t) be a function satisfying k(t), k(t)dt <. If k (t) h(t) for some h(t) L 1 (, ), then k(t) as t. Proof. Supposing the contrary, then we can find A > and a sequence (t j ) j N (1, ) such that t j t j 1 +2, t j as j and k(t j ) A for all j N. On the other hand, by k (t) h(t), we have k(t j τ) k(t j ) tj h(s) ds k(t j ) tj t j τ t j 1 17 h(s) ds (3.18)

18 for all τ (,1). Since h(t) L 1 (, ), we have t j t j 1 h(s) ds as t j and thereby there exists j N such that t j t j 1 h(s) ds A 2 for all j j, which along with (3.18) implies that k(t j τ) k(t j ) A 2 A 2 (3.19) for j j and τ (,1). It follows that t j t j 1 k(t)dt A 2 for all j j, which contradicts k(t)dt < and thus completes the proof of the lemma. Lemma 3.8. If w > 1 1, then there exists a constant C > such that ρ where z = pe ρw. e ρw z x 2 C (3.2) Proof. We know that z x = e ρw p x ρzw x and thus e ρw z x 2 4e ρw p x 2 +4p 2 e ρw ρ 2 w x 2. Integrating over (, ) and taking (2.5) into account, we have e ρw z x 2 4 p x 2 +4ρ 2 sup p(t) L () t 4(1+ρ 2 p x 2 )sup p(t) L ()( + t p Hence by Theorem 1.1 and Lemma 3.5, we get (3.2). p w x 2 p w x 2 ). Lemma 3.9. If w > 1 1, then there exists a constant C > such that ρ d dt e ρw z x e ρw zt 2 3 C( e ρw z x 2 + p w x 2 + c xx 2 + p c x 2 + p w 1 + p(p 1) 2 ) (3.21) with z = pe ρw. Proof. Note that z satisfies z t = e ρw (e ρw z x ) x e ρw ( z xe ρw 1+c c) x +λz(1 ze ρw ) ργe ρw z 2 (1 w). 18

19 Multiplying the above equation by z t e ρw and integrating in the spatial variable, we obtain e ρw zt 2 + e ρw z x z xt = e ρw α z t ( 1+c z xc x + αzρ 1+c w xc x αz (1+c) 2 c x 2 + αz 1+c c xx) (3.22) + e ρw z t (λz(1 ze ρw ) ργe ρw z 2 (1 w)). Notice that e ρw z x z xt = 1 d e ρw z x 2 γρ e 2ρw z(1 w) z x 2 2dt 2 1 d e ρw z x 2 γρ 2dt 2 sup p(t) L () 1 w L () e ρw z x 2, t αe ρw z t 1+c z xc x 1 e ρw zt 2 +c 1 sup c x (t) 2 L 6 () e ρw z x 2, t αze ρw ρ z t 1+c w xc x 1 e ρw zt 2 6 +c 1sup c x (t) 2 L () sup p(t) L () p w x 2, t t αze ρw z t (1+c) 2 c x 2 1 e ρw zt 2 +c 1 sup p(t) L 6 ()sup c x (t) 2 L () p c x 2, t t z t αze ρw 1+c c xx 1 6 e ρw z t (λz(1 ze ρw ) ργe ρw z 2 (1 w)) = λ pz t (1 p) ργ z t p 2 (1 w)) 1 e ρw zt 2 +c 1 sup p(t) L 6 () t e ρw zt 2 +c 1sup p(t) 2 L () t c xx 2, p(1 p) 2 +c 1 sup p(t) 3 L () 1 w L () t p 1 w. Applying Theorem 1.1, (3.5) and inserting the above inequalities into (3.22), we obtain (3.21). Now we focus our attention on the decay properties of the solutions. Indeed, we will show that p(x,t) converges to 1 with respect to the norm in L () as t. Subsequently, we will establish the exponential decay of p(,t) 1 L () with explicit rate. Lemma 3.1. If w > 1 1 ρ, then lim w(,t) 1 L t () =. (3.23) 19

20 Proof. From (3.8), it follows that for any ǫ > w(x,t) 1 w 1 L ()exp{ γ t p(s)ds} where p(t) = 1 p(,t). w 1 L ()exp{γ t p(s) p(s) L ()ds γ t p(s)ds} w 1 L ()exp{ γ ǫ t p(s) p(s) 2 L () ds+ǫγt γ t p(s)ds}, (3.24) On the other hand, by the Poincaré Wirtinger inequality, the Sobolev imbedding theorem in one dimension and (3.9), we have t p(s) p(s) 2 L () ds c 1 p x (s) 2 L 2 () ds c 1 sup p(t) L () t p x (s) 2 ds p(s) (3.25) c 2 for some constant c 2 >. Combining (3.24) with (3.25) yields w(t) 1 L () w 1 L ()exp{ c 2γ ǫ +ǫγt γ t p(s)ds} (3.26) for t. The assertion now follows from the last inequality and the proof is complete. Lemma If w > 1 1 ρ, then Proof. We first show that lim p(,t) 1 L t () =. (3.27) lim z(,t) z(t) L t () = (3.28) where z(t) = 1 z(,t). To this end, we consider the function k(t) defined by k(t) = eρw z x 2 and prove that By Lemmas 3.7, 3.8 and 3.9, it is enough to prove that h(t) := lim k(t) =. (3.29) t e ρw z x 2 + p w x 2 + c xx 2 + p c x 2 + p w 1 + p(p 1) 2 L 1 (, ). 2

21 Noting (3.2), (3.9), (3.3), (3.2) and (2.4), it remains to estimate p(p 1)2. In fact, multiplying the p-equation in (1.1) by p 1, we have 1 d (p 1) 2 = p x 2 p +ρ pw x p x +α 2dt 1+c c xp x λ p(p 1) 2 1 (3.3) p x 2 +C( p 2 w x 2 + p 2 c x 2 ) λ p(p 1) 2. 2 Hence, by the boundedness of p, (3.3) and (3.9), we easily infer that p(p 1)2 C. Furthermore, by the Poincaré Wirtinger inequality and the Sobolev imbedding theorem in one dimension, we have z(t) z(t) L () C p z x (t) L 2 (), (3.31) which along with (3.29) yields (3.28). On the other hand, for any {t j } j N (1, ), there exists a subsequence along which z(,t j ) e ρ a.e. in as j by Lemma 3.6. We apply the dominated convergence theorem along with the uniform majorization z(,t j ) sup z(t j ) L () C to infer that j 1 Hence lim t z(t) e ρ =. (3.32) p(,t) 1 L () = e ρw z 1 L () e ρ(1+ w () ) ( z(,t) z(t) L () + z(t) e ρ )+ e ρw(,t) e ρ L () e ρ(1+ w () ) ( z(,t) z(t) L () + z(t) e ρ )+c 1 w(,t) 1 L () for some c 1 >, which, together with (3.23), (3.28) and (3.32), yields the desired result. Having established that p(x,t) converges to 1 uniformly with respect to x as t, we now go on to establish an explicit exponential convergence rate. Using (3.27), we first look into the decay of w x(t) 2. Lemma Let w > 1 1. Then for any ǫ >, there exists C(ǫ) > such that ρ w x (t) 2 C(ǫ)e 2γ(1 ǫ)t. (3.33) 21

22 Proof. From (3.8), it follows that w x (t) 2 4 w x 2 e 2γ t p(s)ds +4γ 2 w 1 2 e 2γ t p(s)ds ( t 4 w x 2 e 2γ t p(s)ds +4tγ 2 w 1 2 e 2γ t p(s)ds t ) 2 p x (s) ds p x (s) 2 ds. Taking Lemma 3.11 into account, we know that for any ǫ >, there exists t ǫ > 1 such that p(x,t) > 1 ǫ for all x,t > t ǫ. Therefore integrating the above inequality in the space variable yields w x (t) 2 4e 2γ(1 ǫ)(t tǫ) w x 2 +4tγ 2 w 1 2 L ()e 2γ(1 ǫ)(t tǫ) c 1 (ǫ)(1+t)e 2γ(1 ǫ)t ( w x 2 + w 1 2 L ()sup p(t) L () t for all t > t ǫ, which, along with (3.9), implies (3.33). p x 2 p x 2 p ), Now we utilize thedecay properties of c x(t) 2, w x(t) 2 and theuniform convergence of p(x,t) 1 asserted by Lemma 3.11 to establish the decay property of p(,t) 1 L 2 (). Lemma Let w > 1 1. Then for any ǫ (,min{1,γ,λ}), there exists C(ǫ) > ρ such that p(,t) 1 L 2 () C(ǫ)e (min{1,γ,λ} ǫ)t. (3.34) Proof. By (3.27), we know that for any ǫ (,min{1,γ,λ}), there exists t ǫ > 1 such that p(x,t) > 1 ǫ for all x,t > t ǫ. Hence, we multiply the p-equation in (1.1) by p 1 and integrate the result over to get 1 d (p 1) 2 = p x 2 +ρ pw x p x +α 2dt 1 p x 2 +c 1 ( w x p 1+c c xp x λ p(p 1) 2 (3.35) c x 2 ) λ(1 ǫ) (p 1) 2 for all t > t ǫ. Now, applying the Gronwall inequality, (3.4) and Lemma 3.12, we have t (p(t) 1) 2 e 2λ(1 ǫ)(t tǫ) (p(t ǫ ) 1) 2 +c 1 e 2λ(1 ǫ)(t s) ( w x 2 + c x 2 ) c 2 (ǫ)e 2λ(1 ǫ)t +c 3 (ǫ) c 4 (ǫ)e 2min{λ,1,γ}(1 ǫ)t, t e 2λ(1 ǫ)(t s) (e 2γ(1 ǫ)s +e 2(1 ǫ)s ) where c i (ǫ) > (i = 2,3,4) are independent of time t. This completes the proof. 22

23 Moving forward, on the basis of Lemma 3.13, we come to establish the exponential decay of p(,t) p(t) L () by means of a variation-of-constants representation of p, as follows: Lemma Let w > 1 1 ρ. Then for any ǫ (,min{λ 1,1,γ,λ}), there exists C(ǫ) > such that p(,t) p(t) L () C(ǫ)e (min{λ 1,1,γ,λ} ǫ)t. (3.36) Proof. By noting that p t = λp(1 p)(t), applying the variation-of-constants formula to the p-equation in (1.1) yields t p(,t) p(t) = e t (p(,) p()) α ρ t e (t s) (pw x ) x +λ e (t s) p ( 1+c c x) x t Together with (3.4) and Lemma 3.13 and Lemma 2.1, this gives e (t s) (p(1 p) p(1 p)). p(,t) p(t) L () e t (p(,) p()) L () +α +ρ t t e (t s) (pw x ) x L () +λ e (t s) p ( 1+c c x) x L () t e (t s) (p(1 p) p(1 p)) L () t k 1 e λ1t p(,) p() L () +c 1 (1+(t s) 3 4 )e λ 1 (t s) w x L 2 () t +c 1 (1+(t s) 3 4 )e λ 1 (t s) c x L 2 () t +c 1 (1+(t s) 3 4 )e λ 1 (t s) p(1 p) p(1 p) L 2 () k 1 e λ 1t p(,) p() L () +c 2 (ǫ) +c 2 (ǫ) t t (1+(t s) 3 4 )e λ 1 (t s) e (1 ǫ)s (1+(t s) 3 4 )e λ 1 (t s) e γ(1 ǫ)s t +c 1 (1+(t s) 3 4 )e λ 1 (t s) p(1 p) p(1 p) L 2 (). It is observed that (3.37) p(1 p) p(1 p) 2 L 2 () = p(1 p) 2 L 2 () p(1 p) 2 p(1 p) 2 L 2 (). 23

24 Hence from (3.2) and Lemma 3.13, it follows that p(,t) p(t) L () k 1 e λ 1t p(,) p() L () +c 2 (ǫ) +c 2 (ǫ) +c 3 (ǫ) t t t (1+(t s) 3 4 )e λ 1 (t s) e (1 ǫ)s (1+(t s) 3 4 )e λ 1 (t s) e min{1,γ,λ}(1 ǫ)s (1+(t s) 3 4 )e λ 1 (t s) e γ(1 ǫ)s c 4 (ǫ)e min{λ1,1,γ,λ}(1 ǫ)t, (3.38) which implies (3.36). Lemma Let w > 1 1 ρ. Then for any ǫ (,min{λ 1,1,γ,λ}), there exists C(ǫ) > such that p(t) 1 C(ǫ)e (min{2λ 1,2,2γ,λ} ǫ)t. (3.39) Proof. We integrate the p-equation in the spatial variable over to obtain (p 1) t = λ(p p 2 1 (p p) 2 ) = λp(p 1) λ (3.4) p p 2 L 2 (). By (3.12), there exists t ǫ > such that p(t) 1 ǫ for t t ǫ. Hence by (3.36) and (3.4), solving the differential equation entails p(t) 1 p(t ε ) 1 e λ t tε p(s)ds + λ t e λ t s p(σ)dσ p(s) p(s) 2 L 2 () t ε p(t ε ) 1 e λ(1 ǫ)(t tǫ) + λ t e λ(1 ǫ)(t s) p(s) p(s) 2 L 2 () t ε for t t ǫ, which proves (3.39). p(t ε ) 1 e λ(1 ǫ)(t tǫ) +c 1 (ǫ) t e λ(1 ǫ)(t s) e 2min{λ 1,1,γ,λ}(1 ǫ)s p(t ε ) 1 e λ(1 ǫ)(t tǫ) +c 2 (ǫ)e min{2λ 1,2,2γ,λ}(1 ǫ)t c 3 (ǫ)e min{2λ 1,2,2γ,λ}(1 ǫ)t Lemma Let w > 1 1 ρ. Then for any ǫ (,min{λ 1,1,γ,λ}), there exists C(ǫ) > such that p(,t) 1 L () C(ǫ)e (min{λ 1,1,γ,λ} ǫ)t. (3.41) 24

25 Proof. Combining the above two lemmas, we have p(,t) 1 L () p(,t) p(t) L () + p(t) 1 C(ǫ)e (min{λ 1,1,γ,λ} ǫ)t. Proof of Theorem 1.2. (1.5) is the direct consequence of Lemma 3.6 in the previous subsection. As for (1.6) (1.8), we only need to collect (3.4), (3.26), (3.33) and (3.41). Remark 3.2. Incomparisonwith(3.41),by(3.9)and(3.29),wehavesup t w x(t) 2 C and sup t z x(t) 2 C respectively, and thus sup t p(t) W 1,2 () C. Hence an interpolation by means of the Gagliardo Nirenberg inequality in the one-dimensional setting provides p(,t) 1 L () p(,t) 1 2 W 1,2 () p(,t) L 2 () C(ǫ)e (1 2 min{1,γ,λ} ǫ)t, where we have used (3.34). Acknowledgment The authors are grateful to the referee for his illuminating comments. This work is partially supported by the NUS AcRF grant R (PYHP) and by the NNSFC grant (YW). References [1] A.R.A. Anderson, M.A.J. Chaplain, A mathematical model for capillary network formation in the absence of endothelial cell proliferation, Appl. Math. Lett., 11(1998), [2] A.R.A. Anderson, M.A.J. Chaplain, Continuous and discrete mathematical models of tumor-induced angiogenesis, Bull. Math. Biol., 6(1998), [3] N. Bellomo, A. Bellouquid, Y. Tao, M. Winkler, Toward a mathematical theory of Keller Segel models of pattern formation in biological tissues, Math. Models Methods Appl. Sci., 25(215),

26 [4] X. Cao, Global bounded solutions of the higher-dimensional Keller Segel system under smallness conditions in optimal spaces, Discrete Contin. Dyn. Syst. Ser. A 35(215), [5] M.A.J. Chaplain, G. Lolas, Mathematical modelling of cancer cell invasion of tissue: the role of the urokinase plasminogen activation system, Math. Models Methods Appl. Sci., 18(25), [6] M.A.J. Chaplain, G. Lolas, Mathematical modelling of cancer invasion of tissue: dynamic heterogeneity, Net. Hetero. Med., 1(26), [7] M.A.J. Chaplain, A.M. Stuart, A model mechanism for the chemotactic response of endothelial cells to tumor angiogenesis factor, IMA J. Math. Appl. Med. Biol., 1(1993), [8] L. Corrias, B. Perthame, H. Zaag, Global solutions of some chemotaxis and angiogenesis systems in high space dimensions, Milan J. Math., 72(24), [9] T. Hillen, K.J. Painter, M. Winkler, Convergence of a cancer invasion model to a logistic chemotaxis model, Math. Models Methods Appl. Sci., 23(213), [1] D. Horstmann, M. Winkler, Boundedness vs Blow-up in a chemotaxis system, J. Diff. Equations, 215(25), [11] D. Kinderlehrer, G. Stampacchia, An Introduction to Variational Inequalities and Their Applications, Academic Press, 198. [12] R. Kowalczyk, Z. Szymańska, On the global existence of solutions to an aggregation model, J. Math. Anal. Appl., 343 (28), [13] J. Lankeit, Y. Wang, Global existence, boundedness and stabilization in a highdimensional chemotaxis system with consumption, Discrete Contin. Dyn. Syst. Ser. B., 37(217), [14] H.A. Levine, B.D. Sleeman, M. Nilsen-Hamilton, Mathematical modeling of the onset of capillary formation initiating angiogenesis, J. Math. Biol., 42(21), [15] J. Li, Y. Wang, Repulsion effects on boundedness in the higher dimensional fully parabolic attraction repulsion chemotaxis system, J. Math. Anal. Appl., 467(218),

27 [16] Y. Li, J. Lankeit, Boundedness in a chemotaxis haptotaxis model with nonlinear diffusion, Nonlinearity, 29(216), [17] Y. Li. K. Lin, C. Mu, Boundedness and asymptotic behavior of solutions to a chemotaxis haptotaxis model in high dimensions, Appl. Math. Lett., 5(215), [18] G. Liţcanu, C. Morales-Rodrigo, Asymptotic behavior of global solutions to a model of cell invasion, Math. Models Methods Appl. Sci., 2(21), [19] G. Liţcanu, C. Morales-Rodrigo, Global solutions and asymptotic behavior for a parabolic degenerate coupled system arising from biology, Nonlinear Analysis, 72(21), [2] A. Marciniak-Czochra, M. Ptashnyk, Boundedness of solutions of a haptotaxis model, Math. Models Methods Appl. Sci., 2(21), [21] C. Morales-Rodrigo, J.I. Tello, Global existence and asymptotic behavior of a tumor angiogenesis model with chemotaxis and haptotaxis, Math. Models Methods Appl. Sci., 24(214), [22] P.Y.H. Pang, Y. Wang, Global existence of a two-dimensional chemotaxis haptotaxis model with remodeling of non-diffusible attractant, J. Diff. Equations, 263(217), [23] P.Y.H. Pang, Y. Wang, Global boundedness of solutions to a chemotaxis haptotaxis model with tissue remodeling, Math. Models Methods Appl. Sci., in press. [24] N. Paweletz, M. Knierim, Tumor related angiogenesis, Crit. Rev. Oncal. Hematol., 9(1989), [25] M.M. Porzio, V. Vespri, Hölder estimates for local solutions of some doubly nonlinear degenerate parabolic equations, J. Diff. Equations, 13(1993), [26] G. Reyes, J.L. Vázquez, A weighted symmetrization for nonlinear elliptic and parabolic equations in inhomogeneous media, J. Eur. Math. Soc., 8(26), [27] B.D. Sleeman, Mathematical modelling of tumor growth and angiogenesis, Adv. in Exp. Med. Bio., 428(1997), [28] C. Stinner, C. Surulescu, G. Meral, A multiscale model for ph-tactic invasion with time-varying carrying capacities, IMA J. Appl. Math., 8(215),

28 [29] C. Stinner, C. Surulescu, A. Uatay, Global existence for a go-or-grow multiscale model for tumor invasion with therapy, Math. Models Methods Appl. Sci., 26(216), [3] C. Stinner, C. Surulescu, M. Winkler, Global weak solutions in a PDE-ODE system modeling multiscale cancer cell invasion, SIAM J. Math. Anal., 46(214), [31] Y. Tao, Global existence for a haptotaxis model of cancer invasion with tissue remodeling, Nonlinear Anal. Real World Appl., 12(211), [32] Y. Tao, M. Winkler, Eventual smoothness and stabilization of large-data solutions in a three-dimensional chemotaxis system with consumption of chemoattractant, J. Diff. Equations, 252(212), [33] Y. Tao, M. Winkler, Energy-type estimates and global solvability in a two-dimensional chemotaxis haptotaxis model with remodeling of non-diffusible attractant, J. Diff. Equations, 257(214), [34] Y. Tao, M. Winkler, Large time behavior in a multidimensional chemotaxis haptotaxis model with slow signal diffusion, SIAM J. Math. Anal., 47(215), [35] Y. Tao, M. Winkler, A chemotaxis haptotaxis system with haptoattractant remodeling: boundedness enforced by mild saturation of signal production, Commun. Pure Anal. Appl., 18(4)(219), [36] C. Walker, G.F. Webb, Global existence of classical solutions for a haptotaxis model, SIAM J. Math. Anal., 38(27), [37] Y. Wang, Y. Ke, Large time behavior of solution to a fully parabolic chemotaxis haptotaxis model in higher dimensions, J. Diff. Equations, 26(216), [38] M. Winkler, Aggregation vs. global diffusive behavior in the higher-dimensional Keller Segel model, J. Diff. Equations, 248(21), [39] M. Winkler, Singular structure formation in a degenerate haptotaxis model involving myopic diffusion, J. Math. Pures Appl., 112(218),