Search

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

Many systems on our planet are known to shift abruptly and irreversibly from one state to another when they are forced across a "tipping point," such as mass extinctions in ecological networks, cascading failures in infrastructure systems, and social convention changes in human and animal networks. Such a regime shift demonstrates a system's resilience that characterizes the ability of a system to adjust its activity to retain its basic functionality in the face of internal disturbances or external environmental changes. In the past 50 years, attention was almost exclusively given to low dimensional systems and calibration of their resilience functions and indicators of early warning signals without considerations for the interactions between the components. Only in recent years, taking advantages of the network theory and lavish real data sets, network scientists have directed their interest to the real-world complex networked multidimensional systems and their resilience function and early warning indicators. This report is devoted to a comprehensive review of resilience function and regime shift of complex systems in different domains, such as ecology, biology, social systems and infrastructure. We cover the related research about empirical observations, experimental studies, mathematical modeling, and theoretical analysis. We also discuss some ambiguous definitions, such as robustness, resilience, and stability., Comment: Review chapter

We perform short-time Monte Carlo simulations to study the criticality of the isotropic two-state majority-vote model on cubic lattices of volume N = L 3 , with L up to 2048. We obtain the precise location of the critical point by examining the scaling properties of a new auxiliary function Ψ . We perform finite-time scaling analysis to accurately calculate the whole set of critical exponents, including the dynamical critical exponent z = 2 . 027 ( 9 ) , and the initial slip exponent θ = 0 . 1081 ( 1 ) . Our results indicate that the majority-vote model in three dimensions belongs to the same universality class of the three-dimensional Ising model.

Margin trading in which investors purchase shares with money borrowed from brokers is blamed to be a major cause of the 2015 Chinese stock market crash. We propose a cascading failure model and examine how an increase in margin trading increases share price vulnerability. The model is based on a bipartite graph of investors and shares that includes four margin trading factors, (i) initial margin $k$, (ii) minimum maintenance $r$, (iii) volatility $v$, and (iv) diversity $s$. We use our model to simulate margin trading and observe how the share prices are affected by these four factors. The experimental results indicate that a stock market can be either vulnerable or stable. A stock market is vulnerable when an external shock can cause a cascading failure of its share prices. It is stable when its share prices are resilient to external shocks. Furthermore, we investigate how the cascading failure of share price is affected by these four factors, and find that by increasing $v$ and $r$ or decreasing $k$ we increase the probability that the stock market will experience a phase transition from stable to vulnerable. It is also found that increasing $s$ decreases resilience and increases systematic risk. These findings could be useful to regulators supervising margin trading activities., Comment: 14 pages, 9 figures

Detection of power-law behavior and studies of scaling exponents uncover the characteristics of complexity in many real world phenomena. The complexity of financial markets has always presented challenging issues and provided interesting findings, such as the inverse cubic law in the tails of stock price fluctuation distributions. Motivated by the rise of novel digital assets based on blockchain technology, we study the distributions of cryptocurrency price fluctuations. We consider Bitcoin returns over various time intervals and from multiple digital exchanges, in order to investigate the existence of universal scaling behavior in the tails, and ascertain whether the scaling exponent supports the presence of a finite second moment. We provide empirical evidence on slowly decaying tails in the distributions of returns over multiple time intervals and different exchanges, corresponding to a power-law. We estimate the scaling exponent and find an asymptotic power-law behavior with 2 α 2 . 5 suggesting that Bitcoin returns, in addition to being more volatile, also exhibit heavier tails than stocks, which are known to be around 3. Our results also imply the existence of a finite second moment, thus providing a fundamental basis for the usage of standard financial theories and covariance-based techniques in risk management and portfolio optimization scenarios.

Despite persistent efforts in understanding the creativity of scientists over different career stages, little is known about the underlying dynamics of research topic switching that drives innovation. Here, we analyze the publication records of individual scientists, aiming to quantify their topic switching dynamics and its influence. We find that the co-citing network of papers of a scientist exhibits a clear community structure where each major community represents a research topic. Our analysis suggests that scientists have a narrow distribution of number of topics. However, researchers nowadays switch more frequently between topics than those in the early days. We also find that high switching probability in early career is associated with low overall productivity, yet with high overall productivity in latter career. Interestingly, the average citation per paper, however, is in all career stages negatively correlated with the switching probability. We propose a model that can explain the main observed features., How does a scientist’s tendency to explore a variety of topics affect their career? Here, the authors analyze scientific publication data to understand how often scientists switch topics, how topic switching has changed over time, and how it relates to research productivity.

K-core and bootstrap percolation are widely studied models that have been used to represent and understand diverse deactivation and activation processes in natural and social systems. Since these models are considerably similar, it has been suggested in recent years that they could be complementary. In this manuscript we provide a rigorous analysis that shows that for any degree and threshold distributions heterogeneous bootstrap percolation can be mapped into heterogeneous k-core percolation and vice versa, if the functionality thresholds in both processes satisfy a complementary relation. Another interesting problem in bootstrap and k-core percolation is the fraction of nodes belonging to their giant connected components P∞b and P∞c, respectively. We solve this problem analytically for arbitrary randomly connected graphs and arbitrary threshold distributions, and we show that P∞b and P∞c are not complementary. Our theoretical results coincide with computer simulations in the limit of very large graphs. In bootstrap percolation, we show that when using the branching theory to compute the size of the giant component, we must consider two different types of links, which are related to distinct spanning branches of active nodes. Fil: Di Muro, Matias Alberto. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mar del Plata. Instituto de Investigaciones Físicas de Mar del Plata. Universidad Nacional de Mar del Plata. Facultad de Ciencias Exactas y Naturales. Instituto de Investigaciones Físicas de Mar del Plata; Argentina Fil: Valdez, Lucas Daniel. Boston University; Estados Unidos Fil: Stanley, Harry Eugene. Boston University; Estados Unidos Fil: Buldyrev, Sergey V.. Yeshiva University; Estados Unidos Fil: Braunstein, Lidia Adriana. Politecnico di Milano; Italia

The growth of business firms is an example of a system of complex interacting units that resembles complex interacting systems in nature such as earthquakes. Remarkably, work in econophysics has provided evidence that the statistical properties of the growth of business firms follow the same sorts of power laws that characterize physical systems near their critical points. Given how economies change over time, whether these statistical properties are persistent, robust, and universal like those of physical systems remains an open question. Here, we show that the scaling properties of firm growth previously demonstrated for publicly-traded U.S. manufacturing firms from 1974 to 1993 apply to the same sorts of firms from 1993 to 2015, to firms in other broad sectors (such as materials), and to firms in new sectors (such as Internet services). We measure virtually the same scaling exponent for manufacturing for the 1993 to 2015 period as for the 1974 to 1993 period and virtually the same scaling exponent for other sectors as for manufacturing. Furthermore, we show that fluctuations of the growth rate for new industries self-organize into a power law over relatively short time scales., Comment: 15 pages, 7 figures

Although suppressing the spread of a disease is usually achieved by investing in public resources, in the real world only a small percentage of the population have access to government assistance when there is an outbreak, and most must rely on resources from family or friends. We study the dynamics of disease spreading in social-contact multiplex networks when the recovery of infected nodes depends on resources from healthy neighbors in the social layer. We investigate how degree heterogeneity affects the spreading dynamics. Using theoretical analysis and simulations we find that degree heterogeneity promotes disease spreading. The phase transition of the infected density is hybrid and increases smoothly from zero to a finite small value at the first invasion threshold and then suddenly jumps at the second invasion threshold. We also find a hysteresis loop in the transition of the infected density. We further investigate how an overlap in the edges between two layers affects the spreading dynamics. We find that when the amount of overlap is smaller than a critical value the phase transition is hybrid and there is a hysteresis loop, otherwise the phase transition is continuous and the hysteresis loop vanishes. In addition, the edge overlap allows an epidemic outbreak when the transmission rate is below the first invasion threshold, but suppresses any explosive transition when the transmission rate is above the first invasion threshold., Comment: 22 pages, 9 figures

Internet communication channels, e.g., Facebook, Twitter, and email, are multiplex networks that facilitate interaction and information-sharing among individuals. During brief time periods users often use a single communication channel, but then communication channel alteration (CCA) occurs. This means that we must refine our understanding of the dynamics of social contagions. We propose a non-Markovian behavior spreading model in multiplex networks that takes into account the CCA mechanism, and we develop a generalized edge-based compartmental method to describe the spreading dynamics. Through extensive numerical simulations and theoretical analyses we find that the time delays induced by CCA slow the behavior spreading but do not affect the final adoption size. We also find that the CCA suppresses behavior spreading. On two coupled random regular networks, the adoption size exhibits hybrid growth, i.e., it grows first continuously and then discontinuously with the information transmission probability. CCA in ER-SF multiplex networks in which two subnetworks are Erd\H{o}s-R\'{e}nyi (ER) and scale-free (SF) introduces a crossover from continuous to hybrid growth in adoption size versus information transmission probability. Our results extend our understanding of the role of CCA in spreading dynamics, and may elicit further research., Comment: 11 pages, 7 figures, Accepted by Physical Review E

Financial networks have become extremely useful in characterizing the structures of complex financial systems. Meanwhile, the time evolution property of the stock markets can be described by temporal networks. We utilize the temporal network framework to characterize the time-evolving correlation-based networks of stock markets. The market instability can be detected by the evolution of the topology structure of the financial networks. We then employ the temporal centrality as a portfolio selection tool. Those portfolios, which are composed of peripheral stocks with low temporal centrality scores, have consistently better performance under different portfolio optimization frameworks, suggesting that the temporal centrality measure can be used as new portfolio optimization and risk management tool. Our results reveal the importance of the temporal attributes of the stock markets, which should be taken serious consideration in real life applications.

Cascading failures in complex systems have been studied extensively using two different models: $k$-core percolation and interdependent networks. We combine the two models into a general model, solve it analytically and validate our theoretical results through extensive simulations. We also study the complete phase diagram of the percolation transition as we tune the average local $k$-core threshold and the coupling between networks. We find that the phase diagram of the combined processes is very rich and includes novel features that do not appear in the models studying each of the processes separately. For example, the phase diagram consists of first and second-order transition regions separated by two tricritical lines that merge together and enclose a novel two-stage transition region. In the two-stage transition, the size of the giant component undergoes a first-order jump at a certain occupation probability followed by a continuous second-order transition at a lower occupation probability. Furthermore, at certain fixed interdependencies, the percolation transition changes from first-order $\rightarrow$ second-order $\rightarrow$ two-stage $\rightarrow$ first-order as the $k$-core threshold is increased. The analytic equations describing the phase boundaries of the two-stage transition region are set up and the critical exponents for each type of transition are derived analytically.

Being able to predict the occurrence of extreme returns is important in financial risk management. Using the distribution of recurrence intervals---the waiting time between consecutive extremes---we show that these extreme returns are predictable on the short term. Examining a range of different types of returns and thresholds we find that recurrence intervals follow a $q$-exponential distribution, which we then use to theoretically derive the hazard probability $W(\Delta t |t)$. Maximizing the usefulness of extreme forecasts to define an optimized hazard threshold, we indicates a financial extreme occurring within the next day when the hazard probability is greater than the optimized threshold. Both in-sample tests and out-of-sample predictions indicate that these forecasts are more accurate than a benchmark that ignores the predictive signals. This recurrence interval finding deepens our understanding of reoccurring extreme returns and can be applied to forecast extremes in risk management., Comment: 18 pages, 5 figues, 3 tables

Complex systems are composed of mutually interacting components and the output values of these components usually exhibit long-range cross-correlations. Using wavelet analysis, we propose a method of characterizing the joint multifractal nature of these long-range cross correlations, a method we call multifractal cross wavelet analysis (MFXWT). We assess the performance of the MFXWT method by performing extensive numerical experiments on the dual binomial measures with multifractal cross correlations and the bivariate fractional Brownian motions (bFBMs) with monofractal cross correlations. For binomial multifractal measures, we find the empirical joint multifractality of MFXWT to be in approximate agreement with the theoretical formula. For bFBMs, MFXWT may provide spurious multifractality because of the wide spanning range of the multifractal spectrum. We also apply the MFXWT method to stock market indices, and in pairs of index returns and volatilities we find an intriguing joint multifractal behavior. The tests on surrogate series also reveal that the cross correlation behavior, particularly the cross correlation with zero lag, is the main origin of cross multifractality.

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., Comment: 7 figures, typos corrected, references added

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., Comment: 20 pages

Although recent studies have found that the long-term correlations relating to the fat-tailed distribution of inter-event times exist in human activity, and that these correlations indicate the presence of fractality, the property of fractality and its origin have not been analyzed. We use both DFA and MFDFA to analyze the time series in online viewing activity separating from Movielens and Netflix. We find long-term correlations at both the individual and communal level, and that the extent of correlation at the individual level is determined by the activity level. These long-term correlations also indicate that there is fractality in the pattern of online viewing. And, we firstly find a multifractality that results from the combined effect of the fat-tailed distribution of inter-event times (i.e., the times between successive viewing actions of individual) and the long-term correlations in online viewing activity and verify this finding using three synthesized series. Therefore, it can be concluded that the multifractality in online viewing activity is caused by both the fat-tailed distribution of inter-event times and the long-term correlations, and that this enlarges the generic property of human activity to include not just physical space, but also cyberspace., Comment: 10 pages, 13 figures

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., Comment: 16 pages, 7 figures

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.

The robustness of complex networks against node failure and malicious attack has been of interest for decades, while most of the research has focused on random attack or hub-targeted attack. In many real-world scenarios, however, attacks are neither random nor hub-targeted, but localized, where a group of neighboring nodes in a network are attacked and fail. In this paper we develop a percolation framework to analytically and numerically study the robustness of complex networks against such localized attack. In particular, we investigate this robustness in Erd\H{o}s-R\'{e}nyi networks, random-regular networks, and scale-free networks. Our results provide insight into how to better protect networks, enhance cybersecurity, and facilitate the design of more robust infrastructures.

Robustness of network of networks (NON) has been studied only for dependency coupling (Gao et al., 2012) and only for connectivity coupling (Leicht and Souza, arxiv:0907.0894 ). The case of network of networks with both interdependent and interconnected links is more complicated, and also more close to real-life coupled network systems. Understanding the robustness of NON with interdependent and interconnected links is helpful to design resilient infrastructures. Here we develop a framework to study analytically and numerically the robustness of this system with no-feedback and feedback conditions for the case of starlike NON. When assumed that all degree distributions of the connectivity intra- and inter-links are Poissonian, we find that the system undergoes from second order through hybrid order to first order phase transition as coupling strength q increases. Additionally, for both conditions, the results suggest that increasing density of connectivity links (intra-connectivity links or inter-connectivity links) can increase the robustness of the system, while the interdependency links decrease its robustness. Furthermore, by comparing critical attacking strength under same parameters for both conditions, we also find that feedback condition of dependency links, in contrast to no-feedback condition, makes the system extremely vulnerable. Although our detailed analysis is for Poisson degree distribution, the theory can be applied to any degree distribution and other basic topological structure of network such like tree-like and loop-like NON.

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\H{o}s-R\'{e}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., Comment: 31 pages, 12 figures

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 susceptible-infected- susceptible (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 a 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 and 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). We note that, given two isolated networks G1 and G2 with principal eigenvectors x and y, respectively, λ1(A+αB) tends to be higher when nodes i and j with a higher eigenvector component product x iyj are interconnected. This finding suggests essential insights into ways of designing interconnected networks to be robust against epidemics. © 2013 American Physical Society.

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., 4 figures

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 potential, and with predictions for hydrophobic solvation in water using the cavity equation of state and the extended simple point charge 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. This feature of water-like solvation is distinct from the strictly positive and temperature independent enthalpy and entropy of cavity solvation observed in the Lennard-Jones fluid. 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. In addition, the Jagla liquid dewets surfaces of large radii of curvature less readily than the Lennard-Jones liquid, reflecting a greater flexibility or elasticity in the Jagla liquid structure than that of a typical liquid, a behavior also similar to that of water's hydrogen bonding network. 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.

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., Comment: Version 1: 31 pages, 17 figures, 3 tables. Version 2 is streamlined, eliminates substantial material and incorporates referee comments: 19 pages, 14 figures, 3 tables

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., Comment: 36 pages, 20 figures

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., Comment: 11 pages, 2 figures

Studies of liquid water in its supercooled region have helped us better understand the structure and behavior of water. Bulk water freezes at its homogeneous nucleation temperature (approximately 235 K), but protein hydration water avoids this crystallization because each water molecule binds to a protein. Here, we study the dynamics of the hydrogen bond (HB) network of a per-colating layer of water molecules and compare the measurements of a hydrated globular protein with the results of a coarse-grained model that successfully reproduces the properties of hydration water. Using dielectric spectroscopy, we measure the temperature dependence of the relaxation time of proton charge fluctuations. These fluctuations are associated with the dynamics of the HB network of water molecules adsorbed on the protein surface. Using Monte Carlo simulations and mean-field calculations, we study the dynamics and thermodynamics of the model. Both experimental and model analyses are consistent with the interesting possibility of two dynamic crossovers, (i) at approximately 252 K and (ii) at approximately 181 K. Because the experiments agree with the model, we can relate the two crossovers to the presence at ambient pressure of two specific heat maxima. The first is caused by fluctuations in the HB formation, and the second, at a lower temperature, is due to the cooperative reordering of the HB network.

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., Comment: 25 pages in 1 Column Preprint format, 7 Figures, 4 Tables. Version II: changes made in response to referee comments. Note: change in definition of "Paper shares."

We study the volatility time series of 1137 most traded stocks in the US stock markets for the two-year period 2001-02 and analyze their return intervals $\tau$, which are time intervals between volatilities above a given threshold $q$. We explore the probability density function of $\tau$, $P_q(\tau)$, assuming a stretched exponential function, $P_q(\tau) \sim e^{-\tau^\gamma}$. We find that the exponent $\gamma$ depends on the threshold in the range between $q=1$ and 6 standard deviations of the volatility. This finding supports the multiscaling nature of the return interval distribution. To better understand the multiscaling origin, we study how $\gamma$ depends on four essential factors, capitalization, risk, number of trades and return. We show that $\gamma$ depends on the capitalization, risk and return but almost does not depend on the number of trades. This suggests that $\gamma$ relates to the portfolio selection but not on the market activity. To further characterize the multiscaling of individual stocks, we fit the moments of $\tau$, $\mu_m \equiv )^m>^{1/m}$, in the range of $10 < \le 100$ by a power-law, $\mu_m \sim ^\delta$. The exponent $\delta$ is found also to depend on the capitalization, risk and return but not on the number of trades, and its tendency is opposite to that of $\gamma$. Moreover, we show that $\delta$ decreases with $\gamma$ approximately by a linear relation. The return intervals demonstrate the temporal structure of volatilities and our findings suggest that their multiscaling features may be helpful for portfolio optimization., Comment: 16 pages, 6 figures

We investigate the traffic flows of the Korean highway system, which contains both public and private transportation information. We find that the traffic flow T(ij) between city i and j forms a gravity model, the metaphor of physical gravity as described in Newton's law of gravity, P(i)P(j)/r(ij)^2, where P(i) represents the population of city i and r(ij) the distance between cities i and j. It is also shown that the highway network has a heavy tail even though the road network is a rather uniform and homogeneous one. Compared to the highway network, air and public ground transportation establish inhomogeneous systems and have power-law behaviors., Comment: 13 pages

We investigate scaling and memory effects in return intervals between price volatilities above a certain threshold $q$ for the Japanese stock market using daily and intraday data sets. We find that the distribution of return intervals can be approximated by a scaling function that depends only on the ratio between the return interval $\tau$ and its mean $$. We also find memory effects such that a large (or small) return interval follows a large (or small) interval by investigating the conditional distribution and mean return interval. The results are similar to previous studies of other markets and indicate that similar statistical features appear in different financial markets. We also compare our results between the period before and after the big crash at the end of 1989. We find that scaling and memory effects of the return intervals show similar features although the statistical properties of the returns are different., Comment: 11 pages

We propose a numerical method to evaluate the upper critical dimension $d_c$ of random percolation clusters in Erd\H{o}s-R\'{e}nyi networks and in scale-free networks with degree distribution ${\cal P}(k) \sim k^{-\lambda}$, where $k$ is the degree of a node and $\lambda$ is the broadness of the degree distribution. Our results report the theoretical prediction, $d_c = 2(\lambda - 1)/(\lambda - 3)$ for scale-free networks with $3 < \lambda < 4$ and $d_c = 6$ for Erd\H{o}s-R\'{e}nyi networks and scale-free networks with $\lambda > 4$. When the removal of nodes is not random but targeted on removing the highest degree nodes we obtain $d_c = 6$ for all $\lambda > 2$. Our method also yields a better numerical evaluation of the critical percolation threshold, $p_c$, for scale-free networks. Our results suggest that the finite size effects increases when $\lambda$ approaches 3 from above., Comment: 10 pages, 5 figures

We investigate how simultaneously recorded long-range power-law correlated multi-variate signals cross-correlate. To this end we introduce a two-component ARFIMA stochastic process and a two-component FIARCH process to generate coupled fractal signals with long-range power-law correlations which are at the same time long-range cross-correlated. We study how the degree of cross-correlations between these signals depends on the scaling exponents characterizing the fractal correlations in each signal and on the coupling between the signals. Our findings have relevance when studying parallel outputs of multiple-component of physical, physiological and social systems., Comment: 8 pages, 5 figures, elsart.cls

We study a biologically inspired, inherently non-equilibrium model consisting of self-propelled particles. In the model, particles move on a plane with a velocity of constant magnitude; they locally interact with their neighbors by choosing at each time step a velocity direction equal to the average direction of their neighbors. Thus, in the limit of vanishing velocities the model becomes analogous to a Monte-Carlo realization of the classical XY ferromagnet. We show by large-scale numerical simulations that, unlike in the equilibrium XY model, a long-range ordered phase characterized by non-vanishing net flow $\phi$ emerges in this system in a phase space domain bordered by a critical line along which the fluctuations of the order parameter diverge. The corresponding phase diagram as a function of two parameters, the amplitude of noise $\eta$ and the average density of the particles $\varrho$ is calculated and is found to have the form $\eta_c(\varrho)\sim \varrho^{1/2}$. We also find that $\phi$ scales as a function of the external bias $h$ (field or ``wind'') according to a power law $\phi\sim h^{0.9}$. In the ordered phase the system shows long-range correlated fluctuations and $1/f$ noise.

Using molecular dynamics simulations, we study a liquid model which consists of particles interacting via a spherically-symmetric two-scale Jagla ramp potential with both repulsive and attractive ramps. The Jagla potential displays anomalies similar to those found in liquid water, namely expansion upon cooling and an increase of diffusivity upon compression, as well as a liquid-liquid (LL) phase transition in the region of the phase diagram accessible to simulations. The LL coexistence line, unlike in tetrahedrally-coordinated liquids, has a positive slope, because of the Clapeyron relation, corresponding to the fact that the high density phase (HDL) is more ordered than low density phase (LDL). When we cool the system at constant pressure above the critical pressure, the hydrodynamic properties rapidly change from those of LDL-like to those of HDL-like upon crossing the Widom line. The temperature dependence of the diffusivity also changes rapidly in the vicinity of the Widom line, namely the slope of the Arrhenius plot sharply increases upon entering the HDL domain. The properties of the glass transition are different in the two phases, suggesting that the less ordered phase (LDL) is fragile, while the more ordered phase (HDL) is strong, which is consistent with the behavior of tetrahedrally-coordinated liquids such as water silica, silicon, and BeF$_2$., Comment: 9 pages, 18 figures

Coupled nonlinear systems under certain conditions exhibit phase synchronization, which may change for different frequency bands or with presence of additive system noise. In both cases, Fourier filtering is traditionally used to preprocess data. We investigate to what extent the phase synchronization of two coupled R\"{o}ssler oscillators depends on (1) the broadness of their power spectrum, (2) the width of the band-pass filter, and (3) the level of added noise. We find that for identical coupling strengths, oscillators with broader power spectra exhibit weaker synchronization. Further, we find that within a broad band width range, band-pass filtering reduces the effect of noise but can lead to a spurious increase in the degree of synchronization with narrowing band width, even when the coupling between the two oscillators remains the same., Comment: 4 pages,6 figures

We study complex networks with weights w(ij) associated with each link connecting node i and j. The weights are chosen to be correlated with the network topology in the form found in two real world examples: (a) the worldwide airport network and (b) the E. Coli metabolic network. Here w(ij) approximately equals x(ij)(k(i)k(j))alpha, where k(i) and k(j) are the degrees of nodes i and j , x(ij) is a random number, and alpha represents the strength of the correlations. The case alpha0 represents correlation between weights and degree, while alpha0 represents anticorrelation and the case alpha=0 reduces to the case of no correlations. We study the scaling of the lengths of the optimal paths, l(opt), with the system size N in strong disorder for scale-free networks for different alpha. We find two different universality classes for l(opt), in strong disorder depending on alpha: (i) if alpha0 , then for lambda2 the scaling law l(opt) approximately equals N(1/3), where lambda is the power-law exponent of the degree distribution of scale-free networks, and (ii) if alphaor =0 , then l(opt) approximately equals N((nu)(opt)) with nu(opt) identical to its value for the uncorrelated case alpha=0. We calculate the robustness of correlated scale-free networks with different alpha and find the networks with alpha0 to be the most robust networks when compared to the other values of alpha. We propose an analytical method to study percolation phenomena on networks with this kind of correlation, and our numerical results suggest that for scale-free networks with alpha0 , the percolation threshold p(c) is finite for lambda3, which belongs to the same universality class as alpha=0 . We compare our simulation results with the real worldwide airport network, and we find good agreement.

We investigate structural order in glassy water by performing classical molecular dynamics simulations using the extended simple point charge (SPC/E) model of water. We perform isochoric cooling simulations across the glass transition temperature at different cooling rates and densities. We quantify structural order by orientational and translational order metrics. Upon cooling the liquid into the glassy state, both the orientational order parameter $Q$ and translational order parameter $\tau$ increase. At T=0 K, the glasses fall on a line in the $Q$-$\tau$ plane or {\it order map}. The position of this line depends only on density and coincides with the location in the order map of the inherent structures (IS) sampled upon cooling. We evaluate the energy of the IS, $e_{IS}(T)$, and find that both order parameters for the IS are proportional to $e_{IS}$. We also study the structural order during the transformation of low-density amorphous ice (LDA) to high-density amorphous ice (HDA) upon isothermal compression and are able to identify distinct regions in the order map corresponding to these glasses. Comparison of the order parameters for LDA and HDA with those obtained upon isochoric cooling indicates major structural differences between glasses obtained by cooling and glasses obtained by compression. These structural differences are only weakly reflected in the pair correlation function. We also characterize the evolution of structural order upon isobaric annealing, leading at high pressure to very-high density amorphous ice (VHDA)., Comment: submitted

We investigate the relationship between the blood flow velocities (BFV) in the middle cerebral arteries and beat-to-beat blood pressure (BP) recorded from a finger in healthy and post-stroke subjects during the quasi-steady state after perturbation for four different physiologic conditions: supine rest, head-up tilt, hyperventilation and CO_2 rebreathing in upright position. To evaluate whether instantaneous BP changes in the steady state are coupled with instantaneous changes in the BFV, we compare dynamical patterns in the instantaneous phases of these signals, obtained from the Hilbert transform, as a function of time. We find that in post-stroke subjects the instantaneous phase increments of BP and BFV exhibit well pronounced patterns that remain stable in time for all four conditions, while in healthy subjects these patterns are different, less pronounced and more variable. We propose a new approach based on the instantaneous phase increments cross-correlation to quantify the coupling between BP and BFV. We find that the maximum correlation strength is different for the two groups and for the different conditions. For healthy subjects the amplitude of the cross-correlation between the instantaneous phase increments of BP and BFV is small and attenuates within 3-5 heartbeats. In contrast, for post-stroke subjects, this amplitude is significantly larger and cross-correlations persist up to 20 heartbeats. We compare the results of our approach with three complementary methods: direct BP-BFV cross-correlation, transfer function analysis and phase synchronization analysis., Comment: 15 pages, 6 figures; extensively expanded the manuscript including new analyses and results presented in 3 new figures, added more references

We introduce a model of proportional growth to explain the distribution of business firm growth rates. The model predicts that the distribution is exponential in the central part and depicts an asymptotic power-law behavior in the tails with an exponent 3. Because of data limitations, previous studies in this field have been focusing exclusively on the Laplace shape of the body of the distribution. In this article, we test the model at different levels of aggregation in the economy, from products to firms to countries, and we find that the model's predictions agree with empirical growth distributions and size-variance relationships., Comment: 22 pages, 5 Postscript figures, uses revtex4. to be published in Proc. Natl. Acad. Sci. (2005)

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., Comment: Accepted for publication to Physical Review E. 17 Pages, 6 Figures, 1 Table

We study network configurations that provide optimal robustness to random breakdowns for networks with a given number of nodes $N$ and a given cost--which we take as the average number of connections per node $\kav$. We find that the network design that maximizes $f_c$, the fraction of nodes that are randomly removed before global connectivity is lost, consists of $q=[(\kav-1)/\sqrt\kav]\sqrt N$ high degree nodes (``hubs'') of degree $\sqrt{\kav N}$ and $N-q$ nodes of degree 1. Also, we show that $1-f_c$ approaches 0 as $1/\sqrt N$--faster than any other network configuration including scale-free networks. We offer a simple heuristic argument to explain our results.

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., Comment: submitted to Phys. Rev

We respond to the issues discussed by Farmer and Lillo (FL) related to our proposed approach to understanding the origin of power-law distributions in stock price fluctuations. First, we extend our previous analysis to 1000 US stocks and perform a new estimation of market impact that accounts for splitting of large orders and potential autocorrelations in the trade flow. Our new analysis shows clearly that price impact and volume are related by a square-root functional form of market impact for large volumes, in contrast to the claim of FL that this relationship increases as a power law with a smaller exponent. Since large orders are usually executed by splitting into smaller size trades, procedures used by FL give a downward bias for this power law exponent. Second, FL analyze 3 stocks traded on the London Stock Exchange, and solely on this basis they claim that the distribution of transaction volumes do not have a power-law tail for the London Stock Exchange. We perform new empirical analysis on transaction data for the 262 largest stocks listed in the London Stock Exchange, and find that the distribution of volume decays as a power-law with an exponent $\approx 3/2$ -- in sharp contrast to FL's claim that the distribution of transaction volume does not have a power-law tail. Our exponent estimate of $\approx 3/2$ is consistent with our previous results from the New York and Paris Stock Exchanges. We conclude that the available empirical evidence is consistent with our hypothesis on the origin of power-law fluctuations in stock prices., Comment: 6 pages 2 column format

Detrended fluctuation analysis (DFA) and detrended moving average (DMA) are two scaling analysis methods designed to quantify correlations in noisy non-stationary signals. We systematically study the performance of different variants of the DMA method when applied to artificially generated long-range power-law correlated signals with an {\it a-priori} known scaling exponent $\alpha_{0}$ and compare them with the DFA method. We find that the scaling results obtained from different variants of the DMA method strongly depend on the type of the moving average filter. Further, we investigate the optimal scaling regime where the DFA and DMA methods accurately quantify the scaling exponent $\alpha_{0}$, and how this regime depends on the correlations in the signal. Finally, we develop a three-dimensional representation to determine how the stability of the scaling curves obtained from the DFA and DMA methods depends on the scale of analysis, the order of detrending, and the order of the moving average we use, as well as on the type of correlations in the signal., Comment: 15 pages, 16 figures

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^\nu$, where $\nu$ 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 $\xi\sim a^\nu$, below which strong disorder and critical percolation characterize the current path., Comment: 9 pages, 4 figures

We study the optimal distance in networks, $\ell_{\scriptsize opt}$, defined as the length of the path minimizing the total weight, in the presence of disorder. Disorder is introduced by assigning random weights to the links or nodes. For strong disorder, where the maximal weight along the path dominates the sum, we find that $\ell_{\scriptsize opt}\sim N^{1/3}$ in both Erd\H{o}s-R\'enyi (ER) and Watts-Strogatz (WS) networks. For scale free (SF) networks, with degree distribution $P(k) \sim k^{-\lambda}$, we find that $\ell_{\scriptsize opt}$ scales as $N^{(\lambda - 3)/(\lambda - 1)}$ for $3, Comment: 5 pages, 4 figures, accepted for publication in Physical Review Letters

We perform Monte Carlo simulations to determine the average excluded area $$ of randomly oriented squares, randomly oriented widthless sticks and aligned squares in two dimensions. We find significant differences between our results for randomly oriented squares and previous analytical results for the same. The sources of these differences are explained. Using our results for $$ and Monte Carlo simulation results for the percolation threshold, we estimate the mean number of connections per object $B_c$ at the percolation threshold for squares in 2-D. We study systems of squares that are allowed random orientations within a specified angular interval. Our simulations show that the variation in $B_c$ is within 1.6% when the angular interval is varied from 0 to $\pi/2$.