Varying the direction of propagation in reaction-diffusion equations in periodic media

We consider a multidimensional reaction-diffusion equation of either ignition or monostable type, involving periodic heterogeneity, and analyze the dependence of the propagation phenomena on the direction. We prove that the (minimal) speed of the underlying pulsating fronts depends continuously on the direction of propagation, and so does its associated profile provided it is unique up to time shifts. We also prove that the spreading properties \cite{Wein02} are actually uniform with respect to the direction

Nonlinear physics of electrical wave propagation in the heart: a review

The beating of the heart is a synchronized contraction of muscle cells (myocytes) that are triggered by a periodic sequence of electrical waves (action potentials) originating in the sino-atrial node and propagating over the atria and the ventricles. Cardiac arrhythmias like atrial and ventricular fibrillation (AF,VF) or ventricular tachycardia (VT) are caused by disruptions and instabilities of these electrical excitations, that lead to the emergence of rotating waves (VT) and turbulent wave patterns (AF,VF). Numerous simulation and experimental studies during the last 20 years have addressed these topics. In this review we focus on the nonlinear dynamics of wave propagation in the heart with an emphasis on the theory of pulses, spirals and scroll waves and their instabilities in excitable media and their application to cardiac modeling. After an introduction into electrophysiological models for action potential propagation, the modeling and analysis of spatiotemporal alternans, spiral and scroll meandering, spiral breakup and scroll wave instabilities like negative line tension and sproing are reviewed in depth and discussed with emphasis on their impact in cardiac arrhythmias.Peer ReviewedPreprin

Cellular pattern formation during Dictyostelium aggregation

The development of multicellularity in the life cycle of Dictyostelium discoideum provides a paradigm model system for biological pattern formation. Previously, mathematical models have shown how a collective pattern of cell communication by waves of the messenger molecule cyclic adenosine 3′5′-monophosphate (cAMP) arises from excitable local cAMP kinetics and cAMP diffusion. Here we derive a model of the actual cell aggregation process by considering the chemotactic cell response to cAMP and its interplay with the cAMP dynamics. Cell density, which previously has been treated as a spatially homogeneous parameter, is a crucial variable of the aggregation model. We find that the coupled dynamics of cell chemotaxis and cAMP reaction-diffusion lead to the break-up of the initially uniform cell layer and to the formation of the striking cell stream morphology which characterizes the aggregation process in situ. By a combination of stability analysis and two-dimensional simulations of the model equations, we show cell streaming to be the consequence of the growth of a small-amplitude pattern in cell density forced by the large-amplitude cAMP waves, thus representing a novel scenario of spatial patterning in a cell chemotaxis system. The instability mechanism is further analysed by means of an analytic caricature of the model, and the condition for chemotaxis-driven instability is found to be very similar to the one obtained for the standard (non-oscillatory) Keller-Segel system. The growing cell stream pattern feeds back into the cAMP dynamics, which can explain in some detail experimental observations on the time evolution of the cAMP wave pattern, and suggests the characterization of the Dictyostelium aggregation field as a self-organized excitable medium

Front propagation in laminar flows

The problem of front propagation in flowing media is addressed for laminar velocity fields in two dimensions. Three representative cases are discussed: stationary cellular flow, stationary shear flow, and percolating flow. Production terms of Fisher-Kolmogorov-Petrovskii-Piskunov type and of Arrhenius type are considered under the assumption of no feedback of the concentration on the velocity. Numerical simulations of advection-reaction-diffusion equations have been performed by an algorithm based on discrete-time maps. The results show a generic enhancement of the speed of front propagation by the underlying flow. For small molecular diffusivity, the front speed $V_f$ depends on the typical flow velocity $U$ as a power law with an exponent depending on the topological properties of the flow, and on the ratio of reactive and advective time-scales. For open-streamline flows we find always $V_f \sim U$, whereas for cellular flows we observe $V_f \sim U^{1/4}$ for fast advection, and $V_f \sim U^{3/4}$ for slow advection.Comment: Enlarged, revised version, 37 pages, 14 figure

Effect of small-world topology on wave propagation on networks of excitable elements

We study excitation waves on a Newman-Watts small-world network model of coupled excitable elements. Depending on the global coupling strength, we find differing resilience to the added long-range links and different mechanisms of propagation failure. For high coupling strengths, we show agreement between the network and a reaction-diffusion model with additional mean-field term. Employing this approximation, we are able to estimate the critical density of long-range links for propagation failure.Comment: 19 pages, 8 figures and 5 pages supplementary materia