Complex partial synchronization patterns in networks of delay-coupled neurons

We study the spatio-temporal dynamics of a multiplex network of delay-coupled FitzHugh–Nagumo oscillators with non-local and fractal connectivities. Apart from chimera states, a new regime of coexistence of slow and fast oscillations is found. An analytical explanation for the emergence of such coexisting partial synchronization patterns is given. Furthermore, we propose a control scheme for the number of fast and slow neurons in each layer.DFG, 163436311, SFB 910: Kontrolle selbstorganisierender nichtlinearer Systeme: Theoretische Methoden und Anwendungskonzept

Synchronization of coupled neural oscillators with heterogeneous delays

We investigate the effects of heterogeneous delays in the coupling of two excitable neural systems. Depending upon the coupling strengths and the time delays in the mutual and self-coupling, the compound system exhibits different types of synchronized oscillations of variable period. We analyze this synchronization based on the interplay of the different time delays and support the numerical results by analytical findings. In addition, we elaborate on bursting-like dynamics with two competing timescales on the basis of the autocorrelation function.Comment: 18 pages, 14 figure

Leaders do not look back, or do they?

We study the effect of adding to a directed chain of interconnected systems a directed feedback from the last element in the chain to the first. The problem is closely related to the fundamental question of how a change in network topology may influence the behavior of coupled systems. We begin the analysis by investigating a simple linear system. The matrix that specifies the system dynamics is the transpose of the network Laplacian matrix, which codes the connectivity of the network. Our analysis shows that for any nonzero complex eigenvalue $\lambda$ of this matrix, the following inequality holds: $\frac{|\Im \lambda |}{|\Re \lambda |} \leq \cot\frac{\pi}{n}$. This bound is sharp, as it becomes an equality for an eigenvalue of a simple directed cycle with uniform interaction weights. The latter has the slowest decay of oscillations among all other network configurations with the same number of states. The result is generalized to directed rings and chains of identical nonlinear oscillators. For directed rings, a lower bound $\sigma_c$ for the connection strengths that guarantees asymptotic synchronization is found to follow a similar pattern: $\sigma_c=\frac{1}{1-\cos\left( 2\pi /n\right)}$. Numerical analysis revealed that, depending on the network size $n$, multiple dynamic regimes co-exist in the state space of the system. In addition to the fully synchronous state a rotating wave solution occurs. The effect is observed in networks exceeding a certain critical size. The emergence of a rotating wave highlights the importance of long chains and loops in networks of oscillators: the larger the size of chains and loops, the more sensitive the network dynamics becomes to removal or addition of a single connection

Analysis for the hierarchical architecture of the heterogeneous FitzHugh-Nagumo network inducing synchronization

Synchronization is a key topic of research in neuroscience, medicine, and artificial neural networks; however, understanding its principle is difficult, both scientifically and mathematically. Specifically, the synchronization of the FitzHugh-Nagumo network with a hierarchical architecture has previously been studied; however, a mathematical analysis has not been conducted, owing to the network complexity. Therefore, in this paper, we saught to understand synchronization through mathematical analyses. In particular, we consider the most common types of hierarchical architecture and present a condition of the hierarchical architecture to induce synchronization. First, we provide mathematical analyses of a Lyapunov function for each layer, from which we obtain sufficient conditions guaranteeing synchronization and show that the Lyapunov function decreases exponentially. Moreover, we show that the internal connectivity critically affects synchronization in the first layer; however, in the second and subsequent layers, the internal connectivity is not important for synchronization, and the connectivity up to the first layer critically affects synchronization. We expect that the results and mathematical methodology can be applied to study other similar neural models with hierarchical architectures

Detecting partial synchrony in a complex oscillatory network using pseudo-vortices

Partial synchronization is characteristic phase dynamics of coupled oscillators on various natural and artificial networks, which can remain undetected due to the complexity of the systems. With an analogy between pairwise asynchrony of oscillators and topological defects, i.e., vortices, in the two-dimensional XY spin model, we propose a robust and data-driven method to identify the partial synchronization on complex networks. The proposed method is based on an integer matrix whose element is pseudo-vorticity that discretely quantifies asynchronous phase dynamics in every two oscillators, which results in graphical and entropic representations of partial synchrony. As a first trial, we apply our method to 200 FitzHugh-Nagumo neurons on a complex small-world network. Partially synchronized chimera states are revealed by discriminating synchronized states even with phase lags. Such phase lags also appear in partial synchronization in chimera states. Our topological, graphical, and entropic method is implemented solely with measurable phase dynamics data, which will lead to a straightforward application to general oscillatory networks including neural networks in the brain.Comment: 9 pages, 5 figure

Chimeras in physics and biology : Synchronization and desynchronization of rhythms

Rhythmen prägen unser Leben auf vielfältige Weise, z. B. durch Herzschlag und Atmung, oszillierende Gehirnströme, Lebenszyklen und Jahreszeiten, Uhren und Metronome, pulsierende Laser, Übertragung von Datenpaketen, und vieles andere. Die Physik komplexer nichtlinearer Systeme hat Methoden entwickelt, wie periodische Schwingungen und deren Synchronisation in komplexen Netzwerken, die aus vielen Bestandteilen zusammengesetzt sind, beschrieben und analysiert werden können. Synchronisierte Oszillationen, aber auch völlig desynchronisierte, chaotische Oszillationen spielen eine große Rolle in vielen Netzwerken in Natur und Technik. Beispielsweise ist das synchronisierte Feuern aller Neuronen im Gehirn ein pathologischer Zustand, etwa bei Epilepsie oder Parkinson-Erkrankung, und sollte unterdrückt werden, wie auch synchrone mechanische Schwingungen von Brücken. Andererseits ist die Synchronisation erwünscht beim stabilen Betrieb von Stromnetzen oder bei der verschlüsselten Kommunikation mit chaotischen Signalen. In Netzwerken aus identischen Komponenten können sich überraschenderweise auch spontan Hybrid-Zustände („Schimären“) bilden, die aus räumlich koexistierenden synchronisierten und desynchronisierten Bereichen bestehen, welche scheinbar nicht zusammen passen. Diese könnten relevant sein bei der Auslösung oder Beendigung epileptischer Anfälle, oder beim halbseitigen Schlaf einer Gehirnhälfte, der bei bestimmten Zugvögeln oder Säugetieren auftritt, oder beim kaskadenartigen Zusammenbruch des Stromnetzes.Rhythms influence our life in various ways, e.g., through heart beat and respiration, oscillating brain currents, life cycles and seasons, clocks and metronomes, pulsating lasers, transmission of data packets, and many others. The physics of complex nonlinear systems has developed methods to describe and analyze periodic oscillations and their synchronization in complex networks, which are composed of many components. Synchronized oscillations as well as completely asynchronous chaotic oscillations play a major role in many networks in nature and technology. For instance, the synchronous firing of all neurons in the brain represents a pathological state, like in epilepsy or Parkinson’s disease, and should be suppressed, as well as the synchronous mechanical vibration of bridges. On the other hand, synchronization is desirable for the stable operation of power grids or in encrypted communication with chaotic signals. In networks composed of identical components, intriguing hybrid states (“chimeras”) may form spontaneously, which consist of spatially coexisting synchronized and desynchronized domains, i.e., seemingly incongruous parts. This might be of relevance in inducing and terminating epileptic seizures, or in unihemispheric sleep which is found in certain migratory birds and mammals, or in cascading failures of the power grid.DFG, 163436311, Kontrolle selbstorganisierender nichtlinearer Systeme: Theoretische Methoden und AnwendungskonzepteDFG, 308748074, DFG-RSF: Komplexe dynamische Netzwerke: Effekte von heterogenen, adaptiven und zeitverzögerten Kopplunge

The multiplex decomposition: An analytic framework for multilayer dynamical networks

First Published in SIAM Journal on Applied Dynamical Systems in Volume 20, Issue 4 (2021), published by the Society for Industrial and Applied Mathematics (SIAM).Copyright © by SIAM. Unauthorized reproduction of this article is prohibited.Multiplex networks are networks composed of multiple layers such that the number of nodes in all layers is the same and the adjacency matrices between the layers are diagonal. We consider the special class of multiplex networks where the adjacency matrices for each layer are simultaneously triagonalizable. For such networks, we derive the relation between the spectrum of the multiplex network and the eigenvalues of the individual layers. As an application, we propose a generalized master stability approach that allows for a simplified, low-dimensional description of the stability of synchronized solutions in multiplex networks. We illustrate our result with a duplex network of FitzHugh--Nagumo oscillators. In particular, we show how interlayer interaction can lead to stabilization or destabilization of the synchronous state. Finally, we give explicit conditions for the stability of synchronous solutions in duplex networks of linear diffusive systems.DFG, 411803875, Dynamik gekoppelter Systeme mit Zeitverzögerungen und deren AnwendungenDFG, 440145547, Komplexe dynamische Netzwerke: Effekte von Heterogenität, Adaptivität und Topologie der Kopplunge

Neuronal oscillations: from single-unit activity to emergent dynamics and back

L’objectiu principal d’aquesta tesi és avançar en la comprensió del processament d’informació en xarxes neuronals en presència d’oscil lacions subumbrals. La majoria de neurones propaguen la seva activitat elèctrica a través de sinapsis químiques que són activades, exclusivament, quan el corrent elèctric que les travessa supera un cert llindar. És per aquest motiu que les descàrregues ràpides i intenses produïdes al soma neuronal, els anomenats potencials d’acció, són considerades la unitat bàsica d’informació neuronal, és a dir, el senyal mínim i necessari per a iniciar la comunicació entre dues neurones. El codi neuronal és entès, doncs, com un llenguatge binari que expressa qualsevol missatge (estímul sensorial, memòries, etc.) en un tren de potencials d’acció. Tanmateix, cap funció cognitiva rau en la dinàmica d’una única neurona. Circuits de milers de neurones connectades entre sí donen lloc a determinats ritmes, palesos en registres d’activitat colectiva com els electroencefalogrames (EEG) o els potencials de camp local (LFP). Si els potencials d’acció de cada cèl lula, desencadenats per fluctuacions estocàstiques de les corrents sinàptiques, no assolissin un cert grau de sincronia, no apareixeria aquesta periodicitat a nivell de xarxa. Per tal de poder entendre si aquests ritmes intervenen en el codi neuronal hem estudiat tres situacions. Primer, en el Capítol 2, hem mostrat com una cadena oberta de neurones amb un potencial de membrana intrínsecament oscil latori filtra un senyal periòdic arribant per un dels extrems. La resposta de cada neurona (pulsar o no pulsar) depèn de la seva fase, de forma que cada una d’elles rep un missatge filtrat per la precedent. A més, cada potencial d’acció presinàptic provoca un canvi de fase en la neurona postsinàptica que depèn de la seva posició en l’espai de fases. Els períodes d’entrada capaços de sincronitzar les oscil lacions subumbrals són aquells que mantenen la fase d’arribada dels potencials d’acció fixa al llarg de la cadena. Per tal de què el missatge arribi intacte a la darrera neurona cal, a més a més, que aquesta fase permeti la descàrrega del voltatge transmembrana. En segon cas, hem estudiat una xarxa neuronal amb connexions tant a veïns propers com de llarg abast, on les oscil lacions subumbrals emergeixen de l’activitat col lectiva reflectida en els corrents sinàptics (o equivalentment en el LFP). Les neurones inhibidores aporten un ritme a l’excitabilitat de la xarxa, és a dir, que els episodis en què la inhibició és baixa, la probabilitat d’una descàrrega global de la població neuronal és alta. En el Capítol 3 mostrem com aquest ritme implica l’aparició d’una bretxa en la freqüència de descàrrega de les neurones: o bé polsen espaiadament en el temps o bé en ràfegues d’elevada intensitat. La fase del LFP determina l’estat de la xarxa neuronal codificant l’activitat de la població: els mínims indiquen la descàrrega simultània de moltes neurones que, ocasionalment, han superat el llindar d’excitabilitat degut a un decreixement global de la inhibició, mentre que els màxims indiquen la coexistència de ràfegues en diferents punts de la xarxa degut a decreixements locals de la inhibició en estats globals d’excitació. Aquesta dinàmica és possible gràcies al domini de la inhibició sobre l’excitació. En el Capítol 4 considerem acoblament entre dues xarxes neuronals per tal d’estudiar la interacció entre ritmes diferents. Les oscil lacions indiquen recurrència en la sincronització de l’activitat col lectiva, de manera que durant aquestes finestres temporals una població optimitza el seu impacte en una xarxa diana. Quan el ritme de la població receptora i el de l’emissora difereixen significativament, l’eficiència en la comunicació decreix, ja que les fases de màxima resposta de cada senyal LFP no mantenen una diferència constant entre elles. Finalment, en el Capítol 5 hem estudiat com les oscil lacions col lectives pròpies de l’estat de son donen lloc al fenomen de coherència estocàstica. Per a una intensitat òptima del soroll, modulat per l’excitabilitat de la xarxa, el LFP assoleix una regularitat màxima donant lloc a un període refractari de la població neuronal. En resum, aquesta Tesi mostra escenaris d’interacció entre els potencials d’acció, característics de la dinàmica de neurones individuals, i les oscil lacions subumbrals, fruit de l’acoblament entre les cèl lules i ubiqües en la dinàmica de poblacions neuronals. Els resultats obtinguts aporten funcionalitat a aquests ritmes emergents, agents sincronitzadors i moduladors de les descàrregues neuronals i reguladors de la comunicació entre xarxes neuronals.The main objective of this thesis is to better understand information processing in neuronal networks in the presence of subthreshold oscillations. Most neurons propagate their electrical activity via chemical synapses, which are only activated when the electric current that passes through them surpasses a certain threshold. Therefore, fast and intense discharges produced at the neuronal soma (the action potentials or spikes) are considered the basic unit of neuronal information. The neuronal code is understood, then, as a binary language that expresses any message (sensory stimulus, memories, etc.) in a train of action potentials. Circuits of thousands of interconnected neurons give rise to certain rhythms, revealed in collective activity measures such as electroencephalograms (EEG) and local field potentials (LFP). Synchronization of action potentials of each cell, triggered by stochastic fluctuations of the synaptic currents, cause this periodicity at the network level.To understand whether these rhythms are involved in the neuronal code we studied three situations. First, in Chapter 2, we showed how an open chain of neurons with an intrinsically oscillatory membrane potential filters a periodic signal coming from one of its ends. The response of each neuron (to spike or not) depends on its phase, so that each cell receives a message filtered by the preceding one. Each presynaptic action potential causes a phase change in the postsynaptic neuron, which depends on its position in the phase space. Those incoming periods that are able to synchronize the subthreshold oscillations, keep the phase of arrival of action potentials fixed along the chain. The original message reaches intact the last neuron provided that this phase allows the discharge of the transmembrane voltage.I the second case, we studied a neuronal network with connections to both long range and close neighbors, in which the subthreshold oscillations emerge from the collective activity apparent in the synaptic currents. The inhibitory neurons provide a rhythm to the excitability of the network. When inhibition is low, the likelihood of a global discharge of the neuronal population is high. In Chapter 3 we show how this rhythm causes a gap in the discharge frequency of neurons: either they pulse single spikes or they fire bursts of high intensity. The LFP phase determines the state of the neuronal network, coding the activity of the population: its minima indicate the simultaneous discharge of many neurons, while its maxima indicate the coexistence of bursts due to local decreases of inhibition at global states of excitation. In Chapter 4 we consider coupling between two neural networks in order to study the interaction between different rhythms. The oscillations indicate recurrence in the synchronization of collective activity, so that during these time windows a population optimizes its impact on a target network. When the rhythm of the emitter and receiver population differ significantly, the communication efficiency decreases as the phases of maximum response of each LFP signal do not maintain a constant difference between them.Finally, in Chapter 5 we studied how oscillations typical of the collective sleep state give rise to stochastic coherence. For an optimal noise intensity, modulated by the excitability of the network, the LFP reaches a maximal regularity leading to a refractory period of the neuronal population.In summary, this Thesis shows scenarios of interaction between action potentials, characteristics of the dynamics of individual neurons, and the subthreshold oscillations, outcome of the coupling between the cells and ubiquitous in the dynamics of neuronal populations . The results obtained provide functionality to these emerging rhythms, triggers of synchronization and modulator agents of the neuronal discharges and regulators of the communication between neuronal networks

Transition from chimera/solitary states to traveling waves

We study numerically the spatiotemporal dynamics of a ring network of nonlocally coupled nonlinear oscillators, each represented by a two-dimensional discrete-time model of the classical van der Pol oscillator. It is shown that the discretized oscillator exhibits a richer behavior, combining the peculiarities of both the original system and its own dynamics. Moreover, a large variety of spatiotemporal structures is observed in the network of discrete van der Pol oscillators when the discretization parameter and the coupling strength are varied. Such regimes as the coexistence of multichimera state/traveling wave and solitary state are revealed for the first time and studied in detail. It is established that the majority of the observed chimera/solitary states, including the newly found ones, are transient towards the purely traveling wave mode. The peculiarities of the transition process and the lifetime (transient duration) of the chimera structures and the solitary state are analyzed depending on the system parameters, observation time, initial conditions, and influence of external noise