Phase transition in a spatial Lotka-Volterra model

Spatial evolution is investigated in a simulated system of nine competing and mutating bacterium strains, which mimics the biochemical war among bacteria capable of producing two different bacteriocins (toxins) at most. Random sequential dynamics on a square lattice is governed by very symmetrical transition rules for neighborhood invasion of sensitive strains by killers, killers by resistants, and resistants by by sensitives. The community of the nine possible toxicity/resistance types undergoes a critical phase transition as the uniform transmutation rates between the types decreases below a critical value $P_c$ above which all the nine types of strain coexist with equal frequencies. Passing the critical mutation rate from above, the system collapses into one of the three topologically identical states, each consisting of three strain types. Of the three final states each accrues with equal probability and all three maintain themselves in a self-organizing polydomain structure via cyclic invasions. Our Monte Carlo simulations support that this symmetry breaking transition belongs to the universality class of the three-state Potts model.Comment: 4 page

A non trivial extension of the two-dimensional Ising model: the d-dimensional "molecular" model

A recently proposed molecular model is discussed as a non-trivial extension of the Ising model. For d=2 the two models are shown to be equivalent, while for d>2 the molecular model describes a peculiar second order transition from an isotropic high temperature phase to a low-dimensional anisotropic low temperature state. The general mean field analysis is compared with the results achieved by a variational Migdal-Kadanoff real space renormalization group method and by standard Monte Carlo sampling for d=3. By finite size scaling the critical exponent has been found to be 0.44\pm 0.02 thus establishing that the molecular model does not belong to the universality class of the Ising model for d>2.Comment: 25 pages, 5 figure

Defensive alliances in spatial models of cyclical population interactions

As a generalization of the 3-strategy Rock-Scissors-Paper game dynamics in space, cyclical interaction models of six mutating species are studied on a square lattice, in which each species is supposed to have two dominant, two subordinated and a neutral interacting partner. Depending on their interaction topologies, these systems can be classified into four (isomorphic) groups exhibiting significantly different behaviors as a function of mutation rate. On three out of four cases three (or four) species form defensive alliances which maintain themselves in a self-organizing polydomain structure via cyclic invasions. Varying the mutation rate this mechanism results in an ordering phenomenon analogous to that of magnetic Ising model.Comment: 4 pages, 3 figure

Simulation of Potts models with real q and no critical slowing down

A Monte Carlo algorithm is proposed to simulate ferromagnetic q-state Potts model for any real q>0. A single update is a random sequence of disordering and deterministic moves, one for each link of the lattice. A disordering move attributes a random value to the link, regardless of the state of the system, while in a deterministic move this value is a state function. The relative frequency of these moves depends on the two parameters q and beta. The algorithm is not affected by critical slowing down and the dynamical critical exponent z is exactly vanishing. We simulate in this way a 3D Potts model in the range 2<q<3 for estimating the critical value q_c where the thermal transition changes from second-order to first-order, and find q_c=2.620(5).Comment: 5 pages, 3 figures slightly extended version, to appear in Phys. Rev.

Scaling laws for the 2d 8-state Potts model with Fixed Boundary Conditions

We study the effects of frozen boundaries in a Monte Carlo simulation near a first order phase transition. Recent theoretical analysis of the dynamics of first order phase transitions has enabled to state the scaling laws governing the critical regime of the transition. We check these new scaling laws performing a Monte Carlo simulation of the 2d, 8-state spin Potts model. In particular, our results support a pseudo-critical beta finite-size scaling of the form beta(infinity) + a/L + b/L^2, instead of beta(infinity) + c/L^d + d/L^{2d}. Moreover, our value for the latent heat is 0.294(11), which does not coincide with the latent heat analytically derived for the same model if periodic boundary conditions are assumed, which is 0.486358...Comment: 10 pages, 3 postscript figure

Phase transition between synchronous and asynchronous updating algorithms

We update a one-dimensional chain of Ising spins of length $L$ with algorithms which are parameterized by the probability $p$ for a certain site to get updated in one time step. The result of the update event itself is determined by the energy change due to the local change in the configuration. In this way we interpolate between the Metropolis algorithm at zero temperature for $p$ of the order of 1/L and for large $L$, and a synchronous deterministic updating procedure for $p=1$. As function of $p$ we observe a phase transition between the stationary states to which the algorithm drives the system. These are non-absorbing stationary states with antiferromagnetic domains for $p>p_c$, and absorbing states with ferromagnetic domains for $p\leq p_c$. This means that above this transition the stationary states have lost any remnants to the ferromagnetic Ising interaction. A measurement of the critical exponents shows that this transition belongs to the universality class of parity conservation.Comment: 5 pages, 3 figure

Structure Factors and Their Distributions in Driven Two-Species Models

We study spatial correlations and structure factors in a three-state stochastic lattice gas, consisting of holes and two oppositely ``charged'' species of particles, subject to an ``electric'' field at zero total charge. The dynamics consists of two nearest-neighbor exchange processes, occuring on different times scales, namely, particle-hole and particle-particle exchanges. Using both, Langevin equations and Monte Carlo simulations, we study the steady-state structure factors and correlation functions in the disordered phase, where density profiles are homogeneous. In contrast to equilibrium systems, the average structure factors here show a discontinuity singularity at the origin. The associated spatial correlation functions exhibit intricate crossovers between exponential decays and power laws of different kinds. The full probability distributions of the structure factors are universal asymmetric exponential distributions.Comment: RevTex, 18 pages, 4 postscript figures included, mistaken half-empty page correcte

Quantifying the levitation picture of extended states in lattice models

The behavior of extended states is quantitatively analyzed for two dimensional lattice models. A levitation picture is established for both white-noise and correlated disorder potentials. In a continuum limit window of the lattice models we find simple quantitative expressions for the extended states levitation, suggesting an underlying universal behavior. On the other hand, these results point out that the Quantum Hall phase diagrams may be disorder dependent.Comment: 5 pages, submitted to PR

On the Potts model partition function in an external field

We study the partition function of Potts model in an external (magnetic) field, and its connections with the zero-field Potts model partition function. Using a deletion-contraction formulation for the partition function Z for this model, we show that it can be expanded in terms of the zero-field partition function. We also show that Z can be written as a sum over the spanning trees, and the spanning forests, of a graph G. Our results extend to Z the well-known spanning tree expansion for the zero-field partition function that arises though its connections with the Tutte polynomial

Solar Magnetic Carpet I: Simulation of Synthetic Magnetograms

This paper describes a new 2D model for the photospheric evolution of the magnetic carpet. It is the first in a series of papers working towards constructing a realistic 3D non-potential model for the interaction of small-scale solar magnetic fields. In the model, the basic evolution of the magnetic elements is governed by a supergranular flow profile. In addition, magnetic elements may evolve through the processes of emergence, cancellation, coalescence and fragmentation. Model parameters for the emergence of bipoles are based upon the results of observational studies. Using this model, several simulations are considered, where the range of flux with which bipoles may emerge is varied. In all cases the model quickly reaches a steady state where the rates of emergence and cancellation balance. Analysis of the resulting magnetic field shows that we reproduce observed quantities such as the flux distribution, mean field, cancellation rates, photospheric recycle time and a magnetic network. As expected, the simulation matches observations more closely when a larger, and consequently more realistic, range of emerging flux values is allowed (4e16 - 1e19 Mx). The model best reproduces the current observed properties of the magnetic carpet when we take the minimum absolute flux for emerging bipoles to be 4e16 Mx. In future, this 2D model will be used as an evolving photospheric boundary condition for 3D non-potential modeling.Comment: 33 pages, 16 figures, 5 gif movies included: movies may be viewed at http://www-solar.mcs.st-and.ac.uk/~karen/movies_paper1