Topological charge-density-vector method of identifying filaments of scroll waves.

Scroll waves have been found in a variety of three-dimensional excitable media, including physical, chemical, and biological origins. Scroll waves in cardiac tissue are of particular significance as they underlie ventricular fibrillation that can cause sudden death. The behavior of a scroll wave is characterized by a line of phase singularity at its organizing center, known as a filament. A thorough investigation into the filament dynamics is the key to further exploration of the general theory of scroll waves in excitable media and the mechanisms of ventricular fibrillation. In this paper, we propose a method to identify filaments of scroll waves in excitable media. From the definition of the topological charge of filaments, we obtain the discrete expression of the topological charge-density vector, which is useful in calculating the topological charge vectors at each grid in the space directly. The set of starting points of these topological charge vectors represents a set of phase singularities, thereby forming a line of phase singularity, that is, a filament of a scroll wave.

Numerical study of the drift of scroll waves by optical feedback in cardiac tissue.

Nonlinear waves were found in various types of physical, chemical, and biological excitable media, e.g., in heart muscle. They can form three-dimensional (3D) vortices, called scroll waves, that are of particular significance in the heart, as they underlie lethal cardiac arrhythmias. Thus controlling the behavior of scroll waves is interesting and important. Recently, the optical feedback control procedure for two-dimensional vortices, called spiral waves, was developed. It can induce directed linear drift of spiral waves in optogenetically modified cardiac tissue. However, the extension of this methodology to 3D scroll waves is nontrivial, as optogenetic signals only penetrate close to the surface of cardiac tissue. Here we present a study of this extension in a two-variable reaction-diffusion model and in a detailed model of cardiac tissue. We show that the success of the control procedure is determined by the tension of the scroll wave filament. In tissue with positive filament tension the control procedure works in all cases. However, in the case of negative filament tension for a sufficiently large medium, instabilities occur and make drift and control of scroll waves impossible. Because in normal cardiac tissue the filament tension is assumed to be positive, we conclude that the proposed optical feedback scheme can be a robust method in inducing the linear drift of scroll waves that can control their positions in cardiac tissue.

Control of spiral waves in optogenetically modified cardiac tissue by periodic optical stimulation.

Resonant drift of nonlinear spiral waves occurs in various types of excitable media under periodic stimulation. Recently a novel methodology of optogenetics has emerged, which allows to affect excitability of cardiac cells and tissues by optical stimuli. In this paper we study if resonant drift of spiral waves in the heart can be induced by a homogeneous weak periodic optical stimulation of cardiac tissue. We use a two-variable and a detailed model of cardiac tissue and add description of depolarizing and hyperpolarizing optogenetic ionic currents. We show that weak periodic optical stimulation at a frequency equal to the natural rotation frequency of the spiral wave induces resonant drift for both depolarizing and hyperpolarizing optogenetic currents. We quantify these effects and study how the speed of the drift and its direction depend on the initial conditions. We also derive analytical formulas based on the response function theory which correctly predict the drift velocity and its trajectory. We conclude that optogenetic methodology can be used for control of spiral waves in cardiac tissue and discuss its possible applications.

Spiral wave drift under optical feedback in cardiac tissue.

Spiral waves occur in various types of excitable media and their dynamics determine the spatial excitation patterns. An important type of spiral wave dynamics is drift, as it can control the position of a spiral wave or eliminate a spiral wave by forcing it to the boundary. In theoretical and experimental studies of the Belousov-Zhabotinsky reaction, it was shown that the most direct way to induce the controlled drift of spiral waves is by application of an external electric field. Mathematically such drift occurs due to the onset of additional gradient terms in the Laplacian operator describing excitable media. However, this approach does not work for cardiac excitable tissue, where an external electric field does not result in gradient terms. In this paper, we propose a method of how to induce a directed linear drift of spiral waves in cardiac tissue, which can be realized as an optical feedback control in tissue where photosensitive ion channels are expressed. We illustrate our method by using the FitzHugh-Nagumo model for cardiac tissue and the generic model of photosensitive ion channels. We show that our method works for continuous and discrete light sources and can effectively move spiral waves in cardiac tissue, or eliminate them by collisions with the boundary or with another spiral wave. We finally implement our method by using a biophysically motivated photosensitive ion channel model included to the Luo-Rudy model for cardiac cells and show that the proposed feedback control also induces directed linear drift of spiral waves in a wide range of light intensities.

Topological charge-density method of identifying phase singularities in cardiac fibrillation.

Spiral waves represent the key motifs of typical self-sustained dynamical patterns in excitable systems such as cardiac tissue. The motion of phase singularities (PSs) that lies at the center of spiral waves captures many qualitative and, in some cases, quantitative features of their complex dynamics. Recent clinical studies suggested that ablating the tissue at PS locations may cure atrial fibrillation. Here, we propose a different method to determine the location of PSs. Starting from the definition of the topological charge of spiral waves, we obtain the expression of the topological charge density in a discrete case. With this expression, we can calculate the topological charge at each grid in the space directly, so as to accurately identify the position of PSs.

Finding type and location of the source of cardiac arrhythmias from the averaged flow velocity field using the determinant-trace method.

Life threatening cardiac arrhythmias result from abnormal propagation of nonlinear electrical excitation waves in the heart. Finding the locations of the sources of these waves remains a challenging problem. This is mainly due to the low spatial resolution of electrode recordings of these waves. Also, these recordings are subjected to noise. In this paper, we develop a different approach: the AFV-DT method based on an averaged flow velocity (AFV) technique adopted from the analysis of optical flows and the determinant-trace (DT) method used for vector field analysis of dynamical systems. This method can find the location and determine all important types of sources found in excitable media such as focal activity, spiral waves, and waves rotating around obstacles. We test this method on in silico data of various wave excitation patterns obtained using the Luo-Rudy model for cardiac tissue. We show that the method works well for data with low spatial resolutions (up to 8×8) and is stable against noise. Finally, we apply it to two clinical cases and show that it can correctly identify the arrhythmia type and location. We discuss further steps on the development and improvement of this approach.

Removal of pinned scroll waves in cardiac tissues by electric fields in a generic model of three-dimensional excitable media.

Spirals or scroll waves pinned to heterogeneities in cardiac tissues may cause lethal arrhythmias. To unpin these life-threatening spiral waves, methods of wave emission from heterogeneities (WEH) induced by low-voltage pulsed DC electric fields (PDCEFs) and circularly polarized electric fields (CPEFs) have been used in two-dimensional (2D) cardiac tissues. Nevertheless, the unpinning of scroll waves in three-dimensional (3D) cardiac systems is much more difficult than that of spiral waves in 2D cardiac systems, and there are few reports on the removal of pinned scroll waves in 3D cardiac tissues by electric fields. In this article, we investigate in detail the removal of pinned scroll waves in a generic model of 3D excitable media using PDCEF, AC electric field (ACEF) and CPEF, respectively. We find that spherical waves can be induced from the heterogeneities by these electric fields in initially quiescent excitable media. However, only CPEF can induce spherical waves with frequencies higher than that of the pinned scroll wave. Such higher-frequency spherical waves induced by CPEF can be used to drive the pinned scroll wave out of the cardiac systems. We hope this remarkable ability of CPEF can provide a better alternative to terminate arrhythmias caused by pinned scroll waves.

Wave trains induced by circularly polarized electric fields in cardiac tissues.

Clinically, cardiac fibrillation caused by spiral and turbulent waves can be terminated by globally resetting electric activity in cardiac tissues with a single high-voltage electric shock, but it is usually associated with severe side effects. Presently, a promising alternative uses wave emission from heterogeneities induced by a sequence of low-voltage uniform electric field pulses. Nevertheless, this method can only emit waves locally near obstacles in turbulent waves and thereby requires multiple obstacles to globally synchronize myocardium and thus to terminate fibrillation. Here we propose a new approach using wave emission from heterogeneities induced by a low-voltage circularly polarized electric field (i.e., a rotating uniform electric field). We find that, this approach can generate circular wave trains near obstacles and they propagate outwardly. We study the characteristics of such circular wave trains and further find that, the higher-frequency circular wave trains can effectively suppress spiral turbulence.

Core solutions of rigidly rotating spiral waves in highly excitable media.

Analytical spiral wave solutions for reaction-diffusion equations play an important role in studying spiral wave dynamics. In this paper, we focus on such analytical solutions in the case of highly excitable media. We present numerical evidence that, for rigidly rotating spiral waves in highly excitable media, the species values in the spiral core region do harmonic oscillations but not relaxation ones, and their amplitudes grow linearly with the distance from the rotation center. An analytical solution is proposed to describe such spiral wave dynamics, and the quantitative comparisons between the numerical results and the analytical solutions show that the proposed spiral core solution works well in highly excitable media.

Chiralities of spiral waves and their transitions.

The chiralities of spiral waves usually refer to their rotation directions (the turning orientations of the spiral temporal movements as time elapses) and their curl directions (the winding orientations of the spiral spatial geometrical structures themselves). Traditionally, they are the same as each other. Namely, they are both clockwise or both counterclockwise. Moreover, the chiralities are determined by the topological charges of spiral waves, and thus they are conserved quantities. After the inwardly propagating spirals were experimentally observed, the relationship between the chiralities and the one between the chiralities and the topological charges are no longer preserved. The chiralities thus become more complex than ever before. As a result, there is now a desire to further study them. In this paper, the chiralities and their transition properties for all kinds of spiral waves are systemically studied in the framework of the complex Ginzburg-Landau equation, and the general relationships both between the chiralities and between the chiralities and the topological charges are obtained. The investigation of some other models, such as the FitzHugh-Nagumo model, the nonuniform Oregonator model, the modified standard model, etc., is also discussed for comparison.

Electric-field-sustained spiral waves in subexcitable media.

We present numerical evidence that, in the presence of a suitable electric field, an isolated broken plane wave retracting originally in subexcitable media can propagate continuously and eventually evolve into a rotating spiral. Simulations for the FitzHugh-Nagumo, the Barkley, and the Oregonator models are carried out and the same electric-field-sustained spiral phenomena are observed. Semianalytical results in the framework of a kinematic theory are quantitatively consistent with the numerical results.

Influence of impurities on solitons in the nonlinear LC transmission line.

This paper studies the propagation of solitons in the nonlinear LC transmission line (NLCTL) with capacitor impurity. Based on Kirchhoff's laws, the numerical simulation shows that the amplitude of the soliton will be increased or decreased when it is close to the positive or negative impurity. Then, it will be split into reflected and transmitted waves by both the positive and negative impurities. Furthermore, their final amplitude and propagating speeds are almost independent of the impurity polarity. The observations near the impurity can be understood in the physical picture of the linear uncoupled energy absorption. By these results, we find that the impurity-soliton interactions (ISIs) in NLCTLs for both inductance and capacitance impurities, which have been seen to be different before, actually can be unified. They also indicate that the ISIs in NLCTLs essentially can be integrated with those in many other soliton systems, such as the Frenkel-Kontorova model and hydrodynamics. Moreover, the impurity-induced influence on the NLCTL solitons can also be well interpreted in the framework of the nonlinear Schrödinger model with an impurity term derived from the discrete voltage propagation equations by means of the perturbation method.