We consider the problem of ignition of propagating waves in one-dimensional bistable or excitable systems by an instantaneous spatially extended stimulus. Earlier we proposed a method [I. Idris and V. N. Biktashev, Phys. Rev. Lett. 101, 244101 (2008)] for analytical description of the threshold conditions based on an approximation of the (center-)stable manifold of a certain critical solution. Here we generalize this method to address a wider class of excitable systems, such as multicomponent reaction-diffusion systems and systems with non-self-adjoint linearized operators, including systems with moving critical fronts and pulses. We also explore an extension of this method from a linear to a quadratic approximation of the (center-)stable manifold, resulting in some cases in a significant increase in accuracy. The applicability of the approach is demonstrated on five test problems ranging from archetypal examples such as the Zeldovich-Frank-Kamenetsky equation to near realistic examples such as the Beeler-Reuter model of cardiac excitation. While the method is analytical in nature, it is recognized that essential ingredients of the theory can be calculated explicitly only in exceptional cases, so we also describe methods suitable for calculating these ingredients numerically.

We study the problem of initiation of excitation waves in the FitzHugh-Nagumo model. Our approach follows earlier works and is based on the idea of approximating the boundary between basins of attraction of propagating waves and of the resting state as the stable manifold of a critical solution. Here, we obtain analytical expressions for the essential ingredients of the theory by singular perturbation using two small parameters, the separation of time scales of the activator and inhibitor and the threshold in the activator's kinetics. This results in a closed analytical expression for the strength-duration curve.

No abstract

We consider the strength-duration relationship in one-dimensional spatially extended excitable media. In a previous study [1] set out to separate initial (or boundary) conditions leading to propagation wave solutions from those leading to decay solutions, an analytical criterion based on an approximation of the (center-)stable manifold of a certain critical solution was presented. The theoretical prediction in the case of strength-extent curve was later on extended to cover a wider class of excitable systems including multicomponent reaction-diffusion systems, systems with nonself-adjoint linearized operators and in particular, systems with moving critical solutions (critical fronts and critical pulses) [2]. In the present work, we consider extension of the theory to the case of strength-duration curve.