cell invasion is important in many biological processes , particularly in the invasive stage of cancer and in embryonic development [ 23 , 32 ] .
there has recently been some interest in how to mathematically model cell  cell and cell  extracellular matrix adhesions in the context of cancer invasion , since the latter is believed to be characterized by a number of processes that include loss of cell  cell adhesion and enhanced cell 
matrix adhesion in addition to active migration , cell proliferation and the secretion of matrix degrading enzymes .
accurate modelling of cell adhesion is , therefore , important , but it is challenging to do so using continuous mathematical models because in such an approach one uses continuous variables for cell densities .
individual cells are not recognized as such and , therefore , there is no representation of cell boundaries . a modelling approach that aims to address these difficulties was proposed by armstrong et al . .
the ideas have been taken further and applied to various contexts by , for example , [ 3 , 12 , 15 , 23 , 31 ] .
for a cell density p(x , t ) the basic model in one spatial dimension is for x  (  ,  ) , t > 0 with initial condition p(x , 0 ) = 0(x ) .
equation ( 1 ) incorporates a direct representation of cell - cell contact via a non - local flux term , the integral term in equation ( 1 ) .
spatially structured models of tumour invasion that incorporate other aspects such as haptotaxis and age structure were considered in [ 911 ] ; see also [ 34 , 35 ] .
the purpose of the present paper is to formulate and analyse an extension of equation ( 1 ) to the case of n spatial dimensions , considering only cell - cell adhesion . again letting p(x , t )
denote the cell density at ( x , t ) , but now with x   , t > 0 and b denoting the n - dimensional ball centred at 0 and of radius  , we propose the equation with initial condition the function f(p ) models cell loss and gain , and for biological reasons we assume that f(0 ) = 0 and that there is a number p1 > 0 such that f(p ) > 0 for p  ( 0 , p1 ) , and f(p ) < 0 for p > p1 .
these assumptions imply that at lower densities , there is cell gain , but at large enough densities cell loss occurs more rapidly than the generation of new cells via division , due to the effects of crowding .
motion of cells is assumed to be partly due to fickian diffusion , which gives rise to the laplacian term .
note from simulations of their more complex cell  matrix model that , in the absence of a fickian term , and unusually for a continuous model , non - invasive growth is possible for some parameter regimes while invasive growth takes place in others .
the other component to cell motion is the movement of cells due to adhesion as modelled by the cell adhesion term in equation ( 2 ) .
it is an advective flux term and models motion of cells that is caused by mechanical forces attributable to adhesive bonds that cells have formed with their closer neighbours .
cells are able to sense their surroundings over some radius  , which is believed to be of the order of several cell diameters , due to their ability to deform and extend protrusions .
these protrusions cause adhesive bonds to form with other cells , and the resulting forces cause cell movement .
in fact , the cell adhesion term in equation ( 2 ) is an advective flux term of the form   ( p(x , t)u(x , t ) ) , where u(x , t ) is the velocity of the cells at x   at time t. if we imagine the cells as tiny spheres moving through a viscous fluid , then stokes ' law leads us to suppose that they are subject to a resistive force which is directly proportional to the velocity .
this analogy leads us to suppose that the velocity of a cell is proportional to the net adhesive force on it due to the bonds formed with nearby cells .
so , the velocity u(x , t ) can effectively also be thought of as the adhesive force on a cell due to bonds with these other nearby cells within the sensing radius  ( the range over which a cell can detect its surroundings ) . this force is taken to be of the form we think of expression ( 4 ) as a force , but it is really a velocity .
expression ( 4 ) is assumed to contain dimensional constants such as the viscosity of the medium ( perhaps embodied within the functions g and h ) , to ensure that it has the units of velocity .
viewing expression ( 4 ) as a force , we are asserting that the total force on a cell at position x is the sum , over all cells within the sensing radius  , of the local forces attributable to bonds formed with the cells at the nearby points x +  ,   b. the function h :    describes how the magnitude of the force depends on || , and is a positive function since adhesive forces are always directed towards cell centres . the vector  in front of h(|| ) , which is not present in equation ( 1 ) , gives direction to the force on cells at x due to bonds with cells at x + . another formulation would be to make this direction vector a unit vector and write ( /||)(|| ) rather than h(|| ) , as some other authors have done , see for example painter et al . .
there is no real distinction ; we are simply taking h(|| ) = ( ||)/|| .
the magnitude of the adhesive forces from cells at the nearby location x +  will depend on the number of adhesive attachments made by cells at x +  to a cell at x , and hence on the cell density at x + . this explains the term g(p(x +  , t ) ) in the integrand of expression ( 4 ) , and the function g(p ) describes how the forces depend on the local cell density .
we anticipate that g(p ) should increase linearly with p if cell density is not too great , because in this situation a doubling of the number of cells at x +  should roughly double the number of adhesive attachments made to cells at x. however , at higher densities , particularly as close cell packing is approached , the tendency to form attachments drops off and so g(p ) in fact decreases once p has risen above a certain threshold , and is zero for all p above a critical cell density p2 > 0 corresponding to close cell packing .
the shift in balance between cell division and cell loss should occur at an achievable cell density .
. there is unquestionably a need for analytic study of equations of the form ( 2 ) in spatial domains of dimension n
recent interest in this kind of equation has been largely in the areas of cancer cell invasion and wound healing .
gerisch and chaplain   and painter et al .   both consider systems containing a term of the same form as the non - local term in equation ( 2 ) , and with the same interpretation  cell velocity due to adhesive bonds with other cells  with an emphasis on two dimensions so that b becomes a disk of radius . these studies have been largely numerical , which is time - consuming in the 2d case .
in addition to their numerical simulations , gerisch and chaplain   examine the possibility of expanding the non - local term , using taylor series in the case when  is small .
this analytical approach does yield some useful insight , although the approximated equations can allow unrealistic phenomena such as blow up .
two- and three - dimensional studies are extremely important in the cancer modelling context because of the manner in which cancer cells invade tissue and the possibility of phenomena such as fingering at the invasive front ( see figure 5 in ) .
the shape of a tumour - host boundary is an important diagnostic indicator ( , and the references therein ) .
straight and/or sharp boundaries tend to imply non - invasive tumours , whereas if a tumour has a diffuse and/or wavy boundary it is likely to be invasive .
another important reason for considering equations of the form ( 2 ) in higher space dimensions is that the inclusion of another space dimension can be an important test of the stability and robustness of a one - dimensional pattern , and therefore of the plausibility of a type of mathematical model in a situation where many details are unknown .
this was an important issue in armstrong et al .   who considered a system with the same type of adhesive flux term as in equation ( 2 ) as a model for somite formation .
segmentation in a number of organisms proceeds through the formation of somites , which will eventually become repeated structures such as the vertebrae .
somitogenesis is , therefore , seen as a one - dimensional process ( along the anterior - posterior axis of the presomitic mesoderm ) and an important step in the model validation process can be to include a second spatial dimension and investigate whether the model is still capable of forming somites , as was done in . using fractional powers of the diffusion operator , we establish results about the existence and uniqueness of solutions and their positivity and boundedness in spaces of uniformly continuous functions .
we also look at solutions in l spaces and show that under suitable conditions the non - zero uniform state is asymptotically stable in l and that under more stringent conditions all solutions converge to this uniform state .
initially , we will work in the spaces in either case , take the norm ||u||k
= ||k ||du|| , where for i = 1 ,  , n , di : (di ) so di is a closed linear operator .
set d = ( d1 ,  , dn ) , with domain (d ) = (d1 )  
 (dn ) . if  = ( 1 ,  , n ) , d = d1  dn .
we define || = 1 +  + n and say    for n - tuples  and  if the inequality holds component - wise .
note that from   page 33 , u  buc0( ) implies that , for ||  k , du  0 as |x|  . we write   when results apply to either x or x0 and set x = x. we define  : ( )  x  x , by for  > 0 set a =  + i , (a ) = ( ) so that the spectrum of a is contained in the open right half plane .
we rewrite the problem ( 2)(3 ) abstractly as where p(t ) : [ 0 ,  )  x , t  0 , f :    , g :    , and with ith component ki()(x ) .
let t(t ) be the analytic semigroup in x generated by  , so and t(t ) : x0  x0 .
let tk(t ) be the analytic semigroup which is t(t ) restricted to the space   ( k an integer ) ; the corresponding generators k are given by the operator  but restricted to different domains .
it is easy to see that for   x , k > 0 , ||  k so we also have t(t ) : x  x , set ak = k + . we will use the fractional powers ak of ak , and exploit the result that , for suitable  , the operator diak : x  x is bounded , see .
we will write  = (ak ) in x , 0 = (ak ) in the space x0 and   = (ak ) in the space .
thus the mild form of ( 2)(3 ) in x is observe that for   (d ) , g  c( ) , h  l(0 ,  ) , we have k( )  ( ) , note that in our proofs of existence , we will take f , g :    but we will then prove that if 0  0 then p(t )  0 .
thus , provided there is enough smoothness at 0 it is only the properties of f and g on [ 0 ,  ) that are relevant as we can then define f and g appropriately for p < 0
. we also have [ 17 , 18 , 22 , 24 ] if  > 0 , then for every t > 0 the operator aktk(t ) is a bounded linear operator in   and for any   > 0 there exists c1 > 0 such that if 0    1 , then ak :      is a bounded linear operator so there exists a constant c2  1 such that for all u  x , if  <
 < 1 , then diak :      is a bounded linear operator so there exists a constant c3  1 such that for all u  x , the following lemma is proved in very much the same way as lemma 1 in   suppose that for some k  0 , f  c( ) , g  c( ) and h  l(0 ,  ) ,  <  < 1 and f(0 ) = 0 . then g :      and ( a)if ||ak||k  r1 
 there exists a constant k1(r1 ) such that(b)if ||i||k  r2 , i = 1 , 2 , then there exists a constant k2(r2 ) such that(c)if k  2 , and ||||k2  r3 if ||ak||k  r1 
 there exists a constant k1(r1 ) such that if ||i||k  r2 , i = 1 , 2 , then there exists a constant k2(r2 ) such that if k  2 , and ||||k2  r3
suppose that for some k  0 , f  c( ) , g  c( ) , f(0 ) = 0 , that h  l(0 ,  ) , and that 0  q  0 , where  <  < 1 , and either q = 0 or there exists p1 > 0 such that f(p1 ) = 0 and q = p1 . then the problem ( 2 ) and ( 3 ) has a unique mild solution p(t ) such that p(t )  q  c([0 , t0 ] ; 0 ) , for t0 > 0 small enough .
in addition p(t ) is the classical solution on [ 0 , t0 ] of the problem ( 2 ) and ( 3 ) in the sense that and for 0 < t
 t0 , p(t )  (ak ) , and also , under the above conditions , if k  2 m , in particular if k = 2 , p(t )  q  c([0 , t0 ] ; x0 ) and equation ( 23 ) also holds at t = 0 .
let [ 0 , 0 ) be the maximal interval of existence of the solution p(t ) .
now let pn(t ) be the solution of equation ( 2 ) with initial data pn(x , 0 ) = n(x ) where n  q  0 .
suppose that ||ak(n  0)||k  0 as n  . then , given any t1  [ 0 , 0 ) , pn(t ) is defined on [ 0 , t1 ] for large enough
p(t))||k  0 as n   , uniformly on [ 0 , t1 ] . to deal with the case
q = p1 we set p(x , t ) = w(x , t ) + p1 in equations ( 2 ) and ( 3 ) to get where f(w ) = f(w + p1 ) , and g(w ) = g(w + p1 ) .
if we set the abstract form of equation ( 25 ) is if g( ) = d  ( ( a + p1)k(a ) )  f(a ) , then it is easy to see that g satisfies lemma 1 and we can apply the same methods to the resulting equation in w. we can now obtain positivity and boundedness .
suppose that h  c [ 0 ,  ] , f  c( ) , g  c( ) , f(0 ) = 0 , and that 0  q  d0 , where  <  < 1 and q = 0 or q = p1 > 0 where f(p1 ) = 0 . let 0(x )  0 for all x  . if p(t ) is the unique classical solution of equations ( 2 ) and ( 3 ) on [ 0 , t0 ] , then p(x , t )  0 for all 0  t  t0 , x  . proof take m1 > 0 and m2  0 such that and , using the mean value theorem , choose c sufficiently large that where n is the surface area of b1 , and define u(x , t ) = p(x , t ) e , so we now prove that p(x , t )  0 on   [ 0 , t0 ] .
suppose not , then u(x , t ) will be strictly negative somewhere in this set and hence attains a global minimum at some ( x0 , t0 )    ( 0 , t0 ] .
thus for all i in the ith integral in the second term in the right - hand side of equation ( 31 ) , we may replace /xi by /i and therefore an integration by parts shows that the sum at the point ( x0 , t0 ) equals the above quantity is bounded in absolute value by evaluating equation ( 31 ) at ( x0 , t0 ) , and making use of ( 32 ) and the above bound , we obtain using inequality ( 29 ) , the fact that u(x0 , t0 ) < 0 and that c satisfies inequality ( 30 ) .
suppose that h  c[0 ,  ] , h  0 , f  c(),g  c( ) , f(0 ) = 0 and that 0  q  d0 
 where  <  < 1 , and q = 0 or q = p1 > 0 where f(p1 ) = 0
. suppose there exists p2 > p1 
 such that g(p )  0 for p  [ 0 , p2 ] , f(p2 ) < 0 and that let 0  0(x )  p2 
 for all x  . if p(t ) is the unique classical solution of equations ( 2 ) and ( 3 ) on [ 0 , t0 ] , then 0  p(x , t ) 
suppose that h  c [ 0 ,  ] , h  0 , g bounded , f  c( ) , g  c( ) , f(0 ) = 0 , and that 0  q  0 
 where  <  < 1 , and q = 0 or q = p1 > 0 where f(p1 ) = 0
. suppose there exists p3 > p1 
 such that for p > p3 , f(p ) < 0 and let 0  0(x )  p3 
 for all x  . if p(t ) is the unique classical solution of equations ( 2 ) and ( 3 ) on [ 0 , t0 ] , then 0  p(x , t ) 
the method uses proposition 2 or 3 to show first that p(t ) is bounded on bounded intervals of time .
thus by lemma 1(a ) ||g(ap(t))|| is bounded linearly , so we can use the gronwall - type inequality in   1.2.1 to show that ||ap(t)|| is suitably bounded . then we can use lemma 1(c ) to get a linear bound on ||g(a2p(t))||2 , so , using the gronwall - type inequality again , ||a2p(t))||2 is also bounded as required suppose that h  c[0 ,  ] , h  0 , f  c( ) , g  c( ) .
0 
 where  <  < 1 , and q = 0 or q = p1 > 0 where f(p1 ) = 0 and 0(x )  0 for all x  . suppose that either ( a)there exists p2 > p1 
 such that 0(x )  p2 
 for all x   , g(p )  0 for p  [ 0 , p2 ] , f(p2 ) < 0 and inequality ( 34 ) holds .
p1 
 such that 0(x )  p2 
 for all x   , g(p )  0 for p  [ 0 , p2 ] , f(p2 ) < 0 and inequality ( 34 ) holds .
or g is bounded and there exists p3 > p1 
 such that 0(x ) 
p3 
 for all x   , and for p > p3 , f(p ) < 0 and inequality ( 35 ) holds .
p <  define the operator , l : (l )  l( )  l( ) by similarly , we can also define l : (l )  l( )  l( ) , for a suitable domain (l1 ) . for 1  p <  ,  > 0 , set al = l + i , (al ) = (l ) .
it is well known   that l generates a positive analytic semigroup { tl(t)}t0 , with ||tl(t)||  1 .
further , in the case 1 < p <  , from   propositions 4.1.7 and 1.1.4 and example 1.3.9 , if  <  <  < 1 , thus , if we define di : (di )  l( )  l( ) , then inequality ( 18 ) holds in l( ) . to show that inequality ( 18 ) also holds in l( )
, it is sufficient to show that (al )  w( ) . to do this note first that it follows from   example 1.3.9 , example 1.3.11 and remark 1.3.7 that if 1 < s < 2 and 0 <  <
 < 1 , then now , by   proposition 4.1.7 ,   theorem 1.7 and prop 4.8 , and   prop 1.1.4 , for 1 < s < 2 ( the definition and properties of the besov spaces bp , q , and their relationship to the sobolev - slobodeckii spaces w , can be found in [ 1 , 16 , 20 ] . ) as before , from , tl(t ) satisfies equation ( 8) , hence if   x  l( ) , tl(t) = tk(t). equally   (al )  (al ) implies al = ak. we now have the following local existence theorem in l( ) : suppose that f  c( ) , g  c( ) , that h  l(0 ,  ) , f(0 ) = 0 and that 0  q  (al ) , where  <  < 1 , and either q = 0 or there exists p1 > 0 such that f(p1 ) = 0 and q = p1 .
then , for n  3 , the problem ( 2 ) and ( 3 ) has a unique mild solution p(t ) such that p(t )  q  c([0 , t0 ] ; (al ) ) , for t0 > 0 small enough . in addition p(t ) is the classical solution on [ 0 , t0 ] of the problem ( 2 ) and ( 3 ) in the sense that and for 0 < t  t0 , p(t )  (al ) , and also there is continuous dependence on the initial data . let [ 0 , 0 ) be the maximal interval of existence of the solution p(t ) .
now let pn(t ) be the solution of equation ( 2 ) with initial data pn(x , 0 ) = n(x ) where n  q  (al ) .
suppose that ||al(n  0)||l  0 as n  . then , given any t1  [ 0 , 0 ) , pn(t ) is defined on [ 0 , t1 ] for large enough
p(t))||l  0 as n   , uniformly on [ 0 , t1 ] .
finally , if in addition 0  q  (al ) , then also p(t )  q  c([0 , t0 ] ; (al ) ) .
the proof of the first part follows from theorems 3.3.3 and 3.4.1 in   provided we can show that g : l( )  l( ) and g : l( )  l( ) are locally lipschitz .
each product has just one term of the form dial. the rest are functions of al  d(al ) .
but , for  <  <  < 1 , and n  3 ( see   7.34 ) and the local lipschitzness follows . similarly for g. note
further that , also using inequality ( 18 ) in l (  ) , g , g : l (  )  l( )  l (  )  l( ) and are locally lipschitz continuous in l( ) . for 0  (al )  (al ) , we set up the iteration : 0 = al 0 , then the iteration converges in l( ) and l( ) and by considering the restrictions to bounded subsets of  , it can be seen that the limits are equal almost everywhere . similarly for g , and
we now consider the local asymptotic stability of the uniform state p = p1 . for simplicity , we will consider the case where f and g are logistic .
each of f and g in equation ( 2 ) can be of the form rp(1  p / k ) with different r and k. after suitable non - dimensionalization , without loss of generality , they can be taken as f(p )
= p(1  p ) and g(p ) = p(  p ) , with  > 1 .
the uniform steady state is then p = 1 . to examine its stability , set p = w + 1 , 0 = w0 + 1 to obtain the following equation for w : take n  3 .
suppose that h  l(0 ,  ) , and that 0  1  (al ) , where  <  < 1 . for    define then if k( ) < 0 for all    the zero solution of equation ( 40 ) is uniformly asymptotically stable in (al )
. precisely , if  > 0 is such that k( ) <  for all    , then there exist  > 0 , m  1 such that if ||alw0||l   , then moreover , k( ) < 0 for all    holds if conversely if there exists 0 
 such that k(0 ) > 0 , then the stationary solution is unstable .
(al )  l( ) , such that and the operator h : l( )  l( ) such that h( ) = g( )  bal. using the same reasoning as in theorem 3 it can be seen that if ||||l  r , then there exist constants k(r ) and k(r ) such that in the case of stability , we note that this implies that h( ) = o(||||l ) as ||||l  0 .
we have already seen that g : l( )  l( ) is locally lipschitz continuous .
thus theorem 5.1.1 of   holds so that if the solutions of the linearized problem are uniformly asymptotically stable then so too are those of the nonlinear problem .
the right - hand side of equation ( 44 ) is linear and generates an analytic semigroup r(t ) , say .
so , for all u0  l( ) , equation ( 44 ) has a global classical solution u(t ) = r(t)u0 .
first take u0  l( )  l( ) , so that , as in theorem 3 , u(t )  l( )  l( ) , and we may take fourier transforms .
if we denote the fourier transform of u by thus , using the plancherel identity , as required . thus by density ||r(t)u0||l
 ||u0||l2 exp(t ) for all u0  l( ) , and the result follows . the inequality ( 42 ) follows from completing the square in k( ) , while inequality ( 43 ) follows from the fact that |sin(  x)|  |  x|  |||x| . for instability , we apply corollary 5.1.6 of . from inequality ( 45 )
h( ) = o(||||l ) as ||||l  0 , so now we require that there exists  in the spectrum of the generator of r(t ) such that re  > 0 .
suppose not , so that there exists m such that , for all   l( ) , but , using the expression for ( , t ) in ( 46 ) , there exists  > 0 such that if we now choose u0 such that b|0( + 0)| d > 0 , then using the plancherel identity gives a contradiction to inequality ( 47 ) .
the next result follows in a similar fashion to theorem 3 and   theorem 3 .
suppose that f  c( ) , g  c( ) , h  c[0 ,  ] .
let 0  q  0 , where  <  < 1 , and either q = 0 or q = p1 > 0 with f(p1 ) = 0 . suppose also that the hypotheses of proposition 1 and of either proposition 2 or proposition 3 on f , g , h , and 0 
 hold . take also 0  q  (al ) .
q  c([0 ,  ) ; 0 ) be the unique classical solution of equations ( 2 ) and ( 3 ) .
then , for each t  0 , p(t )  q   ( al ) = w( ) . furthermore p(t )  q  c((0 ,  ) ; l( ) ) .
let h  c[0 ,  ] be such that h  0 , and let f(p ) = p(1  p ) and g(p ) = p(  p ) with  > 1 . let 0  1  0  (al(n ) ) , where  <  < 1 , and let p(x , t ) be the solution of the problem ( 2 ) and ( 3 ) .
( a)if 0  0(x )   for all x   , where  satisfiesthen 0  p(x , t )   for all x   , t  0.(b)if   0(x )   for all x   , where 0 < 
  p(x , t )   for all x   , t  0.(c)if   0(x ) <  for all x   , where  > 0 ,  satisfies inequalities ( 48 ) , ( 49 ) , andthen p(x , t )  1 exponentially as t   in the sense that for all t > 0 if 0  0(x )   for all x   , where  satisfies then 0  p(x , t )   for all x   , t  0 . if   0(x )   for all x   , where 0 < 
<  ,  satisfies inequality ( 48 ) , and then   p(x , t )   for all x   , t  0 . if   0(x ) <  for all x   , where  > 0 ,  satisfies inequalities ( 48 ) , ( 49 ) , and then p(x , t )  1 exponentially as t   in the sense that for all t > 0 first , to prove ( a ) , note that inequality ( 48 ) implies inequality ( 34 ) with p2 =  , f(p ) = p(1
 p ) , and thus , solutions remain in the closed interval [ 0 ,  ] for all positive times .
next , to prove ( b ) , suppose that there exists t * > 0 and a corresponding x * with p(x * , t * ) =  , p(x * , t*)/xi = 0 , p(x * , t*)/xi  0 , then note that 0  g(p )  /4 when p  [ 0 ,  ] . therefore by inequality ( 49 ) .
hence p(x , t )   for t > 0 and ( b ) is proved .
last , to prove ( c ) , define w(x , t ) by p(x , t ) = 1 + w(x , t ) , so that w(x , t ) satisfies equation ( 40 ) and   1  w    1 .
multiplying equation ( 40 ) by w(x , t ) and integrating with respect to x over  , and then integrating by parts on the laplacian term and the /xi term ( this is justified by theorem 5 ) , the integrand of the middle term in the right - hand side will now be estimated by using the inequality a  b
this bounds the integral by a sum of two integrals , one of which cancels with the first term on the right - hand side to give since  > 1 , it is easy to show that recalling that 0  1 + w(x , t )   , so using inequality ( 54 ) and theorem iv.15 from , therefore inequality ( 53 ) becomes note that w(x , t )  ( 1  )w(x , t ) , so that inequality ( 56 ) becomes and ( c ) follows from inequality ( 50 ) .
we illustrate in figures 14 the results of the previous two sections with simulations that can be applied to in vitro wound closure experiments . in these experiments ,
cell cultures are scored in a thin line that closes as the cell population proliferates .
a review of mathematical models of wound healing is given in   and other references are found in [ 57 , 13 , 14 , 21 , 2529 , 33 , 36 ] . in previous work , we modelled in vitro wound healing experiments for the model ( 2 ) with n = 1 , and examined the dependence of solutions on the diffusion parameter  and the sensing radius . in   numerical simulations demonstrated that for larger values of the diffusion parameter  and smaller values of the sensing radius  the wound closed completely across the wound opening , but for smaller values of  and larger values of  complex patterns arose across the wound opening with incomplete closing .
either of the alternative conditions ( 42 ) or ( 43 ) is sufficient for the stability of the uniform steady state p = 1 , and note that both of these conditions are largeness conditions on  and smallness conditions on both h( ) and on the sensing radius . similarly conditions ( 48 ) , ( 49 ) and ( 50 ) are sufficient for convergence to the uniform steady state and conditions ( 48 ) and ( 49 ) are smallness conditions on h( ) and  while ( 50 ) is in addition a largeness condition on . in our simulations here , we examine the dependence of the wound healing model with respect to the parameter  , which is the non - dimensionalized cell density corresponding to close cell packing , in equation ( 2 ) with g(p ) = p(
 = 0.1 ,  = 0.2 , f(p ) = p(1.0  p ) , g(p ) = p(  p ) , h(x ) = arctan ( 100.0x)/(10.0x arctan 100.0 ) and initial data 0(x ) = 1.0  0.8 exp((0.1x ) ) . the red plane and lines correspond to  = 0.2 and the green plane and lines correspond to 0.0 .
( a ) graph of p(x , t ) , 0  t < 2.0 . the wound does not close , and instead the solution exhibits an evolving pattern of peaks and valleys which do not converge to 1.0 . ( b ) graphs of p(x , .0 ) , p(x , .5 ) , p(x , 1.0 ) , p(x , 2.0 ) .
simulation with  = 1.0 ,  = 1.5 ( a ) ,  = 17.0 ( b ) ,  = 18.0 ( c ) ,  = 19.0 ( d ) ,  = 1.0 ,  = 0.2 , f(p ) = p(1.0
 p ) , g(p ) = p(  p ) , h(x ) = arctan(100.0x)/(9.0xarctan 100.0 ) , and initial data 0(x ) = 1.0  0.8 exp((0.1x ) ) . for the higher values of  ,
the wound healing simulation wave patterns as in figure 2 with adhesion strength parameters ( a )  = 17.0 , ( b )  = 18.0 , and ( c )  = 19.0 . for  = 17.0
, the solution is shown at times t = 100.0 ( blue ) , t = 95.0 ( red ) and t = 90.0 ( green ) .
the amplitude envelope of the oscillation packet decreases to a minute value for the times shown . for  = 18.0
the solution is shown at times t = 200.0 ( blue ) , t = 150.0 ( red ) and t = 100.0 ( green ) .
for  = 19.0 , the solution is shown at times t = 150.0 ( blue ) , t = 125.0 ( red ) , t = 100.0 ( green ) , t = 75.0 ( yellow ) and t = 50.0 ( cyan ) .
the simulations indicate convergence to the uniform equilibrium  1.0 in ( a )  < 17.62 , but not in ( b ) and ( c ) ,  > 17.62 . wound healing simulations with parameters  = 2.0 ,  = 1.0 ,  = 3.0 , f(p ) = p(1.0
 p ) , h(x ) = 9.0 arctan(100.0x)/(4.0x arctan 100.0 ) , and initial data 0(x ) = 1.0  0.4 exp((0.1x ) ) ( top left ) , 0(x ) = 1.0  0.95 exp((3x ) ) ( top right ) , 0(x ) = 1.0  0.5 exp((0.2x ) ) ( bottom left ) and 0(x ) = 1.0  0.95 exp((0.1x ) ) ( bottom right ) . the wounds top
left and top right both close but the wound bottom right which has the same width as the first and the same depth as the second fails to close .
the widest wound , bottom left , also closes . in figure 1 , n = 1 ,
 p ) ,  = 0.1 ,  = 0.2 ,  = 1.0 , h(x ) = arctan(100.0x)/(10.0x arctan 100.0 ) , and 0(x ) = 1.0  0.8 exp((0.1x ) ) . in this case ,
theorem 6 ( 48 ) is satisfied if  < 18.94 , ( 49 ) is satisfied if  < 1.457 , and ( 50 ) is satisfied if  < 2.823 .
we take  = 19.0 , so that none of the conditions ( 48 ) , ( 49 ) , ( 50 ) of theorem 6 are satisfied .
the simulation in figure 1 shows that for this value of  , the solution does not remain bounded below by  and does not converge to the equilibrium 1.0 . on the other hand ,
if condition ( 48 ) were satisfied , then the solution p would be bounded ( 0  p   ) , and theorem 2 would give global existence of the solution .
indeed , if we take  = 18 , then condition ( 48 ) holds , and we find that the solution behaves qualitatively the same as in figure 1 .
this provides an example where the theory gives global existence of the solution , but the solution does not converge to the equilibrium but instead forms a series of peaks . in figures 2 and 3 , n = 1 ,
 p ) ,  = 1.0 ,  = 0.2 , h(x ) = arctan(100.0x)/(9.0x arctan 100.0 ) , and 0(x ) = 1.0  0.8 exp((0.1x ) ) .
we consider different values of  to illustrate the sensitivity of this parameter for the convergence of the solutions to the equilibrium  1.0 .
the hypotheses of theorem 6 require that  < 16.94 ( 48 ) ,  < 2.683 ( 49 ) , and  < 2.597 ( 50 ) . for  < 2.597
all three conditions are satisfied , and theorem 6 implies the solution remains in the interval [ 0.2 ,  ] and converges to 1.0 for all initial data such that 0.2  0(x )  . on the other hand if one looks at the stability condition from theorem 4 , then the equilibrium  1.0 is asymptotically stable for  < 17.62 with instability if  > 17.62 . in figure 2 ( a )
 = 1.5 and there is convergence to 1.0 ( the solution appears to converge monotonically to 1.0 from below , never rising above 1.0 ) . in figure 2(b)(d )  = 17.0 , 18.0 , 19.0 , respectively , and the solutions develop a series of peaks , symmetric to the right and left of 0.0 , and rising above and falling below 1.0
. further examination of this wave behaviour is given in figure 3 for the cases  = 17.0 ( a ) ,  = 18.0 ( b ) , and  = 19.0 ( c ) .
for  = 17 , the amplitude of the waves decreases as t increases and the solution appears to converge slowly to 1.0 . for  = 18 and  = 19 ,
the wave propagates with amplitude dependent on  ( the greater the value of  , the greater the amplitude ) .
in general , for the roles of the model parameters , we postulate that cells concentrate in interior and exterior regions of the wound in a regular pattern depending on the strength of adhesion  , the sensing radius  , and the diffusion coefficient . in figure 4 , we show that the convergence to the equilibrium  1.0 corresponding to the healed wound depends on the initial conditions . in these simulations ,
 = 3.0 , f(p ) = p(1.0  p ) , g(p ) = p(
 p ) , and h(x ) = 9.0 arctan(100.0x)/(4.0x arctan 100.0 ) . in this figure ,
the only change in the simulations is the initial data  the parameters are the same for all . for initial conditions corresponding to a shallow , wide wound , and to a narrow , deep wound ,
closure occurs , but for an initial condition corresponding to a wound with the same width as the first and the same depth as the second , closure does not occur .
we note that the function k( ) defined in equation ( 41 ) is negative for all    , which by theorem 4 implies that the equilibrium  1.0 is stable in this example .
we have chosen to illustrate the behaviour of solutions for the model using the simple case of an equilibrium corresponding to wound healing experiments , in which wound closure corresponds to convergence to a constant function in the x - direction , perpendicular to the direction of the scoring .
 a model for cell  cell adhesion , based on an equation originally proposed by armstrong et al . 
 [ a continuum approach to modelling cell  cell adhesion , j. theor . 
 biol . 
 243 ( 2006 ) , pp . 
 98113 ] , is considered . 
 the model consists of a nonlinear partial differential equation for the cell density in an n - dimensional infinite domain . 
 it has a non - local flux term which models the component of cell motion attributable to cells having formed bonds with other nearby cells . 
 using the theory of fractional powers of analytic semigroup generators and working in spaces with bounded uniformly continuous derivatives , the local existence of classical solutions is proved . 
 positivity and boundedness of solutions is then established , leading to global existence of solutions . 
 finally , the asymptotic behaviour of solutions about the spatially uniform state is considered . 
 the model is illustrated by simulations that can be applied to in vitro wound closure experiments.ams classifications : 35a01 ; 35b09 ; 35b40 ; 35k57 ; 92c17 