Nonlinear chemical kinetic schemes derived from mechanical and electrical dynamical systems

Chemical Systems with Limit Cycles

The dynamics of a chemical reaction network (CRN) is often modeled under the assumption of mass action kinetics by a system of ordinary differential equations (ODEs) with polynomial right-hand sides that describe the time evolution of concentrations of chemical species involved. Given an arbitrarily large integer [Formula: see text], we show that there exists a CRN such that its ODE model has at least K stable limit cycles. Such a CRN can be constructed with reactions of at most second-order provided that the number of chemical species grows linearly with K. Bounds on the minimal number of chemical species and the minimal number of chemical reactions are presented for CRNs with K stable limit cycles and at most second order or seventh-order kinetics. We also show that CRNs with only two chemical species can have K stable limit cycles, when the order of chemical reactions grows linearly with K.

The onset of dissipative chaos driven by nonequilibrium conditions

Dissipative chaos appears widely in various nonequilibrium systems; however, it is not clear how dissipative chaos originates from nonequilibrium. We discuss a framework based on the potential-flux approach to study chaos from the perspective of nonequilibrium dynamics. In this framework, chaotic systems possess a wide basin on the potential landscape, in which the rotational flux dominates the system dynamics, and chaos occurs with the appearance of this basin. In contrast, the probability flux is particularly associated with the detailed balance-breaking in nonequilibrium systems. This implies that the appearance of dissipative chaos is driven by nonequilibrium conditions.

The dynamic and thermodynamic origin of dissipative chaos: chemical Lorenz system

Chaos appears widely in various chemical and physical systems and is often accompanied by nonequilibrium due to its dissipative nature. However, it is still not clear how dissipative chaos is influenced by nonequilibrium conditions. Here, we study chaos from the perspective of nonequilibrium dynamics by considering a chemical Lorenz system. We found that its nonequilibrium nature can be quantified from the steady-state probability flux in the state space. The dynamic origin for the onset and offset of dissipative chaos was from the sudden appearance and disappearance of such nonequilibrium fluxes. Meanwhile, the dissipation associated with the flux as quantified by the entropy production rate provides the thermodynamic origin of dissipative chaos. Sharp changes in the degree of nonequilibrium also provide alternative quantitative indicators for the onset and offset of dissipative chaos.

Stochastic approach to entropy production in chemical chaos

Methods are presented to evaluate the entropy production rate in stochastic reactive systems. These methods are shown to be consistent with known results from nonequilibrium chemical thermodynamics. Moreover, it is proved that the time average of the entropy production rate can be decomposed into the contributions of the cycles obtained from the stoichiometric matrix in both stochastic processes and deterministic systems. These methods are applied to a complex reaction network constructed on the basis of Rössler's reinjection principle and featuring chemical chaos.

Parametric spatiotemporal oscillation in reaction-diffusion systems

We consider a reaction-diffusion system in a homogeneous stable steady state. On perturbation by a time-dependent sinusoidal forcing of a suitable scaling parameter the system exhibits parametric spatiotemporal instability beyond a critical threshold frequency. We have formulated a general scheme to calculate the threshold condition for oscillation and the range of unstable spatial modes lying within a V-shaped region reminiscent of Arnold's tongue. Full numerical simulations show that depending on the specificity of nonlinearity of the models, the instability may result in time-periodic stationary patterns in the form of standing clusters or spatially localized breathing patterns with characteristic wavelengths. Our theoretical analysis of the parametric oscillation in reaction-diffusion system is corroborated by full numerical simulation of two well-known chemical dynamical models: chlorite-iodine-malonic acid and Briggs-Rauscher reactions.

Potential flux landscapes determine the global stability of a Lorenz chaotic attractor under intrinsic fluctuations

We developed a potential flux landscape theory to investigate the dynamics and the global stability of a chemical Lorenz chaotic strange attractor under intrinsic fluctuations. Landscape was uncovered to have a butterfly shape. For chaotic systems, both landscape and probabilistic flux are crucial to the dynamics of chaotic oscillations. Landscape attracts the system down to the chaotic attractor, while flux drives the coherent motions along the chaotic attractors. Barrier heights from the landscape topography provide a quantitative measure for the robustness of chaotic attractor. We also found that the entropy production rate and phase coherence increase as the molecular numbers increase. Power spectrum analysis of autocorrelation function provides another way to quantify the global stability of chaotic attractor. We further found that limit cycle requires more flux and energy to sustain than the chaotic strange attractor. Finally, by detailed analysis we found that the curl probabilistic flux may provide the origin of the chaotic attractor.

A Chemical Chaotic System Derived from Chua's Circuit

Chua's circuit is converted into a mass action chemical system with Samardzija's nonlinear transformation method.

Lotka-Volterra representation of general nonlinear systems

In this article we elaborate on the structure of the generalized Lotka-Volterra (GLV) form for nonlinear differential equations. We discuss here the algebraic properties of the GLV family, such as the invariance under quasimonomial transformations and the underlying structure of classes of equivalence. Each class possesses a unique representative under the classical quadratic Lotka-Volterra form. We show how other standard modeling forms of biological interest, such as S-systems or mass-action systems, are naturally embedded into the GLV form, which thus provides a formal framework for their comparison and for the establishment of transformation rules. We also focus on the issue of recasting of general nonlinear systems into the GLV format. We present a procedure for doing so and point at possible sources of ambiguity that could make the resulting Lotka-Volterra system dependent on the path followed. We then provide some general theorems that define the operational and algorithmic framework in which this is not the case.