Network Transitivity and Matrix Models

This paper is a step towards a systematic theory of the transitivity (clustering) phenomenon in random networks. A static framework is used, with adjacency matrix playing the role of the dynamical variable. Hence, our model is a matrix model, where matrices are random, but their elements take values 0 and 1 only. Confusion present in some papers where earlier attempts to incorporate transitivity in a similar framework have been made is hopefully dissipated. Inspired by more conventional matrix models, new analytic techniques to develop a static model with non-trivial clustering are introduced. Computer simulations complete the analytic discussion.Comment: 11 pages, 7 eps figures, 2-column revtex format, print bug correcte

Quantum Monte Carlo study of a magnetic-field-driven 2D superconductor-insulator transition

We numerically study the superconductor-insulator phase transition in a model disordered 2D superconductor as a function of applied magnetic field. The calculation involves quantum Monte Carlo calculations of the (2+1)D XY model in the presence of both disorder and magnetic field. The XY coupling is assumed to have the form -J\cos(\theta_i-\theta_j-A_{ij}), where A_{ij} has a mean of zero and a standard deviation \Delta A_{ij}. In a real system, such a model would be approximately realized by a 2D array of small Josephson-coupled grains with slight spatial disorder and a uniform applied magnetic field. The different values \Delta A_{ij} then corresponds to an applied field such that the average number of flux quanta per plaquette has various integer values N: larger N corresponds to larger \Delta A_{ij}. For any value of \Delta A_{ij}, there appears to be a critical coupling constant K_c(\Delta A_{ij})=\sqrt{[J/(2U)]_c}, where U is the charging energy, above which the system is a Mott insulator; there is also a corresponding critical conductivity \sigma^*(\Delta A_{ij}) at the transition. For \Delta A_{ij}=\infty, the order parameter of the transition is a renormalized coupling constant g. Using a numerical technique appropriate for disordered systems, we show that the transition at this value of \Delta A_{ij} takes place from an insulating (I) phase to a Bose glass (BG) phase, and that the dynamical critical exponent characterizing this transition is z \sim 1.3. By contrast, z=1 for this model at \Delta A_{ij}=0. We suggest that the superconductor to insulator transition is actually of this I to BG class at all nonzero \Delta A_{ij}'s, and we support this interpretation by both numerical evidence and an analytical argument based on the Harris criterion.Comment: 17 pages, 23 figures, accepted for publication in Phys. Rev.

Topological enslavement in evolutionary games on correlated multiplex networks

Governments and enterprises strongly rely on incentives to generate favorable outcomes from social and strategic interactions between individuals. The incentives are usually modeled by payoffs in evolutionary games, such as the prisoner's dilemma or the harmony game, with imitation dynamics. Adjusting the incentives by changing the payoff parameters can favor cooperation, as found in the harmony game, over defection, which prevails in the prisoner's dilemma. Here, we show that this is not always the case if individuals engage in strategic interactions in multiple domains. In particular, we investigate evolutionary games on multiplex networks where individuals obtain an aggregate payoff. We explicitly control the strength of degree correlations between nodes in the different layers of the multiplex. We find that if the multiplex is composed of many layers and degree correlations are strong, the topology of the system enslaves the dynamics and the final outcome, cooperation or defection, becomes independent of the payoff parameters. The fate of the system is then determined by the initial conditions

Video Pandemics: Worldwide Viral Spreading of Psy's Gangnam Style Video

Viral videos can reach global penetration traveling through international channels of communication similarly to real diseases starting from a well-localized source. In past centuries, disease fronts propagated in a concentric spatial fashion from the the source of the outbreak via the short range human contact network. The emergence of long-distance air-travel changed these ancient patterns. However, recently, Brockmann and Helbing have shown that concentric propagation waves can be reinstated if propagation time and distance is measured in the flight-time and travel volume weighted underlying air-travel network. Here, we adopt this method for the analysis of viral meme propagation in Twitter messages, and define a similar weighted network distance in the communication network connecting countries and states of the World. We recover a wave-like behavior on average and assess the randomizing effect of non-locality of spreading. We show that similar result can be recovered from Google Trends data as well.Comment: 10 page

Nonlinear Dynamics of Capacitive Charging and Desalination by Porous Electrodes

The rapid and efficient exchange of ions between porous electrodes and aqueous solutions is important in many applications, such as electrical energy storage by super-capacitors, water desalination and purification by capacitive deionization (or desalination), and capacitive extraction of renewable energy from a salinity difference. Here, we present a unified mean-field theory for capacitive charging and desalination by ideally polarizable porous electrodes (without Faradaic reactions or specific adsorption of ions) in the limit of thin double layers (compared to typical pore dimensions). We illustrate the theory in the case of a dilute, symmetric, binary electrolyte using the Gouy-Chapman-Stern (GCS) model of the double layer, for which simple formulae are available for salt adsorption and capacitive charging of the diffuse part of the double layer. We solve the full GCS mean-field theory numerically for realistic parameters in capacitive deionization, and we derive reduced models for two limiting regimes with different time scales: (i) In the "super-capacitor regime" of small voltages and/or early times where the porous electrode acts like a transmission line, governed by a linear diffusion equation for the electrostatic potential, scaled to the RC time of a single pore. (ii) In the "desalination regime" of large voltages and long times, the porous electrode slowly adsorbs neutral salt, governed by coupled, nonlinear diffusion equations for the pore-averaged potential and salt concentration

Scaling Properties of Random Walks on Small-World Networks

Using both numerical simulations and scaling arguments, we study the behavior of a random walker on a one-dimensional small-world network. For the properties we study, we find that the random walk obeys a characteristic scaling form. These properties include the average number of distinct sites visited by the random walker, the mean-square displacement of the walker, and the distribution of first-return times. The scaling form has three characteristic time regimes. At short times, the walker does not see the small-world shortcuts and effectively probes an ordinary Euclidean network in $d$-dimensions. At intermediate times, the properties of the walker shows scaling behavior characteristic of an infinite small-world network. Finally, at long times, the finite size of the network becomes important, and many of the properties of the walker saturate. We propose general analytical forms for the scaling properties in all three regimes, and show that these analytical forms are consistent with our numerical simulations.Comment: 7 pages, 8 figures, two-column format. Submitted to PR

Sign-time distributions for interface growth

We apply the recently introduced distribution of sign-times (DST) to non-equilibrium interface growth dynamics. We are able to treat within a unified picture the persistence properties of a large class of relaxational and noisy linear growth processes, and prove the existence of a non-trivial scaling relation. A new critical dimension is found, relating to the persistence properties of these systems. We also illustrate, by means of numerical simulations, the different types of DST to be expected in both linear and non-linear growth mechanisms.Comment: 4 pages, 5 ps figs, replaced misprint in authors nam

Hawking Radiation of Dirac Particles in a Variable-mass Kerr Space-time

Hawking effect of Dirac particles in a variable-mass Kerr space-time is investigated by using a method called as the generalized tortoise coordinate transformation. The location and the temperature of the event horizon of the non-stationary Kerr black hole are derived. It is shown that the temperature and the shape of the event horizon depend not only on the time but also on the angle. However, the Fermi-Dirac spectrum displays a residual term which is absent from that of Bose-Einstein distribution.Comment: 12 pages in 12pt Revtex, no figure, to appear in Gen. Rel. Grav. Vol.33, No.7 (2001

Properties of Random Graphs with Hidden Color

We investigate in some detail a recently suggested general class of ensembles of sparse undirected random graphs based on a hidden stub-coloring, with or without the restriction to nondegenerate graphs. The calculability of local and global structural properties of graphs from the resulting ensembles is demonstrated. Cluster size statistics are derived with generating function techniques, yielding a well-defined percolation threshold. Explicit rules are derived for the enumeration of small subgraphs. Duality and redundancy is discussed, and subclasses corresponding to commonly studied models are identified.Comment: 14 pages, LaTeX, no figure

Exact eigenvalue spectrum of a class of fractal scale-free networks

The eigenvalue spectrum of the transition matrix of a network encodes important information about its structural and dynamical properties. We study the transition matrix of a family of fractal scale-free networks and analytically determine all the eigenvalues and their degeneracies. We then use these eigenvalues to evaluate the closed-form solution to the eigentime for random walks on the networks under consideration. Through the connection between the spectrum of transition matrix and the number of spanning trees, we corroborate the obtained eigenvalues and their multiplicities.Comment: Definitive version accepted for publication in EPL (Europhysics Letters