{
 "artifacts": null,
 "category": "dynamics",
 "description": "Replace a generic recurrent update by a positive-state continuous-time cell whose interactions are restricted to a quadratic zero-one reaction-network motif with three state variables and six reactions. Select a motif known to possess three positive equilibria, then use the two stable equilibria as binary memory states and the intervening unstable equilibrium as the separatrix. This creates an explicitly multistable RNN module with a bounded attractor count and a measurable stability certificate.",
 "formulas_latex": [
  "$$\\dot{x}=f(x)=N h(x)+u(t),\\qquad h_j(x)=k_j\\prod_{i=1}^{3}x_i^{Y_{ij}},\\quad k_j\u003e0,\\quad Y_{ij}\\in\\{0,1\\},\\quad \\sum_iY_{ij}\\le2.$$",
  "$$J(x)=\\frac{\\partial f}{\\partial x}=N\\,\\operatorname{diag}(h(x))Y^T\\operatorname{diag}(x)^{-1},\\qquad J(z_l)v=\\lambda v.$$",
  "$$\\operatorname{sign}\\det J(z_{l+1})=-\\operatorname{sign}\\det J(z_l),\\qquad l=1,2.$$",
  "$$p_l(\\lambda)=\\det(\\lambda I-J(z_l))=\\lambda^3+a_{1,l}\\lambda^2+a_{2,l}\\lambda+a_{3,l};\\quad z_l\\text{ stable}\\iff a_{1,l}\u003e0,\\ a_{2,l}\u003e0,\\ a_{3,l}\u003e0,\\ a_{1,l}a_{2,l}\u003ea_{3,l}.$$"
 ],
 "id": 1655,
 "implementation": "1. Exact integration point: implement this as a continuous-time recurrent layer inside an RNN or latent world model. Choose one of the paper's 373 three-species, six-reaction motifs, or enumerate candidate integer matrices N and Y and retain only motifs with two stable and one unstable positive equilibria. Parameterize k_j=softplus(theta_j)+epsilon and x_i=softplus(q_i)+epsilon so positivity is preserved. Use an input encoder to produce u(t), with a small leak or bounded input clamp to prevent input-driven escape. Decode the state with a linear layer or small MLP.\n\n2. Pseudocode: initialize x=softplus(q). At every solver step compute h[j]=k[j]*product(x[i]**Y[i,j]), compute dx=N@h+u, and integrate with RK4 or an adaptive ODE solver. During training, periodically solve f(x)=0 from many positive initializations, compute J=N@diag(h)@Y.T@diag(1/x), calculate its eigenvalues and Hurwitz margins, and add a penalty when the two intended memory equilibria violate the Hurwitz inequalities. Train the readout and rates jointly, but freeze the motif topology or perform a discrete motif search before gradient training.\n\n3. Computed from the paper: the allowed quadratic zero-one topology, the six-reaction/full-rank restriction, the alternating determinant-sign rule, and the maximum of three positive equilibria. Estimated empirically: equilibrium locations, basin volumes, eigenvalues, and the parameter value at which a stable/unstable pair disappears. Do not use determinant sign as the sole stability test; verify all eigenvalues or the cubic Hurwitz conditions.\n\n4. First cheap experiment: compare this cell with a three-unit tanh RNN on a two-bit associative-memory task and a delayed XOR sequence task. Use identical parameter counts and train on short sequences, then test recall after long free-running horizons under small perturbations. Sweep one rate multiplier k_1=c k_1^0 and use multistart root finding to measure the number of positive equilibria. The predicted signature is a stable-unstable-stable ordering with opposite signs of det(J) at consecutive roots, at most three positive roots, and a sharp bistability transition at a saddle-node: the smallest Jacobian eigenvalue should approach zero from below at the boundary, while basin occupancy changes sharply. A successful implementation should reproduce the equilibrium pattern across at least 90% of random positive initializations and locate the transition within 20% of a high-resolution parameter sweep.",
 "math_summary": "Let x(t)=(x_1,x_2,x_3)^T in R_{\u003e0}^3 be the cell state, let N in Z^{3x6} be the stoichiometric matrix, Y in {0,1}^{3x6} the reactant matrix, and k_j\u003e0 the six reaction rates. For reaction j, define the mass-action monomial h_j(x)=product_i x_i^{Y_{ij}}, with sum_i Y_{ij}\u003c=2, and let the neural dynamics be f(x)=N h(x)+u(t), where h=(k_1 h_1,...,k_6 h_6)^T and u is a bounded input injection. At a positive equilibrium z_l, f(z_l)=0 when u=0. The Jacobian is J(z_l)=N diag(h(z_l)) Y^T diag(z_l)^{-1}. For the paper's full-dimensional quadratic zero-one class, consecutive positive steady states satisfy sign det J(z_{l+1})=-sign det J(z_l), and the three-state motif has at most three positive equilibria; the bistable case therefore has stable-unstable-stable equilibria. Stability must be checked with the full cubic Hurwitz conditions, not determinant alone: for p(lambda)=lambda^3+a_1 lambda^2+a_2 lambda+a_3, asymptotic stability requires a_1\u003e0, a_2\u003e0, a_3\u003e0, and a_1 a_2\u003ea_3.",
 "math_tags": [
  "dynamical-systems",
  "bifurcations",
  "linear-algebra",
  "control-theory"
 ],
 "ml_areas": [
  "rnn",
  "memory",
  "world-model"
 ],
 "paper": {
  "arxiv_id": "2608.01116",
  "arxiv_url": "https://arxiv.org/abs/2608.01116",
  "summary_what_math_gives_to_ml": "The paper provides a constructive multistability mechanism in three-dimensional quadratic zero-one mass-action networks: among minimal six-reaction motifs, positive steady states occur in an alternating stability pattern, with consecutive Jacobian determinants having opposite signs and exactly half of the states stable. It also establishes a sharp maximum of three positive steady states for the (3,6,3) class, so the minimal bistable motif has two stable equilibria separated by one unstable equilibrium. A transferable neural-network construction is a constrained continuous-time recurrent cell whose vector field is assembled from such quadratic nonnegative reaction terms, yielding certified attractor-based memory rather than relying on a generic RNN to learn a suitable phase portrait.",
  "title": "The Ubiquity of Three Steady States: Minimal Multistable Zero-One Reaction Networks",
  "year": "2026"
 },
 "ratings": {
  "difficulty": 6,
  "novelty": 7,
  "usefulness": 6
 },
 "solves": [
  "stability",
  "memory",
  "accuracy"
 ],
 "title": "Three-Equilibrium Reaction RNN Memory Cell",
 "url": "https://synthcore.org/idea/1655/three-equilibrium-reaction-rnn-memory-cell",
 "verification": {
  "peer_reviewed": false,
  "status": "unverified",
  "status_label": "Unverified",
  "verdict_source": "deterministic test code (paired-seed permutation statistics)",
  "verification_axes": {
   "benchmark_mechanism": {
    "confirmed": null,
    "tested": false
   },
   "practical_benchmark": {
    "beats_baseline": null,
    "tested": false
   },
   "toy_mechanism_gate": {
    "confirmed": null,
    "tested": false
   }
  }
 }
}