1
|
Lohe MA. Generalized N-rotor problems, synchronized subsystems, and associated solitons. CHAOS (WOODBURY, N.Y.) 2024; 34:083114. [PMID: 39141795 DOI: 10.1063/5.0216910] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/02/2024] [Accepted: 07/23/2024] [Indexed: 08/16/2024]
Abstract
We consider systems of N particles interacting on the unit circle through 2π-periodic potentials. An example is the N-rotor problem that arises as the classical limit of coupled Josephson junctions and for various energies is known to have a wide range of behaviors such as global chaos and ergodicity, together with families of periodic solutions and transitions from order to chaos. We focus here on selected initial values for generalized systems in which the second order Euler-Lagrange equations reduce to first order equations, which we show by example can describe an ensemble of oscillators with associated emergent phenomena such as synchronization. A specific case is that of the Kuramoto model with well-known synchronization properties. We further demonstrate the relation of these models to field theories in 1+1 dimensions that allow static kink solitons satisfying first order Bogomolny equations, well-known in soliton physics, which correspond to the first order equations of the generalized N-rotor models. For the nonlinear pendulum, for example, the first order equations define the separatrix in the phase portrait of the system and correspond to kink solitons in the sine-Gordon equation.
Collapse
Affiliation(s)
- M A Lohe
- Department of Physics, The University of Adelaide, Adelaide 5005, Australia
| |
Collapse
|
2
|
Mehta D, Chen T, Tang T, Hauenstein JD. The Loss Surface of Deep Linear Networks Viewed Through the Algebraic Geometry Lens. IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE 2022; 44:5664-5680. [PMID: 33822722 DOI: 10.1109/tpami.2021.3071289] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
By using the viewpoint of modern computational algebraic geometry, we explore properties of the optimization landscapes of deep linear neural network models. After providing clarification on the various definitions of "flat" minima, we show that the geometrically flat minima, which are merely artifacts of residual continuous symmetries of the deep linear networks, can be straightforwardly removed by a generalized L2-regularization. Then, we establish upper bounds on the number of isolated stationary points of these networks with the help of algebraic geometry. Combining these upper bounds with a method in numerical algebraic geometry, we find all stationary points for modest depth and matrix size. We demonstrate that, in the presence of the non-zero regularization, deep linear networks can indeed possess local minima which are not global minima. Finally, we show that even though the number of stationary points increases as the number of neurons (regularization parameters) increases (decreases), higher index saddles are surprisingly rare.
Collapse
|
3
|
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
|
4
|
Kassabov M, Strogatz SH, Townsend A. Sufficiently dense Kuramoto networks are globally synchronizing. CHAOS (WOODBURY, N.Y.) 2021; 31:073135. [PMID: 34340322 DOI: 10.1063/5.0057659] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/23/2021] [Accepted: 07/05/2021] [Indexed: 06/13/2023]
Abstract
Consider any network of n identical Kuramoto oscillators in which each oscillator is coupled bidirectionally with unit strength to at least μ(n-1) other oscillators. There is a critical value of the connectivity, μc, such that whenever μ>μc, the system is guaranteed to converge to the all-in-phase synchronous state for almost all initial conditions, but when μ<μc, there are networks with other stable states. The precise value of the critical connectivity remains unknown, but it has been conjectured to be μc=0.75. In 2020, Lu and Steinerberger proved that μc≤0.7889, and Yoneda, Tatsukawa, and Teramae proved in 2021 that μc>0.6838. This paper proves that μc≤0.75 and explain why this is the best upper bound that one can obtain by a purely linear stability analysis.
Collapse
Affiliation(s)
- Martin Kassabov
- Department of Mathematics, Cornell University, Ithaca, New York 14853, USA
| | - Steven H Strogatz
- Department of Mathematics, Cornell University, Ithaca, New York 14853, USA
| | - Alex Townsend
- Department of Mathematics, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
5
|
Baggio G, Bassett DS, Pasqualetti F. Data-driven control of complex networks. Nat Commun 2021; 12:1429. [PMID: 33658486 PMCID: PMC7930026 DOI: 10.1038/s41467-021-21554-0] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2020] [Accepted: 01/28/2021] [Indexed: 01/31/2023] Open
Abstract
Our ability to manipulate the behavior of complex networks depends on the design of efficient control algorithms and, critically, on the availability of an accurate and tractable model of the network dynamics. While the design of control algorithms for network systems has seen notable advances in the past few years, knowledge of the network dynamics is a ubiquitous assumption that is difficult to satisfy in practice. In this paper we overcome this limitation, and develop a data-driven framework to control a complex network optimally and without any knowledge of the network dynamics. Our optimal controls are constructed using a finite set of data, where the unknown network is stimulated with arbitrary and possibly random inputs. Although our controls are provably correct for networks with linear dynamics, we also characterize their performance against noisy data and in the presence of nonlinear dynamics, as they arise in power grid and brain networks.
Collapse
Affiliation(s)
- Giacomo Baggio
- Department of Information Engineering, University of Padova, Padova, Italy
| | - Danielle S Bassett
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, USA
- Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA, USA
- Department of Electrical and Systems Engineering, University of Pennsylvania, Philadelphia, PA, USA
- Department of Neurology, University of Pennsylvania, Philadelphia, PA, USA
- Department of Psychiatry, University of Pennsylvania, Philadelphia, PA, USA
- Santa Fe Institute, Santa Fe, NM, USA
| | - Fabio Pasqualetti
- Department of Mechanical Engineering, University of California at Riverside, Riverside, CA, USA.
| |
Collapse
|
6
|
DeVille L. Consensus on simplicial complexes: Results on stability and synchronization. CHAOS (WOODBURY, N.Y.) 2021; 31:023137. [PMID: 33653041 DOI: 10.1063/5.0037433] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/13/2020] [Accepted: 02/01/2021] [Indexed: 06/12/2023]
Abstract
We consider a nonlinear flow on simplicial complexes related to the simplicial Laplacian and show that it is a generalization of various consensus and synchronization models commonly studied on networks. In particular, our model allows us to formulate flows on simplices of any dimension so that it includes edge flows, triangle flows, etc. We show that the system can be represented as the gradient flow of an energy functional and use this to deduce the stability of various steady states of the model. Finally, we demonstrate that our model contains higher-dimensional analogs of structures seen in related network models.
Collapse
Affiliation(s)
- Lee DeVille
- Department of Mathematics, University of Illinois, 1409 W. Green St., Urbana, Illinois 61801, USA
| |
Collapse
|
7
|
Reggio A, Delabays R, Jacquod P. Clusterization and phase diagram of the bimodal Kuramoto model with bounded confidence. CHAOS (WOODBURY, N.Y.) 2020; 30:093134. [PMID: 33003922 DOI: 10.1063/5.0020436] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/01/2020] [Accepted: 09/08/2020] [Indexed: 06/11/2023]
Abstract
Inspired by the Deffuant and Hegselmann-Krause models of opinion dynamics, we extend the Kuramoto model to account for confidence bounds, i.e., vanishing interactions between pairs of oscillators when their phases differ by more than a certain value. We focus on Kuramoto oscillators with peaked, bimodal distribution of natural frequencies. We show that, in this case, the fixed-points for the extended model are made of certain numbers of independent clusters of oscillators, depending on the length of the confidence bound-the interaction range-and the distance between the two peaks of the bimodal distribution of natural frequencies. This allows us to construct the phase diagram of attractive fixed-points for the bimodal Kuramoto model with bounded confidence and to analytically explain clusterization in dynamical systems with bounded confidence.
Collapse
Affiliation(s)
- André Reggio
- School of Engineering, University of Applied Sciences of Western Switzerland, CH-1950 Sion, Switzerland
| | - Robin Delabays
- School of Engineering, University of Applied Sciences of Western Switzerland, CH-1950 Sion, Switzerland
| | - Philippe Jacquod
- School of Engineering, University of Applied Sciences of Western Switzerland, CH-1950 Sion, Switzerland
| |
Collapse
|
8
|
Townsend A, Stillman M, Strogatz SH. Dense networks that do not synchronize and sparse ones that do. CHAOS (WOODBURY, N.Y.) 2020; 30:083142. [PMID: 32872810 DOI: 10.1063/5.0018322] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Accepted: 07/24/2020] [Indexed: 06/11/2023]
Abstract
Consider any network of n identical Kuramoto oscillators in which each oscillator is coupled bidirectionally with unit strength to at least μ(n-1) other oscillators. Then, there is a critical value of μ above which the system is guaranteed to converge to the in-phase synchronous state for almost all initial conditions. The precise value of μ remains unknown. In 2018, Ling, Xu, and Bandeira proved that if each oscillator is coupled to at least 79.29% of all the others, global synchrony is ensured. In 2019, Lu and Steinerberger improved this bound to 78.89%. Here, we find clues that the critical connectivity may be exactly 75%. Our methods yield a slight improvement on the best known lower bound on the critical connectivity from 68.18% to 68.28%. We also consider the opposite end of the connectivity spectrum, where the networks are sparse rather than dense. In this regime, we ask how few edges one needs to add to a ring of n oscillators to turn it into a globally synchronizing network. We prove a partial result: all the twisted states in a ring of size n=2m can be destabilized by adding just O(nlog2n) edges. To finish the proof, one needs to rule out all other candidate attractors. We have done this for n≤8 but the problem remains open for larger n. Thus, even for systems as simple as Kuramoto oscillators, much remains to be learned about dense networks that do not globally synchronize and sparse ones that do.
Collapse
Affiliation(s)
- Alex Townsend
- Department of Mathematics, Cornell University, Ithaca, New York 14853, USA
| | - Michael Stillman
- Department of Mathematics, Cornell University, Ithaca, New York 14853, USA
| | - Steven H Strogatz
- Department of Mathematics, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
9
|
Balestra C, Kaiser F, Manik D, Witthaut D. Multistability in lossy power grids and oscillator networks. CHAOS (WOODBURY, N.Y.) 2019; 29:123119. [PMID: 31893663 DOI: 10.1063/1.5122739] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2019] [Accepted: 11/14/2019] [Indexed: 06/10/2023]
Abstract
Networks of phase oscillators are studied in various contexts, in particular, in the modeling of the electric power grid. A functional grid corresponds to a stable steady state such that any bifurcation can have catastrophic consequences up to a blackout. Also, the existence of multiple steady states is undesirable as it can lead to transitions or circulatory flows. Despite the high practical importance there is still no general theory of the existence and uniqueness of steady states in such systems. Analytic results are mostly limited to grids without Ohmic losses. In this article, we introduce a method to systematically construct the solutions of the real power load-flow equations in the presence of Ohmic losses and explicitly compute them for tree and ring networks. We investigate different mechanisms leading to multistability and discuss the impact of Ohmic losses on the existence of solutions.
Collapse
Affiliation(s)
- Chiara Balestra
- Institute for Theoretical Physics, University of Cologne, 50937 Köln, Germany
| | - Franz Kaiser
- Institute for Theoretical Physics, University of Cologne, 50937 Köln, Germany
| | - Debsankha Manik
- Max Planck Institute for Dynamics and Self-Organization, Am Faßberg 17, 37077 Göttingen, Germany
| | - Dirk Witthaut
- Institute for Theoretical Physics, University of Cologne, 50937 Köln, Germany
| |
Collapse
|
10
|
Chen T. Directed acyclic decomposition of Kuramoto equations. CHAOS (WOODBURY, N.Y.) 2019; 29:093101. [PMID: 31575119 DOI: 10.1063/1.5097826] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/28/2019] [Accepted: 08/13/2019] [Indexed: 06/10/2023]
Abstract
The Kuramoto model is one of the most widely studied models for describing synchronization behaviors in a network of coupled oscillators, and it has found a wide range of applications. Finding all possible frequency synchronization configurations in a general nonuniform, heterogeneous, and sparse network is important yet challenging due to complicated nonlinear interactions. From the view point of homotopy deformation, we develop a general framework for decomposing a Kuramoto network into smaller directed acyclic subnetworks, which lays the foundation for a divide-and-conquer approach to studying the configurations of frequency synchronization of large Kuramoto networks.
Collapse
Affiliation(s)
- Tianran Chen
- Department of Mathematics and Computer Science, Auburn University at Montgomery, Montgomery, Alabama 36116, USA and Department of Mathematics Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
11
|
Wolff MF, Lind PG, Maass P. Power grid stability under perturbation of single nodes: Effects of heterogeneity and internal nodes. CHAOS (WOODBURY, N.Y.) 2018; 28:103120. [PMID: 30384670 DOI: 10.1063/1.5040689] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/19/2018] [Accepted: 10/02/2018] [Indexed: 06/08/2023]
Abstract
Power flow dynamics in electricity grids can be described by equations resembling a Kuramoto model of non-linearly coupled oscillators with inertia. The coupling of the oscillators or nodes in a power grid generally exhibits pronounced heterogeneities due to varying features of transmission lines, generators, and loads. In studies aiming at uncovering mechanisms related to failures or malfunction of power systems, these grid heterogeneities are often neglected. However, over-simplification can lead to different results away from reality. We investigate the influence of heterogeneities in power grids on stable grid functioning and show their impact on estimating grid stability. Our conclusions are drawn by comparing the stability of an Institute of Electrical and Electronics Engineers test grid with a homogenized version of this grid.
Collapse
Affiliation(s)
- Matthias F Wolff
- Fachbereich Physik, Universität Osnabrück, Barbarastrasse 7, 49076 Osnabrück, Germany
| | - Pedro G Lind
- Fachbereich Physik, Universität Osnabrück, Barbarastrasse 7, 49076 Osnabrück, Germany
| | - Philipp Maass
- Fachbereich Physik, Universität Osnabrück, Barbarastrasse 7, 49076 Osnabrück, Germany
| |
Collapse
|
12
|
Wieland S, Malerba SB, Aumaitre S, Bercegol H. Mean-field approach for frequency synchronization in complex networks of two oscillator types. Phys Rev E 2018; 97:052310. [PMID: 29906922 DOI: 10.1103/physreve.97.052310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2017] [Indexed: 06/08/2023]
Abstract
Oscillator networks with an asymmetric bipolar distribution of natural frequencies are useful representations of power grids. We propose a mean-field model that captures the onset, form, and linear stability of frequency synchronization in such oscillator networks. The model takes into account a broad class of heterogeneous connection structures and identifies a functional form as well as basic properties that synchronized regimes possess classwide. The framework also captures synchronized regimes with large phase differences that commonly appear just above the critical threshold. Additionally, the accuracy of mean-field assumptions can be gauged internally through two model quantities. With our framework, the impact of local grid structure on frequency synchronization can be systematically explored.
Collapse
Affiliation(s)
- Stefan Wieland
- SPEC, CEA, CNRS UMR 3680, Université Paris-Saclay, F-91191 Gif-sur-Yvette Cedex, France
| | - Simone Blanco Malerba
- SPEC, CEA, CNRS UMR 3680, Université Paris-Saclay, F-91191 Gif-sur-Yvette Cedex, France
| | - Sébastien Aumaitre
- SPEC, CEA, CNRS UMR 3680, Université Paris-Saclay, F-91191 Gif-sur-Yvette Cedex, France
| | - Hervé Bercegol
- SPEC, CEA, CNRS UMR 3680, Université Paris-Saclay, F-91191 Gif-sur-Yvette Cedex, France
| |
Collapse
|
13
|
Manik D, Timme M, Witthaut D. Cycle flows and multistability in oscillatory networks. CHAOS (WOODBURY, N.Y.) 2017; 27:083123. [PMID: 28863499 DOI: 10.1063/1.4994177] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We study multistability in phase locked states in networks of phase oscillators under both Kuramoto dynamics and swing equation dynamics-a popular model for studying coarse-scale dynamics of an electrical AC power grid. We first establish the existence of geometrically frustrated states in such systems-where although a steady state flow pattern exists, no fixed point exists in the dynamical variables of phases due to geometrical constraints. We then describe the stable fixed points of the system with phase differences along each edge not exceeding π/2 in terms of cycle flows-constant flows along each simple cycle-as opposed to phase angles or flows. The cycle flow formalism allows us to compute tight upper and lower bounds to the number of fixed points in ring networks. We show that long elementary cycles, strong edge weights, and spatially homogeneous distribution of natural frequencies (for the Kuramoto model) or power injections (for the oscillator model for power grids) cause such networks to have more fixed points. We generalize some of these bounds to arbitrary planar topologies and derive scaling relations in the limit of large capacity and large cycle lengths, which we show to be quite accurate by numerical computation. Finally, we present an algorithm to compute all phase locked states-both stable and unstable-for planar networks.
Collapse
Affiliation(s)
- Debsankha Manik
- Network Dynamics, Max Planck Institute for Dynamics and Self-Organization (MPIDS), 37077 Göttingen, Germany
| | - Marc Timme
- Network Dynamics, Max Planck Institute for Dynamics and Self-Organization (MPIDS), 37077 Göttingen, Germany
| | - Dirk Witthaut
- Forschungszentrum Jülich, Institute for Energy and Climate Research - Systems Analysis and Technology Evaluation (IEK-STE), 52425 Jülich, Germany
| |
Collapse
|
14
|
Xi K, Dubbeldam JLA, Lin HX. Synchronization of cyclic power grids: Equilibria and stability of the synchronous state. CHAOS (WOODBURY, N.Y.) 2017; 27:013109. [PMID: 28147501 DOI: 10.1063/1.4973770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Synchronization is essential for the proper functioning of power grids; we investigate the synchronous states and their stability for cyclic power grids. We calculate the number of stable equilibria and investigate both the linear and nonlinear stabilities of the synchronous state. The linear stability analysis shows that the stability of the state, determined by the smallest nonzero eigenvalue, is inversely proportional to the size of the network. We use the energy barrier to measure the nonlinear stability and calculate it by comparing the potential energy of the type-1 saddles with that of the stable synchronous state. We find that the energy barrier depends on the network size (N) in a more complicated fashion compared to the linear stability. In particular, when the generators and consumers are evenly distributed in an alternating way, the energy barrier decreases to a constant when N approaches infinity. For a heterogeneous distribution of generators and consumers, the energy barrier decreases with N. The more heterogeneous the distribution is, the stronger the energy barrier depends on N. Finally, we find that by comparing situations with equal line loads in cyclic and tree networks, tree networks exhibit reduced stability. This difference disappears in the limit of N→∞. This finding corroborates previous results reported in the literature and suggests that cyclic (sub)networks may be applied to enhance power transfer while maintaining stable synchronous operation.
Collapse
Affiliation(s)
- Kaihua Xi
- Delft Institute of Applied Mathematics, Delft University of Technology, Mekelweg 4, 2628 CD Delft, The Netherlands
| | - Johan L A Dubbeldam
- Delft Institute of Applied Mathematics, Delft University of Technology, Mekelweg 4, 2628 CD Delft, The Netherlands
| | - Hai Xiang Lin
- Delft Institute of Applied Mathematics, Delft University of Technology, Mekelweg 4, 2628 CD Delft, The Netherlands
| |
Collapse
|
15
|
Abstract
We study the dynamics of coupled phase oscillators on a two-dimensional Kuramoto lattice with periodic boundary conditions. For coupling strengths just below the transition to global phase-locking, we find localized spatiotemporal patterns that we call "frequency spirals." These patterns cannot be seen under time averaging; they become visible only when we examine the spatial variation of the oscillators' instantaneous frequencies, where they manifest themselves as two-armed rotating spirals. In the more familiar phase representation, they appear as wobbly periodic patterns surrounding a phase vortex. Unlike the stationary phase vortices seen in magnetic spin systems, or the rotating spiral waves seen in reaction-diffusion systems, frequency spirals librate: the phases of the oscillators surrounding the central vortex move forward and then backward, executing a periodic motion with zero winding number. We construct the simplest frequency spiral and characterize its properties using analytical and numerical methods. Simulations show that frequency spirals in large lattices behave much like this simple prototype.
Collapse
Affiliation(s)
| | - Steven H Strogatz
- Center for Applied Mathematics, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
16
|
DeVille L, Ermentrout B. Phase-locked patterns of the Kuramoto model on 3-regular graphs. CHAOS (WOODBURY, N.Y.) 2016; 26:094820. [PMID: 27781479 DOI: 10.1063/1.4961064] [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
We consider the existence of non-synchronized fixed points to the Kuramoto model defined on sparse networks: specifically, networks where each vertex has degree exactly three. We show that "most" such networks support multiple attracting phase-locked solutions that are not synchronized and study the depth and width of the basins of attraction of these phase-locked solutions. We also show that it is common in "large enough" graphs to find phase-locked solutions where one or more of the links have angle difference greater than π/2.
Collapse
Affiliation(s)
- Lee DeVille
- Department of Mathematics, University of Illinois, 1409 W. Green St., Urbana, Illinois 61801, USA
| | - Bard Ermentrout
- Department of Mathematics, University of Pittsburgh, 301 Thackeray Hall, Pittsburgh, Pennsylvania 15260, USA
| |
Collapse
|
17
|
Sidky H, Whitmer JK, Mehta D. Reliable mixture critical point computation using polynomial homotopy continuation. AIChE J 2016. [DOI: 10.1002/aic.15319] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Hythem Sidky
- Dept. of Chemical and Biomolecular Engineering; University of Notre Dame; Notre Dame IN 46556
| | - Jonathan K. Whitmer
- Dept. of Chemical and Biomolecular Engineering; University of Notre Dame; Notre Dame IN 46556
| | - Dhagash Mehta
- Dept. of Chemical and Biomolecular Engineering; University of Notre Dame; Notre Dame IN 46556
- Dept. of Applied and Computational Mathematics and Statistics; University of Notre Dame; Notre Dame IN 46556
- Centre for the Subatomic Structure of Matter, Dept. of Physics, School of Physical Sciences; University of Adelaide; Adelaide South Australia 5005 Australia
| |
Collapse
|
18
|
Coletta T, Jacquod P. Linear stability and the Braess paradox in coupled-oscillator networks and electric power grids. Phys Rev E 2016; 93:032222. [PMID: 27078359 DOI: 10.1103/physreve.93.032222] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2015] [Indexed: 06/05/2023]
Abstract
We investigate the influence that adding a new coupling has on the linear stability of the synchronous state in coupled-oscillator networks. Using a simple model, we show that, depending on its location, the new coupling can lead to enhanced or reduced stability. We extend these results to electric power grids where a new line can lead to four different scenarios corresponding to enhanced or reduced grid stability as well as increased or decreased power flows. Our analysis shows that the Braess paradox may occur in any complex coupled system, where the synchronous state may be weakened and sometimes even destroyed by additional couplings.
Collapse
Affiliation(s)
- Tommaso Coletta
- School of Engineering, University of Applied Sciences of Western Switzerland, CH-1951 Sion, Switzerland
| | - Philippe Jacquod
- School of Engineering, University of Applied Sciences of Western Switzerland, CH-1951 Sion, Switzerland
| |
Collapse
|
19
|
Metz FL, Stariolo DA. Index statistical properties of sparse random graphs. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:042153. [PMID: 26565214 DOI: 10.1103/physreve.92.042153] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/24/2015] [Indexed: 06/05/2023]
Abstract
Using the replica method, we develop an analytical approach to compute the characteristic function for the probability P(N)(K,λ) that a large N×N adjacency matrix of sparse random graphs has K eigenvalues below a threshold λ. The method allows to determine, in principle, all moments of P(N)(K,λ), from which the typical sample-to-sample fluctuations can be fully characterized. For random graph models with localized eigenvectors, we show that the index variance scales linearly with N≫1 for |λ|>0, with a model-dependent prefactor that can be exactly calculated. Explicit results are discussed for Erdös-Rényi and regular random graphs, both exhibiting a prefactor with a nonmonotonic behavior as a function of λ. These results contrast with rotationally invariant random matrices, where the index variance scales only as lnN, with an universal prefactor that is independent of λ. Numerical diagonalization results confirm the exactness of our approach and, in addition, strongly support the Gaussian nature of the index fluctuations.
Collapse
Affiliation(s)
- F L Metz
- Departamento de Física, Universidade Federal de Santa Maria, 97105-900 Santa Maria, Brazil
- Departamento de Física, Universidade Federal do Rio Grande do Sul, 91501-970 Porto Alegre, Brazil
| | - Daniel A Stariolo
- Departamento de Física, Universidade Federal do Rio Grande do Sul, 91501-970 Porto Alegre, Brazil
| |
Collapse
|