Network-driven synchronisation of phase-coupled oscillators
|
|
- Liliana Evans
- 6 years ago
- Views:
Transcription
1 Network-driven synchronisation of phase-coupled oscillators University of Exeter CEMPS Exeter, 4 th Feb 2014
2 Outline Motivation The Classical Kuramoto Model
3 A normal EEG Fp1 AVG Fp2 AVG F3 AVG F4 AVG C3 AVG C4 AVG P3 AVG P4 AVG O1 AVG O2 AVG F7 AVG F8 AVG T3 AVG T4 AVG T5 AVG T6 AVG Fz AVG Cz AVG Pz AVG Scale
4 EEG recording of a generalised seizure Fp1 Fp2 F7 F3 Fz F4 F8 A1 T3 C3 Cz Electrode C4 T4 A2 T5 P3 Pz P4 T6 O1 O2 Poly1 Poly2 Poly Time (ms) x 10 4
5 EEG recording of a generalised seizure Fp1 Fp2 F7 F3 Fz F4 F8 A1 T3 C3 Cz Electrode C4 T4 A2 T5 P3 Pz P4 T6 O1 O2 Poly1 Poly2 Poly Time (ms) x 10 4
6 Motivation Epilepsy - A Network Problem!
7 phase and moto net- 5). tatively in isotropic ted by the anatomializing on the ampliity in the hem to the nnectivity the funcronization entioned nizability chard and scopically tural cons of averes appears tradition ilson and an, 1975; 8). Neural a rhythm, transition t al., 1974; 999; David of two or dam et al., 01; David networks al., 2010). ass model scriptions ral models token, the where N e and N i are the numbers c of excitatory and inhibitory Motivation neurons. By this averaging, E and The I represent Classical the mean Kuramoto firing rates Model of dt N l = 1 d h EE wk = vk + sin( wl wk) all excitatory and inhibitory neurons, Networks respectively, of Kuramoto of the neural Models For the sake of legibility, however, we here refer to (2) also as population in question, i.e., that at node Application k. to EEG Data the Kuramoto model. Within that population, P every E k k neuron receives Cinput ki E I from all As mentioned above, the effect of increasing h in the isotropic other neurons of the population. Furthermore, the excitatory units case is to increase the phase synchrony amongst the oscillators. individually receive constant external inputs p n, whose average is 1 N Suppose the coupling is weak (i.e., smaller than the critical value, e given by P = N p e n = 1 n. The c EI sum of all cinputs IE is (instantaneously) or h h integrated in time when it exceeds some threshold u. This thresholding is realized by means of a sigmoid function S. Without loss of c c, then the oscillators phases disperse, whereas for strong coupling h h the oscillators become synchronous, i.e., the phases are locked at fixed differences. In the intermediate case h h generality we here chose S[x] = (1 + e x ) 1 ; we note that, in general, c, clusters of synchronous oscillators may emerge. However, many other I k the thresholds may differ between excitatory and inhibitory units 1. oscillators, whose natural ORIGINAL frequencies RESEARCH ARTICLE are at the tails of the distribution, are not locked into a cluster. doi: /fninf In other words, as h increases, the In consequence, the mean NEUROINFORMATICS firing rates of the neural populations can published: 15 July 2011 be cast in the following dynamical system c II interaction functions overcome the dispersion of natural frequencies d On the influence of amplitude on the connectivity between dt E E S a c E phases c I P v n resulting in a transition from incoherence, to partial and then = + E( EE IE ue + ) FIGURE 1 Network of Wilson Cowan models. At each node k a neural full synchronization (Acebron et al., 2005; Breakspear et al., 2010). population containing excitatory and inhibitory units (E d Andreas Daffertshofer* and Bernadette C. M. van Wijk dt I I S a c E c I k and I k, respectively) = yields + self-sustained I( oscillations. EI II uother I) nodes are connected to the excitatory Linking neural mass models to phase oscillators unit by means of ΣC kl E l. Note that this (mean-field) coupling is scaled by a Research Institute MOVE, VU University Amsterdam, Amsterdam, Netherlands When deriving the Kuramoto network from the Wilson Cowan oscillator network, the major ingredient is to average every oscillator over scalar h see Eq. 1 for details. The characteristics of this dynamical system range from a mere Edited by: In recent studies, functional connectivities have been reported to display characteristics of fixed-point relaxation to Olaf limit Sporns, Indiana cycle University, oscillations USA (self-sustained one cycle when assuming that its amplitude and phase change slowly complex networks that have been suggested concur with those of the underlying structural, oscillations) depending on Reviewed parameter by: settings (Wilson i.e., anatomical, and Cowan, as compared to the oscillator s frequency. That is, time-dependent oscillators. The Kuramoto networks. Do functional networks always agree with structural ones? In all Joana R. B. model Cabral, Universitat and its Pompeu link to the here-discussed 1972), in particular on the Fabra, choice Spain of the external input generality, P. That this input question can amplitude be answered and with phase no : are for fixed, instance, the a system fully synchronized is integrated state over one period network of Wilson Cowan models will be briefly sketched in the Juergen Kurths, Humboldt Universität, would imply isotropic homogeneous functional connections irrespective of the real underlying is usually chosen at random. In the current study, we restrict all to remove all harmonic oscillations, and, subsequently, amplitude and following two sub-sections. Germany structure. A proper inference of structure from function and vice versa requires more than a parameter values to the regime *Correspondence: within which the dynamics sole focus on displays phase are again considered to be time-dependent (Guckenheimer and phase synchronization. We show that functional connectivity critically depends Andreas Daffertshofer, Research self-sustained oscillations; see Appendix. on amplitude variations, which Holmes, implies 1990) that, in general, this procedure phase patterns is also should referred be analyzed to as a combination of Kuramoto network of phase Institute MOVE, oscillators VU University To combine Wilson Cowan Amsterdam, models Van der Boechorststraat in a network, 9, conjunction different with populations are now connected a.daffertshofer@vu.nl via their excitatory units of phase by oscillators. virtue of mation (Haken, 1974). As shown in more detail in the Appendix, the the corresponding rotating amplitude. wave approximation We discuss this issue and slowly by comparing varying the amplitude phase approxi- The collective behavior 1081 of BT Amsterdam, a network Netherlands. of oscillators, synchronization whose states patterns are of interconnected Wilson Cowan models vis-à-vis Kuramoto networks captured by a single scalar phase w For the interconnected Wilson Cowan models we derive analytically how the sum of all E l in the dynamics of E k each, can, in first approximation, k (see Figure 1). connectivity The dynamics phase dynamics of the system (1) can in this way be approximated as between phases explicitly depends on the generating oscillators amplitudes. be represented by the set of N coupled differential Eq. at node k then becomes In consequence, the link between neurophysiological N studies and computational models d h always requires the incorporation of the amplitude dynamics. Supplementing synchronization Ncharacteristics by amplitude patterns, as captured by, e.g., spectral power in M/EEG recordings, d will certainly aid our understanding of the relation between structural and functional organizations dt E E S a c E c I P h N CE dt N CaS Rl wk = vk + ( 0) N kl E x Ek, sin wl wk d h ( ) w 2 l = 1 Rk k = vk + Dkl sin( wl wk) (2) dt k = k + NE k k E + k + kl l i = 1 EE IE u in neural networks at large. l N h = 1 Keywords: connectivity, phase synchronization, Kuramoto N Ca 3 ( 0) kl ES 16 network, Wilson Cowan x Ek, model, RkRl( ( cee + 3c IE ) sin( wl wk) amplitude dependency d That is, the k-th oscillator, with natural frequency v l = 1 k, adjusts its Iphase k = Iaccording k + S ai( ceiek ciiik I) dt to input from other u (1) INTRODUCTION oscillators through a pair-wise + 2c discriminating between frequency and phase-locked activity EEcIE cos( w (Baker l wk)) phase interaction function The interplay sin(w between l wstructural k ). The and connectivity functional brain matrix = networks has et al., 1999; Mima and Hallett, 1999; Salinas and Sejnowski, 2001; In D kl words, is again all Wilson Cowan scaled by become overall a oscillators, popular coupling topic of located research strength, at in recent nodes h. years. As l in will It the is currently be with Fries, S 2005; and Womelsdorf S referring et al., to 2007). the It first is usually and assumed third derivative that of the network drive the change believed of the that firing the topologies rate of structural the excitatory and functional units networks in sigmoid amplitude function or power S. variations The parameter take place on long x ( 0) sketched below, h serves as a bifurcation parameter in that small Ek, time is given scales when by various empirical systems may disagree (Sporns and Kötter, 2004) compared to the phase dynamics and are therefore considered negligible. The coupling that does, or does not, yield synchrony N between E k. The values connectivity of h yield is a given network but systematic by the behavior analyses real-valued tackling that essentially matrix this issue are C kl agrees few that and has far with between. 0 h vanishing diagonal elements, In a combined i.e., Cneural mass and graph theoretical model of electroencephalographic signals, it was found that patterns of functional whether this assumption is valid, and by this, tackle if a sole focus oscillators hence 0exclusively 0 depends on the phase. Here we ask kk = 0. That connectivity matrix ( ) ( ) ( ) xek ae c Ek c Ik ue Pk l is scaled via the overall coupling strength h. It is important to note N CE ( 0 the entirely uncoupled case (i.e., the phases are not synchronized), ), = + + EE IE kl l = 1 whereas h larger than a certain critical value h connectivity are influenced by but c causes the phases not identical to those of on phase really covers all functional characteristics of networks. In that to the synchronize. C connectivity The the frequencies matrix corresponding is here v k structural are always distributed Helmut level identified (Ponten according Schmidt et al., as 2010). the to In this the present ( 0) ( study 0) we describe the dynamics of neural populations Ubiquitousness of the Kuramoto Model
8 The Kuramoto Model The model: θ j = ω j + K N N sin(θ n θ j ), j = {1,..., N} n=1 The order parameter r defines the degree of synchronicity : re iψ = 1 N N e iθn = θ j = ω j + rk sin(ψ θ j ), n=1 Thermodynamic limit (N ), g(ω) = 1 π e (Ω ω)2 : r = π rk ( ) ( )) /2 (rk) 2 (rk) 2 (I 2 e (rk)2 0 + I 1 = F (rk) 2 2
9 The Kuramoto Model Synchronisation sets in at critical coupling strength K c : r, f(rk) K=2 K=1.5 1 r K=1 K= r K C K 3
10 The Kuramoto model
11 Directed Networks of Kuramoto Populations :
12 Directed Networks of Kuramoto Populations : The full model with adjacency matrix Â: θ p j = ω p j + K p N N P sin(θn p θ p j ) + n=1 s=1 A s,p N N sin(θn s θ p j ), n=1 define local order parameters: = r p e iψp = 1 N N e iθp n, n=1 θ p j P = ω p j + r p K p sin(ψ p θ p j ) + r s A s,p sin(ψ s θ p j ). s=1
13 Directed Networks of Kuramoto Populations : under the assumption that ψ p = ψ s p, s we obtain r p = F ( K p r p + ) P A s,p r s, p = {1,..., P}. s=1 linearisation about (r 1,..., r P ) = (0,..., 0) yields ( ) Kp 1 r p + 1 K c K c P m=1 A s,p r s = 0. nontrivial solutions exist if ( ) Â det + ˆK = 0, K p,p = K p 1, K c K c K s,p = 0 if s p
14 Directed Networks of Kuramoto Populations : rewriting  = C ˆρ yields a critical value for the global coupling parameter C: ( ) Cc det ˆρ + ˆK = 0. K c if ˆρ is triangular, ( ) Cc ( ) det ˆρ + ˆK = det ˆK 0 K c
15 2 populations
16 cycles
17 a loop of 3 populations plus one extra connection r 1 = F (r 1 K + r 2 C + r 3 C) r 2 = F (r 2 K + r 1 C) r 3 = F (r 3 K + r 2 C) 1 r K = C 1.8 2
18 a hierarchical network of 3: r 1 = F (r 1 K + r 2 C + r 3 C) r 2 = F (r 2 K) r 3 = F (r 3 K + r 2 C) 1 r K = C 1.8 2
19 Larger Networks
20 Larger Networks Network-driven synchronisation of phase-coupled oscillators
21 Larger Networks Helmut Schmidt
22 Motivation The Classical Kuramoto Model Larger Networks
23 Directed Networks from EEG Data To compute (delayed) cross-correlation between time-series x i and x j we follow Bialonski & Lehnertz (2013): with ρ i,j = max τ ξ(x i, x j )(τ) = { } ξ(x i, x j )(τ), ξ(xi, x i )(0)ξ(x j, x j )(0) { T τ t=1 x i(t + τ)x j (t) if τ 0 T t=1+τ x i(t + τ)x j (t) if τ < 0. Fp1 AVG Fp2 AVG F3 AVG F4 AVG C3 AVG C4 AVG P3 AVG P4 AVG O1 AVG O2 AVG F7 AVG F8 AVG T3 AVG T4 AVG T5 AVG T6 AVG Fz AVG Cz AVG Pz AVG Scale
24 Motivation The Classical Kuramoto Model Networks of Kuramoto Models Directed Networks from EEG Data
25 Critical Coupling Constant C c critical coupling constant controls patients ** * 0.2 low high delta theta beta gamma alpha alpha
26 Seizure Onset Zones P 2 P 3 P 1 P 4 P 5 P 6 P 7
27 Seizure Onset Zones
Phase Synchronization
Phase Synchronization Lecture by: Zhibin Guo Notes by: Xiang Fan May 10, 2016 1 Introduction For any mode or fluctuation, we always have where S(x, t) is phase. If a mode amplitude satisfies ϕ k = ϕ k
More informationSynchronization Transitions in Complex Networks
Synchronization Transitions in Complex Networks Y. Moreno 1,2,3 1 Institute for Biocomputation and Physics of Complex Systems (BIFI) University of Zaragoza, Zaragoza 50018, Spain 2 Department of Theoretical
More informationConsider the following spike trains from two different neurons N1 and N2:
About synchrony and oscillations So far, our discussions have assumed that we are either observing a single neuron at a, or that neurons fire independent of each other. This assumption may be correct in
More informationOscillator synchronization in complex networks with non-uniform time delays
Oscillator synchronization in complex networks with non-uniform time delays Jens Wilting 12 and Tim S. Evans 13 1 Networks and Complexity Programme, Imperial College London, London SW7 2AZ, United Kingdom
More informationHysteretic Transitions in the Kuramoto Model with Inertia
Rostock 4 p. Hysteretic Transitions in the uramoto Model with Inertia A. Torcini, S. Olmi, A. Navas, S. Boccaletti http://neuro.fi.isc.cnr.it/ Istituto dei Sistemi Complessi - CNR - Firenze, Italy Istituto
More informationChimera states in networks of biological neurons and coupled damped pendulums
in neural models in networks of pendulum-like elements in networks of biological neurons and coupled damped pendulums J. Hizanidis 1, V. Kanas 2, A. Bezerianos 3, and T. Bountis 4 1 National Center for
More informationarxiv: v1 [math.ds] 22 Jan 2019
Noname manuscript No. (will be inserted by the editor The uncoupling limit of identical Hopf bifurcations with an application to perceptual bistability Alberto Pérez-Cervera Peter Ashwin Gemma Huguet Tere
More informationSaturation of Information Exchange in Locally Connected Pulse-Coupled Oscillators
Saturation of Information Exchange in Locally Connected Pulse-Coupled Oscillators Will Wagstaff School of Computer Science, Georgia Institute of Technology, Atlanta, Georgia 30332, USA (Dated: 13 December
More informationSynchronization in delaycoupled bipartite networks
Synchronization in delaycoupled bipartite networks Ram Ramaswamy School of Physical Sciences Jawaharlal Nehru University, New Delhi February 20, 2015 Outline Ø Bipartite networks and delay-coupled phase
More informationAn analysis of how coupling parameters influence nonlinear oscillator synchronization
An analysis of how coupling parameters influence nonlinear oscillator synchronization Morris Huang, 1 Ben McInroe, 2 Mark Kingsbury, 2 and Will Wagstaff 3 1) School of Mechanical Engineering, Georgia Institute
More informationInfluence of Criticality on 1/f α Spectral Characteristics of Cortical Neuron Populations
Influence of Criticality on 1/f α Spectral Characteristics of Cortical Neuron Populations Robert Kozma rkozma@memphis.edu Computational Neurodynamics Laboratory, Department of Computer Science 373 Dunn
More informationSynchrony in Stochastic Pulse-coupled Neuronal Network Models
Synchrony in Stochastic Pulse-coupled Neuronal Network Models Katie Newhall Gregor Kovačič and Peter Kramer Aaditya Rangan and David Cai 2 Rensselaer Polytechnic Institute, Troy, New York 2 Courant Institute,
More informationFrom Baseline to Epileptiform Activity: A Path to Synchronized Rhythmicity in Large-Scale Neural Networks
From Baseline to Epileptiform Activity: A Path to Synchronized Rhythmicity in Large-Scale Neural Networks Vladimir Shusterman and William C. Troy Cardiovascular Institute and Department of Mathematics
More informationarxiv:nlin/ v1 [nlin.cd] 4 Oct 2005
Synchronization of Coupled Chaotic Dynamics on Networks R. E. Amritkar and Sarika Jalan Physical Research Laboratory, Navrangpura, Ahmedabad 380 009, India. arxiv:nlin/0510008v1 [nlin.cd] 4 Oct 2005 Abstract
More informationA Multivariate Time-Frequency Based Phase Synchrony Measure for Quantifying Functional Connectivity in the Brain
A Multivariate Time-Frequency Based Phase Synchrony Measure for Quantifying Functional Connectivity in the Brain Dr. Ali Yener Mutlu Department of Electrical and Electronics Engineering, Izmir Katip Celebi
More informationFast neural network simulations with population density methods
Fast neural network simulations with population density methods Duane Q. Nykamp a,1 Daniel Tranchina b,a,c,2 a Courant Institute of Mathematical Science b Department of Biology c Center for Neural Science
More informationNeural Networks 1 Synchronization in Spiking Neural Networks
CS 790R Seminar Modeling & Simulation Neural Networks 1 Synchronization in Spiking Neural Networks René Doursat Department of Computer Science & Engineering University of Nevada, Reno Spring 2006 Synchronization
More informationNETWORKS Lecture 2: Synchronization Fundamentals
Short-course: Complex Systems Beyond the Metaphor UNSW, February 2007 NETWORKS Lecture 2: Synchronization Fundamentals David J Hill Research School of Information Sciences and Engineering The ANU 8/2/2007
More informationkg meter ii) Note the dimensions of ρ τ are kg 2 velocity 2 meter = 1 sec 2 We will interpret this velocity in upcoming slides.
II. Generalizing the 1-dimensional wave equation First generalize the notation. i) "q" has meant transverse deflection of the string. Replace q Ψ, where Ψ may indicate other properties of the medium that
More informationSynchronization of Limit Cycle Oscillators by Telegraph Noise. arxiv: v1 [cond-mat.stat-mech] 5 Aug 2014
Synchronization of Limit Cycle Oscillators by Telegraph Noise Denis S. Goldobin arxiv:148.135v1 [cond-mat.stat-mech] 5 Aug 214 Department of Physics, University of Potsdam, Postfach 61553, D-14415 Potsdam,
More informationSynchrony in Neural Systems: a very brief, biased, basic view
Synchrony in Neural Systems: a very brief, biased, basic view Tim Lewis UC Davis NIMBIOS Workshop on Synchrony April 11, 2011 components of neuronal networks neurons synapses connectivity cell type - intrinsic
More informationPhase Locking. 1 of of 10. The PRC's amplitude determines which frequencies a neuron locks to. The PRC's slope determines if locking is stable
Printed from the Mathematica Help Browser 1 1 of 10 Phase Locking A neuron phase-locks to a periodic input it spikes at a fixed delay [Izhikevich07]. The PRC's amplitude determines which frequencies a
More informationUniversity of Bristol - Explore Bristol Research. Early version, also known as pre-print
Chizhov, Anto, V., Rodrigues, S., & Terry, JR. (6). A comparative analysis of an EEG model and a conductance-based. https://doi.org/1.116/j.physleta.7.4.6 Early version, also known as pre-print Link to
More informationHierarchical Clustering Identifies Hub Nodes in a Model of Resting-State Brain Activity
WCCI 22 IEEE World Congress on Computational Intelligence June, -5, 22 - Brisbane, Australia IJCNN Hierarchical Clustering Identifies Hub Nodes in a Model of Resting-State Brain Activity Mark Wildie and
More informationBifurcations and global stability of synchronized stationary states in the Kuramoto model for oscillator populations
PHYSICAL REVIEW E, VOLUME 64, 16218 Bifurcations and global stability of synchronized stationary states in the Kuramoto model for oscillator populations J. A. Acebrón, 1,2, * A. Perales, 3 and R. Spigler
More informationInternal and external synchronization of self-oscillating oscillators with non-identical control parameters
Internal and external synchronization of self-oscillating oscillators with non-identical control parameters Emelianova Yu.P., Kuznetsov A.P. Department of Nonlinear Processes, Saratov State University,
More informationFrom Baseline to Epileptiform Activity: A Path to Synchronized Rhythmicity in Large-Scale Neural Networks
From Baseline to Epileptiform Activity: A Path to Synchronized Rhythmicity in Large-Scale Neural Networks Vladimir Shusterman and William C. Troy Cardiovascular Institute and Department of Mathematics
More informationSelf-organized Criticality and Synchronization in a Pulse-coupled Integrate-and-Fire Neuron Model Based on Small World Networks
Commun. Theor. Phys. (Beijing, China) 43 (2005) pp. 466 470 c International Academic Publishers Vol. 43, No. 3, March 15, 2005 Self-organized Criticality and Synchronization in a Pulse-coupled Integrate-and-Fire
More informationNetworks: Lectures 9 & 10 Random graphs
Networks: Lectures 9 & 10 Random graphs Heather A Harrington Mathematical Institute University of Oxford HT 2017 What you re in for Week 1: Introduction and basic concepts Week 2: Small worlds Week 3:
More informationarxiv:cond-mat/ v1 [cond-mat.dis-nn] 6 Oct 2006
Revealing Network Connectivity From Dynamics Marc Timme arxiv:cond-mat/0610188v1 [cond-mat.dis-nn] 6 Oct 2006 Center for Applied Mathematics, Theoretical & Applied Mechanics, Kimball Hall, Cornell University,
More informationPrinciples of DCM. Will Penny. 26th May Principles of DCM. Will Penny. Introduction. Differential Equations. Bayesian Estimation.
26th May 2011 Dynamic Causal Modelling Dynamic Causal Modelling is a framework studying large scale brain connectivity by fitting differential equation models to brain imaging data. DCMs differ in their
More informationarxiv: v3 [physics.soc-ph] 12 Dec 2008
Synchronization in complex networks Alex Arenas Departament d Enginyeria Informàtica i Matemàtiques, Universitat Rovira i Virgili, 43007 Tarragona, Spain Institute for Biocomputation and Physics of Complex
More informationCommon noise vs Coupling in Oscillator Populations
Common noise vs Coupling in Oscillator Populations A. Pimenova, D. Goldobin, M. Rosenblum, and A. Pikovsky Institute of Continuous Media Mechanics UB RAS, Perm, Russia Institut for Physics and Astronomy,
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/319/5869/1543/dc1 Supporting Online Material for Synaptic Theory of Working Memory Gianluigi Mongillo, Omri Barak, Misha Tsodyks* *To whom correspondence should be addressed.
More informationUniversity of Bristol - Explore Bristol Research. Peer reviewed version. Link to published version (if available): /j.physleta
Rodrigues, S., Terry, JR., & Breakspear, M. (2006). On the genesis of spikewave oscillations in a mean-field model of human thalamic and corticothalamic dynamics. Physics letters A, 355, 352-357. DOI:
More informationNonlinear and Collective Effects in Mesoscopic Mechanical Oscillators
Dynamics Days Asia-Pacific: Singapore, 2004 1 Nonlinear and Collective Effects in Mesoscopic Mechanical Oscillators Alexander Zumdieck (Max Planck, Dresden), Ron Lifshitz (Tel Aviv), Jeff Rogers (HRL,
More informationSynchronization and Phase Oscillators
1 Synchronization and Phase Oscillators Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 Synchronization
More informationDynamics of delayed-coupled chaotic logistic maps: Influence of network topology, connectivity and delay times
PRAMANA c Indian Academy of Sciences Vol. 7, No. 6 journal of June 28 physics pp. 1 9 Dynamics of delayed-coupled chaotic logistic maps: Influence of network topology, connectivity and delay times ARTURO
More informationEEG- Signal Processing
Fatemeh Hadaeghi EEG- Signal Processing Lecture Notes for BSP, Chapter 5 Master Program Data Engineering 1 5 Introduction The complex patterns of neural activity, both in presence and absence of external
More informationPHASE-LOCKED SOLUTIONS IN A HUB CONNECTED OSCILLATOR RING NETWORK
Copyright c 29 by ABCM PHASE-LOCKED SOLUTIONS IN A HUB CONNECTED OSCILLATOR RING NETWORK Jacqueline Bridge, Jacqueline.Bridge@sta.uwi.edu Department of Mechanical Engineering, The University of the West
More informationSpiral Wave Chimeras
Spiral Wave Chimeras Carlo R. Laing IIMS, Massey University, Auckland, New Zealand Erik Martens MPI of Dynamics and Self-Organization, Göttingen, Germany Steve Strogatz Dept of Maths, Cornell University,
More informationCISC 3250 Systems Neuroscience
CISC 3250 Systems Neuroscience Systems Neuroscience How the nervous system performs computations How groups of neurons work together to achieve intelligence Professor Daniel Leeds dleeds@fordham.edu JMH
More informationEntrainment of coupled oscillators on regular networks by pacemakers
PHYSICAL REVIEW E 73, 036218 2006 Entrainment of coupled oscillators on regular networks by pacemakers Filippo Radicchi* and Hildegard Meyer-Ortmanns School of Engineering and Science, International University
More informationDynamical modelling of systems of coupled oscillators
Dynamical modelling of systems of coupled oscillators Mathematical Neuroscience Network Training Workshop Edinburgh Peter Ashwin University of Exeter 22nd March 2009 Peter Ashwin (University of Exeter)
More information7 Rate-Based Recurrent Networks of Threshold Neurons: Basis for Associative Memory
Physics 178/278 - David Kleinfeld - Fall 2005; Revised for Winter 2017 7 Rate-Based Recurrent etworks of Threshold eurons: Basis for Associative Memory 7.1 A recurrent network with threshold elements The
More informationEffects of Noisy Drive on Rhythms in Networks of Excitatory and Inhibitory Neurons
Effects of Noisy Drive on Rhythms in Networks of Excitatory and Inhibitory Neurons Christoph Börgers 1 and Nancy Kopell 2 1 Department of Mathematics, Tufts University, Medford, MA 2155 2 Department of
More informationCHAPTER 4 RELAXATION OSCILLATORS WITH TIME DELAY COUPLING
CHAPTER 4 RELAXATION OSCILLATORS WITH TIME DELAY COUPLING 4.1 Introduction Thus far we have studied synchronization in networks of locally coupled neurobiologically based oscillator models. One additional
More informationDynamical systems in neuroscience. Pacific Northwest Computational Neuroscience Connection October 1-2, 2010
Dynamical systems in neuroscience Pacific Northwest Computational Neuroscience Connection October 1-2, 2010 What do I mean by a dynamical system? Set of state variables Law that governs evolution of state
More informationSpectral Graph Theory for. Dynamic Processes on Networks
Spectral Graph Theory for Dynamic Processes on etworks Piet Van Mieghem in collaboration with Huijuan Wang, Dragan Stevanovic, Fernando Kuipers, Stojan Trajanovski, Dajie Liu, Cong Li, Javier Martin-Hernandez,
More informationON THE INVERSE KINEMATICS OF A FRAGMENT OF PROTEIN BACKBONE
ON THE INVERSE KINEMATIS OF A FRAGMENT OF PROTEIN BAKBONE Guanfeng Liu, R.J. Milgram, A. Dhanik, and J.. Latombe Department of omputer Science, Stanford University Department of Mathematics, Stanford University
More informationARTIFICIAL NEURAL NETWORK PART I HANIEH BORHANAZAD
ARTIFICIAL NEURAL NETWORK PART I HANIEH BORHANAZAD WHAT IS A NEURAL NETWORK? The simplest definition of a neural network, more properly referred to as an 'artificial' neural network (ANN), is provided
More information2 1. Introduction. Neuronal networks often exhibit a rich variety of oscillatory behavior. The dynamics of even a single cell may be quite complicated
GEOMETRIC ANALYSIS OF POPULATION RHYTHMS IN SYNAPTICALLY COUPLED NEURONAL NETWORKS J. Rubin and D. Terman Dept. of Mathematics; Ohio State University; Columbus, Ohio 43210 Abstract We develop geometric
More informationIDENTIFICATION OF TARGET IMAGE REGIONS BASED ON BIFURCATIONS OF A FIXED POINT IN A DISCRETE-TIME OSCILLATOR NETWORK
International Journal of Innovative Computing, Information and Control ICIC International c 203 ISSN 349-498 Volume 9, Number, January 203 pp. 355 363 IDENTIFICATION OF TARGET IMAGE REGIONS BASED ON BIFURCATIONS
More informationBursting Oscillations of Neurons and Synchronization
Bursting Oscillations of Neurons and Synchronization Milan Stork Applied Electronics and Telecommunications, Faculty of Electrical Engineering/RICE University of West Bohemia, CZ Univerzitni 8, 3064 Plzen
More informationProceedings of Neural, Parallel, and Scientific Computations 4 (2010) xx-xx PHASE OSCILLATOR NETWORK WITH PIECEWISE-LINEAR DYNAMICS
Proceedings of Neural, Parallel, and Scientific Computations 4 (2010) xx-xx PHASE OSCILLATOR NETWORK WITH PIECEWISE-LINEAR DYNAMICS WALTER GALL, YING ZHOU, AND JOSEPH SALISBURY Department of Mathematics
More informationHSND-2015, IPR. Department of Physics, University of Burdwan, Burdwan, West Bengal.
New kind of deaths: Oscillation Death and Chimera Death HSND-2015, IPR Dr. Tanmoy Banerjee Department of Physics, University of Burdwan, Burdwan, West Bengal. Points to be discussed Oscillation suppression
More informationCausality and communities in neural networks
Causality and communities in neural networks Leonardo Angelini, Daniele Marinazzo, Mario Pellicoro, Sebastiano Stramaglia TIRES-Center for Signal Detection and Processing - Università di Bari, Bari, Italy
More informationSpontaneous Synchronization in Complex Networks
B.P. Zeegers Spontaneous Synchronization in Complex Networks Bachelor s thesis Supervisors: dr. D. Garlaschelli (LION) prof. dr. W.Th.F. den Hollander (MI) August 2, 25 Leiden Institute of Physics (LION)
More informationCharacterising Probability Distributions via Entropies
1 Characterising Probability Distributions via Entropies Satyajit Thakor, Terence Chan and Alex Grant Indian Institute of Technology Mandi University of South Australia Myriota Pty Ltd arxiv:1602.03618v2
More informationStructure of Correlations in Neuronal Networks
Structure of Correlations in Neuronal Networks Krešimir Josić University of Houston James Trousdale (UH), Yu Hu (UW) Eric Shea-Brown (UW) The Connectome Van J. Wedeen, MGH/Harvard The Connectome Van J.
More informationWave Propagation in Uniaxial Media. Reflection and Transmission at Interfaces
Lecture 5: Crystal Optics Outline 1 Homogeneous, Anisotropic Media 2 Crystals 3 Plane Waves in Anisotropic Media 4 Wave Propagation in Uniaxial Media 5 Reflection and Transmission at Interfaces Christoph
More informationPhase Oscillators. and at r, Hence, the limit cycle at r = r is stable if and only if Λ (r ) < 0.
1 Phase Oscillators Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 2 Phase Oscillators Oscillations
More informationSpectral Clustering. Guokun Lai 2016/10
Spectral Clustering Guokun Lai 2016/10 1 / 37 Organization Graph Cut Fundamental Limitations of Spectral Clustering Ng 2002 paper (if we have time) 2 / 37 Notation We define a undirected weighted graph
More informationCHAPTER 3 SYNCHRONY IN RELAXATION OSCILLATORS
CHAPTER 3 SYNCHRONY IN RELAXATION OSCILLATORS 3.1 Introduction The phrase relaxation oscillations was coined by van der Pol in 1926 in his analysis of a triode circuit [van der Pol, 1926]. A relaxation
More informationAn Inverse Mass Expansion for Entanglement Entropy. Free Massive Scalar Field Theory
in Free Massive Scalar Field Theory NCSR Demokritos National Technical University of Athens based on arxiv:1711.02618 [hep-th] in collaboration with Dimitris Katsinis March 28 2018 Entanglement and Entanglement
More information7 Recurrent Networks of Threshold (Binary) Neurons: Basis for Associative Memory
Physics 178/278 - David Kleinfeld - Winter 2019 7 Recurrent etworks of Threshold (Binary) eurons: Basis for Associative Memory 7.1 The network The basic challenge in associative networks, also referred
More information1 Balanced networks: Trading speed for noise
Physics 178/278 - David leinfeld - Winter 2017 (Corrected yet incomplete notes) 1 Balanced networks: Trading speed for noise 1.1 Scaling of neuronal inputs An interesting observation is that the subthresold
More informationSTUDY OF SYNCHRONIZED MOTIONS IN A ONE-DIMENSIONAL ARRAY OF COUPLED CHAOTIC CIRCUITS
International Journal of Bifurcation and Chaos, Vol 9, No 11 (1999) 19 4 c World Scientific Publishing Company STUDY OF SYNCHRONIZED MOTIONS IN A ONE-DIMENSIONAL ARRAY OF COUPLED CHAOTIC CIRCUITS ZBIGNIEW
More informationSeveral ways to solve the MSO problem
Several ways to solve the MSO problem J. J. Steil - Bielefeld University - Neuroinformatics Group P.O.-Box 0 0 3, D-3350 Bielefeld - Germany Abstract. The so called MSO-problem, a simple superposition
More informationDesigning Information Devices and Systems II Fall 2015 Note 5
EE 16B Designing Information Devices and Systems II Fall 01 Note Lecture given by Babak Ayazifar (9/10) Notes by: Ankit Mathur Spectral Leakage Example Compute the length-8 and length-6 DFT for the following
More informationSynchronization in complex networks
Synchronization in complex networks Alex Arenas, 1, 2,3 Albert Díaz-Guilera, 4,2 Jurgen Kurths, 5 Yamir Moreno, 2,6 and Changsong Zhou 7 1 Departament d Enginyeria Informàtica i Matemàtiques, Universitat
More informationMATH 23 Exam 2 Review Solutions
MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution
More informationECE662: Pattern Recognition and Decision Making Processes: HW TWO
ECE662: Pattern Recognition and Decision Making Processes: HW TWO Purdue University Department of Electrical and Computer Engineering West Lafayette, INDIANA, USA Abstract. In this report experiments are
More informationDerivation of amplitude equations for nonlinear oscillators subject to arbitrary forcing
PHYSICAL REVIEW E 69, 066141 (2004) Derivation of amplitude equations for nonlinear oscillators subject to arbitrary forcing Catalina Mayol, Raúl Toral, and Claudio R. Mirasso Department de Física, Universitat
More informationCommon noise induces clustering in populations of globally coupled oscillators
OFFPRINT Common noise induces clustering in populations of globally coupled oscillators S. Gil, Y. Kuramoto and A. S. Mikhailov EPL, 88 (29) 65 Please visit the new website www.epljournal.org TARGET YOUR
More informationEE C245 ME C218 Introduction to MEMS Design
EE C45 ME C18 Introduction to MEMS Design Fall 008 Prof. Clark T.-C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 9470 Lecture 6: Output
More informationNORMAL MODES, WAVE MOTION AND THE WAVE EQUATION. Professor G.G.Ross. Oxford University Hilary Term 2009
NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION Professor G.G.Ross Oxford University Hilary Term 009 This course of twelve lectures covers material for the paper CP4: Differential Equations, Waves and
More informationThe Phase Response Curve of Reciprocally Inhibitory Model Neurons Exhibiting Anti-Phase Rhythms
The Phase Response Curve of Reciprocally Inhibitory Model Neurons Exhibiting Anti-Phase Rhythms Jiawei Zhang Timothy J. Lewis Department of Mathematics, University of California, Davis Davis, CA 9566,
More informationIntroduction to Neural Networks
Introduction to Neural Networks What are (Artificial) Neural Networks? Models of the brain and nervous system Highly parallel Process information much more like the brain than a serial computer Learning
More informationSynchronization in Phase-Coupled Kuramoto Oscillator Networks with Axonal Delay and Synaptic Plasticity
Dickinson College Dickinson Scholar Faculty and Staff Publications By Year Faculty and Staff Publications 3-214 Synchronization in Phase-Coupled Kuramoto Oscillator Networks with Axonal Delay and Synaptic
More informationarxiv: v1 [nlin.ao] 23 Sep 2015
One node driving synchronisation Chengwei Wang 1,*, Celso Grebogi 1, and Murilo S. Baptista 1 1 Institute for Comple Systems and Mathematical Biology, King s College, University of Aberdeen, Aberdeen,
More informationEffects of Interactive Function Forms in a Self-Organized Critical Model Based on Neural Networks
Commun. Theor. Phys. (Beijing, China) 40 (2003) pp. 607 613 c International Academic Publishers Vol. 40, No. 5, November 15, 2003 Effects of Interactive Function Forms in a Self-Organized Critical Model
More informationDigital Control: Summary # 7
Digital Control: Summary # 7 Proportional, integral and derivative control where K i is controller parameter (gain). It defines the ratio of the control change to the control error. Note that e(k) 0 u(k)
More informationYi Ming Lai. Rachel Nicks. Pedro Soares. Dynamical Systems with Noise from Renewal Processes. Synchrony Branching Lemma for Regular Networks
Yi Ming Lai Dynamical Systems with Noise from Renewal Processes One important area in the area of dynamics on networks is the study of temporal patterns of edges, for example, when edges appear or disappear
More informationLinear second-order differential equations with constant coefficients and nonzero right-hand side
Linear second-order differential equations with constant coefficients and nonzero right-hand side We return to the damped, driven simple harmonic oscillator d 2 y dy + 2b dt2 dt + ω2 0y = F sin ωt We note
More informationPatterns, Memory and Periodicity in Two-Neuron Delayed Recurrent Inhibitory Loops
Math. Model. Nat. Phenom. Vol. 5, No. 2, 2010, pp. 67-99 DOI: 10.1051/mmnp/20105203 Patterns, Memory and Periodicity in Two-Neuron Delayed Recurrent Inhibitory Loops J. Ma 1 and J. Wu 2 1 Department of
More informationarxiv: v1 [q-bio.nc] 16 May 2017
Optimized brute-force algorithms for the bifurcation analysis of a spin-glass-like neural network model Diego Fasoli 1,2,, Stefano Panzeri 1 arxiv:1705.05647v1 [q-bio.nc] 16 May 2017 1 Laboratory of Neural
More informationNeural Nets and Symbolic Reasoning Hopfield Networks
Neural Nets and Symbolic Reasoning Hopfield Networks Outline The idea of pattern completion The fast dynamics of Hopfield networks Learning with Hopfield networks Emerging properties of Hopfield networks
More informationDate: 1 April (1) The only reference material you may use is one 8½x11 crib sheet and a calculator.
PH1140: Oscillations and Waves Name: Solutions Conference: Date: 1 April 2005 EXAM #1: D2005 INSTRUCTIONS: (1) The only reference material you may use is one 8½x11 crib sheet and a calculator. (2) Show
More informationCollective and Stochastic Effects in Arrays of Submicron Oscillators
DYNAMICS DAYS: Long Beach, 2005 1 Collective and Stochastic Effects in Arrays of Submicron Oscillators Ron Lifshitz (Tel Aviv), Jeff Rogers (HRL, Malibu), Oleg Kogan (Caltech), Yaron Bromberg (Tel Aviv),
More informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of n-dof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:
More information5.3 METABOLIC NETWORKS 193. P (x i P a (x i )) (5.30) i=1
5.3 METABOLIC NETWORKS 193 5.3 Metabolic Networks 5.4 Bayesian Networks Let G = (V, E) be a directed acyclic graph. We assume that the vertices i V (1 i n) represent for example genes and correspond to
More informationData analysis methods in the neuroscience Zoltán Somogyvári
Data analysis methods in the neuroscience Zoltán Somogyvári Wigner Research Centre for Physics of the Hungarian Academy of Sciences Methods applicable to one time series The Fourier transformation The
More informationPhase synchronization of an ensemble of weakly coupled oscillators: A paradigm of sensor fusion
Phase synchronization of an ensemble of weakly coupled oscillators: A paradigm of sensor fusion Ariën J. van der Wal Netherlands Defence Academy (NLDA) 15 th ICCRTS Santa Monica, CA, June 22-24, 2010 Copyright
More informationDate: 31 March (1) The only reference material you may use is one 8½x11 crib sheet and a calculator.
PH1140: Oscillations and Waves Name: SOLUTIONS AT END Conference: Date: 31 March 2005 EXAM #1: D2006 INSTRUCTIONS: (1) The only reference material you may use is one 8½x11 crib sheet and a calculator.
More information1. Synchronization Phenomena
1. Synchronization Phenomena In nature synchronization happens all the time. In mechanical systems, in biological systems, in epidemiology, basically everywhere. When we talk about synchronization we usually
More information11/17/10. Chapter 14. Oscillations. Chapter 14. Oscillations Topics: Simple Harmonic Motion. Simple Harmonic Motion
11/17/10 Chapter 14. Oscillations This striking computergenerated image demonstrates an important type of motion: oscillatory motion. Examples of oscillatory motion include a car bouncing up and down,
More informationThe homogeneous Poisson process
The homogeneous Poisson process during very short time interval Δt there is a fixed probability of an event (spike) occurring independent of what happened previously if r is the rate of the Poisson process,
More informationResponse Selection Using Neural Phase Oscillators
Response Selection Using Neural Phase Oscillators J. Acacio de Barros 1 A symposium on the occasion of Patrick Suppes 90th birthday March 10, 2012 1 With Pat Suppes and Gary Oas (Suppes et al., 2012).
More information,, rectilinear,, spherical,, cylindrical. (6.1)
Lecture 6 Review of Vectors Physics in more than one dimension (See Chapter 3 in Boas, but we try to take a more general approach and in a slightly different order) Recall that in the previous two lectures
More information