1
|
Chauhan K, Neiman AB, Tass PA. Synaptic reorganization of synchronized neuronal networks with synaptic weight and structural plasticity. PLoS Comput Biol 2024; 20:e1012261. [PMID: 38980898 DOI: 10.1371/journal.pcbi.1012261] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2023] [Accepted: 06/20/2024] [Indexed: 07/11/2024] Open
Abstract
Abnormally strong neural synchronization may impair brain function, as observed in several brain disorders. We computationally study how neuronal dynamics, synaptic weights, and network structure co-emerge, in particular, during (de)synchronization processes and how they are affected by external perturbation. To investigate the impact of different types of plasticity mechanisms, we combine a network of excitatory integrate-and-fire neurons with different synaptic weight and/or structural plasticity mechanisms: (i) only spike-timing-dependent plasticity (STDP), (ii) only homeostatic structural plasticity (hSP), i.e., without weight-dependent pruning and without STDP, (iii) a combination of STDP and hSP, i.e., without weight-dependent pruning, and (iv) a combination of STDP and structural plasticity (SP) that includes hSP and weight-dependent pruning. To accommodate the diverse time scales of neuronal firing, STDP, and SP, we introduce a simple stochastic SP model, enabling detailed numerical analyses. With tools from network theory, we reveal that structural reorganization may remarkably enhance the network's level of synchrony. When weaker contacts are preferentially eliminated by weight-dependent pruning, synchrony is achieved with significantly sparser connections than in randomly structured networks in the STDP-only model. In particular, the strengthening of contacts from neurons with higher natural firing rates to those with lower rates and the weakening of contacts in the opposite direction, followed by selective removal of weak contacts, allows for strong synchrony with fewer connections. This activity-led network reorganization results in the emergence of degree-frequency, degree-degree correlations, and a mixture of degree assortativity. We compare the stimulation-induced desynchronization of synchronized states in the STDP-only model (i) with the desynchronization of models (iii) and (iv). The latter require stimuli of significantly higher intensity to achieve long-term desynchronization. These findings may inform future pre-clinical and clinical studies with invasive or non-invasive stimulus modalities aiming at inducing long-lasting relief of symptoms, e.g., in Parkinson's disease.
Collapse
Affiliation(s)
- Kanishk Chauhan
- Department of Physics and Astronomy, Ohio University, Athens, Ohio, United States of America
- Neuroscience Program, Ohio University, Athens, Ohio, United States of America
| | - Alexander B Neiman
- Department of Physics and Astronomy, Ohio University, Athens, Ohio, United States of America
- Neuroscience Program, Ohio University, Athens, Ohio, United States of America
| | - Peter A Tass
- Department of Neurosurgery, Stanford University, Stanford, California, United States of America
| |
Collapse
|
2
|
Bolotov MI, Munyayev VO, Smirnov LA, Osipov GV, Belykh I. Breathing and switching cyclops states in Kuramoto networks with higher-mode coupling. Phys Rev E 2024; 109:054202. [PMID: 38907462 DOI: 10.1103/physreve.109.054202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2024] [Accepted: 04/10/2024] [Indexed: 06/24/2024]
Abstract
Cyclops states are intriguing cluster patterns observed in oscillator networks, including neuronal ensembles. The concept of cyclops states formed by two distinct, coherent clusters and a solitary oscillator was introduced by Munyaev et al. [Phys. Rev. Lett. 130, 107201 (2023)0031-900710.1103/PhysRevLett.130.107201], where we explored the surprising prevalence of such states in repulsive Kuramoto networks of rotators with higher-mode harmonics in the coupling. This paper extends our analysis to understand the mechanisms responsible for destroying the cyclops' states and inducing dynamical patterns called breathing and switching cyclops states. We first analytically study the existence and stability of cyclops states in the Kuramoto-Sakaguchi networks of two-dimensional oscillators with inertia as a function of the second coupling harmonic. We then describe two bifurcation scenarios that give birth to breathing and switching cyclops states. We demonstrate that these states and their hybrids are prevalent across a wide coupling range and are robust against a relatively large intrinsic frequency detuning. Beyond the Kuramoto networks, breathing and switching cyclops states promise to strongly manifest in other physical and biological networks, including coupled theta neurons.
Collapse
Affiliation(s)
- Maxim I Bolotov
- Department of Control Theory, Lobachevsky State University of Nizhny Novgorod, 23 Gagarin Avenue, Nizhny Novgorod, 603022, Russia
| | - Vyacheslav O Munyayev
- Department of Control Theory, Lobachevsky State University of Nizhny Novgorod, 23 Gagarin Avenue, Nizhny Novgorod, 603022, Russia
| | - Lev A Smirnov
- Department of Control Theory, Lobachevsky State University of Nizhny Novgorod, 23 Gagarin Avenue, Nizhny Novgorod, 603022, Russia
| | - Grigory V Osipov
- Department of Control Theory, Lobachevsky State University of Nizhny Novgorod, 23 Gagarin Avenue, Nizhny Novgorod, 603022, Russia
| | - Igor Belykh
- Department of Mathematics and Statistics and Neuroscience Institute, Georgia State University, P.O. Box 4110, Atlanta, Georgia 30302-410, USA
| |
Collapse
|
3
|
Papo D, Buldú JM. Does the brain behave like a (complex) network? I. Dynamics. Phys Life Rev 2024; 48:47-98. [PMID: 38145591 DOI: 10.1016/j.plrev.2023.12.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2023] [Accepted: 12/10/2023] [Indexed: 12/27/2023]
Abstract
Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.
Collapse
Affiliation(s)
- D Papo
- Department of Neuroscience and Rehabilitation, Section of Physiology, University of Ferrara, Ferrara, Italy; Center for Translational Neurophysiology, Fondazione Istituto Italiano di Tecnologia, Ferrara, Italy.
| | - J M Buldú
- Complex Systems Group & G.I.S.C., Universidad Rey Juan Carlos, Madrid, Spain
| |
Collapse
|
4
|
Bröhl T, Rings T, Pukropski J, von Wrede R, Lehnertz K. The time-evolving epileptic brain network: concepts, definitions, accomplishments, perspectives. FRONTIERS IN NETWORK PHYSIOLOGY 2024; 3:1338864. [PMID: 38293249 PMCID: PMC10825060 DOI: 10.3389/fnetp.2023.1338864] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/15/2023] [Accepted: 12/19/2023] [Indexed: 02/01/2024]
Abstract
Epilepsy is now considered a network disease that affects the brain across multiple levels of spatial and temporal scales. The paradigm shift from an epileptic focus-a discrete cortical area from which seizures originate-to a widespread epileptic network-spanning lobes and hemispheres-considerably advanced our understanding of epilepsy and continues to influence both research and clinical treatment of this multi-faceted high-impact neurological disorder. The epileptic network, however, is not static but evolves in time which requires novel approaches for an in-depth characterization. In this review, we discuss conceptual basics of network theory and critically examine state-of-the-art recording techniques and analysis tools used to assess and characterize a time-evolving human epileptic brain network. We give an account on current shortcomings and highlight potential developments towards an improved clinical management of epilepsy.
Collapse
Affiliation(s)
- Timo Bröhl
- Department of Epileptology, University of Bonn Medical Centre, Bonn, Germany
- Helmholtz Institute for Radiation and Nuclear Physics, University of Bonn, Bonn, Germany
| | - Thorsten Rings
- Department of Epileptology, University of Bonn Medical Centre, Bonn, Germany
- Helmholtz Institute for Radiation and Nuclear Physics, University of Bonn, Bonn, Germany
| | - Jan Pukropski
- Department of Epileptology, University of Bonn Medical Centre, Bonn, Germany
| | - Randi von Wrede
- Department of Epileptology, University of Bonn Medical Centre, Bonn, Germany
| | - Klaus Lehnertz
- Department of Epileptology, University of Bonn Medical Centre, Bonn, Germany
- Helmholtz Institute for Radiation and Nuclear Physics, University of Bonn, Bonn, Germany
- Interdisciplinary Center for Complex Systems, University of Bonn, Bonn, Germany
| |
Collapse
|
5
|
Lizier JT, Bauer F, Atay FM, Jost J. Analytic relationship of relative synchronizability to network structure and motifs. Proc Natl Acad Sci U S A 2023; 120:e2303332120. [PMID: 37669393 PMCID: PMC10500263 DOI: 10.1073/pnas.2303332120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Accepted: 08/07/2023] [Indexed: 09/07/2023] Open
Abstract
Synchronization phenomena on networks have attracted much attention in studies of neural, social, economic, and biological systems, yet we still lack a systematic understanding of how relative synchronizability relates to underlying network structure. Indeed, this question is of central importance to the key theme of how dynamics on networks relate to their structure more generally. We present an analytic technique to directly measure the relative synchronizability of noise-driven time-series processes on networks, in terms of the directed network structure. We consider both discrete-time autoregressive processes and continuous-time Ornstein-Uhlenbeck dynamics on networks, which can represent linearizations of nonlinear systems. Our technique builds on computation of the network covariance matrix in the space orthogonal to the synchronized state, enabling it to be more general than previous work in not requiring either symmetric (undirected) or diagonalizable connectivity matrices and allowing arbitrary self-link weights. More importantly, our approach quantifies the relative synchronization specifically in terms of the contribution of process motif (walk) structures. We demonstrate that in general the relative abundance of process motifs with convergent directed walks (including feedback and feedforward loops) hinders synchronizability. We also reveal subtle differences between the motifs involved for discrete or continuous-time dynamics. Our insights analytically explain several known general results regarding synchronizability of networks, including that small-world and regular networks are less synchronizable than random networks.
Collapse
Affiliation(s)
- Joseph T. Lizier
- School of Computer Science and Centre for Complex Systems, Faculty of Engineering, The University of Sydney, Sydney, NSW2006, Australia
- Max Planck Institute for Mathematics in the Sciences, Leipzig04103, Germany
| | - Frank Bauer
- Max Planck Institute for Mathematics in the Sciences, Leipzig04103, Germany
- Department of Mathematics, Harvard University, Cambridge, MA02138
| | - Fatihcan M. Atay
- Max Planck Institute for Mathematics in the Sciences, Leipzig04103, Germany
- Department of Mathematics, Bilkent University, Ankara06800, Turkey
| | - Jürgen Jost
- Max Planck Institute for Mathematics in the Sciences, Leipzig04103, Germany
- Santa Fe Institute, Santa Fe, NM87501
| |
Collapse
|
6
|
Ruiz-Herrera A. The role of the spatial topology in trophic metacommunities: Species with reduced mobility and total population size. J Theor Biol 2023; 566:111479. [PMID: 37075827 DOI: 10.1016/j.jtbi.2023.111479] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2022] [Revised: 02/10/2023] [Accepted: 03/24/2023] [Indexed: 04/21/2023]
Abstract
A central question in ecology is understanding the influence of the spatial topology on the dynamics of a metacommunity. This is not an easy task, as most fragmented ecosystems have trophic interactions involving many species and patches. Recent attempts to solve this challenge have introduced certain simplifying assumptions or focused on a limited set of examples. These simplifications make the models mathematically tractable but keep away from real-world problems. In this paper, we provide a novel methodology to describe the influence of the spatial topology on the total population size of the species when the dispersal rates are small. The main conclusion is that the influence of the spatial topology is the result of the influence of each path in isolation. Here, a path refers to a pairwise connection between two patches. Our framework can be readily used with any metacommunity, and therefore represents a unification of biological insights. We also discuss several applications regarding the construction of ecological corridors.
Collapse
|
7
|
Munyayev VO, Bolotov MI, Smirnov LA, Osipov GV, Belykh I. Cyclops States in Repulsive Kuramoto Networks: The Role of Higher-Order Coupling. PHYSICAL REVIEW LETTERS 2023; 130:107201. [PMID: 36962033 DOI: 10.1103/physrevlett.130.107201] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Accepted: 02/09/2023] [Indexed: 06/18/2023]
Abstract
Repulsive oscillator networks can exhibit multiple cooperative rhythms, including chimera and cluster splay states. Yet, understanding which rhythm prevails remains challenging. Here, we address this fundamental question in the context of Kuramoto-Sakaguchi networks of rotators with higher-order Fourier modes in the coupling. Through analysis and numerics, we show that three-cluster splay states with two distinct coherent clusters and a solitary oscillator are the prevalent rhythms in networks with an odd number of units. We denote such tripod patterns cyclops states with the solitary oscillator reminiscent of the Cyclops' eye. As their mythological counterparts, the cyclops states are giants that dominate the system's phase space in weakly repulsive networks with first-order coupling. Astonishingly, the addition of the second or third harmonics to the Kuramoto coupling function makes the cyclops states global attractors practically across the full range of coupling's repulsion. Beyond the Kuramoto oscillators, we show that this effect is robustly present in networks of canonical theta neurons with adaptive coupling. At a more general level, our results suggest clues for finding dominant rhythms in repulsive physical and biological networks.
Collapse
Affiliation(s)
- Vyacheslav O Munyayev
- Department of Control Theory, Lobachevsky State University of Nizhny Novgorod, 23 Gagarin Avenue, Nizhny Novgorod 603022, Russia
| | - Maxim I Bolotov
- Department of Control Theory, Lobachevsky State University of Nizhny Novgorod, 23 Gagarin Avenue, Nizhny Novgorod 603022, Russia
| | - Lev A Smirnov
- Department of Control Theory, Lobachevsky State University of Nizhny Novgorod, 23 Gagarin Avenue, Nizhny Novgorod 603022, Russia
| | - Grigory V Osipov
- Department of Control Theory, Lobachevsky State University of Nizhny Novgorod, 23 Gagarin Avenue, Nizhny Novgorod 603022, Russia
| | - Igor Belykh
- Department of Mathematics and Statistics and Neuroscience Institute, Georgia State University, P.O. Box 4110, Atlanta, Georgia 30302-410, USA
| |
Collapse
|
8
|
Khramenkov VA, Dmitrichev AS, Nekorkin VI. A new scenario for Braess's paradox in power grids. CHAOS (WOODBURY, N.Y.) 2022; 32:113116. [PMID: 36456330 DOI: 10.1063/5.0093980] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Accepted: 09/30/2022] [Indexed: 06/17/2023]
Abstract
We consider several topologies of power grids and analyze how the addition of transmission lines affects their dynamics. The main example we are dealing with is a power grid that has a tree-like three-element motif at the periphery. We establish conditions where the addition of a transmission line in the motif enhances its stability or induces Braess's paradox and reduces stability of the entire grid. By using bifurcation theory and nonlocal stability analysis, we show that two scenarios for Braess's paradox are realized in the grid. The first scenario is well described and is associated with the disappearance of the synchronous mode. The second scenario has not been previously described and is associated with the reduction of nonlocal stability of the synchronous mode due to the appearance of asynchronous modes. The necessary conditions for stable operation of the grid, under the addition of a line, are derived. It is proved that the new scenario for Braess's paradox is realized in the grids with more complex topologies even when several lines are added in their bulks.
Collapse
Affiliation(s)
- V A Khramenkov
- Department of Nonlinear Dynamics, Institute of Applied Physics of RAS, 46 Ulyanov Str., 603950 Nizhny Novgorod, Russia
| | - A S Dmitrichev
- Department of Nonlinear Dynamics, Institute of Applied Physics of RAS, 46 Ulyanov Str., 603950 Nizhny Novgorod, Russia
| | - V I Nekorkin
- Department of Nonlinear Dynamics, Institute of Applied Physics of RAS, 46 Ulyanov Str., 603950 Nizhny Novgorod, Russia
| |
Collapse
|
9
|
Functional control of oscillator networks. Nat Commun 2022; 13:4721. [PMID: 35953467 PMCID: PMC9372149 DOI: 10.1038/s41467-022-31733-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Accepted: 06/28/2022] [Indexed: 11/11/2022] Open
Abstract
Oscillatory activity is ubiquitous in natural and engineered network systems. The interaction scheme underlying interdependent oscillatory components governs the emergence of network-wide patterns of synchrony that regulate and enable complex functions. Yet, understanding, and ultimately harnessing, the structure-function relationship in oscillator networks remains an outstanding challenge of modern science. Here, we address this challenge by presenting a principled method to prescribe exact and robust functional configurations from local network interactions through optimal tuning of the oscillators’ parameters. To quantify the behavioral synchrony between coupled oscillators, we introduce the notion of functional pattern, which encodes the pairwise relationships between the oscillators’ phases. Our procedure is computationally efficient and provably correct, accounts for constrained interaction types, and allows to concurrently assign multiple desired functional patterns. Further, we derive algebraic and graph-theoretic conditions to guarantee the feasibility and stability of target functional patterns. These conditions provide an interpretable mapping between the structural constraints and their functional implications in oscillator networks. As a proof of concept, we apply the proposed method to replicate empirically recorded functional relationships from cortical oscillations in a human brain, and to redistribute the active power flow in different models of electrical grids. In network systems governed by oscillatory activity, such as brain networks or power grids, configurations of synchrony may define network functions. The authors introduce a control approach for the formation of desired synchrony patterns through optimal interventions on the network parameters.
Collapse
|
10
|
Tang L, Smith K, Daley K, Belykh I. When multilayer links exchange their roles in synchronization. Phys Rev E 2022; 106:024214. [PMID: 36109922 DOI: 10.1103/physreve.106.024214] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Accepted: 08/09/2022] [Indexed: 06/15/2023]
Abstract
Real world networks contain multiple layers of links whose interactions can lead to extraordinary collective dynamics, including synchronization. The fundamental problem of assessing how network topology controls synchronization in multilayer networks remains open due to serious limitations of the existing stability methods. Towards removing this obstacle, we propose an approximation method which significantly enhances the predictive power of the master stability function for stable synchronization in multilayer networks. For a class of saddle-focus oscillators, including Rössler and piecewise linear systems, our method reduces the complex stability analysis to simply solving a set of linear algebraic equations. Using the method, we analytically predict surprising effects due to multilayer coupling. In particular, we prove that two coupling layers-one of which would alone hamper synchronization and the other would foster it-reverse their roles when used in a multilayer network. We also analytically demonstrate that increasing the size of a globally coupled layer, that in isolation would induce stable synchronization, makes the multilayer network unsynchronizable.
Collapse
Affiliation(s)
- Longkun Tang
- School of Mathematical Sciences, Huaqiao University, Quanzhou 362021, China
- Department of Mathematics and Statistics, Georgia State University, P.O. Box 4110, Atlanta, Georgia 30302-410, USA
| | - Kelley Smith
- Department of Mathematics and Statistics, Georgia State University, P.O. Box 4110, Atlanta, Georgia 30302-410, USA
| | - Kevin Daley
- Department of Mathematics and Statistics, Georgia State University, P.O. Box 4110, Atlanta, Georgia 30302-410, USA
| | - Igor Belykh
- Department of Mathematics and Statistics, Georgia State University, P.O. Box 4110, Atlanta, Georgia 30302-410, USA
| |
Collapse
|
11
|
Wang B, Zhang S, Yu X, Niu Y, Niu J, Li D, Zhang S, Xiang J, Yan T, Yang J, Wu J, Liu M. Alterations in white matter network dynamics in patients with schizophrenia and bipolar disorder. Hum Brain Mapp 2022; 43:3909-3922. [PMID: 35567336 PMCID: PMC9374889 DOI: 10.1002/hbm.25892] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2022] [Revised: 03/17/2022] [Accepted: 04/08/2022] [Indexed: 11/12/2022] Open
Abstract
Emerging evidence suggests white matter network abnormalities in patients with schizophrenia (SZ) and bipolar disorder (BD), but the alterations in dynamics of the white matter network in patients with SZ and BD are largely unknown. The white matter network of patients with SZ (n = 45) and BD (n = 47) and that of healthy controls (HC, n = 105) were constructed. We used dynamics network control theory to quantify the dynamics metrics of the network, including controllability and synchronizability, to measure the ability to transfer between different states. Experiments show that the patients with SZ and BD showed decreasing modal controllability and synchronizability and increasing average controllability. The correlations between the average controllability and synchronizability of patients were broken, especially for those with SZ. The patients also showed alterations in brain regions with supercontroller roles and their distribution in the cognitive system. Finally, we were able to accurately discriminate and predict patients with SZ and BD. Our findings provide novel dynamic metrics evidence that patients with SZ and BD are characterized by a selective disruption of brain network controllability, potentially leading to reduced brain state transfer capacity, and offer new guidance for the clinical diagnosis of mental illness.
Collapse
Affiliation(s)
- Bin Wang
- Department of Information and Computer, Taiyuan University of Technology, Taiyuan, China
| | - Shanshan Zhang
- Department of Information and Computer, Taiyuan University of Technology, Taiyuan, China
| | - Xuexue Yu
- Department of Information and Computer, Taiyuan University of Technology, Taiyuan, China
| | - Yan Niu
- Department of Information and Computer, Taiyuan University of Technology, Taiyuan, China
| | - Jinliang Niu
- Department of Medical Imaging, The Second Hospital of Shanxi Medical University, Taiyuan, China
| | - Dandan Li
- Department of Information and Computer, Taiyuan University of Technology, Taiyuan, China
| | - Shan Zhang
- Department of Information and Computer, Taiyuan University of Technology, Taiyuan, China
| | - Jie Xiang
- Department of Information and Computer, Taiyuan University of Technology, Taiyuan, China
| | - Ting Yan
- Teranslational Medicine Research Center, Shanxi Medical University, Taiyuan, China
| | - Jiajia Yang
- Graduate School of Interdisciplinary Science and Engineering in Health Systems, Okayama, Japan
| | - Jinglong Wu
- Graduate School of Interdisciplinary Science and Engineering in Health Systems, Okayama, Japan
| | - Miaomiao Liu
- School of Psychology, Shenzhen University, Shenzhen, China
| |
Collapse
|
12
|
Barabash NV, Belykh VN, Osipov GV, Belykh IV. Partial synchronization in the second-order Kuramoto model: An auxiliary system method. CHAOS (WOODBURY, N.Y.) 2021; 31:113113. [PMID: 34881584 DOI: 10.1063/5.0066663] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/11/2021] [Accepted: 10/11/2021] [Indexed: 06/13/2023]
Abstract
Partial synchronization emerges in an oscillator network when the network splits into clusters of coherent and incoherent oscillators. Here, we analyze the stability of partial synchronization in the second-order finite-dimensional Kuramoto model of heterogeneous oscillators with inertia. Toward this goal, we develop an auxiliary system method that is based on the analysis of a two-dimensional piecewise-smooth system whose trajectories govern oscillating dynamics of phase differences between oscillators in the coherent cluster. Through a qualitative bifurcation analysis of the auxiliary system, we derive explicit bounds that relate the maximum natural frequency mismatch, inertia, and the network size that can support stable partial synchronization. In particular, we predict threshold-like stability loss of partial synchronization caused by increasing inertia. Our auxiliary system method is potentially applicable to cluster synchronization with multiple coherent clusters and more complex network topology.
Collapse
Affiliation(s)
- Nikita V Barabash
- Department of Mathematics, Volga State University of Water Transport, 5A, Nesterov str., Nizhny Novgorod 603950, Russia
| | - Vladimir N Belykh
- Department of Mathematics, Volga State University of Water Transport, 5A, Nesterov str., Nizhny Novgorod 603950, Russia
| | - Grigory V Osipov
- Department of Control Theory, Lobachevsky State University of Nizhny Novgorod, 23, Gagarin Ave., 603950 Nizhny Novgorod, Russia
| | - Igor V Belykh
- Department of Control Theory, Lobachevsky State University of Nizhny Novgorod, 23, Gagarin Ave., 603950 Nizhny Novgorod, Russia
| |
Collapse
|
13
|
Designing temporal networks that synchronize under resource constraints. Nat Commun 2021; 12:3273. [PMID: 34075037 PMCID: PMC8169648 DOI: 10.1038/s41467-021-23446-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Accepted: 04/22/2021] [Indexed: 11/10/2022] Open
Abstract
Being fundamentally a non-equilibrium process, synchronization comes with unavoidable energy costs and has to be maintained under the constraint of limited resources. Such resource constraints are often reflected as a finite coupling budget available in a network to facilitate interaction and communication. Here, we show that introducing temporal variation in the network structure can lead to efficient synchronization even when stable synchrony is impossible in any static network under the given budget, thereby demonstrating a fundamental advantage of temporal networks. The temporal networks generated by our open-loop design are versatile in the sense of promoting synchronization for systems with vastly different dynamics, including periodic and chaotic dynamics in both discrete-time and continuous-time models. Furthermore, we link the dynamic stabilization effect of the changing topology to the curvature of the master stability function, which provides analytical insights into synchronization on temporal networks in general. In particular, our results shed light on the effect of network switching rate and explain why certain temporal networks synchronize only for intermediate switching rate. The ability of complex networks to synchronize themselves is limited by available coupling resources. Zhang and Strogatz show that allowing temporal variation in the network structure can lead to synchronization even when stable synchrony is impossible in any static network under the fixed budget.
Collapse
|
14
|
Sugitani Y, Zhang Y, Motter AE. Synchronizing Chaos with Imperfections. PHYSICAL REVIEW LETTERS 2021; 126:164101. [PMID: 33961469 DOI: 10.1103/physrevlett.126.164101] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Revised: 01/17/2021] [Accepted: 02/26/2021] [Indexed: 06/12/2023]
Abstract
Previous research on nonlinear oscillator networks has shown that chaos synchronization is attainable for identical oscillators but deteriorates in the presence of parameter mismatches. Here, we identify regimes for which the opposite occurs and show that oscillator heterogeneity can synchronize chaos for conditions under which identical oscillators cannot. This effect is not limited to small mismatches and is observed for random oscillator heterogeneity on both homogeneous and heterogeneous network structures. The results are demonstrated experimentally using networks of Chua's oscillators and are further supported by numerical simulations and theoretical analysis. In particular, we propose a general mechanism based on heterogeneity-induced mode mixing that provides insights into the observed phenomenon. Since individual differences are ubiquitous and often unavoidable in real systems, it follows that such imperfections can be an unexpected source of synchronization stability.
Collapse
Affiliation(s)
- Yoshiki Sugitani
- Department of Electrical and Electronic System Engineering, Ibaraki University, 4-12-1 Nakanarusawa, Hitachi, Ibaraki 316-8511, Japan
| | - Yuanzhao Zhang
- Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208, USA
- Center for Applied Mathematics, Cornell University, Ithaca, New York 14853, USA
| | - Adilson E Motter
- Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208, USA
- Northwestern Institute on Complex Systems, Northwestern University, Evanston, Illinois 60208, USA
| |
Collapse
|
15
|
Bou Assi E, Zerouali Y, Robert M, Lesage F, Pouliot P, Nguyen DK. Large-Scale Desynchronization During Interictal Epileptic Discharges Recorded With Intracranial EEG. Front Neurol 2020; 11:529460. [PMID: 33424733 PMCID: PMC7785800 DOI: 10.3389/fneur.2020.529460] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2020] [Accepted: 11/27/2020] [Indexed: 11/13/2022] Open
Abstract
It is increasingly recognized that deep understanding of epileptic seizures requires both localizing and characterizing the functional network of the region where they are initiated, i. e., the epileptic focus. Previous investigations of the epileptogenic focus' functional connectivity have yielded contrasting results, reporting both pathological increases and decreases during resting periods and seizures. In this study, we shifted paradigm to investigate the time course of connectivity in relation to interictal epileptiform discharges. We recruited 35 epileptic patients undergoing intracranial EEG (iEEG) investigation as part of their presurgical evaluation. For each patient, 50 interictal epileptic discharges (IEDs) were marked and iEEG signals were epoched around those markers. Signals were narrow-band filtered and time resolved phase-locking values were computed to track the dynamics of functional connectivity during IEDs. Results show that IEDs are associated with a transient decrease in global functional connectivity, time-locked to the peak of the discharge and specific to the high range of the gamma frequency band. Disruption of the long-range connectivity between the epileptic focus and other brain areas might be an important process for the generation of epileptic activity. Transient desynchronization could be a potential biomarker of the epileptogenic focus since 1) the functional connectivity involving the focus decreases significantly more than the connectivity outside the focus and 2) patients with good surgical outcome appear to have a significantly more disconnected focus than patients with bad outcomes.
Collapse
Affiliation(s)
- Elie Bou Assi
- University of Montreal Hospital Research Center (CRCHUM), University of Montreal, Montreal, QC, Canada.,Department of Neuroscience, University of Montreal, Montreal, QC, Canada
| | - Younes Zerouali
- University of Montreal Hospital Research Center (CRCHUM), University of Montreal, Montreal, QC, Canada.,Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada
| | - Manon Robert
- University of Montreal Hospital Research Center (CRCHUM), University of Montreal, Montreal, QC, Canada
| | - Frederic Lesage
- Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada.,Montreal Heart Institute, Montreal, QC, Canada
| | - Philippe Pouliot
- Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, QC, Canada.,Montreal Heart Institute, Montreal, QC, Canada
| | - Dang K Nguyen
- University of Montreal Hospital Research Center (CRCHUM), University of Montreal, Montreal, QC, Canada.,Department of Neuroscience, University of Montreal, Montreal, QC, Canada
| |
Collapse
|
16
|
|
17
|
Daley K, Zhao K, Belykh IV. Synchronizability of directed networks: The power of non-existent ties. CHAOS (WOODBURY, N.Y.) 2020; 30:043102. [PMID: 32357666 DOI: 10.1063/1.5134920] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/03/2019] [Accepted: 03/16/2020] [Indexed: 06/11/2023]
Abstract
The understanding of how synchronization in directed networks is influenced by structural changes in network topology is far from complete. While the addition of an edge always promotes synchronization in a wide class of undirected networks, this addition may impede synchronization in directed networks. In this paper, we develop the augmented graph stability method, which allows for explicitly connecting the stability of synchronization to changes in network topology. The transformation of a directed network into a symmetrized-and-augmented undirected network is the central component of this new method. This transformation is executed by symmetrizing and weighting the underlying connection graph and adding new undirected edges with consideration made for the mean degree imbalance of each pair of nodes. These new edges represent "non-existent ties" in the original directed network and often control the location of critical nodes whose directed connections can be altered to manipulate the stability of synchronization in a desired way. In particular, we show that the addition of small-world shortcuts to directed networks, which makes "non-existent ties" disappear, can worsen the synchronizability, thereby revealing a destructive role of small-world connections in directed networks. An extension of our method may open the door to studying synchronization in directed multilayer networks, which cannot be effectively handled by the eigenvalue-based methods.
Collapse
Affiliation(s)
- Kevin Daley
- Department of Mathematics and Statistics, Georgia State University, 30 Pryor Street, Atlanta, Georgia 30303, USA
| | - Kun Zhao
- Department of Mathematics and Statistics, Georgia State University, 30 Pryor Street, Atlanta, Georgia 30303, USA
| | - Igor V Belykh
- Department of Mathematics and Statistics, Georgia State University, 30 Pryor Street, Atlanta, Georgia 30303, USA
| |
Collapse
|
18
|
Topology Effects on Sparse Control of Complex Networks with Laplacian Dynamics. Sci Rep 2019; 9:9034. [PMID: 31227756 PMCID: PMC6588614 DOI: 10.1038/s41598-019-45476-6] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2018] [Accepted: 06/10/2019] [Indexed: 12/23/2022] Open
Abstract
Ease of control of complex networks has been assessed extensively in terms of structural controllability and observability, and minimum control energy criteria. Here we adopt a sparsity-promoting feedback control framework for undirected networks with Laplacian dynamics and distinct topological features. The control objective considered is to minimize the effect of disturbance signals, magnitude of control signals and cost of feedback channels. We show that depending on the cost of feedback channels, different complex network structures become the least expensive option to control. Specifically, increased cost of feedback channels favors organized topological complexity such as modularity and centralization. Thus, although sparse and heterogeneous undirected networks may require larger numbers of actuators and sensors for structural controllability, networks with Laplacian dynamics are shown to be easier to control when accounting for the cost of feedback channels.
Collapse
|
19
|
He Z, Yao C, Yu J, Zhan M. Perturbation analysis and comparison of network synchronization methods. Phys Rev E 2019; 99:052207. [PMID: 31212531 DOI: 10.1103/physreve.99.052207] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2018] [Indexed: 11/07/2022]
Abstract
In many networked systems, synchronization is important and useful, and how to enhance synchronizability is an interesting problem. Based on the matrix perturbation theory, we analyze five methods of network synchronization enhancement, including the link removal, node removal, dividing hub node, pull control, and pinning control methods, and obtain explicit expressions for eigenvalue changes. By these comparisons, we find that, among all these methods, the pull control method is remarkable, as it can extend the synchronization (coupling strength) region from both the left and right sides, for any controlled node. Extensive simulation results are given to support the accuracy of the perturbation-based analysis.
Collapse
Affiliation(s)
- Zhiwei He
- Department of Mathematics, Shaoxing University, Shaoxing 312000, China
| | - Chenggui Yao
- Department of Mathematics, Shaoxing University, Shaoxing 312000, China
| | - Jun Yu
- Institute of Nonlinear Science, Shaoxing University, Shaoxing 312000, China
| | - Meng Zhan
- State Key Laboratory of Advanced Electromagnetic Engineering and Technology, School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
| |
Collapse
|
20
|
Papo D, Buldú JM. Brain synchronizability, a false friend. Neuroimage 2019; 196:195-199. [PMID: 30986500 DOI: 10.1016/j.neuroimage.2019.04.029] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2018] [Revised: 03/28/2019] [Accepted: 04/08/2019] [Indexed: 01/20/2023] Open
Abstract
Synchronization plays a fundamental role in healthy cognitive and motor function. However, how synchronization depends on the interplay between local dynamics, coupling and topology and how prone to synchronization a network is, given its topological organization, are still poorly understood issues. To investigate the synchronizability of both anatomical and functional brain networks various studies resorted to the Master Stability Function (MSF) formalism, an elegant tool which allows analysing the stability of synchronous states in a dynamical system consisting of many coupled oscillators. Here, we argue that brain dynamics does not fulfil the formal criteria under which synchronizability is usually quantified and, perhaps more importantly, this measure refers to a global dynamical condition that never holds in the brain (not even in the most pathological conditions), and therefore no neurophysiological conclusions should be drawn based on it. We discuss the meaning of synchronizability and its applicability to neuroscience and propose alternative ways to quantify brain networks synchronization.
Collapse
Affiliation(s)
- D Papo
- SCALab UMR CNRS 9193, Université de Lille, Villeneuve d'Ascq, France.
| | - J M Buldú
- Laboratory of Biological Networks, Center for Biomedical Technology (UPM), 28223, Pozuelo de Alarcón, Madrid, Spain; Complex Systems Group & G.I.S.C., Universidad Rey Juan Carlos, 28933, Móstoles, Madrid, Spain
| |
Collapse
|
21
|
Hart JD, Zhang Y, Roy R, Motter AE. Topological Control of Synchronization Patterns: Trading Symmetry for Stability. PHYSICAL REVIEW LETTERS 2019; 122:058301. [PMID: 30822003 DOI: 10.1103/physrevlett.122.058301] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/18/2018] [Revised: 01/08/2019] [Indexed: 06/09/2023]
Abstract
Symmetries are ubiquitous in network systems and have profound impacts on the observable dynamics. At the most fundamental level, many synchronization patterns are induced by underlying network symmetry, and a high degree of symmetry is believed to enhance the stability of identical synchronization. Yet, here we show that the synchronizability of almost any symmetry cluster in a network of identical nodes can be enhanced precisely by breaking its structural symmetry. This counterintuitive effect holds for generic node dynamics and arbitrary network structure and is, moreover, robust against noise and imperfections typical of real systems, which we demonstrate by implementing a state-of-the-art optoelectronic experiment. These results lead to new possibilities for the topological control of synchronization patterns, which we substantiate by presenting an algorithm that optimizes the structure of individual clusters under various constraints.
Collapse
Affiliation(s)
- Joseph D Hart
- Institute for Research in Electronics and Applied Physics, University of Maryland, College Park, Maryland 20742, USA
- Department of Physics, University of Maryland, College Park, Maryland 20742, USA
| | - Yuanzhao Zhang
- Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208, USA
| | - Rajarshi Roy
- Institute for Research in Electronics and Applied Physics, University of Maryland, College Park, Maryland 20742, USA
- Department of Physics, University of Maryland, College Park, Maryland 20742, USA
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
| | - Adilson E Motter
- Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208, USA
- Northwestern Institute on Complex Systems, Northwestern University, Evanston, Illinois 60208, USA
| |
Collapse
|
22
|
Jeter R, Porfiri M, Belykh I. Overcoming network resilience to synchronization through non-fast stochastic broadcasting. CHAOS (WOODBURY, N.Y.) 2018; 28:071104. [PMID: 30070517 DOI: 10.1063/1.5044420] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/12/2018] [Accepted: 07/06/2018] [Indexed: 06/08/2023]
Abstract
Stochastic broadcasting is an important and understudied paradigm for controlling networks. In this paper, we examine the feasibility of on-off broadcasting from a single reference node to induce synchronization in a target network with connections from the reference node that stochastically switch in time with an arbitrary switching period. Internal connections within the target network are static and promote the network's resilience to externally induced synchronization. Through rigorous mathematical analysis, we uncover a complex interplay between the network topology and the switching period of stochastic broadcasting, fostering or hindering synchronization to the reference node. We derive a criterion which reveals an explicit dependence of induced synchronization on the properties of the network (the Laplacian spectrum) and the switching process (strength of broadcasting, switching period, and switching probabilities). With coupled chaotic tent maps as our test-bed, we prove the emergence of "windows of opportunity" where only non-fast switching periods are favorable to synchronization. The size of these windows of opportunity is shaped by the Laplacian spectrum such that the switching period needs to be manipulated accordingly to induce synchronization. Surprisingly, only the zero and the largest eigenvalues of the Laplacian matrix control these windows of opportunities for tent maps within a wide parameter region.
Collapse
Affiliation(s)
- Russell Jeter
- Department of Mathematics and Statistics, Georgia State University, P.O. Box 4110, Atlanta, Georgia 30302-4110, USA
| | - Maurizio Porfiri
- Department of Mechanical and Aerospace Engineering, New York University Tandon School of Engineering, Brooklyn, New York 11201, USA
| | - Igor Belykh
- Department of Mathematics and Statistics, Georgia State University, P.O. Box 4110, Atlanta, Georgia 30302-4110, USA
| |
Collapse
|
23
|
Arola-Fernández L, Díaz-Guilera A, Arenas A. Synchronization invariance under network structural transformations. Phys Rev E 2018; 97:060301. [PMID: 30011485 DOI: 10.1103/physreve.97.060301] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2018] [Indexed: 06/08/2023]
Abstract
Synchronization processes are ubiquitous despite the many connectivity patterns that complex systems can show. Usually, the emergence of synchrony is a macroscopic observable; however, the microscopic details of the system, as, e.g., the underlying network of interactions, is many times partially or totally unknown. We already know that different interaction structures can give rise to a common functionality, understood as a common macroscopic observable. Building upon this fact, here we propose network transformations that keep the collective behavior of a large system of Kuramoto oscillators invariant. We derive a method based on information theory principles, that allows us to adjust the weights of the structural interactions to map random homogeneous in-degree networks into random heterogeneous networks and vice versa, keeping synchronization values invariant. The results of the proposed transformations reveal an interesting principle; heterogeneous networks can be mapped to homogeneous ones with local information, but the reverse process needs to exploit higher-order information. The formalism provides analytical insight to tackle real complex scenarios when dealing with uncertainty in the measurements of the underlying connectivity structure.
Collapse
Affiliation(s)
- Lluís Arola-Fernández
- Departament d'Enginyeria Informàtica i Matemàtiques, Universitat Rovira i Virgili, 43007 Tarragona, Spain
| | - Albert Díaz-Guilera
- Departament de Física Fonamental, Universitat de Barcelona, Martí i Franquès 1, 08028 Barcelona, Spain
- Universitat de Barcelona Institute for Complex Systems (UBICS), Barcelona, Spain
| | - Alex Arenas
- Departament d'Enginyeria Informàtica i Matemàtiques, Universitat Rovira i Virgili, 43007 Tarragona, Spain
| |
Collapse
|
24
|
Motter AE, Timme M. Antagonistic Phenomena in Network Dynamics. ANNUAL REVIEW OF CONDENSED MATTER PHYSICS 2018; 9:463-484. [PMID: 30116502 PMCID: PMC6089548 DOI: 10.1146/annurev-conmatphys-033117-054054] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Recent research on the network modeling of complex systems has led to a convenient representation of numerous natural, social, and engineered systems that are now recognized as networks of interacting parts. Such systems can exhibit a wealth of phenomena that not only cannot be anticipated from merely examining their parts, as per the textbook definition of complexity, but also challenge intuition even when considered in the context of what is now known in network science. Here we review the recent literature on two major classes of such phenomena that have far-reaching implications: (i) antagonistic responses to changes of states or parameters and (ii) coexistence of seemingly incongruous behaviors or properties-both deriving from the collective and inherently decentralized nature of the dynamics. They include effects as diverse as negative compressibility in engineered materials, rescue interactions in biological networks, negative resistance in fluid networks, and the Braess paradox occurring across transport and supply networks. They also include remote synchronization, chimera states and the converse of symmetry breaking in brain, power-grid and oscillator networks as well as remote control in biological and bio-inspired systems. By offering a unified view of these various scenarios, we suggest that they are representative of a yet broader class of unprecedented network phenomena that ought to be revealed and explained by future research.
Collapse
Affiliation(s)
- Adilson E Motter
- Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208, USA
- Northwestern Institute on Complex Systems, Northwestern University, Evanston, Illinois 60208, USA
| | - Marc Timme
- Chair of Network Dynamics, Institute for Theoretical Physics and Center for Advancing Electronics (cfaed), Technical University of Dresden, 01062 Dresden, Germany
- Network Dynamics, Max Planck Institute for Dynamics and Self-Organization, 37077 Göttingen, Germany
| |
Collapse
|
25
|
Connection adaption for control of networked mobile chaotic agents. Sci Rep 2017; 7:16069. [PMID: 29167510 PMCID: PMC5700208 DOI: 10.1038/s41598-017-16235-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2017] [Accepted: 11/09/2017] [Indexed: 11/25/2022] Open
Abstract
In this paper, we propose a strategy for the control of mobile chaotic oscillators by adaptively rewiring connections between nearby agents with local information. In contrast to the dominant adaptive control schemes where coupling strength is adjusted continuously according to the states of the oscillators, our method does not request adaption of coupling strength. As the resulting interaction structure generated by this proposed strategy is strongly related to unidirectional chains, by investigating synchronization property of unidirectional chains, we reveal that there exists a certain coupling range in which the agents could be controlled regardless of the length of the chain. This feature enables the adaptive strategy to control the mobile oscillators regardless of their moving speed. Compared with existing adaptive control strategies for networked mobile agents, our proposed strategy is simpler for implementation where the resulting interaction networks are kept unweighted at all time.
Collapse
|
26
|
Belykh I, Jeter R, Belykh V. Foot force models of crowd dynamics on a wobbly bridge. SCIENCE ADVANCES 2017; 3:e1701512. [PMID: 29296679 PMCID: PMC5688767 DOI: 10.1126/sciadv.1701512] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/09/2017] [Accepted: 10/18/2017] [Indexed: 06/07/2023]
Abstract
Modern pedestrian and suspension bridges are designed using industry standard packages, yet disastrous resonant vibrations are observed, necessitating multimillion dollar repairs. Recent examples include pedestrian-induced vibrations during the opening of the Solférino Bridge in Paris in 1999 and the increased bouncing of the Squibb Park Bridge in Brooklyn in 2014. The most prominent example of an unstable lively bridge is the London Millennium Bridge, which started wobbling as a result of pedestrian-bridge interactions. Pedestrian phase locking due to footstep phase adjustment is suspected to be the main cause of its large lateral vibrations; however, its role in the initiation of wobbling was debated. We develop foot force models of pedestrians' response to bridge motion and detailed, yet analytically tractable, models of crowd phase locking. We use biomechanically inspired models of crowd lateral movement to investigate to what degree pedestrian synchrony must be present for a bridge to wobble significantly and what is a critical crowd size. Our results can be used as a safety guideline for designing pedestrian bridges or limiting the maximum occupancy of an existing bridge. The pedestrian models can be used as "crash test dummies" when numerically probing a specific bridge design. This is particularly important because the U.S. code for designing pedestrian bridges does not contain explicit guidelines that account for the collective pedestrian behavior.
Collapse
Affiliation(s)
- Igor Belykh
- Department of Mathematics and Statistics and Neuroscience Institute, Georgia State University, P. O. Box 4110, Atlanta, GA 30302–4110, USA
| | - Russell Jeter
- Department of Mathematics and Statistics and Neuroscience Institute, Georgia State University, P. O. Box 4110, Atlanta, GA 30302–4110, USA
| | - Vladimir Belykh
- Department of Mathematics, Volga State University of Water Transport, 5A Nesterov Street, Nizhny Novgorod 603951, Russia
- Department of Control Theory, Lobachevsky State University of Nizhny Novgorod, 23 Gagarin Avenue, Nizhny Novgorod 603950, Russia
| |
Collapse
|
27
|
Tang E, Giusti C, Baum GL, Gu S, Pollock E, Kahn AE, Roalf DR, Moore TM, Ruparel K, Gur RC, Gur RE, Satterthwaite TD, Bassett DS. Developmental increases in white matter network controllability support a growing diversity of brain dynamics. Nat Commun 2017; 8:1252. [PMID: 29093441 PMCID: PMC5665937 DOI: 10.1038/s41467-017-01254-4] [Citation(s) in RCA: 96] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2017] [Accepted: 09/01/2017] [Indexed: 11/17/2022] Open
Abstract
As the human brain develops, it increasingly supports coordinated control of neural activity. The mechanism by which white matter evolves to support this coordination is not well understood. Here we use a network representation of diffusion imaging data from 882 youth ages 8-22 to show that white matter connectivity becomes increasingly optimized for a diverse range of predicted dynamics in development. Notably, stable controllers in subcortical areas are negatively related to cognitive performance. Investigating structural mechanisms supporting these changes, we simulate network evolution with a set of growth rules. We find that all brain networks are structured in a manner highly optimized for network control, with distinct control mechanisms predicted in child vs. older youth. We demonstrate that our results cannot be explained by changes in network modularity. This work reveals a possible mechanism of human brain development that preferentially optimizes dynamic network control over static network architecture.
Collapse
Affiliation(s)
- Evelyn Tang
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Chad Giusti
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Graham L Baum
- Brain Behavior Laboratory, Department of Psychiatry, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Shi Gu
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Eli Pollock
- Department of Physics, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Ari E Kahn
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - David R Roalf
- Brain Behavior Laboratory, Department of Psychiatry, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Tyler M Moore
- Brain Behavior Laboratory, Department of Psychiatry, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Kosha Ruparel
- Brain Behavior Laboratory, Department of Psychiatry, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Ruben C Gur
- Brain Behavior Laboratory, Department of Psychiatry, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Raquel E Gur
- Brain Behavior Laboratory, Department of Psychiatry, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Theodore D Satterthwaite
- Brain Behavior Laboratory, Department of Psychiatry, University of Pennsylvania, Philadelphia, PA, 19104, USA
| | - Danielle S Bassett
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, 19104, USA.
- Department of Electrical and Systems Engineering, University of Pennsylvania, Philadelphia, PA, 19104, USA.
| |
Collapse
|
28
|
Shirasaka S, Watanabe N, Kawamura Y, Nakao H. Optimizing stability of mutual synchronization between a pair of limit-cycle oscillators with weak cross coupling. Phys Rev E 2017; 96:012223. [PMID: 29347076 DOI: 10.1103/physreve.96.012223] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2017] [Indexed: 12/22/2022]
Abstract
We consider optimization of the linear stability of synchronized states between a pair of weakly coupled limit-cycle oscillators with cross coupling, where different components of state variables of the oscillators are allowed to interact. On the basis of the phase reduction theory, we derive the coupling matrix between different components of the oscillator states that maximizes the linear stability of the synchronized state under given constraints on the overall coupling intensity and the stationary phase difference. The improvement in the linear stability is illustrated by using several types of limit-cycle oscillators as examples.
Collapse
Affiliation(s)
- Sho Shirasaka
- Research Center for Advanced Science and Technology, University of Tokyo, Tokyo 153-8904, Japan
| | - Nobuhiro Watanabe
- Department of Systems and Control Engineering, Tokyo Institute of Technology, Tokyo 152-8552, Japan
| | - Yoji Kawamura
- Department of Mathematical Science and Advanced Technology, Japan Agency for Marine-Earth Science and Technology, Yokohama 236-0001, Japan
| | - Hiroya Nakao
- Department of Systems and Control Engineering, Tokyo Institute of Technology, Tokyo 152-8552, Japan.,Department of Mechanical Engineering, University of California, Santa Barbara, California 93106, USA
| |
Collapse
|
29
|
Saha S, Mishra A, Padmanaban E, Bhowmick SK, Roy PK, Dam B, Dana SK. Coupling conditions for globally stable and robust synchrony of chaotic systems. Phys Rev E 2017; 95:062204. [PMID: 28709232 DOI: 10.1103/physreve.95.062204] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2016] [Indexed: 11/07/2022]
Abstract
We propose a set of general coupling conditions to select a coupling profile (a set of coupling matrices) from the linear flow matrix of dynamical systems for realizing global stability of complete synchronization (CS) in identical systems and robustness to parameter perturbation. The coupling matrices define the coupling links between any two oscillators in a network that consists of a conventional diffusive coupling link (self-coupling link) as well as a cross-coupling link. The addition of a selective cross-coupling link in particular plays constructive roles that ensure the global stability of synchrony and furthermore enables robustness of synchrony against small to nonsmall parameter perturbation. We elaborate the general conditions for the selection of coupling profiles for two coupled systems, three- and four-node network motifs analytically as well as numerically using benchmark models, the Lorenz system, the Hindmarsh-Rose neuron model, the Shimizu-Morioka laser model, the Rössler system, and a Sprott system. The role of the cross-coupling link is, particularly, exemplified with an example of a larger network, where it saves the network from a breakdown of synchrony against large parameter perturbation in any node. The perturbed node in the network transits from CS to generalized synchronization (GS) when all the other nodes remain in CS. The GS is manifested by an amplified response of the perturbed node in a coherent state.
Collapse
Affiliation(s)
- Suman Saha
- Department of Instrumentation and Electronics Engineering, Jadavpur University, Kolkata 700090, India.,Department of Electronics, Asutosh College, Kolkata 700026, India.,Dumkal Institute of Engineering and Technology, Murshidabad 742406, India
| | - Arindam Mishra
- Department of Physics, Jadavpur University, Kolkata 700032, India
| | - E Padmanaban
- CSIR-Indian Institute of Chemical Biology, Kolkata 700032, India.,Center for Complex System Research Kolkata, Kolkata 700094, India
| | - Sourav K Bhowmick
- Department of Electronics, Asutosh College, Kolkata 700026, India.,Center for Complex System Research Kolkata, Kolkata 700094, India
| | - Prodyot K Roy
- Center for Complex System Research Kolkata, Kolkata 700094, India.,Department of Mathematics, Presidency University, Kolkata 700073, India
| | - Bivas Dam
- Department of Instrumentation and Electronics Engineering, Jadavpur University, Kolkata 700090, India
| | - Syamal K Dana
- Center for Complex System Research Kolkata, Kolkata 700094, India.,Department of Mathematics, Jadavpur University, Kolkata 700032, India
| |
Collapse
|
30
|
Zhang Y, Nishikawa T, Motter AE. Asymmetry-induced synchronization in oscillator networks. Phys Rev E 2017; 95:062215. [PMID: 28709352 DOI: 10.1103/physreve.95.062215] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2017] [Indexed: 06/07/2023]
Abstract
A scenario has recently been reported in which in order to stabilize complete synchronization of an oscillator network-a symmetric state-the symmetry of the system itself has to be broken by making the oscillators nonidentical. But how often does such behavior-which we term asymmetry-induced synchronization (AISync)-occur in oscillator networks? Here we present the first general scheme for constructing AISync systems and demonstrate that this behavior is the norm rather than the exception in a wide class of physical systems that can be seen as multilayer networks. Since a symmetric network in complete synchrony is the basic building block of cluster synchronization in more general networks, AISync should be common also in facilitating cluster synchronization by breaking the symmetry of the cluster subnetworks.
Collapse
Affiliation(s)
- Yuanzhao Zhang
- Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208, USA
| | - Takashi Nishikawa
- Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208, USA
- Northwestern Institute on Complex Systems, Northwestern University, Evanston, Illinois 60208, USA
| | - Adilson E Motter
- Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208, USA
- Northwestern Institute on Complex Systems, Northwestern University, Evanston, Illinois 60208, USA
| |
Collapse
|
31
|
Schröder M, Chakraborty S, Witthaut D, Nagler J, Timme M. Interaction Control to Synchronize Non-synchronizable Networks. Sci Rep 2016; 6:37142. [PMID: 27853266 PMCID: PMC5112558 DOI: 10.1038/srep37142] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2016] [Accepted: 10/24/2016] [Indexed: 11/18/2022] Open
Abstract
Synchronization constitutes one of the most fundamental collective dynamics across networked systems and often underlies their function. Whether a system may synchronize depends on the internal unit dynamics as well as the topology and strength of their interactions. For chaotic units with certain interaction topologies synchronization might be impossible across all interaction strengths, meaning that these networks are non-synchronizable. Here we propose the concept of interaction control, generalizing transient uncoupling, to induce desired collective dynamics in complex networks and apply it to synchronize even such non-synchronizable systems. After highlighting that non-synchronizability prevails for a wide range of networks of arbitrary size, we explain how a simple binary control may localize interactions in state space and thereby synchronize networks. Intriguingly, localizing interactions by a fixed control scheme enables stable synchronization across all connected networks regardless of topological constraints. Interaction control may thus ease the design of desired collective dynamics even without knowledge of the networks’ exact interaction topology and consequently have implications for biological and self-organizing technical systems.
Collapse
Affiliation(s)
- Malte Schröder
- Network Dynamics, Max Planck Institute for Dynamics and Self-Organization (MPIDS), 37077 Göttingen, Germany
| | - Sagar Chakraborty
- Department of Physics, Indian Institute of Technology Kanpur, U.P. 208016, India
| | - Dirk Witthaut
- Forschungszentrum Jülich, Institute for Energy and Climate Research - Systems Analysis and Technology Evaluation (IEK-STE), 52428 Jülich, Germany.,Institute for Theoretical Physics, University of Cologne, 50937 Köln, Germany
| | - Jan Nagler
- Computational Physics, IfB, ETH Zurich, Wolfgang-Pauli-Strasse 27, 8093 Zurich, Switzerland
| | - Marc Timme
- Network Dynamics, Max Planck Institute for Dynamics and Self-Organization (MPIDS), 37077 Göttingen, Germany.,Department of Physics, Technical University of Darmstadt, 64289 Darmstadt, Germany
| |
Collapse
|
32
|
del Genio CI, Gómez-Gardeñes J, Bonamassa I, Boccaletti S. Synchronization in networks with multiple interaction layers. SCIENCE ADVANCES 2016; 2:e1601679. [PMID: 28138540 PMCID: PMC5262445 DOI: 10.1126/sciadv.1601679] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/20/2016] [Accepted: 10/13/2016] [Indexed: 05/17/2023]
Abstract
The structure of many real-world systems is best captured by networks consisting of several interaction layers. Understanding how a multilayered structure of connections affects the synchronization properties of dynamical systems evolving on top of it is a highly relevant endeavor in mathematics and physics and has potential applications in several socially relevant topics, such as power grid engineering and neural dynamics. We propose a general framework to assess the stability of the synchronized state in networks with multiple interaction layers, deriving a necessary condition that generalizes the master stability function approach. We validate our method by applying it to a network of Rössler oscillators with a double layer of interactions and show that highly rich phenomenology emerges from this. This includes cases where the stability of synchronization can be induced even if both layers would have individually induced unstable synchrony, an effect genuinely arising from the true multilayer structure of the interactions among the units in the network.
Collapse
Affiliation(s)
| | - Jesús Gómez-Gardeñes
- Departamento de Física de la Materia Condensada, University of Zaragoza, 50009 Zaragoza, Spain
- Institute for Biocomputation and Physics of Complex Systems, University of Zaragoza, 50018 Zaragoza, Spain
| | - Ivan Bonamassa
- Department of Physics, Bar-Ilan University, 52900 Ramat Gan, Israel
| | - Stefano Boccaletti
- CNR–Istituto dei Sistemi Complessi, Via Madonna del Piano, 10, 50019 Sesto Fiorentino, Italy
- Embassy of Italy in Israel, 25 Hamered Street, 68125 Tel Aviv, Israel
| |
Collapse
|
33
|
Zhou J, Zou Y, Guan S, Liu Z, Boccaletti S. Synchronization in slowly switching networks of coupled oscillators. Sci Rep 2016; 6:35979. [PMID: 27779253 PMCID: PMC5078792 DOI: 10.1038/srep35979] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2016] [Accepted: 10/07/2016] [Indexed: 11/17/2022] Open
Abstract
Networks whose structure of connections evolves in time constitute a big challenge in the study of synchronization, in particular when the time scales for the evolution of the graph topology are comparable with (or even longer than) those pertinent to the units’ dynamics. We here focus on networks with a slow-switching structure, and show that the necessary conditions for synchronization, i.e. the conditions for which synchronization is locally stable, are determined by the time average of the largest Lyapunov exponents of transverse modes of the switching topologies. Comparison between fast- and slow-switching networks allows elucidating that slow-switching processes prompt synchronization in the cases where the Master Stability Function is concave, whereas fast-switching schemes facilitate synchronization for convex curves. Moreover, the condition of slow-switching enables the introduction of a control strategy for inducing synchronization in networks with arbitrary structure and coupling strength, which is of evident relevance for broad applications in real world systems.
Collapse
Affiliation(s)
- Jie Zhou
- Department of Physics, East China Normal University, Shanghai 200241, China
| | - Yong Zou
- Department of Physics, East China Normal University, Shanghai 200241, China
| | - Shuguang Guan
- Department of Physics, East China Normal University, Shanghai 200241, China
| | - Zonghua Liu
- Department of Physics, East China Normal University, Shanghai 200241, China
| | - S Boccaletti
- CNR-Institute of Complex Systems, Via Madonna del Piano, 10, 50019 Sesto Fiorentino, Florence, Italy.,The Embassy of Italy in Tel Aviv, 25 Hamered street, 68125 Tel Aviv, Israel
| |
Collapse
|
34
|
Nishikawa T, Motter AE. Network-complement transitions, symmetries, and cluster synchronization. CHAOS (WOODBURY, N.Y.) 2016; 26:094818. [PMID: 27781466 DOI: 10.1063/1.4960617] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Synchronization in networks of coupled oscillators is known to be largely determined by the spectral and symmetry properties of the interaction network. Here, we leverage this relation to study a class of networks for which the threshold coupling strength for global synchronization is the lowest among all networks with the same number of nodes and links. These networks, defined as being uniform, complete, and multi-partite (UCM), appear at each of an infinite sequence of network-complement transitions in a larger class of networks characterized by having near-optimal thresholds for global synchronization. We show that the distinct symmetry structure of the UCM networks, which by design are optimized for global synchronizability, often leads to formation of clusters of synchronous oscillators, and that such states can coexist with the state of global synchronization.
Collapse
Affiliation(s)
- Takashi Nishikawa
- Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208, USA
| | - Adilson E Motter
- Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208, USA
| |
Collapse
|
35
|
Belykh IV, Brister BN, Belykh VN. Bistability of patterns of synchrony in Kuramoto oscillators with inertia. CHAOS (WOODBURY, N.Y.) 2016; 26:094822. [PMID: 27781476 DOI: 10.1063/1.4961435] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We study the co-existence of stable patterns of synchrony in two coupled populations of identical Kuramoto oscillators with inertia. The two populations have different sizes and can split into two clusters where the oscillators synchronize within a cluster while there is a phase shift between the dynamics of the two clusters. Due to the presence of inertia, which increases the dimensionality of the oscillator dynamics, this phase shift can oscillate, inducing a breathing cluster pattern. We derive analytical conditions for the co-existence of stable two-cluster patterns with constant and oscillating phase shifts. We demonstrate that the dynamics, that governs the bistability of the phase shifts, is described by a driven pendulum equation. We also discuss the implications of our stability results to the stability of chimeras.
Collapse
Affiliation(s)
- Igor V Belykh
- Department of Mathematics and Statistics and Neuroscience Institute, Georgia State University, 30 Pryor Street, Atlanta, Georgia 30303, USA
| | - Barrett N Brister
- Department of Mathematics and Statistics and Neuroscience Institute, Georgia State University, 30 Pryor Street, Atlanta, Georgia 30303, USA
| | - Vladimir N Belykh
- Department of Control Theory, Lobachevsky State University of Nizhny Novgorod, 23, Gagarin Ave., 603950 Nizhny Novgorod, Russia
| |
Collapse
|
36
|
Skardal PS, Taylor D, Sun J. Optimal synchronization of directed complex networks. CHAOS (WOODBURY, N.Y.) 2016; 26:094807. [PMID: 27781463 PMCID: PMC4920812 DOI: 10.1063/1.4954221] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/19/2016] [Accepted: 05/16/2016] [Indexed: 05/29/2023]
Abstract
We study optimal synchronization of networks of coupled phase oscillators. We extend previous theory for optimizing the synchronization properties of undirected networks to the important case of directed networks. We derive a generalized synchrony alignment function that encodes the interplay between the network structure and the oscillators' natural frequencies and serves as an objective measure for the network's degree of synchronization. Using the generalized synchrony alignment function, we show that a network's synchronization properties can be systematically optimized. This framework also allows us to study the properties of synchrony-optimized networks, and in particular, investigate the role of directed network properties such as nodal in- and out-degrees. For instance, we find that in optimally rewired networks, the heterogeneity of the in-degree distribution roughly matches the heterogeneity of the natural frequency distribution, but no such relationship emerges for out-degrees. We also observe that a network's synchronization properties are promoted by a strong correlation between the nodal in-degrees and the natural frequencies of oscillators, whereas the relationship between the nodal out-degrees and the natural frequencies has comparatively little effect. This result is supported by our theory, which indicates that synchronization is promoted by a strong alignment of the natural frequencies with the left singular vectors corresponding to the largest singular values of the Laplacian matrix.
Collapse
Affiliation(s)
| | - Dane Taylor
- Department of Mathematics, Carolina Center for Interdisciplinary Applied Mathematics, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| | - Jie Sun
- Department of Mathematics, Clarkson University, Potsdam, New York 13699, USA
| |
Collapse
|
37
|
Girón A, Saiz H, Bacelar FS, Andrade RFS, Gómez-Gardeñes J. Synchronization unveils the organization of ecological networks with positive and negative interactions. CHAOS (WOODBURY, N.Y.) 2016; 26:065302. [PMID: 27368792 DOI: 10.1063/1.4952960] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Network science has helped to understand the organization principles of the interactions among the constituents of large complex systems. However, recently, the high resolution of the data sets collected has allowed to capture the different types of interactions coexisting within the same system. A particularly important example is that of systems with positive and negative interactions, a usual feature appearing in social, neural, and ecological systems. The interplay of links of opposite sign presents natural difficulties for generalizing typical concepts and tools applied to unsigned networks and, moreover, poses some questions intrinsic to the signed nature of the network, such as how are negative interactions balanced by positive ones so to allow the coexistence and survival of competitors/foes within the same system? Here, we show that synchronization phenomenon is an ideal benchmark for uncovering such balance and, as a byproduct, to assess which nodes play a critical role in the overall organization of the system. We illustrate our findings with the analysis of synthetic and real ecological networks in which facilitation and competitive interactions coexist.
Collapse
Affiliation(s)
- Andrea Girón
- Department of Condensed Matter Physics, University of Zaragoza, E-50009 Zaragoza, Spain
| | - Hugo Saiz
- UMR CNRS 6553 Ecosystems-Biodiversity-Evolution, University of Rennes 1, Campus de Beaulieu, Bâtiment 14A, 35042 Rennes Cedex, France
| | - Flora S Bacelar
- Instituto de Fisica, Universidade Federal da Bahia, 40210-340 Salvador, Brazil
| | - Roberto F S Andrade
- Instituto de Fisica, Universidade Federal da Bahia, 40210-340 Salvador, Brazil
| | - Jesús Gómez-Gardeñes
- Department of Condensed Matter Physics, University of Zaragoza, E-50009 Zaragoza, Spain
| |
Collapse
|
38
|
Huddy SR, Sun J. Master stability islands for amplitude death in networks of delay-coupled oscillators. Phys Rev E 2016; 93:052209. [PMID: 27300882 DOI: 10.1103/physreve.93.052209] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2016] [Indexed: 06/06/2023]
Abstract
This paper presents a master stability function (MSF) approach for analyzing the stability of amplitude death (AD) in networks of delay-coupled oscillators. Unlike the familiar MSFs for instantaneously coupled networks, which typically have a single input encoding for the effects of the eigenvalues of the network Laplacian matrix, for delay-coupled networks we show that such MSFs generally require two additional inputs: the time delay and the coupling strength. To utilize the MSF for determining the stability of AD of general networks for a chosen nonlinear system (node dynamics) and coupling function, we introduce the concept of master stability islands (MSIs), which are two-dimensional stability islands of the delay-coupling parameter space together with a third dimension ("altitude") encoding for eigenvalues that result in stable AD. We numerically compute the MSFs and visualize the corresponding MSIs for several common chaotic systems including the Rössler, the Lorenz, and Chen's system and find that it is generally possible to achieve AD and that a nonzero time delay is necessary for the stabilization of the AD states.
Collapse
Affiliation(s)
- Stanley R Huddy
- Department of Computer Sciences and Engineering, Fairleigh Dickinson University, Teaneck, New Jersey 07666, USA
| | - Jie Sun
- Department of Mathematics, Clarkson University, Potsdam, New York 13699, USA and Department of Physics, Clarkson University, Potsdam, New York 13699, USA
| |
Collapse
|
39
|
Diwadkar A, Vaidya U. Limitations and tradeoffs in synchronization of large-scale networks with uncertain links. Sci Rep 2016; 6:21157. [PMID: 27067994 PMCID: PMC4828643 DOI: 10.1038/srep21157] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2014] [Accepted: 12/21/2015] [Indexed: 11/25/2022] Open
Abstract
The synchronization of nonlinear systems connected over large-scale networks has gained popularity in a variety of applications, such as power grids, sensor networks, and biology. Stochastic uncertainty in the interconnections is a ubiquitous phenomenon observed in these physical and biological networks. We provide a size-independent network sufficient condition for the synchronization of scalar nonlinear systems with stochastic linear interactions over large-scale networks. This sufficient condition, expressed in terms of nonlinear dynamics, the Laplacian eigenvalues of the nominal interconnections, and the variance and location of the stochastic uncertainty, allows us to define a synchronization margin. We provide an analytical characterization of important trade-offs between the internal nonlinear dynamics, network topology, and uncertainty in synchronization. For nearest neighbour networks, the existence of an optimal number of neighbours with a maximum synchronization margin is demonstrated. An analytical formula for the optimal gain that produces the maximum synchronization margin allows us to compare the synchronization properties of various complex network topologies.
Collapse
Affiliation(s)
- Amit Diwadkar
- Electrical and Computer Engineering, Iowa State University Coover Hall, Ames, IA, USA 50011
| | - Umesh Vaidya
- Electrical and Computer Engineering, Iowa State University Coover Hall, Ames, IA, USA 50011
| |
Collapse
|
40
|
Taylor D, Skardal PS, Sun J. SYNCHRONIZATION OF HETEROGENEOUS OSCILLATORS UNDER NETWORK MODIFICATIONS: PERTURBATION AND OPTIMIZATION OF THE SYNCHRONY ALIGNMENT FUNCTION. SIAM JOURNAL ON APPLIED MATHEMATICS 2016; 76:1984-2008. [PMID: 27872501 PMCID: PMC5115605 DOI: 10.1137/16m1075181] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
Abstract
Synchronization is central to many complex systems in engineering physics (e.g., the power-grid, Josephson junction circuits, and electro-chemical oscillators) and biology (e.g., neuronal, circadian, and cardiac rhythms). Despite these widespread applications-for which proper functionality depends sensitively on the extent of synchronization-there remains a lack of understanding for how systems can best evolve and adapt to enhance or inhibit synchronization. We study how network modifications affect the synchronization properties of network-coupled dynamical systems that have heterogeneous node dynamics (e.g., phase oscillators with non-identical frequencies), which is often the case for real-world systems. Our approach relies on a synchrony alignment function (SAF) that quantifies the interplay between heterogeneity of the network and of the oscillators and provides an objective measure for a system's ability to synchronize. We conduct a spectral perturbation analysis of the SAF for structural network modifications including the addition and removal of edges, which subsequently ranks the edges according to their importance to synchronization. Based on this analysis, we develop gradient-descent algorithms to efficiently solve optimization problems that aim to maximize phase synchronization via network modifications. We support these and other results with numerical experiments.
Collapse
Affiliation(s)
- Dane Taylor
- Carolina Center for Interdisciplinary Applied Mathematics, Department of Mathematics, University of North Carolina, Chapel Hill, NC 27599, USA; and Statistical and Applied Mathematical Sciences Institute (SAMSI), Research Triangle Park, NC, 27709, USA
| | | | - Jie Sun
- Department of Mathematics, Clarkson University, Potsdam, NY, 13699, USA; Department of Physics, Potsdam, NY, 13699, USA; Department of Computer Science, Potsdam, NY, 13699, USA
| |
Collapse
|
41
|
Choudhary A, Sinha S. Balance of Interactions Determines Optimal Survival in Multi-Species Communities. PLoS One 2015; 10:e0145278. [PMID: 26710077 PMCID: PMC4699223 DOI: 10.1371/journal.pone.0145278] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2015] [Accepted: 12/02/2015] [Indexed: 12/02/2022] Open
Abstract
We consider a multi-species community modelled as a complex network of populations, where the links are given by a random asymmetric connectivity matrix J, with fraction 1 - C of zero entries, where C reflects the over-all connectivity of the system. The non-zero elements of J are drawn from a Gaussian distribution with mean μ and standard deviation σ. The signs of the elements Jij reflect the nature of density-dependent interactions, such as predatory-prey, mutualism or competition, and their magnitudes reflect the strength of the interaction. In this study we try to uncover the broad features of the inter-species interactions that determine the global robustness of this network, as indicated by the average number of active nodes (i.e. non-extinct species) in the network, and the total population, reflecting the biomass yield. We find that the network transitions from a completely extinct system to one where all nodes are active, as the mean interaction strength goes from negative to positive, with the transition getting sharper for increasing C and decreasing σ. We also find that the total population, displays distinct non-monotonic scaling behaviour with respect to the product μC, implying that survival is dependent not merely on the number of links, but rather on the combination of the sparseness of the connectivity matrix and the net interaction strength. Interestingly, in an intermediate window of positive μC, the total population is maximal, indicating that too little or too much positive interactions is detrimental to survival. Rather, the total population levels are optimal when the network has intermediate net positive connection strengths. At the local level we observe marked qualitative changes in dynamical patterns, ranging from anti-phase clusters of period 2 cycles and chaotic bands, to fixed points, under the variation of mean μ of the interaction strengths. We also study the correlation between synchronization and survival, and find that synchronization does not necessarily lead to extinction. Lastly, we propose an effective low dimensional map to capture the behavior of the entire network, and this provides a broad understanding of the interplay of the local dynamical patterns and the global robustness trends in the network.
Collapse
Affiliation(s)
- Anshul Choudhary
- Department of Physical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Knowledge City, SAS Nagar, Sector 81, Manauli PO 140 306, Punjab, India
| | - Sudeshna Sinha
- Department of Physical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Knowledge City, SAS Nagar, Sector 81, Manauli PO 140 306, Punjab, India
| |
Collapse
|
42
|
Allahverdyan AE, Steeg GV, Galstyan A. Memory-induced mechanism for self-sustaining activity in networks. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:062824. [PMID: 26764761 DOI: 10.1103/physreve.92.062824] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/20/2015] [Indexed: 06/05/2023]
Abstract
We study a mechanism of activity sustaining on networks inspired by a well-known model of neuronal dynamics. Our primary focus is the emergence of self-sustaining collective activity patterns, where no single node can stay active by itself, but the activity provided initially is sustained within the collective of interacting agents. In contrast to existing models of self-sustaining activity that are caused by (long) loops present in the network, here we focus on treelike structures and examine activation mechanisms that are due to temporal memory of the nodes. This approach is motivated by applications in social media, where long network loops are rare or absent. Our results suggest that under a weak behavioral noise, the nodes robustly split into several clusters, with partial synchronization of nodes within each cluster. We also study the randomly weighted version of the models where the nodes are allowed to change their connection strength (this can model attention redistribution) and show that it does facilitate the self-sustained activity.
Collapse
Affiliation(s)
- A E Allahverdyan
- Yerevan Physics Institute, Alikhanian Brothers Street 2, Yerevan 375036, Armenia
| | - G Ver Steeg
- USC Information Sciences Institute, 4676 Admiralty Way, Marina del Rey, California 90292, USA
| | - A Galstyan
- USC Information Sciences Institute, 4676 Admiralty Way, Marina del Rey, California 90292, USA
| |
Collapse
|
43
|
Witthaut D, Timme M. Nonlocal effects and countermeasures in cascading failures. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:032809. [PMID: 26465530 DOI: 10.1103/physreve.92.032809] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/17/2015] [Indexed: 06/05/2023]
Abstract
We study the propagation of cascading failures in complex supply networks with a focus on nonlocal effects occurring far away from the initial failure. It is shown that a high clustering and a small average path length of a network generally suppress nonlocal overloads. These properties are typical for many real-world networks, often called small-world networks, such that cascades propagate mostly locally in these networks. Furthermore, we analyze the spatial aspects of countermeasures based on the intentional removal of additional edges. Nonlocal actions are generally required in networks that have a low redundancy and are thus especially vulnerable to cascades.
Collapse
Affiliation(s)
- Dirk Witthaut
- Forschungszentrum Jülich, Institute for Energy and Climate Research - Systems Analysis and Technology Evaluation (IEK-STE), 52428 Jülich, Germany
- Institute for Theoretical Physics, University of Cologne, 50937 Köln, Germany
| | - Marc Timme
- Network Dynamics, Max Planck Institute for Dynamics and Self-Organization (MPIDS), Am Fassberg 17, 37077 Göttingen, Germany
- Faculty of Physics, Georg August University Göttingen, Germany
| |
Collapse
|
44
|
Werner S, Lehnertz K. Transitions between dynamical behaviors of oscillator networks induced by diversity of nodes and edges. CHAOS (WOODBURY, N.Y.) 2015; 25:073101. [PMID: 26232952 DOI: 10.1063/1.4922836] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We study the impact of dynamical and structural heterogeneity on the collective dynamics of large small-world networks of pulse-coupled integrate-and-fire oscillators endowed with refractory periods and time delay. Depending on the choice of homogeneous control parameters (here, refractoriness and coupling strength), these networks exhibit a large spectrum of dynamical behaviors, including asynchronous, partially synchronous, and fully synchronous states. Networks exhibit transitions between these dynamical behaviors upon introducing heterogeneity. We show that the probability for a network to exhibit a certain dynamical behavior (network susceptibility) is affected differently by dynamical and structural heterogeneity and depends on the respective homogeneous dynamics.
Collapse
Affiliation(s)
- Sebastian Werner
- Department of Epileptology, University of Bonn, Sigmund-Freud-Straße 25, 53105 Bonn, Germany
| | - Klaus Lehnertz
- Department of Epileptology, University of Bonn, Sigmund-Freud-Straße 25, 53105 Bonn, Germany
| |
Collapse
|
45
|
Belykh I, Reimbayev R, Zhao K. Synergistic effect of repulsive inhibition in synchronization of excitatory networks. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:062919. [PMID: 26172784 DOI: 10.1103/physreve.91.062919] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/09/2014] [Indexed: 06/04/2023]
Abstract
We show that the addition of pairwise repulsive inhibition to excitatory networks of bursting neurons induces synchrony, in contrast to one's expectations. Through stability analysis, we reveal the mechanism underlying this purely synergistic phenomenon and demonstrate that it originates from the transition between different types of bursting, caused by excitatory-inhibitory synaptic coupling. This effect is generic and observed in different models of bursting neurons and fast synaptic interactions. We also find a universal scaling law for the synchronization stability condition for large networks in terms of the number of excitatory and inhibitory inputs each neuron receives, regardless of the network size and topology. This general law is in sharp contrast with linearly coupled networks with positive (attractive) and negative (repulsive) coupling where the placement and structure of negative connections heavily affect synchronization.
Collapse
Affiliation(s)
- Igor Belykh
- Department of Mathematics and Statistics and Neuroscience Institute, Georgia State University, 30 Pryor Street, Atlanta, Georgia 30303, USA
| | - Reimbay Reimbayev
- Department of Mathematics and Statistics and Neuroscience Institute, Georgia State University, 30 Pryor Street, Atlanta, Georgia 30303, USA
| | - Kun Zhao
- Department of Mathematics and Statistics and Neuroscience Institute, Georgia State University, 30 Pryor Street, Atlanta, Georgia 30303, USA
| |
Collapse
|
46
|
Psarras AI, Karafyllidis IG. Simulation of the Dynamics of Bacterial Quorum Sensing. IEEE Trans Nanobioscience 2015; 14:440-446. [PMID: 25594975 DOI: 10.1109/tnb.2014.2385109] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Quorum sensing (QS) is a signaling mechanism that pathogenic bacteria use to communicate and synchronize the production of exofactors to attack their hosts. Understanding and controlling QS is an important step towards a possible solution to the growing problem of antibiotic resistance. QS is a cooperative effort of a bacterial population in which some of the bacteria do not participate. This phenomenon is usually studied using game theory and the non-participating bacteria are modeled as cheaters that exploit the production of common goods (exofactors) by other bacteria. Here, we take a different approach to study the QS dynamics of a growing bacterial population. We model the bacterial population as a growing graph and use spectral graph theory to compute the evolution of its synchronizability. We also treat each bacterium as a source of signaling molecules and use the diffusion equation to compute the signaling molecule distribution. We formulate a cost function based on Lagrangian dynamics that combines the time-like synchronization with the space-like diffusion of signaling molecules. Our results show that the presence of non-participating bacteria improves the homogeneity of the signaling molecule distribution preventing thus an early onset of exofactor production and has a positive effect on the optimization of QS signaling and on attack synchronization.
Collapse
|
47
|
Jabr-Hamdan A, Sun J, Ben-Avraham D. Growing networks with superjoiners. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:052812. [PMID: 25493839 DOI: 10.1103/physreve.90.052812] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/27/2014] [Indexed: 06/04/2023]
Abstract
We study the Krapivsky-Redner (KR) network growth model, but where new nodes can connect to any number of existing nodes, m, picked from a power-law distribution p(m)∼m^{-α}. Each of the m new connections is still carried out as in the KR model with probability redirection r (corresponding to degree exponent γ_{KR}=1+1/r in the original KR model). The possibility to connect to any number of nodes resembles a more realistic type of growth in several settings, such as social networks, routers networks, and networks of citations. Here we focus on the in-, out-, and total-degree distributions and on the potential tension between the degree exponent α, characterizing new connections (outgoing links), and the degree exponent γ_{KR}(r) dictated by the redirection mechanism.
Collapse
Affiliation(s)
- Ameerah Jabr-Hamdan
- Department of Physics, Clarkson University, Potsdam, New York 13699-5820, USA
| | - Jie Sun
- Department of Mathematics & Computer Science, Clarkson University, Potsdam, New York 13699-5815, USA
| | - Daniel Ben-Avraham
- Department of Physics, Clarkson University, Potsdam, New York 13699-5820, USA and Department of Mathematics & Computer Science, Clarkson University, Potsdam, New York 13699-5815, USA
| |
Collapse
|
48
|
Skardal PS, Taylor D, Sun J. Optimal synchronization of complex networks. PHYSICAL REVIEW LETTERS 2014; 113:144101. [PMID: 25325646 DOI: 10.1103/physrevlett.113.144101] [Citation(s) in RCA: 76] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2014] [Indexed: 05/16/2023]
Abstract
We study optimal synchronization in networks of heterogeneous phase oscillators. Our main result is the derivation of a synchrony alignment function that encodes the interplay between network structure and oscillators' frequencies and that can be readily optimized. We highlight its utility in two general problems: constrained frequency allocation and network design. In general, we find that synchronization is promoted by strong alignments between frequencies and the dominant Laplacian eigenvectors, as well as a matching between the heterogeneity of frequencies and network structure.
Collapse
Affiliation(s)
- Per Sebastian Skardal
- Departament d'Enginyeria Informatica i Matemátiques, Universitat Rovira i Virgili, 43007 Tarragona, Spain and Department of Applied Mathematics, University of Colorado at Boulder, Boulder, Colorado 80309, USA
| | - Dane Taylor
- Department of Applied Mathematics, University of Colorado at Boulder, Boulder, Colorado 80309, USA and Statistical and Applied Mathematical Sciences Institute, Research Triangle Park, North Carolina 27709, USA and Department of Mathematics, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| | - Jie Sun
- Department of Mathematics, Clarkson University, Potsdam, New York 13699, USA
| |
Collapse
|
49
|
Kohar V, Ji P, Choudhary A, Sinha S, Kurths J. Synchronization in time-varying networks. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:022812. [PMID: 25215786 DOI: 10.1103/physreve.90.022812] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/06/2014] [Indexed: 06/03/2023]
Abstract
We study the stability of the synchronized state in time-varying complex networks using the concept of basin stability, which is a nonlocal and nonlinear measure of stability that can be easily applied to high-dimensional systems [P. J. Menck, J. Heitzig, N. Marwan, and J. Kurths, Nature Phys. 9, 89 (2013)]. The time-varying character is included by stochastically rewiring each link with the average frequency f. We find that the time taken to reach synchronization is lowered and the stability range of the synchronized state increases considerably in dynamic networks. Further we uncover that small-world networks are much more sensitive to link changes than random ones, with the time-varying character of the network having a significant effect at much lower rewiring frequencies. At very high rewiring frequencies, random networks perform better than small-world networks and the synchronized state is stable over a much wider window of coupling strengths. Lastly we show that the stability range of the synchronized state may be quite different for small and large perturbations, and so the linear stability analysis and the basin stability criterion provide complementary indicators of stability.
Collapse
Affiliation(s)
- Vivek Kohar
- Potsdam Institute for Climate Impact Research (PIK), 14473 Potsdam, Germany and Indian Institute of Science Education and Research (IISER) Mohali, Knowledge City, SAS Nagar, Sector 81, Manauli PO 140 306, Punjab, India
| | - Peng Ji
- Potsdam Institute for Climate Impact Research (PIK), 14473 Potsdam, Germany and Department of Physics, Humboldt University, 12489 Berlin, Germany
| | - Anshul Choudhary
- Indian Institute of Science Education and Research (IISER) Mohali, Knowledge City, SAS Nagar, Sector 81, Manauli PO 140 306, Punjab, India
| | - Sudeshna Sinha
- Indian Institute of Science Education and Research (IISER) Mohali, Knowledge City, SAS Nagar, Sector 81, Manauli PO 140 306, Punjab, India
| | - Jüergen Kurths
- Potsdam Institute for Climate Impact Research (PIK), 14473 Potsdam, Germany and Department of Physics, Humboldt University, 12489 Berlin, Germany and Institute for Complex Systems and Mathematical Biology, University of Aberdeen, Aberdeen AB24 3UE, United Kingdom and Department of Control Theory, Nizhny Novgorod State University, Gagarin Avenue 23, 606950 Nizhny Novgorod, Russia
| |
Collapse
|
50
|
San Martín J, Porter MA. Convergence time towards periodic orbits in discrete dynamical systems. PLoS One 2014; 9:e92652. [PMID: 24736594 PMCID: PMC3988007 DOI: 10.1371/journal.pone.0092652] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2013] [Accepted: 02/03/2014] [Indexed: 12/05/2022] Open
Abstract
We investigate the convergence towards periodic orbits in discrete dynamical systems. We examine the probability that a randomly chosen point converges to a particular neighborhood of a periodic orbit in a fixed number of iterations, and we use linearized equations to examine the evolution near that neighborhood. The underlying idea is that points of stable periodic orbit are associated with intervals. We state and prove a theorem that details what regions of phase space are mapped into these intervals (once they are known) and how many iterations are required to get there. We also construct algorithms that allow our theoretical results to be implemented successfully in practice.
Collapse
Affiliation(s)
- Jesús San Martín
- Escuela Técnica Superior de Ingeniería y Diseño Industrial (ETSIDI), Universidad Politécnica de Madrid, Madrid, Spain
| | - Mason A. Porter
- Oxford Centre for Industrial and Applied Mathematics, Mathematical Institute, University of Oxford, Oxford, United Kingdom
| |
Collapse
|