Source author record

Hernan A. Makse

Hernan A. Makse appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

43works
18topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

43 published item(s)

preprint2022arXiv

Radiologist-level Performance by Using Deep Learning for Segmentation of Breast Cancers on MRI Scans

Purpose: To develop a deep network architecture that would achieve fully automated radiologist-level segmentation of cancers at breast MRI. Materials and Methods: In this retrospective study, 38229 examinations (composed of 64063 individual breast scans from 14475 patients) were performed in female patients (age range, 12-94 years; mean age, 52 years +/- 10 [standard deviation]) who presented between 2002 and 2014 at a single clinical site. A total of 2555 breast cancers were selected that had been segmented on two-dimensional (2D) images by radiologists, as well as 60108 benign breasts that served as examples of noncancerous tissue; all these were used for model training. For testing, an additional 250 breast cancers were segmented independently on 2D images by four radiologists. Authors selected among several three-dimensional (3D) deep convolutional neural network architectures, input modalities, and harmonization methods. The outcome measure was the Dice score for 2D segmentation, which was compared between the network and radiologists by using the Wilcoxon signed rank test and the two one-sided test procedure. Results: The highest-performing network on the training set was a 3D U-Net with dynamic contrast-enhanced MRI as input and with intensity normalized for each examination. In the test set, the median Dice score of this network was 0.77 (interquartile range, 0.26). The performance of the network was equivalent to that of the radiologists (two one-sided test procedures with radiologist performance of 0.69-0.84 as equivalence bounds, P <= .001 for both; n = 250). Conclusion: When trained on a sufficiently large dataset, the developed 3D U-Net performed as well as fellowship-trained radiologists in detailed 2D segmentation of breast cancers at routine clinical MRI.

preprint2021arXiv

Fast algorithm to identify cluster synchrony through fibration symmetries in large information-processing networks

Recent studies revealed an important interplay between the detailed structure of fibration symmetric circuits and the functionality of biological and non-biological networks within which they have be identified. The presence of these circuits in complex networks are directed related to the phenomenon of cluster synchronization, which produces patterns of synchronized group of nodes. Here we present a fast, and memory efficient, algorithm to identify fibration symmetries over information-processing networks. This algorithm is specially suitable for large and sparse networks since it has runtime of complexity $O(M\log N)$ and requires $O(M+N)$ of memory resources, where $N$ and $M$ are the number of nodes and edges in the network, respectively. We propose a modification on the so-called refinement paradigm to identify circuits symmetrical to information flow (i.e., fibers) by finding the coarsest refinement partition over the network. Finally, we show that the presented algorithm provides an optimal procedure for identifying fibers, overcoming the current approaches used in the literature.

preprint2021arXiv

Why polls fail to predict elections

In the past decade we have witnessed the failure of traditional polls in predicting presidential election outcomes across the world. To understand the reasons behind these failures we analyze the raw data of a trusted pollster which failed to predict, along with the rest of the pollsters, the surprising 2019 presidential election in Argentina which has led to a major market collapse in that country. Analysis of the raw and re-weighted data from longitudinal surveys performed before and after the elections reveals clear biases (beyond well-known low-response rates) related to mis-representation of the population and, most importantly, to social-desirability biases, i.e., the tendency of respondents to hide their intention to vote for controversial candidates. We then propose a longitudinal opinion tracking method based on big-data analytics from social media, machine learning, and network theory that overcomes the limits of traditional polls. The model achieves accurate results in the 2019 Argentina elections predicting the overwhelming victory of the candidate Alberto Fernández over the president Mauricio Macri; a result that none of the traditional pollsters in the country was able to predict. Beyond predicting political elections, the framework we propose is more general and can be used to discover trends in society; for instance, what people think about economics, education or climate change.

preprint2020arXiv

Circuits with broken fibration symmetries perform core logic computations in biological networks

We show that logic computational circuits in gene regulatory networks arise from a fibration symmetry breaking in the network structure. From this idea we implement a constructive procedure that reveals a hierarchy of genetic circuits, ubiquitous across species, that are surprising analogues to the emblematic circuits of solid-state electronics: starting from the transistor and progressing to ring oscillators, current-mirror circuits to toggle switches and flip-flops. These canonical variants serve fundamental operations of synchronization and clocks (in their symmetric states) and memory storage (in their broken symmetry states). These conclusions introduce a theoretically principled strategy to search for computational building blocks in biological networks, and present a systematic route to design synthetic biological circuits.

preprint2020arXiv

Fibration symmetries uncover the building blocks of biological networks

A major ambition of systems science is to uncover the building blocks of any biological network to decipher how cellular function emerges from their interactions. Here, we introduce a graph representation of the information flow in these networks as a set of input trees, one for each node, which contains all pathways along which information can be transmitted in the network. In this representation, we find remarkable symmetries in the input trees that deconstruct the network into functional building blocks called fibers. Nodes in a fiber have isomorphic input trees and thus process equivalent dynamics and synchronize their activity. Each fiber can then be collapsed into a single representative base node through an information-preserving transformation called 'symmetry fibration', introduced by Grothendieck in the context of algebraic geometry. We exemplify the symmetry fibrations in gene regulatory networks and then show that they universally apply across species and domains from biology to social and infrastructure networks. The building blocks are classified into topological classes of input trees characterized by integer branching ratios and fractal golden ratios of Fibonacci sequences representing cycles of information. Thus, symmetry fibrations describe how complex networks are built from the bottom up to process information through the synchronization of their constitutive building blocks.

preprint2019arXiv

Symmetry group factorization reveals the structure-function relation in the neural connectome of Caenorhabditis elegans

The neural connectome of the nematode Caenorhabditis elegans has been completely mapped, yet in spite of being one of the smallest connectomes (302 neurons), the design principles that explain how the connectome structure determines its function remain unknown. Here, we find symmetries in the locomotion neural circuit of C. elegans, each characterized by its own symmetry group which can be factorized into the direct product of normal subgroups. The action of these normal subgroups partitions the connectome into sectors of neurons that match broad functional categories. Furthermore, symmetry principles predict the existence of novel finer structures inside these normal subgroups forming feedforward and recurrent networks made of blocks of imprimitivity. These blocks constitute structures made of circulant matrices nested in a hierarchy of block-circulant matrices, whose functionality is understood in terms of neural processing filters responsible for fast processing of information.

preprint2016arXiv

A first-order phase transition at the random close packing of hard spheres

Randomly packing spheres of equal size into a container consistently results in a static configuration with a density of ~64%. The ubiquity of random close packing (RCP) rather than the optimal crystalline array at 74% begs the question of the physical law behind this empirically deduced state. Indeed, there is no signature of any macroscopic quantity with a discontinuity associated with the observed packing limit. Here we show that RCP can be interpreted as a manifestation of a thermodynamic singularity, which defines it as the "freezing point" in a first-order phase transition between ordered and disordered packing phases. Despite the athermal nature of granular matter, we show the thermodynamic character of the transition in that it is accompanied by sharp discontinuities in volume and entropy. This occurs at a critical compactivity, which is the intensive variable that plays the role of temperature in granular matter. Our results predict the experimental conditions necessary for the formation of a jammed crystal by calculating an analogue of the "entropy of fusion". This approach is useful since it maps out-of-equilibrium problems in complex systems onto simpler established frameworks in statistical mechanics.

preprint2016arXiv

Angoricity and compactivity describe the jamming transition in soft particulate matter

The application of concepts from equilibrium statistical mechanics to out-of-equilibrium systems has a long history of describing diverse systems ranging from glasses to granular materials. For dissipative jammed systems-- particulate grains or droplets-- a key concept is to replace the energy ensemble describing conservative systems by the volume-stress ensemble. Here, we test the applicability of the volume-stress ensemble to describe the jamming transition by comparing the jammed configurations obtained by dynamics with those averaged over the ensemble as a probe of ergodicity. Agreement between both methods suggests the idea of "thermalization" at a given angoricity and compactivity. We elucidate the thermodynamic order of the jamming transition by showing the absence of critical fluctuations in static observables like pressure and volume. The approach allows to calculate observables such as the entropy, volume, pressure, coordination number and distribution of forces to characterize the scaling laws near the jamming transition from a statistical mechanics viewpoint.

preprint2016arXiv

Collective behavior in the spatial spreading of obesity

Non-communicable diseases like diabetes, obesity and certain forms of cancer have been increasing in many countries at alarming levels. A difficulty in the conception of policies to reverse these trends is the identification of the drivers behind the global epidemics. Here, we implement a spatial spreading analysis to investigate whether diabetes, obesity and cancer show spatial correlations revealing the effect of collective and global factors acting above individual choices. We adapt a theoretical framework for critical physical systems displaying collective behavior to decipher the laws of spatial spreading of diseases. We find a regularity in the spatial fluctuations of their prevalence revealed by a pattern of scale-free long-range correlations. The fluctuations are anomalous, deviating in a fundamental way from the weaker correlations found in the underlying population distribution. This collective behavior indicates that the spreading dynamics of obesity, diabetes and some forms of cancer like lung cancer are analogous to a critical point of fluctuations, just as a physical system in a second-order phase transition. According to this notion, individual interactions and habits may have negligible influence in shaping the global patterns of spreading. Thus, obesity turns out to be a global problem where local details are of little importance. Interestingly, we find the same critical fluctuations in obesity and diabetes, and in the activities of economic sectors associated with food production such as supermarkets, food and beverage stores--- which cluster in a different universality class than other generic sectors of the economy. These results motivate future interventions to investigate the causality of this relation providing guidance for the implementation of preventive health policies.

preprint2016arXiv

Collective Influence Algorithm to find influencers via optimal percolation in massively large social media

We elaborate on a linear time implementation of the Collective Influence (CI) algorithm introduced by Morone, Makse, Nature 524, 65 (2015) to find the minimal set of influencers in a network via optimal percolation. We show that the computational complexity of CI is O(N log N) when removing nodes one-by-one, with N the number of nodes. This is made possible by using an appropriate data structure to process the CI values, and by the finite radius l of the CI sphere. Furthermore, we introduce a simple extension of CI when l is infinite, the CI propagation (CI_P) algorithm, that considers the global optimization of influence via message passing in the whole network and identifies a slightly smaller fraction of influencers than CI. Remarkably, CI_P is able to reproduce the exact analytical optimal percolation threshold obtained by Bau, Wormald, Random Struct. Alg. 21, 397 (2002) for cubic random regular graphs, leaving little improvement left for random graphs. We also introduce the Collective Immunization Belief Propagation algorithm (CI_BP), a belief-propagation (BP) variant of CI based on optimal immunization, which has the same performance as CI_P. However, this small augmented performance of the order of 1-2 % in the low influencers tail comes at the expense of increasing the computational complexity from O(N log N) to O(N^2 log N), rendering both, CI_P and CI_BP, prohibitive for finding influencers in modern-day big-data. The same nonlinear running time drawback pertains to a recently introduced BP-decimation (BPD) algorithm by Mugisha, Zhou, arXiv:1603.05781. For instance, we show that for big-data social networks of typically 200 million users (eg, active Twitter users sending 500 million tweets per day), CI finds the influencers in less than 3 hours running on a single CPU, while the BP algorithms (CI_P, CI_BP and BDP) would take more than 3,000 years to accomplish the same task.

preprint2016arXiv

Dynamic effective mass of granular media

We develop the concept of frequency dependent effective mass, M(omega), of jammed granular materials which occupy a rigid cavity to a filling fraction of 48%, the remaining volume being air of normal room condition or controlled humidity. The dominant features of M(omega) provide signatures of the dissipation of acoustic modes, elasticity and aging effects in the granular medium. We perform humidity controlled experiments and interpret the data in terms of a continuum model and a "trap" model of thermally activated capillary bridges at the contact points. The results suggest that attenuation in the granular materials is influenced significantly by the kinetics of capillary condensation between the asperities at the contacts.

preprint2016arXiv

Frequency-dependent attenuation and elasticity in unconsolidated earth materials: effect of damping

We use the Discrete Element Method (DEM) to understand the underlying attenuation mechanism in granular media, with special applicability to the measurements of the so-called effective mass developed earlier. We consider that the particles interact via Hertz-Mindlin elastic contact forces and that the damping is describable as a force proportional to the velocity difference of contacting grains. We determine the behavior of the complex-valued normal mode frequencies using 1) DEM, 2) direct diagonalization of the relevant matrix, and 3) a numerical search for the zeros of the relevant determinant. All three methods are in strong agreement with each other. The real and the imaginary parts of each normal mode frequency characterize the elastic and the dissipative properties, respectively, of the granular medium. We demonstrate that, as the interparticle damping, $ξ$, increases, the normal modes exhibit nearly circular trajectories in the complex frequency plane and that for a given value of $ξ$ they all lie on or near a circle of radius $R$ centered on the point $-iR$ in the complex plane, where $R\propto 1/ξ$. We show that each normal mode becomes critically damped at a value of the damping parameter $ξ\approx 1/ω_n^0$, where $ω_n^0$ is the (real-valued) frequency when there is no damping. The strong indication is that these conclusions carry over to the properties of real granular media whose dissipation is dominated by the relative motion of contacting grains. For example, compressional or shear waves in unconsolidated dry sediments can be expected to become overdamped beyond a critical frequency, depending upon the strength of the intergranular damping constant.

preprint2016arXiv

Modularity map of the network of human cell differentiation

Cell differentiation in multicellular organisms is a complex process whose mechanism can be understood by a reductionist approach, in which the individual processes that control the generation of different cell types are identified. Alternatively, a large scale approach in search of different organizational features of the growth stages promises to reveal its modular global structure with the goal of discovering previously unknown relations between cell types. Here we sort and analyze a large set of scattered data to construct the network of human cell differentiation (NHCD) based on cell types (nodes) and differentiation steps (links) from the fertilized egg to a crying baby. We discover a dynamical law of critical branching, which reveals a fractal regularity in the modular organization of the network, and allows us to observe the network at different scales. The emerging picture clearly identifies clusters of cell types following a hierarchical organization, ranging from sub-modules to super-modules of specialized tissues and organs on varying scales. This discovery will allow one to treat the development of a particular cell function in the context of the complex network of human development as a whole. Our results point to an integrated large-scale view of the network of cell types systematically revealing ties between previously unrelated domains in organ functions.

preprint2016arXiv

The Area and Population of Cities: New Insights from a Different Perspective on Cities

The distribution of the population of cities has attracted a great deal of attention, in part because it sharply constrains models of local growth. However, to this day, there is no consensus on the distribution below the very upper tail, because available data need to rely on the "legal" rather than "economic" definition of cities for medium and small cities. To remedy this difficulty, in this work we construct cities "from the bottom up" by clustering populated areas obtained from high-resolution data. This method allows us to investigate the population and area of cities for urban agglomerations of all sizes using clustering methods from percolation theory. We find that Zipf's law (a power law with exponent close to 1) for population holds for cities as small as 12,000 inhabitants in the USA and 5,000 inhabitants in Great Britain. In addition the distribution of city areas is also close to a Zipf's law. We provide a parsimonious model with endogenous city area that is consistent with those findings.

preprint2016arXiv

Theory of random packings

We review a recently proposed theory of random packings. We describe the volume fluctuations in jammed matter through a volume function, amenable to analytical and numerical calculations. We combine an extended statistical mechanics approach 'a la Edwards' (where the role traditionally played by the energy and temperature in thermal systems is substituted by the volume and compactivity) with a constraint on mechanical stability imposed by the isostatic condition. We show how such approaches can bring results that can be compared to experiments and allow for an exploitation of the statistical mechanics framework. The key result is the use of a relation between the local Voronoi volume of the constituent grains and the number of neighbors in contact that permits a simple combination of the two approaches to develop a theory of random packings. We predict the density of random loose packing (RLP) and random close packing (RCP) in close agreement with experiments and develop a phase diagram of jammed matter that provides a unifying view of the disordered hard sphere packing problem and further shedding light on a diverse spectrum of data, including the RLP state. Theoretical results are well reproduced by numerical simulations that confirm the essential role played by friction in determining both the RLP and RCP limits. Finally we present an extended discussion on the existence of geometrical and mechanical coordination numbers and how to measure both quantities in experiments and computer simulations.

preprint2015arXiv

Influence maximization in complex networks through optimal percolation

The whole frame of interconnections in complex networks hinges on a specific set of structural nodes, much smaller than the total size, which, if activated, would cause the spread of information to the whole network [1]; or, if immunized, would prevent the diffusion of a large scale epidemic [2,3]. Localizing this optimal, i.e. minimal, set of structural nodes, called influencers, is one of the most important problems in network science [4,5]. Despite the vast use of heuristic strategies to identify influential spreaders [6-14], the problem remains unsolved. Here, we map the problem onto optimal percolation in random networks to identify the minimal set of influencers, which arises by minimizing the energy of a many-body system, where the form of the interactions is fixed by the non-backtracking matrix [15] of the network. Big data analyses reveal that the set of optimal influencers is much smaller than the one predicted by previous heuristic centralities. Remarkably, a large number of previously neglected weakly-connected nodes emerges among the optimal influencers. These are topologically tagged as low-degree nodes surrounded by hierarchical coronas of hubs, and are uncovered only through the optimal collective interplay of all the influencers in the network. Eventually, the present theoretical framework may hold a larger degree of universality, being applicable to other hard optimization problems exhibiting a continuous transition from a known phase [16].

preprint2014arXiv

Cavity method for force transmission in jammed disordered packings of hard particles

The force distribution of jammed disordered packings has always been considered a central object in the physics of granular materials. However, many of its features are poorly understood. In particular, analytic relations to other key macroscopic properties of jammed matter, such as the contact network and its coordination number, are still lacking. Here we develop a mean-field theory for this problem, based on the consideration of the contact network as a random graph where the force transmission becomes a constraint optimization problem. We can thus use the cavity method developed in the last decades within the statistical physics of spin glasses and hard computer science problems. This method allows us to compute the force distribution $\text P(f)$ for random packings of hard particles of any shape, with or without friction. We find a new signature of jamming in the small force behavior $\text P(f) \sim f^θ$, whose exponent has attracted recent active interest: we find a finite value for $\text P(f=0)$, along with $θ=0$. Furthermore, we relate the force distribution to a lower bound of the average coordination number $\, {\bar z}_{\rm c}^{\rm min}(μ)$ of jammed packings of frictional spheres with coefficient $μ$. This bridges the gap between the two known isostatic limits $\, {\bar z}_{\rm c}(μ=0)=2D$ (in dimension $D$) and $\, {\bar z}_{\rm c}(μ\to \infty)=D+1$ by extending the naive Maxwell's counting argument to frictional spheres. The theoretical framework describes different types of systems, such as non-spherical objects in arbitrary dimensions, providing a common mean-field scenario to investigate force transmission, contact networks and coordination numbers of jammed disordered packings.

preprint2014arXiv

Searching for superspreaders of information in real-world social media

A number of predictors have been suggested to detect the most influential spreaders of information in online social media across various domains such as Twitter or Facebook. In particular, degree, PageRank, k-core and other centralities have been adopted to rank the spreading capability of users in information dissemination media. So far, validation of the proposed predictors has been done by simulating the spreading dynamics rather than following real information flow in social networks. Consequently, only model-dependent contradictory results have been achieved so far for the best predictor. Here, we address this issue directly. We search for influential spreaders by following the real spreading dynamics in a wide range of networks. We find that the widely-used degree and PageRank fail in ranking users' influence. We find that the best spreaders are consistently located in the k-core across dissimilar social platforms such as Twitter, Facebook, Livejournal and scientific publishing in the American Physical Society. Furthermore, when the complete global network structure is unavailable, we find that the sum of the nearest neighbors' degree is a reliable local proxy for user's influence. Our analysis provides practical instructions for optimal design of strategies for "viral" information dissemination in relevant applications.

preprint2014arXiv

Statistical theory of correlations in random packings of hard particles

A random packing of hard particles represents a fundamental model for granular matter. Despite its importance, analytical modeling of random packings remains difficult due to the existence of strong correlations which preclude the development of a simple theory. Here, we take inspiration from liquid theories for the $n$-particle angular correlation function to develop a formalism of random packings of hard particles from the bottom-up. A progressive expansion into a shell of particles converges in the large layer limit under a Kirkwood-like approximation of higher-order correlations. We apply the formalism to hard disks and predict the density of two-dimensional random close packing (RCP), $ϕ_{\rm rcp} = 0.85\pm0.01$, and random loose packing (RLP), $ϕ_{\rm rlp} = 0.67\pm0.01$. Our theory also predicts a phase diagram and angular correlation functions that are in good agreement with experimental and numerical data.

preprint2014arXiv

Structural properties of dense hard spheres near random close packing

We numerically study structural properties of mechanically stable packings of hard spheres (HS), in a wide range of packing fractions $0.53 \le ϕ\le 0.72$. Detailed structural information is obtained from the analysis of orientational order parameters, which clearly reveals a disorder-order phase transition at the random close packing (RCP) density, $ϕ_{\rm c} \simeq 0.64$. Above $ϕ_{\rm c}$ the crystalline nuclei form 3D-like clusters, which upon further desification, transform into alternating planar-like layers. We also find that particles with icosahedral symmetry survive only in a narrow density range in the vicinity of the RCP transition.

preprint2013arXiv

Calculation of the Voronoi boundary for lens-shaped particles and spherocylinders

We have recently developed a mean-field theory to estimate the packing fraction of non-spherical particles [A. Baule et al., Nature Commun. (2013)]. The central quantity in this framework is the Voronoi excluded volume, which generalizes the standard hard-core excluded volume appearing in Onsager's theory. The Voronoi excluded volume is defined from an exclusion condition for the Voronoi boundary between two particles, which is usually not tractable analytically. Here, we show how the technical difficulties in calculating the Voronoi boundary can be overcome for lens-shaped particles and spherocylinders, two standard prolate and oblate shapes with rotational symmetry. By decomposing these shapes into unions and intersections of spheres analytical expressions can be obtained.

preprint2013arXiv

IMDB network revisited: unveiling fractal and modular properties from a typical small-world network

We study a subset of the movie collaboration network, imdb.com, where only adult movies are included. We show that there are many benefits in using such a network, which can serve as a prototype for studying social interactions. We find that the strength of links, i.e., how many times two actors have collaborated with each other, is an important factor that can significantly influence the network topology. We see that when we link all actors in the same movie with each other, the network becomes small-world, lacking a proper modular structure. On the other hand, by imposing a threshold on the minimum number of links two actors should have to be in our studied subset, the network topology becomes naturally fractal. This occurs due to a large number of meaningless links, namely, links connecting actors that did not actually interact. We focus our analysis on the fractal and modular properties of this resulting network, and show that the renormalization group analysis can characterize the self-similar structure of these networks.

preprint2013arXiv

Mean-field theory of random close packings of axisymmetric particles

Finding the optimal random packing of non-spherical particles is an open problem with great significance in a broad range of scientific and engineering fields. So far, this search has been performed only empirically on a case-by-case basis, in particular, for shapes like dimers, spherocylinders and ellipsoids of revolution. Here, we present a mean-field formalism to estimate the packing density of axisymmetric non-spherical particles. We derive an analytic continuation from the sphere that provides a phase diagram predicting that, for the same coordination number, the density of monodisperse random packings follows the sequence of increasing packing fractions: spheres < oblate ellipsoids < prolate ellipsoids < dimers < spherocylinders. We find the maximal packing densities of 73.1% for spherocylinders and 70.7% for dimers, in good agreement with the largest densities found in simulations. Moreover, we find a packing density of 73.6% for lens-shaped particles, representing the densest random packing of the axisymmetric objects studied so far.

preprint2013arXiv

Origins of power-law degree distribution in the heterogeneity of human activity in social networks

The probability distribution of number of ties of an individual in a social network follows a scale-free power-law. However, how this distribution arises has not been conclusively demonstrated in direct analyses of people's actions in social networks. Here, we perform a causal inference analysis and find an underlying cause for this phenomenon. Our analysis indicates that heavy-tailed degree distribution is causally determined by similarly skewed distribution of human activity. Specifically, the degree of an individual is entirely random - following a "maximum entropy attachment" model - except for its mean value which depends deterministically on the volume of the users' activity. This relation cannot be explained by interactive models, like preferential attachment, since the observed actions are not likely to be caused by interactions with other people.

preprint2013arXiv

Spreading dynamics in complex networks

Searching for influential spreaders in complex networks is an issue of great significance for applications across various domains, ranging from the epidemic control, innovation diffusion, viral marketing, social movement to idea propagation. In this paper, we first display some of the most important theoretical models that describe spreading processes, and then discuss the problem of locating both the individual and multiple influential spreaders respectively. Recent approaches in these two topics are presented. For the identification of privileged single spreaders, we summarize several widely used centralities, such as degree, betweenness centrality, PageRank, k-shell, etc. We investigate the empirical diffusion data in a large scale online social community -- LiveJournal. With this extensive dataset, we find that various measures can convey very distinct information of nodes. Of all the users in LiveJournal social network, only a small fraction of them involve in spreading. For the spreading processes in LiveJournal, while degree can locate nodes participating in information diffusion with higher probability, k-shell is more effective in finding nodes with large influence. Our results should provide useful information for designing efficient spreading strategies in reality.

preprint2013arXiv

The evolutionary dynamics of protein-protein interaction networks inferred from the reconstruction of ancient networks

Cellular functions are based on the complex interplay of proteins, therefore the structure and dynamics of these protein-protein interaction (PPI) networks are the key to the functional understanding of cells. In the last years, large-scale PPI networks of several model organisms were investigated. Methodological improvements now allow the analysis of PPI networks of multiple organisms simultaneously as well as the direct modeling of ancestral networks. This provides the opportunity to challenge existing assumptions on network evolution. We utilized present-day PPI networks from integrated datasets of seven model organisms and developed a theoretical and bioinformatic framework for studying the evolutionary dynamics of PPI networks. A novel filtering approach using percolation analysis was developed to remove low confidence interactions based on topological constraints. We then reconstructed the ancient PPI networks of different ancestors, for which the ancestral proteomes, as well as the ancestral interactions, were inferred. Ancestral proteins were reconstructed using orthologous groups on different evolutionary levels. A stochastic approach, using the duplication-divergence model, was developed for estimating the probabilities of ancient interactions from today's PPI networks. The growth rates for nodes, edges, sizes and modularities of the networks indicate multiplicative growth and are consistent with the results from independent static analysis. Our results support the duplication-divergence model of evolution and indicate fractality and multiplicative growth as general properties of the PPI network structure and dynamics.

preprint2012arXiv

A small-world of weak ties provides optimal global integration of self-similar modules in functional brain networks

The human brain is organized in functional modules. Such an organization presents a basic conundrum: modules ought to be sufficiently independent to guarantee functional specialization and sufficiently connected to bind multiple processors for efficient information transfer. It is commonly accepted that small-world architecture of short lengths and large local clustering may solve this problem. However, there is intrinsic tension between shortcuts generating small-worlds and the persistence of modularity; a global property unrelated to local clustering. Here, we present a possible solution to this puzzle. We first show that a modified percolation theory can define a set of hierarchically organized modules made of strong links in functional brain networks. These modules are "large-world" self-similar structures and, therefore, are far from being small-world. However, incorporating weaker ties to the network converts it into a small-world preserving an underlying backbone of well-defined modules. Remarkably, weak ties are precisely organized as predicted by theory maximizing information transfer with minimal wiring cost. This trade-off architecture is reminiscent of the "strength of weak ties" crucial concept of social networks. Such a design suggests a natural solution to the paradox of efficient information flow in the highly modular structure of the brain.

preprint2012arXiv

Communication activity in a social network: relation between long-term correlations and inter-event clustering

The timing patterns of human communication in social networks is not random. On the contrary, communication is dominated by emergent statistical laws such as non-trivial correlations and clustering. Recently, we found long-term correlations in the user's activity in social communities. Here, we extend this work to study collective behavior of the whole community. The goal is to understand the origin of clustering and long-term persistence. At the individual level, we find that the correlations in activity are a byproduct of the clustering expressed in the power-law distribution of inter-event times of single users. On the contrary, the activity of the whole community presents long-term correlations that are a true emergent property of the system, i.e. they are not related to the distribution of inter-event times. This result suggests the existence of collective behavior, possible arising from nontrivial communication patterns through the embedding social network.

preprint2012arXiv

The conundrum of functional brain networks: small-world efficiency or fractal modularity

The human brain has been studied at multiple scales, from neurons, circuits, areas with well defined anatomical and functional boundaries, to large-scale functional networks which mediate coherent cognition. In a recent work, we addressed the problem of the hierarchical organization in the brain through network analysis. Our analysis identified functional brain modules of fractal structure that were inter-connected in a small-world topology. Here, we provide more details on the use of network science tools to elaborate on this behavior. We indicate the importance of using percolation theory to highlight the modular character of the functional brain network. These modules present a fractal, self-similar topology, identified through fractal network methods. When we lower the threshold of correlations to include weaker ties, the network as a whole assumes a small-world character. These weak ties are organized precisely as predicted by theory maximizing information transfer with minimal wiring costs.

preprint2011arXiv

Communication activity in social networks: growth and correlations

We investigate the timing of messages sent in two online communities with respect to growth fluctuations and long-term correlations. We find that the timing of sending and receiving messages comprises pronounced long-term persistence. Considering the activity of the community members as growing entities, i.e. the cumulative number of messages sent (or received) by the individuals, we identify non-trivial scaling in the growth fluctuations which we relate to the long-term correlations. We find a connection between the scaling exponents of the growth and the long-term correlations which is supported by numerical simulations based on peaks over threshold. In addition, we find that the activity on directed links between pairs of members exhibits long-term correlations, indicating that communication activity with the most liked partners may be responsible for the long-term persistence in the timing of messages. Finally, we show that the number of messages, $M$, and the number of communication partners, $K$, of the individual members are correlated following a power-law, $K\sim M^λ$, with exponent $λ\approx 3/4$.

preprint2011arXiv

How people interact in evolving online affiliation networks

The study of human interactions is of central importance for understanding the behavior of individuals, groups and societies. Here, we observe the formation and evolution of networks by monitoring the addition of all new links and we analyze quantitatively the tendencies used to create ties in these evolving online affiliation networks. We first show that an accurate estimation of these probabilistic tendencies can only be achieved by following the time evolution of the network. For example, actions that are attributed to the usual friend of a friend mechanism through a static snapshot of the network are overestimated by a factor of two. A detailed analysis of the dynamic network evolution shows that half of those triangles were generated through other mechanisms, in spite of the characteristic static pattern. We start by characterizing every single link when the tie was established in the network. This allows us to describe the probabilistic tendencies of tie formation and extract sociological conclusions as follows. The tendencies to add new links differ significantly from what we would expect if they were not affected by the individuals' structural position in the network, i.e., from random link formation. We also find significant differences in behavioral traits among individuals according to their degree of activity, gender, age, popularity and other attributes. For instance, in the particular datasets analyzed here, we find that women reciprocate connections three times as much as men and this difference increases with age. Men tend to connect with the most popular people more often than women across all ages. On the other hand, triangular ties tendencies are similar and independent of gender. Our findings can be useful to build models of realistic social network structures and discover the underlying laws that govern establishment of ties in evolving social networks.

preprint2011arXiv

Identification of influential spreaders in complex networks

Networks portray a multitude of interactions through which people meet, ideas are spread, and infectious diseases propagate within a society. Identifying the most efficient "spreaders" in a network is an important step to optimize the use of available resources and ensure the more efficient spread of information. Here we show that, in contrast to common belief, the most influential spreaders in a social network do not correspond to the best connected people or to the most central people (high betweenness centrality). Instead, we find: (i) The most efficient spreaders are those located within the core of the network as identified by the k-shell decomposition analysis. (ii) When multiple spreaders are considered simultaneously, the distance between them becomes the crucial parameter that determines the extend of the spreading. Furthermore, we find that-- in the case of infections that do not confer immunity on recovered individuals-- the infection persists in the high k-shell layers of the network under conditions where hubs may not be able to preserve the infection. Our analysis provides a plausible route for an optimal design of efficient dissemination strategies.

preprint2010arXiv

Distribution of volumes and coordination number in jammed matter: mesoscopic ensemble

We investigate the distribution of the volume and coordination number associated to each particle in a jammed packing of monodisperse hard sphere using the mesoscopic ensemble developed in Nature 453, 606 (2008). Theory predicts an exponential distribution of the orientational volumes for random close packings and random loose packings. A comparison with computer generated packings reveals deviations from the theoretical prediction in the volume distribution, which can be better modeled by a compressed exponential function. On the other hand, the average of the volumes is well reproduced by the theory leading to good predictions of the limiting densities of RCP and RLP. We discuss a more exact theory to capture the volume distribution in its entire range. The available data suggests a plausible order/disorder transition defining random close packings. Finally, we consider an extended ensemble to calculate the coordination number distribution which is shown to be of an exponential and inverse exponential form for coordinations larger and smaller than the average, respectively, in reasonable agreement with the simulated data.

preprint2010arXiv

Irreversible Incremental Behavior in a Granular Material

We test the elasticity of granular aggregates using increments of shear and volume strain in a numerical simulation. We find that the increment in volume strain is almost reversible, but the increment in shear strain is not. The strength of this irreversibility increases as the average number of contacts per particle (the coordination number) decreases. For increments of volume strain, an elastic model that includes both average and fluctuating motions between contacting particles reproduces well the numerical results over the entire range of coordination numbers. For increments of shear strain, the theory and simulations agree quite well for high values of the coordination number.

preprint2010arXiv

Jamming I: A volume function for jammed matter

We introduce a "Hamiltonian"-like function, called the volume function, indispensable to describe the ensemble of jammed matter such as granular materials and emulsions from a geometrical point of view. The volume function represents the available volume of each particle in the jammed systems. At the microscopic level, we show that the volume function is the Voronoi volume associated to each particle and in turn we provide an analytical formula for the Voronoi volume in terms of the contact network, valid for any dimension. We then develop a statistical theory for the probability distribution of the volumes in 3d to calculate an average volume function coarse-grained at a mesoscopic level. The salient result is the discovery of a mesoscopic volume function inversely proportional to the coordination number. Our analysis is the first step toward the calculation of macroscopic observables and equations of state using the statistical mechanics of jammed matter, when supplemented by the condition of mechanical equilibrium of jamming that properly defines jammed matter at the ensemble level.

preprint2010arXiv

Jamming II: Edwards' statistical mechanics of random packings of hard spheres

The problem of finding the most efficient way to pack spheres has an illustrious history, dating back to the crystalline arrays conjectured by Kepler and the random geometries explored by Bernal in the 60's. This problem finds applications spanning from the mathematician's pencil, the processing of granular materials, the jamming and glass transitions, all the way to fruit packing in every grocery. There are presently numerous experiments showing that the loosest way to pack spheres gives a density of ~55% (RLP) while filling all the loose voids results in a maximum density of ~63-64% (RCP). While those values seem robustly true, to this date there is no physical explanation or theoretical prediction for them. Here we show that random packings of monodisperse hard spheres in 3d can pack between the densities 4/(4 + 2 \sqrt 3) or 53.6% and 6/(6 + 2 \sqrt 3) or 63.4%, defining RLP and RCP, respectively. The reason for these limits arises from a statistical picture of jammed states in which the RCP can be interpreted as the ground state of the ensemble of jammed matter with zero compactivity, while the RLP arises in the infinite compactivity limit. We combine an extended statistical mechanics approach 'a la Edwards' (where the role traditionally played by the energy and temperature in thermal systems is substituted by the volume and compactivity) with a constraint on mechanical stability imposed by the isostatic condition. Ultimately, our results lead to a phase diagram that provides a unifying view of the disordered hard sphere packing problem.

preprint2010arXiv

Model of random packings of different size balls

We develop a model to describe the properties of random assemblies of polydisperse hard spheres. We show that the key features to describe the system are (i) the dependence between the free volume of a sphere and the various coordination numbers between the species, and (ii) the dependence of the coordination numbers with the concentration of species; quantities that are calculated analytically. The model predicts the density of random close packing and random loose packing of polydisperse systems for a given distribution of ball size and describes packings for any interparticle friction coefficient. The formalism allows to determine the optimal packing over different distributions and may help to treat packing problems of non-spherical particles which are notoriously difficult to solve.

preprint2009arXiv

Entropy of jammed matter

We investigate the nature of randomness in disordered packings of frictional spheres. We calculate the entropy of 3D packings through the force and volume ensemble of jammed matter, a mesoscopic ensemble and numerical simulations using volume fluctuation analysis and graph theoretical methods. Equations of state are obtained relating entropy, volume fraction and compactivity characterizing the different states of jammed matter. At the mesoscopic level the entropy vanishes at random close packing. The entropy of the jammed system reveals that the random loose packings are more disordered than random close packings, allowing for an unambiguous interpretation of both limits.

preprint2009arXiv

Jamming III: Characterizing Randomness via the Entropy of Jammed Matter

The nature of randomness in disordered packings of frictional and frictionless spheres is investigated using theory and simulations of identical spherical grains. The entropy of the packings is defined through the force and volume ensemble of jammed matter and shown difficult to calculate analytically. A mesoscopic ensemble of isostatic states is then utilized in an effort to predict the entropy through the defnition of a volume function dependent on the coordination number. Equations of state are obtained relating entropy, volume fraction and compactivity characterizing the different states of jammed matter, and elucidating the phase diagram for jammed granular matter. Analytical calculations are compared to numerical simulations using volume fluctuation analysis and graph theoretical methods, with reasonable agreement. The entropy of the jammed system reveals that the random loose packings are more disordered than random close packings, allowing for an unambiguous interpretation of both limits. Ensemble calculations show that the entropy vanishes at random close packing (RCP), while numerical simulations show that a finite entropy remains in the microscopic states at RCP. The notion of a negative compactivity, that explores states with volume fractions below those achievable by existing simulation protocols, is also explored, expanding the equations of state. We discuss possible extensions to the present mesoscopic approach describing packings from RLP to RCP to the ordered branch of the equation of state in an effort to understand the entropy of jammed matter in the full range of densities from RLP to FCC.

preprint2009arXiv

Scaling laws of human interaction activity

Even though people in our contemporary, technological society are depending on communication, our understanding of the underlying laws of human communicational behavior continues to be poorly understood. Here we investigate the communication patterns in two social Internet communities in search of statistical laws in human interaction activity. This research reveals that human communication networks dynamically follow scaling laws that may also explain the observed trends in economic growth. Specifically, we identify a generalized version of Gibrat's law of social activity expressed as a scaling law between the fluctuations in the number of messages sent by members and their level of activity. Gibrat's law has been essential in understanding economic growth patterns, yet without an underlying general principle for its origin. We attribute this scaling law to long-term correlation patterns in human activity, which surprisingly span from days to the entire period of the available data of more than one year. Further, we provide a mathematical framework that relates the generalized version of Gibrat's law to the long-term correlated dynamics, which suggests that the same underlying mechanism could be the source of Gibrat's law in economics, ranging from large firms, research and development expenditures, gross domestic product of countries, to city population growth. These findings are also of importance for designing communication networks and for the understanding of the dynamics of social systems in which communication plays a role, such as economic markets and political systems.

preprint2008arXiv

From force distribution to average coordination number in frictional granular matter

We study the joint probability distribution of normal and tangential frictional forces in jammed granular media, $P_μ(f_t, f_n)$, for various friction coefficient $μ$, especially when $μ= \infty$. A universal scaling law is found to collapse the data for $μ=0$ to $\infty$ demonstrating a link between force distribution $P_μ(f_t, f_n)$ and average coordination number, $z^μ_c$. The results determine $z_c^μ$ for a finite friction coefficient, extending the constraints counting argument of isostatic granular packing to finite frictional packings.

preprint2006arXiv

Dynamic particle tracking reveals the aging temperature of a colloidal glass

Understanding glasses is considered to be one of the most fundamental problems in statistical physics. A theoretical approach to unravel their universal properties is to consider the validity of equilibrium concepts such as temperature and thermalization in these out-of-equilibrium systems. Here we investigate the autocorrelation and response function to monitor the aging of a colloidal glass. At equilibrium, all the observables are stationary while in the out-of-equilibrium glassy state they have an explicit dependence on the age of the system. We find that the transport coefficients scale with the aging-time as a power-law, a signature of the slow relaxation. Nevertheless, our analysis reveals that the glassy system has thermalized at a constant temperature independent of the age and larger than the bath, reflecting the structural rearrangements of cage-dynamics. Furthermore, a universal scaling law is found to describe the global and local fluctuations of the observables.

preprint1998arXiv

Spontaneous Stratification in Granular Mixtures

Granular materials size segregate when exposed to external periodic perturbations such as vibrations. Moreover, mixtures of grains of different sizes spontaneously segregate in the absence of external perturbations: when a mixture is simply poured onto a pile, the large grains are more likely to be found near the base, while the small grains are more likely to be near the top. Here, we report a spontaneous phenomenon arising when we pour a mixture between two vertical plates: the mixture spontaneously stratifies into alternating layers of small and large grains whenever the large grains are rougher than the small grains. In contrast, we find only spontaneous segregation when the large grains are more rounded than the small grains. The stratification is related to the occurrence of avalanches; during each avalanche the grains comprising the avalanche spontaneously stratify into a pair of layers through a "kink" mechanism, with the small grains forming a sublayer underneath the layer of large grains.