Multistability and Hysteresis

Topological Determinism and Biological Switches in Chemical Networks

Far-from-equilibrium chemical networks governed by non-linear feedback mechanisms frequently exhibit multistability, a dynamical regime where a single parameter set supports multiple simultaneously stable non-equilibrium stationary states. Within the framework of Chemical Reaction Network Theory, the topological structure of the reaction manifold—specifically its deficiency and the presence of complex stoichiometric cycles—dictates the network’s capacity for multiple steady-state behaviors, independent of the magnitude of phenomenological rate constants.

This structural determinism is critical in biological signaling cascades, where bistability functions as a switch-like threshold between divergent cellular trajectories, filtering noise from sustained environmental stimuli 1. The chemical network acts as a physical switch; crossing a specific concentration of signaling molecules forces the cell to choose one response, such as antiviral defense versus programmed cell death, over another. The emergence of these discrete operational branches necessitates a departure from local linear stability analysis, requiring a global topological assessment of the thermodynamic state space to capture the multiplicity of dissipative structures and their respective basins of attraction.

Quantifying Attractor Stability Using Quasi-Potential Landscapes

To quantify transition probabilities and the relative stabilities of these coexisting stationary states, we construct a generalized quasi-potential landscape, extending equilibrium thermodynamic potentials into the far-from-equilibrium regime. This approach maps the invariant probability density of homogeneous steady states, rendering the kinetic network as a topological surface characterized by landscape peaks and deep basins of attraction 2. Mapping this potential requires multidimensional thematic overlays of state variables to isolate regions of high thermodynamic stability from the unstable saddle points that mediate stochastic transitions 3.

The features of this landscape correspond directly to thermodynamic properties:

• Basins of attraction: The depth of each basin correlates inversely with the local entropy production rate of the corresponding non-equilibrium stationary state.
• Landscape peaks: These represent kinetic activation barriers defined by excess entropy production along the transition path.

When system parameters traverse a critical threshold, the topological deformation of this quasi-potential landscape induces a catastrophic shift in the probability distribution, annihilating one local minimum and forcing the system to relax into the surviving macroscopic attractor. We will formally quantify the stochastic fluctuations driving these transitions across saddle points when we introduce the Fokker-Planck formalism in subsequent analyses.

Hysteresis and Chemical Memory under Modulated Excitation

The coexistence of several attractors within the quasi-potential landscape gives rise to hysteresis, a phenomenon where the instantaneous state of the chemical network exhibits an irreversible dependence on its kinetic history. When highly nonlinear chemical systems undergo amplitude-modulated excitation, external periodic forcing interacts with the intrinsic relaxation timescales of the reaction network, generating dynamical behaviors including:

• Vibrational resonance
• Multistability
• Deterministic chaos 4

The modulation frequency determines the efficiency of the control force; analytical investigations utilizing Lyapunov exponents and Poincaré sections reveal that unmodulated cases, or those where the frequency ratio avoids integer proportions, significantly alter the structural stability of coexisting attractors 4. Consequently, traversing the hysteresis loop via the cyclic variation of a primary bifurcation parameter does not simply reverse the forward kinetic trajectory. Instead, it forces the system to remain on its current thermodynamic branch until the boundary of the basin of attraction is compromised via a saddle-node bifurcation. This fundamental asymmetry in forward and reverse kinetic pathways underscores the macroscopic irreversibility of dissipative structures, providing a robust thermodynamic mechanism for chemical memory and symmetry breaking in self-organizing manifolds.