Parameter identification in networks of dynamical systems

Mathematical models of real systems allow to simulate their behavior in conditions that are not easily or affordably reproducible in real life. Defining accurate models, however, is far from trivial and there is no one-size-fits-all solution. This thesis focuses on parameter identification in models of networks of dynamical systems, considering three case studies that fall under this umbrella: two of them are related to neural networks and one to power grids. The first case study is concerned with central pattern generators, i.e. small neural networks involved in animal locomotion. In this case, a design strategy for optimal tuning of biologically-plausible model parameters is developed, resulting in network models able to reproduce key characteristics of animal locomotion. The second case study is in the context of brain networks. In this case, a method to derive the weights of the connections between brain areas is proposed, utilizing both imaging data and nonlinear dynamics principles. The third and last case study deals with a method for the estimation of the inertia constant, a key parameter in determining the frequency stability in power grids. In this case, the method is customized to different challenging scenarios involving renewable energy sources, resulting in accurate estimations of this parameter

Locomotor Network Dynamics Governed By Feedback Control In Crayfish Posture And Walking

Sensorimotor circuits integrate biomechanical feedback with ongoing motor activity to produce behaviors that adapt to unpredictable environments. Reflexes are critical in modulating motor output by facilitating rapid responses. During posture, resistance reflexes generate negative feedback that opposes perturbations to stabilize a body. During walking, assistance reflexes produce positive feedback that facilitates fast transitions between swing and stance of each step cycle. Until recently, sensorimotor networks have been studied using biomechanical feedback based on external perturbations in the presence or absence of intrinsic motor activity. Experiments in which biomechanical feedback driven by intrinsic motor activity is studied in the absence of perturbation have been limited. Thus, it is unclear whether feedback plays a role in facilitating transitions between behavioral states or mediating different features of network activity independent of perturbation. These properties are important to understand because they can elucidate how a circuit coordinates with other neural networks or contributes to adaptable motor output. Computational simulations and mathematical models have been used extensively to characterize interactions of negative and positive feedback with nonlinear oscillators. For example, neuronal action potentials are generated by positive and negative feedback of ionic currents via a membrane potential. While simulations enable manipulation of system parameters that are inaccessible through biological experiments, mathematical models ascertain mechanisms that help to generate biological hypotheses and can be translated across different systems. Here, a three-tiered approach was employed to determine the role of sensory feedback in a crayfish locomotor circuit involved in posture and walking. In vitro experiments using a brain-machine interface illustrated that unperturbed motor output of the circuit was changed by closing the sensory feedback loop. Then, neuromechanical simulations of the in vitro experiments reproduced a similar range of network activity and showed that the balance of sensory feedback determined how the network behaved. Finally, a reduced mathematical model was designed to generate waveforms that emulated simulation results and demonstrated how sensory feedback can control the output of a sensorimotor circuit. Together, these results showed how the strengths of different approaches can complement each other to facilitate an understanding of the mechanisms that mediate sensorimotor integration

Understanding spiking and bursting electrical activity through piece-wise linear systems

In recent years there has been an increased interest in working with piece-wise linear caricatures of nonlinear models. Such models are often preferred over more detailed conductance based models for their small number of parameters and low computational overhead. Moreover, their piece-wise linear (PWL) form, allow the construction of action potential shapes in closed form as well as the calculation of phase response curves (PRC). With the inclusion of PWL adaptive currents they can also support bursting behaviour, though remain amenable to mathematical analysis at both the single neuron and network level. In fact, PWL models caricaturing conductance based models such as that of Morris-Lecar or McKean have also been studied for some time now and are known to be mathematically tractable at the network level. In this work we proceed to analyse PWL neuron models of conductance type. In particular we focus on PWL models of the FitzHugh-Nagumo type and describe in detail the mechanism for a canard explosion. This model is further explored at the network level in the presence of gap junction coupling. The study moves to a different area where excitable cells (pancreatic beta-cells) are used to explain insulin secretion phenomena. Here, Ca2+ signals obtained from pancreatic beta-cells of mice are extracted from image data and analysed using signal processing techniques. Both synchrony and functional connectivity analyses are performed. As regards to PWL bursting models we focus on a variant of the adaptive absolute IF model that can support bursting. We investigate the bursting electrical activity of such models with an emphasis on pancreatic beta-cells

Mean field modelling of human EEG: application to epilepsy

Aggregated electrical activity from brain regions recorded via an electroencephalogram (EEG), reveal that the brain is never at rest, producing a spectrum of ongoing oscillations that change as a result of different behavioural states and neurological conditions. In particular, this thesis focusses on pathological oscillations associated with absence seizures that typically affect 2–16 year old children. Investigation of the cellular and network mechanisms for absence seizures studies have implicated an abnormality in the cortical and thalamic activity in the generation of absence seizures, which have provided much insight to the potential cause of this disease. A number of competing hypotheses have been suggested, however the precise cause has yet to be determined. This work attempts to provide an explanation of these abnormal rhythms by considering a physiologically based, macroscopic continuum mean-field model of the brain's electrical activity. The methodology taken in this thesis is to assume that many of the physiological details of the involved brain structures can be aggregated into continuum state variables and parameters. The methodology has the advantage to indirectly encapsulate into state variables and parameters, many known physiological mechanisms underlying the genesis of epilepsy, which permits a reduction of the complexity of the problem. That is, a macroscopic description of the involved brain structures involved in epilepsy is taken and then by scanning the parameters of the model, identification of state changes in the system are made possible. Thus, this work demonstrates how changes in brain state as determined in EEG can be understood via dynamical state changes in the model providing an explanation of absence seizures. Furthermore, key observations from both the model and EEG data motivates a number of model reductions. These reductions provide approximate solutions of seizure oscillations and a better understanding of periodic oscillations arising from the involved brain regions. Local analysis of oscillations are performed by employing dynamical systems theory which provide necessary and sufficient conditions for their appearance. Finally local and global stability is then proved for the reduced model, for a reduced region in the parameter space. The results obtained in this thesis can be extended and suggestions are provided for future progress in this area

Mathematical frameworks for oscillatory network dynamics in neuroscience

The tools of weakly coupled phase oscillator theory have had a profound impact on the neuroscience community, providing insight into a variety of network behaviours ranging from central pattern generation to synchronisation, as well as predicting novel network states such as chimeras. However, there are many instances where this theory is expected to break down, say in the presence of strong coupling, or must be carefully interpreted, as in the presence of stochastic forcing. There are also surprises in the dynamical complexity of the attractors that can robustly appear—for example, heteroclinic network attractors. In this review we present a set of mathemat- ical tools that are suitable for addressing the dynamics of oscillatory neural networks, broadening from a standard phase oscillator perspective to provide a practical frame- work for further successful applications of mathematics to understanding network dynamics in neuroscience

Generalized reconfigurable memristive dynamical system (MDS) for neuromorphic applications

This study firstly presents (i) a novel general cellular mapping scheme for two dimensional neuromorphic dynamical systems such as bio-inspired neuron models, and (ii) an efficient mixed analog-digital circuit, which can be conveniently implemented on a hybrid memristor-crossbar/CMOS platform, for hardware implementation of the scheme. This approach employs 4n memristors and no switch for implementing an n-cell system in comparison with 2n2 memristors and 2n switches of a Cellular Memristive Dynamical System (CMDS). Moreover, this approach allows for dynamical variables with both analog and one-hot digital values opening a wide range of choices for interconnections and networking schemes. Dynamical response analyses show that this circuit exhibits various responses based on the underlying bifurcation scenarios which determine the main characteristics of the neuromorphic dynamical systems. Due to high programmability of the circuit, it can be applied to a variety of learning systems, real-time applications, and analytically indescribable dynamical systems. We simulate the FitzHugh-Nagumo (FHN), Adaptive Exponential (AdEx) integrate and fire, and Izhikevich neuron models on our platform, and investigate the dynamical behaviors of these circuits as case studies. Moreover, error analysis shows that our approach is suitably accurate. We also develop a simple hardware prototype for experimental demonstration of our approach.Unión Europea H2020 ECOMODE project under grant agreement 604102Unión Europea HBP project under grant number FP7-ICT-2013-FET-F-60410