New approach to Dynamical Monte Carlo Methods: application to an Epidemic Model

A new approach to Dynamical Monte Carlo Methods is introduced to simulate markovian processes. We apply this approach to formulate and study an epidemic Generalized SIRS model. The results are in excellent agreement with the forth order Runge-Kutta Method in a region of deterministic solution. We also demonstrate that purely local interactions reproduce a poissonian-like process at mesoscopic level. The simulations for this case are checked self-consistently using a stochastic version of the Euler Method.Comment: Written with Scientific WorkPlace 3.51 in REVTex4 format, 11 pages with 2 figures in postscript forma

First-order transitions and triple point on a random p-spin interaction model

The effects of competing quadrupolar- and spin-glass orderings are investigated on a spin-1 Ising model with infinite-range random $p$-spin interactions. The model is studied through the replica approach and a phase diagram is obtained in the limit $p\to\infty$. The phase diagram, obtained within replica-symmetry breaking, exhibits a very unusual feature in magnetic models: three first-order transition lines meeting at a commom triple point, where all phases of the model coexist.Comment: 9 pages, 2 ps figures include

Fermionic van Hemmen Spin Glass Model with a Transverse Field

In the present work it is studied the fermionic van Hemmen model for the spin glass (SG) with a transverse magnetic field $\Gamma$. In this model, the spin operators are written as a bilinear combination of fermionic operators, which allows the analysis of the interplay between charge and spin fluctuations in the presence of a quantum spin flipping mechanism given by $\Gamma$. The problem is expressed in the fermionic path integral formalism. As results, magnetic phase diagrams of temperature versus the ferromagnetic interaction are obtained for several values of chemical potential $\mu$ and $\Gamma$. The $\Gamma$ field suppresses the magnetic orders. The increase of $\mu$ alters the average occupation per site that affects the magnetic phases. For instance, the SG and the mixed SG+ferromagnetic phases are also suppressed by $\mu$. In addition, $\mu$ can change the nature of the phase boundaries introducing a first order transition.Comment: 9 pages, 4 figures, accepted for publication in Phys. Lett.

Sandpiles with height restrictions

We study stochastic sandpile models with a height restriction in one and two dimensions. A site can topple if it has a height of two, as in Manna's model, but, in contrast to previously studied sandpiles, here the height (or number of particles per site), cannot exceed two. This yields a considerable simplification over the unrestricted case, in which the number of states per site is unbounded. Two toppling rules are considered: in one, the particles are redistributed independently, while the other involves some cooperativity. We study the fixed-energy system (no input or loss of particles) using cluster approximations and extensive simulations, and find that it exhibits a continuous phase transition to an absorbing state at a critical value zeta_c of the particle density. The critical exponents agree with those of the unrestricted Manna sandpile.Comment: 10 pages, 14 figure

Selforganized 3-band structure of the doped fermionic Ising spin glass

The fermionic Ising spin glass is analyzed for arbitrary filling and for all temperatures. A selforganized 3-band structure of the model is obtained in the magnetically ordered phase. Deviation from half filling generates a central nonmagnetic band, which becomes sharply separated at T=0 by (pseudo)gaps from upper and lower magnetic bands. Replica symmetry breaking effects are derived for several observables and correlations. They determine the shape of the 3-band DoS, and, for given chemical potential, influence the fermion filling strongly in the low temperature regime.Comment: 13 page

A p-Spin Interaction Ashkin-Teller Spin-Glass Model

A p-spin interaction Ashkin-Teller spin glass, with three independent Gaussian probability distributions for the exchange interactions, is studied by means of the replica method. A simple phase diagram is obtained within the replica-symmetric approximation, presenting an instability of the paramagnetic solution at low temperatures. The replica-symmetry-breaking procedure is implemented and a rich phase diagram is obtained; besides the paramagnetic phase, three distinct spin-glass phases appear. Three first-order critical frontiers are found and they all meet at a triple point; among such lines, two of them present discontinuities in the order parameters, but no latent heat, whereas the other one exhibits both discontinuities in the order parameters and a finite latent heat.Comment: 17 pages, 2 figures, submitted to Physica

Curie Temperatures for Three-Dimensional Binary Ising Ferromagnets

Using the Swendsen and Wang algorithm, high accuracy Monte Carlo simulations were performed to study the concentration dependence of the Curie temperature in binary, ferromagnetic Ising systems on the simple-cubic lattice. Our results are in good agreement with known mean-field like approaches. Based on former theoretical formulas we propose a new way of estimating the Curie temperature of these systems.Comment: nr. of pages:13, LATEX. Version 2.09, Scientific Report :02/1994 (Univ. of Bergen, Norway), 7 figures upon reques

Phase transition of meshwork models for spherical membranes

We have studied two types of meshwork models by using the canonical Monte Carlo simulation technique. The first meshwork model has elastic junctions, which are composed of vertices, bonds, and triangles, while the second model has rigid junctions, which are hexagonal (or pentagonal) rigid plates. Two-dimensional elasticity is assumed only at the elastic junctions in the first model, and no two-dimensional bending elasticity is assumed in the second model. Both of the meshworks are of spherical topology. We find that both models undergo a first-order collapsing transition between the smooth spherical phase and the collapsed phase. The Hausdorff dimension of the smooth phase is H\simeq 2 in both models as expected. It is also found that H\simeq 2 in the collapsed phase of the second model, and that H is relatively larger than 2 in the collapsed phase of the first model, but it remains in the physical bound, i.e., H<3. Moreover, the first model undergoes a discontinuous surface fluctuation transition at the same transition point as that of the collapsing transition, while the second model undergoes a continuous transition of surface fluctuation. This indicates that the phase structure of the meshwork model is weakly dependent on the elasticity at the junctions.Comment: 21 pages, 12 figure

Generating a checking sequence with a minimum number of reset transitions

Given a finite state machine M, a checking sequence is an input sequence that is guaranteed to lead to a failure if the implementation under test is faulty and has no more states than M. There has been much interest in the automated generation of a short checking sequence from a finite state machine. However, such sequences can contain reset transitions whose use can adversely affect both the cost of applying the checking sequence and the effectiveness of the checking sequence. Thus, we sometimes want a checking sequence with a minimum number of reset transitions rather than a shortest checking sequence. This paper describes a new algorithm for generating a checking sequence, based on a distinguishing sequence, that minimises the number of reset transitions used.This work was supported in part by Leverhulme Trust grant number F/00275/D, Testing State Based Systems, Natural Sciences and Engineering Research Council (NSERC) of Canada grant number RGPIN 976, and Engineering and Physical Sciences Research Council grant number GR/R43150, Formal Methods and Testing (FORTEST)

Transport Properties of Random Walks on Scale-Free/Regular-Lattice Hybrid Networks

We study numerically the mean access times for random walks on hybrid disordered structures formed by embedding scale-free networks into regular lattices, considering different transition rates for steps across lattice bonds ($F$) and across network shortcuts ($f$). For fast shortcuts ($f/F\gg 1$) and low shortcut densities, traversal time data collapse onto an universal curve, while a crossover behavior that can be related to the percolation threshold of the scale-free network component is identified at higher shortcut densities, in analogy to similar observations reported recently in Newman-Watts small-world networks. Furthermore, we observe that random walk traversal times are larger for networks with a higher degree of inhomogeneity in their shortcut distribution, and we discuss access time distributions as functions of the initial and final node degrees. These findings are relevant, in particular, when considering the optimization of existing information networks by the addition of a small number of fast shortcut connections.Comment: 8 pages, 6 figures; expanded discussions, added figures and references. To appear in J Stat Phy