Source author record

Sergey V. Buldyrev

Sergey V. Buldyrev 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

48works
19topics
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

48 published item(s)

preprint2022arXiv

Cascading traffic jamming in a two-dimensional Motter and Lai model

We study the cascading traffic jamming on a two-dimensional random geometric graph using the Motter and Lai model. The traffic jam is caused by a localized attack incapacitating circular region or a line of a certain size, as well as a dispersed attack on an equal number of randomly selected nodes. We investigate if there is a critical size of the attack above which the network becomes completely jammed due to cascading jamming, and how this critical size depends on the average degree $\langle k\rangle$ of the graph, on the number of nodes $N$ in the system, and the tolerance parameter $α$ of the Motter and Lai model.

preprint2022arXiv

Formation of Dissipative Structures in Microscopic Models of Mixtures with Species Interconversion

The separation of substances into different phases is ubiquitous in nature and important scientifically and technologically. This phenomenon may become drastically different if the species involved, whether molecules or supramolecular assemblies, interconvert. In the presence of an external force large enough to overcome energetic differences between the interconvertible species (forced interconversion), the two alternative species will be present in equal amounts, and the striking phenomenon of steady-state, restricted phase separation into mesoscales is observed. Such microphase separation is one of the simplest examples of dissipative structures in condensed matter. In this work, we investigate the formation of such mesoscale steady-state structures through Monte Carlo and Molecular Dynamics simulations of three physically distinct microscopic models of binary mixtures that exhibit both equilibrium (natural) interconversion and a nonequilibrium source of forced interconversion. We show that this source can be introduced through an internal imbalance of intermolecular forces or an external flux of energy that promotes molecular interconversion, possible manifestations of which could include the internal nonequilibrium environment of living cells or a flux of photons. The main trends and observations from the simulations are well captured by a non-equilibrium thermodynamic theory of phase transitions affected by interconversion. We show how a nonequilibrium bicontinuous microemulsion or a spatially modulated state may be generated depending on the interplay between diffusion, natural interconversion, and forced interconversion.

preprint2022arXiv

Modeling Fluid Polyamorphism Through a Maximum-Valence Approach

We suggest a simple model to describe polyamorphism in single-component fluids using a maximum-valence approach. The model contains three types of interactions: i) atoms attract each other by van der Waals forces that generate a liquid-gas transition at low pressures, ii) atoms may form covalent bonds that induce association, and iii) bonded atoms attract or repel each other stronger than non-bonded atoms, thus generating liquid-liquid separation. As an example, we qualitatively compare this model with the behavior of liquid sulfur and show that condition (iii) generates a liquid-liquid phase transition in addition to the liquid-gas phase transition.

preprint2022arXiv

Percolation in heterogeneous spatial networks with long-range interactions

We study the emergence of a giant component in a spatial network where the distribution of the metric distances between the nodes is scale-invariant, and the interaction between the nodes has a long-range power-law behavior. The nodes are positioned in the metric space using a Levy flight procedure, with an associated scale-invariant step probability density function, and is then followed by a process of connecting each pair of nodes with a probability function that depends on the distance between them. A natural way to analyze the system is to consider the total probability for an edge between steps in term of their indexes, by summing over their possible positions. By doing so, a correspondence is found between this model and a model of percolation in a one-dimensional lattice with long-range interactions, which allows the identification of the conditions for which a percolation transition is possible. We find that the emergence of a giant component and percolation transitions is determined by a complicated phase diagram, that exhibits a transition from weak long-range interactions to strong long-range interactions.

preprint2020arXiv

Cascading Failures in Complex Networks

Cascading failure is a potentially devastating process that spreads on real-world complex networks and can impact the integrity of wide-ranging infrastructures, natural systems, and societal cohesiveness. One of the essential features that create complex network vulnerability to failure propagation is the dependency among their components, exposing entire systems to significant risks from destabilizing hazards such as human attacks, natural disasters or internal breakdowns. Developing realistic models for cascading failures as well as strategies to halt and mitigate the failure propagation can point to new approaches to restoring and strengthening real-world networks. In this review, we summarize recent progress on models developed based on physics and complex network science to understand the mechanisms, dynamics and overall impact of cascading failures. We present models for cascading failures in single networks and interdependent networks and explain how different dynamic propagation mechanisms can lead to an abrupt collapse and a rich dynamic behavior. Finally, we close the review with novel emerging strategies for containing cascades of failures and discuss open questions that remain to be addressed.

preprint2020arXiv

Efficient network immunization under limited knowledge

Targeted immunization or attacks of large-scale networks has attracted significant attention by the scientific community. However, in real-world scenarios, knowledge and observations of the network may be limited thereby precluding a full assessment of the optimal nodes to immunize (or remove) in order to avoid epidemic spreading such as that of current COVID-19 epidemic. Here, we study a novel immunization strategy where only $n$ nodes are observed at a time and the most central between these $n$ nodes is immunized (or attacked). This process is continued repeatedly until $1-p$ fraction of nodes are immunized (or attacked). We develop an analytical framework for this approach and determine the critical percolation threshold $p_c$ and the size of the giant component $P_{\infty}$ for networks with arbitrary degree distributions $P(k)$. In the limit of $n\to\infty$ we recover prior work on targeted attack, whereas for $n=1$ we recover the known case of random failure. Between these two extremes, we observe that as $n$ increases, $p_c$ increases quickly towards its optimal value under targeted immunization (attack) with complete information. In particular, we find a new scaling relationship between $|p_c(\infty)-p_c(n)|$ and $n$ as $|p_c(\infty)-p_c(n)|\sim n^{-1}\exp(-αn)$. For Scale-free (SF) networks, where $P(k)\sim k^{-γ}, 2<γ<3$, we find that $p_c$ has a transition from zero to non-zero when $n$ increases from $n=1$ to order of $\log N$ ($N$ is the size of network). Thus, for SF networks, knowledge of order of $\log N$ nodes and immunizing them can reduce dramatically an epidemics.

preprint2020arXiv

Two transitions in spatial modular networks

Understanding the resilience of infrastructures such as transportation network has significant importance for our daily life. Recently, a homogeneous spatial network model was developed for studying spatial embedded networks with characteristic link length such as power-grids and the brain. However, although many real-world networks are spatially embedded and their links have characteristics length such as pipelines, power lines or ground transportation lines they are not homogeneous but rather heterogeneous. For example, density of links within cities are significantly higher than between cities. Here we present and study numerically and analytically a similar realistic heterogeneous spatial modular model using percolation process to better understand the effect of heterogeneity on such networks. The model assumes that inside a city there are many lines connecting different locations, while long lines between the cities are sparse and usually directly connecting only a few nearest neighbours cities in a two dimensional plane. We find that this model experiences two distinct continues transitions, one when the cities disconnect from each other and the second when each city breaks apart. Although the critical threshold for site percolation in 2D grid remains an open question we analytically find the critical threshold for site percolation in this model. In addition, while the homogeneous model experience a single transition having a unique phenomenon called \textit{critical stretching} where a geometric crossover from random to spatial structure in different scales found to stretch non-linearly with the characteristic length at criticality. Here we show that the heterogeneous model does not experience such a phenomenon indicating that critical stretching strongly depends on the network structure.

preprint2019arXiv

Faster calculation of the percolation correlation length on spatial networks

The divergence of the correlation length $ξ$ at criticality is an important phenomenon of percolation in two-dimensional systems. Substantial speed-ups to the calculation of the percolation threshold and component distribution have been achieved by utilizing disjoint sets, but existing algorithms of this sort cannot measure the correlation length. Here, we utilize the parallel axis theorem to track the correlation length as nodes are added to the system, allowing us to utilize disjoint sets to measure $ξ$ for the entire percolation process with arbitrary precision in a single sweep. This algorithm enables direct measurement of the correlation length in lattices as well as spatial network topologies, and provides an important tool for understanding critical phenomena in spatial systems.

preprint2016arXiv

A study of cascading failures in real and synthetic power grid topologies using DC power flows

Using the linearized DC power flow model, we study cascading failures and their spatial and temporal properties in the US Western Interconnect (USWI) power grid. We also introduce the preferential Degree And Distance Attachment (DADA) model, with similar degree distributions, resistances, and currents to the USWI. We investigate the behavior of both grids resulting from the failure of a single line. We find that the DADA model and the USWI model react very similarly to that failure, and that their blackout characteristics resemble each other. In many cases, the failure of a single line can cause cascading failures, which impact the entire grid. We characterize the resilience of the grid by three parameters, the most important of which is tolerance $α$, which is the ratio of the maximal load a line can carry to its initial load. We characterize a blackout by its yield, which we define as the ratio of the final to the initial consumed currents. We find that if $α\leq2$, the probability of a large blackout occurring is very small. By contrast, in a broad range of $1 < α < 2$, the initial failure of a single line can result, with a high probability, in cascading failures leading to a massive blackout with final yield less than 80%. The yield has a bimodal distribution typical of a first-order transition, i.e., the failure of a randomly selected line leads either to an insignificant current reduction or to a major blackout. We find that there is a latent period in the development of major blackouts during which few lines are overloaded, and the yield remains high. The duration of this latent period is proportional to the tolerance. The existence of the latent period suggests that intervention during early time steps of a cascade can significantly reduce the risk of a major blackout.

preprint2016arXiv

Interdependent Lattice Networks in High Dimensions

We study the mutual percolation of two interdependent lattice networks ranging from two to seven dimensions, denoted as $D$. We impose that the length of interdependent links connecting nodes in the two lattices be less than or equal to a certain value, $r$. For each value of $D$ and $r$, we find the mutual percolation threshold, $p_c[D,r]$ below which the system completely collapses through a cascade of failures following an initial destruction of a fraction $ (1-p)$ of the nodes in one of the lattices. We find that for each dimension, $D<6$, there is a value of $r=r_I>1$ such that for $r\geq r_I$ the cascading failures occur as a discontinuous first order transition, while for $r<r_I$ the system undergoes a continuous second order transition, as in the classical percolation theory. Remarkably, for $D=6$, $r_I=1$ which is the same as in random regular (RR) graphs with the same degree (coordination number) of nodes. We also find that in all dimensions, the interdependent lattices reach maximal vulnerability (maximal $p_c[D,r]$) at a distance $r=r_{max}>r_I$, and for $r>r_{max}$ the vulnerability starts to decrease as $r\to\infty$. However the decrease becomes less significant as $D$ increases and $p_c[D,r_{max}]-p_c[D,\infty]$ decreases exponentially with $D$. We also investigate the dependence of $p_c[D,r]$ on the system size as well as how the nature of the transition changes as the number of lattice sites, $N\to\infty$.

preprint2014arXiv

Classical mechanics of economic networks

Financial networks are dynamic. To assess their systemic importance to the world-wide economic network and avert losses we need models that take the time variations of the links and nodes into account. Using the methodology of classical mechanics and Laplacian determinism we develop a model that can predict the response of the financial network to a shock. We also propose a way of measuring the systemic importance of the banks, which we call BankRank. Using European Bank Authority 2011 stress test exposure data, we apply our model to the bipartite network connecting the largest institutional debt holders of the troubled European countries (Greece, Italy, Portugal, Spain, and Ireland). From simulating our model we can determine whether a network is in a "stable" state in which shocks do not cause major losses, or a "unstable" state in which devastating damages occur. Fitting the parameters of the model, which play the role of physical coupling constants, to Eurozone crisis data shows that before the Eurozone crisis the system was mostly in a "stable" regime, and that during the crisis it transitioned into an "unstable" regime. The numerical solutions produced by our model match closely the actual time-line of events of the crisis. We also find that, while the largest holders are usually more important, in the unstable regime smaller holders also exhibit systemic importance. Our model also proves useful for determining the vulnerability of banks and assets to shocks. This suggests that our model may be a useful tool for simulating the response dynamics of shared portfolio networks.

preprint2014arXiv

Inverse melting in a two-dimensional off-lattice model

We carry out computer simulations of a simple, two-dimensional off-lattice model that exhibits inverse melting. The monodisperse system comprises core-softened disks interacting through a repulsive square shoulder located inside an attractive square well. By systematically varying the potential parameters, we increase the pressure range over which the liquid freezes to a crystal upon isobaric heating. The effect is largely controlled by the extent of the shoulder. Despite occurring in two dimensions, the melting transition is first order and to a liquid, rather than to a hexatic or quasicrystal phase. We also provide comment on a commonly employed correlation function used to determine the degree of translational ordering in a system.

preprint2014arXiv

Search for a liquid-liquid critical point in models of silica

Previous research has indicated the possible existence of a liquid-liquid critical point (LLCP) in models of silica at high pressure. To clarify this interesting question we run extended molecular dynamics simulations of two different silica models (WAC and BKS) and perform a detailed analysis of the liquid at temperatures much lower than those previously simulated. We find no LLCP in either model within the accessible temperature range, although it is closely approached in the case of the WAC potential near 4000 K and 5 GPa. Comparing our results with those obtained for other tetrahedral liquids, and relating the average Si-O-Si bond angle and liquid density at the model glass temperature to those of the ice-like beta-cristobalite structure, we conclude that the absence of a critical point can be attributed to insufficient "stiffness" in the bond angle. We hypothesize that a modification of the potential function to mildly favor larger average bond angles will generate a LLCP in a temperature range that is accessible to simulation. The tendency to crystallize in these models is extremely weak in the pressure range studied, although this tendency will undoubtedly increase with increasing stiffness.

preprint2014arXiv

The Italian primary school-size distribution and the city-size: a complex nexus

We characterize the statistical law according to which Italian primary school-size distributes. We find that the school-size can be approximated by a log-normal distribution, with a fat lower tail that collects a large number of very small schools. The upper tail of the school-size distribution decreases exponentially and the growth rates are distributed with a Laplace PDF. These distributions are similar to those observed for firms and are consistent with a Bose-Einstein preferential attachment process. The body of the distribution features a bimodal shape suggesting some source of heterogeneity in the school organization that we uncover by an in-depth analysis of the relation between schools-size and city-size. We propose a novel cluster methodology and a new spatial interaction approach among schools which outline the variety of policies implemented in Italy. Different regional policies are also discussed shedding lights on the relation between policy and geographical features.

preprint2013arXiv

Cascading Failures in Networks with Proximate Dependent Nodes

We study the mutual percolation of a system composed of two interdependent random regular networks. We introduce a notion of distance to explore the effects of the proximity of interdependent nodes on the cascade of failures after an initial attack. We find a non-trivial relation between the nature of the transition through which the networks disintegrate and the parameters of the system, which are the degree of the nodes and the maximum distance between interdependent nodes. We explain this relation by solving the problem analytically for the relevant set of cases.

preprint2013arXiv

Finite-size scaling investigation of the liquid-liquid critical point in ST2 water and its stability with respect to crystallization

The liquid-liquid critical point scenario of water hypothesizes the existence of two metastable liquid phases---low-density liquid (LDL) and high-density liquid (HDL)---deep within the supercooled region. The hypothesis originates from computer simulations of the ST2 water model, but the stability of the LDL phase with respect to the crystal is still being debated. We simulate supercooled ST2 water at constant pressure, constant temperature and constant number of molecules N for N<=729 and times up to 1000 ns. We observe clear differences between the two liquids, both structural and dynamical. Using several methods, including finite-size scaling, we confirm the presence of a liquid-liquid phase transition ending in a critical point. We find that the LDL is stable with respect to the crystal in 98% of our runs (we perform 372 runs for LDL or LDL-like states), and in 100% of our runs for the two largest system sizes (N=512 and 729, for which we perform 136 runs for LDL or LDL-like states). In all these runs tiny crystallites grow and then melt within 1000 ns. Only for N<=343 we observe six events (over 236 runs for LDL or LDL-like states) of spontaneous crystallization after crystallites reach an estimated critical size of about 70+/-10 molecules.

preprint2013arXiv

Percolation of a general network of networks

Percolation theory is an approach to study vulnerability of a system. We develop analytical framework and analyze percolation properties of a network composed of interdependent networks (NetONet). Typically, percolation of a single network shows that the damage in the network due to a failure is a continuous function of the fraction of failed nodes. In sharp contrast, in NetONet, due to the cascading failures, the percolation transition may be discontinuous and even a single node failure may lead to abrupt collapse of the system. We demonstrate our general framework for a NetONet composed of $n$ classic Erdős-Rényi (ER) networks, where each network depends on the same number $m$ of other networks, i.e., a random regular network of interdependent ER networks. In contrast to a \emph{treelike} NetONet in which the size of the largest connected cluster (mutual component) depends on $n$, the loops in the RR NetONet cause the largest connected cluster to depend only on $m$. We also analyzed the extremely vulnerable feedback condition of coupling. In the case of ER networks, the NetONet only exhibits two phases, a second order phase transition and collapse, and there is no first phase transition regime unlike the no feedback condition. In the case of NetONet composed of RR networks, there exists a first order phase transition when $q$ is large and second order phase transition when $q$ is small. Our results can help in designing robust interdependent systems.

preprint2012arXiv

Cascading Failures in Interdependent Lattice Networks: The Critical Role of the Length of Dependency Links

We study the cascading failures in a system composed of two interdependent square lattice networks A and B placed on the same Cartesian plane, where each node in network A depends on a node in network B randomly chosen within a certain distance $r$ from the corresponding node in network A and vice versa. Our results suggest that percolation for small $r$ below $r_{\rm max}\approx 8$ (lattice units) is a second-order transition, and for larger $r$ is a first-order transition. For $r<r_{\rm max}$, the critical threshold increases linearly with $r$ from 0.593 at $r=0$ and reaches a maximum, 0.738 for $r=r_{\rm max}$ and then gradually decreases to 0.683 for $r=\infty$. Our analytical considerations are in good agreement with simulations. Our study suggests that interdependent infrastructures embedded in Euclidean space become most vulnerable when the distance between interdependent nodes is in the intermediate range, which is much smaller than the size of the system.

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

Nanoparticle Confinement in Anomalous Liquids

We investigate using molecular dynamics the effect of the structure of nanoconfinement for liquids with water-like anomalies and liquid-liquid phase transition (LLPT). We find that if the confinement is in an ordered matrix of nanoparticles (NPs) the anomalies are preserved, although the LLPT shifts to lower temperatures, higher pressures and higher densities with respect to bulk. On the contrary, if the NPs matrix is disordered, we find a drastically different phase diagram: the LLPT occurs only in a reduced interval of densities and the anomalies are washed out. To understand this effect we calculate the changes in the system at the microscopic level. In all the different confinements considered here we observe a dramatic increase of density of liquid near the confining NPs. In the disordered case the confinement induces larger heterogeneity in the local density, responsible for the weakening of the LLPT and the disappearance of anomalies.

preprint2012arXiv

Phase diagram of a two-dimensional system with anomalous liquid properties

Using Monte Carlo simulation techniques, we calculate the phase diagram for a square shoulder-square well potential in two dimensions that has been previously shown to exhibit liquid anomalies consistent with a metastable liquid-liquid critical point. We consider the liquid, gas and five crystal phases, and find that all the melting lines are first order, despite a small range of metastability. One melting line exhibits a temperature maximum, as well as a pressure maximum that implies inverse melting over a small range in pressure.

preprint2012arXiv

Physics of the Jagla Model as the Liquid-Liquid Coexistence Line Approaches Horizontal

The slope of the coexistence line of the liquid-liquid phase transition (LLPT) can be positive, negative, or zero. All three possibilities have been found in Monte-Carlo simulations of a modified spherically symmetric two-scale Jagla model. Since the liquid-liquid critical point (LLCP) frequently lies in a region of the phase diagram that is difficult to access experimentally, it is of great interest to study critical phenomena in the supercritical region. We therefore study the properties of the Widom line, which is defined in the one-phase region above the critical point as the locus of maximum correlation length as function of the ordering field at constant thermal field. Asymptotically close to the critical point, the Widom line coincides with the loci of the response function extrema, because all response functions can be asymptotically expressed as functions of the diverging correlation length. We find that the method of identifying the Widom line as the loci of heat capacity maxima becomes unfruitful when the slope of the coexistence line approaches zero in the $T$-$P$ plane. In this case the specific heat displays no maximum in the one-phase region because for a horizontal phase coexistence line, according to the Clapeyron equation, the enthalpy difference between the coexisting phases is zero, and thus there can be no contribution to enthalpy fluctuations from the critical fluctuations. The extension of the coexistence line beyond the critical point into the one-phase region must in this case be performed using density fluctuations; the line of compressibility maxima remains well defined, though it bifurcates into a symmetrical pair of lines. These findings agree well with the linear scaling theory of the LLCP by Anisimov and collaborators.

preprint2012arXiv

Temperature and Length Scale Dependence of Solvophobic Solvation in a Single-site Water-like Liquid

The temperature and length scale dependence of solvation properties of spherical hard solvophobic solutes is investigated in the Jagla liquid, a simple liquid that consists of particles interacting via a spherically symmetric potential combining a hard core repulsion and a longer ranged soft core interaction, yet exhibits water-like anomalies. The results are compared with equivalent calculations for a model of a typical atomic liquid, the Lennard-Jones (LJ) potential, and with predictions for hydrophobic solvation in water using the cavity equation of state and the extended simple point charge (SPC/E) model. We find that the Jagla liquid captures the qualitative thermodynamic behavior of hydrophobic hydration as a function of temperature for both small and large length scale solutes. In particular, for both the Jagla liquid and water, we observe temperature-dependent enthalpy and entropy of solvation for all solute sizes as well as a negative solvation entropy for sufficiently small solutes at low temperature. The results suggest that, compared to a simple liquid, it is the presence of a second thermally accessible repulsive energy scale, acting to increasingly favor larger separations for decreasing temperature, that is the essential characteristic of a liquid that favors low-density, open structures and models hydrophobic hydration, and that it is the presence of this second energy scale that leads to the similarity in the behavior of water and the Jagla liquid. The implications of the temperature and length scale dependence of solvation free energies in water-like liquids are explored with a simple model for the aggregation of solvophobic solutes. We show how aggregate stability depends upon the size of the aggregate and the size of its constituent solutes, and we relate this dependence to cold-induced destabilization phenomena such as the cold-induced denaturation of proteins.

preprint2012arXiv

The extreme vulnerability of interdependent spatially embedded networks

Recent studies show that in interdependent networks a very small failure in one network may lead to catastrophic consequences. Above a critical fraction of interdependent nodes, even a single node failure can invoke cascading failures that may abruptly fragment the system, while below this "critical dependency" (CD) a failure of few nodes leads only to small damage to the system. So far, the research has been focused on interdependent random networks without space limitations. However, many real systems, such as power grids and the Internet, are not random but are spatially embedded. Here we analytically and numerically analyze the stability of systems consisting of interdependent spatially embedded networks modeled as lattice networks. Surprisingly, we find that in lattice systems, in contrast to non-embedded systems, there is no CD and \textit{any} small fraction of interdependent nodes leads to an abrupt collapse. We show that this extreme vulnerability of very weakly coupled lattices is a consequence of the critical exponent describing the percolation transition of a single lattice. Our results are important for understanding the vulnerabilities and for designing robust interdependent spatial embedded networks.

preprint2012arXiv

The robustness of interdependent clustered networks

It was recently found that cascading failures can cause the abrupt breakdown of a system of interdependent networks. Using the percolation method developed for single clustered networks by Newman [Phys. Rev. Lett. {\bf 103}, 058701 (2009)], we develop an analytical method for studying how clustering within the networks of a system of interdependent networks affects the system's robustness. We find that clustering significantly increases the vulnerability of the system, which is represented by the increased value of the percolation threshold $p_c$ in interdependent networks.

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

Modeling Amphiphilic Solutes in a Jagla Solvent

Methanol is an amphiphilic solute whose aqueous solutions exhibit distinctive physical properties. The volume change upon mixing, for example, is negative across the entire composition range, indicating strong association. We explore the corresponding behavior of a Jagla solvent, which has been previously shown to exhibit many of the anomalous properties of water. We consider two models of an amphiphilic solute: (i) a "dimer" model, which consists of one hydrophobic hard sphere linked to a Jagla particle with a permanent bond, and (ii) a "monomer" model, which is a limiting case of the dimer, formed by concentrically overlapping a hard sphere and a Jagla particle. Using discrete molecular dynamics, we calculate the thermodynamic properties of the resulting solutions. We systematically vary the set of parameters of the dimer and monomer models and find that one can readily reproduce the experimental behavior of the excess volume of the methanolwater system as a function of methanol volume fraction. We compare the pressure and temperature dependence of the excess volume and the excess enthalpy of both models with experimental data on methanol-water solutions and find qualitative agreement in most cases. We also investigate the solute effect on the temperature of maximum density and find that the effect of concentration isorders of magnitude stronger than measured experimentally.

preprint2010arXiv

Cascade of failures in coupled network systems with multiple support-dependent relations

We study, both analytically and numerically, the cascade of failures in two coupled network systems A and B, where multiple support-dependent relations are randomly built between nodes of networks A and B. In our model we assume that each node in one network can function only if it has at least a single support node in the other network. If both networks A and B are Erdős-Rényi networks, A and B, with (i) sizes $N^A$ and $N^B$, (ii) average degrees $a$ and $b$, and (iii) $c^{AB}_0N^B$ support links from network A to B and $c^{BA}_0N^B$ support links from network B to A, we find that under random attack with removal of fractions $(1-R^A)N^A$ and $(1-R^B)N^B$ nodes respectively, the percolating giant components of both networks at the end of the cascading failures, $μ^A_\infty$ and $μ^B_\infty$, are given by the percolation laws $μ^A_\infty = R^A [1-\exp{({-c^{BA}_0μ^B_\infty})}] [1-\exp{({-aμ^A_\infty})}]$ and $μ^B_\infty = R^B [1-\exp{({-c^{AB}_0μ^A_\infty})}] [1-\exp{({-bμ^B_\infty})}]$. In the limit of $c^{BA}_0 \to \infty$ and $c^{AB}_0 \to \infty$, both networks become independent, and the giant components are equivalent to a random attack on a single Erdős-Rényi network. We also test our theory on two coupled scale-free networks, and find good agreement with the simulations.

preprint2010arXiv

Dynamics of the contact between a ruthenium surface with a single nanoasperity and a flat ruthenium surface: Molecular dynamics simulations

We study the dynamics of the contact between a pair of surfaces (with properties designed to mimic ruthenium) via molecular dynamics simulations. In particular, we study the contact between a ruthenium surface with a single nanoasperity and a flat ruthenium surface. The results of such simulations suggest that contact behavior is highly variable. The goal of this study is to investigate the source and degree of this variability. We find that during compression, the behavior of the contact force displacement curves is reproducible, while during contact separation, the behavior is highly variable. Examination of the contact surfaces suggest that two separation mechanism are in operation and give rise to this variability. One mechanism corresponds to the formation of a bridge between the two surfaces that plastically stretches as the surfaces are drawn apart and eventually separates in shear. This leads to a morphology after separation in which there are opposing asperities on the two surfaces. This plastic separation/bridge formation mechanism leads to a large work of separation. The other mechanism is a more brittle-like mode in which a crack propagates across the base of the asperity (slightly below the asperity/substrate junction) leading to most of the asperity on one surface or the other after separation and a slight depression facing this asperity on the opposing surface. This failure mode corresponds to a smaller work of separation. those in which a single mechanism operates. Furthermore, contacts made from materials that exhibit predominantly brittle-like behavior will tend to require lower work of separation than those made from ductile-like contact materials.

preprint2010arXiv

Effect of hydrophobic solutes on the liquid-liquid critical point

Jagla ramp particles, interacting through a ramp potential with two characteristic length scales, are known to show in their bulk phase thermodynamic and dynamic anomalies, similar to what is found in water. Jagla particles also exhibit a line of phase transitions separating a low density liquid phase and a high density liquid phase, terminating in a liquid-liquid critical point in a region of the phase diagram that can be studied by simulations. Employing molecular dynamics computer simulations, we study the thermodynamics and the dynamics of solutions of hard spheres (HS) in a solvent formed by Jagla ramp particles. We consider the cases of HS mole fraction x = 0.10, 0.15 and 0.20, and also the case x = 0.50 (a 1:1 mixture of HS and Jagla particles). We find a liquid-liquid critical point, up to the highest HS mole fraction; its position shifts to higher pressures and lower temperatures upon increasing x. We also find that the diffusion coefficient anomalies appear to be preserved for all the mole fractions studied.

preprint2010arXiv

Interdependent networks with correlated degrees of mutually dependent nodes

We study a problem of failure of two interdependent networks in the case of correlated degrees of mutually dependent nodes. We assume that both networks (A and B) have the same number of nodes $N$ connected by the bidirectional dependency links establishing a one-to-one correspondence between the nodes of the two networks in a such a way that the mutually dependent nodes have the same number of connectivity links, i.e. their degrees coincide. This implies that both networks have the same degree distribution $P(k)$. We call such networks correspondently coupled networks (CCN). We assume that the nodes in each network are randomly connected. We define the mutually connected clusters and the mutual giant component as in earlier works on randomly coupled interdependent networks and assume that only the nodes which belong to the mutual giant component remain functional. We assume that initially a $1-p$ fraction of nodes are randomly removed due to an attack or failure and find analytically, for an arbitrary $P(k)$, the fraction of nodes $μ(p)$ which belong to the mutual giant component. We find that the system undergoes a percolation transition at certain fraction $p=p_c$ which is always smaller than the $p_c$ for randomly coupled networks with the same $P(k)$. We also find that the system undergoes a first order transition at $p_c>0$ if $P(k)$ has a finite second moment. For the case of scale free networks with $2<λ\leq 3$, the transition becomes a second order transition. Moreover, if $λ<3$ we find $p_c=0$ as in percolation of a single network. For $λ=3$ we find an exact analytical expression for $p_c>0$. Finally, we find that the robustness of CCN increases with the broadness of their degree distribution.

preprint2010arXiv

Interdependent networks: Reducing the coupling strength leads to a change from a first to second order percolation transition

We study a system composed from two interdependent networks A and B, where a fraction of the nodes in network A depends on the nodes of network B and a fraction of the nodes in network B depends on the nodes of network A. Due to the coupling between the networks when nodes in one network fail they cause dependent nodes in the other network to also fail. This invokes an iterative cascade of failures in both networks. When a critical fraction of nodes fail the iterative process results in a percolation phase transition that completely fragments both networks. We show both analytically and numerically that reducing the coupling between the networks leads to a change from a first order percolation phase transition to a second order percolation transition at a critical point. The scaling of the percolation order parameter near the critical point is characterized by the critical exponent beta=1.

preprint2010arXiv

Robustness of a Network of Networks

Almost all network research has been focused on the properties of a single network that does not interact and depends on other networks. In reality, many real-world networks interact with other networks. Here we develop an analytical framework for studying interacting networks and present an exact percolation law for a network of $n$ interdependent networks. In particular, we find that for $n$ Erdős-Rényi networks each of average degree $k$, the giant component, $P_{\infty}$, is given by $P_{\infty}=p[1-\exp(-kP_{\infty})]^n$ where $1-p$ is the initial fraction of removed nodes. Our general result coincides for $n=1$ with the known Erdős-Rényi second-order phase transition for a single network. For any $n \geq 2$ cascading failures occur and the transition becomes a first-order percolation transition. The new law for $P_{\infty}$ shows that percolation theory that is extensively studied in physics and mathematics is a limiting case ($n=1$) of a more general general and different percolation law for interdependent networks.

preprint2010arXiv

Robustness of interdependent networks under targeted attack

When an initial failure of nodes occurs in interdependent networks, a cascade of failure between the networks occurs. Earlier studies focused on random initial failures. Here we study the robustness of interdependent networks under targeted attack on high or low degree nodes. We introduce a general technique and show that the {\it targeted-attack} problem in interdependent networks can be mapped to the {\it random-attack} problem in a transformed pair of interdependent networks. We find that when the highly connected nodes are protected and have lower probability to fail, in contrast to single scale free (SF) networks where the percolation threshold $p_c=0$, coupled SF networks are significantly more vulnerable with $p_c$ significantly larger than zero. The result implies that interdependent networks are difficult to defend by strategies such as protecting the high degree nodes that have been found useful to significantly improve robustness of single networks.

preprint2010arXiv

The critical effect of dependency groups on the function of networks

Current network models assume one type of links to define the relations between the network entities. However, many real networks can only be correctly described using two different types of relations. Connectivity links that enable the nodes to function cooperatively as a network and dependency links that bind the failure of one network element to the failure of other network elements. Here we present for the first time an analytical framework for studying the robustness of networks that include both connectivity and dependency links. We show that the synergy between the two types of failures leads to an iterative process of cascading failures that has a devastating effect on the network stability and completely alters the known assumptions regarding the robustness of networks. We present exact analytical results for the dramatic change in the network behavior when introducing dependency links. For a high density of dependency links the network disintegrates in a form of a first order phase transition while for a low density of dependency links the network disintegrates in a second order transition. Moreover, opposed to networks containing only connectivity links where a broader degree distribution results in a more robust network, when both types of links are present a broad degree distribution leads to higher vulnerability.

preprint2009arXiv

Catastrophic cascade of failures in interdependent networks

Many systems, ranging from engineering to medical to societal, can only be properly characterized by multiple interdependent networks whose normal functioning depends on one another. Failure of a fraction of nodes in one network may lead to a failure in another network. This in turn may cause further malfunction of additional nodes in the first network and so on. Such a cascade of failures, triggered by a failure of a small faction of nodes in only one network, may lead to the complete fragmentation of all networks. We introduce a model and an analytical framework for studying interdependent networks. We obtain interesting and surprising results that should significantly effect the design of robust real-world networks. For two interdependent Erdos-Renyi (ER) networks, we find that the critical average degree below which both networks collapse is <k_c>=2.445, compared to <k_c>=1 for a single ER network. Furthermore, while for a single network a broader degree distribution of the network nodes results in higher robustness to random failure, for interdependent networks, the broader the distribution is, the more vulnerable the networks become to random failure.

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.

preprint2009arXiv

Structure of shells in complex networks

In a network, we define shell $\ell$ as the set of nodes at distance $\ell$ with respect to a given node and define $r_\ell$ as the fraction of nodes outside shell $\ell$. In a transport process, information or disease usually diffuses from a random node and reach nodes shell after shell. Thus, understanding the shell structure is crucial for the study of the transport property of networks. For a randomly connected network with given degree distribution, we derive analytically the degree distribution and average degree of the nodes residing outside shell $\ell$ as a function of $r_\ell$. Further, we find that $r_\ell$ follows an iterative functional form $r_\ell=ϕ(r_{\ell-1})$, where $ϕ$ is expressed in terms of the generating function of the original degree distribution of the network. Our results can explain the power-law distribution of the number of nodes $B_\ell$ found in shells with $\ell$ larger than the network diameter $d$, which is the average distance between all pairs of nodes. For real world networks the theoretical prediction of $r_\ell$ deviates from the empirical $r_\ell$. We introduce a network correlation function $c(r_\ell)\equiv r_{\ell+1}/ϕ(r_\ell)$ to characterize the correlations in the network, where $r_{\ell+1}$ is the empirical value and $ϕ(r_\ell)$ is the theoretical prediction. $c(r_\ell)=1$ indicates perfect agreement between empirical results and theory. We apply $c(r_\ell)$ to several model and real world networks. We find that the networks fall into two distinct classes: (i) a class of {\it poorly-connected} networks with $c(r_\ell)>1$, which have larger average distances compared with randomly connected networks with the same degree distributions; and (ii) a class of {\it well-connected} networks with $c(r_\ell)<1$.

preprint2008arXiv

Fractal Boundaries of Complex Networks

We introduce the concept of boundaries of a complex network as the set of nodes at distance larger than the mean distance from a given node in the network. We study the statistical properties of the boundaries nodes of complex networks. We find that for both Erdös-Rényi and scale-free model networks, as well as for several real networks, the boundaries have fractal properties. In particular, the number of boundaries nodes {\it B} follows a power-law probability density function which scales as $B^{-2}$. The clusters formed by the boundary nodes are fractals with a fractal dimension $d_{f} \approx 2$. We present analytical and numerical evidence supporting these results for a broad class of networks. Our findings imply potential applications for epidemic spreading.

preprint2006arXiv

Anomalous electrical and frictionless flow conductance in complex networks

We study transport properties such as electrical and frictionless flow conductance on scale-free and Erdos-Renyi networks. We consider the conductance G between two arbitrarily chosen nodes where each link has the same unit resistance. Our theoretical analysis for scale-free networks predicts a broad range of values of G, with a power-law tail distribution Φ_{SF}(G) \sim G^{g_G}, where g_G = 2λ- 1, where λis the decay exponent for the scale-free network degree distribution. We confirm our predictions by simulations of scale-free networks solving the Kirchhoff equations for the conductance between a pair of nodes. The power-law tail in Φ_{SF}(G) leads to large values of G, thereby significantly improving the transport in scale-free networks, compared to Erdos-Renyi networks where the tail of the conductivity distribution decays exponentially. Based on a simple physical 'transport backbone' picture we suggest that the conductances of scale-free and Erdos-Renyi networks can be approximated by ck_Ak_B/(k_A+k_B) for any pair of nodes A and B with degrees k_A and k_B. Thus, a single quantity c, which depends on the average degree <k> of the network, characterizes transport on both scale-free and Erdos-Renyi networks. We determine that c tends to 1 for increasing <k>, and it is larger for scale-free networks. We compare the electrical results with a model for frictionless transport, where conductance is defined as the number of link-independent paths between A and B, and find that a similar picture holds. The effects of distance on the value of conductance are considered for both models, and some differences emerge. Finally, we use a recent data set for the AS (autonomous system) level of the Internet and confirm that our results are valid in this real-world example.

preprint2005arXiv

Current Flow in Random Resistor Networks: The Role of Percolation in Weak and Strong Disorder

We study the current flow paths between two edges in a random resistor network on a $L\times L$ square lattice. Each resistor has resistance $e^{ax}$, where $x$ is a uniformly-distributed random variable and $a$ controls the broadness of the distribution. We find (a) the scaled variable $u\equiv L/a^ν$, where $ν$ is the percolation connectedness exponent, fully determines the distribution of the current path length $\ell$ for all values of $u$. For $u\gg 1$, the behavior corresponds to the weak disorder limit and $\ell$ scales as $\ell\sim L$, while for $u\ll 1$, the behavior corresponds to the strong disorder limit with $\ell\sim L^{d_{\scriptsize opt}}$, where $d_{\scriptsize opt} = 1.22\pm0.01$ is the optimal path exponent. (b) In the weak disorder regime, there is a length scale $ξ\sim a^ν$, below which strong disorder and critical percolation characterize the current path.

preprint2005arXiv

Possible Connection between the Optimal Path and Flow in Percolation Clusters

We study the behavior of the optimal path between two sites separated by a distance $r$ on a $d$-dimensional lattice of linear size $L$ with weight assigned to each site. We focus on the strong disorder limit, i.e., when the weight of a single site dominates the sum of the weights along each path. We calculate the probability distribution $P(\ell_{\rm opt}|r,L)$ of the optimal path length $\ell_{\rm opt}$, and find for $r\ll L$ a power law decay with $\ell_{\rm opt}$, characterized by exponent $g_{\rm opt}$. We determine the scaling form of $P(\ell_{\rm opt}|r,L)$ in two- and three-dimensional lattices. To test the conjecture that the optimal paths in strong disorder and flow in percolation clusters belong to the same universality class, we study the tracer path length $\ell_{\rm tr}$ of tracers inside percolation through their probability distribution $P(\ell_{\rm tr}|r,L)$. We find that, because the optimal path is not constrained to belong to a percolation cluster, the two problems are different. However, by constraining the optimal paths to remain inside the percolation clusters in analogy to tracers in percolation, the two problems exhibit similar scaling properties.

preprint2005arXiv

Scale-Free Networks Emerging from Weighted Random Graphs

We study Erdös-Rényi random graphs with random weights associated with each link. We generate a new ``Supernode network'' by merging all nodes connected by links having weights below the percolation threshold (percolation clusters) into a single node. We show that this network is scale-free, i.e., the degree distribution is $P(k)\sim k^{-λ}$ with $λ=2.5$. Our results imply that the minimum spanning tree (MST) in random graphs is composed of percolation clusters, which are interconnected by a set of links that create a scale-free tree with $λ=2.5$. We show that optimization causes the percolation threshold to emerge spontaneously, thus creating naturally a scale-free ``supernode network''. We discuss the possibility that this phenomenon is related to the evolution of several real world scale-free networks.

preprint2004arXiv

Anomalous Transport in Complex Networks

To study transport properties of complex networks, we analyze the equivalent conductance $G$ between two arbitrarily chosen nodes of random scale-free networks with degree distribution $P(k)\sim k^{-λ}$ in which each link has the same unit resistance. We predict a broad range of values of $G$, with a power-law tail distribution $Φ_{\rm SF}(G)\sim G^{-g_G}$, where $g_G=2λ-1$, and confirm our predictions by simulations. The power-law tail in $Φ_{\rm SF}(G)$ leads to large values of $G$, thereby significantly improving the transport in scale-free networks, compared to Erdős-Rényi random graphs where the tail of the conductivity distribution decays exponentially. Based on a simple physical ``transport backbone'' picture we show that the conductances are well approximated by $ck_Ak_B/(k_A+k_B)$ for any pair of nodes $A$ and $B$ with degrees $k_A$ and $k_B$. Thus, a single parameter $c$ characterizes transport on scale-free networks.

preprint2004arXiv

Folding Trp-cage to NMR resolution native structure using a coarse-grained model

We develop a coarse-grained protein model with a simplified amino acid interaction potential. We perform discrete molecular dynamics folding simulations of a small 20 residue protein - Trp-cage - from a fully extended conformation. We demonstrate the ability of the Trp-cage model to consistently reach conformations within 2angstrom backbone root-mean-square distance (RMSD) from the corresponding NMR structures. The minimum RMSD of Trp-cage conformations in the simulation can be smaller than 1.00angstrom. Our findings suggest that, at least for the case of Trp-cage, a detailed all-atom protein model with a physical molecular mechanics force field is not necessary to reach the native state of a protein. Our results also suggest that the success folding Trp-cage in our simulations and in the reported all-atom molecular mechanics simulations studies may be mainly due to the special stabilizing features specific to this miniprotein.

preprint2002arXiv

Connection between Adam-Gibbs Theory and Spatially Heterogeneous Dynamics

We investigate the spatially heterogeneous dynamics in the SPC/E model of water by using molecular dynamics simulations. We relate the average mass $n^*$ of mobile particle clusters to the diffusion constant and the configurational entropy. Hence, $n^*$ can be interpreted as the mass of the ``cooperatively rearranging regions'' that form the basis of the Adam-Gibbs theory of the dynamics of supercooled liquids. Finally, we examine the time and temperature dependence of these transient clusters.