On Damage Spreading Transitions

We study the damage spreading transition in a generic one-dimensional stochastic cellular automata with two inputs (Domany-Kinzel model) Using an original formalism for the description of the microscopic dynamics of the model, we are able to show analitically that the evolution of the damage between two systems driven by the same noise has the same structure of a directed percolation problem. By means of a mean field approximation, we map the density phase transition into the damage phase transition, obtaining a reliable phase diagram. We extend this analysis to all symmetric cellular automata with two inputs, including the Ising model with heath-bath dynamics.Comment: 12 pages LaTeX, 2 PostScript figures, tar+gzip+u

Synchronization and Maximum Lyapunov Exponent in Cellular Automata

We study the synchronization of totalistic one dimensional cellular automata (CA). The CA with a non zero synchronization threshold exhibit complex non periodic space time patterns and conversely. This synchronization transition is related to directed percolation. We study also the maximum Lyapunov exponent for CA, defined in analogy with continuous dynamical systems as the exponential rate of expansion of the linear map induced by the evolution rule of CA, constructed with the aid of the Boolean derivatives. The synchronization threshold is strongly correlated to the maximum Lyapunov exponent and we propose approximate relations between these quantities. The value of this threshold can be used to parametrize the space time complexity of CA.Comment: 8 pages and 2 figure

A Self-Organized Method for Computing the Epidemic Threshold in Computer Networks

In many cases, tainted information in a computer network can spread in a way similar to an epidemics in the human world. On the other had, information processing paths are often redundant, so a single infection occurrence can be easily "reabsorbed". Randomly checking the information with a central server is equivalent to lowering the infection probability but with a certain cost (for instance processing time), so it is important to quickly evaluate the epidemic threshold for each node. We present a method for getting such information without resorting to repeated simulations. As for human epidemics, the local information about the infection level (risk perception) can be an important factor, and we show that our method can be applied to this case, too. Finally, when the process to be monitored is more complex and includes "disruptive interference", one has to use actual simulations, which however can be carried out "in parallel" for many possible infection probabilities

Small world effects in evolution

For asexual organisms point mutations correspond to local displacements in the genotypic space, while other genotypic rearrangements represent long-range jumps. We investigate the spreading properties of an initially homogeneous population in a flat fitness landscape, and the equilibrium properties on a smooth fitness landscape. We show that a small-world effect is present: even a small fraction of quenched long-range jumps makes the results indistinguishable from those obtained by assuming all mutations equiprobable. Moreover, we find that the equilibrium distribution is a Boltzmann one, in which the fitness plays the role of an energy, and mutations that of a temperature.Comment: 13 pages and 5 figures. New revised versio

A population study of type II bursts in the Rapid Burster

Type II bursts are thought to arise from instabilities in the accretion flow onto a neutron star in an X-ray binary. Despite having been known for almost 40 years, no model can yet satisfactorily account for all their properties. To shed light on the nature of this phenomenon and provide a reference for future theoretical work, we study the entire sample of Rossi X-ray Timing Explorer data of type II bursts from the Rapid Burster (MXB 1730-335). We find that type II bursts are Eddington-limited in flux, that a larger amount of energy goes in the bursts than in the persistent emission, that type II bursts can be as short as 0.130 s, and that the distribution of recurrence times drops abruptly below 15-18 s. We highlight the complicated feedback between type II bursts and the NS surface thermonuclear explosions known as type I bursts, and between type II bursts and the persistent emission. We review a number of models for type II bursts. While no model can reproduce all the observed burst properties and explain the source uniqueness, models involving a gating role for the magnetic field come closest to matching the properties of our sample. The uniqueness of the source may be explained by a special combination of magnetic field strength, stellar spin period and alignment between the magnetic field and the spin axis.Comment: Accepted 2015 February 12. Received 2015 February 10; in original form 2014 December 1

Nature of phase transitions in a probabilistic cellular automaton with two absorbing states

We present a probabilistic cellular automaton with two absorbing states, which can be considered a natural extension of the Domany-Kinzel model. Despite its simplicity, it shows a very rich phase diagram, with two second-order and one first-order transition lines that meet at a tricritical point. We study the phase transitions and the critical behavior of the model using mean field approximations, direct numerical simulations and field theory. A closed form for the dynamics of the kinks between the two absorbing phases near the tricritical point is obtained, providing an exact correspondence between the presence of conserved quantities and the symmetry of absorbing states. The second-order critical curves and the kink critical dynamics are found to be in the directed percolation and parity conservation universality classes, respectively. The first order phase transition is put in evidence by examining the hysteresis cycle. We also study the "chaotic" phase, in which two replicas evolving with the same noise diverge, using mean field and numerical techniques. Finally, we show how the shape of the potential of the field-theoretic formulation of the problem can be obtained by direct numerical simulations.Comment: 19 pages with 7 figure

Fast vectorized algorithm for the Monte Carlo Simulation of the Random Field Ising Model

An algoritm for the simulation of the 3--dimensional random field Ising model with a binary distribution of the random fields is presented. It uses multi-spin coding and simulates 64 physically different systems simultaneously. On one processor of a Cray YMP it reaches a speed of 184 Million spin updates per second. For smaller field strength we present a version of the algorithm that can perform 242 Million spin updates per second on the same machine.Comment: 13 pp., HLRZ 53/9

Quasi-stationary distributions for the Domany-Kinzel stochastic cellular automaton

We construct the {\it quasi-stationary} (QS) probability distribution for the Domany-Kinzel stochastic cellular automaton (DKCA), a discrete-time Markov process with an absorbing state. QS distributions are derived at both the one- and two-site levels. We characterize the distribuitions by their mean, and various moment ratios, and analyze the lifetime of the QS state, and the relaxation time to attain this state. Of particular interest are the scaling properties of the QS state along the critical line separating the active and absorbing phases. These exhibit a high degree of similarity to the contact process and the Malthus-Verhulst process (the closest continuous-time analogs of the DKCA), which extends to the scaling form of the QS distribution.Comment: 15 pages, 9 figures, submited to PR