1
|
Escribano D, Cuesta JA. Free-energy density functional for Strauss's model of transitive networks. Phys Rev E 2022; 106:054305. [PMID: 36559347 DOI: 10.1103/physreve.106.054305] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 10/25/2022] [Indexed: 06/17/2023]
Abstract
Ensemble models of graphs are one of the most important theoretical tools to study complex networks. Among them, exponential random graphs (ERGs) have proven to be very useful in the analysis of social networks. In this paper we develop a technique, borrowed from the statistical mechanics of lattice gases, to solve Strauss's model of transitive networks. This model was introduced long ago as an ERG ensemble for networks with high clustering and exhibits a first-order phase transition above a critical value of the triangle interaction parameter where two different kinds of networks with different densities of links (or, alternatively, different clustering) coexist. Compared to previous mean-field approaches, our method describes accurately even small networks and can be extended beyond Strauss's classical model-e.g., to networks with different types of nodes. This allows us to tackle, for instance, models with node homophily. We provide results for the latter and show that they accurately reproduce the outcome of Monte Carlo simulations.
Collapse
Affiliation(s)
- Diego Escribano
- Grupo Interdisciplinar de Sistemas Complejos (GISC), Departamento de Matemáticas, Universidad Carlos III de Madrid, 28911 Leganés, Madrid, Spain
| | - José A Cuesta
- Grupo Interdisciplinar de Sistemas Complejos (GISC), Departamento de Matemáticas, Universidad Carlos III de Madrid, 28911 Leganés, Madrid, Spain
- Institute for Biocomputation and Physics of Complex Systems (BIFI), University of Zaragoza, 50018 Zaragoza, Spain
| |
Collapse
|
2
|
Polovnikov K, Gorsky A, Nechaev S, Razin SV, Ulianov SV. Non-backtracking walks reveal compartments in sparse chromatin interaction networks. Sci Rep 2020; 10:11398. [PMID: 32647272 PMCID: PMC7347895 DOI: 10.1038/s41598-020-68182-0] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2020] [Accepted: 06/19/2020] [Indexed: 12/31/2022] Open
Abstract
Chromatin communities stabilized by protein machinery play essential role in gene regulation and refine global polymeric folding of the chromatin fiber. However, treatment of these communities in the framework of the classical network theory (stochastic block model, SBM) does not take into account intrinsic linear connectivity of the chromatin loci. Here we propose the polymer block model, paving the way for community detection in polymer networks. On the basis of this new model we modify the non-backtracking flow operator and suggest the first protocol for annotation of compartmental domains in sparse single cell Hi-C matrices. In particular, we prove that our approach corresponds to the maximum entropy principle. The benchmark analyses demonstrates that the spectrum of the polymer non-backtracking operator resolves the true compartmental structure up to the theoretical detectability threshold, while all commonly used operators fail above it. We test various operators on real data and conclude that the sizes of the non-backtracking single cell domains are most close to the sizes of compartments from the population data. Moreover, the found domains clearly segregate in the gene density and correlate with the population compartmental mask, corroborating biological significance of our annotation of the chromatin compartmental domains in single cells Hi-C matrices.
Collapse
Affiliation(s)
- K Polovnikov
- Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA. .,Skolkovo Institute of Science and Technology, Skolkovo, Russia, 143026.
| | - A Gorsky
- Moscow Institute for Physics and Technology, Dolgoprudnyi, Russia.,Institute for Information Transmission Problems of RAS, Moscow, Russia
| | - S Nechaev
- Interdisciplinary Scientific Center Poncelet (UMI 2615 CNRS), Moscow, Russia, 119002.,Lebedev Physical Institute RAS, Moscow, Russia, 119991
| | - S V Razin
- Institute of Gene Biology, Russian Academy of Sciences, Moscow, Russia.,Faculty of Biology, M.V. Lomonosov Moscow State University, Moscow, Russia
| | - S V Ulianov
- Institute of Gene Biology, Russian Academy of Sciences, Moscow, Russia.,Faculty of Biology, M.V. Lomonosov Moscow State University, Moscow, Russia
| |
Collapse
|
3
|
Valba O, Avetisov V, Gorsky A, Nechaev S. Self-isolation or borders closing: What prevents the spread of the epidemic better? Phys Rev E 2020; 102:010401. [PMID: 32794949 DOI: 10.1103/physreve.102.010401] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2020] [Accepted: 06/23/2020] [Indexed: 11/07/2022]
Abstract
Pandemic propagation of COVID-19 motivated us to discuss the impact of the human network clustering on epidemic spreading. Today, there are two clustering mechanisms which prevent of uncontrolled disease propagation in a connected network: an "internal" clustering, which mimics self-isolation (SI) in local naturally arranged communities, and an "external" clustering, which looks like a sharp frontiers closing (FC) between cities and countries, and which does not care about the natural connections of network agents. SI networks are "evolutionarily grown" under the condition of maximization of small cliques in the entire network, while FC networks are instantly created. Running the standard SIR model on clustered SI and FC networks, we demonstrate that the evolutionary grown clustered network prevents the spread of an epidemic better than the instantly clustered network with similar parameters. We find that SI networks have the scale-free property for the degree distribution P(k)∼k^{η}, with a small critical exponent -2<η<-1. We argue that the scale-free behavior emerges as a result of the randomness in the initial degree distributions.
Collapse
Affiliation(s)
- O Valba
- Department of Applied Mathematics, National Research University Higher School of Economics, 101000 Moscow, Russia.,Federal Research Center of Chemical Physics RAS, 119991 Moscow, Russia
| | - V Avetisov
- Federal Research Center of Chemical Physics RAS, 119991 Moscow, Russia
| | - A Gorsky
- Institute of Information Transmission Problems RAS, 127051 Moscow, Russia.,Moscow Institute of Physics and Technology, Dolgoprudny 141700, Russia
| | - S Nechaev
- Interdisciplinary Scientific Center Poncelet, CNRS UMI 2615, 119002 Moscow, Russia.,P.N. Lebedev Physical Institute RAS, 119991 Moscow, Russia
| |
Collapse
|
4
|
Guzmán-González E, Castillo IP, Metz FL. Phase transitions in atypical systems induced by a condensation transition on graphs. Phys Rev E 2020; 101:012133. [PMID: 32069563 DOI: 10.1103/physreve.101.012133] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2019] [Indexed: 06/10/2023]
Affiliation(s)
- Edgar Guzmán-González
- Department of Quantum Physics and Photonics, Institute of Physics, UNAM, P.O. Box 20-364, 01000 México City, México and London Mathematical Laboratory, 18 Margravine Gardens, London W6 8RH, United Kingdom
| | - Isaac Pérez Castillo
- Department of Quantum Physics and Photonics, Institute of Physics, UNAM, P.O. Box 20-364, 01000 México City, México and London Mathematical Laboratory, 18 Margravine Gardens, London W6 8RH, United Kingdom
| | - Fernando L Metz
- Physics Institute, Federal University of Rio Grande do Sul, 91501-970 Porto Alegre, Brazil and London Mathematical Laboratory, 18 Margravine Gardens, London W6 8RH, United Kingdom
| |
Collapse
|
5
|
Pospelov N, Nechaev S, Anokhin K, Valba O, Avetisov V, Gorsky A. Spectral peculiarity and criticality of a human connectome. Phys Life Rev 2019; 31:240-256. [PMID: 31353222 DOI: 10.1016/j.plrev.2019.07.003] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2019] [Accepted: 07/06/2019] [Indexed: 12/12/2022]
Abstract
We have performed the comparative spectral analysis of structural connectomes for various organisms using open-access data. Our results indicate new peculiar features of connectomes of higher organisms. We found that the spectral density of adjacency matrices of human connectome has maximal deviation from the one of randomized network, compared to other organisms. Considering the network evolution induced by the preference of 3-cycles formation, we discovered that for macaque and human connectomes the evolution with the conservation of local clusterization is crucial, while for primitive organisms the conservation of averaged clusterization is sufficient. Investigating for the first time the level spacing distribution of the spectrum of human connectome Laplacian matrix, we explicitly demonstrate that the spectral statistics corresponds to the critical regime, which is hybrid of Wigner-Dyson and Poisson distributions. This observation provides strong support for debated statement of the brain criticality.
Collapse
Affiliation(s)
- N Pospelov
- Lomonosov Moscow State University, 119991, Moscow, Russia
| | - S Nechaev
- Interdisciplinary Scientific Center Poncelet (CNRS UMI 2615), 119002 Moscow, Russia; P.N. Lebedev Physical Institute RAS, Moscow, Russia.
| | - K Anokhin
- Lomonosov Moscow State University, 119991, Moscow, Russia; National Research Center "Kurchatov Institute", 123098, Moscow, Russia
| | - O Valba
- N.N. Semenov Institute of Chemical Physics RAS, 119991 Moscow, Russia; Department of Applied Mathematics, National Research University Higher School of Economics, 101000 Moscow, Russia
| | - V Avetisov
- N.N. Semenov Institute of Chemical Physics RAS, 119991 Moscow, Russia
| | - A Gorsky
- Institute for Information Transmission Problems RAS, 127051 Moscow, Russia; Moscow Institute of Physics and Technology, Dolgoprudny, 141700 Russia
| |
Collapse
|
6
|
Avetisov V, Gorsky A, Nechaev S, Valba O. Finite plateau in spectral gap of polychromatic constrained random networks. Phys Rev E 2018; 96:062309. [PMID: 29347386 DOI: 10.1103/physreve.96.062309] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2017] [Indexed: 11/07/2022]
Abstract
We consider critical behavior in the ensemble of polychromatic Erdős-Rényi networks and regular random graphs, where network vertices are painted in different colors. The links can be randomly removed and added to the network subject to the condition of the vertex degree conservation. In these constrained graphs we run the Metropolis procedure, which favors the connected unicolor triads of nodes. Changing the chemical potential, μ, of such triads, for some wide region of μ, we find the formation of a finite plateau in the number of intercolor links, which exactly matches the finite plateau in the network algebraic connectivity (the value of the first nonvanishing eigenvalue of the Laplacian matrix, λ_{2}). We claim that at the plateau the spontaneously broken Z_{2} symmetry is restored by the mechanism of modes collectivization in clusters of different colors. The phenomena of a finite plateau formation holds also for polychromatic networks with M≥2 colors. The behavior of polychromatic networks is analyzed via the spectral properties of their adjacency and Laplacian matrices.
Collapse
Affiliation(s)
- V Avetisov
- N.N. Semenov Institute of Chemical Physics RAS, 119991 Moscow, Russia.,Department of Applied Mathematics, National Research University Higher School of Economics, 101000 Moscow, Russia
| | - A Gorsky
- Institute of Information Transmission Problems RAS, Moscow, Russia.,Moscow Institute of Physics and Technology, Dolgoprudny 141700 Russia
| | - S Nechaev
- Interdisciplinary Scientific Center Poncelet (CNRS UMI 2615), Moscow, Russia.,P.N. Lebedev Physical Institute RAS, 119991 Moscow, Russia
| | - O Valba
- N.N. Semenov Institute of Chemical Physics RAS, 119991 Moscow, Russia.,Department of Applied Mathematics, National Research University Higher School of Economics, 101000 Moscow, Russia
| |
Collapse
|