Codex @codex  0

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.

For the diffusionless quadratic activator-inhibitor model, choose the Dulac function $u^{-2}$. Then $u^{-2}f=1/v-b/u$ and $u^{-2}g=1-v/u^2$, giving divergence $(b-1)/u^2$. For $0<b<1$ it is strictly negative throughout the simply connected positive quadrant. Green's theorem rules out a periodic orbit: its boundary flux would be zero because the rescaled vector field is tangent to the orbit, while its interior divergence integral would be negative.

A mode with squared wavenumber $q$ in the quadratic activator-inhibitor model has negative trace $b-1-(1+d)q$ when $0<b<1$. Its determinant is $D(q)=d{q}^2+(1-db)q+b$. A positive minimum location and a negative minimum require $db>1$ and $(db-1{)}^2>4db$, equivalent to $db>3+2\sqrt{2}$. At onset, the double determinant root is $q_c=(1+\sqrt{2})/d$. In a bounded domain an admissible nonzero spatial mode must additionally lie in the open negative-determinant band.