1
|
Deng S, Ódor G. Chimera-like states in neural networks and power systems. Chaos 2024; 34:033135. [PMID: 38526980 DOI: 10.1063/5.0154581] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Accepted: 02/27/2024] [Indexed: 03/27/2024]
Abstract
Partial, frustrated synchronization, and chimera-like states are expected to occur in Kuramoto-like models if the spectral dimension of the underlying graph is low: ds<4. We provide numerical evidence that this really happens in the case of the high-voltage power grid of Europe (ds<2), a large human connectome (KKI113) and in the case of the largest, exactly known brain network corresponding to the fruit-fly (FF) connectome (ds<4), even though their graph dimensions are much higher, i.e., dgEU≃2.6(1) and dgFF≃5.4(1), dgKKI113≃3.4(1). We provide local synchronization results of the first- and second-order (Shinomoto) Kuramoto models by numerical solutions on the FF and the European power-grid graphs, respectively, and show the emergence of chimera-like patterns on the graph community level as well as by the local order parameters.
Collapse
Affiliation(s)
- Shengfeng Deng
- School of Physics and Information Technology, Shaanxi Normal University, Xi'an 710062, China
| | - Géza Ódor
- Institute of Technical Physics and Materials Science, HUN-REN Centre for Energy Research, P.O. Box 49, Budapest H-1525, Hungary
| |
Collapse
|
2
|
Pál G, Danku Z, Batool A, Kádár V, Yoshioka N, Ito N, Ódor G, Kun F. Scaling laws of failure dynamics on complex networks. Sci Rep 2023; 13:19733. [PMID: 37957302 PMCID: PMC10643452 DOI: 10.1038/s41598-023-47152-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2023] [Accepted: 11/09/2023] [Indexed: 11/15/2023] Open
Abstract
The topology of the network of load transmitting connections plays an essential role in the cascading failure dynamics of complex systems driven by the redistribution of load after local breakdown events. In particular, as the network structure is gradually tuned from regular to completely random a transition occurs from the localized to mean field behavior of failure spreading. Based on finite size scaling in the fiber bundle model of failure phenomena, here we demonstrate that outside the localized regime, the load bearing capacity and damage tolerance on the macro-scale, and the statistics of clusters of failed nodes on the micro-scale obey scaling laws with exponents which depend on the topology of the load transmission network and on the degree of disorder of the strength of nodes. Most notably, we show that the spatial structure of damage governs the emergence of the localized to mean field transition: as the network gets gradually randomized failed clusters formed on locally regular patches merge through long range links generating a percolation like transition which reduces the load concentration on the network. The results may help to design network structures with an improved robustness against cascading failure.
Collapse
Affiliation(s)
- Gergő Pál
- Department of Theoretical Physics, Faculty of Science and Technology, Doctoral School of Physics, University of Debrecen, P.O.Box: 400, Debrecen, H-4002, Hungary
| | - Zsuzsa Danku
- Department of Theoretical Physics, Faculty of Science and Technology, Doctoral School of Physics, University of Debrecen, P.O.Box: 400, Debrecen, H-4002, Hungary
| | - Attia Batool
- Department of Theoretical Physics, Faculty of Science and Technology, Doctoral School of Physics, University of Debrecen, P.O.Box: 400, Debrecen, H-4002, Hungary
| | - Viktória Kádár
- Department of Theoretical Physics, Faculty of Science and Technology, Doctoral School of Physics, University of Debrecen, P.O.Box: 400, Debrecen, H-4002, Hungary
| | - Naoki Yoshioka
- RIKEN Center for Computational Science, 7-1-26 Minatojima-minami-machi, Chuo-ku, Kobe, Hyogo, 650-0047, Japan
| | - Nobuyasu Ito
- RIKEN Center for Computational Science, 7-1-26 Minatojima-minami-machi, Chuo-ku, Kobe, Hyogo, 650-0047, Japan
| | - Géza Ódor
- Centre for Energy Research, Institute of Technical Physics and Materials Science, P.O. Box 49, H-1525, Budapest, Hungary
| | - Ferenc Kun
- Department of Theoretical Physics, Faculty of Science and Technology, Doctoral School of Physics, University of Debrecen, P.O.Box: 400, Debrecen, H-4002, Hungary.
- Institute for Nuclear Research (Atomki), P.O. Box 51, Debrecen, H-4001, Hungary.
| |
Collapse
|
3
|
Ódor G, Deng S. Synchronization Transition of the Second-Order Kuramoto Model on Lattices. Entropy (Basel) 2023; 25:164. [PMID: 36673304 PMCID: PMC9857586 DOI: 10.3390/e25010164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/02/2022] [Revised: 01/04/2023] [Accepted: 01/12/2023] [Indexed: 06/17/2023]
Abstract
The second-order Kuramoto equation describes the synchronization of coupled oscillators with inertia, which occur, for example, in power grids. On the contrary to the first-order Kuramoto equation, its synchronization transition behavior is significantly less known. In the case of Gaussian self-frequencies, it is discontinuous, in contrast to the continuous transition for the first-order Kuramoto equation. Herein, we investigate this transition on large 2D and 3D lattices and provide numerical evidence of hybrid phase transitions, whereby the oscillator phases θi exhibit a crossover, while the frequency is spread over a real phase transition in 3D. Thus, a lower critical dimension dlO=2 is expected for the frequencies and dlR=4 for phases such as that in the massless case. We provide numerical estimates for the critical exponents, finding that the frequency spread decays as ∼t-d/2 in the case of an aligned initial state of the phases in agreement with the linear approximation. In 3D, however, in the case of the initially random distribution of θi, we find a faster decay, characterized by ∼t-1.8(1) as the consequence of enhanced nonlinearities which appear by the random phase fluctuations.
Collapse
|
4
|
Deng S, Ódor G. Critical behavior of the diffusive susceptible-infected-recovered model. Phys Rev E 2023; 107:014303. [PMID: 36797889 DOI: 10.1103/physreve.107.014303] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2022] [Accepted: 01/04/2023] [Indexed: 06/18/2023]
Abstract
The critical behavior of the nondiffusive susceptible-infected-recovered model on lattices had been well established in virtue of its duality symmetry. By performing simulations and scaling analyses for the diffusive variant on the two-dimensional lattice, we show that diffusion for all agents, while rendering this symmetry destroyed, constitutes a singular perturbation that induces asymptotically distinct dynamical and stationary critical behavior from the nondiffusive model. In particular, the manifested crossover behavior in the effective mean-square radius exponents reveals that slow crossover behavior in general diffusive multispecies reaction systems may be ascribed to the interference of multiple length scales and timescales at early times.
Collapse
Affiliation(s)
- Shengfeng Deng
- Institute of Technical Physics and Materials Science, Centre for Energy Research, P.O. Box 49, H-1525 Budapest, Hungary
| | - Géza Ódor
- Institute of Technical Physics and Materials Science, Centre for Energy Research, P.O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|
5
|
Ódor G, Deng S, Hartmann B, Kelling J. Synchronization dynamics on power grids in Europe and the United States. Phys Rev E 2022; 106:034311. [PMID: 36266845 DOI: 10.1103/physreve.106.034311] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Accepted: 08/16/2022] [Indexed: 06/16/2023]
Abstract
Dynamical simulation of the cascade failures on the Europe and United States (U.S.) high-voltage power grids has been done via solving the second-order Kuramoto equation. We show that synchronization transition happens by increasing the global coupling parameter K with metasatble states depending on the initial conditions so that hysteresis loops occur. We provide analytic results for the time dependence of frequency spread in the large-K approximation and by comparing it with numerics of d=2,3 lattices, we find agreement in the case of ordered initial conditions. However, different power-law (PL) tails occur, when the fluctuations are strong. After thermalizing the systems we allow a single line cut failure and follow the subsequent overloads with respect to threshold values T. The PDFs p(N_{f}) of the cascade failures exhibit PL tails near the synchronization transition point K_{c}. Near K_{c} the exponents of the PLs for the U.S. power grid vary with T as 1.4≤τ≤2.1, in agreement with the empirical blackout statistics, while on the Europe power grid we find somewhat steeper PLs characterized by 1.4≤τ≤2.4. Below K_{c}, we find signatures of T-dependent PLs, caused by frustrated synchronization, reminiscent of Griffiths effects. Here we also observe stability growth following the blackout cascades, similar to intentional islanding, but for K>K_{c} this does not happen. For T<T_{c}, bumps appear in the PDFs with large mean values, known as "dragon king" blackout events. We also analyze the delaying or stabilizing effects of instantaneous feedback or increased dissipation and show how local synchronization behaves on geographic maps.
Collapse
Affiliation(s)
- Géza Ódor
- Centre for Energy Research, Institute of Technical Physics and Materials Science, H-1525 Budapest, Hungary
| | - Shengfeng Deng
- Centre for Energy Research, Institute of Technical Physics and Materials Science, H-1525 Budapest, Hungary
| | - Bálint Hartmann
- Centre for Energy Research, Institute for Energy Security and Environmental Safety, H-1525 Budapest, Hungary
| | - Jeffrey Kelling
- Faculty of Natural Sciences, Technische Universität Chemnitz, 09111 Chemnitz, Germany
- Department of Information Services and Computing, Helmholtz-Zentrum Dresden-Rossendorf, 01314 Dresden, Germany
| |
Collapse
|
6
|
Ódor G. Erratum: Nonuniversal power-law dynamics of susceptible infected recovered models on hierarchical modular networks [Phys. Rev. E 103, 062112 (2021)]. Phys Rev E 2022; 106:029901. [PMID: 36110027 DOI: 10.1103/physreve.106.029901] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2022] [Indexed: 06/15/2023]
Abstract
This corrects the article DOI: 10.1103/PhysRevE.103.062112.
Collapse
|
7
|
Ódor G. Nonuniversal power-law dynamics of susceptible infected recovered models on hierarchical modular networks. Phys Rev E 2021; 103:062112. [PMID: 34271752 DOI: 10.1103/physreve.103.062112] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Accepted: 05/25/2021] [Indexed: 11/07/2022]
Abstract
Power-law (PL) time-dependent infection growth has been reported in many COVID-19 statistics. In simple susceptible infected recovered (SIR) models, the number of infections grows at the outbreak as I(t)∝t^{d-1} on d-dimensional Euclidean lattices in the endemic phase, or it follows a slower universal PL at the critical point, until finite sizes cause immunity and a crossover to an exponential decay. Heterogeneity may alter the dynamics of spreading models, and spatially inhomogeneous infection rates can cause slower decays, posing a threat of a long recovery from a pandemic. COVID-19 statistics have also provided epidemic size distributions with PL tails in several countries. Here I investigate SIR-like models on hierarchical modular networks, embedded in 2d lattices with the addition of long-range links. I show that if the topological dimension of the network is finite, average degree-dependent PL growth of prevalence emerges. Supercritically, the same exponents as those of regular graphs occur, but the topological disorder alters the critical behavior. This is also true for the epidemic size distributions. Mobility of individuals does not affect the form of the scaling behavior, except for the d=2 lattice, but it increases the magnitude of the epidemic. The addition of a superspreader hot spot also does not change the growth exponent and the exponential decay in the herd immunity regime.
Collapse
Affiliation(s)
- Géza Ódor
- Institute of Technical Physics and Materials Science, Center for Energy Research, P.O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|
8
|
Ódor G, Hartmann B. Power-Law Distributions of Dynamic Cascade Failures in Power-Grid Models. Entropy (Basel) 2020; 22:e22060666. [PMID: 33286438 PMCID: PMC7517205 DOI: 10.3390/e22060666] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/07/2020] [Revised: 06/10/2020] [Accepted: 06/11/2020] [Indexed: 11/16/2022]
Abstract
Power-law distributed cascade failures are well known in power-grid systems. Understanding this phenomena has been done by various DC threshold models, self-tuned at their critical point. Here, we attempt to describe it using an AC threshold model, with a second-order Kuramoto type equation of motion of the power-flow. We have focused on the exploration of network heterogeneity effects, starting from homogeneous two-dimensional (2D) square lattices to the US power-grid, possessing identical nodes and links, to a realistic electric power-grid obtained from the Hungarian electrical database. The last one exhibits node dependent parameters, topologically marginally on the verge of robust networks. We show that too weak quenched heterogeneity, coming solely from the probabilistic self-frequencies of nodes (2D square lattice), is not sufficient for finding power-law distributed cascades. On the other hand, too strong heterogeneity destroys the synchronization of the system. We found agreement with the empirically observed power-law failure size distributions on the US grid, as well as on the Hungarian networks near the synchronization transition point. We have also investigated the consequence of replacing the usual Gaussian self-frequencies to exponential distributed ones, describing renewable energy sources. We found a drop in the steady state synchronization averages, but the cascade size distribution, both for the US and Hungarian systems, remained insensitive and have kept the universal tails, being characterized by the exponent τ≃1.8. We have also investigated the effect of an instantaneous feedback mechanism in case of the Hungarian power-grid.
Collapse
|
9
|
Ódor G, Kelling J. Critical synchronization dynamics of the Kuramoto model on connectome and small world graphs. Sci Rep 2019; 9:19621. [PMID: 31873076 PMCID: PMC6928153 DOI: 10.1038/s41598-019-54769-9] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2019] [Accepted: 11/15/2019] [Indexed: 11/19/2022] Open
Abstract
The hypothesis, that cortical dynamics operates near criticality also suggests, that it exhibits universal critical exponents which marks the Kuramoto equation, a fundamental model for synchronization, as a prime candidate for an underlying universal model. Here, we determined the synchronization behavior of this model by solving it numerically on a large, weighted human connectome network, containing 836733 nodes, in an assumed homeostatic state. Since this graph has a topological dimension d < 4, a real synchronization phase transition is not possible in the thermodynamic limit, still we could locate a transition between partially synchronized and desynchronized states. At this crossover point we observe power-law–tailed synchronization durations, with τt ≃ 1.2(1), away from experimental values for the brain. For comparison, on a large two-dimensional lattice, having additional random, long-range links, we obtain a mean-field value: τt ≃ 1.6(1). However, below the transition of the connectome we found global coupling control-parameter dependent exponents 1 < τt ≤ 2, overlapping with the range of human brain experiments. We also studied the effects of random flipping of a small portion of link weights, mimicking a network with inhibitory interactions, and found similar results. The control-parameter dependent exponent suggests extended dynamical criticality below the transition point.
Collapse
Affiliation(s)
- Géza Ódor
- Institute of Technical Physics and Materials Science, Centre for Energy Research, P.O.Box 49, H-1525, Budapest, Hungary
| | - Jeffrey Kelling
- Department of Information Services and Computing, Helmholtz-Zentrum Dresden - Rossendorf, P.O.Box 51 01 19, 01314, Dresden, Germany.
| |
Collapse
|
10
|
Abstract
I provide numerical evidence for the robustness of the Griffiths phase (GP) reported previously in dynamical threshold model simulations on a large human brain network with N=836733 connected nodes. The model, with equalized network sensitivity, is extended in two ways: introduction of refractory states or by randomized time-dependent thresholds. The nonuniversal power-law dynamics in an extended control parameter region survives these modifications for a short refractory state and weak disorder. In case of temporal disorder the GP shrinks and for stronger heterogeneity disappears, leaving behind a mean-field type of critical transition. Activity avalanche size distributions below the critical point decay faster than in the original model, but the addition of inhibitory interactions sets it back to the range of experimental values.
Collapse
Affiliation(s)
- Géza Ódor
- Research Institute for Technical Physics and Materials Science, Centre for Energy Research of the Hungarian Academy of Sciences, P.O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|
11
|
Abstract
We have compared the phase synchronization transition of the second-order Kuramoto model on two-dimensional (2D) lattices and on large, synthetic power grid networks, generated from real data. The latter are weighted, hierarchical modular networks. Due to the inertia the synchronization transitions are of first-order type, characterized by fast relaxation and hysteresis by varying the global coupling parameter K. Finite-size scaling analysis shows that there is no real phase transition in the thermodynamic limit, unlike in the mean-field model. The order parameter and its fluctuations depend on the network size without any real singular behavior. In case of power grids the phase synchronization breaks down at lower global couplings, than in case of 2D lattices of the same sizes, but the hysteresis is much narrower or negligible due to the low connectivity of the graphs. The temporal behavior of desynchronization avalanches after a sudden quench to low K values has been followed and duration distributions with power-law tails have been detected. This suggests rare region effects, caused by frozen disorder, resulting in heavy-tailed distributions, even without a self-organization mechanism as a consequence of a catastrophic drop event in the couplings.
Collapse
Affiliation(s)
- Géza Ódor
- Centre for Energy Research of the Hungarian Academy of Sciences, P. O. Box 49, H-1525 Budapest, Hungary
| | - Bálint Hartmann
- Centre for Energy Research of the Hungarian Academy of Sciences, P. O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|
12
|
Abstract
Griffiths phases (GPs), generated by the heterogeneities on modular networks, have recently been suggested to provide a mechanism, rid of fine parameter tuning, to explain the critical behavior of complex systems. One conjectured requirement for systems with modular structures was that the network of modules must be hierarchically organized and possess finite dimension. We investigate the dynamical behavior of an activity spreading model, evolving on heterogeneous random networks with highly modular structure and organized non-hierarchically. We observe that loosely coupled modules act as effective rare-regions, slowing down the extinction of activation. As a consequence, we find extended control parameter regions with continuously changing dynamical exponents for single network realizations, preserved after finite size analyses, as in a real GP. The avalanche size distributions of spreading events exhibit robust power-law tails. Our findings relax the requirement of hierarchical organization of the modular structure, which can help to rationalize the criticality of modular systems in the framework of GPs.
Collapse
Affiliation(s)
- Wesley Cota
- Departamento de Física, Universidade Federal de Viçosa, 36570-000, Viçosa, Minas Gerais, Brazil.
| | - Géza Ódor
- MTA-EK-MFA, Centre for Energy Research of the Hungarian Academy of Sciences, H-1121, P.O. Box 49, Budapest, Hungary
| | - Silvio C Ferreira
- Departamento de Física, Universidade Federal de Viçosa, 36570-000, Viçosa, Minas Gerais, Brazil.,National Institute of Science and Technology for Complex Systems, Rio de Janeiro, Brazil
| |
Collapse
|
13
|
Abstract
Extended numerical simulations of threshold models have been performed on a human brain network with N=836733 connected nodes available from the Open Connectome Project. While in the case of simple threshold models a sharp discontinuous phase transition without any critical dynamics arises, variable threshold models exhibit extended power-law scaling regions. This is attributed to fact that Griffiths effects, stemming from the topological or interaction heterogeneity of the network, can become relevant if the input sensitivity of nodes is equalized. I have studied the effects of link directness, as well as the consequence of inhibitory connections. Nonuniversal power-law avalanche size and time distributions have been found with exponents agreeing with the values obtained in electrode experiments of the human brain. The dynamical critical region occurs in an extended control parameter space without the assumption of self-organized criticality.
Collapse
Affiliation(s)
- Géza Ódor
- Institute of Technical Physics and Materials Science, Centre for Energy Research of the Hungarian Academy of Sciences, P.O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|
14
|
Abstract
Extensive dynamical simulations of restricted solid-on-solid models in D=2+1 dimensions have been done using parallel multisurface algorithms implemented on graphics cards. Numerical evidence is presented that these models exhibit Kardar-Parisi-Zhang surface growth scaling, irrespective of the step heights N. We show that by increasing N the corrections to scaling increase, thus smaller step-sized models describe better the asymptotic, long-wave-scaling behavior.
Collapse
Affiliation(s)
- Jeffrey Kelling
- Department of Information Services and Computing, Helmholtz-Zentrum Dresden-Rossendorf, P. O. Box 51 01 19, 01314 Dresden, Germany.,Institute of Ion Beam Physics and Materials Research, Helmholtz-Zentrum Dresden-Rossendorf, P. O. Box 51 01 19, 01314 Dresden, Germany
| | - Géza Ódor
- Institute of Technical Physics and Materials Science, Centre for Energy Research of the Hungarian Academy of Sciences, P. O. Box 49, H-1525 Budapest, Hungary
| | - Sibylle Gemming
- Institute of Ion Beam Physics and Materials Research, Helmholtz-Zentrum Dresden-Rossendorf, P. O. Box 51 01 19, 01314 Dresden, Germany.,Institute of Physics, TU Chemnitz, 09107 Chemnitz, Germany
| |
Collapse
|
15
|
Cota W, Ferreira SC, Ódor G. Griffiths effects of the susceptible-infected-susceptible epidemic model on random power-law networks. Phys Rev E 2016; 93:032322. [PMID: 27078381 DOI: 10.1103/physreve.93.032322] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2015] [Indexed: 06/05/2023]
Abstract
We provide numerical evidence for slow dynamics of the susceptible-infected-susceptible model evolving on finite-size random networks with power-law degree distributions. Extensive simulations were done by averaging the activity density over many realizations of networks. We investigated the effects of outliers in both highly fluctuating (natural cutoff) and nonfluctuating (hard cutoff) most connected vertices. Logarithmic and power-law decays in time were found for natural and hard cutoffs, respectively. This happens in extended regions of the control parameter space λ(1)<λ<λ(2), suggesting Griffiths effects, induced by the topological inhomogeneities. Optimal fluctuation theory considering sample-to-sample fluctuations of the pseudothresholds is presented to explain the observed slow dynamics. A quasistationary analysis shows that response functions remain bounded at λ(2). We argue these to be signals of a smeared transition. However, in the thermodynamic limit the Griffiths effects loose their relevancy and have a conventional critical point at λ(c)=0. Since many real networks are composed by heterogeneous and weakly connected modules, the slow dynamics found in our analysis of independent and finite networks can play an important role for the deeper understanding of such systems.
Collapse
Affiliation(s)
- Wesley Cota
- Departamento de Física, Universidade Federal de Viçosa, 36570-000, Viçosa, MG, Brazil
| | - Silvio C Ferreira
- Departamento de Física, Universidade Federal de Viçosa, 36570-000, Viçosa, MG, Brazil
| | - Géza Ódor
- MTA-MFA-EK Research Institute for Technical Physics and Materials Science, H-1121 Budapest, P. O. Box 49, Hungary
| |
Collapse
|
16
|
Abstract
We study variants of hierarchical modular network models suggested by Kaiser and Hilgetag [ Front. in Neuroinform., 4 (2010) 8] to model functional brain connectivity, using extensive simulations and quenched mean-field theory (QMF), focusing on structures with a connection probability that decays exponentially with the level index. Such networks can be embedded in two-dimensional Euclidean space. We explore the dynamic behavior of the contact process (CP) and threshold models on networks of this kind, including hierarchical trees. While in the small-world networks originally proposed to model brain connectivity, the topological heterogeneities are not strong enough to induce deviations from mean-field behavior, we show that a Griffiths phase can emerge under reduced connection probabilities, approaching the percolation threshold. In this case the topological dimension of the networks is finite, and extended regions of bursty, power-law dynamics are observed. Localization in the steady state is also shown via QMF. We investigate the effects of link asymmetry and coupling disorder, and show that localization can occur even in small-world networks with high connectivity in case of link disorder.
Collapse
Affiliation(s)
- Géza Ódor
- MTA-MFA-EK Research Institute for Technical Physics and Materials Science, H-1121 Budapest, P.O. Box 49, Hungary
| | - Ronald Dickman
- Departamento de Fisica and National Institute of Science and Technology of Complex Systems, ICEx, Universidade Federal de Minas Gerais, Caixa Postal 702, 30161-970, Belo Horizonte - Minas Gerais, Brazil
| | - Gergely Ódor
- Massachusetts Institute of Technology, 77 Massachusetts Avenue Cambridge, MA 02139-4307, USA
| |
Collapse
|
17
|
Ódor G. Localization transition, Lifschitz tails, and rare-region effects in network models. Phys Rev E Stat Nonlin Soft Matter Phys 2014; 90:032110. [PMID: 25314398 DOI: 10.1103/physreve.90.032110] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/16/2014] [Indexed: 06/04/2023]
Abstract
Effects of heterogeneity in the suspected-infected-susceptible model on networks are investigated using quenched mean-field theory. The emergence of localization is described by the distributions of the inverse participation ratio and compared with the rare-region effects appearing in simulations and in the Lifschitz tails. The latter, in the linear approximation, is related to the spectral density of the Laplacian matrix and to the time dependent order parameter. I show that these approximations indicate correctly Griffiths phases both on regular one-dimensional lattices and on small-world networks exhibiting purely topological disorder. I discuss the localization transition that occurs on scale-free networks at γ=3 degree exponent.
Collapse
Affiliation(s)
- Géza Ódor
- Research Center for Natural Sciences, Hungarian Academy of Sciences, P.O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|
18
|
Ódor G. Slow, bursty dynamics as a consequence of quenched network topologies. Phys Rev E Stat Nonlin Soft Matter Phys 2014; 89:042102. [PMID: 24827188 DOI: 10.1103/physreve.89.042102] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/27/2014] [Indexed: 06/03/2023]
Abstract
Bursty dynamics of agents is shown to appear at criticality or in extended Griffiths phases, even in case of Poisson processes. I provide numerical evidence for a power-law type of intercommunication time distributions by simulating the contact process and the susceptible-infected-susceptible model. This observation suggests that in the case of nonstationary bursty systems, the observed non-Poissonian behavior can emerge as a consequence of an underlying hidden Poissonian network process, which is either critical or exhibits strong rare-region effects. On the contrary, in time-varying networks, rare-region effects do not cause deviation from the mean-field behavior, and heterogeneity-induced burstyness is absent.
Collapse
Affiliation(s)
- Géza Ódor
- Research Center for Natural Sciences, Hungarian Academy of Sciences, MTA TTK MFA, P.O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|
19
|
Ódor G, Kelling J, Gemming S. Aging of the (2+1)-dimensional Kardar-Parisi-Zhang model. Phys Rev E Stat Nonlin Soft Matter Phys 2014; 89:032146. [PMID: 24730828 DOI: 10.1103/physreve.89.032146] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/30/2013] [Indexed: 06/03/2023]
Abstract
Extended dynamical simulations have been performed on a (2+1)-dimensional driven dimer lattice-gas model to estimate aging properties. The autocorrelation and the autoresponse functions are determined and the corresponding scaling exponents are tabulated. Since this model can be mapped onto the (2+1)-dimensional Kardar-Parisi-Zhang surface growth model, our results contribute to the understanding of the universality class of that basic system.
Collapse
Affiliation(s)
- Géza Ódor
- MTA TTK MFA Research Institute for Natural Sciences, P. O. Box 49, H-1525 Budapest, Hungary
| | - Jeffrey Kelling
- Institute of Ion Beam Physics and Materials Research Helmholtz-Zentrum, Dresden-Rossendorf, P. O. Box 51 01 19, 01314 Dresden, Germany and Institute of Physics, TU Chemnitz 09107 Chemnitz, Germany
| | - Sibylle Gemming
- Institute of Ion Beam Physics and Materials Research Helmholtz-Zentrum, Dresden-Rossendorf, P. O. Box 51 01 19, 01314 Dresden, Germany and Institute of Physics, TU Chemnitz 09107 Chemnitz, Germany
| |
Collapse
|
20
|
Ódor G. Rare regions of the susceptible-infected-susceptible model on Barabási-Albert networks. Phys Rev E Stat Nonlin Soft Matter Phys 2013; 87:042132. [PMID: 23679396 DOI: 10.1103/physreve.87.042132] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2013] [Revised: 04/05/2013] [Indexed: 06/02/2023]
Abstract
I extend a previous work to susceptible-infected-susceptible (SIS) models on weighted Barabási-Albert scale-free networks. Numerical evidence is provided that phases with slow, power-law dynamics emerge as the consequence of quenched disorder and tree topologies studied previously with the contact process. I compare simulation results with spectral analysis of the networks and show that the quenched mean-field (QMF) approximation provides a reliable, relatively fast method to explore activity clustering. This suggests that QMF can be used for describing rare-region effects due to network inhomogeneities. Finite-size study of the QMF shows the expected disappearance of the epidemic threshold λ(c) in the thermodynamic limit and an inverse participation ratio ~0.25, meaning localization in case of disassortative weight scheme. Contrarily, for the multiplicative weights and the unweighted trees, this value vanishes in the thermodynamic limit, suggesting only weak rare-region effects in agreement with the dynamical simulations. Strong corrections to the mean-field behavior in case of disassortative weights explains the concave shape of the order parameter ρ(λ) at the transition point. Application of this method to other models may reveal interesting rare-region effects, Griffiths phases as the consequence of quenched topological heterogeneities.
Collapse
Affiliation(s)
- Géza Ódor
- Research Centre for Natural Sciences, Hungarian Academy of Sciences, MTA TTK MFA, P.O. Box 49, H-1525 Budapest, Hungary
| |
Collapse
|
21
|
|
22
|
Juhász R, Ódor G, Castellano C, Muñoz MA. Rare-region effects in the contact process on networks. Phys Rev E Stat Nonlin Soft Matter Phys 2012; 85:066125. [PMID: 23005180 DOI: 10.1103/physreve.85.066125] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/22/2011] [Indexed: 06/01/2023]
Abstract
Networks and dynamical processes occurring on them have become a paradigmatic representation of complex systems. Studying the role of quenched disorder, both intrinsic to nodes and topological, is a key challenge. With this in mind, here we analyze the contact process (i.e., the simplest model for propagation phenomena) with node-dependent infection rates (i.e., intrinsic quenched disorder) on complex networks. We find Griffiths phases and other rare-region effects, leading rather generically to anomalously slow (algebraic, logarithmic, etc.) relaxation, on Erdős-Rényi networks. We predict similar effects to exist for other topologies as long as a nonvanishing percolation threshold exists. More strikingly, we find that Griffiths phases can also emerge--even with constant epidemic rates--as a consequence of mere topological heterogeneity. In particular, we find Griffiths phases in finite-dimensional networks as, for instance, a family of generalized small-world networks. These results have a broad spectrum of implications for propagation phenomena and other dynamical processes on networks, and are relevant for the analysis of both models and empirical data.
Collapse
Affiliation(s)
- Róbert Juhász
- Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, H-1525 Budapest, P.O. Box 49, Hungary
| | | | | | | |
Collapse
|
23
|
Menyhárd N, Ódor G. Non-Markovian persistence at the parity conserving point of a one-dimensional nonequilibrium kinetic Ising model. ACTA ACUST UNITED AC 1999. [DOI: 10.1088/0305-4470/30/24/015] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
24
|
Menyhárd N, Ódor G. Phase transitions and critical behaviour in one-dimensional non-equilibrium kinetic Ising models with branching annihilating random walk of kinks. ACTA ACUST UNITED AC 1999. [DOI: 10.1088/0305-4470/29/23/030] [Citation(s) in RCA: 46] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
25
|
Ódor G, Szolnoki A. Directed-percolation conjecture for cellular automata. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics 1996; 53:2231-2238. [PMID: 9964504 DOI: 10.1103/physreve.53.2231] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
26
|
Ódor G. Estimation of the order-parameter exponent of critical cellular automata using the enhanced coherent anomaly method. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics 1995; 51:6261-6262. [PMID: 9963370 DOI: 10.1103/physreve.51.6261] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
27
|
Ódor G, Szabó G. Universality change in stochastic cellular automaton with applied site exchange. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics 1994; 49:R3555-R3557. [PMID: 9961782 DOI: 10.1103/physreve.49.r3555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
28
|
Szabó G, Ódor G. Extended mean-field study of a stochastic cellular automaton. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics 1994; 49:2764-2768. [PMID: 9961541 DOI: 10.1103/physreve.49.2764] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
29
|
Ódor G, Boccara N, Szabó G. Phase-transition study of a one-dimensional probabilistic site-exchange cellular automaton. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics 1993; 48:3168-3171. [PMID: 9960955 DOI: 10.1103/physreve.48.3168] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
|
30
|
|