Asymptotics of conduction velocity restitution in models of electrical excitation in the heart.

We extend a non-Tikhonov asymptotic embedding, proposed earlier, for calculation of conduction velocity restitution curves in ionic models of cardiac excitability. Conduction velocity restitution is the simplest non-trivial spatially extended problem in excitable media, and in the case of cardiac tissue it is an important tool for prediction of cardiac arrhythmias and fibrillation. An idealized conduction velocity restitution curve requires solving a non-linear eigenvalue problem with periodic boundary conditions, which in the cardiac case is very stiff and calls for the use of asymptotic methods. We compare asymptotics of restitution curves in four examples, two generic excitable media models, and two ionic cardiac models. The generic models include the classical FitzHugh-Nagumo model and its variation by Barkley. They are treated with standard singular perturbation techniques. The ionic models include a simplified "caricature" of Noble (J. Physiol. Lond. 160:317-352, 1962) model and Beeler and Reuter (J. Physiol. Lond. 268:177-210, 1977) model, which lead to non-Tikhonov problems where known asymptotic results do not apply. The Caricature Noble model is considered with particular care to demonstrate the well-posedness of the corresponding boundary-value problem. The developed method for calculation of conduction velocity restitution is then applied to the Beeler-Reuter model. We discuss new mathematical features appearing in cardiac ionic models and possible applications of the developed method.

Semianalytical approach to criteria for ignition of excitation waves.

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.

Asymptotic analysis and analytical solutions of a model of cardiac excitation.

We describe an asymptotic approach to gated ionic models of single-cell cardiac excitability. It has a form essentially different from the Tikhonov fast-slow form assumed in standard asymptotic reductions of excitable systems. This is of interest since the standard approaches have been previously found inadequate to describe phenomena such as the dissipation of cardiac wave fronts and the shape of action potential at repolarization. The proposed asymptotic description overcomes these deficiencies by allowing, among other non-Tikhonov features, that a dynamical variable may change its character from fast to slow within a single solution. The general asymptotic approach is best demonstrated on an example which should be both simple and generic. The classical model of Purkinje fibers (Noble in J. Physiol. 160:317-352, 1962) has the simplest functional form of all cardiac models but according to the current understanding it assigns a physiologically incorrect role to the Na current. This leads us to suggest an "Archetypal Model" with the simplicity of the Noble model but with a structure more typical to contemporary cardiac models. We demonstrate that the Archetypal Model admits a complete asymptotic solution in quadratures. To validate our asymptotic approach, we proceed to consider an exactly solvable "caricature" of the Archetypal Model and demonstrate that the asymptotic of its exact solution coincides with the solutions obtained by substituting the "caricature" right-hand sides into the asymptotic solution of the generic Archetypal Model. This is necessary, because, unlike in standard asymptotic descriptions, no general results exist which can guarantee the proximity of the non-Tikhonov asymptotic solutions to the solutions of the corresponding detailed ionic model.

Asymptotic properties of mathematical models of excitability.

We analyse small parameters in selected models of biological excitability, including Hodgkin-Huxley (Hodgkin & Huxley 1952 J. Physiol.117, 500-544) model of nerve axon, Noble (Noble 1962 J. Physiol.160, 317-352) model of heart Purkinje fibres and Courtemanche et al. (Courtemanche et al. 1998 Am. J. Physiol.275, H301-H321) model of human atrial cells. Some of the small parameters are responsible for differences in the characteristic time-scales of dynamic variables, as in the traditional singular perturbation approaches. Others appear in a way which makes the standard approaches inapplicable. We apply this analysis to study the behaviour of fronts of excitation waves in spatially extended cardiac models. Suppressing the excitability of the tissue leads to a decrease in the propagation speed, but only to a certain limit; further suppression blocks active propagation and leads to a passive diffusive spread of voltage. Such a dissipation may happen if a front propagates into a tissue recovering after a previous wave, e.g. re-entry. A dissipated front does not recover even when the excitability restores. This has no analogy in FitzHugh-Nagumo model and its variants, where fronts can stop and then start again. In two spatial dimensions, dissipation accounts for breakups and self-termination of re-entrant waves in excitable media with Courtemanche et al. kinetics.