Synchronization Transitions in Complex Networks
|
|
- Jade Morgan
- 5 years ago
- Views:
Transcription
1 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 Physics University of Zaragoza, Zaragoza 50009, Spain. 3 ISI Foundation, Turin, Italy Thematic School on Complex Networks Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 1 / 40
2 Outline A. Arenas, A. Diaz-Guilera,J.Kurths, Y. Moreno, and C. Zhou, Synchronization in Complex Networks, Physics Reports, 469, (2008). 1 Kuramoto model 2 Explosive Synchronization: Set up and Results 3 Master Stability Approach Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 2 / 40
3 Synchronization One of the most paradigmatic example of collective behavior Emergence of coherent rythms in the dynamics of coupled agents We will model these coupled agents as phase-oscillators... But first lets discuss a very useful approach. Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 3 / 40
4 The Kuramoto model We have a collection of N agents, being the dynamics of each of them described by a phase θ i [0, 2π): θ i = ω i + K N sin(θ j θ i ), (1) N j=1 Synchronization is measured via the Kuramoto order parameter: r = 1 N N e iθ j (2) j r 0 dynamical incoherence r 1 full synchronization. There exist a phase transition when increasing λ = K/N. At some critical coupling λ c = 2/[πg(0)] the incoherent solution becomes unstable and r starts to increase with λ. Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 4 / 40
5 The Synchronization transition Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 5 / 40
6 The Kuramoto model in Complex Networks We have a network of N agents, coupled as dictated by the adjacency matrix A. The Kuramoto dynamics of each agent is: N θ i = ω i + λ A ij sin(θ j θ i ), (3) j=1 Synchronization can be measured via the Kuramoto order parameter: r = 1 N N e iθ j, (4) and by the degree of synchronization between pairs of connected nodes: D ij = lim A 1 τ+t ij T e i[θ i (t) θ j (t)] dt. (5) T D ij 0 local incoherence τ j D ij = 1 local synchronization. Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 6 / 40
7 Synchronization transitions in Networks Two remarkable examples of network architecture: Scale-free networks: P(k) k γ Erdös-Rényi graphs: P(k) e k k k /k! Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 7 / 40
8 Synchronization transitions in Networks Synchronization diagram for SF Networks: evolution of the Kuramoto order parameter as a function of the coupling r(λ). The onset of synchronization occurs at a non-zero value of λ, i.e., λ c > 0 even when N R N=10 3 N=5x10 3 N=10 4 N=2x10 4 N=3x10 4 N=4x10 4 N=5x ! Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 8 / 40
9 r Synchronization transitions in Networks Finite Size Scaling Analysis The idea is that: if λ < λ c (subcritical regime), then r falls of as N 1/2. if λ > λ c (supercritical regime), then r const, though with O(N 1/2 ) fluctuations. at λ = λ c r falls as a power law. Our best estimate gives λ c = 0.05 ± 0.01 Besides, when both N, t we get r (λ λ c ) β, with β = !1 10!2 10!3!=0.02!=0.04!=0.06!= ! N Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 9 / 40
10 Synchronization transitions in Networks Stability of the Fully Synchronized state τ : Average time it takes for a node to be again in the synchronized state after being perturbed <!> 10 0 N=10 3, k max =91 τ k ν 10! k The more connected a node is, the more stable it is. Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 10 / 40
11 Synchronization transitions in Networks Stability of the Fully Synchronized state As we are perturbing a single node and this perturbation, ξ i, is small, we can consider that it only affects the first neighbors of the perturbed node. The stability analysis can be locally reduced to the problem of how such a perturbation relaxes in a star-like topology (perturbed node attached to k >> 1 nodes). Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 11 / 40
12 Synchronization transitions in Networks Stability of the Fully Synchronized state For a start graph we have: η i = 1 for i = 1,..., N 2 η N 1 = η hub = N = k hub 1 The fastest relaxation rate corresponds to the hub and goes like: 1 k hub for k hub 1 Finally, the superposition of many perturbations, each one corresponding to a k hub + 1 star, leads to different contributions of 1 k hub with k hub = k min,..., k max. Thus, τ k 1. Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 12 / 40
13 Synchronization transitions in Networks Synchronization diagrams in SF and ER networks: evolution of the Kuramoto order parameter as a function of the coupling r(λ) Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 13 / 40
14 Synchronization transitions in Networks Looking at the microscopic patterns with D ij : How nodes and links are incorporated into the Giant Synchronized Component? Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 14 / 40
15 Synchronization transitions in Networks Different paths towards synchronization: The importance of hubs in the entrainment of oscillators Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 15 / 40
16 Explosive Percolation In percolation processes the addition of a new rule for link creation turns the percolation transition explosive: D. Achlioptas, R.M. D Souza, and J. Spencer, Science 323, 1453 (2009). F. Radicchi, and S. Fortunato, Phys. Rev. Lett. 103, (2009). Y.S. Cho et al., Phys. Rev. Lett. 103, (2009). J. Nagler, A. Levina, and M. Timme, to appear in Nature Phyics (2011). Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 16 / 40
17 Explosive Synchronization Can we change the order of the synchronization transition from microscopic dynamics? Yes! Hubs promote synchronization: by integrating a mean-field dynamics they are able to collect many nodes in the synchronized component. We will destroy this topological ability dynamically: by assigning hubs a natural frequency far from the average one Ω. To this aim we propose this setting: This condition automatically implies: ω i = k i with i = 1,..., N (6) Ω = ω i = k i (7) g(ω) = P(k) (8) Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 17 / 40
18 Explosive Synchronization Going smoothly from ER graphs to SF networks: We confirm that for SF networks (i.e. only when hubs are present) the synchronization transition is explosive displaying a strong histeresis effect. Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 18 / 40
19 Explosive Synchronization We compute for each value of λ the effective frequency of each node: ω eff i = 1 T t+t t θ i (τ) dτ, with T 1. (9) Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 19 / 40
20 Explosive Synchronization What happens with a general model of SF networks (configurational ensembles with P(k) k γ ): The rule ω i = k i changes the order of the transition in SF networks Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 20 / 40
21 Explosive Synchronization In SF networks P(k) k γ so (when ω i = k i and thus P(k) = g(ω)) we have g(ω) ω γ. Is this broad frequency distribution behind the first-order transition? No. We do need the correlation structure-dynamics to obtain the explosive synchronization! Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 21 / 40
22 The star graph The star graph is a first proxy of a SF network capturing the neiborhood of a hub: One central hub connecting with K leaves. We set the natural frequency of the leaves identical with value ω The hub naturally pace at frequency ω h. The average frequency is Ω = (Kω + ω h )/(K + 1) Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 22 / 40
23 The star graph The equations of motion for the hub and the leaves read: θ h = ω h + λ K sin(θ j θ h ), (10) j=1 θ j = ω + λ sin(θ h θ j ), with j = 1,..., K. (11) Hub motion: We settle in a rotating frame with frequency Ω so that the hub motion is described as: φ h = (ω h Ω) + λ(k + 1)r sin(φ h ), (12) Imposing that the phase of the hub is locked, φ h = 0, we obtain: sin φ h = (ω h Ω) λ(k + 1)r. (13) Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 23 / 40
24 The star graph Leaves motion: Again, in a rotating frame with frequency Ω so that the motion of leaf is described as: φ j = (ω Ω) + λ sin(φ h φ j ), with j = 1,..., K. (14) Imposing that the phase of the hub is locked, φ h = 0, we obtain: [1 cos φ j = (Ω ω) sin φ h ± sin 2 ] φ h [λ 2 (Ω ω) 2 ]. (15) λ This equation implies that locking is lost at λ c = (Ω ω). The order parameter is r = cosφ. Thus, at λ c, r takes the value: r c = 1 K (Ω ω) (ω h Ω) λ c(k+1)r c + 1 (ω h Ω) 2 K + 1 λ c λ 2 c(k + 1) 2 rc 2. (16) Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 24 / 40
25 The star graph In our case (ω i = k i ): ω = 1, ω h = K, Ω = 2K/(K + 1) thus: λ c = K 1 K + 1 and r c = K K + 1. (17) Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 25 / 40
26 Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 26 / 40
27 Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 27 / 40
28 Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 28 / 40
29 Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 29 / 40
30 Master Stability Function formalism Derivation follows the paper Dynamical Systems on Networks: a tutorial a highly recommendable review by MAP & JPG, arxiv preprint Let s suppose that each node i is associated with a single variable x i. We use x to denote the vector of variables. Consider the continuous dynamical system ẋ i = f i (x i ) + N A ij g ij (x i, x j ), i {1,..., N}, (18) j=1 A ij = adjacency matrix g ij (x i, x j ) effects coupling to neighbors. the equilibrium points Eq. (18) satisfy x i = 0 for all nodes i. local stability of these points via linear stability analysis: let x i = xi (where ɛ i 1) and take a Taylor expansion. + ɛ i Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 30 / 40
31 For each i, we obtain Master Stability Function formalism N ẋ i = ɛ i = f i (x i + ɛ i ) + A ij g ij (x i + ɛ i, x j + ɛ j ), j=1 N = f i (x i ) + A ij g ij (x i, x j ) + ɛ i f N i xi =x + ɛ i A ij j=1 i j=1 g ij x i xi =x i,x j =x j N + ɛ j A ij j=1 g ij x j xi =x i,x j =x j +... = α 1 + α 2 + α 3 + α 4 + α , (19) where f i := df i dx i and the α l terms are defined in order. Because x is an equilibrium, it follows that α 1 + α 2 = 0. The terms α 3 and α 4 are linear in ɛ i, and α 5 is linear in ɛ j. We now neglect all higher-order terms. Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 31 / 40
32 Master Stability Function formalism To simplify notation, we define a i := ɛ i f i b ij := g ij x i c ij := g ij x j xi =x i, xi =x i,x j =x j xi =x i,x j =x j,. (20) That is, ɛ = Mɛ +..., (21) where M = [M ij ] and [ M ij = δ ij a i + ] b ik A ik + c ij A ij. (22) k Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 32 / 40
33 Master Stability Function formalism Assume that M has N unique eigenvectors (not always the case, although it will usually be if away from a bifurcation point) and is thus diagonalizable, we obtain ɛ = N α r (t)v r, (23) r=1 where v r (with corresponding eigenvalue µ r ) is the rth (right) eigenvector of the matrix M. It follows that ɛ = N r=1 N α r v r = M α r (t)v r = r=1 r=1 N µ r α r (t)v r. (24) Separately equating the linearly independent terms in equation (24) then yields α r = µ r α r, which in turn implies that α r (t) = α r (0) exp(µ r t). Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 33 / 40
34 Master Stability Function formalism As usual for dynamical systems, we obtain local asymptotic stability if Re(µ r ) < 0 for all r, instability if any Re(µ r ) > 0, and a marginal stability (for which one needs to examine nonlinear terms) if any Re(µ r ) = 0 for some r and none of the eigenvalues have a positive real part. Example, let s consider a (significantly) simplified situation in which every node has the same equilibrium location: i.e., xi = x for all nodes i. (This arises, for example, in the SI model of a biological contagion.) We will also assume that f i f for all nodes and g ij g for all node pairs. Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 34 / 40
35 Master Stability Function formalism With these simplifications, it follows that f (x ) + N A ij g(x, x ) = f (x ) + k i g(x, x ) = 0, (25) j=1 where we recall that k i is the degree of node i. Equation (25) implies that either all nodes have the same degree (i.e., that our graph z-regular ) or that g(x, x ) = 0. We do not wish to restrict the network structure severely, so we will suppose that the latter condition holds. It follows that f (x ) = 0, so the equilibria of the coupled equation (18) in this simplified situation are necessarily the same as the equilibria of the intrinsic dynamics that are satisfied by individual (uncoupled) nodes. This yields a simplified version of the notation from equation (20): a i a := f x=x, b ij b := g, x i xi =x j =x c ij c := g. x j xi =x j =x (26) Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 35 / 40
36 Master Stability Function formalism We thus obtain ɛ i = (a + bk i )ɛ i + c N A ij ɛ j, i {1,..., N}. (27) j=1 If we assume that g(x i, x j ) = g(x j ), which is yet another major simplifying assumption, we obtain Consequently, b = 0 and ẋ i = f (x i ) + where I is the N N identity matrix. N A ij g(x j ). (28) j=1 ɛ = (ai + ca)ɛ, (29) Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 36 / 40
37 Master Stability Function formalism An equilibrium of (29) is (locally) asymptotically stable if and only if all of the eigenvalues of P := ai + ca = P T are negative. (The matrix P is symmetric, so all of its eigenvalues are guaranteed to be real.) Letting w r denote an eigenvector of A with corresponding eigenvalue λ r. It follows that (ai + ca)w r = (a + cλ r )w r (30) for all r (where there are at most N eigenvectors and there will be exactly N of them if we are able to diagonalize A), so w r is also an eigenvector of the matrix P, with corresponding eigenvalue (a + cλ r ). For (local) asymptotic stability, we thus need a + cλ r < 0 to hold for all λ r. This, in turn, implies that we need a < 0 because the adjacency matrix A is guaranteed to have both positive and negative eigenvalues. Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 37 / 40
38 Master Stability Function formalism We thus need: (i) λ r < a/c for c > 0 and (ii) λ r > a/c for c < 0. If (i) is satisfied for the most positive eigenvalue λ 1 of A, then it (obviously) must be satisfied for all eigenvalues of A. If (ii) is satisfied for the most negative eigenvalue λ N of A, then it (obviously) must be satisfied for all eigenvalues of A. It follows that 1 λ N < c a < 1 λ 1, (31) which is more insightful when we insert the definitions of a and c. This yields g xi 1 x j =x < j =x λ N f x=x < 1 λ 1. (32) Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 38 / 40
39 Master Stability Function formalism The left and right terms in equation (32), which is called a master stability condition, depend only on the structure of the network, and the central term depends only on the nature (i.e., functional forms of the individual dynamics and of the coupling terms) of the dynamics. It also illustrates that the eigenvalues of adjacency matrices have important ramifications for dynamical behavior when studying dynamical systems on networks. Indeed, investigations of the spectra (i.e., set of eigenvales) of adjacency matrices can yield crucial insights about dynamical systems on networks. Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 39 / 40
40 Master Stability Function formalism N ẋ = F (x i ) + σ L ij H(x j ), j=1 i = 1,..., N where L ij is the Laplacian matrix: L ij = 1 if i and j are connected, L ii = k i, and 0 otherwise. It can be shown that the linear stability of the synchronized state is determined by the eigenvalues of L, effectively decoupling the dynamics and the topology. In particular, if λ 1 = 0 λ 2... λ N, are the eigenvalues of L, then it follows that: The larger the ratio λ N /λ 2 is, the more stable the synchronized state is and vice versa. Y. Moreno (BIFI) Synchronization in Complex Networks Les Houches, 7-18/04/14 40 / 40
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 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 informationSolutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)
Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddle-node bifurcation and show that the first order system dx dt = r x e x
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 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 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 informationMATH6142 Complex Networks Exam 2016
MATH642 Complex Networks Exam 26 Solution Problem a) The network is directed because the adjacency matrix is not symmetric. b) The network is shown in Figure. (4 marks)[unseen] 4 3 Figure : The directed
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
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 informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationDynamical systems tutorial. Gregor Schöner, INI, RUB
Dynamical systems tutorial Gregor Schöner, INI, RUB Dynamical systems: Tutorial the word dynamics time-varying measures range of a quantity forces causing/accounting for movement => dynamical systems dynamical
More informationarxiv: v1 [cond-mat.stat-mech] 2 Apr 2013
Link-disorder fluctuation effects on synchronization in random networks Hyunsuk Hong, Jaegon Um, 2 and Hyunggyu Park 2 Department of Physics and Research Institute of Physics and Chemistry, Chonbuk National
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 information8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)
8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)
More informationFINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.
FINAL EXAM MATH0 Theory of Ordinary Differential Equations There are 5 problems on 2 pages. Spring 2009. 25 points Consider the linear plane autonomous system x + y x y. Find a fundamental matrix of the
More information5.2.2 Planar Andronov-Hopf bifurcation
138 CHAPTER 5. LOCAL BIFURCATION THEORY 5.. Planar Andronov-Hopf bifurcation What happens if a planar system has an equilibrium x = x 0 at some parameter value α = α 0 with eigenvalues λ 1, = ±iω 0, ω
More informationSolutions 2: Simple Harmonic Oscillator and General Oscillations
Massachusetts Institute of Technology MITES 2017 Physics III Solutions 2: Simple Harmonic Oscillator and General Oscillations Due Wednesday June 21, at 9AM under Rene García s door Preface: This problem
More informationPhase Desynchronization as a Mechanism for Transitions to High-Dimensional Chaos
Commun. Theor. Phys. (Beijing, China) 35 (2001) pp. 682 688 c International Academic Publishers Vol. 35, No. 6, June 15, 2001 Phase Desynchronization as a Mechanism for Transitions to High-Dimensional
More informationMAT 22B - Lecture Notes
MAT 22B - Lecture Notes 4 September 205 Solving Systems of ODE Last time we talked a bit about how systems of ODE arise and why they are nice for visualization. Now we'll talk about the basics of how to
More informationLong-wave Instability in Anisotropic Double-Diffusion
Long-wave Instability in Anisotropic Double-Diffusion Jean-Luc Thiffeault Institute for Fusion Studies and Department of Physics University of Texas at Austin and Neil J. Balmforth Department of Theoretical
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationEmergence of synchronization in complex networks of interacting dynamical systems
Physica D 224 (2006) 114 122 www.elsevier.com/locate/physd Emergence of synchronization in complex networks of interacting dynamical systems Juan G. Restrepo a,,1, Edward Ott a,b, Brian R. Hunt c a Institute
More informationLINEAR RESPONSE THEORY
MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior
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 informationRecent progress on the classical and quantum Kuramoto synchronization
Recent progress on the classical and quantum Kuramoto synchronization Seung Yeal Ha Department of Mathematical Sciences Seoul National University March 23rd, 2017 Outline PROLOGUE What is synchronization?
More informationExplosive percolation in graphs
Home Search Collections Journals About Contact us My IOPscience Explosive percolation in graphs This article has been downloaded from IOPscience. Please scroll down to see the full text article. 2011 J.
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
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 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 informationModels Involving Interactions between Predator and Prey Populations
Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predator-prey models are used to show the intricate
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
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 information= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :
1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change
More informationME 680- Spring Geometrical Analysis of 1-D Dynamical Systems
ME 680- Spring 2014 Geometrical Analysis of 1-D Dynamical Systems 1 Geometrical Analysis of 1-D Dynamical Systems Logistic equation: n = rn(1 n) velocity function Equilibria or fied points : initial conditions
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for COMPUTATIONAL BIOLOGY A. COURSE CODES: FFR 110, FIM740GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for COMPUTATIONAL BIOLOGY A COURSE CODES: FFR 110, FIM740GU, PhD Time: Place: Teachers: Allowed material: Not allowed: June 8, 2018, at 08 30 12 30 Johanneberg Kristian
More informationStructural Consensus Controllability of Singular Multi-agent Linear Dynamic Systems
Structural Consensus Controllability of Singular Multi-agent Linear Dynamic Systems M. ISAL GARCÍA-PLANAS Universitat Politècnica de Catalunya Departament de Matèmatiques Minería 1, sc. C, 1-3, 08038 arcelona
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More information7 Two-dimensional bifurcations
7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed
More informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty
More informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More informationSynchronization between different motifs
Synchronization between different motifs Li Ying( ) a) and Liu Zeng-Rong( ) b) a) College of Information Technology, Shanghai Ocean University, Shanghai 201306, China b) Institute of Systems Biology, Shanghai
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More informationEpidemics in Complex Networks and Phase Transitions
Master M2 Sciences de la Matière ENS de Lyon 2015-2016 Phase Transitions and Critical Phenomena Epidemics in Complex Networks and Phase Transitions Jordan Cambe January 13, 2016 Abstract Spreading phenomena
More informationAlgebraic Representation of Networks
Algebraic Representation of Networks 0 1 2 1 1 0 0 1 2 0 0 1 1 1 1 1 Hiroki Sayama sayama@binghamton.edu Describing networks with matrices (1) Adjacency matrix A matrix with rows and columns labeled by
More information35. RESISTIVE INSTABILITIES: CLOSING REMARKS
35. RESISTIVE INSTABILITIES: CLOSING REMARKS In Section 34 we presented a heuristic discussion of the tearing mode. The tearing mode is centered about x = 0, where F x that B 0 ( ) = k! B = 0. Note that
More informationDispersion relations, stability and linearization
Dispersion relations, stability and linearization 1 Dispersion relations Suppose that u(x, t) is a function with domain { < x 0}, and it satisfies a linear, constant coefficient partial differential
More informationData Mining and Analysis: Fundamental Concepts and Algorithms
Data Mining and Analysis: Fundamental Concepts and Algorithms dataminingbook.info Mohammed J. Zaki 1 Wagner Meira Jr. 2 1 Department of Computer Science Rensselaer Polytechnic Institute, Troy, NY, USA
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 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 informationPhysics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension
Physics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension In these notes we examine Bloch s theorem and band structure in problems with periodic potentials, as a part of our survey
More informationSection 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d
Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d July 6, 2009 Today s Session Today s Session A Summary of This Session: Today s Session A Summary of This
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
More informationPhase Synchronization of Coupled Rossler Oscillators: Amplitude Effect
Commun. Theor. Phys. (Beijing, China) 47 (2007) pp. 265 269 c International Academic Publishers Vol. 47, No. 2, February 15, 2007 Phase Synchronization of Coupled Rossler Oscillators: Amplitude Effect
More informationECS 289 F / MAE 298, Lecture 15 May 20, Diffusion, Cascades and Influence
ECS 289 F / MAE 298, Lecture 15 May 20, 2014 Diffusion, Cascades and Influence Diffusion and cascades in networks (Nodes in one of two states) Viruses (human and computer) contact processes epidemic thresholds
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 informationGeneralized correlated states in a ring of coupled nonlinear oscillators with a local injection
PHYSICAL REVIEW E 66, 6621 22 Generalized correlated states in a ring of coupled nonlinear oscillators with a local injection Y. Chembo ouomou 1,2 and P. Woafo 1, * 1 Laboratoire de Mécanique, Faculté
More informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More information6.207/14.15: Networks Lectures 4, 5 & 6: Linear Dynamics, Markov Chains, Centralities
6.207/14.15: Networks Lectures 4, 5 & 6: Linear Dynamics, Markov Chains, Centralities 1 Outline Outline Dynamical systems. Linear and Non-linear. Convergence. Linear algebra and Lyapunov functions. Markov
More informationKrein-Rutman Theorem and the Principal Eigenvalue
Chapter 1 Krein-Rutman Theorem and the Principal Eigenvalue The Krein-Rutman theorem plays a very important role in nonlinear partial differential equations, as it provides the abstract basis for the proof
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 informationExplosive Percolation in Erdős-Rényi-Like Random Graph Processes
Explosive Percolation in Erdős-Rényi-Like Random Graph Processes Konstantinos Panagiotou a,1 Reto Spöhel a,1,4 Angelika Steger b,2 Henning Thomas b,2,3 a Max Planck Institute for Informatics, 66123 Saarbrücken,
More informationPractice Final Exam Solutions for Calculus II, Math 1502, December 5, 2013
Practice Final Exam Solutions for Calculus II, Math 5, December 5, 3 Name: Section: Name of TA: This test is to be taken without calculators and notes of any sorts. The allowed time is hours and 5 minutes.
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationIn these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.
1 2 Linear Systems In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation 21 Matrix ODEs Let and is a scalar A linear function satisfies Linear superposition ) Linear
More informationWaves and characteristics: Overview 5-1
Waves and characteristics: Overview 5-1 Chapter 5: Waves and characteristics Overview Physics and accounting: use example of sound waves to illustrate method of linearization and counting of variables
More informationSeminar 6: COUPLED HARMONIC OSCILLATORS
Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached
More information7 Pendulum. Part II: More complicated situations
MATH 35, by T. Lakoba, University of Vermont 60 7 Pendulum. Part II: More complicated situations In this Lecture, we will pursue two main goals. First, we will take a glimpse at a method of Classical Mechanics
More informationRobot Control Basics CS 685
Robot Control Basics CS 685 Control basics Use some concepts from control theory to understand and learn how to control robots Control Theory general field studies control and understanding of behavior
More informationThe Hamiltonian Mean Field Model: Effect of Network Structure on Synchronization Dynamics. Yogesh Virkar University of Colorado, Boulder.
The Hamiltonian Mean Field Model: Effect of Network Structure on Synchronization Dynamics Yogesh Virkar University of Colorado, Boulder. 1 Collaborators Prof. Juan G. Restrepo Prof. James D. Meiss Department
More informationStochastic models, patterns formation and diffusion
Stochastic models, patterns formation and diffusion Duccio Fanelli Francesca Di Patti, Tommaso Biancalani Dipartimento di Energetica, Università degli Studi di Firenze CSDC Centro Interdipartimentale per
More informationPattern formation and Turing instability
Pattern formation and Turing instability. Gurarie Topics: - Pattern formation through symmetry breaing and loss of stability - Activator-inhibitor systems with diffusion Turing proposed a mechanism for
More informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More informationOn the topology of synchrony optimized networks of a Kuramoto-model with non-identical oscillators
On the topology of synchrony optimized networks of a Kuramoto-model with non-identical oscillators David Kelly and Georg A. Gottwald School of Mathematics and Statistics University of Sydney NSW 26, Australia
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationRandom Networks. Complex Networks, CSYS/MATH 303, Spring, Prof. Peter Dodds
Complex Networks, CSYS/MATH 303, Spring, 2010 Prof. Peter Dodds Department of Mathematics & Statistics Center for Complex Systems Vermont Advanced Computing Center University of Vermont Licensed under
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationGeometric Series and the Ratio and Root Test
Geometric Series and the Ratio and Root Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 5, 2018 Outline 1 Geometric Series
More information16 Period doubling route to chaos
16 Period doubling route to chaos We now study the routes or scenarios towards chaos. We ask: How does the transition from periodic to strange attractor occur? The question is analogous to the study of
More informationAverage-Consensus of Multi-Agent Systems with Direct Topology Based on Event-Triggered Control
Outline Background Preliminaries Consensus Numerical simulations Conclusions Average-Consensus of Multi-Agent Systems with Direct Topology Based on Event-Triggered Control Email: lzhx@nankai.edu.cn, chenzq@nankai.edu.cn
More informationConsensus Stabilizability and Exact Consensus Controllability of Multi-agent Linear Systems
Consensus Stabilizability and Exact Consensus Controllability of Multi-agent Linear Systems M. ISABEL GARCÍA-PLANAS Universitat Politècnica de Catalunya Departament de Matèmatiques Minería 1, Esc. C, 1-3,
More information8.385 MIT (Rosales) Hopf Bifurcations. 2 Contents Hopf bifurcation for second order scalar equations. 3. Reduction of general phase plane case to seco
8.385 MIT Hopf Bifurcations. Rodolfo R. Rosales Department of Mathematics Massachusetts Institute of Technology Cambridge, Massachusetts MA 239 September 25, 999 Abstract In two dimensions a Hopf bifurcation
More informationExact Consensus Controllability of Multi-agent Linear Systems
Exact Consensus Controllability of Multi-agent Linear Systems M. ISAEL GARCÍA-PLANAS Universitat Politècnica de Catalunya Departament de Matèmatiques Minería 1, Esc. C, 1-3, 08038 arcelona SPAIN maria.isabel.garcia@upc.edu
More informationNonlinear Control Lecture 2:Phase Plane Analysis
Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 r. Farzaneh Abdollahi Nonlinear Control Lecture 2 1/53
More informationThe Kuramoto Model. Gerald Cooray. U.U.D.M. Project Report 2008:23. Department of Mathematics Uppsala University
U.U.D.M. Project Report 008:3 The Kuramoto Model Gerald Cooray Examensarbete i matematik, 30 hp Handledare och examinator: David Sumpter September 008 Department of Mathematics Uppsala University THE
More informationSummary of topics relevant for the final. p. 1
Summary of topics relevant for the final p. 1 Outline Scalar difference equations General theory of ODEs Linear ODEs Linear maps Analysis near fixed points (linearization) Bifurcations How to analyze a
More informationarxiv: v1 [nlin.ao] 8 Dec 2017
Nonlinearity 31(1):R1-R23 (2018) Identical synchronization of nonidentical oscillators: when only birds of different feathers flock together arxiv:1712.03245v1 [nlin.ao] 8 Dec 2017 Yuanzhao Zhang 1 and
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 informationPhonons and lattice dynamics
Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed
More informationVirgili, Tarragona (Spain) Roma (Italy) Zaragoza, Zaragoza (Spain)
Int.J.Complex Systems in Science vol. 1 (2011), pp. 47 54 Probabilistic framework for epidemic spreading in complex networks Sergio Gómez 1,, Alex Arenas 1, Javier Borge-Holthoefer 1, Sandro Meloni 2,3
More informationBifurcations of mutually coupled equations in random graphs
Bifurcations of mutually coupled equations in random graphs Eduardo Garibaldi UNICAMP Department of Mathematics 13083-859 Campinas SP, Brazil garibaldi@ime.unicamp.br Tiago Pereira Department of Mathematics
More informationLecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems
University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 45 Lecture Notes 6: Dynamic Equations Part C: Linear Difference Equation Systems Peter J. Hammond latest revision 2017 September
More informationRandom Networks. Complex Networks CSYS/MATH 303, Spring, Prof. Peter Dodds
Complex Networks CSYS/MATH 303, Spring, 2011 Prof. Peter Dodds Department of Mathematics & Statistics Center for Complex Systems Vermont Advanced Computing Center University of Vermont Licensed under the
More information8 Ecosystem stability
8 Ecosystem stability References: May [47], Strogatz [48]. In these lectures we consider models of populations, with an emphasis on the conditions for stability and instability. 8.1 Dynamics of a single
More informationDYNAMICS OF THREE COUPLED VAN DER POL OSCILLATORS WITH APPLICATION TO CIRCADIAN RHYTHMS
Proceedings of IDETC/CIE 2005 ASME 2005 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference September 24-28, 2005, Long Beach, California USA DETC2005-84017
More informationON THE CRITICAL COUPLING FOR KURAMOTO OSCILLATORS
ON THE CRITICAL COUPLING FOR KURAMOTO OSCILLATORS FLORIAN DÖRFLER AND FRANCESCO BULLO Abstract. The celebrated Kuramoto model captures various synchronization phenomena in biological and man-made dynamical
More informationarxiv: v2 [nlin.cd] 8 Nov 2017
Non-identical multiplexing promotes chimera states Saptarshi Ghosh a, Anna Zakharova b, Sarika Jalan a a Complex Systems Lab, Discipline of Physics, Indian Institute of Technology Indore, Simrol, Indore
More informationLotka Volterra Predator-Prey Model with a Predating Scavenger
Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
More information