Synchronous behavior in networks of coupled systems : with applications to neuronal dynamics

Synchronization in networks of interacting dynamical systems is an interesting phenomenon that arises in nature, science and engineering. Examples include the simultaneous flashing of thousands of fireflies, the synchronous firing of action potentials by groups of neurons, cooperative behavior of robots and synchronization of chaotic systems with applications to secure communication. How is it possible that systems in a network synchronize? A key ingredient is that the systems in the network "communicate" information about their state to the systems they are connected to. This exchange of information ultimately results in synchronization of the systems in the network. The question is how the systems in the network should be connected and respond to the received information to achieve synchronization. In other words, which network structures and what kind of coupling functions lead to synchronization of the systems? In addition, since the exchange of information is likely to take some time, can systems in networks show synchronous behavior in presence of time-delays? The first part of this thesis focusses on synchronization of identical systems that interact via diffusive coupling, that is a coupling defined through the weighted difference of the output signals of the systems. The coupling might contain timedelays. In particular, two types of diffusive time-delay coupling are considered: coupling type I is diffusive coupling in which only the transmitted signals contain a time-delay, and coupling type II is diffusive coupling in which every signal is timedelayed. It is proven that networks of diffusive time-delay coupled systems that satisfy a strict semipassivity property have solutions that are ultimately bounded. This means that the solutions of the interconnected systems always enter some compact set in finite time and remain in that set ever after. Moreover, it is proven that nonlinear minimum-phase strictly semipassive systems that interact via diffusive coupling always synchronize provided the interaction is sufficiently strong. If the coupling functions contain time-delays, then these systems synchronize if, in addition to the sufficiently strong interaction, the product of the time-delay and the coupling strength is sufficiently small. Next, the specific role of the topology of the network in relation to synchronization is discussed. First, using symmetries in the network, linear invariant manifolds for networks of the diffusively time-delayed coupled systems are identified. If such a linear invariant manifold is also attracting, then the network possibly shows partial synchronization. Partial synchronization is the phenomenon that some, at least two, systems in the network synchronize with each other but not with every system in the network. It is proven that a linear invariant manifold defined by a symmetry in a network of strictly semipassive systems is attracting if the coupling strength is sufficiently large and the product of the coupling strength and the time-delay is sufficiently small. The network shows partial synchronization if the values of the coupling strength and time-delay for which this manifold is attracting differ from those for which all systems in the network synchronize. Next, for systems that interact via symmetric coupling type II, it is shown that the values of the coupling strength and time-delay for which any network synchronizes can be determined from the structure of that network and the values of the coupling strength and time-delay for which two systems synchronize. In the second part of the thesis the theory presented in the first part is used to explain synchronization in networks of neurons that interact via electrical synapses. In particular, it is proven that four important models for neuronal activity, namely the Hodgkin-Huxley model, the Morris-Lecar model, the Hindmarsh-Rose model and the FitzHugh-Nagumo model, all have the semipassivity property. Since electrical synapses can be modeled by diffusive coupling, and all these neuronal models are nonlinear minimum-phase, synchronization in networks of these neurons happens if the interaction is sufficiently strong and the product of the time-delay and the coupling strength is sufficiently small. Numerical simulations with various networks of Hindmarsh-Rose neurons support this result. In addition to the results of numerical simulations, synchronization and partial synchronization is witnessed in an experimental setup with type II coupled electronic realizations of Hindmarsh-Rose neurons. These experimental results can be fully explained by the theoretical findings that are presented in the first part of the thesis. The thesis continues with a study of a network of pancreatic -cells. There is evidence that these beta-cells are diffusively coupled and that the synchronous bursting activity of the network is related to the secretion of insulin. However, if the network consists of active (oscillatory) beta-cells and inactive (dead) beta-cells, it might happen that, due to the interaction between the active and inactive cells, the activity of the network dies out which results in a inhibition of the insulin secretion. This problem is related to Diabetes Mellitus type 1. Whether the activity dies out or not depends on the number of cells that are active relative to the number of inactive cells. A bifurcation analysis gives estimates of the number of active cells relative to the number of inactive cells for which the network remains active. At last the controlled synchronization problem for all-to-all coupled strictly semipassive systems is considered. In particular, a systematic design procedure is presented which gives (nonlinear) coupling functions that guarantee synchronization of the systems. The coupling functions have the form of a definite integral of a scalar weight function on a interval defined by the outputs of the systems. The advantage of these coupling functions over linear diffusive coupling is that they provide high gain only when necessary, i.e. at those parts of the state space of the network where nonlinearities need to be suppressed. Numerical simulations in networks of Hindmarsh-Rose neurons support the theoretical results

A stability-theory perspective to synchronisation of heterogeneous networks

Dans ce mémoire, nous faisons une présentation de nos recherches dans le domaine de la synchronisation des systèmes dynamiques interconnectés en réseau. Une des originalités de nos travaux est qu'ils portent sur les réseaux hétérogènes, c'est à dire, des systèmes à dynamiques diverses. Au centre du cadre d'analyse que nous proposons, nous introduisons le concept de dynamique émergente. Il s'agit d'une dynamique "moyennée'' propre au réseau lui-même. Sous l'hypothèse qu'il existe un attracteur pour cette dynamique, nous montrons que le problème de synchronisation se divise en deux problèmes duaux : la stabilité de l'attracteur et la convergence des trajectoires de chaque système vers celles générées par la dynamique émergente. Nous étudions aussi le cas particulier des oscillateurs de Stuart-Landau

Network Identification: A Passivity and Network Optimization Approach

The theory of network identification, namely identifying the interaction topology among a known number of agents, has been widely developed for linear agents over recent years. However, the theory for nonlinear agents remains less extensive. We use the notion maximal equilibrium-independent passivity (MEIP) and network optimization theory to present a network identification method for nonlinear agents.We do so by introducing a specially designed exogenous input, and exploiting the properties of networked MEIP systems. We then specialize on LTI agents, showing that the method gives a distributed cubic-time algorithm for network reconstruction in that case. We also discuss different methods of choosing the exogenous input, and provide an example on a neural network model.Comment: 8 Pages, 3 Figure

Synchronization of complex dynamical networks with fractional order

Complex dynamical networks (CDN) can be applied to many areas in real world, from medicine, biology, Internet to sociology. Study on CDNs has drawn great attention in recent years. Nodes in a CDN can be modelled as systems represented by differential equations. Study has shown that fractional order differential equations (DF) can better represent some real world systems than integer-order DFs. This research work focuses on synchronization in fractional CDNs.  A literature review on CDNs with fractional order has summarized the latest works in this area.  Fractional chaotic systems are studied in our initial investigation.  Fractional calculus is introduced and the relevant fundamentals to model, describe and analyse dynamical networks are presented. It is shown that the structure and topological characteristics of a network can have a big impact on its synchronizability. Synchronizability and its various interpretations in dynamical networks are studied. To synchronize a CDN efficiently, controllers are generally needed. Controller design is one of the main tasks in this research. Our first design is a new sliding mode control to synchronize a dynamical network with two nodes. Its stability has been proven and verified by simulations.  Its convergence speed outperforms Vaidyanathan's scheme, a well-recognized scheme in this area. The design can be generalized to CDNs with more nodes.  As many applications can be modelled as CDNs with node clustering, a different sliding mode control is designed for cluster synchronization of a CDN with fractional order. Its stability is proven by using Lyapunov method. Its convergence and efficiency is shown in a simulation. Besides these nonlinear methods mentioned, linear control is also studied intensively for the synchronization.  A novel linear method for synchronization of fractional CDNs using a new fractional Proportional-Integral (PI) pinning control is proposed.  Its stability is proven and the synchronization criteria are obtained. The criteria have been simplified using two corollaries so the right value for the variables can be easily assigned. The proposed method is compared with the conventional linear method which uses Proportional (P) controller. In the comparison, the mean squared error function is used. The function measures the average of the squared errors and it is an instant indicator of the synchronization efficiency. A numerical simulation is repeated 100 times to obtain the averages over these runs. Each simulation has different random initial values for both controllers. The average of the errors in all the 100 simulations is obtained and the area under the function curve is defined as an overall performance index (OPI), which indicates the controller's overall performance. In control, small overshoot is always desired. In our work, the error variation is also used as a measure.  The maximum variation from the average of 100 simulations is calculated and compared for both methods. With all the statistical comparisons, it is clear that with the same power consumption, the proposed method outperforms the conventional one and achieves faster and smoother synchronization. Communication constraints exist in most real world CDNs. Communication constraints and their impact on control and synchronization of CDNs with fractional order are investigated in our study. A new adaptive method for synchronizing fractional CDN with disturbance and uncertainty is designed. Its stability is proven and its synchronization criteria are obtained for both fractional CDN with known and unknown parameters. Random disturbance is also included in both cases. Our results show that the new method is efficient in synchronizing CDNs with presence of both disturbance and uncertainty

Putting reaction-diffusion systems into port-Hamiltonian framework

Reaction-diffusion systems model the evolution of the constituents distributed in space under the influence of chemical reactions and diffusion [6], [10]. These systems arise naturally in chemistry [5], but can also be used to model dynamical processes beyond the realm of chemistry such as biology, ecology, geology, and physics. In this paper, by adopting the viewpoint of port-controlled Hamiltonian systems [7] we cast reaction-diffusion systems into the portHamiltonian framework. Aside from offering conceptually a clear geometric interpretation formalized by a Stokes-Dirac structure [8], a port-Hamiltonian perspective allows to treat these dissipative systems as interconnected and thus makes their analysis, both quantitative and qualitative, more accessible from a modern dynamical systems and control theory point of view. This modeling approach permits us to draw immediately some conclusions regarding passivity and stability of reaction-diffusion systems. It is well-known that adding diffusion to the reaction system can generate behaviors absent in the ode case. This primarily pertains to the problem of diffusion-driven instability which constitutes the basis of Turing’s mechanism for pattern formation [11], [5]. Here the treatment of reaction-diffusion systems as dissipative distributed portHamiltonian systems could prove to be instrumental in supply of the results on absorbing sets, the existence of the maximal attractor and stability analysis. Furthermore, by adopting a discrete differential geometrybased approach [9] and discretizing the reaction-diffusion system in port-Hamiltonian form, apart from preserving a geometric structure, a compartmental model analogous to the standard one [1], [2] is obtaine