1
|
Papo D, Buldú JM. Does the brain behave like a (complex) network? I. Dynamics. Phys Life Rev 2024; 48:47-98. [PMID: 38145591 DOI: 10.1016/j.plrev.2023.12.006] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2023] [Accepted: 12/10/2023] [Indexed: 12/27/2023]
Abstract
Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.
Collapse
Affiliation(s)
- D Papo
- Department of Neuroscience and Rehabilitation, Section of Physiology, University of Ferrara, Ferrara, Italy; Center for Translational Neurophysiology, Fondazione Istituto Italiano di Tecnologia, Ferrara, Italy.
| | - J M Buldú
- Complex Systems Group & G.I.S.C., Universidad Rey Juan Carlos, Madrid, Spain
| |
Collapse
|
2
|
Safari A, Moretti P, Diez I, Cortes JM, Muñoz MA. Persistence of hierarchical network organization and emergent topologies in models of functional connectivity. Neurocomputing 2021. [DOI: 10.1016/j.neucom.2021.02.096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
3
|
St-Onge G, Thibeault V, Allard A, Dubé LJ, Hébert-Dufresne L. Master equation analysis of mesoscopic localization in contagion dynamics on higher-order networks. Phys Rev E 2021; 103:032301. [PMID: 33862710 DOI: 10.1103/physreve.103.032301] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Accepted: 01/04/2021] [Indexed: 12/19/2022]
Abstract
Simple models of infectious diseases tend to assume random mixing of individuals, but real interactions are not random pairwise encounters: they occur within various types of gatherings such as workplaces, households, schools, and concerts, best described by a higher-order network structure. We model contagions on higher-order networks using group-based approximate master equations, in which we track all states and interactions within a group of nodes and assume a mean-field coupling between them. Using the susceptible-infected-susceptible dynamics, our approach reveals the existence of a mesoscopic localization regime, where a disease can concentrate and self-sustain only around large groups in the network overall organization. In this regime, the phase transition is smeared, characterized by an inhomogeneous activation of the groups. At the mesoscopic level, we observe that the distribution of infected nodes within groups of the same size can be very dispersed, even bimodal. When considering heterogeneous networks, both at the level of nodes and at the level of groups, we characterize analytically the region associated with mesoscopic localization in the structural parameter space. We put in perspective this phenomenon with eigenvector localization and discuss how a focus on higher-order structures is needed to discern the more subtle localization at the mesoscopic level. Finally, we discuss how mesoscopic localization affects the response to structural interventions and how this framework could provide important insights for a broad range of dynamics.
Collapse
Affiliation(s)
- Guillaume St-Onge
- Département de physique, de génie physique et d'optique, Université Laval, Québec (Québec), Canada G1V 0A6.,Centre interdisciplinaire en modélisation mathématique, Université Laval, Québec (Québec), Canada G1V 0A6
| | - Vincent Thibeault
- Département de physique, de génie physique et d'optique, Université Laval, Québec (Québec), Canada G1V 0A6.,Centre interdisciplinaire en modélisation mathématique, Université Laval, Québec (Québec), Canada G1V 0A6
| | - Antoine Allard
- Département de physique, de génie physique et d'optique, Université Laval, Québec (Québec), Canada G1V 0A6.,Centre interdisciplinaire en modélisation mathématique, Université Laval, Québec (Québec), Canada G1V 0A6
| | - Louis J Dubé
- Département de physique, de génie physique et d'optique, Université Laval, Québec (Québec), Canada G1V 0A6.,Centre interdisciplinaire en modélisation mathématique, Université Laval, Québec (Québec), Canada G1V 0A6
| | - Laurent Hébert-Dufresne
- Département de physique, de génie physique et d'optique, Université Laval, Québec (Québec), Canada G1V 0A6.,Vermont Complex Systems Center, University of Vermont, Burlington, Vermont 05405, USA.,Department of Computer Science, University of Vermont, Burlington, Vermont 05405, USA
| |
Collapse
|
4
|
Mata AS. An overview of epidemic models with phase transitions to absorbing states running on top of complex networks. CHAOS (WOODBURY, N.Y.) 2021; 31:012101. [PMID: 33754778 DOI: 10.1063/5.0033130] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Accepted: 12/18/2020] [Indexed: 06/12/2023]
Abstract
Dynamical systems running on the top of complex networks have been extensively investigated for decades. But this topic still remains among the most relevant issues in complex network theory due to its range of applicability. The contact process (CP) and the susceptible-infected-susceptible (SIS) model are used quite often to describe epidemic dynamics. Despite their simplicity, these models are robust to predict the kernel of real situations. In this work, we review concisely both processes that are well-known and very applied examples of models that exhibit absorbing-state phase transitions. In the epidemic scenario, individuals can be infected or susceptible. A phase transition between a disease-free (absorbing) state and an active stationary phase (where a fraction of the population is infected) are separated by an epidemic threshold. For the SIS model, the central issue is to determine this epidemic threshold on heterogeneous networks. For the CP model, the main interest is to relate critical exponents with statistical properties of the network.
Collapse
Affiliation(s)
- Angélica S Mata
- Departamento de Física, Universidade Federal de Lavras, Caixa postal 3037, CEP:37200-900, Lavras, Minas Gerais, Brazil
| |
Collapse
|
5
|
St-Onge G, Young JG, Hébert-Dufresne L, Dubé LJ. Efficient sampling of spreading processes on complex networks using a composition and rejection algorithm. COMPUTER PHYSICS COMMUNICATIONS 2019; 240:30-37. [PMID: 31708586 PMCID: PMC6839824 DOI: 10.1016/j.cpc.2019.02.008] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Efficient stochastic simulation algorithms are of paramount importance to the study of spreading phenomena on complex networks. Using insights and analytical results from network science, we discuss how the structure of contacts affects the efficiency of current algorithms. We show that algorithms believed to require O ( log N ) or even O ( 1 ) operations per update-where N is the number of nodes-display instead a polynomial scaling for networks that are either dense or sparse and heterogeneous. This significantly affects the required computation time for simulations on large networks. To circumvent the issue, we propose a node-based method combined with a composition and rejection algorithm, a sampling scheme that has an average-case complexity of O [ log ( log N ) ] per update for general networks. This systematic approach is first set-up for Markovian dynamics, but can also be adapted to a number of non-Markovian processes and can enhance considerably the study of a wide range of dynamics on networks.
Collapse
Affiliation(s)
- Guillaume St-Onge
- Département de Physique, de Génie Physique, et d'Optique, Université Laval, Québec (Québec), Canada, G1V 0A6
- Centre interdisciplinaire de modélisation mathématique de l'Université Laval, Québec (Québec), Canada, G1V 0A6
| | - Jean-Gabriel Young
- Département de Physique, de Génie Physique, et d'Optique, Université Laval, Québec (Québec), Canada, G1V 0A6
- Centre interdisciplinaire de modélisation mathématique de l'Université Laval, Québec (Québec), Canada, G1V 0A6
| | - Laurent Hébert-Dufresne
- Département de Physique, de Génie Physique, et d'Optique, Université Laval, Québec (Québec), Canada, G1V 0A6
- Department of Computer Science and Vermont Complex Systems Center, University of Vermont, Burlington, VT 05401, USA
| | - Louis J Dubé
- Département de Physique, de Génie Physique, et d'Optique, Université Laval, Québec (Québec), Canada, G1V 0A6
- Centre interdisciplinaire de modélisation mathématique de l'Université Laval, Québec (Québec), Canada, G1V 0A6
| |
Collapse
|
6
|
St-Onge G, Young JG, Laurence E, Murphy C, Dubé LJ. Phase transition of the susceptible-infected-susceptible dynamics on time-varying configuration model networks. Phys Rev E 2018; 97:022305. [PMID: 29548152 DOI: 10.1103/physreve.97.022305] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2017] [Indexed: 06/08/2023]
Abstract
We present a degree-based theoretical framework to study the susceptible-infected-susceptible (SIS) dynamics on time-varying (rewired) configuration model networks. Using this framework on a given degree distribution, we provide a detailed analysis of the stationary state using the rewiring rate to explore the whole range of the time variation of the structure relative to that of the SIS process. This analysis is suitable for the characterization of the phase transition and leads to three main contributions: (1) We obtain a self-consistent expression for the absorbing-state threshold, able to capture both collective and hub activation. (2) We recover the predictions of a number of existing approaches as limiting cases of our analysis, providing thereby a unifying point of view for the SIS dynamics on random networks. (3) We obtain bounds for the critical exponents of a number of quantities in the stationary state. This allows us to reinterpret the concept of hub-dominated phase transition. Within our framework, it appears as a heterogeneous critical phenomenon: observables for different degree classes have a different scaling with the infection rate. This phenomenon is followed by the successive activation of the degree classes beyond the epidemic threshold.
Collapse
Affiliation(s)
- Guillaume St-Onge
- Département de Physique, de Génie Physique, et d'Optique, Université Laval, Québec City, Québec, Canada, G1V 0A6
| | - Jean-Gabriel Young
- Département de Physique, de Génie Physique, et d'Optique, Université Laval, Québec City, Québec, Canada, G1V 0A6
| | - Edward Laurence
- Département de Physique, de Génie Physique, et d'Optique, Université Laval, Québec City, Québec, Canada, G1V 0A6
| | - Charles Murphy
- Département de Physique, de Génie Physique, et d'Optique, Université Laval, Québec City, Québec, Canada, G1V 0A6
| | - Louis J Dubé
- Département de Physique, de Génie Physique, et d'Optique, Université Laval, Québec City, Québec, Canada, G1V 0A6
| |
Collapse
|
7
|
Slanina F. Localization in random bipartite graphs: Numerical and empirical study. Phys Rev E 2017; 95:052149. [PMID: 28618645 DOI: 10.1103/physreve.95.052149] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2016] [Indexed: 06/07/2023]
Abstract
We investigate adjacency matrices of bipartite graphs with a power-law degree distribution. Motivation for this study is twofold: first, vibrational states in granular matter and jammed sphere packings; second, graphs encoding social interaction, especially electronic commerce. We establish the position of the mobility edge and show that it strongly depends on the power in the degree distribution and on the ratio of the sizes of the two parts of the bipartite graph. At the jamming threshold, where the two parts have the same size, localization vanishes. We found that the multifractal spectrum is nontrivial in the delocalized phase, but still near the mobility edge. We also study an empirical bipartite graph, namely, the Amazon reviewer-item network. We found that in this specific graph the mobility edge disappears, and we draw a conclusion from this fact regarding earlier empirical studies of the Amazon network.
Collapse
Affiliation(s)
- František Slanina
- Institute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, CZ-18221 Praha, Czech Republic
| |
Collapse
|
8
|
Li S. Griffiths phase on hierarchical modular networks with small-world edges. Phys Rev E 2017; 95:032306. [PMID: 28415342 PMCID: PMC7217519 DOI: 10.1103/physreve.95.032306] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2016] [Revised: 12/30/2016] [Indexed: 12/21/2022]
Abstract
The Griffiths phase has been proposed to induce a stretched critical regime that facilitates self-organizing of brain networks for optimal function. This phase stems from the intrinsic structural heterogeneity of brain networks, i.e., the hierarchical modular structure. In this work, the Griffiths phase is studied in modified hierarchical networks with small-world connections based on the 3-regular Hanoi network. Through extensive simulations, the hierarchical level-dependent inter-module wiring probabilities are identified to determine the emergence of the Griffiths phase. Numerical results and the complementary spectral analysis of the relevant networks can be helpful for a deeper understanding of the essential structural characteristics of finite-dimensional networks to support the Griffiths phase.
Collapse
Affiliation(s)
- Shanshan Li
- Department of Physics, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
9
|
Ferreira RS, da Costa RA, Dorogovtsev SN, Mendes JFF. Metastable localization of diseases in complex networks. Phys Rev E 2016; 94:062305. [PMID: 28085335 DOI: 10.1103/physreve.94.062305] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2016] [Indexed: 11/07/2022]
Abstract
We describe the phenomenon of localization in the epidemic susceptible-infective-susceptible model on highly heterogeneous networks in which strongly connected nodes (hubs) play the role of centers of localization. We find that in this model the localized states below the epidemic threshold are metastable. The longevity and scale of the metastable outbreaks do not show a sharp localization transition; instead there is a smooth crossover from localized to delocalized states as we approach the epidemic threshold from below. Analyzing these long-lasting local outbreaks for a random regular graph with a hub, we show how this localization can be detected from the shape of the distribution of the number of infective nodes.
Collapse
Affiliation(s)
- R S Ferreira
- Departamento de Ciências Exatas e Aplicadas, Universidade Federal de Ouro Preto, 35931-008 João Monlevade, Brazil
| | - R A da Costa
- Department of Physics & I3N, University of Aveiro, 3810-193 Aveiro, Portugal
| | - S N Dorogovtsev
- Department of Physics & I3N, University of Aveiro, 3810-193 Aveiro, Portugal.,A. F. Ioffe Physico-Technical Institute, 194021 St. Petersburg, Russia
| | - J F F Mendes
- Department of Physics & I3N, University of Aveiro, 3810-193 Aveiro, Portugal
| |
Collapse
|
10
|
Sander RS, Costa GS, Ferreira SC. Sampling methods for the quasistationary regime of epidemic processes on regular and complex networks. Phys Rev E 2016; 94:042308. [PMID: 27841497 PMCID: PMC7217511 DOI: 10.1103/physreve.94.042308] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2016] [Indexed: 11/24/2022]
Abstract
A major hurdle in the simulation of the steady state of epidemic processes is that the system will unavoidably visit an absorbing, disease-free state at sufficiently long times due to the finite size of the networks where epidemics evolves. In the present work, we compare different quasistationary (QS) simulation methods where the absorbing states are suitably handled and the thermodynamical limit of the original dynamics can be achieved. We analyze the standard QS (SQS) method, where the sampling is constrained to active configurations, the reflecting boundary condition (RBC), where the dynamics returns to the pre-absorbing configuration, and hub reactivation (HR), where the most connected vertex of the network is reactivated after a visit to an absorbing state. We apply the methods to the contact process (CP) and susceptible-infected-susceptible (SIS) models on regular and scale free networks. The investigated methods yield the same epidemic threshold for both models. For CP, that undergoes a standard collective phase transition, the methods are equivalent. For SIS, whose phase transition is ruled by the hub mutual reactivation, the SQS and HR methods are able to capture localized epidemic phases while RBC is not. We also apply the autocorrelation time as a tool to characterize the phase transition and observe that this analysis provides the same finite-size scaling exponents for the critical relaxation time for the investigated methods. Finally, we verify the equivalence between RBC method and a weak external field for epidemics on networks.
Collapse
Affiliation(s)
- Renan S Sander
- Departamento de Física, Universidade Federal de Viçosa, 36570-900, Viçosa, Minas Gerais, Brazil
| | - Guilherme S Costa
- Departamento de Física, Universidade Federal de Viçosa, 36570-900, Viçosa, Minas Gerais, Brazil
| | - Silvio C Ferreira
- Departamento de Física, Universidade Federal de Viçosa, 36570-900, Viçosa, Minas Gerais, Brazil
| |
Collapse
|
11
|
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.3] [Reference Citation Analysis] [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
|
12
|
Mata AS, Ferreira SC. Multiple transitions of the susceptible-infected-susceptible epidemic model on complex networks. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:012816. [PMID: 25679666 DOI: 10.1103/physreve.91.012816] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/09/2014] [Indexed: 05/04/2023]
Abstract
The epidemic threshold of the susceptible-infected-susceptible (SIS) dynamics on random networks having a power law degree distribution with exponent γ>3 has been investigated using different mean-field approaches, which predict different outcomes. We performed extensive simulations in the quasistationary state for a comparison with these mean-field theories. We observed concomitant multiple transitions in individual networks presenting large gaps in the degree distribution and the obtained multiple epidemic thresholds are well described by different mean-field theories. We observed that the transitions involving thresholds which vanish at the thermodynamic limit involve localized states, in which a vanishing fraction of the network effectively contributes to epidemic activity, whereas an endemic state, with a finite density of infected vertices, occurs at a finite threshold. The multiple transitions are related to the activations of distinct subdomains of the network, which are not directly connected.
Collapse
Affiliation(s)
- Angélica S Mata
- Departamento de Física, Universidade Federal de Viçosa, 36570-000 Viçosa, Minas Gerais, Brazil
| | - Silvio C Ferreira
- Departamento de Física, Universidade Federal de Viçosa, 36570-000 Viçosa, Minas Gerais, Brazil
| |
Collapse
|