Paper 3, Section II, B

Mathematical Biology
Part II, 2011

The number density of a population of amoebae is n(x,t)n(\mathbf{x}, t). The amoebae exhibit chemotaxis and are attracted to high concentrations of a chemical which has concentration a(x,t)a(\mathbf{x}, t). The equations governing nn and aa are

nt=αn(n02n2)+2n(χ(n)na)at=βnγa+D2a\begin{aligned} &\frac{\partial n}{\partial t}=\alpha n\left(n_{0}^{2}-n^{2}\right)+\nabla^{2} n-\nabla \cdot(\chi(n) n \nabla a) \\ &\frac{\partial a}{\partial t}=\beta n-\gamma a+D \nabla^{2} a \end{aligned}

where the constants n0,α,β,γn_{0}, \alpha, \beta, \gamma and DD are all positive.

(i) Give a biological interpretation of each term in these equations and discuss the sign of χ(n)\chi(n).

(ii) Show that there is a non-trivial (i.e. a0,n0a \neq 0, n \neq 0 ) steady-state solution for nn and aa, independent of x\mathbf{x}, and show further that it is stable to small disturbances that are also independent of x\mathbf{x}.

(iii) Consider small spatially varying disturbances to the steady state, with spatial structure such that 2ψ=k2ψ\nabla^{2} \psi=-k^{2} \psi, where ψ\psi is any disturbance quantity. Show that if such disturbances also satisfy ψ/t=pψ\partial \psi / \partial t=p \psi, where pp is a constant, then pp satisfies a quadratic equation, to be derived. By considering the conditions required for p=0p=0 to be a possible solution of this quadratic equation, or otherwise, deduce that instability is possible if

βχ0n0>2αn02D+γ+2(2Dαn02γ)1/2\beta \chi_{0} n_{0}>2 \alpha n_{0}^{2} D+\gamma+2\left(2 D \alpha n_{0}^{2} \gamma\right)^{1 / 2}

where χ0=χ(n0)\chi_{0}=\chi\left(n_{0}\right).

(iv) Explain briefly how your conclusions might change if an additional geometric constraint implied that k2>k02k^{2}>k_{0}^{2}, where k0k_{0} is a given constant.