Extreme multistability in a chemical model system

Introduction to Focus Issue: Chimera states: From theory and experiments to technology and living systems

One of the pillars of modern science is the concept of symmetries. Spontaneously breaking such symmetries gives rise to non-trivial states, which can explain a variety of phenomena around us. Chimera states, characterized by the coexistence of localized synchronized and unsynchronized dynamics, are a novel example. This Focus Issue covers recent developments in the study of chimera states, from both theoretical and experimental points of view, including an emphasis on prospective practical realization for application in technology and living systems.

Exactly solvable Stuart-Landau models in arbitrary dimensions

We use Clifford's geometric algebra to extend the Stuart-Landau system to dimensions D>2 and give an exact solution of the oscillator equations in the general case. At the supercritical Hopf bifurcation marked by a transition from stable fixed-point dynamics to oscillatory motion, the Jacobian matrix evaluated at the fixed point has N=⌊D/2⌋ pairs of complex conjugate eigenvalues which cross the imaginary axis simultaneously. For odd D there is an additional purely real eigenvalue that does the same. Oscillatory dynamics is asymptotically confined to a hypersphere S^{D-1} and is characterised by extreme multistability, namely the coexistence of an infinite number of limiting orbits, each of which has the geometry of a torus T^{N} on which the motion is either periodic or quasiperiodic. We also comment on similar Clifford extensions of other limit cycle oscillator systems and their generalizations.

Multistable ghost attractors in a switching laser system

This paper studies the effects of a switching parameter on the dynamics of a multistable laser model. The laser model represents multistability in distinct ranges of parameters. We assume that the system's parameter switches periodically between different values. Since the system is multistable, the presence of a ghost attractor is also dependent on the initial condition. It is shown that when the composing subsystems are chaotic, a periodic ghost attractor can emerge and vice versa, depending on the initial conditions. In contrast to the previous studies in which the attractor of the fast blinking systems approximates the average attractor, here, the blinking attractor differs from the average in some cases. It is shown that when the switching parameter values are distant from their average, the blinking and the average attractors are different, and as they approach, the blinking attractor approaches the average attractor too.

Extreme Multistability and Its Incremental Integral Reconstruction in a Non-Autonomous Memcapacitive Oscillator

Extreme multistability has frequently been reported in autonomous circuits involving memory-circuit elements, since these circuits possess line/plane equilibrium sets. However, this special phenomenon has rarely been discovered in non-autonomous circuits. Luckily, extreme multistability is found in a simple non-autonomous memcapacitive oscillator in this paper. The oscillator only contains a memcapacitor, a linear resistor, a linear inductor, and a sinusoidal voltage source, which are connected in series. The memcapacitive system model is firstly built for further study. The equilibrium points of the memcapacitive system evolve between a no equilibrium point and a line equilibrium set with the change in time. This gives rise to the emergence of extreme multistability, but the forming mechanism is not clear. Thus, the incremental integral method is employed to reconstruct the memcapacitive system. In the newly reconstructed system, the number and stability of the equilibrium points have complex time-varying characteristics due to the presence of fold bifurcation. Furthermore, the forming mechanism of the extreme multistability is further explained. Note that the initial conditions of the original memcapacitive system are mapped onto the controlling parameters of the newly reconstructed system. This makes it possible to achieve precise control of the extreme multistability. Furthermore, an analog circuit is designed for the reconstructed system, and then PSIM circuit simulations are performed to verify the numerical results.

Oscillatory microcells connected on a ring by chemical waves

The dynamics of four coupled microcells with the oscillatory Belousov-Zhabotinsky (BZ) reaction in them is analyzed with the aid of partial differential equations. Identical BZ microcells are coupled in a circle via identical narrow channels containing all the components of the BZ reaction, which is in the stationary excitable state in the channels. Spikes in the BZ microcells generate unidirectional chemical waves in the channels. A thin filter is put in between the end of the channel and the cell. To make coupling between neighboring cells of the inhibitory type, hydrophobic filters are used, which let only Br₂ molecules, the inhibitor of the BZ reaction, go through the filter. To simulate excitatory coupling, we use a hypothetical filter that let only HBrO₂ molecules, the activator of the BZ reaction, go through it. New dynamic modes found in the described system are compared with the "old" dynamic modes found earlier in the analogous system of the "single point" BZ oscillators coupled in a circle by pulses with time delay. The "new" and "old" dynamic modes found for inhibitory coupling match well, the only difference being much broader regions of multi-rhythmicity in the "new" dynamic modes. For the excitatory type of coupling, in addition to four symmetrical modes of the "old" type, many new asymmetrical modes coexisting with the symmetrical ones have been found. Asymmetrical modes are characterized by the spikes occurring any time within some finite time intervals.

Synchronization within synchronization: transients and intermittency in ecological networks

Transients are fundamental to ecological systems with significant implications to management, conservation and biological control. We uncover a type of transient synchronization behavior in spatial ecological networks whose local dynamics are of the chaotic, predator–prey type. In the parameter regime where there is phase synchronization among all the patches, complete synchronization (i.e. synchronization in both phase and amplitude) can arise in certain pairs of patches as determined by the network symmetry—henceforth the phenomenon of  ‘synchronization within synchronization.’ Distinct patterns of complete synchronization coexist but, due to intrinsic instability or noise, each pattern is a transient and there is random, intermittent switching among the patterns in the course of time evolution. The probability distribution of the transient time is found to follow an algebraic scaling law with a divergent average transient lifetime. Based on symmetry considerations, we develop a stability analysis to understand these phenomena. The general principle of symmetry can also be exploited to explain previously discovered, counterintuitive synchronization behaviors in ecological networks.

Birhythmic Analog Circuit Maze: A Nonlinear Neurostimulation Testbed

Brain dynamics can exhibit narrow-band nonlinear oscillations and multistability. For a subset of disorders of consciousness and motor control, we hypothesized that some symptoms originate from the inability to spontaneously transition from one attractor to another. Using external perturbations, such as electrical pulses delivered by deep brain stimulation devices, it may be possible to induce such transition out of the pathological attractors. However, the induction of transition may be non-trivial, rendering the current open-loop stimulation strategies insufficient. In order to develop next-generation neural stimulators that can intelligently learn to induce attractor transitions, we require a platform to test the efficacy of such systems. To this end, we designed an analog circuit as a model for the multistable brain dynamics. The circuit spontaneously oscillates stably on two periods as an instantiation of a 3-dimensional continuous-time gated recurrent neural network. To discourage simple perturbation strategies, such as constant or random stimulation patterns from easily inducing transition between the stable limit cycles, we designed a state-dependent nonlinear circuit interface for external perturbation. We demonstrate the existence of nontrivial solutions to the transition problem in our circuit implementation.

Noise-enabled species recovery in the aftermath of a tipping point

The beneficial role of noise in promoting species coexistence and preventing extinction has been recognized in theoretical ecology, but previous studies were mostly concerned with low-dimensional systems. We investigate the interplay between noise and nonlinear dynamics in real-world complex mutualistic networks with a focus on species recovery in the aftermath of a tipping point. Particularly, as a critical parameter such as the mutualistic interaction strength passes through a tipping point, the system collapses and approaches an extinction state through a dramatic reduction in the species populations to near-zero values. We demonstrate the striking effect of noise: when the direction of parameter change is reversed through the tipping point, noise enables species recovery which otherwise would not be possible. We uncover an algebraic scaling law between the noise amplitude and the parameter distance from the tipping point to the recovery point and provide a physical understanding through analyzing the nonlinear dynamics based on an effective, reduced-dimension model. Noise, in the form of small population fluctuations, can thus play a positive role in protecting high-dimensional, complex ecological networks.

Multistability and basin stability in coupled pendulum clocks

In this paper, we investigate the phenomenon of multistability and the concept of basin stability in two coupled pendula with escapement mechanisms, suspended on horizontally oscillating beam. The dynamics of a single pendulum clock is studied and described, showing possible responses of the unit. The basin stability maps are discussed in two-parameters plane, where we vary both the system's stiffness as well as the damping. The possible attractors for the investigated clocks are discussed, showing that different patterns of synchronization and desynchronization can occur. The oscillators may completely synchronize in one of the three possible combinations (including inphase and antiphase ones), practically synchronize with some fluctuations or stay in the irregular pattern, which includes chaotic motion. The transitions between solutions are studied, uncovering that the road from one type of dynamics into another may become very complex. Moreover, we examine the multistability property of our model using the bifurcation diagrams and the basins of attraction maps, discussing possible scenarios in which the states co-exist. The analysis of attractors' basins uncovers complicated structure of the latter ones, exhibiting that the final behavior of investigated model may be hard to determine and trace. Our results are discussed for the cases of identical and nonidentical pendula, as well as light and heavy beam, showing that depending on considered scenario, various patterns of behaviors and transitions may be observed. The research described in this paper proves that the mechanical properties of the system's suspension may play a crucial role in the possibility of the appearance of different types of attractors and that the basin stabilities of states strictly depend on the values of considered parameters.

Extreme multistability in a Josephson-junction-based circuit

We design and report an electrical circuit using a Josephson junction under periodic forcing that reveals extreme multistability. Its overall state equations surprisingly recall those of a well-known model of Josephson junction initially introduced in our circuit. The final circuit is characterized by the presence of two new and different current sources in parallel with the nonlinear internal current source sin[ϕ(t)] of the Josephson junction single electronic component. Furthermore, the model presents an interesting extreme multistability which is justified by a very large number of different attractors (chaotic or not) when slightly changing the initial conditions.

Complex Dynamics in a Memcapacitor-Based Circuit

In this paper, a new memcapacitor model and its corresponding circuit emulator are proposed, based on which, a chaotic oscillator is designed and the system dynamic characteristics are investigated, both analytically and experimentally. Extreme multistability and coexisting attractors are observed in this complex system. The basins of attraction, multistability, bifurcations, Lyapunov exponents, and initial-condition-triggered similar bifurcation are analyzed. Finally, the memcapacitor-based chaotic oscillator is realized via circuit implementation with experimental results presented.

A Chaotic System with Infinite Equilibria and Its S-Box Constructing Application

Systems with many equilibrium points have attracted considerable interest recently. A chaotic system with a line equilibrium has been studied in this work. The system has infinite equilibria and exhibits coexisting chaotic attractors. The system with an infinite number of equilibria has been realized by an electronic circuit, which confirms the feasibility of the system. Based on such a system, we have developed a new S-Box generation algorithm. With the developed algorithm, two new S-Boxes are produced. Performance tests of S-Boxes are performed. The tests have shown that proposed S-Boxes have good performance results.

New topological tool for multistable dynamical systems

We introduce a new method for investigation of dynamical systems which allows us to extract as much information as possible about potential system dynamics, based only on the form of equations describing it. The discussed tool of critical surfaces, defined by the zero velocity (and/or) acceleration field for particular variables of the system is related to the geometry of the attractors. Particularly, the developed method provides a new and simple procedure allowing to localize hidden oscillations. Our approach is based on the dimension reduction of the searched area in the phase space and has an advantage (in terms of complexity) over standard procedures for investigating full-dimensional space. The two approaches have been compared using typical examples of oscillators with hidden states. Our topological tool allows us not only to develop alternate ways of extracting information from the equations of motion of the dynamical system, but also provides a better understanding of attractors geometry and their capturing in complex cases, especially including multistable and hidden attractors. We believe that the introduced method can be widely used in the studies of dynamical systems and their applications in science and engineering.

Emergence, evolution, and control of multistability in a hybrid topological quantum/classical system

We present a novel class of nonlinear dynamical systems-a hybrid of relativistic quantum and classical systems and demonstrate that multistability is ubiquitous. A representative setting is coupled systems of a topological insulator and an insulating ferromagnet, where the former possesses an insulating bulk with topologically protected, dissipationless, and conducting surface electronic states governed by the relativistic quantum Dirac Hamiltonian and the latter is described by the nonlinear classical evolution of its magnetization vector. The interactions between the two are essentially the spin transfer torque from the topological insulator to the ferromagnet and the local proximity induced exchange coupling in the opposite direction. The hybrid system exhibits a rich variety of nonlinear dynamical phenomena besides multistability such as bifurcations, chaos, and phase synchronization. The degree of multistability can be controlled by an external voltage. In the case of two coexisting states, the system is effectively binary, opening a door to exploitation for developing spintronic memory devices. Because of the dissipationless and spin-momentum locking nature of the surface currents of the topological insulator, little power is needed for generating a significant current, making the system appealing for potential applications in next generation of low power memory devices.

Multistability and tipping: From mathematics and physics to climate and brain-Minireview and preface to the focus issue

Multistability refers to the coexistence of different stable states in nonlinear dynamical systems. This phenomenon has been observed in laboratory experiments and in nature. In this introduction, we briefly introduce the classes of dynamical systems in which this phenomenon has been found and discuss the extension to new system classes. Furthermore, we introduce the concept of critical transitions and discuss approaches to distinguish them according to their characteristics. Finally, we present some specific applications in physics, neuroscience, biology, ecology, and climate science.

Sensitivity analysis of the noise-induced oscillatory multistability in Higgins model of glycolysis

A phenomenon of the noise-induced oscillatory multistability in glycolysis is studied. As a basic deterministic skeleton, we consider the two-dimensional Higgins model. The noise-induced generation of mixed-mode stochastic oscillations is studied in various parametric zones. Probabilistic mechanisms of the stochastic excitability of equilibria and noise-induced splitting of randomly forced cycles are analysed by the stochastic sensitivity function technique. A parametric zone of supersensitive Canard-type cycles is localized and studied in detail. It is shown that the generation of mixed-mode stochastic oscillations is accompanied by the noise-induced transitions from order to chaos.

Coexisting multiple attractors and riddled basins of a memristive system

In this paper, a new memristor-based chaotic system is designed, analyzed, and implemented. Multistability, multiple attractors, and complex riddled basins are observed from the system, which are investigated along with other dynamical behaviors such as equilibrium points and their stabilities, symmetrical bifurcation diagrams, and sustained chaotic states. With different sets of system parameters, the system can also generate various multi-scroll attractors. Finally, the system is realized by experimental circuits.

Constructing an autonomous system with infinitely many chaotic attractors

Some classical chaotic systems such as the Lorenz system and Chua system have finite numbers of chaotic attractors. This letter develops a simple, effective method for constructing lower-dimensional autonomous systems with infinitely many chaotic attractors. As an application, a Lorenz-type system and a Rössler-type system with infinitely many chaotic attractors are constructed with bifurcation analysis, and with an extension to the fractional-order setting.

Intermittent feedback induces attractor selection

We present a method for attractor selection in multistable dynamical systems. It involves a feedback term that is active only when the dynamics of the system is in a particular fraction of state space of the attractor. We implement this method first on a simplest symmetric chaotic flow and then on a bistable neuronal system. We find that adding this space-dependent feedback term to the dynamical equations of these systems will drive the dynamics to the desired attractor by annihilating the other. We further demonstrate that the attractor selection due to this feedback term can be used in construction of logic gates, which is one of the practical applications of the proposed method.

Bistability and Bifurcation in Minimal Self-Replication and Nonenzymatic Catalytic Networks

Bistability and bifurcation, found in a wide range of biochemical networks, are central to the proper function of living systems. We investigate herein recent model systems that show bistable behavior based on nonenzymatic self-replication reactions. Such models were used before to investigate catalytic growth, chemical logic operations, and additional processes of self-organization leading to complexification. By solving for their steady-state solutions by using various analytical and numerical methods, we analyze how and when these systems yield bistability and bifurcation and discover specific cases and conditions producing bistability. We demonstrate that the onset of bistability requires at least second-order catalysis and results from a mismatch between the various forward and reverse processes. Our findings may have far-reaching implications in understanding early evolutionary processes of complexification, emergence, and potentially the origin of life.

Extreme multistability in a memristor-based multi-scroll hyper-chaotic system

In this paper, a new memristor-based multi-scroll hyper-chaotic system is designed. The proposed memristor-based system possesses multiple complex dynamic behaviors compared with other chaotic systems. Various coexisting attractors and hidden coexisting attractors are observed in this system, which means extreme multistability arises. Besides, by adjusting parameters of the system, this chaotic system can perform single-scroll attractors, double-scroll attractors, and four-scroll attractors. Basic dynamic characteristics of the system are investigated, including equilibrium points and stability, bifurcation diagrams, Lyapunov exponents, and so on. In addition, the presented system is also realized by an analog circuit to confirm the correction of the numerical simulations.

Linear Augmentation for Stabilizing Stationary Solutions: Potential Pitfalls and Their Application

Linear augmentation has recently been shown to be effective in targeting desired stationary solutions, suppressing bistablity, in regulating the dynamics of drive response systems and in controlling the dynamics of hidden attractors. The simplicity of the procedure is the main highlight of this scheme but questions related to its general applicability still need to be addressed. Focusing on the issue of targeting stationary solutions, this work demonstrates instances where the scheme fails to stabilize the required solutions and leads to other complicated dynamical scenarios. Examples from conservative as well as dissipative systems are presented in this regard and important applications in dissipative predator-prey systems are discussed, which include preventative measures to avoid potentially catastrophic dynamical transitions in these systems.

Selective monostability in multi-stable systems

We propose a robust method that allows a periodic or a chaotic multi-stable system to be transformed to a monostable system at an orbit with dominant frequency of any of the coexisting attractors. Our approach implies the selection of a particular attractor by periodic external modulation with frequency close to the dominant frequency in the power spectrum of a desired orbit and simultaneous annihilation of all other coexisting states by positive feedback, both applied to one of the system parameters. The method does not require any preliminary knowledge of the system dynamics and the phase space structure. The efficiency of the method is demonstrated in both a non-autonomous multi-stable laser with coexisting periodic orbits and an autonomous Rössler-like oscillator with coexisting chaotic attractors. The experiments with an erbium-doped fibre laser provide evidence for the robustness of the proposed method in making the system monostable at an orbit with dominant frequency of any preselected attractor.

Delayed feedback induced multirhythmicity in the oscillatory electrodissolution of copper

Occurrence of bi- and trirhythmicities (coexistence of two or three stable limit cycles, respectively, with distinctly different periods) has been studied experimentally by applying delayed feedback control to the copper-phosphoric acid electrochemical system oscillating close to a Hopf bifurcation point under potentiostatic condition. The oscillating electrode potential is delayed by τ and the difference between the present and delayed values is fed back to the circuit potential with a feedback gain K. The experiments were performed by determining the period of current oscillations T as a function of (both increasing and decreasing) τ at several fixed values of K. With small delay times, the period exhibits a sinusoidal type dependence on τ. However, with relatively large delays (typically τ ≫ T) for each feedback gain K, there exists a critical delay τcrit above which birhythmicity emerges. The experiments show that for weak feedback, Kτcrit is approximately constant. At very large delays, the dynamics becomes even more complex, and trirhythmicity could be observed. Results of numerical simulations based on a general kinetic model for metal electrodissolution were consistent with the experimental observations. The experimental and numerical results are also interpreted by using a phase model; the model parameters can be obtained from experimental data measured at small delay times. Analytical solutions to the phase model quantitatively predict the parameter regions for the appearance of birhythmicity in the experiments, and explain the almost constant value of Kτcrit for weak feedback.

Extreme multistability: Attractor manipulation and robustness

The coexistence of infinitely many attractors is called extreme multistability in dynamical systems. In coupled systems, this phenomenon is closely related to partial synchrony and characterized by the emergence of a conserved quantity. We propose a general design of coupling that leads to partial synchronization, which may be a partial complete synchronization or partial antisynchronization and even a mixed state of complete synchronization and antisynchronization in two coupled systems and, thereby reveal the emergence of extreme multistability. The proposed design of coupling has wider options and allows amplification or attenuation of the amplitude of the attractors whenever it is necessary. We demonstrate that this phenomenon is robust to parameter mismatch of the coupled oscillators.

Comment on "how to obtain extreme multistability in coupled dynamical systems"

We note that extreme multistability of the type described in the referenced paper can be achieved in virtually any dynamical system by adding extraneous variables and using their initial conditions in place of the existing parameters or as additional parameters. We show several simple examples of this and show how the referenced examples are similar.

Reply to "comment on 'how to obtain extreme multistability in coupled dynamical systems' "

In this Reply we answer the two major issues raised by the Comment. First, we point out that the idea of constructing extreme multistability in simple dynamical systems is not new and has been demonstrated previously by other authors. Furthermore, we emphasize the importance of the concept of a conserved quantity and its consequences for the dynamics, which applies to all the examples in the Comment. Second, we show that the design of controllers to achieve extreme multistability in coupled systems is as general as described in Phys. Rev. E 85, 035202(R) (2012) by providing two examples which do not lead to a master-slave dynamics.

Experimental observation of extreme multistability in an electronic system of two coupled Rössler oscillators

We report the first experimental observation of extreme multistability in a controlled laboratory investigation. Extreme multistability arises when infinitely many attractors coexist for the same set of system parameters. The behavior was predicted earlier on theoretical grounds, supported by numerical studies of models of two coupled identical or nearly identical systems. We construct and couple two analog circuits based on a modified coupled Rössler system and demonstrate the occurrence of extreme multistability through a controlled switching to different attractor states purely through a change in initial conditions for a fixed set of system parameters. Numerical studies of the coupled model equations are in agreement with our experimental findings.

Three-body interactions with time delay

This work focuses on the dynamics of globally coupled phase oscillators with three-body interaction and time delay. Analytic estimates regarding the stability of the incoherent solution are presented. Expressions for the phase synchronization frequencies and their stability are also derived. These theoretical results are supplemented with appropriate numerical computations. Numerical results regarding the fluctuations observed in the synchronization order parameter are then discussed. Some comparative results for phase synchronization in two-body, three-body, and mixed-coupled systems for different coupling combinations are also presented.

How to obtain extreme multistability in coupled dynamical systems

We present a method for designing an appropriate coupling scheme for two dynamical systems in order to realize extreme multistability. We achieve the coexistence of infinitely many attractors for a given set of parameters by using the concept of partial synchronization based on Lyapunov function stability. We show that the method is very general and allows a great flexibility in choosing the coupling. Furthermore, we demonstrate its applicability in different models, such as the Rössler system and a chemical oscillator. Finally we show that extreme multistability is robust with respect to parameter mismatch and, hence, a very general phenomenon in coupled systems.

Structures of chaos in open reaction systems

By numerically simulating the Bray-Liebhafsky (BL) reaction (the hydrogen peroxide decomposition in the presence of hydrogen and iodate ions) in a continuously fed well stirred tank reactor (CSTR), we find "structured" types of chaos emerging in regular order with respect to flow rate as the control parameter. These chaotic "structures" appear between each two successive periodic states, and have forms and evolution resembling to the neighboring periodic dynamics. More precisely, in the transition from period-doubling route to chaos to the arising periodic mixture of different mixed-mode oscillations, we are able to recognize and qualitatively and quantitatively distinguish the sequence of "period-doubling" chaos and chaos consisted of mixed-mode oscillations (the "mixed-mode structured" chaos), both appearing in regular order between succeeding periodic states. Additionally, between these types of chaos, the chaos without such recognizable "structures" ("unstructured" chaos) is also distinguished. Furthermore, all transitions between two successive periodic states are realized through bifurcation of chaotic states. This scenario is a universal feature throughout the whole mixed-mode region, as well as throughout other mixed-mode regions obtained under different initial conditions.