-
Synchronization of phase oscillators on complex hypergraphs
Authors:
Sabina Adhikari,
Juan G. Restrepo,
Per Sebastian Skardal
Abstract:
We study the effect of structured higher-order interactions on the collective behavior of coupled phase oscillators. By combining a hypergraph generative model with dimensionality reduction techniques, we obtain a reduced system of differential equations for the system's order parameters. We illustrate our framework with the example of a hypergraph with hyperedges of sizes 2 (links) and 3 (triangl…
▽ More
We study the effect of structured higher-order interactions on the collective behavior of coupled phase oscillators. By combining a hypergraph generative model with dimensionality reduction techniques, we obtain a reduced system of differential equations for the system's order parameters. We illustrate our framework with the example of a hypergraph with hyperedges of sizes 2 (links) and 3 (triangles). For this case, we obtain a set of 2 coupled nonlinear algebraic equations for the order parameters. For strong values of coupling via triangles, the system exhibits bistability and explosive synchronization transitions. We find conditions that lead to bistability in terms of hypergraph properties and validate our predictions with numerical simulations. Our results provide a general framework to study synchronization of phase oscillators in hypergraphs, and they can be extended to hypergraphs with hyperedges of arbitrary sizes, dynamic-structural correlations, and other features.
△ Less
Submitted 31 January, 2023; v1 submitted 1 August, 2022;
originally announced August 2022.
-
Stochastic and deterministic dynamics in networks with excitable nodes
Authors:
Milad Rahimi-Majd,
Juan G. Restrepo,
Morteza Nattagh-Najafi
Abstract:
The analysis of the dynamics of a large class of excitable systems on locally tree-like networks leads to the conclusion that at $λ=1$ a continuous phase transition takes place, where $λ$ is the largest eigenvalue of the adjacency matrix of the network. This paper is devoted to evaluate this claim for a more general case where the assumption of the linearity of the dynamical transfer function is v…
▽ More
The analysis of the dynamics of a large class of excitable systems on locally tree-like networks leads to the conclusion that at $λ=1$ a continuous phase transition takes place, where $λ$ is the largest eigenvalue of the adjacency matrix of the network. This paper is devoted to evaluate this claim for a more general case where the assumption of the linearity of the dynamical transfer function is violated with a non-linearity parameter $β$ which interpolates between stochastic ($β=0$) and deterministic ($β\rightarrow\infty$) dynamics. Our model shows a rich phase diagram with an absorbing state and extended critical and oscillatory regimes separated by transition and bifurcation lines which depend on the initial state. We test initial states with ($\mathbb{I}$) only one initial excited node, ($\mathbb{II}$) a fixed fraction ($10\%$) of excited nodes, for all of which the transition is of first order for $β>0$ with a hysteresis effect and a gap function. For the case ($\mathbb{I}$) in the thermodynamic limit the absorbing state in the only phase for all $λ$ values and $β>0$. We further develop mean-field theories for cases ($\mathbb{I}$) and ($\mathbb{II}$). For case ($\mathbb{II}$) we obtain an analytic one-dimensional map which explains the essential properties of the model, including the hysteresis diagrams and fixed points of the dynamics.
△ Less
Submitted 24 November, 2021;
originally announced December 2021.
-
Dirac synchronization is rhythmic and explosive
Authors:
Lucille Calmon,
Juan G. Restrepo,
Joaquín J. Torres,
Ginestra Bianconi
Abstract:
Topological signals defined on nodes, links and higher dimensional simplices define the dynamical state of a network or of a simplicial complex. As such, topological signals are attracting increasing attention in network theory, dynamical systems, signal processing and machine learning. Topological signals defined on the nodes are typically studied in network dynamics, while topological signals de…
▽ More
Topological signals defined on nodes, links and higher dimensional simplices define the dynamical state of a network or of a simplicial complex. As such, topological signals are attracting increasing attention in network theory, dynamical systems, signal processing and machine learning. Topological signals defined on the nodes are typically studied in network dynamics, while topological signals defined on links are much less explored. Here we investigate Dirac synchronization, describing locally coupled topological signals defined on the nodes and on the links of a network, and treated using the topological Dirac operator. The dynamics of signals defined on the nodes is affected by a phase lag depending on the dynamical state of nearby links and vice versa. We show that Dirac synchronization on a fully connected network is explosive with a hysteresis loop characterized by a discontinuous forward transition and a continuous backward transition. The analytical investigation of the phase diagram provides a theoretical understanding of this topological explosive synchronization. The model also displays an exotic coherent synchronized phase, also called rhythmic phase, characterized by non-stationary order parameters which can shed light on topological mechanisms for the emergence of brain rhythms.
△ Less
Submitted 3 September, 2022; v1 submitted 11 July, 2021;
originally announced July 2021.
-
Geometry, Topology and Simplicial Synchronization
Authors:
Ana Paula Millán,
Juan G. Restrepo,
Joaquín J. Torres,
Ginestra Bianconi
Abstract:
Simplicial synchronization reveals the role that topology and geometry have in determining the dynamical properties of simplicial complexes. Simplicial network geometry and topology are naturally encoded in the spectral properties of the graph Laplacian and of the higher-order Laplacians of simplicial complexes. Here we show how the geometry of simplicial complexes induces spectral dimensions of t…
▽ More
Simplicial synchronization reveals the role that topology and geometry have in determining the dynamical properties of simplicial complexes. Simplicial network geometry and topology are naturally encoded in the spectral properties of the graph Laplacian and of the higher-order Laplacians of simplicial complexes. Here we show how the geometry of simplicial complexes induces spectral dimensions of the simplicial complex Laplacians that are responsible for changing the phase diagram of the Kuramoto model. In particular, simplicial complexes displaying a non-trivial simplicial network geometry cannot sustain a synchronized state in the infinite network limit if their spectral dimension is smaller or equal to four. This theoretical result is here verified on the Network Geometry with Flavor simplicial complex generative model displaying emergent hyperbolic geometry. On its turn simplicial topology is shown to determine the dynamical properties of the higher-order Kuramoto model. The higher-orderKuramoto model describes synchronization of topological signals, i.e. phases not only associated to the nodes of a simplicial complexes but associated also to higher-order simplices, including links, triangles and so on. This model displays discontinuous synchronization transitions when topological signals of different dimension and/or their solenoidal and irrotational projections are coupled in an adaptive way.
△ Less
Submitted 23 February, 2022; v1 submitted 3 May, 2021;
originally announced May 2021.
-
Higher-order simplicial synchronization of coupled topological signals
Authors:
Reza Ghorbanchian,
Juan G. Restrepo,
Joaquín J. Torres,
Ginestra Bianconi
Abstract:
Simplicial complexes capture the underlying network topology and geometry of complex systems ranging from the brain to social networks. Here we show that algebraic topology is a fundamental tool to capture the higher-order dynamics of simplicial complexes. In particular we consider topological signals, i.e., dynamical signals defined on simplices of different dimension, here taken to be nodes and…
▽ More
Simplicial complexes capture the underlying network topology and geometry of complex systems ranging from the brain to social networks. Here we show that algebraic topology is a fundamental tool to capture the higher-order dynamics of simplicial complexes. In particular we consider topological signals, i.e., dynamical signals defined on simplices of different dimension, here taken to be nodes and links for simplicity. We show that coupling between signals defined on nodes and links leads to explosive topological synchronization in which phases defined on nodes synchronize simultaneously to phases defined on links at a discontinuous phase transition. We study the model on real connectomes and on simplicial complexes and network models. Finally, we provide a comprehensive theoretical approach that captures this transition on fully connected networks and on random networks treated within the annealed approximation, establishing the conditions for observing a closed hysteresis loop in the large network limit.
△ Less
Submitted 13 March, 2021; v1 submitted 2 November, 2020;
originally announced November 2020.
-
Using Machine Learning to Assess Short Term Causal Dependence and Infer Network Links
Authors:
Amitava Banerjee,
Jaideep Pathak,
Rajarshi Roy,
Juan G. Restrepo,
Edward Ott
Abstract:
We introduce and test a general machine-learning-based technique for the inference of short term causal dependence between state variables of an unknown dynamical system from time series measurements of its state variables. Our technique leverages the results of a machine learning process for short time prediction to achieve our goal. The basic idea is to use the machine learning to estimate the e…
▽ More
We introduce and test a general machine-learning-based technique for the inference of short term causal dependence between state variables of an unknown dynamical system from time series measurements of its state variables. Our technique leverages the results of a machine learning process for short time prediction to achieve our goal. The basic idea is to use the machine learning to estimate the elements of the Jacobian matrix of the dynamical flow along an orbit. The type of machine learning that we employ is reservoir computing. We present numerical tests on link inference of a network of interacting dynamical nodes. It is seen that dynamical noise can greatly enhance the effectiveness of our technique, while observational noise degrades the effectiveness. We believe that the competition between these two opposing types of noise will be the key factor determining the success of causal inference in many of the most important application situations.
△ Less
Submitted 5 December, 2019;
originally announced December 2019.
-
Shattered Time: Can a Dissipative Time Crystal Survive Many-Body Correlations?
Authors:
Kristopher Tucker,
Bihui Zhu,
Robert J. Lewis-Swan,
Jamir Marino,
Felix Jimenez,
Juan G. Restrepo,
Ana Maria Rey
Abstract:
We investigate the emergence of a time crystal in a driven-dissipative many-body spin array. In this system the interplay between incoherent spin pumping and collective emission stabilizes a synchronized non-equilibrium steady state which in the thermodynamic limit features a self-generated time-periodic pattern imposed by collective elastic interactions. In contrast to prior realizations where th…
▽ More
We investigate the emergence of a time crystal in a driven-dissipative many-body spin array. In this system the interplay between incoherent spin pumping and collective emission stabilizes a synchronized non-equilibrium steady state which in the thermodynamic limit features a self-generated time-periodic pattern imposed by collective elastic interactions. In contrast to prior realizations where the time symmetry is already broken by an external drive, here it is only spontaneously broken by the elastic exchange interactions and manifest in the two-time correlation spectrum. Employing a combination of exact numerical calculations and a second-order cumulant expansion, we investigate the impact of many-body correlations on the time crystal formation and establish a connection between the regime where it is stable and a slow growth rate of the mutual information, signalling that the time crystal studied here is an emergent semi-classical out-of-equilibrium state of matter. We also confirm the rigidity of the time crystal to single-particle dephasing. Finally, we discuss an experimental implementation using long-lived dipoles in an optical cavity.
△ Less
Submitted 18 January, 2019; v1 submitted 8 May, 2018;
originally announced May 2018.
-
Robust entropy requires strong and balanced excitatory and inhibitory synapses
Authors:
Vidit Agrawal,
Andrew B. Cowley,
Qusay Alfaori,
Juan G. Restrepo,
Daniel B. Larremore,
Woodrow L. Shew
Abstract:
It is widely appreciated that well-balanced excitation and inhibition are necessary for proper function in neural networks. However, in principle, such balance could be achieved by many possible configurations of excitatory and inhibitory strengths, and relative numbers of excitatory and inhibitory neurons. For instance, a given level of excitation could be balanced by either numerous inhibitory n…
▽ More
It is widely appreciated that well-balanced excitation and inhibition are necessary for proper function in neural networks. However, in principle, such balance could be achieved by many possible configurations of excitatory and inhibitory strengths, and relative numbers of excitatory and inhibitory neurons. For instance, a given level of excitation could be balanced by either numerous inhibitory neurons with weak synapses, or few inhibitory neurons with strong synapses. Among the continuum of different but balanced configurations, why should any particular configuration be favored? Here we address this question in the context of the entropy of network dynamics by studying an analytically tractable network of binary neurons. We find that entropy is highest at the boundary between excitation-dominant and inhibition-dominant regimes. Entropy also varies along this boundary with a trade-off between high and robust entropy: weak synapse strengths yield high network entropy which is fragile to parameter variations, while strong synapse strengths yield a lower, but more robust, network entropy. In the case where inhibitory and excitatory synapses are constrained to have similar strength, we find that a small, but non-zero fraction of inhibitory neurons, like that seen in mammalian cortex, results in robust and relatively high entropy.
△ Less
Submitted 14 April, 2018;
originally announced April 2018.
-
Dynamic regulation of resource transport induces criticality in interdependent networks of excitable units
Authors:
Yogesh S. Virkar,
Juan G. Restrepo,
Woodrow L. Shew,
Edward Ott
Abstract:
Various functions of a network of excitable units can be enhanced if the network is in the `critical regime', where excitations are, on average, neither damped nor amplified. An important question is how can such networks self-organize to operate in the critical regime. Previously it was shown that regulation via resource transport on a secondary network can robustly maintain the primary network d…
▽ More
Various functions of a network of excitable units can be enhanced if the network is in the `critical regime', where excitations are, on average, neither damped nor amplified. An important question is how can such networks self-organize to operate in the critical regime. Previously it was shown that regulation via resource transport on a secondary network can robustly maintain the primary network dynamics in a balanced state where activity doesn't grow or decay. Here we show that this inter-network regulation process robustly produces a power-law distribution of activity avalanches, as observed in experiments, over ranges of model parameters spanning orders of magnitude. We also show that the resource transport over the secondary network protects the system against the destabilizing effect of local variations in parameters and heterogeneity in network structure. For homogeneous networks, we derive a reduced 3-dimensional map which reproduces the behavior of the full system.
△ Less
Submitted 20 October, 2019; v1 submitted 6 February, 2018;
originally announced February 2018.
-
Ensemble-based estimates of eigenvector error for empirical covariance matrices
Authors:
Dane Taylor,
Juan G. Restrepo,
Francois G. Meyer
Abstract:
Covariance matrices are fundamental to the analysis and forecast of economic, physical and biological systems. Although the eigenvalues $\{λ_i\}$ and eigenvectors $\{{\bf u}_i\}$ of a covariance matrix are central to such endeavors, in practice one must inevitably approximate the covariance matrix based on data with finite sample size $n$ to obtain empirical eigenvalues $\{\tildeλ_i\}$ and eigenve…
▽ More
Covariance matrices are fundamental to the analysis and forecast of economic, physical and biological systems. Although the eigenvalues $\{λ_i\}$ and eigenvectors $\{{\bf u}_i\}$ of a covariance matrix are central to such endeavors, in practice one must inevitably approximate the covariance matrix based on data with finite sample size $n$ to obtain empirical eigenvalues $\{\tildeλ_i\}$ and eigenvectors $\{\tilde{\bf u}_i\}$, and therefore understanding the error so introduced is of central importance. We analyze eigenvector error $\|{\bf u}_i - \tilde{\bf u}_i \|^2$ while leveraging the assumption that the true covariance matrix having size $p$ is drawn from a matrix ensemble with known spectral properties---particularly, we assume the distribution of population eigenvalues weakly converges as $p\to\infty$ to a spectral density $ρ(λ)$ and that the spacing between population eigenvalues is similar to that for the Gaussian orthogonal ensemble. Our approach complements previous analyses of eigenvector error that require the full set of eigenvalues to be known, which can be computationally infeasible when $p$ is large. To provide a scalable approach for uncertainty quantification of eigenvector error, we consider a fixed eigenvalue $λ$ and approximate the distribution of the expected square error $r= \mathbb{E}\left[\| {\bf u}_i - \tilde{\bf u}_i \|^2\right]$ across the matrix ensemble for all ${\bf u}_i$ associated with $λ_i=λ$. We find, for example, that for sufficiently large matrix size $p$ and sample size $n>p$, the probability density of $r$ scales as $1/nr^2$. This power-law scaling implies that eigenvector error is extremely heterogeneous---even if $r$ is very small for most eigenvectors, it can be large for others with non-negligible probability. We support this and further results with numerical experiments.
△ Less
Submitted 28 February, 2018; v1 submitted 28 December, 2016;
originally announced December 2016.
-
Metabolite transport through glial networks stabilizes the dynamics of learning
Authors:
Yogesh S. Virkar,
Woodrow L. Shew,
Juan G. Restrepo,
Edward Ott
Abstract:
Learning and memory are acquired through long-lasting changes in synapses. In the simplest models, such synaptic potentiation typically leads to runaway excitation, but in reality there must exist processes that robustly preserve overall stability of the neural system dynamics. How is this accomplished? Various approaches to this basic question have been considered. Here we propose a particularly…
▽ More
Learning and memory are acquired through long-lasting changes in synapses. In the simplest models, such synaptic potentiation typically leads to runaway excitation, but in reality there must exist processes that robustly preserve overall stability of the neural system dynamics. How is this accomplished? Various approaches to this basic question have been considered. Here we propose a particularly compelling and natural mechanism for preserving stability of learning neural systems. This mechanism is based on the global processes by which metabolic resources are distributed to the neurons by glial cells. Specifically, we introduce and study a model comprised of two interacting networks: a model neural network interconnected by synapses which undergo spike-timing dependent plasticity (STDP); and a model glial network interconnected by gap junctions which diffusively transport metabolic resources among the glia and, ultimately, to neural synapses where they are consumed. Our main result is that the biophysical constraints imposed by diffusive transport of metabolic resources through the glial network can prevent runaway growth of synaptic strength, both during ongoing activity and during learning. Our findings suggest a previously unappreciated role for glial transport of metabolites in the feedback control stabilization of neural network dynamics during learning.
△ Less
Submitted 10 May, 2016;
originally announced May 2016.
-
The Hamiltonian Mean Field model: effect of network structure on synchronization dynamics
Authors:
Yogesh S. Virkar,
Juan G. Restrepo,
James D. Meiss
Abstract:
The Hamiltonian Mean Field (HMF) model of coupled inertial, Hamiltonian rotors is a prototype for conservative dynamics in systems with long-range interactions. We consider the case where the interactions between the rotors are governed by a network described by a weighted adjacency matrix. By studying the linear stability of the incoherent state, we find that the transition to synchrony occurs at…
▽ More
The Hamiltonian Mean Field (HMF) model of coupled inertial, Hamiltonian rotors is a prototype for conservative dynamics in systems with long-range interactions. We consider the case where the interactions between the rotors are governed by a network described by a weighted adjacency matrix. By studying the linear stability of the incoherent state, we find that the transition to synchrony occurs at a coupling constant $K$ inversely proportional to the largest eigenvalue of the adjacency matrix. We derive a closed system of equations for a set of local order parameters and use these equations to study the effect of network heterogeneity on the synchronization of the rotors. We find that for values of $K$ just beyond the transition to synchronization the degree of synchronization is highly dependent on the network's heterogeneity, but that for large values of $K$ the degree of synchronization is robust to changes in the heterogeneity of the network's degree distribution. Our results are illustrated with numerical simulations on Erdös-Renyi networks and networks with power-law degree distributions.
△ Less
Submitted 16 March, 2015;
originally announced March 2015.
-
Synchronization of Interacting Quantum Dipoles
Authors:
Bihui Zhu,
Johannes Schachenmayer,
Minghui Xu,
F. Herrera,
Juan G. Restrepo,
Murray J. Holland,
Ana Maria Rey
Abstract:
Macroscopic ensembles of radiating dipoles are ubiquitous in the physical and natural sciences. In the classical limit the dipoles can be described as damped-driven oscillators, which are able to spontaneously synchronize and collectively lock their phases. Here we investigate the correspond- ing phenomenon in the quantum regime with arrays of quantized two-level systems coupled via long-range and…
▽ More
Macroscopic ensembles of radiating dipoles are ubiquitous in the physical and natural sciences. In the classical limit the dipoles can be described as damped-driven oscillators, which are able to spontaneously synchronize and collectively lock their phases. Here we investigate the correspond- ing phenomenon in the quantum regime with arrays of quantized two-level systems coupled via long-range and anisotropic dipolar interactions. Our calculations demonstrate that the dipoles may overcome the decoherence induced by quantum fluctuations and inhomogeneous couplings and evolve to a synchronized steady-state. This steady-state bears much similarity to that observed in classical systems, and yet also exhibits genuine quantum properties such as quantum correlations and quan- tum phase diffusion (reminiscent of lasing). Our predictions could be relevant for the development of better atomic clocks and a variety of noise tolerant quantum devices.
△ Less
Submitted 17 July, 2015; v1 submitted 20 February, 2015;
originally announced February 2015.
-
Spatiotemporal Dynamics of Calcium-Driven Cardiac Alternans
Authors:
Per Sebastian Skardal,
Alain Karma,
Juan G. Restrepo
Abstract:
We investigate the dynamics of spatially discordant alternans (SDA) driven by an instability of intracellular calcium cycling using both amplitude equations [P. S. Skardal, A. Karma, and J. G. Restrepo, Phys. Rev. Lett. 108, 108103 (2012)] and ionic model simulations. We find that, close to the alternans bifurcation, SDA is manifested as a smooth wavy modulation of the amplitudes of both repolariz…
▽ More
We investigate the dynamics of spatially discordant alternans (SDA) driven by an instability of intracellular calcium cycling using both amplitude equations [P. S. Skardal, A. Karma, and J. G. Restrepo, Phys. Rev. Lett. 108, 108103 (2012)] and ionic model simulations. We find that, close to the alternans bifurcation, SDA is manifested as a smooth wavy modulation of the amplitudes of both repolarization and calcium transient (CaT) alternans, similarly to the well-studied case of voltage-driven alternans. In contrast, further away from the bifurcation, the amplitude of CaT alternans jumps discontinuously at the nodes separating out-of-phase regions, while the amplitude of repolarization alternans remains smooth. We show that node motion of discontinuous SDA patterns is strongly hysteretic even in homogeneous tissue due to the novel phenomenon of "unidirectional pinning": node movement can only be induced towards, but not away from, the pacing site in response to a change of pacing rate or physiological parameter. In addition, we show that the wavelength of discontinuous SDA patterns scales linearly with the conduction velocity restitution length scale, in contrast to the wavelength of smooth patterns that scales sub-linearly with this length scale. Those results are also shown to be robust against cell-to-cell fluctuations owing to the property that unidirectional node motion collapses multiple jumps accumulating in nodal regions into a single jump. Amplitude equation predictions are in good overall agreement with ionic model simulations. Finally, we briefly discuss physiological implications of our findings. In particular, we suggest that due to the tendency of conduction blocks to form near nodes, the presence of unidirectional pinning makes calcium-driven alternans potentially more arrhythmogenic than voltage-driven alternans.
△ Less
Submitted 26 April, 2014; v1 submitted 26 March, 2014;
originally announced March 2014.
-
Onset of Synchronization in the Disordered Hamiltonian Mean Field Model
Authors:
Juan G. Restrepo,
James D. Meiss
Abstract:
We study the Hamiltonian Mean Field (HMF) model of coupled Hamiltonian rotors with a heterogeneous distribution of moments of inertia and coupling strengths. We show that when the parameters of the rotors are heterogeneous, finite size fluctuations can greatly modify the coupling strength at which the incoherent state loses stability by inducing correlations between the momenta and parameters of t…
▽ More
We study the Hamiltonian Mean Field (HMF) model of coupled Hamiltonian rotors with a heterogeneous distribution of moments of inertia and coupling strengths. We show that when the parameters of the rotors are heterogeneous, finite size fluctuations can greatly modify the coupling strength at which the incoherent state loses stability by inducing correlations between the momenta and parameters of the rotors. When the distribution of initial frequencies of the oscillators is sufficiently narrow, an analytical expression for the modification in critical coupling strength is obtained that confirms numerical simulations. We find that heterogeneity in the moments of inertia tends to stabilize the incoherent state, while heterogeneity in the coupling strengths tends to destabilize the incoherent state. Numerical simulations show that these effects disappear for a wide, bimodal frequency distribution.
△ Less
Submitted 11 May, 2014; v1 submitted 13 January, 2014;
originally announced January 2014.
-
Inhibition causes ceaseless dynamics in networks of excitable nodes
Authors:
Daniel B. Larremore,
Woodrow L. Shew,
Edward Ott,
Francesco Sorrentino,
Juan G. Restrepo
Abstract:
The collective dynamics of a network of excitable nodes changes dramatically when inhibitory nodes are introduced. We consider inhibitory nodes which may be activated just like excitatory nodes but, upon activating, decrease the probability of activation of network neighbors. We show that, although the direct effect of inhibitory nodes is to decrease activity, the collective dynamics becomes self-…
▽ More
The collective dynamics of a network of excitable nodes changes dramatically when inhibitory nodes are introduced. We consider inhibitory nodes which may be activated just like excitatory nodes but, upon activating, decrease the probability of activation of network neighbors. We show that, although the direct effect of inhibitory nodes is to decrease activity, the collective dynamics becomes self-sustaining. We explain this counterintuitive result by defining and analyzing a "branching function" which may be thought of as an activity-dependent branching ratio. The shape of the branching function implies that for a range of global coupling parameters dynamics are self-sustaining. Within the self-sustaining region of parameter space lies a critical line along which dynamics take the form of avalanches with universal scaling of size and duration, embedded in ceaseless timeseries of activity. Our analyses, confirmed by numerical simulation, suggest that inhibition may play a counterintuitive role in excitable networks.
△ Less
Submitted 21 February, 2014; v1 submitted 29 July, 2013;
originally announced July 2013.
-
Statistical Properties of Avalanches in Networks
Authors:
Daniel B. Larremore,
Marshall Y. Carpenter,
Edward Ott,
Juan G. Restrepo
Abstract:
We characterize the distributions of size and duration of avalanches propagating in complex networks. By an avalanche we mean the sequence of events initiated by the externally stimulated `excitation' of a network node, which may, with some probability, then stimulate subsequent firings of the nodes to which it is connected, resulting in a cascade of firings. This type of process is relevant to a…
▽ More
We characterize the distributions of size and duration of avalanches propagating in complex networks. By an avalanche we mean the sequence of events initiated by the externally stimulated `excitation' of a network node, which may, with some probability, then stimulate subsequent firings of the nodes to which it is connected, resulting in a cascade of firings. This type of process is relevant to a wide variety of situations, including neuroscience, cascading failures on electrical power grids, and epidemology. We find that the statistics of avalanches can be characterized in terms of the largest eigenvalue and corresponding eigenvector of an appropriate adjacency matrix which encodes the structure of the network. By using mean-field analyses, previous studies of avalanches in networks have not considered the effect of network structure on the distribution of size and duration of avalanches. Our results apply to individual networks (rather than network ensembles) and provide expressions for the distributions of size and duration of avalanches starting at particular nodes in the network. These findings might find application in the analysis of branching processes in networks, such as cascading power grid failures and critical brain dynamics. In particular, our results show that some experimental signatures of critical brain dynamics (i.e., power-law distributions of size and duration of neuronal avalanches), are robust to complex underlying network topologies.
△ Less
Submitted 17 April, 2012;
originally announced April 2012.
-
A network-specific approach to percolation in networks with bidirectional links
Authors:
Dane Taylor,
Juan G. Restrepo
Abstract:
Methods for determining the percolation threshold usually study the behavior of network ensembles and are often restricted to a particular type of probabilistic node/link removal strategy. We propose a network-specific method to determine the connectivity of nodes below the percolation threshold and offer an estimate to the percolation threshold in networks with bidirectional links. Our analysis d…
▽ More
Methods for determining the percolation threshold usually study the behavior of network ensembles and are often restricted to a particular type of probabilistic node/link removal strategy. We propose a network-specific method to determine the connectivity of nodes below the percolation threshold and offer an estimate to the percolation threshold in networks with bidirectional links. Our analysis does not require the assumption that a network belongs to a specific ensemble and can at the same time easily handle arbitrary removal strategies (previously an open problem for undirected networks). In validating our analysis, we find that it predicts the effects of many known complex structures (e.g., degree correlations) and may be used to study both probabilistic and deterministic attacks.
△ Less
Submitted 7 March, 2012; v1 submitted 23 August, 2011;
originally announced August 2011.
-
Network connectivity during mergers and growth: optimizing the addition of a module
Authors:
Dane Taylor,
Juan G. Restrepo
Abstract:
The principal eigenvalue $λ$ of a network's adjacency matrix often determines dynamics on the network (e.g., in synchronization and spreading processes) and some of its structural properties (e.g., robustness against failure or attack) and is therefore a good indicator for how ``strongly'' a network is connected. We study how $λ$ is modified by the addition of a module, or community, which has bro…
▽ More
The principal eigenvalue $λ$ of a network's adjacency matrix often determines dynamics on the network (e.g., in synchronization and spreading processes) and some of its structural properties (e.g., robustness against failure or attack) and is therefore a good indicator for how ``strongly'' a network is connected. We study how $λ$ is modified by the addition of a module, or community, which has broad applications, ranging from those involving a single modification (e.g., introduction of a drug into a biological process) to those involving repeated additions (e.g., power-grid and transit development). We describe how to optimally connect the module to the network to either maximize or minimize the shift in $λ$, noting several applications of directing dynamics on networks.
△ Less
Submitted 6 June, 2011; v1 submitted 23 February, 2011;
originally announced February 2011.
-
Effects of network topology, transmission delays, and refractoriness on the response of coupled excitable systems to a stochastic stimulus
Authors:
Daniel B. Larremore,
Woodrow L. Shew,
Edward Ott,
Juan G. Restrepo
Abstract:
We study the effects of network topology on the response of networks of coupled discrete excitable systems to an external stochastic stimulus. We extend recent results that characterize the response in terms of spectral properties of the adjacency matrix by allowing distributions in the transmission delays and in the number of refractory states, and by developing a nonperturbative approximation to…
▽ More
We study the effects of network topology on the response of networks of coupled discrete excitable systems to an external stochastic stimulus. We extend recent results that characterize the response in terms of spectral properties of the adjacency matrix by allowing distributions in the transmission delays and in the number of refractory states, and by developing a nonperturbative approximation to the steady state network response. We confirm our theoretical results with numerical simulations. We find that the steady state response amplitude is inversely proportional to the duration of refractoriness, which reduces the maximum attainable dynamic range. We also find that transmission delays alter the time required to reach steady state. Importantly, neither delays nor refractoriness impact the general prediction that criticality and maximum dynamic range occur when the largest eigenvalue of the adjacency matrix is unity.
△ Less
Submitted 20 May, 2011; v1 submitted 14 February, 2011;
originally announced February 2011.
-
Predicting criticality and dynamic range in complex networks: effects of topology
Authors:
Daniel B. Larremore,
Woodrow L. Shew,
Juan G. Restrepo
Abstract:
The collective dynamics of a network of coupled excitable systems in response to an external stimulus depends on the topology of the connections in the network. Here we develop a general theoretical approach to study the effects of network topology on dynamic range, which quantifies the range of stimulus intensities resulting in distinguishable network responses. We find that the largest eigenvalu…
▽ More
The collective dynamics of a network of coupled excitable systems in response to an external stimulus depends on the topology of the connections in the network. Here we develop a general theoretical approach to study the effects of network topology on dynamic range, which quantifies the range of stimulus intensities resulting in distinguishable network responses. We find that the largest eigenvalue of the weighted network adjacency matrix governs the network dynamic range. Specifically, a largest eigenvalue equal to one corresponds to a critical regime with maximum dynamic range. We gain deeper insight on the effects of network topology using a nonlinear analysis in terms of additional spectral properties of the adjacency matrix. We find that homogeneous networks can reach a higher dynamic range than those with heterogeneous topology. Our analysis, confirmed by numerical simulations, generalizes previous studies in terms of the largest eigenvalue of the adjacency matrix.
△ Less
Submitted 15 October, 2010; v1 submitted 30 July, 2010;
originally announced August 2010.
-
Spontaneous synchronization of coupled oscillator systems with frequency adaptation
Authors:
Dane Taylor,
Edward Ott,
Juan G. Restrepo
Abstract:
We study the synchronization of Kuramoto oscillators with all-to-all coupling in the presence of slow, noisy frequency adaptation. In this paper we develop a new model for oscillators which adapt both their phases and frequencies. It is found that this model naturally reproduces some observed phenomena that are not qualitatively produced by the standard Kuramoto model, such as long waiting times…
▽ More
We study the synchronization of Kuramoto oscillators with all-to-all coupling in the presence of slow, noisy frequency adaptation. In this paper we develop a new model for oscillators which adapt both their phases and frequencies. It is found that this model naturally reproduces some observed phenomena that are not qualitatively produced by the standard Kuramoto model, such as long waiting times before the synchronization of clapping audiences. By assuming a self-consistent steady state solution, we find three stability regimes for the coupling constant k, separated by critical points k1 and k2: (i) for k<k1, only the stable incoherent state exists; (ii) for k>k2, the incoherent state becomes unstable and only the synchronized state exists; (iii) for k1<k<k2, both the incoherent and synchronized states are stable. In the bistable regime spontaneous transitions between the incoherent and synchronized states are observed for finite ensembles. These transitions are well described as a stochastic process on the order parameter r undergoing fluctuations due to the system's finite size, leading to the following conclusions: (a) in the bistable regime, the average waiting time of an incoherent-to-coherent transition can be predicted by using Kramer's escape time formula and grows exponentially with the number of oscillators; (b) when the incoherent state is unstable (k>k2), the average waiting time grows logarithmically with the number of oscillators.
△ Less
Submitted 31 January, 2010; v1 submitted 11 June, 2009;
originally announced June 2009.
-
Approximating the largest eigenvalue of network adjacency matrices
Authors:
Juan G. Restrepo,
Edward Ott,
Brian R. Hunt
Abstract:
The largest eigenvalue of the adjacency matrix of a network plays an important role in several network processes (e.g., synchronization of oscillators, percolation on directed networks, linear stability of equilibria of network coupled systems, etc.). In this paper we develop approximations to the largest eigenvalue of adjacency matrices and discuss the relationships between these approximations…
▽ More
The largest eigenvalue of the adjacency matrix of a network plays an important role in several network processes (e.g., synchronization of oscillators, percolation on directed networks, linear stability of equilibria of network coupled systems, etc.). In this paper we develop approximations to the largest eigenvalue of adjacency matrices and discuss the relationships between these approximations. Numerical experiments on simulated networks are used to test our results.
△ Less
Submitted 30 May, 2007;
originally announced May 2007.
-
Weighted percolation on directed networks
Authors:
Juan G. Restrepo,
Edward Ott,
Brian R. Hunt
Abstract:
We present an analysis of the percolation transition for general node removal strategies valid for locally tree-like directed networks. On the basis of heuristic arguments we predict that, if the probability of removing node $i$ is $p_i$, the network disintegrates if $p_i$ is such that the largest eigenvalue of the matrix with entries $A_{ij}(1-p_i)$ is less than 1, where $A$ is the adjacency ma…
▽ More
We present an analysis of the percolation transition for general node removal strategies valid for locally tree-like directed networks. On the basis of heuristic arguments we predict that, if the probability of removing node $i$ is $p_i$, the network disintegrates if $p_i$ is such that the largest eigenvalue of the matrix with entries $A_{ij}(1-p_i)$ is less than 1, where $A$ is the adjacency matrix of the network. The knowledge or applicability of a Markov network model is not required by our theory, thus making it applicable to situations not covered by previous works. We test our predicted percolation criterion against numerical results for different networks and node removal strategies.
△ Less
Submitted 4 April, 2007;
originally announced April 2007.
-
Characterizing the dynamical importance of network nodes and links
Authors:
J. G. Restrepo,
E. Ott,
B. R. Hunt
Abstract:
The largest eigenvalue of the adjacency matrix of the networks is a key quantity determining several important dynamical processes on complex networks. Based on this fact, we present a quantitative, objective characterization of the dynamical importance of network nodes and links in terms of their effect on the largest eigenvalue. We show how our characterization of the dynamical importance of n…
▽ More
The largest eigenvalue of the adjacency matrix of the networks is a key quantity determining several important dynamical processes on complex networks. Based on this fact, we present a quantitative, objective characterization of the dynamical importance of network nodes and links in terms of their effect on the largest eigenvalue. We show how our characterization of the dynamical importance of nodes can be affected by degree-degree correlations and network community structure. We discuss how our characterization can be used to optimize techniques for controlling certain network dynamical processes and apply our results to real networks.
△ Less
Submitted 5 June, 2006;
originally announced June 2006.
-
The emergence of coherence in complex networks of heterogeneous dynamical systems
Authors:
Juan G. Restrepo,
Edward Ott,
Brian R. Hunt
Abstract:
We present a general theory for the onset of coherence in collections of heterogeneous maps interacting via a complex connection network. Our method allows the dynamics of the individual uncoupled systems to be either chaotic or periodic, and applies generally to networks for which the number of connections per node is large. We find that the critical coupling strength at which a transition to s…
▽ More
We present a general theory for the onset of coherence in collections of heterogeneous maps interacting via a complex connection network. Our method allows the dynamics of the individual uncoupled systems to be either chaotic or periodic, and applies generally to networks for which the number of connections per node is large. We find that the critical coupling strength at which a transition to synchrony takes place depends separately on the dynamics of the individual uncoupled systems and on the largest eigenvalue of the adjacency matrix of the coupling network. Our theory directly generalizes the Kuramoto model of equal strength, all-to-all coupled phase oscillators to the case of oscillators with more realistic dynamics coupled via a large heterogeneous network.
△ Less
Submitted 26 July, 2006; v1 submitted 27 January, 2006;
originally announced January 2006.
-
Synchronization in large directed networks of coupled phase oscillators
Authors:
Juan G. Restrepo,
Edward Ott,
Brian R. Hunt
Abstract:
We extend recent theoretical approximations describing the transition to synchronization in large undirected networks of coupled phase oscillators to the case of directed networks. We also consider extensions to networks with mixed positive/negative coupling strengths. We compare our theory with numerical simulations and find good agreement.
We extend recent theoretical approximations describing the transition to synchronization in large undirected networks of coupled phase oscillators to the case of directed networks. We also consider extensions to networks with mixed positive/negative coupling strengths. We compare our theory with numerical simulations and find good agreement.
△ Less
Submitted 9 September, 2005;
originally announced September 2005.
-
The onset of synchronization in large networks of coupled oscillators
Authors:
Juan G. Restrepo,
Brian R. Hunt,
Edward Ott
Abstract:
We study the transition from incoherence to coherence in large networks of coupled phase oscillators. We present various approximations that describe the behavior of an appropriately defined order parameter past the transition, and generalize recent results for the critical coupling strength. We find that, under appropriate conditions, the coupling strength at which the transition occurs is dete…
▽ More
We study the transition from incoherence to coherence in large networks of coupled phase oscillators. We present various approximations that describe the behavior of an appropriately defined order parameter past the transition, and generalize recent results for the critical coupling strength. We find that, under appropriate conditions, the coupling strength at which the transition occurs is determined by the largest eigenvalue of the adjacency matrix. We show how, with an additional assumption, a mean field approximation recently proposed is recovered from our results. We test our theory with numerical simulations, and find that it describes the transition when our assumptions are satisfied. We find that our theory describes the transition well in situations in which the mean field approximation fails. We study the finite size effects caused by nodes with small degree and find that they cause the critical coupling strength to increase.
△ Less
Submitted 31 January, 2005; v1 submitted 8 November, 2004;
originally announced November 2004.
-
Desynchronization waves and localized instabilities in oscillator arrays
Authors:
Juan G. Restrepo,
Edward Ott,
Brian R. Hunt
Abstract:
We consider a ring of identical or near identical coupled periodic oscillators in which the connections have randomly heterogeneous strength. We use the master stability function method to determine the possible patterns at the desynchronization transition that occurs as the coupling strengths are increased. We demonstrate Anderson localization of the modes of instability, and show that such loc…
▽ More
We consider a ring of identical or near identical coupled periodic oscillators in which the connections have randomly heterogeneous strength. We use the master stability function method to determine the possible patterns at the desynchronization transition that occurs as the coupling strengths are increased. We demonstrate Anderson localization of the modes of instability, and show that such localized instability generates waves of desynchronization that spread to the whole array. Similar results should apply to other networks with regular topology and heterogeneous connection strengths.
△ Less
Submitted 7 April, 2004;
originally announced April 2004.