BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

H. Eugene Stanley

H. Eugene Stanley appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.soc-phcond-mat.stat-mechphysics.data-anSocial and Information Networksq-fin.STcond-mat.softcond-mat.dis-nnphysics.pop-phq-fin.GNq-fin.RMphysics.geo-phDigital Librariesq-bioq-fin.TRComputation and Languagecond-mat.mtrl-scicond-mat.othercs.CYnlin.CDphysics.med-phApplicationsBiological Physicscond-matcond-mat.str-elCryptography and Securityecon.GNHuman-Computer InteractionInformation RetrievalMachine Learningmath-phmath.COmath.MPNetworking and Internet ArchitectureNeurons and Cognitionnlin.AOphysics.comp-phq-fin.ECquant-ph

Catalog footprint

What is connected

122works
38topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

122 published item(s)

preprint2022arXiv

Emergence of scaling in dockless bike-sharing systems

Fundamental laws of human mobility have been extensively studied, yet we are still lacking a comprehensive understanding of the mobility patterns of sharing conveyances. Since travellers would highly probably no longer possess their own conveyances in the near future, the interplay between travellers and sharing bikes is a central question for developing more sustainable transportation. Dockless bike-sharing systems that record detailed information of every trip provide us a unique opportunity for revealing the hidden patterns behind riding activities. By treating each bike as an individual entity, we reveal that distributions of mobility indicators of bikes are quite different from humans; and mobility patterns are even inconsistent across cities. All above discrepancies can be well explained by a choice model that is characterized by a universal scaling. Our model unveils that instead of choosing among the newest bikes, the distribution of rank values of selected bikes on usage condition manifests a truncated power-law and is quite stable across several cities despite various diversities. Our framework would have broad implications in sharing economy and contribute towards developing a greener, healthier, and more sustainable future city.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multiple-scale integro-differential perturbation method for generic non-Markovian environments

Non-Markovianity may significantly speed up quantum dynamics when the system interacts strongly with an infinite large reservoir, of which the coupling spectrum should be fine-tuned. The potential benefits are evident in many dynamics schemes, especially the continuous-time quantum walk. Difficulty exists, however, in producing closed-form solutions with controllable accuracy against the complexity of memory kernels. Here, we introduce a new multiple-scale perturbation method that works on integro-differential equations for general study of memory effects in dynamical systems. We propose an open-system model in which a continuous-time quantum walk is enclosed in a non-Markovian reservoir, that naturally corresponds to an error correction algorithm scheme. By applying the multiple-scale method we show how emergence of different time scales is related to transition of system dynamics into the non-Markovian regime. We find that up to two long-term modes and two short-term modes exist in regular networks, limited by their intrinsic symmetries. In addition to the effective approximation by our perturbation method on general forms of reservoirs, the speed-up of quantum walks assisted by non-Markovianity is also confirmed, revealing the advantage of reservoir engineering in designing time-sensitive quantum algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Matrix Factorization Model for Hellinger-based Trust Management in Social Internet of Things

The Social Internet of Things (SIoT), integration of the Internet of Things and Social Networks paradigms, has been introduced to build a network of smart nodes that are capable of establishing social links. In order to deal with misbehaving service provider nodes, service requestor nodes must evaluate their trustworthiness levels. In this paper, we propose a novel trust management mechanism in the SIoT to predict the most reliable service providers for each service requestor, which leads to reduce the risk of being exposed to malicious nodes. We model the SIoT with a flexible bipartite graph (containing two sets of nodes: service providers and service requestors), then build a social network among the service requestor nodes, using the Hellinger distance. Afterward, we develop a social trust model using nodes' centrality and similarity measures to extract trust behaviors among the social network nodes. Finally, a matrix factorization technique is designed to extract latent features of SIoT nodes, find trustworthy nodes, and mitigate the data sparsity and cold start problems. We analyze the effect of parameters in the proposed trust prediction mechanism on prediction accuracy. The results indicate that feedbacks from the neighboring nodes of a specific service requestor with high Hellinger similarity in our mechanism outperforms the best existing methods. We also show that utilizing the social trust model, which only considers a similarity measure, significantly improves the accuracy of the prediction mechanism. Furthermore, we evaluate the effectiveness of the proposed trust management system through a real-world SIoT use case. Our results demonstrate that the proposed mechanism is resilient to different types of network attacks, and it can accurately find the most proper and trustworthy service provider.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A new look at calendar anomalies: Multifractality and day of the week effect

Stock markets can become inefficient due to calendar anomalies known as day-of-the-week effect. Calendar anomalies are well-known in financial literature, but the phenomena remain to be explored in econophysics. In this paper we use multifractal analysis to evaluate if the temporal dynamics of market returns also exhibits calendar anomalies such as day-of-the-week effects. We apply the multifractal detrended fluctuation analysis (MF-DFA) to daily returns of market indices around the world for each day of the week. Our results indicate that individual days of the week are characterized by distinct multifractal properties. Monday returns tend to exhibit more persistent behavior and richer multifractal structures than other day-resolved returns. Shuffling the series reveals that multifractality arises both from a broad probability density function and from long-term correlations. From the time-dependent multifractal analysis we find that multifractal spectra for Monday returns are much wider than for other days during periods of financial crises. The presence of day-of-the-week effects in multifractal dynamics of market returns motivates further research on calendar anomalies from an econophysics perspective.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Impacts of export restrictions on the global personal protective equipment trade network during COVID-19

The COVID-19 pandemic has caused a dramatic surge in demand for personal protective equipment (PPE) worldwide. Many countries have imposed export restrictions on PPE to ensure the sufficient domestic supply. The surging demand and export restrictions cause shortage contagions on the global PPE trade network. Here, we develop an integrated network model, which integrates a metapopulation model and a threshold model, to investigate the shortage contagion patterns. The metapopulation model captures disease contagion across countries. The threshold model captures the shortage contagion on the global PPE trade network. Results show that, the shortage contagion patterns are mainly decided by top exporters. Export restrictions exacerbate the shortages of PPE and cause the shortage contagion to transmit even faster than the disease contagion. Besides, export restrictions lead to ineffective and inefficient allocation of PPE around the world, which has no benefits for the world to fight against the pandemic.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An Analytical Solution to the $k$-core Pruning Process

$k$-core decomposition is widely used to identify the center of a large network, it is a pruning process in which the nodes with degrees less than $k$ are recursively removed. Although the simplicity and effectiveness of this method facilitate its implementation on broad applications across many scientific fields, it produces few analytical results. We here simplify the existing theoretical framework to a simple iterative relationship and obtain the exact analytical solutions of the $k$-core pruning process on large uncorrelated networks. From these solutions we obtain such statistical properties as the degree distribution and the size of the remaining subgraph in each of the pruning steps. Our theoretical results resolve the long-lasting puzzle of the $k$-core pruning dynamics and provide an intuitive description of the dynamic process.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Assessing the attraction of cities on venture capital from a scaling law perspective

Cities are centers for the integration of capital and incubators of invention, and attracting venture capital (VC) is of great importance for cities to advance in innovative technology and business models towards a sustainable and prosperous future. Yet we still lack a quantitative understanding of the relationship between urban characteristics and VC activities. In this paper, we find a clear nonlinear scaling relationship between VC activities and the urban population of Chinese cities. In such nonlinear systems, the widely applied linear per capita indicators would be either biased to larger cities or smaller cities depends on whether it is superlinear or sublinear, while the residual of cities relative to the prediction of scaling law is a more objective and scale-invariant metric. %(i.e., independent of the city size). Such a metric can distinguish the effects of local dynamics and scaled growth induced by the change of population size. The spatiotemporal evolution of such metrics on VC activities reveals three distinct groups of cities, two of which stand out with increasing and decreasing trends, respectively. And the taxonomy results together with spatial analysis also signify different development modes between large urban agglomeration regions. Besides, we notice the evolution of scaling exponents on VC activities are of much larger fluctuations than on socioeconomic output of cities, and a conceptual model that focuses on the growth dynamics of different sized cities can well explain it, which we assume would be general to other scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

History-dependent percolation on multiplex networks

The structure of interconnected systems and its impact on the system dynamics is a much-studied cross-disciplinary topic. Although various critical phenomena have been found in different models, the study on the connections between different percolation transitions is still lacking. Here we propose a unified framework to study the origins of the discontinuous transitions of the percolation process on interacting networks. The model evolves in generations with the result of the present percolation depending on the previous state and thus is history-dependent. Both theoretical analysis and Monte Carlo simulations reveal that the nature of the transition remains the same at finite generations but exhibits an abrupt change for the infinite generation. We use brain functional correlation and morphological similarity data to show that our model also provides a general method to explore the network structure and can contribute to many practical applications, such as detecting the abnormal structures of human brain networks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Universal Urban Spreading Pattern of COVID-19 and Its Underlying Mechanism

Currently, the global situation of COVID-19 is aggravating, pressingly calling for efficient control and prevention measures. Understanding spreading pattern of COVID-19 has been widely recognized as a vital step for implementing non-pharmaceutical measures. Previous studies investigated such an issue in large-scale (e.g., inter-country or inter-state) scenarios while urban spreading pattern still remains an open issue. Here, we fill this gap by leveraging the trajectory data of 197,808 smartphone users (including 17,808 anonymous confirmed cases) in 9 cities in China. We find a universal spreading pattern in all cities: the spatial distribution of confirmed cases follows a power-law-like model and the spreading centroid is time-invariant. Moreover, we reveal that human mobility in a city drives the spatialtemporal spreading process: long average travelling distance results in a high growth rate of spreading radius and wide spatial diffusion of confirmed cases. With such insight, we adopt Kendall model to simulate urban spreading of COVID-19 that can well fit the real spreading process. Our results unveil the underlying mechanism behind the spatial-temporal urban evolution of COVID-19, and can be used to evaluate the performance of mobility restriction policies implemented by many governments and to estimate the evolving spreading situation of COVID-19.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unveiling the Physics of the Mutual Interactions in Paramagnets

In real paramagnets, there is always a subtle many-body contribution to the system's energy, which can be regarded as a small effective local magnetic field $B_{loc}$. Usually, it is neglected, since it is very small when compared with thermal fluctuations and/or external magnetic fields $B$. Nevertheless, as both the temperature $T \rightarrow$ 0K and $B \rightarrow$ 0T, such many-body contributions become ubiquitous. Here, employing the magnetic Grüneisen parameter $Γ_{mag}$ and entropy arguments, we report on the pivotal role played by the mutual interactions in the regime of ultra-low-$T$ and vanishing $B$. Our key results are: $i$) absence of a genuine zero-field quantum phase transition due to the presence of $B_{loc}$; $ii$) connection between the canonical definition of temperature and $Γ_{mag}$; and $iii$) possibility of performing adiabatic magnetization by only manipulating the mutual interactions. Our findings unveil unprecedented aspects emerging from the mutual interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Network analysis of the worldwide footballer transfer market

The transfer of football players is an important part in football games. Most studies on the transfer of football players focus on the transfer system and transfer fees but not on the transfer behavior itself. Based on the 470,792 transfer records from 1990 to 2016 among 23,605 football clubs in 206 countries and regions, we construct a directed footballer transfer network (FTN), where the nodes are the football clubs and the links correspond to the footballer transfers. A systemic analysis is conduced on the topological properties of the FTN. We find that the in-degrees, out-degrees, in-strengths and out-strengths of nodes follow bimodal distributions (a power law with exponential decay), while the distribution of link weights has a power-law tail. We further figure out the correlations between node degrees, node strengths and link weights. We also investigate the general characteristics of different measures of network centrality. Our network analysis of the global footballer transfer market sheds new lights into the investigation of the characteristics of transfer activities.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Model of Brain Activation Predicts the Neural Collective Influence Map of the Brain

Efficient complex systems have a modular structure, but modularity does not guarantee robustness, because efficiency also requires an ingenious interplay of the interacting modular components. The human brain is the elemental paradigm of an efficient robust modular system interconnected as a network of networks (NoN). Understanding the emergence of robustness in such modular architectures from the interconnections of its parts is a long-standing challenge that has concerned many scientists. Current models of dependencies in NoN inspired by the power grid express interactions among modules with fragile couplings that amplify even small shocks, thus preventing functionality. Therefore, we introduce a model of NoN to shape the pattern of brain activations to form a modular environment that is robust. The model predicts the map of neural collective influencers (NCIs) in the brain, through the optimization of the influence of the minimal set of essential nodes responsible for broadcasting information to the whole-brain NoN. Our results suggest new intervention protocols to control brain activity by targeting influential neural nodes predicted by network theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

$k$-core percolation on complex networks: Comparing random, localized and targeted attacks

The type of malicious attack inflicting on networks greatly influences their stability under ordinary percolation in which a node fails when it becomes disconnected from the giant component. Here we study its generalization, $k$-core percolation, in which a node fails when it loses connection to a threshold $k$ number of neighbors. We study and compare analytically and by numerical simulations of $k$-core percolation the stability of networks under random attacks (RA), localized attacks (LA) and targeted attacks (TA), respectively. By mapping a network under LA or TA into an equivalent network under RA, we find that in both single and interdependent networks, TA exerts the greatest damage to the core structure of a network. We also find that for Erdős-Rényi (ER) networks, LA and RA exert equal damage to the core structure whereas for scale-free (SF) networks, LA exerts much more damage than RA does to the core structure.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How fear of future outcomes affects social dynamics

Mutualistic relationships among the different species are ubiquitous in nature. To prevent mutualism from slipping into antagonism, a host often invokes a "carrot and stick" approach towards symbionts with a stabilizing effect on their symbiosis. In open human societies, a mutualistic relationship arises when a native insider population attracts outsiders with benevolent incentives in hope that the additional labor will improve the standard of all. A lingering question, however, is the extent to which insiders are willing to tolerate outsiders before mutualism slips into antagonism. To test the assertion by Karl Popper that unlimited tolerance leads to the demise of tolerance, we model a society under a growing incursion from the outside. Guided by their traditions of maintaining the social fabric and prizing tolerance, the insiders reduce their benevolence toward the growing subpopulation of outsiders but do not invoke punishment. This reduction of benevolence intensifies as less tolerant insiders (e.g., "radicals") openly renounce benevolence. Although more tolerant insiders maintain some level of benevolence, they may also tacitly support radicals out of fear for the future. If radicals and their tacit supporters achieve a critical majority, herd behavior ensues and the relation between the insider and outsider subpopulations turns antagonistic. To control the risk of unwanted social dynamics, we map the parameter space within which the tolerance of insiders is in balance with the assimilation of outsiders, the tolerant insiders maintain a sustainable majority, and any reduction in benevolence occurs smoothly. We also identify the circumstances that cause the relations between insiders and outsiders to collapse or that lead to the dominance of the outsiders.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Large Slowdown of Water Dynamics at Stacked Phospholipid Membranes for Increasing Hydration Level: All-Atoms Molecular Dynamics

Water hydrating phospholipid membranes determine their stability and function, as well as their interaction with other molecules. In this article we study, using all-atom molecular dynamics simulations, the rotational and translational dynamical properties of water molecules confined in stacked phospholipid membranes at different levels of hydration, from poorly hydrated to a completely hydrated membrane. We find that both the translational and the reorientation dynamics of water are dramatically slowed down as the hydration is reduced. Analyzing in details the structure and dynamics of the hydrogen bond at the interface, we show that both those among water molecules and those between water and lipids slow down by decreasing the hydration, however the latter are always slower than the former. By increasing hydration, water saturates all the possible hydrogen bonds with the lipids and, by further increase of hydration, the hydrogen bonds among waters becomes the majority. However, the dynamics of the water-lipids hydrogen bonds becomes up to one order of magnitude slower than that of the water-water hydrogen bonds, inducing a large slowing down of the dynamics of the entire system even at large hydration level.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Locating the source of diffusion in complex networks by time-reversal backward spreading

Locating the source that triggers a dynamical process is a fundamental but challenging problem in complex networks, ranging from epidemic spreading in society and on the Internet to cancer metastasis in the human body. An accurate localization of the source is inherently limited by our ability to simultaneously access the information of all nodes in a large-scale complex network. This thus raises two critical questions: how do we locate the source from incomplete information and can we achieve full localization of sources at any possible location from a given set of observable nodes. Here we develop a time-reversal backward spreading algorithm to locate the source of a diffusion-like process efficiently and propose a general locatability condition. We test the algorithm by employing epidemic spreading and consensus dynamics as typical dynamical processes and apply it to the H1N1 pandemic in China. We find that the sources can be precisely located in arbitrary networks insofar as the locatability condition is assured. Our tools greatly improve our ability to locate the source of diffusion in complex networks based on limited accessibility of nodal information. Moreover, they have implications for controlling a variety of dynamical processes taking place on complex networks, such as inhibiting epidemics, slowing the spread of rumors, pollution control and environmental protection.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Percolation of networks with directed dependency links

The self-consistent probabilistic approach has proven itself powerful in studying the percolation behavior of interdependent or multiplex networks without tracking the percolation process through each cascading step. In order to understand how directed dependency links impact criticality, we employ this approach to study the percolation properties of networks with both undirected connectivity links and directed dependency links. We find that when a random network with a given degree distribution undergoes a second-order phase transition, the critical point and the unstable regime surrounding the second-order phase transition regime are determined by the proportion of nodes that do not depend on any other nodes. Moreover, we also find that the triple point and the boundary between first- and second-order transitions are determined by the proportion of nodes that depend on no more than one node. This implies that it is maybe general for multiplex network systems, some important properties of phase transitions can be determined only by a few parameters. We illustrate our findings using Erdos-Renyi (ER) networks.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Predicting the rise of right-wing populism in response to unbalanced immigration

Among the central tenets of globalization is free migration of labor. Although much has been written about its benefits, little is known about the limitations of globalization, including how immigration affects the anti-globalist sentiment. Analyzing polls data, we find that over the last three years in a group of EU countries affected by the recent migrant crisis, the percentage of right-wing (RW) populist voters in a given country depends on the prevalence of immigrants in this country's population and the total immigration inflow into the entire EU. The latter is likely due to the EU resembling a supranational state, where the lack of inner borders causes that "somebody else's problem" easily turns into "my problem". We further find that the increase in the percentage of RW voters substantially surpasses the immigration inflow, implying that if this process continues, RW populism may democratically prevail and eventually lead to a demise of globalization. We present evidence for tipping points in relation to the rise of RW populism. Finally, we model these empirical findings using a complex network framework wherein the success of globalization largely rests on the balance between immigration and immigrant integration.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stochastic model of financial markets reproducing scaling and memory in volatility return intervals

We investigate the volatility return intervals in the NYSE and FOREX markets. We explain previous empirical findings using a model based on the interacting agent hypothesis instead of the widely-used efficient market hypothesis. We derive macroscopic equations based on the microscopic herding interactions of agents and find that they are able to reproduce various stylized facts of different markets and different assets with the same set of model parameters. We show that the power-law properties and the scaling of return intervals and other financial variables have a similar origin and could be a result of a general class of non-linear stochastic differential equations derived from a master equation of an agent system that is coupled by herding interactions. Specifically, we find that this approach enables us to recover the volatility return interval statistics as well as volatility probability and spectral densities for the NYSE and FOREX markets, for different assets, and for different time-scales. We find also that the historical S\&P500 monthly series exhibits the same volatility return interval properties recovered by our proposed model. Our statistical results suggest that human herding is so strong that it persists even when other evolving fluctuations perturbate the financial system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Suppressing disease spreading by using information diffusion on multiplex networks

Although there is always an interplay between the dynamics of information diffusion and disease spreading, the empirical research on the systemic coevolution mechanisms connecting these two spreading dynamics is still lacking. Here we investigate the coevolution mechanisms and dynamics between information and disease spreading by utilizing real data and a proposed spreading model on multiplex network. Our empirical analysis finds asymmetrical interactions between the information and disease spreading dynamics. Our results obtained from both the theoretical framework and extensive stochastic numerical simulations suggest that an information outbreak can be triggered in a communication network by its own spreading dynamics or by a disease outbreak on a contact network, but that the disease threshold is not affected by information spreading. Our key finding is that there is an optimal information transmission rate that markedly suppresses the disease spreading. We find that the time evolution of the dynamics in the proposed model qualitatively agrees with the real-world spreading processes at the optimal information transmission rate.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Detrended partial cross-correlation analysis of two nonstationary time series influenced by common external forces

When common factors strongly influence two power-law cross-correlated time series recorded in complex natural or social systems, using classic detrended cross-correlation analysis (DCCA) without considering these common factors will bias the results. We use detrended partial cross-correlation analysis (DPXA) to uncover the intrinsic power-law cross-correlations between two simultaneously recorded time series in the presence of nonstationarity after removing the effects of other time series acting as common forces. The DPXA method is a generalization of the detrended cross-correlation analysis that takes into account partial correlation analysis. We demonstrate the method by using bivariate fractional Brownian motions contaminated with a fractional Brownian motion. We find that the DPXA is able to recover the analytical cross Hurst indices, and thus the multi-scale DPXA coefficients are a viable alternative to the conventional cross-correlation coefficient. We demonstrate the advantage of the DPXA coefficients over the DCCA coefficients by analyzing contaminated bivariate fractional Brownian motions. We calculate the DPXA coefficients and use them to extract the intrinsic cross-correlation between crude oil and gold futures by taking into consideration the impact of the US dollar index. We develop the multifractal DPXA (MF-DPXA) method in order to generalize the DPXA method and investigate multifractal time series. We analyze multifractal binomial measures masked with strong white noises and find that the MF-DPXA method quantifies the hidden multifractal nature while the MF-DCCA method fails.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Division of labor, skill complementarity, and heterophily in socioeconomic networks

Constituents of complex systems interact with each other and self-organize to form complex networks. Empirical results show that the link formation process of many real networks follows either the global principle of popularity or the local principle of similarity or a tradeoff between the two. In particular, it has been shown that in social networks individuals exhibit significant homophily when choosing their collaborators. We demonstrate, however, that in populations in which there is a division of labor, skill complementarity is an important factor in the formation of socioeconomic networks and an individual's choice of collaborators is strongly affected by heterophily. We analyze 124 evolving virtual worlds of a popular "massively multiplayer online role-playing game" (MMORPG) in which people belong to three different professions and are allowed to work and interact with each other in a somewhat realistic manner. We find evidence of heterophily in the formation of collaboration networks, where people prefer to forge social ties with people who have professions different from their own. We then construct an economic model to quantify the heterophily by assuming that individuals in socioeconomic systems choose collaborators that are of maximum utility. The results of model calibration confirm the presence of heterophily. Both empirical analysis and model calibration show that the heterophilous feature is persistent along the evolution of virtual worlds. We also find that the degree of complementarity in virtual societies is positively correlated with their economic output. Our work sheds new light on the scientific research utility of virtual worlds for studying human behaviors in complex socioeconomic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Early warning of large volatilities based on recurrence interval analysis in Chinese stock markets

Being able to forcast extreme volatility is a central issue in financial risk management. We present a large volatility predicting method based on the distribution of recurrence intervals between volatilities exceeding a certain threshold $Q$ for a fixed expected recurrence time $τ_Q$. We find that the recurrence intervals are well approximated by the $q$-exponential distribution for all stocks and all $τ_Q$ values. Thus a analytical formula for determining the hazard probability $W(Δt |t)$ that a volatility above $Q$ will occur within a short interval $Δt$ if the last volatility exceeding $Q$ happened $t$ periods ago can be directly derived from the $q$-exponential distribution, which is found to be in good agreement with the empirical hazard probability from real stock data. Using these results, we adopt a decision-making algorithm for triggering the alarm of the occurrence of the next volatility above $Q$ based on the hazard probability. Using a "receiver operator characteristic" (ROC) analysis, we find that this predicting method efficiently forecasts the occurrance of large volatility events in real stock data. Our analysis may help us better understand reoccurring large volatilities and more accurately quantify financial risks in stock markets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Echo chambers in the age of misinformation

The wide availability of user-provided content in online social media facilitates the aggregation of people around common interests, worldviews, and narratives. Despite the enthusiastic rhetoric on the part of some that this process generates "collective intelligence", the WWW also allows the rapid dissemination of unsubstantiated conspiracy theories that often elicite rapid, large, but naive social responses such as the recent case of Jade Helm 15 -- where a simple military exercise turned out to be perceived as the beginning of the civil war in the US. We study how Facebook users consume information related to two different kinds of narrative: scientific and conspiracy news. We find that although consumers of scientific and conspiracy stories present similar consumption patterns with respect to content, the sizes of the spreading cascades differ. Homogeneity appears to be the primary driver for the diffusion of contents, but each echo chamber has its own cascade dynamics. To mimic these dynamics, we introduce a data-driven percolation model on signed networks.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiple Tipping Points and Optimal Repairing in Interacting Networks

Systems that comprise many interacting dynamical networks, such as the human body with its biological networks or the global economic network consisting of regional clusters, often exhibit complicated collective dynamics. To understand the collective behavior of such systems, we investigate a model of interacting networks exhibiting the fundamental processes of failure, damage spread, and recovery. We find a very rich phase diagram that becomes exponentially more complex as the number of networks is increased. In the simplest example of $n=2$ interacting networks we find two critical points, 4 triple points, 10 allowed transitions, and two "forbidden" transitions, as well as complex hysteresis loops. Remarkably, we find that triple points play the dominant role in constructing the optimal repairing strategy in damaged interacting systems. To support our model, we analyze an example of real interacting financial networks and find evidence of rapid dynamical transitions between well-defined states, in agreement with the predictions of our model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Predicting the epidemic threshold of the susceptible-infected-recovered model

Researchers have developed several theoretical methods for predicting epidemic thresholds, including the mean-field like (MFL) method, the quenched mean-field (QMF) method, and the dynamical message passing (DMP) method. When these methods are applied to predict epidemic threshold they often produce differing results and their relative levels of accuracy are still unknown. We systematically analyze these two issues---relationships among differing results and levels of accuracy---by studying the susceptible-infected-recovered (SIR) model on uncorrelated configuration networks and a group of 56 real-world networks. In uncorrelated configuration networks the MFL and DMP methods yield identical predictions that are larger and more accurate than the prediction generated by the QMF method. When compared to the 56 real-world networks, the epidemic threshold obtained by the DMP method is closer to the actual epidemic threshold because it incorporates full network topology information and some dynamical correlations. We find that in some scenarios---such as networks with positive degree-degree correlations, with an eigenvector localized on the high $k$-core nodes, or with a high level of clustering---the epidemic threshold predicted by the MFL method, which uses the degree distribution as the only input parameter, performs better than the other two methods. We also find that the performances of the three predictions are irregular versus modularity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The cost of attack in competing networks

Real-world attacks can be interpreted as the result of competitive interactions between networks, ranging from predator-prey networks to networks of countries under economic sanctions. Although the purpose of an attack is to damage a target network, it also curtails the ability of the attacker, which must choose the duration and magnitude of an attack to avoid negative impacts on its own functioning. Nevertheless, despite the large number of studies on interconnected networks, the consequences of initiating an attack have never been studied. Here, we address this issue by introducing a model of network competition where a resilient network is willing to partially weaken its own resilience in order to more severely damage a less resilient competitor. The attacking network can take over the competitor nodes after their long inactivity. However, due to a feedback mechanism the takeovers weaken the resilience of the attacking network. We define a conservation law that relates the feedback mechanism to the resilience dynamics for two competing networks. Within this formalism, we determine the cost and optimal duration of an attack, allowing a network to evaluate the risk of initiating hostilities.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The influence of the broadness of the degree distribution on network's robustness: comparing localized attack and random attack

The stability of networks is greatly influenced by their degree distributions and in particular by their broadness. Networks with broader degree distributions are usually more robust to random failures but less robust to localized attacks. To better understand the effect of the broadness of the degree distribution we study here two models where the broadness is controlled and compare their robustness against localized attacks (LA) and random attacks (RA). We study analytically and by numerical simulations the cases where the degrees in the networks follow a Bi-Poisson distribution $P(k)=αe^{-λ_1}\frac{λ_1^k}{k!}+(1-α) e^{-λ_2}\frac{λ_2^k}{k!},α\in[0,1]$, and a Gaussian distribution $P(k)=A \cdot exp{(-\frac{(k-μ)^2}{2σ^2})}$ with a normalization constant $A$ where $k\geq 0$. In the Bi-Poisson distribution the broadness is controlled by the values of $α$, $λ_1$ and $λ_2$, while in the Gaussian distribution it is controlled by the standard deviation, $σ$. We find that only for $α=0$ or $α=1$, namely degrees obeying a pure Poisson distribution, LA and RA are the same but for all other cases networks are more vulnerable under LA compared to RA. For Gaussian distribution, with an average degree $μ$ fixed, we find that when $σ^2$ is smaller than $μ$ the network is more vulnerable against random attack. However, when $σ^2$ is larger than $μ$ the network becomes more vulnerable against localized attack. Similar qualitative results are also shown for interdependent networks.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Universality of market superstatistics

We use a continuous-time random walk (CTRW) to model market fluctuation data from times when traders experience excessive losses or excessive profits. We analytically derive "superstatistics" that accurately model empirical market activity data (supplied by Bogachev, Ludescher, Tsallis, and Bunde)that exhibit transition thresholds. We measure the interevent times between excessive losses and excessive profits, and use the mean interevent time as a control variable to derive a universal description of empirical data collapse. Our superstatistic value is a weighted sum of two components, (i) a powerlaw corrected by the lower incomplete gamma function, which asymptotically tends toward robustness but initially gives an exponential, and (ii) a powerlaw damped by the upper incomplete gamma function, which tends toward the power-law only during short interevent times. We find that the scaling shape exponents that drive both components subordinate themselves and a "superscaling" configuration emerges. We use superstatistics to describe the hierarchical activity when component (i) reproduces the negative feedback and component (ii) reproduces the stylized fact of volatility clustering. Our results indicate that there is a functional (but not literal) balance between excessive profits and excessive losses that can be described using the same body of superstatistics, but different calibration values and driving parameters.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Competing for Attention in Social Media under Information Overload Conditions

Although the many forms of modern social media have become major channels for the dissemination of information, they are becoming overloaded because of the rapidly-expanding number of information feeds. We analyze the expanding user-generated content in Sina Weibo, the largest micro-blog site in China, and find evidence that popular messages often follow a mechanism that differs from that found in the spread of disease, in contrast to common believe. In this mechanism, an individual with more friends needs more repeated exposures to spread further the information. Moreover, our data suggest that in contrast to epidemics, for certain messages the chance of an individual to share the message is proportional to the fraction of its neighbours who shared it with him/her. Thus the greater the number of friends an individual has the greater the number of repeated contacts needed to spread the message, which is a result of competition for attention. We model this process using a fractional susceptible infected recovered (FSIR) model, where the infection probability of a node is proportional to its fraction of infected neighbors. Our findings have dramatic implications for information contagion. For example, using the FSIR model we find that real-world social networks have a finite epidemic threshold. This is in contrast to the zero threshold that conventional wisdom derives from disease epidemic models. This means that when individuals are overloaded with excess information feeds, the information either reaches out the population if it is above the critical epidemic threshold, or it would never be well received, leading to only a handful of information contents that can be widely spread throughout the population.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Correlation between centrality metrics and their application to the opinion model

In recent decades, a number of centrality metrics describing network properties of nodes have been proposed to rank the importance of nodes. In order to understand the correlations between centrality metrics and to approximate a high-complexity centrality metric by a strongly correlated low-complexity metric, we first study the correlation between centrality metrics in terms of their Pearson correlation coefficient and their similarity in ranking of nodes. In addition to considering the widely used centrality metrics, we introduce a new centrality measure, the degree mass. The m order degree mass of a node is the sum of the weighted degree of the node and its neighbors no further than m hops away. We find that the B_{n}, the closeness, and the components of x_{1} are strongly correlated with the degree, the 1st-order degree mass and the 2nd-order degree mass, respectively, in both network models and real-world networks. We then theoretically prove that the Pearson correlation coefficient between x_{1} and the 2nd-order degree mass is larger than that between x_{1} and a lower order degree mass. Finally, we investigate the effect of the inflexible antagonists selected based on different centrality metrics in helping one opinion to compete with another in the inflexible antagonists opinion model. Interestingly, we find that selecting the inflexible antagonists based on the leverage, the B_{n}, or the degree is more effective in opinion-competition than using other centrality metrics in all types of networks. This observation is supported by our previous observations, i.e., that there is a strong linear correlation between the degree and the B_{n}, as well as a high centrality similarity between the leverage and the degree.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Emergence of statistically validated financial intraday lead-lag relationships

According to the leading models in modern finance, the presence of intraday lead-lag relationships between financial assets is negligible in efficient markets. With the advance of technology, however, markets have become more sophisticated. To determine whether this has resulted in an improved market efficiency, we investigate whether statistically significant lagged correlation relationships exist in financial markets. We introduce a numerical method to statistically validate links in correlation-based networks, and employ our method to study lagged correlation networks of equity returns in financial markets. Crucially, our statistical validation of lead-lag relationships accounts for multiple hypothesis testing over all stock pairs. In an analysis of intraday transaction data from the periods 2002--2003 and 2011--2012, we find a striking growth in the networks as we increase the frequency with which we sample returns. We compute how the number of validated links and the magnitude of correlations change with increasing sampling frequency, and compare the results between the two data sets. Finally, we compare topological properties of the directed correlation-based networks from the two periods using the in-degree and out-degree distributions and an analysis of three-node motifs. Our analysis suggests a growth in both the efficiency and instability of financial markets over the past decade.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

NetSci High: Bringing Network Science Research to High Schools

We present NetSci High, our NSF-funded educational outreach program that connects high school students who are underrepresented in STEM (Science Technology Engineering and Mathematics), and their teachers, with regional university research labs and provides them with the opportunity to work with researchers and graduate students on team-based, year-long network science research projects, culminating in a formal presentation at a network science conference. This short paper reports the content and materials that we have developed to date, including lesson plans and tools for introducing high school students and teachers to network science; empirical evaluation data on the effect of participation on students' motivation and interest in pursuing STEM careers; the application of professional development materials for teachers that are intended to encourage them to use network science concepts in their lesson plans and curriculum; promoting district-level interest and engagement; best practices gained from our experiences; and the future goals for this project and its subsequent outgrowth.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-consensus opinion model on directed networks

Dynamic social opinion models have been widely studied on undirected networks, and most of them are based on spin interaction models that produce a consensus. In reality, however, many networks such as Twitter and the World Wide Web are directed and are composed of both unidirectional and bidirectional links. Moreover, from choosing a coffee brand to deciding who to vote for in an election, two or more competing opinions often coexist. In response to this ubiquity of directed networks and the coexistence of two or more opinions in decision-making situations, we study a non-consensus opinion model introduced by Shao et al. \cite{shao2009dynamic} on directed networks. We define directionality $ξ$ as the percentage of unidirectional links in a network, and we use the linear correlation coefficient $ρ$ between the indegree and outdegree of a node to quantify the relation between the indegree and outdegree. We introduce two degree-preserving rewiring approaches which allow us to construct directed networks that can have a broad range of possible combinations of directionality $ξ$ and linear correlation coefficient $ρ$ and to study how $ξ$ and $ρ$ impact opinion competitions. We find that, as the directionality $ξ$ or the indegree and outdegree correlation $ρ$ increases, the majority opinion becomes more dominant and the minority opinion's ability to survive is lowered.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Partial correlation analysis: Applications for financial markets

The presence of significant cross-correlations between the synchronous time evolution of a pair of equity returns is a well-known empirical fact. The Pearson correlation is commonly used to indicate the level of similarity in the price changes for a given pair of stocks, but it does not measure whether other stocks influence the relationship between them. To explore the influence of a third stock on the relationship between two stocks, we use a partial correlation measurement to determine the underlying relationships between financial assets. Building on previous work, we present a statistically robust approach to extract the underlying relationships between stocks from four different financial markets: the United States, the United Kingdom, Japan, and India. This methodology provides new insights into financial market dynamics and uncovers implicit influences in play between stocks. To demonstrate the capabilities of this methodology, we (i) quantify the influence of different companies and, by studying market similarity across time, present new insights into market structure and market stability, and (ii) we present a practical application, which provides information on the how a company is influenced by different economic sectors, and how the sectors interact with each other. These examples demonstrate the effectiveness of this methodology in uncovering information valuable for a range of individuals, including not only investors and traders but also regulators and policy makers.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Predicting market instability: New dynamics between volume and volatility

Econophysics and econometrics agree that there is a correlation between volume and volatility in a time series. Using empirical data and their distributions, we further investigate this correlation and discover new ways that volatility and volume interact, particularly when the levels of both are high. We find that the distribution of the volume-conditional volatility is well fit by a power-law function with an exponential cutoff. We find that the volume-conditional volatility distribution scales with volume, and collapses these distributions to a single curve. We exploit the characteristics of the volume-volatility scatter plot to find a strong correlation between logarithmic volume and a quantity we define as local maximum volatility (LMV), which indicates the largest volatility observed in a given range of trading volumes. This finding supports our empirical analysis showing that volume is an excellent predictor of the maximum value of volatility for both same-day and near-future time periods. We also use a joint conditional probability that includes both volatility and volume to demonstrate that invoking both allows us to better predict the largest next-day volatility than invoking either one alone.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ranking the Economic Importance of Countries and Industries

In the current era of worldwide stock market interdependencies, the global financial village has become increasingly vulnerable to systemic collapse. The recent global financial crisis has highlighted the necessity of understanding and quantifying interdependencies among the world's economies, developing new effective approaches to risk evaluation, and providing mitigating solutions. We present a methodological framework for quantifying interdependencies in the global market and for evaluating risk levels in the world-wide financial network. The resulting information will enable policy and decision makers to better measure, understand, and maintain financial stability. We use the methodology to rank the economic importance of each industry and country according to the global damage that would result from their failure. Our quantitative results shed new light on China's increasing economic dominance over other economies, including that of the USA, to the global economy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Reputation and Impact in Academic Careers

Reputation is an important social construct in science, which enables informed quality assessments of both publications and careers of scientists in the absence of complete systemic information. However, the relation between reputation and career growth of an individual remains poorly understood, despite recent proliferation of quantitative research evaluation methods. Here we develop an original framework for measuring how a publication's citation rate $Δc$ depends on the reputation of its central author $i$, in addition to its net citation count $c$. To estimate the strength of the reputation effect, we perform a longitudinal analysis on the careers of 450 highly-cited scientists, using the total citations $C_{i}$ of each scientist as his/her reputation measure. We find a citation crossover $c_{\times}$ which distinguishes the strength of the reputation effect. For publications with $c < c_{\times}$, the author's reputation is found to dominate the annual citation rate. Hence, a new publication may gain a significant early advantage corresponding to roughly a 66% increase in the citation rate for each tenfold increase in $C_{i}$. However, the reputation effect becomes negligible for highly cited publications meaning that for $c\geq c_{\times}$ the citation rate measures scientific impact more transparently. In addition we have developed a stochastic reputation model, which is found to reproduce numerous statistical observations for real careers, thus providing insight into the microscopic mechanisms underlying cumulative advantage in science.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Calling patterns in human communication dynamics

Modern technologies not only provide a variety of communication modes, e.g., texting, cellphone conversation, and online instant messaging, but they also provide detailed electronic traces of these communications between individuals. These electronic traces indicate that the interactions occur in temporal bursts. Here, we study the inter-call durations of the 100,000 most-active cellphone users of a Chinese mobile phone operator. We confirm that the inter-call durations follow a power-law distribution with an exponential cutoff at the population level but find differences when focusing on individual users. We apply statistical tests at the individual level and find that the inter-call durations follow a power-law distribution for only 3460 individuals (3.46%). The inter-call durations for the majority (73.34%) follow a Weibull distribution. We quantify individual users using three measures: out-degree, percentage of outgoing calls, and communication diversity. We find that the cellphone users with a power-law duration distribution fall into three anomalous clusters: robot-based callers, telecom frauds, and telephone sales. This information is of interest to both academics and practitioners, mobile telecom operator in particular. In contrast, the individual users with a Weibull duration distribution form the fourth cluster of ordinary cellphone users. We also discover more information about the calling patterns of these four clusters, e.g., the probability that a user will call the $c_r$-th most contact and the probability distribution of burst sizes. Our findings may enable a more detailed analysis of the huge body of data contained in the logs of massive users.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Carbon-dioxide emissions trading and hierarchical structure in worldwide finance and commodities markets

In a highly interdependent economic world, the nature of relationships between financial entities is becoming an increasingly important area of study. Recently, many studies have shown the usefulness of minimal spanning trees (MST) in extracting interactions between financial entities. Here, we propose a modified MST network whose metric distance is defined in terms of cross-correlation coefficient absolute values, enabling the connections between anticorrelated entities to manifest properly. We investigate 69 daily time series, comprising three types of financial assets: 28 stock market indicators, 21 currency futures, and 20 commodity futures. We show that though the resulting MST network evolves over time, the financial assets of similar type tend to have connections which are stable over time. In addition, we find a characteristic time lag between the volatility time series of the stock market indicators and those of the EU CO2 emission allowance (EUA) and crude oil futures (WTI). This time lag is given by the peak of the cross-correlation function of the volatility time series EUA (or WTI) with that of the stock market indicators, and is markedly different (>20 days) from 0, showing that the volatility of stock market indicators today can predict the volatility of EU emissions allowances and of crude oil in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cascading Failures in Bi-partite Graphs: Model for Systemic Risk Propagation

As economic entities become increasingly interconnected, a shock in a financial network can provoke significant cascading failures throughout the system. To study the systemic risk of financial systems, we create a bi-partite banking network model composed of banks and bank assets and propose a cascading failure model to describe the risk propagation process during crises. We empirically test the model with 2007 US commercial banks balance sheet data and compare the model prediction of the failed banks with the real failed banks after 2007. We find that our model efficiently identifies a significant portion of the actual failed banks reported by Federal Deposit Insurance Corporation. The results suggest that this model could be useful for systemic risk stress testing for financial systems. The model also identifies that commercial rather than residential real estate assets are major culprits for the failure of over 350 US commercial banks during 2008-2011.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effect of the Interconnected Network Structure on the Epidemic Threshold

Most real-world networks are not isolated. In order to function fully, they are interconnected with other networks, and this interconnection influences their dynamic processes. For example, when the spread of a disease involves two species, the dynamics of the spread within each species (the contact network) differs from that of the spread between the two species (the interconnected network). We model two generic interconnected networks using two adjacency matrices, A and B, in which A is a 2N*2N matrix that depicts the connectivity within each of two networks of size N, and B a 2N*2N matrix that depicts the interconnections between the two. Using an N-intertwined mean-field approximation, we determine that a critical susceptable-infected-susceptable (SIS) epidemic threshold in two interconnected networks is 1/λ1(A+αB), where the infection rate is βwithin each of the two individual networks and αβin the interconnected links between the two networks and λ1(A+αB) is the largest eigenvalue of the matrix A+αB. In order to determine how the epidemic threshold is dependent upon the structure of interconnected networks, we analytically derive λ1(A+αB) using perturbation approximation for small and large α, the lower and upper bound for any αas a function of the adjacency matrix of the two individual networks, and the interconnections between the two and their largest eigenvalues/eigenvectors. We verify these approximation and boundary values for λ1(A+αB) using numerical simulations, and determine how component network features affect λ1(A+αB).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Percolation of Partially Interdependent Scale-free Networks

We study the percolation behavior of two interdependent scale-free (SF) networks under random failure of 1-$p$ fraction of nodes. Our results are based on numerical solutions of analytical expressions and simulations. We find that as the coupling strength between the two networks $q$ reduces from 1 (fully coupled) to 0 (no coupling), there exist two critical coupling strengths $q_1$ and $q_2$, which separate three different regions with different behavior of the giant component as a function of $p$. (i) For $q \geq q_1$, an abrupt collapse transition occurs at $p=p_c$. (ii) For $q_2<q<q_1$, the giant component has a hybrid transition combined of both, abrupt decrease at a certain $p=p^{jump}_c$ followed by a smooth decrease to zero for $p < p^{jump}_c$ as $p$ decreases to zero. (iii) For $q \leq q_2$, the giant component has a continuous second-order transition (at $p=p_c$). We find that $(a)$ for $λ\leq 3$, $q_1 \equiv 1$; and for $λ> 3$, $q_1$ decreases with increasing $λ$. $(b)$ In the hybrid transition, at the $q_2 < q < q_1$ region, the mutual giant component $P_{\infty}$ jumps discontinuously at $p=p^{jump}_c$ to a very small but non-zero value, and when reducing $p$, $P_{\infty}$ continuously approaches to 0 at $p_c = 0$ for $λ< 3$ and at $p_c > 0$ for $λ> 3$. Thus, the known theoretical $p_c=0$ for a single network with $λ\leqslant 3$ is expected to be valid also for strictly partial interdependent networks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Percolation on interacting networks with feedback-dependency links

When real networks are considered, coupled networks with connectivity and feedback-dependency links are not rare but more general. Here we develop a mathematical framework and study numerically and analytically percolation of interacting networks with feedback-dependency links. We find that when nodes of between networks are lowly connected, the system undergoes from second order transition through hybrid order transition to first order transition as coupling strength increases. And, as average degree of each inter-network increases, first order region becomes smaller and second-order region becomes larger but hybrid order region almost keep constant. Especially, the results implies that average degree \bar{k} between intra-networks has a little influence on robustness of system for weak coupling strength, but for strong coupling strength corresponding to first order transition system become robust as \bar{k} increases. However, when average degree k of inter-network is increased, the system become robust for all coupling strength. Additionally, when nodes of between networks are highly connected, the hybrid order region disappears and the system first order region becomes larger and secondorder region becomes smaller. Moreover, we find that the existence of feedback dependency links between interconnecting networks makes the system extremely vulnerable by comparing non-feedback condition for the same parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Robustness of Network of Networks with Interdependent and Interconnected links

Robustness of network of networks (NON) has been studied only for dependency coupling (J.X. Gao et. al., Nature Physics, 2012) and only for connectivity coupling (E.A. Leicht and R.M. D Souza, arxiv:0907.0894). The case of network of n networks with both interdependent and interconnected links is more complicated, and also more closely to real-life coupled network systems. Here we develop a framework to study analytically and numerically the robustness of this system. For the case of starlike network of n ER networks, we find that the system undergoes from second order to first order phase transition as coupling strength q increases. We find that increasing intra-connectivity links or inter-connectivity links can increase the robustness of the system, while the interdependency links decrease its robustness. Especially when q=1, we find exact analytical solutions of the giant component and the first order transition point. Understanding the robustness of network of networks with interdependent and interconnected links is helpful to design resilient infrastructures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Robustness of partially interdependent network formed of clustered networks

Clustering, or transitivity has been observed in real networks and its effects on their structure and function has been discussed extensively. The focus of these studies has been on clustering of single networks while the effect of clustering on the robustness of coupled networks received very little attention. Only the case of a pair of fully coupled networks with clustering has been studied recently. Here we generalize the study of clustering of a fully coupled pair of networks to the study of partially interdependent network of networks with clustering within the network components. We show both analytically and numerically, how clustering within the networks, affects the percolation properties of interdependent networks, including percolation threshold, size of giant component and critical coupling point where first order phase transition changes to second order phase transition as the coupling between the networks reduces. We study two types of clustering: one type proposed by Newman where the average degree is kept constant while changing the clustering and the other proposed by Hackett $et$ $al.$ where the degree distribution is kept constant. The first type of clustering is treated both analytically and numerically while the second one is treated only numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Systemic risk and spatiotemporal dynamics of the US housing market

Housing markets play a crucial role in economies and the collapse of a real-estate bubble usually destabilizes the financial system and causes economic recessions. We investigate the systemic risk and spatiotemporal dynamics of the US housing market (1975-2011) at the state level based on the Random Matrix Theory (RMT). We identify rich economic information in the largest eigenvalues deviating from RMT predictions and unveil that the component signs of the eigenvectors contain either geographical information or the extent of differences in house price growth rates or both. Our results show that the US housing market experienced six different regimes, which is consistent with the evolution of state clusters identified by the box clustering algorithm and the consensus clustering algorithm on the partial correlation matrices. Our analysis uncovers that dramatic increases in the systemic risk are usually accompanied with regime shifts, which provides a means of early detection of housing bubbles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Assortativity Decreases the Robustness of Interdependent Networks

It was recently recognized that interdependencies among different networks can play a crucial role in triggering cascading failures and hence system-wide disasters. A recent model shows how pairs of interdependent networks can exhibit an abrupt percolation transition as failures accumulate. We report on the effects of topology on failure propagation for a model system consisting of two interdependent networks. We find that the internal node correlations in each of the two interdependent networks significantly changes the critical density of failures that triggers the total disruption of the two-network system. Specifically, we find that the assortativity (i.e. the likelihood of nodes with similar degree to be connected) within a single network decreases the robustness of the entire system. The results of this study on the influence of assortativity may provide insights into ways of improving the robustness of network architecture, and thus enhances the level of protection of critical infrastructures.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Earthquake networks based on similar activity patterns

Earthquakes are a complex spatiotemporal phenomenon, the underlying mechanism for which is still not fully understood despite decades of research and analysis. We propose and develop a network approach to earthquake events. In this network, a node represents a spatial location while a link between two nodes represents similar activity patterns in the two different locations. The strength of a link is proportional to the strength of the cross-correlation in activities of two nodes joined by the link. We apply our network approach to a Japanese earthquake catalog spanning the 14-year period 1985-1998. We find strong links representing large correlations between patterns in locations separated by more than 1000 km, corroborating prior observations that earthquake interactions have no characteristic length scale. We find network characteristics not attributable to chance alone, including a large number of network links, high node assortativity, and strong stability over time.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effect of pressure on the anomalous response functions of a confined water monolayer at low temperature

We study a coarse-grained model for a water monolayer that cannot crystallize due to the presence of confining interfaces, such as protein powders or inorganic surfaces. Using both Monte Carlo simulations and mean field calculations, we calculate three response functions: the isobaric specific heat $C_P$, the isothermal compressibility $K_T$, and the isobaric thermal expansivity $α_P$. At low temperature $T$, we find two distinct maxima in $C_P$, $K_T$ and $|α_P|$, all converging toward a liquid-liquid critical point (LLCP) with increasing pressure $P$. We show that the maximum in $C_P$ at higher $T$ is due to the fluctuations of hydrogen (H) bond formation and that the second maximum at lower $T$ is due to the cooperativity among the H bonds. We discuss a similar effect in $K_T$ and $|α_P|$. If this cooperativity were not taken into account, both the lower-$T$ maximum and the LLCP would disappear. However, comparison with recent experiments on water hydrating protein powders provides evidence for the existence of the lower-$T$ maximum, supporting the hypothesized LLCP at positive $P$ and finite $T$. The model also predicts that when $P$ moves closer to the critical $P$ the $C_P$ maxima move closer in $T$ until they merge at the LLCP. Considering that other scenarios for water are thermodynamically possible, we discuss how an experimental measurement of the changing separation in $T$ between the two maxima of $C_P$ as $P$ increases could determine the best scenario for describing water.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Languages cool as they expand: Allometric scaling and the decreasing need for new words

We analyze the occurrence frequencies of over 15 million words recorded in millions of books published during the past two centuries in seven different languages. For all languages and chronological subsets of the data we confirm that two scaling regimes characterize the word frequency distributions, with only the more common words obeying the classic Zipf law. Using corpora of unprecedented size, we test the allometric scaling relation between the corpus size and the vocabulary size of growing languages to demonstrate a decreasing marginal need for new words, a feature that is likely related to the underlying correlations between words. We calculate the annual growth fluctuations of word use which has a decreasing trend as the corpus size increases, indicating a slowdown in linguistic evolution following language expansion. This "cooling pattern" forms the basis of a third statistical regularity, which unlike the Zipf and the Heaps law, is dynamical in nature.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nanoscale Dynamics of Phase Flipping in Water near its Hypothesized Liquid-Liquid Critical Point

Achieving a coherent understanding of the many thermodynamic and dynamic anomalies of water is among the most important unsolved puzzles in physics, chemistry, and biology. One hypothesized explanation imagines the existence of a line of first order phase transitions separating two liquid phases and terminating at a novel "liquid-liquid" critical point in a region of low temperature ($T \approx 250 \rm{K}$) and high pressure ($P \approx 200 \rm{MPa}$). Here we analyze a common model of water, the ST2 model, and find that the entire system flips between liquid states of high and low density. Further, we find that in the critical region crystallites melt on a time scale of nanoseconds. We perform a finite-size scaling analysis that accurately locates both the liquid-liquid coexistence line and its associated liquid-liquid critical point.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-consensus opinion models on complex networks

We focus on non-consensus opinion models in which above a certain threshold two opinions coexist in a stable relationship. We revisit and extend the non-consensus opinion (NCO) model introduced by Shao. We generalize the NCO model by adding a weight factor W to individual's own opinion when determining its future opinion (NCOW model). We find that as W increases the minority opinion holders tend to form stable clusters with a smaller initial minority fraction compared to the NCO model. We also revisit another non-consensus opinion, the inflexible contrarian opinion (ICO) model, which introduces inflexible contrarians to model a competition between two opinions in the steady state. In the ICO model, the inflexible contrarians effectively decrease the size of the largest cluster of the rival opinion. All of the above models have previously been explored in terms of a single network. However opinions propagate not only within single networks but also between networks, we study here the opinion dynamics in coupled networks. We apply the NCO rule on each individual network and the global majority rule on interdependent pairs. We find that the interdependent links effectively force the system from a second order phase transition, which is characteristic of the NCO model on a single network, to a hybrid phase transition, i.e., a mix of second-order and abrupt jump-like transitions that ultimately becomes, as we increase the percentage of interdependent agents, a pure abrupt transition. We conclude that for the NCO model on coupled networks, interactions through interdependent links could push the non-consensus opinion type model to a consensus opinion type model, which mimics the reality that increased mass communication causes people to hold opinions that are increasingly similar.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Persistence and Uncertainty in the Academic Career

Understanding how institutional changes within academia may affect the overall potential of science requires a better quantitative representation of how careers evolve over time. Since knowledge spillovers, cumulative advantage, competition, and collaboration are distinctive features of the academic profession, both the employment relationship and the procedures for assigning recognition and allocating funding should be designed to account for these factors. We study the annual production n_{i}(t) of a given scientist i by analyzing longitudinal career data for 200 leading scientists and 100 assistant professors from the physics community. We compare our results with 21,156 sports careers. Our empirical analysis of individual productivity dynamics shows that (i) there are increasing returns for the top individuals within the competitive cohort, and that (ii) the distribution of production growth is a leptokurtic "tent-shaped" distribution that is remarkably symmetric. Our methodology is general, and we speculate that similar features appear in other disciplines where academic publication is essential and collaboration is a key feature. We introduce a model of proportional growth which reproduces these two observations, and additionally accounts for the significantly right-skewed distributions of career longevity and achievement in science. Using this theoretical model, we show that short-term contracts can amplify the effects of competition and uncertainty making careers more vulnerable to early termination, not necessarily due to lack of individual talent and persistence, but because of random negative production shocks. We show that fluctuations in scientific production are quantitatively related to a scientist's collaboration radius and team efficiency.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Preferential Attachment in the Interaction between Dynamically Generated Interdependent Networks

We generalize the scale-free network model of Barabàsi and Albert [Science 286, 509 (1999)] by proposing a class of stochastic models for scale-free interdependent networks in which interdependent nodes are not randomly connected but rather are connected via preferential attachment (PA). Each network grows through the continuous addition of new nodes, and new nodes in each network attach preferentially and simultaneously to (a) well-connected nodes within the same network and (b) well-connected nodes in other networks. We present analytic solutions for the power-law exponents as functions of the number of links both between networks and within networks. We show that a cross-clustering coefficient vs. size of network $N$ follows a power law. We illustrate the models using selected examples from the Internet and finance.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Robustness of onion-like correlated networks against targeted attacks

Recently, it was found by Schneider et al. [Proc. Natl. Acad. Sci. USA, 108, 3838 (2011)], using simulations, that scale-free networks with "onion structure" are very robust against targeted high degree attacks. The onion structure is a network where nodes with almost the same degree are connected. Motivated by this work, we propose and analyze, based on analytical considerations, an onion-like candidate for a nearly optimal structure against simultaneous random and targeted high degree node attacks. The nearly optimal structure can be viewed as a hierarchically interconnected random regular graphs, the degrees and populations of which are specified by the degree distribution. This network structure exhibits an extremely assortative degree-degree correlation and has a close relationship to the "onion structure." After deriving a set of exact expressions that enable us to calculate the critical percolation threshold and the giant component of a correlated network for an arbitrary type of node removal, we apply the theory to the cases of random scale-free networks that are highly vulnerable against targeted high degree node removal. Our results show that this vulnerability can be significantly reduced by implementing this onion-like type of degree-degree correlation without much undermining the almost complete robustness against random node removal. We also investigate in detail the robustness enhancement due to assortative degree-degree correlation by introducing a joint degree-degree probability matrix that interpolates between an uncorrelated network structure and the onion-like structure proposed here by tuning a single control parameter. The optimal values of the control parameter that maximize the robustness against simultaneous random and targeted attacks are also determined. Our analytical calculations are supported by numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Statistical Laws Governing Fluctuations in Word Use from Word Birth to Word Death

We analyze the dynamic properties of 10^7 words recorded in English, Spanish and Hebrew over the period 1800--2008 in order to gain insight into the coevolution of language and culture. We report language independent patterns useful as benchmarks for theoretical models of language evolution. A significantly decreasing (increasing) trend in the birth (death) rate of words indicates a recent shift in the selection laws governing word use. For new words, we observe a peak in the growth-rate fluctuations around 40 years after introduction, consistent with the typical entry time into standard dictionaries and the human generational timescale. Pronounced changes in the dynamics of language during periods of war shows that word correlations, occurring across time and between words, are largely influenced by coevolutionary social, technological, and political factors. We quantify cultural memory by analyzing the long-term correlations in the use of individual words using detrended fluctuation analysis.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The competitiveness versus the wealth of a country

Politicians world-wide frequently promise a better life for their citizens. We find that the probability that a country will increase its {\it per capita} GDP ({\it gdp}) rank within a decade follows an exponential distribution with decay constant $λ= 0.12$. We use the Corruption Perceptions Index (CPI) and the Global Competitiveness Index (GCI) and find that the distribution of change in CPI (GCI) rank follows exponential functions with approximately the same exponent as $λ$, suggesting that the dynamics of {\it gdp}, CPI, and GCI may share the same origin. Using the GCI, we develop a new measure, which we call relative competitiveness, to evaluate an economy's competitiveness relative to its {\it gdp}. For all European and EU countries during the 2008-2011 economic downturn we find that the drop in {\it gdp} in more competitive countries relative to {\it gdp} was substantially smaller than in relatively less competitive countries, which is valuable information for policymakers.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Detrending career statistics in professional baseball: Accounting for the steroids era and beyond

There is a long standing debate over how to objectively compare the career achievements of professional athletes from different historical eras. Developing an objective approach will be of particular importance over the next decade as Major League Baseball (MLB) players from the "steroids era" become eligible for Hall of Fame induction. Here we address this issue, as well as the general problem of comparing statistics from distinct eras, by detrending the seasonal statistics of professional baseball players. We detrend player statistics by normalizing achievements to seasonal averages, which accounts for changes in relative player ability resulting from both exogenous and endogenous factors, such as talent dilution from expansion, equipment and training improvements, as well as performance enhancing drugs (PED). In this paper we compare the probability density function (pdf) of detrended career statistics to the pdf of raw career statistics for five statistical categories -- hits (H), home runs (HR), runs batted in (RBI), wins (W) and strikeouts (K) -- over the 90-year period 1920-2009. We find that the functional form of these pdfs are stationary under detrending. This stationarity implies that the statistical regularity observed in the right-skewed distributions for longevity and success in professional baseball arises from both the wide range of intrinsic talent among athletes and the underlying nature of competition. Using this simple detrending technique, we examine the top 50 all-time careers for H, HR, RBI, W and K. We fit the pdfs for career success by the Gamma distribution in order to calculate objective benchmarks based on extreme statistics which can be used for the identification of extraordinary careers.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Financial factor influence on scaling and memory of trading volume in stock market

We study the daily trading volume volatility of 17,197 stocks in the U.S. stock markets during the period 1989--2008 and analyze the time return intervals $τ$ between volume volatilities above a given threshold q. For different thresholds q, the probability density function P_q(τ) scales with mean interval <τ> as P_q(τ)=<τ>^{-1}f(τ/<τ>) and the tails of the scaling function can be well approximated by a power-law f(x)~x^{-γ}. We also study the relation between the form of the distribution function P_q(τ) and several financial factors: stock lifetime, market capitalization, volume, and trading value. We find a systematic tendency of P_q(τ) associated with these factors, suggesting a multi-scaling feature in the volume return intervals. We analyze the conditional probability P_q(τ|τ_0) for $τ$ following a certain interval τ_0, and find that P_q(τ|τ_0) depends on τ_0 such that immediately following a short/long return interval a second short/long return interval tends to occur. We also find indications that there is a long-term correlation in the daily volume volatility. We compare our results to those found earlier for price volatility.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Large decrease of fluctuations for supercooled water in hydrophobic nanoconfinement

We study by Monte Carlo simulations a coarse-grained model of a water layer confined in a fixed disordered matrix of hydrophobic nanoparticles at different particle concentrations c. For c=0 (bulk water) we find a first order liquid-liquid phase transition (LLPT). For c>0 our simulations are consistent with a LLPT line ending in two critical points at low and high pressure P. For c=25% we find a dramatic decrease of compressibility K_T, thermal expansion coefficient α_P, specific heat C_P at high P and low temperature T. Surprisingly, the effect is present also for c as low as 2.4%. We conclude that the higher the confinement, the lower the thermodynamic fluctuations and the more difficult it will be to detect the LLPT.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

More than one dynamic crossover in protein hydration water

Studies of liquid water in its supercooled region have led to many insights into the structure and behavior of water. While bulk water freezes at its homogeneous nucleation temperature of approximately 235 K, for protein hydration water, the binding of water molecules to the protein avoids crystallization. Here we study the dynamics of the hydrogen bond (HB) network of a percolating layer of water molecules, comparing measurements of a hydrated globular protein with the results of a coarse-grained model that has been shown to successfully reproduce the properties of hydration water. With dielectric spectroscopy we measure the temperature dependence of the relaxation time of protons charge fluctuations. These fluctuations are associated to the dynamics of the HB network of water molecules adsorbed on the protein surface. With Monte Carlo (MC) simulations and mean--field (MF) calculations we study the dynamics and thermodynamics of the model. In both experimental and model analyses we find two dynamic crossovers: (i) one at about 252 K, and (ii) one at about 181 K. The agreement of the experiments with the model allows us to relate the two crossovers to the presence of two specific heat maxima at ambient pressure. The first is due to fluctuations in the HB formation, and the second, at lower temperature, is due to the cooperative reordering of the HB network.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantifying and Modeling Long-Range Cross-Correlations in Multiple Time Series with Applications to World Stock Indices

We propose a modified time lag random matrix theory in order to study time lag cross-correlations in multiple time series. We apply the method to 48 world indices, one for each of 48 different countries. We find long-range power-law cross-correlations in the absolute values of returns that quantify risk, and find that they decay much more slowly than cross-correlations between the returns. The magnitude of the cross-correlations constitute "bad news" for international investment managers who may believe that risk is reduced by diversifying across countries. We find that when a market shock is transmitted around the world, the risk decays very slowly. We explain these time lag cross-correlations by introducing a global factor model (GFM) in which all index returns fluctuate in response to a single global factor. For each pair of individual time series of returns, the cross-correlations between returns (or magnitudes) can be modeled with the auto-correlations of the global factor returns (or magnitudes). We estimate the global factor using principal component analysis, which minimizes the variance of the residuals after removing the global trend. Using random matrix theory, a significant fraction of the world index cross-correlations can be explained by the global factor, which supports the utility of the GFM. We demonstrate applications of the GFM in forecasting risks at the world level, and in finding uncorrelated individual indices. We find 10 indices are practically uncorrelated with the global factor and with the remainder of the world indices, which is relevant information for world managers in reducing their portfolio risk. Finally, we argue that this general method can be applied to a wide range of phenomena in which time series are measured, ranging from seismology and physiology to atmospheric geophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantitative and empirical demonstration of the Matthew effect in a study of career longevity

The Matthew effect refers to the adage written some two-thousand years ago in the Gospel of St. Matthew: "For to all those who have, more will be given." Even two millennia later, this idiom is used by sociologists to qualitatively describe the dynamics of individual progress and the interplay between status and reward. Quantitative studies of professional careers are traditionally limited by the difficulty in measuring progress and the lack of data on individual careers. However, in some professions, there are well-defined metrics that quantify career longevity, success, and prowess, which together contribute to the overall success rating for an individual employee. Here we demonstrate testable evidence of the age-old Matthew "rich get richer" effect, wherein the longevity and past success of an individual lead to a cumulative advantage in further developing his/her career. We develop an exactly solvable stochastic career progress model that quantitatively incorporates the Matthew effect, and validate our model predictions for several competitive professions. We test our model on the careers of 400,000 scientists using data from six high-impact journals, and further confirm our findings by testing the model on the careers of more than 20,000 athletes in four sports leagues. Our model highlights the importance of early career development, showing that many careers are stunted by the relative disadvantage associated with inexperience.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Scaling of Seismic Memory with Earthquake Size

It has been observed that the earthquake events possess short-term memory, i.e. that events occurring in a particular location are dependent on the short history of that location. We conduct an analysis to see whether real-time earthquake data also possess long-term memory and, if so, whether such autocorrelations depend on the size of earthquakes within close spatiotemporal proximity. We analyze the seismic waveform database recorded by 64 stations in Japan, including the 2011 "Great East Japan Earthquake", one of the five most powerful earthquakes ever recorded which resulted in a tsunami and devastating nuclear accidents. We explore the question of seismic memory through use of mean conditional intervals and detrended fluctuation analysis (DFA). We find that the waveform sign series show long-range power-law anticorrelations while the interval series show long-range power-law correlations. We find size-dependence in earthquake auto-correlations---as earthquake size increases, both of these correlation behaviors strengthen. We also find that the DFA scaling exponent $α$ has no dependence on earthquake hypocenter depth or epicentral distance.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Statistical regularities in the rank-citation profile of scientists

Recent "science of science" research shows that scientific impact measures for journals and individual articles have quantifiable regularities across both time and discipline. However, little is known about the scientific impact distribution at the scale of an individual scientist. We analyze the aggregate scientific production and impact of individual careers using the rank-citation profile c_{i}(r) of 200 distinguished professors and 100 assistant professors. For the entire range of paper rank r, we fit each c_{i}(r) to a common distribution function that is parameterized by two scaling exponents. Since two scientists with equivalent Hirsch h-index can have significantly different c_{i}(r) profiles, our results demonstrate the utility of the β_{i} scaling parameter in conjunction with h_{i} for quantifying individual publication impact. We show that the total number of citations C_{i} tallied from a scientist's N_{i} papers scales as C_{i} \sim h_{i}^{1+β_{i}}. Such statistical regularities in the input-output patterns of scientists can be used as benchmarks for theoretical models of career progress.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Strategy of Competition between Two Groups based on a Contrarian Opinion Model

We introduce a contrarian opinion (CO) model in which a fraction p of contrarians within a group holds a strong opinion opposite to the opinion held by the rest of the group. At the initial stage, stable clusters of two opinions, A and B exist. Then we introduce contrarians which hold a strong B opinion into the opinion A group. Through their interactions, the contrarians are able to decrease the size of the largest A opinion cluster, and even destroy it. We see this kind of method in operation, e.g when companies send free new products to potential customers in order to convince them to adopt the product and influence others. We study the CO model, using two different strategies, on both ER and scale-free networks. In strategy I, the contrarians are positioned at random. In strategy II, the contrarians are chosen to be the highest degrees nodes. We find that for both strategies the size of the largest A cluster decreases to zero as p increases as in a phase transition. At a critical threshold value p_c the system undergoes a second-order phase transition that belongs to the same universality class of mean field percolation. We find that even for an ER type model, where the degrees of the nodes are not so distinct, strategy II is significantly more effctive in reducing the size of the largest A opinion cluster and, at very small values of p, the largest A opinion cluster is destroyed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Random Quadratic Assignment Problem

Optimal assignment of classes to classrooms \cite{dickey}, design of DNA microarrays \cite{carvalho}, cross species gene analysis \cite{kolar}, creation of hospital layouts cite{elshafei}, and assignment of components to locations on circuit boards \cite{steinberg} are a few of the many problems which have been formulated as a quadratic assignment problem (QAP). Originally formulated in 1957, the QAP is one of the most difficult of all combinatorial optimization problems. Here, we use statistical mechanical methods to study the asymptotic behavior of problems in which the entries of at least one of the two matrices that specify the problem are chosen from a random distribution $P$. Surprisingly, this case has not been studied before using statistical methods despite the fact that the QAP was first proposed over 50 years ago \cite{Koopmans}. We find simple forms for $C_{\rm min}$ and $C_{\rm max}$, the costs of the minimal and maximum solutions respectively. Notable features of our results are the symmetry of the results for $C_{\rm min}$ and $C_{\rm max}$ and the dependence on $P$ only through its mean and standard deviation, independent of the details of $P$. After the asymptotic cost is determined for a given QAP problem, one can straightforwardly calculate the asymptotic cost of a QAP problem specified with a different random distribution $P$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Bankruptcy risk model and empirical tests

We analyze the size dependence and temporal stability of firm bankruptcy risk in the US economy by applying Zipf scaling techniques. We focus on a single risk factor-the debt-to-asset ratio R-in order to study the stability of the Zipf distribution of R over time. We find that the Zipf exponent increases during market crashes, implying that firms go bankrupt with larger values of R. Based on the Zipf analysis, we employ Bayes's theorem and relate the conditional probability that a bankrupt firm has a ratio R with the conditional probability of bankruptcy for a firm with a given R value. For 2,737 bankrupt firms, we demonstrate size dependence in assets change during the bankruptcy proceedings. Prepetition firm assets and petition firm assets follow Zipf distributions but with different exponents, meaning that firms with smaller assets adjust their assets more than firms with larger assets during the bankruptcy process. We compare bankrupt firms with nonbankrupt firms by analyzing the assets and liabilities of two large subsets of the US economy: 2,545 Nasdaq members and 1,680 New York Stock Exchange (NYSE) members. We find that both assets and liabilities follow a Pareto distribution. The finding is not a trivial consequence of the Zipf scaling relationship of firm size quantified by employees-although the market capitalization of Nasdaq stocks follows a Pareto distribution, the same distribution does not describe NYSE stocks. We propose a coupled Simon model that simultaneously evolves both assets and debt with the possibility of bankruptcy, and we also consider the possibility of firm mergers.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cross-correlations between volume change and price change

In finance, one usually deals not with prices but with growth rates $R$, defined as the difference in logarithm between two consecutive prices. Here we consider not the trading volume, but rather the volume growth rate $\tilde R$, the difference in logarithm between two consecutive values of trading volume. To this end, we use several methods to analyze the properties of volume changes $|\tilde R|$, and their relationship to price changes $|R|$. We analyze $14,981$ daily recordings of the S\&P 500 index over the 59-year period 1950--2009, and find power-law {\it cross-correlations\/} between $|R|$ and $|\tilde R|$ using detrended cross-correlation analysis (DCCA). We introduce a joint stochastic process that models these cross-correlations. Motivated by the relationship between $| R|$ and $|\tilde R|$, we estimate the tail exponent ${\tildeα}$ of the probability density function $P(|\tilde R|) \sim |\tilde R|^{-1 -\tildeα}$ for both the S\&P 500 index as well as the collection of 1819 constituents of the New York Stock Exchange Composite index on 17 July 2009. As a new method to estimate $\tildeα$, we calculate the time intervals $τ_q$ between events where $\tilde R>q$. We demonstrate that $\barτ_q$, the average of $τ_q$, obeys $\bar τ_q \sim q^{\tildeα}$. We find $\tilde α\approx 3$. Furthermore, by aggregating all $τ_q$ values of 28 global financial indices, we also observe an approximate inverse cubic law.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Market dynamics immediately before and after financial shocks: quantifying the Omori, productivity and Bath laws

We study the cascading dynamics immediately before and immediately after 219 market shocks. We define the time of a market shock T_{c} to be the time for which the market volatility V(T_{c}) has a peak that exceeds a predetermined threshold. The cascade of high volatility "aftershocks" triggered by the "main shock" is quantitatively similar to earthquakes and solar flares, which have been described by three empirical laws --- the Omori law, the productivity law, and the Bath law. We analyze the most traded 531 stocks in U.S. markets during the two-year period 2001-2002 at the 1-minute time resolution. We find quantitative relations between (i) the "main shock" magnitude M \equiv \log V(T_{c}) occurring at the time T_{c} of each of the 219 "volatility quakes" analyzed, and (ii) the parameters quantifying the decay of volatility aftershocks as well as the volatility preshocks. We also find that stocks with larger trading activity react more strongly and more quickly to market shocks than stocks with smaller trading activity. Our findings characterize the typical volatility response conditional on M, both at the market and the individual stock scale. We argue that there is potential utility in these three statistical quantitative relations with applications in option pricing and volatility trading.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Methods for measuring the citations and productivity of scientists across time and discipline

Publication statistics are ubiquitous in the ratings of scientific achievement, with citation counts and paper tallies factoring into an individual's consideration for postdoctoral positions, junior faculty, tenure, and even visa status for international scientists. Citation statistics are designed to quantify individual career achievement, both at the level of a single publication, and over an individual's entire career. While some academic careers are defined by a few significant papers (possibly out of many), other academic careers are defined by the cumulative contribution made by the author's publications to the body of science. Several metrics have been formulated to quantify an individual's publication career, yet none of these metrics account for the dependence of citation counts and journal size on time. In this paper, we normalize publication metrics across both time and discipline in order to achieve a universal framework for analyzing and comparing scientific achievement. We study the publication careers of individual authors over the 50-year period 1958-2008 within six high-impact journals: CELL, the New England Journal of Medicine (NEJM), Nature, the Proceedings of the National Academy of Science (PNAS), Physical Review Letters (PRL), and Science. In comparing the achievement of authors within each journal, we uncover quantifiable statistical regularity in the probability density function (pdf) of scientific achievement across both time and discipline. The universal distribution of career success within these arenas for publication raises the possibility that a fundamental driving force underlying scientific achievement is the competitive nature of scientific advancement.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantitative law describing market dynamics before and after interest-rate change

We study the behavior of U.S. markets both before and after U.S. Federal Open Market Committee (FOMC) meetings, and show that the announcement of a U.S. Federal Reserve rate change causes a financial shock, where the dynamics after the announcement is described by an analogue of the Omori earthquake law. We quantify the rate n(t) of aftershocks following an interest rate change at time T, and find power-law decay which scales as n(t-T) (t-T)^-$Ω$, with $Ω$ positive. Surprisingly, we find that the same law describes the rate n'(|t-T|) of "pre-shocks" before the interest rate change at time T. This is the first study to quantitatively relate the size of the market response to the news which caused the shock and to uncover the presence of quantifiable preshocks. We demonstrate that the news associated with interest rate change is responsible for causing both the anticipation before the announcement and the surprise after the announcement. We estimate the magnitude of financial news using the relative difference between the U. S. Treasury Bill and the Federal Funds Effective rate. Our results are consistent with the "sign effect," in which "bad news" has a larger impact than "good news." Furthermore, we observe significant volatility aftershocks, confirming a "market underreaction" that lasts at least 1 trading day.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Statistical Properties of Cross-Correlation in the Korean Stock Market

We investigate the statistical properties of the correlation matrix between individual stocks traded in the Korean stock market using the random matrix theory (RMT) and observe how these affect the portfolio weights in the Markowitz portfolio theory. We find that the distribution of the correlation matrix is positively skewed and changes over time. We find that the eigenvalue distribution of original correlation matrix deviates from the eigenvalues predicted by the RMT, and the largest eigenvalue is 52 times larger than the maximum value among the eigenvalues predicted by the RMT. The $β_{473}$ coefficient, which reflect the largest eigenvalue property, is 0.8, while one of the eigenvalues in the RMT is approximately zero. Notably, we show that the entropy function $E(σ)$ with the portfolio risk $σ$ for the original and filtered correlation matrices are consistent with a power-law function, $E(σ) \sim σ^{-γ}$, with the exponent $γ\sim 2.92$ and those for Asian currency crisis decreases significantly.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Disconnected Glass-Glass Transitions and Diffusion Anomalies in a Model with two Repulsive Length Scales

Building on mode-coupling-theory calculations, we report a novel scenario for multiple glass transitions in a purely repulsive spherical potential: the square-shoulder. The liquid-glass transition lines exhibit both melting by cooling and melting by compression as well as associated diffusion anomalies, similar to the ones observed in water. Differently from all previously investigated models, here for small shoulder widths a glass-glass line is found that is disconnected from the liquid phase. Upon increasing the shoulder width such a glass-glass line merges with the liquid-glass transition lines, featuring two distinct endpoint singularities that give rise to logarithmic decays in the dynamics. These findings can be explained analytically by the interplay of different repulsive length scales.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamic networks and directed percolation

We introduce a model for dynamic networks, where the links or the strengths of the links change over time. We solve the model by mapping dynamic networks to the problem of directed percolation, where the direction corresponds to the evolution of the network in time. We show that the dynamic network undergoes a percolation phase transition at a critical concentration $p_c$, which decreases with the rate $r$ at which the network links are changed. The behavior near criticality is universal and independent of $r$. We find fundamental network laws are changed. (i) For Erdős-Rényi networks we find that the size of the giant component at criticality scales with the network size $N$ for all values of $r$, rather than as $N^{2/3}$. (ii) In the presence of a broad distribution of disorder, the optimal path length between two nodes in a dynamic network scales as $N^{1/2}$, compared to $N^{1/3}$ in a static network.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Communication Bottlenecks in Scale-Free Networks

We consider the effects of network topology on the optimality of packet routing quantified by $γ_c$, the rate of packet insertion beyond which congestion and queue growth occurs. The key result of this paper is to show that for any network, there exists an absolute upper bound, expressed in terms of vertex separators, for the scaling of $γ_c$ with network size $N$, irrespective of the routing algorithm used. We then derive an estimate to this upper bound for scale-free networks, and introduce a novel static routing protocol which is superior to shortest path routing under intense packet insertion rates.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Relationship between Fragility, Diffusive Directions and Energy Barriers in a Supercooled Liquid

An analysis of diffusion in a supercooled liquid based solely in the density of diffusive directions and the value of energy barriers shows how the potential energy landscape (PEL) approach is capable of explaining the $α$ and $β$ relaxations and the fragility of a glassy system. We find that the $β$ relaxation is directly related to the search for diffusive directions. Our analysis shows how in strong liquids diffusion is mainly energy activated, and how in fragile liquids the diffusion is governed by the density of diffusive directions. We describe the fragile-to-strong crossover as a change in the topography of the PEL sampled by the system at a certain crossover temperature $T_\times$.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Transport of multiple users in complex networks

We study the transport properties of model networks such as scale-free and Erdős-Rényi networks as well as a real network. 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 $Φ_{\rm SF}(G)\sim G^{-g_G}$, where $g_G=2λ-1$, and $λ$ is the decay exponent for the scale-free network degree distribution. We confirm our predictions by large scale 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 networks where the tail of the conductivity distribution decays exponentially. We develop a simple physical picture of the transport to account for the results. We study another model for transport, the \emph{max-flow} model, where conductance is defined as the number of link-independent paths between the two nodes, 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. We then extend our study to the case of multiple sources, where the transport is define between two \emph{groups} of nodes. We find a fundamental difference between the two forms of flow when considering the quality of the transport with respect to the number of sources, and find an optimal number of sources, or users, for the max-flow case. A qualitative (and partially quantitative) explanation is also given.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Universal behavior of optimal paths in weighted networks with general disorder

We study the statistics of the optimal path in both random and scale free networks, where weights $w$ are taken from a general distribution $P(w)$. We find that different types of disorder lead to the same universal behavior. Specifically, we find that a single parameter ($S \equiv AL^{-1/ν}$ for $d$-dimensional lattices, and $S\equiv AN^{-1/3}$ for random networks) determines the distributions of the optimal path length, including both strong and weak disorder regimes. Here $ν$ is the percolation connectivity exponent, and $A$ depends on the percolation threshold and $P(w)$. For $P(w)$ uniform, Poisson or Gaussian the crossover from weak to strong does not occur, and only weak disorder exists.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Relation between Rotational and Translational Dynamic Heterogeneities in Water

We use molecular dynamics simulations to probe the rotational dynamics of the SPC/E model of water for a range of temperatures down to 200 K, 13 K above to the mode coupling temperature. We find that rotational dynamics is spatially heterogeneous, i.e., there are clusters of molecules that rotate significantly more than the average for a given time interval, and we study the size and the temporal behavior of these clusters. We find that the position of a rotational heterogeneity is strongly correlated with the position of a translational heterogeneity, and that the fraction of molecules belonging to both kinds of heterogeneities increases with decreasing temperature. We further find that although the two types of heterogeneities are not identical, they are related to the same physical picture.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Levels of Complexity in Scale-Invariant Neural Signals

Many physical and physiological signals exhibit complex scale-invariant features characterized by $1/f$ scaling and long-range power-law correlations, suggesting a possibly common control mechanism. Specifically, it has been suggested that dynamical processes influenced by inputs and feedback on multiple time scales may be sufficient to give rise to $1/f$ scaling and scale invariance. Two examples of physiologic signals that are the output of hierarchical, multi-scale physiologic systems under neural control are the human heartbeat and human gait. Here we show that while both cardiac interbeat interval and gait interstride interval time series under healthy conditions have comparable $1/f$ scaling, they still may belong to different complexity classes. Our analysis of the magnitude series correlations and multifractal scaling exponents of the fluctuations in these two signals demonstrates that in contrast with the nonlinear multifractal behavior found in healthy heartbeat dynamics, gait time series exhibit less complex, close to monofractal behavior and a low degree of nonlinearity. These findings are of interest because they underscore the limitations of traditional two-point correlation methods in fully characterizing physiologic and physical dynamics. In addition, these results suggest that different mechanisms of control may be responsible for varying levels of complexity observed in physiological systems under neural regulation and in physical systems that possess similar $1/f$ scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The Glass Transition Temperature of Water: A Simulation Study

We report a computer simulation study of the glass transition for water. To mimic the difference between standard and hyperquenched glass, we generate glassy configurations with different cooling rates and calculate the $T$ dependence of the specific heat on heating. The absence of crystallization phenomena allows us, for properly annealed samples, to detect in the specific heat the simultaneous presence of a weak pre-peak (``shadow transition''), and an intense glass transition peak at higher temperature. We discuss the implications for the currently debated value of the glass transition temperature of water. We also compare our simulation results with the Tool-Narayanaswamy-Moynihan phenomenological model.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Intra-molecular coupling as a mechanism for a liquid-liquid phase transition

We study a model for water with a tunable intra-molecular interaction $J_σ$, using mean field theory and off-lattice Monte Carlo simulations. For all $J_σ\geq 0$, the model displays a temperature of maximum density.For a finite intra-molecular interaction $J_σ> 0$,our calculations support the presence of a liquid-liquid phase transition with a possible liquid-liquid critical point for water, likely pre-empted by inevitable freezing. For J=0 the liquid-liquid critical point disappears at T=0.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Noise Effects on the Complex Patterns of Abnormal Heartbeats

Patients at high risk for sudden death often exhibit complex heart rhythms in which abnormal heartbeats are interspersed with normal heartbeats. We analyze such a complex rhythm in a single patient over a 12-hour period and show that the rhythm can be described by a theoretical model consisting of two interacting oscillators with stochastic elements. By varying the magnitude of the noise, we show that for an intermediate level of noise, the model gives best agreement with key statistical features of the dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Classes of behavior of small-world networks

Small-world networks are the focus of recent interest because they appear to circumvent many of the limitations of either random networks or regular lattices as frameworks for the study of interaction networks of complex systems. Here, we report an empirical study of the statistical properties of a variety of diverse real-world networks. We present evidence of the occurrence of three classes of small-world networks: (a) scale-free networks, characterized by a vertex connectivity distribution that decays as a power law; (b) broad-scale networks, characterized by a connectivity distribution that has a power-law regime followed by a sharp cut-off; (c) single-scale networks, characterized by a connectivity distribution with a fast decaying tail. Moreover, we note for the classes of broad-scale and single-scale networks that there are constraints limiting the addition of new links. Our results suggest that the nature of such constraints may be the controlling factor for the emergence of different classes of networks.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Configurational Entropy and Diffusivity of Supercooled Water

We calculate the configurational entropy S_conf for the SPC/E model of water for state points covering a large region of the (T,rho) plane. We find that (i) the (T,rho) dependence of S_conf correlates with the diffusion constant and (ii) that the line of maxima in S_conf tracks the line of density maxima. Our simulation data indicate that the dynamics are strongly influenced by S_conf even above the mode-coupling temperature T_MCT(rho).

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Ivory Tower Universities and Competitive Business Firms

There is nowadays considerable interest on ways to quantify the dynamics of research activities, in part due to recent changes in research and development (R&D) funding. Here, we seek to quantify and analyze university research activities, and compare their growth dynamics with those of business firms. Specifically, we analyze five distinct databases, the largest of which is a National Science Foundation database of the R&D expenditures for science and engineering of 719 United States (US) universities for the 17-year period 1979--1995. We find that the distribution of growth rates displays a ``universal'' form that does not depend on the size of the university or on the measure of size used, and that the width of this distribution decays with size as a power law. Our findings are quantitatively similar to those independently uncovered for business firms, and consistent with the hypothesis that the growth dynamics of complex organizations may be governed by universal mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Multifractality in Human Heartbeat Dynamics

Recent evidence suggests that physiological signals under healthy conditions may have a fractal temporal structure. We investigate the possibility that time series generated by certain physiological control systems may be members of a special class of complex processes, termed multifractal, which require a large number of exponents to characterize their scaling properties. We report on evidence for multifractality in a biological dynamical system --- the healthy human heartbeat. Further, we show that the multifractal character and nonlinear properties of the healthy heart rate are encoded in the Fourier phases. We uncover a loss of multifractality for a life-threatening condition, congestive heart failure.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Sleep-Wake Differences in Scaling Behavior of the Human Heartbeat: Analysis of Terrestrial and Long-Term Space Flight Data

We compare scaling properties of the cardiac dynamics during sleep and wake periods for healthy individuals, cosmonauts during orbital flight, and subjects with severe heart disease. For all three groups, we find a greater degree of anticorrelation in the heartbeat fluctuations during sleep compared to wake periods. The sleep-wake difference in the scaling exponents for the three groups is comparable to the difference between healthy and diseased individuals. The observed scaling differences are not accounted for simply by different levels of activity, but appear related to intrinsic changes in the neuroautonomic control of the heartbeat.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Liquid State Anomalies for the Stell-Hemmer Core-Softened Potential

We study the Stell-Hemmer potential using both analytic (exact $1d$ and approximate $2d$) solutions and numerical $2d$ simulations. We observe in the liquid phase an anomalous decrease in specific volume and isothermal compressibility upon heating, and an anomalous increase in the diffusion coefficient with pressure. We relate the anomalies to the existence of two different local structures in the liquid phase. Our results are consistent with the possibility of a low temperature/high pressure liquid-liquid phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Modeling of Financial Data: Comparison of the Truncated Lévy Flight and the ARCH(1) and GARCH(1,1) processes

We compare our results on empirical analysis of financial data with simulations of two stochastic models of the dynamics of stock market prices. The two models are (i) the truncated Lévy flight recently introduced by us and (ii) the ARCH(1) and GARCH(1,1) processes. We find that the TLF well describes the scaling and its breakdown observed in empirical data, while it is not able to properly describe the fluctuations of volatility empirically detected. The ARCH(1) and GARCH(1,1) models are able to describe the probability density function of price changes at a given time horizon, but both fail to describe the scaling properties of the PDFs for short time horizons.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Correlations in Economic Time Series

The correlation function of a financial index of the New York stock exchange, the S&P 500, is analyzed at 1 min intervals over the 13-year period, Jan 84 -- Dec 96. We quantify the correlations of the absolute values of the index increment. We find that these correlations can be described by two different power laws with a crossover time t_\times\approx 600 min. Detrended fluctuation analysis gives exponents $α_1=0.66$ and $α_2=0.93$ for $t<t_\times$ and $t>t_\times$ respectively. Power spectrum analysis gives corresponding exponents $β_1=0.31$ and $β_2=0.90$ for $f>f_\times$ and $f< f_\times$ respectively.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Sandpile Model on Sierpinski Gasket Fractal

We investigate the sandpile model on the two--dimensional Sierpinski gasket fractal. We find that the model displays novel critical behavior, and we analyze the distribution functions of avalanche sizes, lifetimes and topplings and calculate the associated critical exponents $τ= 1.51 \pm 0.04$, $α= 1.63 \pm 0.04$ and $μ= 1.36 \pm 0.04$. The avalanche size distribution shows power law behavior modulated by logarithmic oscillations which can be related to the discrete scale invariance of the underlying lattice. Such a distribution can be formally described by introducing a complex scaling exponent $τ^{*} \equiv τ+ i δ$, where the real part $τ$ corresponds to the power law and the imaginary part $δ$ is related to the period of the logarithmic oscillations.

0 citations0 reviewsNo code linkedOpen access