Geometric Properties of Isostables and Basins of Attraction of Monotone Systems

In this paper, we study geometric properties of basins of attraction of monotone systems. Our results are based on a combination of monotone systems theory and spectral operator theory. We exploit the framework of the Koopman operator, which provides a linear infinite-dimensional description of nonlinear dynamical systems and spectral operator-theoretic notions such as eigenvalues and eigenfunctions. The sublevel sets of the dominant eigenfunction form a family of nested forward-invariant sets and the basin of attraction is the largest of these sets. The boundaries of these sets, called isostables, allow studying temporal properties of the system. Our first observation is that the dominant eigenfunction is increasing in every variable in the case of monotone systems. This is a strong geometric property which simplifies the computation of isostables. We also show how variations in basins of attraction can be bounded under parametric uncertainty in the vector field of monotone systems. Finally, we study the properties of the parameter set for which a monotone system is multistable. Our results are illustrated on several systems of two to four dimensions.Comment: 12 pages, to appear in IEEE Transaction on Automatic Contro

Computing an Inner and an Outer Approximation of the Viability Kernel

International audienceThe viability kernel corresponds to the set of all state vectors of a controlled dynamic system that are viable, i.e., such that there exists an input such that the system will not enter inside a forbidden zone. In this paper, we propose a method which computes an inner and an outer approximation of the viability kernel in a guaranteed way. Our method is based on interval analysis and uses the notions of V-viability and capture basin. We illustrate our approach on the car on the hill problem. A software package has been developed to solve any 2D-problem

Boundary effects on localized structures in spatially extended systems

We present a general method of analyzing the influence of finite size and boundary effects on the dynamics of localized solutions of non-linear spatially extended systems. The dynamics of localized structures in infinite systems involve solvability conditions that require projection onto a Goldstone mode. Our method works by extending the solvability conditions to finite sized systems, by incorporating the finite sized modifications of the Goldstone mode and associated nonzero eigenvalue. We apply this method to the special case of non-equilibrium domain walls under the influence of Dirichlet boundary conditions in a parametrically forced complex Ginzburg Landau equation, where we examine exotic nonuniform domain wall motion due to the influence of boundary conditions.Comment: 9 pages, 5 figures, submitted to Physical Review

Metastability in a stochastic neural network modeled as a velocity jump Markov process

One of the major challenges in neuroscience is to determine how noise that is present at the molecular and cellular levels affects dynamics and information processing at the macroscopic level of synaptically coupled neuronal populations. Often noise is incorprated into deterministic network models using extrinsic noise sources. An alternative approach is to assume that noise arises intrinsically as a collective population effect, which has led to a master equation formulation of stochastic neural networks. In this paper we extend the master equation formulation by introducing a stochastic model of neural population dynamics in the form of a velocity jump Markov process. The latter has the advantage of keeping track of synaptic processing as well as spiking activity, and reduces to the neural master equation in a particular limit. The population synaptic variables evolve according to piecewise deterministic dynamics, which depends on population spiking activity. The latter is characterised by a set of discrete stochastic variables evolving according to a jump Markov process, with transition rates that depend on the synaptic variables. We consider the particular problem of rare transitions between metastable states of a network operating in a bistable regime in the deterministic limit. Assuming that the synaptic dynamics is much slower than the transitions between discrete spiking states, we use a WKB approximation and singular perturbation theory to determine the mean first passage time to cross the separatrix between the two metastable states. Such an analysis can also be applied to other velocity jump Markov processes, including stochastic voltage-gated ion channels and stochastic gene networks

Effective Fokker-Planck Equation for Birhythmic Modified van der Pol Oscillator

We present an explicit solution based on the phase-amplitude approximation of the Fokker-Planck equation associated with the Langevin equation of the birhythmic modified van der Pol system. The solution enables us to derive probability distributions analytically as well as the activation energies associated to switching between the coexisting different attractors that characterize the birhythmic system. Comparing analytical and numerical results we find good agreement when the frequencies of both attractors are equal, while the predictions of the analytic estimates deteriorate when the two frequencies depart. Under the effect of noise the two states that characterize the birhythmic system can merge, inasmuch as the parameter plane of the birhythmic solutions is found to shrink when the noise intensity increases. The solution of the Fokker-Planck equation shows that in the birhythmic region, the two attractors are characterized by very different probabilities of finding the system in such a state. The probability becomes comparable only for a narrow range of the control parameters, thus the two limit cycles have properties in close analogy with the thermodynamic phases

On Norm-Based Estimations for Domains of Attraction in Nonlinear Time-Delay Systems

For nonlinear time-delay systems, domains of attraction are rarely studied despite their importance for technological applications. The present paper provides methodological hints for the determination of an upper bound on the radius of attraction by numerical means. Thereby, the respective Banach space for initial functions has to be selected and primary initial functions have to be chosen. The latter are used in time-forward simulations to determine a first upper bound on the radius of attraction. Thereafter, this upper bound is refined by secondary initial functions, which result a posteriori from the preceding simulations. Additionally, a bifurcation analysis should be undertaken. This analysis results in a possible improvement of the previous estimation. An example of a time-delayed swing equation demonstrates the various aspects.Comment: 33 pages, 8 figures, "This is a pre-print of an article published in 'Nonlinear Dynamics'. The final authenticated version is available online at https://doi.org/10.1007/s11071-020-05620-8

Investigating a simple model of cutaneous wound healing angiogenesis

A simple model of wound healing angiogenesis is presented, and investigated using numerical and asymptotic techniques. The model captures many key qualitative features of the wound healing angiogenic response, such as the propagation of a structural unit into the wound centre. A detailed perturbative study is pursued, and is shown to capture all features of the model. This enables one to show that the level of the angiogenic response predicted by the model is governed to a good approximation by a small number of parameter groupings. Further investigation leads to predictions concerning how one should select between potential optimal means of stimulating cell proliferation in order to increase the level of the angiogenic response

Capture basin approximation using interval analysis

This paper proposes a new approach for computing the capture basin C of a target T. The capture basin corresponds to the set of initial state vectors such that the target could be reached in finite time via an appropriate control input, before possibly leaving the target. Whereas classical capture basin characterization does not provide any guarantee on the set of state vectors that belong to the capture basin, interval analysis and guaranteed numerical integration allow us to avoid any indetermination. We present an algorithm that is able to provide guaranteed approximation of the inner C− and the outer C+ of the capture basin, such that C−⊆C⊂C+. In order to illustrate the principle and the efficiency of the approach, a testcase on the ‘car on the hill’ problem is provided. Copyright © 2010 John Wiley & Sons, Ltd