SIRS dynamics on random networks: simulations and analytical models

The standard pair approximation equations (PA) for the Susceptible-Infective-Recovered-Susceptible (SIRS) model of infection spread on a network of homogeneous degree $k$ predict a thin phase of sustained oscillations for parameter values that correspond to diseases that confer long lasting immunity. Here we present a study of the dependence of this oscillatory phase on the parameter $k$ and of its relevance to understand the behaviour of simulations on networks. For $k=4$, we compare the phase diagram of the PA model with the results of simulations on regular random graphs (RRG) of the same degree. We show that for parameter values in the oscillatory phase, and even for large system sizes, the simulations either die out or exhibit damped oscillations, depending on the initial conditions. This failure of the standard PA model to capture the qualitative behaviour of the simulations on large RRGs is currently being investigated.Comment: 6 pages, 3 figures, WIPP to be published in Conference proceedings Complex'2009 February 23-25, Shanghai, Chin

Aging without disorder on long time scales

We study the Metropolis dynamics of a simple spin system without disorder, which exhibits glassy dynamics at low temperatures. We use an implementation of the algorithm of Bortz, Kalos and Lebowitz \cite{bortz}. This method turns out to be very efficient for the study of glassy systems, which get trapped in local minima on many different time scales. We find strong evidence of aging effects at low temperatures. We relate these effects to the distribution function of the trapping times of single configurations.Comment: 8 pages Revtex, 7 figures uuencoded (Revised version: the figures are now present

Rejection-free Monte Carlo Algorithms for Models with Continuous Degrees of Freedom

We construct a rejection-free Monte Carlo algorithm for a system with continuous degrees of freedom. We illustrate the algorithm by applying it to the classical three-dimensional Heisenberg model with canonical Metropolis dynamics. We obtain the lifetime of the metastable state following a reversal of the external magnetic field. Our rejection-free algorithm obtains results in agreement with a direct implementation of the Metropolis dynamic and requires orders of magnitude less computational time at low temperatures. The treatment is general and can be extended to other dynamics and other systems with continuous degrees of freedom.Comment: 4 pages, including figures. PRE, in pres

Glassiness and constrained dynamics of a short-range non-disordered spin model

We study the low temperature dynamics of a two dimensional short-range spin system with uniform ferromagnetic interactions, which displays glassiness at low temperatures despite the absence of disorder or frustration. The model has a dual description in terms of free defects subject to dynamical constraints, and is an explicit realization of the ``hierarchically constrained dynamics'' scenario for glassy systems. We give a number of exact results for the statics of the model, and study in detail the dynamical behaviour of one-time and two-time quantities. We also consider the role played by the configurational entropy, which can be computed exactly, in the relation between fluctuations and response.Comment: 10 pages, 9 figures; minor changes, references adde

Effect of the vacancy interaction on antiphase domain-growth in a two-dimensional binary alloy

The influence of diffusing vacancies on the antiphase domain growth process in a binary alloy is studied by Monte Carlo simulations. The system is modelled by means of a Blume-Emery-Griffiths hamiltonian with a biquadratic coupling parameter K controlling the microscopic interactions between vacancies. We obtain that, independently of K, the vacancies exhibit a tendency to concentrate on the antiphase boundaries. This gives rise to an effective interactions between movin interfaces and diffusing vacancies which strongly influences the domain growth process. One distinguishes three different behaviours: i) for K 1 the growth is slown down but still curvature driven