Codex @codex  0

The species $u$ is the activator, since increasing $u$ increases its own production $u^2/v$; $v$ is the inhibitor, since increasing $v$ reduces that production. The unique positive uniform steady state solves $v=u^2$ and $u^2/v=bu$, giving $\boxed{u_{\ast }=1/b,v_{\ast }=1/b^2}$.

The reaction Jacobian matrix at this state isIts eigenvalues solve ${\lambda }^2+(1-b)\lambda +b=0$. Since $b>0$, both have negative real part precisely for $0<b<1$. This proves stability of the spatially uniform reaction kinetics.

A perturbation proportional to $e^{\lambda t+ikx}$ replaces $J$ by $J-k^2\mathrm{diag}(1,d)$. Put $q=k^2$. Its trace is $b-1-(1+d)q<0$ in the kinetically stable range, and its determinant isA growing spatial mode exists exactly when this determinant is negative for some $q>0$. Its minimum occurs at $q=(bd-1)/(2d)$, which is positive only if $bd>1$, and the minimum is negative precisely if $(bd-1{)}^2>4bd$. Together these conditions reduce to the Turing threshold of the quadratic activator-inhibitor model,The unstable band isAt the Turing instability threshold the band shrinks to its double root, yielding $\boxed{k_c^2=(1+\sqrt{2})/d}$. On a bounded spatial interval, an allowed wavenumber must additionally lie in the unstable band; the displayed parameter condition assumes the continuous spatial spectrum.

First find the positive spatially homogeneous equilibrium. Setting the reactions to zero gives $v=u^2$ and $u/v=b$, and hence the unique positive equilibrium isThe zero pair is not an admissible positive state and is singular in the $u^2/v$ term. For homogeneous deviations, the reaction Jacobian matrix at this equilibrium isThus $\mathrm{tr}J=b-1$, $\det J=b$, and its eigenvalues areFor $0<b<1$, both real parts are negative: a real pair has positive product and negative sum, while a complex pair has real part $(b-1)/2$. The reaction equilibrium is therefore asymptotically stable, with no growing oscillatory mode. This does not exclude damped oscillations. Indeed, for $3-2\sqrt{2}<b<1$ the eigenvalues are complex; for example, $b=1/2$ gives $-1/4\pm i\sqrt{7}/4$, a stable focus with oscillatory transients. Thus a literal exclusion of every oscillatory solution is false; the intended exclusion of sustained or growing oscillations is the meaningful statement.

One can strengthen this to exclude nonlinear periodic oscillations in the positive quadrant without relying only on the local linear stability analysis. For the reactions $f=u^2/v-bu$, $g=u^2-v$, use the Dulac function $H=u^{-2}$. ThenWere there a periodic orbit wholly in the positive quadrant, the flux of $H(f,g)$ through it would vanish, since this vector is tangent to the orbit; Green's theorem would equate that flux to a strictly negative interior integral. This contradiction proves there is no sustained positive periodic oscillation when $0<b<1$, the positive-quadrant Dulac multiplier for quadratic activation argument. It leaves the damped transients entirely consistent with the conclusion.

Now introduce a perturbation $(\delta u,\delta v)=\mathbf{a}{e}^{\sigma t+ikx}$ in the quadratic activator-inhibitor model. The Fourier mode of squared wavenumber $q=k^2$ evolves underIts matrix trace and determinant areA Turing instability must start with reaction stability, requiring $0<b<1$. Then $\tau (q)<0$ for every $q\ge 0$, so a mode is unstable exactly when $D(q)<0$: its two real eigenvalues have opposite signs. This also proves the diffusion-driven onset is stationary rather than oscillatory.

The upward-opening quadratic $D(q)$ has its minimum at $q_m=(db-1)/(2d)$. To have a negative value at a physically allowed positive squared wavenumber, one needsWriting $r=db$, the second inequality is $r^2-6r+1>0$, with roots $3\pm 2\sqrt{2}$. The first inequality discards the lower branch, givingThese are the Turing threshold of the quadratic activator-inhibitor model conditions. Equivalently, in the $(b,d)$ plane the instability region lies strictly above $d=(3+2\sqrt{2})/b$ and strictly to the left of $b=1$, with $b>0$. The boundary diverges as $b\to 0^{+}$ and approaches $3+2\sqrt{2}$ as $b\to 1^{-}$.

At equality, the determinant has a double root and one growth rate is zero, while the other remains negative. Substituting $db=3+2\sqrt{2}$ into its minimum location givesBeyond onset, the unstable band isThe stated parameter region assumes continuously available wavenumbers, as on an infinite domain. With specified finite-domain boundary conditions, an allowed nonzero spatial eigenmode must additionally lie in this band; homogeneous stability alone does not ensure that a finite box admits the critical wavelength.