Terminal Transient Phase of Chaotic Transients

Annihilation dynamics during spiral defect chaos revealed by particle models

Pair-annihilation events are ubiquitous in a variety of spatially extended systems and are often studied using computationally expensive simulations. Here, we develop an approach in which we simulate the pair-annihilation of spiral wave tips in cardiac models using a computationally efficient particle model. Spiral wave tips are represented as particles with dynamics governed by diffusive behavior and short-ranged attraction. The parameters for diffusion and attraction are obtained by comparing particle motion to the trajectories of spiral wave tips in cardiac models during spiral defect chaos. The particle model reproduces the annihilation rates of the cardiac models and can determine the statistics of spiral wave dynamics, including its mean termination time. We show that increasing the attraction coefficient sharply decreases the mean termination time, making it a possible target for pharmaceutical intervention.

Dissolution of spiral wave's core using cardiac optogenetics

Rotating spiral waves in the heart are associated with life-threatening cardiac arrhythmias such as ventricular tachycardia and fibrillation. These arrhythmias are treated by a process called defibrillation, which forces electrical resynchronization of the heart tissue by delivering a single global high-voltage shock directly to the heart. This method leads to immediate termination of spiral waves. However, this may not be the only mechanism underlying successful defibrillation, as certain scenarios have also been reported, where the arrhythmia terminated slowly, over a finite period of time. Here, we investigate the slow termination dynamics of an arrhythmia in optogenetically modified murine cardiac tissue both in silico and ex vivo during global illumination at low light intensities. Optical imaging of an intact mouse heart during a ventricular arrhythmia shows slow termination of the arrhythmia, which is due to action potential prolongation observed during the last rotation of the wave. Our numerical studies show that when the core of a spiral is illuminated, it begins to expand, pushing the spiral arm towards the inexcitable boundary of the domain, leading to termination of the spiral wave. We believe that these fundamental findings lead to a better understanding of arrhythmia dynamics during slow termination, which in turn has implications for the improvement and development of new cardiac defibrillation techniques.

Optimising low-energy defibrillation in 2D cardiac tissue with a genetic algorithm

Sequences of low-energy electrical pulses can effectively terminate ventricular fibrillation (VF) and avoid the side effects of conventional high-energy electrical defibrillation shocks, including tissue damage, traumatic pain, and worsening of prognosis. However, the systematic optimisation of sequences of low-energy pulses remains a major challenge. Using 2D simulations of homogeneous cardiac tissue and a genetic algorithm, we demonstrate the optimisation of sequences with non-uniform pulse energies and time intervals between consecutive pulses for efficient VF termination. We further identify model-dependent reductions of total pacing energy ranging from ∼4% to ∼80% compared to reference adaptive-deceleration pacing (ADP) protocols of equal success rate (100%).

Multimodal distribution of transient time of predator extinction in a three-species food chain

The transient dynamics capture the time history in the behavior of a system before reaching an attractor. This paper deals with the statistics of transient dynamics in a classic tri-trophic food chain with bistability. The species of the food chain model either coexist or undergo a partial extinction with predator death after a transient time depending upon the initial population density. The distribution of transient time to predator extinction shows interesting patterns of inhomogeneity and anisotropy in the basin of the predator-free state. More precisely, the distribution shows a multimodal character when the initial points are located near a basin boundary and a unimodal character when chosen from a location far away from the boundary. The distribution is also anisotropic because the number of modes depends on the direction of the local of initial points. We define two new metrics, viz., homogeneity index and local isotropic index, to characterize the distinctive features of the distribution. We explain the origin of such multimodal distributions and try to present their ecological implications.

The physics of heart rhythm disorders

The global burden caused by cardiovascular disease is substantial, with heart disease representing the most common cause of death around the world. There remains a need to develop better mechanistic models of cardiac function in order to combat this health concern. Heart rhythm disorders, or arrhythmias, are one particular type of disease which has been amenable to quantitative investigation. Here we review the application of quantitative methodologies to explore dynamical questions pertaining to arrhythmias. We begin by describing single-cell models of cardiac myocytes, from which two and three dimensional models can be constructed. Special focus is placed on results relating to pattern formation across these spatially-distributed systems, especially the formation of spiral waves of activation. Next, we discuss mechanisms which can lead to the initiation of arrhythmias, focusing on the dynamical state of spatially discordant alternans, and outline proposed mechanisms perpetuating arrhythmias such as fibrillation. We then review experimental and clinical results related to the spatio-temporal mapping of heart rhythm disorders. Finally, we describe treatment options for heart rhythm disorders and demonstrate how statistical physics tools can provide insights into the dynamics of heart rhythm disorders.

Controlled generation of self-sustained oscillations in complex artificial neural networks

Spatially distinct, self-sustained oscillations in artificial neural networks are fundamental to information encoding, storage, and processing in these systems. Here, we develop a method to induce a large variety of self-sustained oscillatory patterns in artificial neural networks and a controlling strategy to switch between different patterns. The basic principle is that, given a complex network, one can find a set of nodes-the minimum feedback vertex set (mFVS), whose removal or inhibition will result in a tree-like network without any loop structure. Reintroducing a few or even a single mFVS node into the tree-like artificial neural network can recover one or a few of the loops and lead to self-sustained oscillation patterns based on these loops. Reactivating various mFVS nodes or their combinations can then generate a large number of distinct neuronal firing patterns with a broad distribution of the oscillation period. When the system is near a critical state, chaos can arise, providing a natural platform for pattern switching with remarkable flexibility. With mFVS guided control, complex networks of artificial neurons can thus be exploited as potential prototypes for local, analog type of processing paradigms.

Mitigating long transient time in deterministic systems by resetting

How long does a trajectory take to reach a stable equilibrium point in the basin of attraction of a dynamical system? This is a question of quite general interest and has stimulated a lot of activities in dynamical and stochastic systems where the metric of this estimation is often known as the transient or first passage time. In nonlinear systems, one often experiences long transients due to their underlying dynamics. We apply resetting or restart, an emerging concept in statistical physics and stochastic process, to mitigate the detrimental effects of prolonged transients in deterministic dynamical systems. We show that resetting the intrinsic dynamics intermittently to a spatial control line that passes through the equilibrium point can dramatically expedite its completion, resulting in a huge reduction in mean transient time and fluctuations around it. Moreover, our study reveals the emergence of an optimal restart time that globally minimizes the mean transient time. We corroborate the results with detailed numerical studies on two canonical setups in deterministic dynamical systems, namely, the Stuart-Landau oscillator and the Lorenz system. The key features-expedition of transient time-are found to be very generic under different resetting strategies. Our analysis opens up a door to control the mean and fluctuations in transient time by unifying the original dynamics with an external stochastic or periodic timer and poses open questions on the optimal way to harness transients in dynamical systems.

Susceptibility of transient chimera states

Chaotic dynamics of a dynamical system is not necessarily persistent. If there is (without any active intervention from outside) a transition towards a (possibly nonchaotic) attractor, this phenomenon is called transient chaos, which can be observed in a variety of systems, e.g., in chemical reactions, population dynamics, neuronal activity, or cardiac dynamics. Also, chimera states, which show coherent and incoherent dynamics in spatially distinct regions of the system, are often chaotic transients. In many practical cases, the control of the chaotic dynamics (either the termination or the preservation of the chaotic dynamics) is desired. Although the self-termination typically occurs quite abruptly and can so far in general not be properly predicted, previous studies showed that in many systems a 'terminal transient phase" (TTP) prior to the self-termination existed, where the system was less susceptible against small but finite perturbations in different directions in state space. In this study, we show that, in the specific case of chimera states, these susceptible directions can be related to the structure of the chimera, which we divide into the coherent part, the incoherent part and the boundary in between. That means, in practice, if self-termination is close we can identify the direction of perturbation which is likely to maintain the chaotic dynamics (the chimera state). This finding improves the general understanding of the state space structure during the TTP, and could contribute also to practical applications like future control strategies of epileptic seizures which have been recently related to the collapse of chimera states.

Terminating transient chaos in spatially extended systems

In many real-life systems, transient chaotic dynamics plays a major role. For instance, the chaotic spiral or scroll wave dynamics of electrical excitation waves during life-threatening cardiac arrhythmias can terminate by itself. Epileptic seizures have recently been related to the collapse of transient chimera states. Controlling chaotic transients, either by maintaining the chaotic dynamics or by terminating it as quickly as possible, is often desired and sometimes even vital (as in the case of cardiac arrhythmias). We discuss in this study that the difference of the underlying structures in state space between a chaotic attractor (persistent chaos) and a chaotic saddle (transient chaos) may have significant implications for efficient control strategies in real life systems. In particular, we demonstrate that in the latter case, chaotic dynamics in spatially extended systems can be terminated via a relatively low number of (spatially and temporally) localized perturbations. We demonstrate as a proof of principle that control and targeting of high-dimensional systems exhibiting transient chaos can be achieved with exceptionally small interactions with the system. This insight may impact future control strategies in real-life systems like cardiac arrhythmias.

State-dependent vulnerability of synchronization

A state-dependent vulnerability of synchronization is shown to exist in a complex network composed of numerically simulated electronic circuits. We demonstrate that disturbances to the local dynamics of network units can produce different outcomes to synchronization depending on the current state of its trajectory. We address such state dependence by systematically perturbing the synchronized system at states equally distributed along its trajectory. We find the states at which the perturbation desynchronizes the network to be complicatedly mixed with the ones that restore synchronization. Additionally, we characterize perturbation sets obtained for consecutive states by defining a safety index between them. Finally, we demonstrate that the observed vulnerability is due to the existence of an unstable chaotic set in the system's state space.

Spontaneous termination of chaotic spiral wave dynamics in human cardiac ion channel models

Chaotic spiral or scroll wave dynamics can be found in diverse systems. In cardiac dynamics, spiral or scroll waves of electrical excitation determine the dynamics during life-threatening arrhythmias like ventricular fibrillation. In numerical studies it was found that chaotic episodes of spiral and scroll waves can be transient, thus they terminate spontaneously. We show in this study that this behavior can also be observed using models which describe the ion channel dynamics of human cardiomyocytes (Bueno-Orovio-Cherry-Fenton model and the Ten Tusscher-Noble-Noble-Panfilov model). For both models we find that the average lifetime of the chaotic transients grows exponentially with the system size. With this behavior, we classify the systems into the group of type-II supertransients. We observe a significant difference of the breakup behavior between the models, which results in a distinct dynamics during the final phase just before the termination. The observation of a (temporally) stable single-spiral state affects the prevailing description of the dynamics of type-II supertransients as being “quasi-stationary” and also the feasibility of predicting the spontaneous termination of the spiral wave dynamics. In the long term, the relation between the breakup behavior of spiral waves and properties of chaotic transients like predictability or average transient lifetime may contribute to an improved understanding and classification of cardiac arrhythmias.

Scaling behavior of the terminal transient phase

Transient chaos can emerge in a variety of diverse systems, e.g., in chemical reactions, population dynamics, neuronal activity, or cardiac dynamics. The end of the chaotic episode can either be desired or not, depending on the specific system and application. In both cases, however, a prediction of the end of the chaotic dynamics is required. Despite the general challenges of reliably predicting chaotic dynamics for a long time period, the recent observation of a "terminal transient phase" of chaotic transients provides new insights into the transition from chaos to the subsequent (nonchaotic) regime. In spatially extended systems and also low-dimensional maps it was shown that the structure of the state space changes already a significant amount of time before the actual end of the chaotic dynamics. In this way, the terminal transient phase provides the conceptual foundation for a possible prediction of the upcoming end of the chaotic episode a significant amount of time in advance. In this study, we strengthen the general validity of the terminal transient phase by verifying its existence in another spatially extended model (Gray-Scott model) and the Hénon map, where in the latter case the underlying mechanisms can be understood in an intuitive way. Furthermore, we show that the temporal length of the terminal transient phase remains approximately constant, when changing the system size (Gray-Scott) or parameters (Hénon map) of the investigated models, although the average lifetime of the observed chaotic transients sensitively depends on these variations. Since the timescale of the terminal transient phase is in this sense relatively robust, this insight might be essential for possible applications, where the ratio between the length of the terminal transient phase and the relevant timescale of the dynamics may probably be crucial when a reasonable prediction (thus a sufficient time before) the end of the chaotic episode is required.