TOPICS
Search

Schnakenberg System


The Schnakenberg system is a planar mass-action system for two chemical species (Schnakenberg 1979). A dimensionless form for the concentrations x and y is

x^'=b-x+x^2y
(1)
y^'=a-x^2y,
(2)

where a,b>0 and primes denote derivatives of x(t) and y(t) with respect to time t. The opposite signs of x^2y represent conversion of one unit of y into one of x by autocatalysis. Adding the two equations gives (x+y)^'=a+b-x. Writing s=a+b, the unique positive equilibrium point is (s,a/s^2). The Jacobian at this point has determinant s^2 and trace -1+2a/s-s^2, so its two eigenvalues are purely imaginary when the trace vanishes. If a=mus and b=(1-mu)s with 0<s<1, varying mu through (1+s^2)/2 makes the pair of eigenvalues cross the imaginary axis. The resulting supercritical Hopf bifurcation produces a small attracting limit cycle on the side where mu>(1+s^2)/2. Adding diffusion terms gives a spatially dependent system of partial differential equations for the two concentrations.


See also

Brusselator Equations, Mass-Action System

Explore with Wolfram|Alpha

References

Schnakenberg, J. "Simple Chemical Reaction Systems with Limit Cycle Behaviour." J. Theoret. Biol. 81, 389-400, 1979. https://doi.org/10.1016/0022-5193(79)90042-0.

Cite this as:

Weisstein, Eric W. "Schnakenberg System." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/SchnakenbergSystem.html

Subject classifications