Feedback control of canards

We present a control mechanism for tuning a fast-slow dynamical system undergoing a supercritical Hopf bifurcation to be in the canard regime, the tiny parameter window between small and large periodic behavior. Our control strategy uses continuous feedback control via a slow control variable to cause the system to drift on average toward canard orbits. We apply this to tune the FitzHugh-Nagumo model to produce maximal canard orbits. When the controller is improperly configured, periodic or chaotic mixed-mode oscillations are found. We also investigate the effects of noise on this control mechanism. Finally, we demonstrate that a sensor tuned in this way to operate near the canard regime can detect tiny changes in system parameters.     

Spiking behavior in a noise-driven system combining oscillatory and excitatory properties

We show that firing activity (spiking) can be regularized by noise in a FitzHugh-Nagumo (FHN) neuron model when operating slightly beyond the supercritical Hopf bifurcation (in the "canard" region). We also provide the conditions for imperfect phase locking between interspike intervals and low amplitude quasiharmonic oscillations. For the imperfect phase locking no need exists of an external signal as it follows from the FHN intrinsic dynamics.     

An elementary model of torus canards

We study the recently observed phenomena of torus canards. These are a higher-dimensional generalization of the classical canard orbits familiar from planar systems and arise in fast-slow systems of ordinary differential equations in which the fast subsystem contains a saddle-node bifurcation of limit cycles. Torus canards are trajectories that pass near the saddle-node and subsequently spend long times near a repelling branch of slowly varying limit cycles. In this article, we carry out a study of torus canards in an elementary third-order system that consists of a rotated planar system of van der Pol type in which the rotational symmetry is broken by including a phase-dependent term in the slow component of the vector field. In the regime of fast rotation, the torus canards behave much like their planar counterparts. In the regime of slow rotation, the phase dependence creates rich torus canard dynamics and dynamics of mixed mode type. The results of this elementary model provide insight into the torus canards observed in a higher-dimensional neuroscience model.     

Interpretation of gas oscillations in multicylinder fluid machinery manifolds by using lumped parameter descriptions

Fluid pressure oscillations are complicated in multicylinder positive displacement machinery manifolds because of the cylinder interactions due to (i) the kinematic arrangements between cylinders and (ii) the inter-connected manifold elements. In practice, manifolds possess irregular shapes and geometries. These are difficult to analyze by using a continuous parameters approach; however, the manifold components can be discretized easily and described by acoustic lumped parameters. Additionally, this approach is more suitable for developing design guidelines, which is the primary aim of the present paper. The multicylinder interaction problem is formulated in the frequency domain. Mass flow rates, with proper crank phase relationships, are considered excitation sources. Manifold components are described as acoustic elements. The equations of acoustic motion are solved for the eigenvalues and the natural modes of gas oscillation, and these are used to provide the forced acoustic response. The formulations are then applied to a two-cylinder compressor discharge system. The computed results compare well with measurements.     

Emergence of chaotic attractor and anti-synchronization for two coupled monostable neurons

The dynamics of two coupled piece-wise linear one-dimensional monostable maps is investigated. The single map is associated with Poincare section of the FitzHugh-Nagumo neuron model. It is found that a diffusive coupling leads to the appearance of chaotic attractor. The attractor exists in an invariant region of phase space bounded by the manifolds of the saddle fixed point and the saddle periodic point. The oscillations from the chaotic attractor have a spike-burst shape with anti-phase synchronized spiking.     

A general algorithm for improving ILDMs

The method of intrinsic low-dimensional manifolds (ILDMs) has proven to be an efficient tool for the simplification of chemical kinetics. Nevertheless, there are still some open questions with respect to an efficient calculation and implementation of ILDMs. In this paper, we focus on the efficient calculation of ILDMs and present a refinement method for ILDMs. The method is based on an evolution equation of the manifolds towards a steady state solution which then represents the slow manifold. It has the property that it is continuous, differentiable, and in addition, it is an inertial manifold which represents the slow dynamics of the chemical system. In this way, many of the problems associated with the concept of ILDM are overcome.     

Bifurcation analysis of the travelling waveform of FitzHugh-Nagumo nerve conduction model equation

The FitzHugh-Nagumo model for travelling wave type neuron excitation is studied in detail. Carrying out a linear stability analysis near the equilibrium point, we bring out various interesting bifurcations which the system admits when a specific Z(2) symmetry is present and when it is not. Based on a center manifold reduction and normal form analysis, the Hopf normal form is deduced. The condition for the onset of limit cycle oscillations is found to agree well with the numerical results. We further demonstrate numerically that the system admits a period doubling route to chaos both in the presence as well as in the absence of constant external stimuli. (c) 1997 American Institute of Physics.     

Bifurcations of mixed-mode oscillations in a stellate cell model

Experimental recordings of the membrane potential of stellate cells within the entorhinal cortex show a transition from subthreshold oscillations (STOs) via mixed-mode oscillations (MMOs) to relaxation oscillations under increased injection of depolarizing current. Acker et al. introduced a 7D conductance based model which reproduces many features of the oscillatory patterns observed in these experiments. For the first time, we present a comprehensive bifurcation analysis of this model by using the software package AUTO. In particular, we calculate the stable MMO branches within the bifurcation diagram of this model, as well as other MMO patterns which are unstable. We then use geometric singular perturbation theory to demonstrate how the bifurcations are governed by a 3D reduced model introduced by Rotstein et al. We extend their analysis to explain all observed MMO patterns within the bifurcation diagram. A key role in this bifurcation analysis is played by a novel homoclinic bifurcation structure connecting to a saddle equilibrium on the unstable branch of the corresponding critical manifold. This type of homoclinic connection is possible due to canards of folded node (folded saddle-node) type.     

Instanton solutions from Abelian sinh-Gordon and Tzitzeica vortices

We study the Abelian Higgs vortex solutions to the sinh-Gordon equation and the elliptic Tzitzeica equation. Starting from these particular vortices, we construct solutions to the Taubes equation with higher vortex number, on surfaces with conical singularities.We then, analyse more general properties of vortices on such singular surfaces and propose a method to obtain vortices on conifolds from vortices on surfaces of revolution. We apply our method to construct explicit vortex solutions on the Poincaré disk with a conical singularity in the centre, to which we refer as the “hyperbolic cone”.We uplift the Abelian sinh-Gordon and Tzitzeica vortex solutions to four dimensions and construct cylindrically symmetric, self-dual Yang–Mills instantons on a non-self-dual (nor anti-self-dual) 4-dimensional Kähler manifold with non-vanishing scalar curvature. The instantons we construct in this way cannot be obtained via a twistorial approach.     

Minimal curvature trajectories: Riemannian geometry concepts for slow manifold computation in chemical kinetics

In dissipative ordinary differential equation systems different time scales cause anisotropic phase volume contraction along solution trajectories. Model reduction methods exploit this for simplifying chemical kinetics via a time scale separation into fast and slow modes. The aim is to approximate the system dynamics with a dimension-reduced model after eliminating the fast modes by enslaving them to the slow ones via computation of a slow attracting manifold. We present a novel method for computing approximations of such manifolds using trajectory-based optimization. We discuss Riemannian geometry concepts as a basis for suitable optimization criteria characterizing trajectories near slow attracting manifolds and thus provide insight into fundamental geometric properties of multiple time scale chemical kinetics. The optimization criteria correspond to a suitable mathematical formulation of “minimal relaxation” of chemical forces along reaction trajectories under given constraints. We present various geometrically motivated criteria and the results of their application to four test case reaction mechanisms serving as examples. We demonstrate that accurate numerical approximations of slow invariant manifolds can be obtained.     

Slow Motion and Metastability for a Nonlocal Evolution Equation

In this paper we consider a nonlocal evolution equation in one dimension, which describes the dynamics of a ferromagnetic system in the mean field approximation. In the presence of a small magnetic field, it admits two stationary and homogeneous solutions, representing the stable and metastable phases of the physical system. We prove the existence of an invariant, one dimensional manifold connecting the stable and metastable phases. This is the unstable manifold of a distinguished, spatially nonhomogeneous, stationary solution, called the critical droplet.^{(4, 10)} We show that the points on the manifold are droplets longer or shorter than the critical one, and that their motion is very slow in agreement with the theory of metastable patterns. We also obtain a new proof of the existence of the critical droplet, which is supplied with a local uniqueness result.     

Energy-Momentum and Equivalence Principle in Non-Riemannian Geometries

We introduce an energy-momentum density vector which is independent of the affine structure of the manifold and whose conservation is linked to observers. Integrating this quantity over time-like surfaces we can define Hamiltonian and momentum for the system which coincide with the corresponding ADM definitions for the case of irrotational Riemannian manifolds. As a consequence of our formalism, a Weak Equivalence Principle version for manifolds with torsion appears as the natural extension to non-Riemannian geometries from the Equivalence Principle of General Relativity.     

On the superpropagator of fields with exponential coupling

We define the vacuum expectation value of the time-ordered product of two exponentials of free fields as a distribution using minimal singularity as a criterion. The implication of this definition for an exponentially self-coupled scalar field is studied in second order of a perturbation expansion.     

Food chain chaos with canard explosion

The "tea-cup" attractor of a classical prey-predator-superpredator food chain model is studied analytically. Under the assumption that each species has its own time scale, ranging from fast for the prey to intermediate for the predator and to slow for the superpredator, the model is transformed into a singular perturbed system. It is demonstrated that the singular limit of the attractor contains a canard singularity. Singular return maps are constructed for which some subdynamics are shown to be equivalent to chaotic shift maps. Parameter regions in which the described chaotic dynamics exist are explicitly given.     

Bifurcation, amplitude death and oscillation patterns in a system of three coupled van der Pol oscillators with diffusively delayed velocity coupling

In this paper, we study a system of three coupled van der Pol oscillators that are coupled through the damping terms. Hopf bifurcations and amplitude death induced by the coupling time delay are first investigated by analyzing the related characteristic equation. Then the oscillation patterns of these bifurcating periodic oscillations are determined and we find that there are two kinds of critical values of the coupling time delay: one is related to the synchronous periodic oscillations, the other is related to eight branches of asynchronous periodic solutions bifurcating simultaneously from the zero solution. The stability of these bifurcating periodic solutions are also explicitly determined by calculating the normal forms on center manifolds, and the stable synchronous and stable phase-locked periodic solutions are found. Finally, some numerical simulations are employed to illustrate and extend our obtained theoretical results and numerical studies also describe the switches of stable synchronous and phase-locked periodic oscillations.     

On precise center stable manifold theorems for certain reaction–diffusion and Klein–Gordon equations

We consider positive, radial and exponentially decaying steady state solutions of the general reaction–diffusion and Klein–Gordon type equations and present an explicit construction of infinite-dimensional invariant manifolds in the vicinity of these solutions. The result is a precise stable manifold theorem for the reaction–diffusion equation and a precise center-stable manifold theorem for the Klein–Gordon equation, which include the co-dimension of the manifolds and the decay rates for even perturbations.     

Invariant manifold methods for metabolic model reduction

After the decay of transients, the behavior of a set of differential equations modeling a chemical or biochemical system generally rests on a low-dimensional surface which is an invariant manifold of the flow. If an equation for such a manifold can be obtained, the model has effectively been reduced to a smaller system of differential equations. Using perturbation methods, we show that the distinction between rapidly decaying and long-lived (slow) modes has a rigorous basis. We show how equations for attracting invariant (slow) manifolds can be constructed by a geometric approach based on functional equations derived directly from the differential equations. We apply these methods to two simple metabolic models. (c) 2001 American Institute of Physics.     

Synchronization of multi-frequency noise-induced oscillations

Using a model system of FitzHugh-Nagumo type in the excitable regime, the similarity between synchronization of self-sustained and noise-induced oscillations is studied for the case of more than one main frequency in the spectrum. It is shown that this excitable system undergoes the same frequency lockings as a self-sustained quasiperiodic oscillator. The presence of noise-induced both stable and unstable limit cycles and tori, as well as their tangential bifurcations, are discussed. As the FitzHugh-Nagumo oscillator represents one of the basic neural models, the obtained results are of high importance for neuroscience.     

Understanding anomalous delays in a model of intracellular calcium dynamics

In many cell types, oscillations in the concentration of free intracellular calcium ions are used to control a variety of cellular functions. It has been suggested [J. Sneyd et al., "A method for determining the dependence of calcium oscillations on inositol trisphosphate oscillations," Proc. Natl. Acad. Sci. U.S.A. 103, 1675-1680 (2006)] that the mechanisms underlying the generation and control of such oscillations can be determined by means of a simple experiment, whereby a single exogenous pulse of inositol trisphosphate (IP(3)) is applied to the cell. However, more detailed mathematical investigations [M. Domijan et al., "Dynamical probing of the mechanisms underlying calcium oscillations," J. Nonlinear Sci. 16, 483-506 (2006)] have shown that this is not necessarily always true, and that the experimental data are more difficult to interpret than first thought. Here, we use geometric singular perturbation techniques to study the dynamics of models that make different assumptions about the mechanisms underlying the calcium oscillations. In particular, we show how recently developed canard theory for singularly perturbed systems with three or more slow variables [M. Wechselberger, "A propos de canards (Apropos canards)," Preprint, 2010] applies to these calcium models and how the presence of a curve of folded singularities and corresponding canards can result in anomalous delays in the response of these models to a pulse of IP(3).     

On the number of solutions for the two-point boundary value problem on Riemannian manifolds

We study the solutions joining two fixed points of a time-independent dynamical system on a Riemannian manifold (M,g) from an enumerative point of view. We prove a finiteness result for solutions joining two points p,q∈M that are non-conjugate in a suitable sense, under the assumption that (M,g) admits a non-trivial convex function. We discuss in some detail the notion of conjugacy induced by a general dynamical system on a Riemannian manifold. Using techniques of infinite dimensional Morse theory on Hilbert manifolds we also prove that, under generic circumstances, the number of solutions joining two fixed points is odd. We present some examples where our theory applies.