Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2019 Jun;99(6-1):062409.
doi: 10.1103/PhysRevE.99.062409.

Chaotic tip trajectories of a single spiral wave in the presence of heterogeneities

Daniel M Lombardo  1 , Wouter-Jan Rappel  1
Affiliations

Chaotic tip trajectories of a single spiral wave in the presence of heterogeneities

Daniel M Lombardo et al. Phys Rev E. 2019 Jun.

Abstract

Spiral waves have been observed in a variety of physical, chemical, and biological systems. They play a major role in cardiac arrhythmias, including fibrillation, where the observed irregular activation patterns are generally thought to arise from the continuous breakup of multiple unstable spiral waves. Using spatially extended simulations of different electrophysiological models of cardiac tissue, we show that a single spiral wave in the presence of heterogeneities can display chaotic tip trajectories, consistent with fibrillation. We also show that the simulated spiral tip dynamics, including chaotic trajectories, can be captured by a simple particle model which only describes the dynamics of the spiral tip. This shows that spiral wave breakup, or interactions with other waves, are not necessary to initiate chaos in spiral waves.

PubMed Disclaimer

Figures

FIG. 1.
FIG. 1.
A, C: Snapshot of a counter-clockwise rotating spiral wave for two different parameter sets of the FK model (voltage color-coded ranging from high (red) to low (blue) values, tip trajectory shown in white, A: Set I, C: Set II). B, D: Phase diagrams in the spacing-size space for the FK model. In this, and all other figures, blue area indicates regular trajectories that wrap around both heterogeneities, purple region represents tip trajectories that circle either one of the heterogeneities, and yellow region corresponds to chaotic trajectories. Displayed spiral tip trajectories correspond to the white dots in the phase diagram. Red X’s mark the locations of the circular regions with decreased excitability. Lyapunov exponents (λ) for the chaotic trajectories are given in units of bits/second (B: Set I, D: Set II).
FIG. 2.
FIG. 2.
A: Power spectrum of the x-component of the tip trajectory for the different electrophysiological models and parameter sets used in the main text. The spectra have either one or two peaks, corresponding to the frequencies ω1 = 2π/T1 and ω2 = 2π/T2 (left panel: FK Set I, middle panel: FK Set II, right panel: KKT). B: The spiral tip trajectory for set II. The maximum and minimum value of the spatial extent of the tip trajectory are denoted by A1 and A2, respectively.
FIG. 3.
FIG. 3.
A: Snapshot of a counter-clockwise rotating spiral wave in the homogeneous KKT model. B: Sample tip trajectories of the KKT model, for different spacings of heterogeneities with radius 0.5cm. C: Phase diagrams of tip trajectories for the SP model. All scale bars are 1cm. See also Fig. 1.
FIG. 4.
FIG. 4.
A: Comparison between the full models (FK Set I right panel, FK Set II middle panel, and KKT right panel) and SP model for homogeneous media. Scale bars are 0.5cm. B, C: Cumulative error (measured in cm, and plotted using a color scheme) in the tip trajectory from the SP model for Set I (B) and II (C) as compared to the trajectory of the full model. The error was computed for every ms and for an entire rotation of the full model. D: Example of a trajectory of the SP model for the non-optimal parameters corresponding to the white symbol in C (A1 = 1.075 cm, A2 = 0.2 cm, ω1 = 0.0058 ms−1 , and ω2 = 0.0439 ms−1).
FIG. 5.
FIG. 5.
Schematic representation of a particle moving in a potential landscape with two wells, representing tissue heterogeneities.
FIG. 6.
FIG. 6.
A, B: Phase diagrams of tip trajectories for set I (A) and set II (B) in the SP model, with tip trajectories and Lyapunov exponents. See also Fig. 1. C: Phase diagram for Set I, illustrating the effect of different well depths g. The circles correspond to the parameter values of the circles in (A).
FIG. 7.
FIG. 7.
Spiral trajectories for the FK (left column) and SP models (right column) in the presence of six randomly placed heterogeneities. Regular (A, B) and chaotic (C, D) trajectories with matching placement of the heterogeneities. The radius of the heterogeneities in the FK model was 0.3cm, and the s parameter in the SP model was 0.275cm. All scale bars are 1cm.
FIG. 8.
FIG. 8.
Example trajectories for the homogeneous (left column) and heterogeneous (right column) case with added anisotropy. In the FK model, the diffusion constant in the y (vertical) direction was chosen to be 50% of the value for the homogeneous case (i.e., Dy = 0.0005 cm2/ms). In the SP model, the forcing amplitude in the x and y equation can be directly determined from the tip trajectory in the homogeneous FK model. The tip trajectories in the SP model capture the FK trajectories for both the homogeneous (left column) and heterogeneous case (right colum). Specifically, the chaotic regime is still present in both models. All scale bars are 1cm.
See this image and copyright information in PMC

Similar articles

See all similar articles

Cited by

References

    1. Lee KJ, Cox EC, and Goldstein RE, Physical Review Letters 76, 1174 (1996). - PubMed
    1. Gorelova N and Bureš J, Developmental Neurobiology 14, 353 (1983). - PubMed
    1. Jakubith S, Rotermund H, Engel W, Von Oertzen A, and Ertl G, Physical Review Letters 65, 3013 (1990). - PubMed
    1. Epstein IR and Showalter K, The Journal of Physical Chemistry 100, 13132 (1996).
    1. Winfree AT, Science 175, 634 (1972). - PubMed

LinkOut - more resources

-