Harmonic analysis of oscillators through standard numerical continuation tools

In this paper, we describe a numerical continuation method that enables harmonic analysis of nonlinear periodic oscillators. This method is formulated as a boundary value problem that can be readily implemented by resorting to a standard continuation package - without modification - such as AUTO, which we used. Our technique works for any kind of oscillator, including electronic, mechanical and biochemical systems. We provide two case studies. The first study concerns itself with the autonomous electronic oscillator known as the Colpitts oscillator, and the second one with a nonlinear damped oscillator, a non-autonomous mechanical oscillator. As shown in the case studies, the proposed technique can aid both the analysis and the design of the oscillators, by following curves for which a certain constraint, related to harmonic analysis, is fulfilled.Comment: 20 pages, 4 figure

Phase analysis method for burst onset prediction

The response of bursting neurons to fluctuating inputs is usually hard to predict, due to their strong nonlinearity. For the same reason, decoding the injected stimulus from the activity of a bursting neuron is generally difficult. In this paper we propose a method describing (for neuron models) a mechanism of phase coding relating the burst onsets with the phase profile of the input current. This relation suggests that burst onset may provide a way for postsynaptic neurons to track the input phase. Moreover, we define a method of phase decoding to solve the inverse problem and estimate the likelihood of burst onset given the input state. Both methods are presented here in a unified framework, describing a complete coding-decoding procedure. This procedure is tested by using different neuron models, stimulated with different inputs (stochastic, sinusoidal, up, and down states). The results obtained show the efficacy and broad range of application of the proposed methods. Possible applications range from the study of sensory information processing, in which phase-of-firing codes are known to play a crucial role, to clinical applications such as deep brain stimulation, helping to design stimuli in order to trigger or prevent neural bursting

One-way dependent clusters and stability of cluster synchronization in directed networks

Cluster synchronization in networks of coupled oscillators is the subject of broad interest from the scientific community, with applications ranging from neural to social and animal networks and technological systems. Most of these networks are directed, with flows of information or energy that propagate unidirectionally from given nodes to other nodes. Nevertheless, most of the work on cluster synchronization has focused on undirected networks. Here we characterize cluster synchronization in general directed networks. Our first observation is that, in directed networks, a cluster A of nodes might be one-way dependent on another cluster B: in this case, A may remain synchronized provided that B is stable, but the opposite does not hold. The main contribution of this paper is a method to transform the cluster stability problem in an irreducible form. In this way, we decompose the original problem into subproblems of the lowest dimension, which allows us to immediately detect inter-dependencies among clusters. We apply our analysis to two examples of interest, a human network of violin players executing a musical piece for which directed interactions may be either activated or deactivated by the musicians, and a multilayer neural network with directed layer-to-layer connections.Comment: This is a preprint of an article published in Nature Communications. The final authenticated version is available online at: https://doi.org/10.1038/s41467-021-24363-7 or https://rdcu.be/cnya

Codimension-two homoclinic bifurcations underlying spike adding in the Hindmarsh-Rose burster

The well-studied Hindmarsh-Rose model of neural action potential is revisited from the point of view of global bifurcation analysis. This slow-fast system of three paremeterised differential equations is arguably the simplest reduction of Hodgkin-Huxley models capable of exhibiting all qualitatively important distinct kinds of spiking and bursting behaviour. First, keeping the singular perturbation parameter fixed, a comprehensive two-parameter bifurcation diagram is computed by brute force. Of particular concern is the parameter regime where lobe-shaped regions of irregular bursting undergo a transition to stripe-shaped regions of periodic bursting. The boundary of each stripe represents a fold bifurcation that causes a smooth spike-adding transition where the number of spikes in each burst is increased by one. Next, numerical continuation studies reveal that the global structure is organised by various curves of homoclinic bifurcations. In particular the lobe to stripe transition is organised by a sequence of codimension-two orbit- and inclination-flip points that occur along {\em each} homoclinic branch. Each branch undergoes a sharp turning point and hence approximately has a double-cover of the same curve in parameter space. The sharp turn is explained in terms of the interaction between a two-dimensional unstable manifold and a one-dimensional slow manifold in the singular limit. Finally, a new local analysis is undertaken using approximate Poincar\'{e} maps to show that the turning point on each homoclinic branch in turn induces an inclination flip that gives birth to the fold curve that organises the spike-adding transition. Implications of this mechanism for explaining spike-adding behaviour in other excitable systems are discussed.Comment: 32 pages, 18 figures, submitted to SIAM Journal on Applied Dynamical System

Dimensional reduction in networks of non- Markovian spiking neurons: Equivalence of synaptic filtering and heterogeneous propagation delays

Understanding the collective behavior of the intricate web of neurons composing a brain is one of the most challenging and complex tasks of modern neuroscience. Part of this complexity resides in the distributed nature of the interactions between the network components: for instance, the neurons transmit their messages (through spikes) with delays, which are due to different axonal lengths (i.e., communication distances) and/or noninstantaneous synaptic transmission. In developing effective network models, both of these aspects have to be taken into account. In addition, a satisfactory description level must be chosen as a compromise between simplicity and faithfulness in reproducing the system behavior. Here we propose a method to derive an effective theoretical description - validated through network simulations at microscopic level - of the neuron population dynamics in many different working conditions and parameter settings, valid for any synaptic time scale. In doing this we assume relatively small instantaneous fluctuations of the input synaptic current. As a by-product of this theoretical derivation, we prove analytically that a network with non-instantaneous synaptic transmission with fixed spike delivery delay is equivalent to a network characterized by a suited distribution of spike delays and instantaneous synaptic transmission, the latter being easier to treat

Co-Design of a Controller and Its Digital Implementation: The MOBY-DIC2 Toolbox for Embedded Model Predictive Control

Several software tools are available in the literature for the design and embedded implementation of linear model predictive control (MPC), both in its implicit and explicit (either exact or approximate) forms. Most of them generate C code for easy implementation on a microcontroller, and the others can convert the C code into hardware description language code for implementation on a field programmable gate array (FPGA). However, a unified tool allowing one to generate efficient embedded MPC for an FPGA, starting from the definition of the plant and its constraints, was still missing. The MOBY-DIC2 toolbox described in this brief bridges this gap. To illustrate its functionalities, the tool is exploited to embed the controller and observer for a real buck power converter in an FPGA. This implementation achieves a latency of about 30 μs with the implicit controller and 240 ns with the approximate explicit controlle