Nonextensivity of the cyclic Lattice Lotka Volterra model

We numerically show that the Lattice Lotka-Volterra model, when realized on a square lattice support, gives rise to a {\it finite} production, per unit time, of the nonextensive entropy $S_q= \frac{1- \sum_ip_i^q}{q-1}$ $(S_1=-\sum_i p_i \ln p_i)$. This finiteness only occurs for $q=0.5$ for the $d=2$ growth mode (growing droplet), and for $q=0$ for the $d=1$ one (growing stripe). This strong evidence of nonextensivity is consistent with the spontaneous emergence of local domains of identical particles with fractal boundaries and competing interactions. Such direct evidence is for the first time exhibited for a many-body system which, at the mean field level, is conservative.Comment: Latex, 6 pages, 5 figure

Co-existence in the two-dimensional May-Leonard model with random rates

We employ Monte Carlo simulations to numerically study the temporal evolution and transient oscillations of the population densities, the associated frequency power spectra, and the spatial correlation functions in the (quasi-)steady state in two-dimensional stochastic May--Leonard models of mobile individuals, allowing for particle exchanges with nearest-neighbors and hopping onto empty sites. We therefore consider a class of four-state three-species cyclic predator-prey models whose total particle number is not conserved. We demonstrate that quenched disorder in either the reaction or in the mobility rates hardly impacts the dynamical evolution, the emergence and structure of spiral patterns, or the mean extinction time in this system. We also show that direct particle pair exchange processes promote the formation of regular spiral structures. Moreover, upon increasing the rates of mobility, we observe a remarkable change in the extinction properties in the May--Leonard system (for small system sizes): (1) As the mobility rate exceeds a threshold that separates a species coexistence (quasi-)steady state from an absorbing state, the mean extinction time as function of system size N crosses over from a functional form ~ e^{cN} / N (where c is a constant) to a linear dependence; (2) the measured histogram of extinction times displays a corresponding crossover from an (approximately) exponential to a Gaussian distribution. The latter results are found to hold true also when the mobility rates are randomly distributed.Comment: 9 pages, 4 figures; to appear in Eur. Phys. J. B (2011