 Solution

At a positive spatially homogeneous equilibrium point, $v=u^2$ and $2+u-v=0$. Hence $(u,v)=(2,4)$. The reaction Jacobian matrix there isThe two-dimensional Routh-Hurwitz stability criterion gives strict linear asymptotic stability for $a>2$, and instability for $0<a<2$.

For a spatial mode with $s=k^2$, the reaction-diffusion linearization is $J_s=J-\mathrm{diag}(s,ds)$. When $a>2$, its trace is negative for every $s\ge 0$, whileA growing mode exists precisely when this quadratic is negative at some $s>0$. Its minimum must occur at positive $s$, and must be strictly negative, givingWith $z=d/a$, the latter inequality becomes $4z^2-28z+1>0$. Selecting the root compatible with $2d>a$ yields the Turing instability thresholdThe unstable band isOn a bounded domain a boundary condition must permit a mode in this band. Equality at the threshold gives a neutral mode, not positive growth.

At $a=2$, the homogeneous linearization has purely imaginary eigenvalues, so it is not strictly linearly asymptotically stable. If “stable” is read as nonlinear stability, the boundary deserves attention. Put $x=u-2$, $y=(v-4-x)/\sqrt{3}$, and $x=r\cos \theta$, $y=r\sin \theta$. The exact equations give $\dot{r}=r^{2}F(\theta )$, $\dot{\theta }=\omega +rG(\theta )$, where $\omega =2\sqrt{3}$,Expanding $dr/d\theta$ through cubic order, ${\int }_0^{2\pi }F=0$ cancels the quadratic displacement. The correction from the quadratic change in $r$ also integrates to zero, because it is proportional to $(\int F{)}^2$. Since ${\int }_0^{2\pi }FG=\pi /(2\sqrt{3})$, the Poincaré return map isThis cubic return-map stability of a weak focus shows that $a=2$ is a weak attracting focus. The strict linear range used in the Turing instability calculation is $a>2$; the positive homogeneous equilibrium is also nonlinearly asymptotically stable at the boundary $a=2$. For that boundary, $\mathrm{tr}{J}_s=-(1+d)s<0$ for every nonzero spatial mode. The same strict Turing instability threshold and unstable band therefore also apply there, with the homogeneous mode understood through nonlinear stability.