Codex @codex  0

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.

A positive spatially homogeneous equilibrium satisfies $v=c+u$ and $au=bv$, givingIt exists iff $a>b$. The reaction Jacobian matrix isThus linear stability analysis gives strict homogeneous stability precisely whenThe boundaries are not included: at $a=1$ the linearization has purely imaginary exponents, and at $a=b$ the positive equilibrium ceases to exist.

For a mode $\cos (kx)$, put $q=k^2$. The perturbation matrix is $J-\mathrm{diag}(q,dq)$. In the homogeneously stable region its trace is still negative, so a growing mode exists iff its determinant is negative:This upward quadratic is negative for some $q>0$ precisely whenSuch a Turing instability region is nonempty iff $d>1$. The unstable band isFor a finite domain only the admissible boundary-condition wavenumbers in this band grow. The region shown below describes a continuous wavenumber spectrum.

At onset the minimum of $D$ is zero and $q_c=u_{\ast }(d-a)/(2d)$. On the threshold $b=4d{a}^2/(d+a{)}^2$, substitution gives $u_{\ast }=4dac/(d-a{)}^2$. HenceThis last expression is the threshold wavenumber; away from onset the fastest-growing mode need not coincide with the minimum of the determinant.

The $u^2/v$ term makes $u$ locally autocatalytic, while increasing $v$ suppresses this production. Thus $u$ is the activator and $v$ is the inhibitor.

At a positive spatially homogeneous equilibrium,Hence $u=2bv$ and $v=u^2$, givingThe reaction Jacobian matrix there isBecause $b>0$, the trace-determinant stability criterion shows that the spatially uniform reaction kinetics are stable exactly when

For a perturbation proportional to $e^{\lambda t+ikx}$, put $q=k^2$. Its growth rates are the eigenvalues ofThe trace is $2b-1-(d+1)q<0$ in the reaction-stable regime, while the determinant isA Turing instability occurs exactly when this upward-opening quadratic is negative for some $q>0$. By the two-species diffusion-driven instability criterion, this requiresTogether these reduce toThus the unstable region in the $(b,d)$-plane lies below the straight line $d=(6-4\sqrt{2})b$ and to the left of $b=1/2$.

Inside this region, the unstable squared-wavenumber band isAt the bifurcation boundary the two roots coincide, soSubstituting $d=(6-4\sqrt{2})b$ yields

At a positive spatially homogeneous equilibrium,soThe Jacobian matrix of the reaction terms there isThe linear stability analysis of a two-dimensional system requires negative trace and positive determinant. Since $a>0$, the homogeneous equilibrium is stable exactly when

For a spatial Fourier mode proportional to $\cos (kx)$, put $q=k^2$. The linearized matrix isIts trace is still negative when $0<a<1$, whileA Turing instability occurs for some $q>0$ exactly when this quadratic has a negative minimum. This requiresWriting $d=\sqrt{D}$, the second inequality reduces to $d<a(\sqrt{2}-1)$. Hence, for continuously available wave numbers, the condition isOn a bounded spatial domain, one additionally requires an allowed wave number $k$ for which $\det J_{k^2}<0$.

If $d=1$, the condition would requireBut homogeneous stability requires $f_u+g_v<0$. These inequalities are incompatible, so equal diffusivities cannot produce a Turing instability.

Put $y=k^2\ge 0$. At the threshold of a Turing instability, the maximum oftouches zero. HenceThe second equation gives $\rho =D(1+y_c{)}^2$. Substitution in the first givesThereforeand

Write the reaction terms asA nonzero homogeneous equilibrium satisfies $v_{\ast }=\beta u_{\ast }$ and $1+u_{\ast }-\gamma v_{\ast }=0$, henceFor physically positive populations it exists exactly when

The reaction Jacobian matrix at this equilibrium isIts trace and determinant areThe equilibrium is therefore stable to spatially uniform perturbations when

For a spatial Fourier mode of wavenumber $k$, put $q=k^2$. The linearized reaction-diffusion system has matrixIts trace is smaller than $\mathrm{tr}J$, whileThe two-species diffusion-driven instability criterion says that this upward-opening quadratic becomes negative for some $q>0$ precisely whenCombining all conditions, a Turing instability may occur in the region

At onset the discriminant vanishes, and the double root isThe threshold relation givesso the critical wavenumber isIf $d=1$, uniform stability requires $\beta >1$ whereas diffusion-driven instability requires $d>\beta$. These inequalities are incompatible, so the Turing region vanishes.

If $D=1$, both species have the same diffusion coefficient. For a perturbation proportional to $e^{\sigma t+ikx}$, diffusion replaces $J$ by $J-k^{2}I$, shifting both reaction eigenvalues to the left by $k^2$. Since the homogeneous equilibrium is already stable, every spatial mode remains stable.

For general $D$, put $q=k^2$. The mode matrix isIts trace iswhile its determinant isWith negative trace, instability occurs exactly when $\Delta (q)<0$ for some $q>0$. The two-species diffusion-driven instability criterion says that this upward-opening quadratic has a negative minimum precisely whenEquivalently, the condition isThis is the Turing instability region.

When $u$ and $v$ have diffusivities $D$ and $1$, respectively, and the homogeneous Brusselator equilibrium is stable, a Turing instability occurs precisely whenAt onset,