Mixed-mode oscillations in complex-plasma instabilities

Bifurcations and mixed mode oscillations in a bi-stable plasma model with slow parametric excitation

In this study, considering a bi-stable plasma model with slow parametric excitation, the bifurcation of periodic and chaotic responses as well as the resulting fast-slow motions is discussed analytically and numerically. For a nonautonomous fast sub-system, the generalized harmonic balancing method is utilized to obtain an averaged system. Bifurcation analysis about the averaged system shows that the critical manifolds form a S-shape structure. Meanwhile, supercritical and subcritical period doubling (PD) occurs on the upper branch simultaneously. As the frequency of the external excitation changes, bifurcation points on the limit cycle manifolds can present different relative locations. Moreover, an additional bi-stable structure induced by Cusp bifurcation emanates from the upper branch. On the other hand, the existence of a chaotic attractor and the corresponding boundary crisis phenomenon are verified using the Melnikov method and the basin of attraction. The structures of the numerical bifurcation diagram show good agreements with the analytical results. Considering two cases of low-frequency excitation, the corresponding fast-slow dynamics are discussed. It is found that, when the fast-slow flow passing the subcritical PD point, a low frequency with different magnitudes will lead to two patterns of bifurcation delay, i.e., the typical one and the excessive delay, which suppress the PD. As for the boundary crisis point, the slow passage effects show no distinct influence. Thus, three transition mechanisms based on two cases of the bifurcation structure are explained, including "fold of cycle-fold of cycle" type, "fold of cycle-delayed subcritical PD" type, and "fold of cycle-boundary crisis" type.

Mixed-mode oscillations and extreme events in fractional-order Bonhoeffer-van der Pol oscillator

In the present study, we investigate the dynamic behavior of the fractional-order Bonhoeffer-van der Pol (BVP) oscillator. Previous studies on the integer-order BVP have shown that it exhibits mixed-mode oscillations (MMOs) with respect to the frequency of external forcing. We explore the effect of fractional-order on these MMOs and observe interesting phenomena. For fractional-order q1, we find that as we vary the frequency of external forcing, the system exhibits increasingly small amplitude oscillations. Eventually, as q1 decreases, the MMOs disappear entirely, indicating that lower fractional orders eliminate the presence of MMOs in the BVP oscillator. On the other hand, for the fractional-order q2, we observe more complex MMOs compared to q1. However, we find that the elimination of MMOs occurs with less variation from the integer order 1. Intriguingly, as we change q2, the fractional-order BVP oscillator undergoes a phenomenon known as a crisis, where the attractor expands and extreme events occur. Overall, our study highlights the rich dynamics of the fractional-order BVP oscillator and its ability to display various modes of oscillations and crises as the order is changed.

Heartbeat instability as auto-oscillation between dim and bright void regimes

We investigated the self-excited as well as optogalvanically stimulated heartbeat instability in RF discharge complex plasma. Three video cameras measured the motion of the microparticles, the plasma emission, and the laser-induced fluorescence simultaneously. Comprehensive studies of the optogalvanic control of the heartbeat instability revealed that the microparticle suspension can be stabilized by a continuous laser, whereas a modulated laser beam induces the void contraction either transiently or resonantly. The resonance occurred when the laser modulation frequency coincided with the frequency of small breathing oscillations of the microparticle suspension, which are known to be a prerequisite to the heartbeat instability. Based on the experimental results we suggest that the void contraction during the instability is caused by an abrupt void transition from the dim to the bright regime [Pikalev et al., Plasma Sources Sci. Technol. 30, 035014 (2021)PSTEEU0963-025210.1088/1361-6595/abe0a2]. In the bright regime, a time-averaged electric field at the void boundary heats the electrons causing bright plasma emission inside the void. The dim void has much lower electric field at the boundary and exhibits therefore no emission feature associated with it.

Why pacing frequency affects the production of early afterdepolarizations in cardiomyocytes: An explanation revealed by slow-fast analysis of a minimal model

Early afterdepolarizations (EADs) are pathological voltage oscillations in cardiomyocytes that have been observed in response to a number of pharmacological agents and disease conditions. Phase-2 EADs consist of small voltage fluctuations during the plateau of an action potential, typically under conditions in which the action potential is elongated. Although a single-cell behavior, EADs can lead to tissue-level arrhythmias. Much is currently known about the biophysical mechanisms (i.e., the roles of ion channels and intracellular Ca^{2+} stores) for the various forms of EADs, due partially to the development and analysis of mathematical models. This includes the application of slow-fast analysis, which takes advantage of timescale separation inherent in the system to simplify its analysis. We take this further, using a minimal three-dimensional model to demonstrate that phase-2 EADs are canards formed in the neighborhood of a folded node singularity. This allows us to predict the number of EADs that can be produced for a given parameter set, and provides guidance on parameter changes that facilitate or inhibit EAD production. With this approach, we demonstrate why periodic stimulation, as occurs in intact heart, preferentially facilitates EAD production when applied at low frequencies. We also explain the origin of complex alternan dynamics that can occur with intermediate-frequency stimulation, in which varying numbers of EADs are produced with each pulse. These revelations fall out naturally from an understanding of folded node singularities, but are difficult to glean from knowledge of the biophysical mechanism for EADs alone. Therefore, understanding the canard mechanism is a useful complement to understanding of the biophysical mechanism that has been developed over years of experimental and computational investigations.

Bursting and mixed mode oscillations during the transition to limit cycle oscillations in a matrix burner

We investigate the route to self-excited thermoacoustic instability in a laminar flow multiple flame matrix burner. With an increase in the equivalence ratio, the thermoacoustic system that is initially quiet (stable operation) transitions to limit cycle oscillations through two distinct dynamical states, namely, bursting oscillations and mixed mode oscillations. The acoustic pressure oscillations transition from quiescence to large amplitudes during bursting oscillations. Such high amplitude bursting oscillations that occur well ahead of the onset of limit cycle oscillations can potentially cause structural damage. The thermoacoustic system exhibits hysteresis. The transition to limit cycle oscillations is replicated in a phenomenological model containing slow-fast time scales.

Multi-timescale systems and fast-slow analysis

Mathematical models of biological systems often have components that vary on different timescales. This multi-timescale character can lead to problems when doing computer simulations, which can require a great deal of computer time so that the components that change on the fastest time scale can be resolved. Mathematical analysis of these multi-timescale systems can be greatly simplified by partitioning them into subsystems that evolve on different time scales. The subsystems are then analyzed semi-independently, using a technique called fast-slow analysis. In this review we describe the fast-slow analysis technique and apply it to relaxation oscillations, neuronal bursting oscillations, canard oscillations, and mixed-mode oscillations. Although these examples all involve neural systems, the technique can and has been applied to other biological, chemical, and physical systems. It is a powerful analysis method that will become even more useful in the future as new experimental techniques push forward the complexity of biological models.

Multiparameter bifurcations and mixed-mode oscillations in Q-switched CO2 lasers

We study the origin of mixed-mode oscillations and related bifurcations in a fully molecular laser model that describes CO2 monomode lasers with a slow saturable absorber. Our study indicates that the presence of isolas of periodic mixed-mode oscillations, as the pump parameter and the cavity-frequency detuning change, is inherent to Q-switched CO2 monomode lasers. We compare this model, known as the dual four-level model, to the more conventional 3:2 model and to a CO2 laser model for fast saturable absorbers. In these models, we find similarities as well as qualitative differences, such as the different nature of the homoclinic tangency to a relevant unstable periodic orbit, where the Gavrilov-Shilnikov theory and its extensions may hold. We also show that there are isolas of periodic mixed-mode oscillations in a model for CO2 lasers with modulated losses, as the pump parameter varies. The coarse-grained bifurcation diagrams of the periodic mixed-mode oscillations in these models suggest that these oscillations belong to similar classes.

Merging and splitting of plasma spheroids in a dusty plasma

Dust particle growth in a plasma is a strongly disturbing phenomenon for the plasma equilibrium. It can induce many different types of low-frequency instabilities that can be experimentally observed, especially using high-speed imaging. A spectacular case has been observed in a krypton plasma where a huge density of dust particles is grown by material sputtering. The instability consists of well-defined regions of enhanced optical emission that emerge from the electrode vicinity and propagate towards the discharge center. These plasma spheroids have complex motions resulting from their mutual interaction that can also lead to the merging of two plasma spheroids into a single one. The reverse situation is also observed with the splitting of a plasma spheroid into two parts. These results are presented for the first time and reveal new behaviors in dusty plasmas.

Mixed-mode oscillations via canard explosions in light-emitting diodes with optoelectronic feedback

Chaotically spiking attractors in semiconductor lasers with optoelectronic feedback have been recently observed to be the result of canard phenomena in three-dimensional phase space (incomplete homoclinic scenarios). Since light-emitting diodes display the same dynamics and are much more easily controllable, we use one of these systems to complete the attractor analysis demonstrating experimentally and theoretically the occurrence of complex sequences of periodic mixed-mode oscillations. In particular, we investigate the transition between periodic and chaotic mixed-mode states and analyze the effects of the unavoidable experimental noise on these transitions.

Threshold phenomena in a throbbing complex plasma

In complex plasmas, the trapped dust particle cloud is often characterized by a central dust-free region ("void"). The void induces a spatial inhomogeneity of the dust particle distribution and is at the origin of many intricate unstable phenomena. One type of this kind of behavior is the so-called heartbeat instability consisting of successive contractions and expansions of the void. This instability is characterized by a strong nonlinear dynamics which can reveal the occurrence of incomplete sequences corresponding to failed contractions. Experimental results based on high-speed imaging are presented for the first time and underline this threshold effect in both the dust cloud motion and the evolution of the plasma light emission.