Nonlinear q-voter model with inflexible zealots

We study the dynamics of the nonlinear q-voter model with inflexible zealots in a finite well-mixed population. In this system, each individual supports of one of two parties and is either a susceptible voter or an inflexible zealot. At each time step, a susceptible adopts the opinion of a neighbor if this belongs to a group of q ≥ 2 neighbors all in the same state, whereas inflexible zealots never change their opinion. In the presence of zealots of both parties the model is characterized by a fluctuating stationary state and, below a zealotry density threshold, the distribution of opinions is bimodal. After a characteristic time, most susceptibles become supporters of the party having more zealots and the opinion distribution is asymmetric. When the number of zealots of both parties is the same, the opinion distribution is symmetric and, in the long run, susceptibles endlessly swing from the state where they all support one party to the opposite state. Above the zealotry density threshold, when there is an unequal number of zealots of each type, the probability distribution is single-peaked and non-Gaussian. These properties are investigated analytically and with stochastic simulations. We also study the mean time to reach a consensus when zealots support only one party

Nonuniform autonomous one-dimensional exclusion nearest-neighbor reaction-diffusion models

The most general nonuniform reaction-diffusion models on a one-dimensional lattice with boundaries, for which the time evolution equations of corre- lation functions are closed, are considered. A transfer matrix method is used to find the static solution. It is seen that this transfer matrix can be obtained in a closed form, if the reaction rates satisfy certain conditions. We call such models superautonomous. Possible static phase transitions of such models are investigated. At the end, as an example of superau- tonomous models, a nonuniform voter model is introduced, and solved explicitly.Comment: 14 page

Influence of local carrying capacity restrictions on stochastic predator-prey models

We study a stochastic lattice predator-prey system by means of Monte Carlo simulations that do not impose any restrictions on the number of particles per site, and discuss the similarities and differences of our results with those obtained for site-restricted model variants. In accord with the classic Lotka-Volterra mean-field description, both species always coexist in two dimensions. Yet competing activity fronts generate complex, correlated spatio-temporal structures. As a consequence, finite systems display transient erratic population oscillations with characteristic frequencies that are renormalized by fluctuations. For large reaction rates, when the processes are rendered more local, these oscillations are suppressed. In contrast with site-restricted predator-prey model, we observe species coexistence also in one dimension. In addition, we report results on the steady-state prey age distribution.Comment: Latex, IOP style, 17 pages, 9 figures included, related movies available at http://www.phys.vt.edu/~tauber/PredatorPrey/movies

Spatial Patterns Emerging from a Stochastic Process Near Criticality

There is mounting empirical evidence that many communities of living organisms display key features which closely resemble those of physical systems at criticality. We here introduce a minimal model framework for the dynamics of a community of individuals which undergoes local birth-death, immigration, and local jumps on a regular lattice. We study its properties when the system is close to its critical point. Even if this model violates detailed balance, within a physically relevant regime dominated by fluctuations, it is possible to calculate analytically the probability density function of the number of individuals living in a given volume, which captures the close-to-critical behavior of the community across spatial scales. We find that the resulting distribution satisfies an equation where spatial effects are encoded in appropriate functions of space, which we calculate explicitly. The validity of the analytical formulae is confirmed by simulations in the expected regimes. We finally discuss how this model in the critical-like regime is in agreement with several biodiversity patterns observed in tropical rain forests

Survival behavior in the cyclic Lotka-Volterra model with a randomly switching reaction rate

We study the influence of a randomly switching reproduction-predation rate on the survival behavior of the nonspatial cyclic Lotka-Volterra model, also known as the zero-sum rock-paper-scissors game, used to metaphorically describe the cyclic competition between three species. In large and finite populations, demographic fluctuations (internal noise) drive two species to extinction in a finite time, while the species with the smallest reproduction-predation rate is the most likely to be the surviving one (law of the weakest). Here we model environmental (external) noise by assuming that the reproduction-predation rate of the strongest species (the fastest to reproduce and predate) in a given static environment randomly switches between two values corresponding to more and less favorable external conditions. We study the joint effect of environmental and demographic noise on the species survival probabilities and on the mean extinction time. In particular, we investigate whether the survival probabilities follow the law of the weakest and analyze their dependence on the external noise intensity and switching rate. Remarkably, when, on average, there is a finite number of switches prior to extinction, the survival probability of the predator of the species whose reaction rate switches typically varies nonmonotonically with the external noise intensity (with optimal survival about a critical noise strength). We also outline the relationship with the case where all reaction rates switch on markedly different time scales

Solution of a one-dimensional stochastic model with branching and coagulation reactions

We solve an one-dimensional stochastic model of interacting particles on a chain. Particles can have branching and coagulation reactions, they can also appear on an empty site and disappear spontaneously. This model which can be viewed as an epidemic model and/or as a generalization of the {\it voter} model, is treated analytically beyond the {\it conventional} solvable situations. With help of a suitably chosen {\it string function}, which is simply related to the density and the non-instantaneous two-point correlation functions of the particles, exact expressions of the density and of the non-instantaneous two-point correlation functions, as well as the relaxation spectrum are obtained on a finite and periodic lattice.Comment: 5 pages, no figure. To appear as a Rapid Communication in Physical Review E (September 2001

Complete Solution of the Kinetics in a Far-from-equilibrium Ising Chain

The one-dimensional Ising model is easily generalized to a \textit{genuinely nonequilibrium} system by coupling alternating spins to two thermal baths at different temperatures. Here, we investigate the full time dependence of this system. In particular, we obtain the evolution of the magnetisation, starting with arbitrary initial conditions. For slightly less general initial conditions, we compute the time dependence of all correlation functions, and so, the probability distribution. Novel properties, such as oscillatory decays into the steady state, are presented. Finally, we comment on the relationship to a reaction-diffusion model with pair annihilation and creation.Comment: Submitted to J. Phys. A (Letter to the editor

When does cyclic dominance lead to stable spiral waves?

Species diversity in ecosystems is often accompanied by characteristic spatio-temporal patterns. Here, we consider a generic two-dimensional population model and study the spiraling patterns arising from the combined effects of cyclic dominance of three species, mutation, pair-exchange and individual hopping. The dynamics is characterized by nonlinear mobility and a Hopf bifurcation around which the system's four-phase state diagram is inferred from a complex Ginzburg-Landau equation derived using a perturbative multiscale expansion. While the dynamics is generally characterized by spiraling patterns, we show that spiral waves are stable in only one of the four phases. Furthermore, we characterize a phase where nonlinearity leads to the annihilation of spirals and to the spatially uniform dominance of each species in turn. Away from the Hopf bifurcation, when the coexistence fixed point is unstable, the spiraling patterns are also affected by the nonlinear diffusion

Coupled environmental and demographic fluctuations shape the evolution of cooperative antimicrobial resistance

There is a pressing need to better understand how microbial populations respond to antimicrobial drugs, and to find mechanisms to possibly eradicate antimicrobial-resistant cells. The inactivation of antimicrobials by resistant microbes can often be viewed as a cooperative behavior leading to the coexistence of resistant and sensitive cells in large populations and static environments. This picture is however greatly altered by the fluctuations arising in volatile environments, in which microbial communities commonly evolve. Here, we study the eco-evolutionary dynamics of a population consisting of an antimicrobial resistant strain and microbes sensitive to antimicrobial drugs in a time-fluctuating environment, modeled by a carrying capacity randomly switching between states of abundance and scarcity. We assume that antimicrobial resistance is a shared public good when the number of resistant cells exceeds a certain threshold. Eco-evolutionary dynamics is thus characterized by demographic noise (birth and death events) coupled to environmental fluctuations which can cause population bottlenecks. By combining analytical and computational means, we determine the environmental conditions for the long-lived coexistence and fixation of both strains, and characterize a fluctuation-driven antimicrobial resistance eradication mechanism, where resistant microbes experience bottlenecks leading to extinction. We also discuss the possible applications of our findings to laboratory-controlled experiments.Comment: 19+7 pages, 4+1 figures. Simulation data and codes for all figures are electronically available from the University of Leeds Data Repository. DOI: https://doi.org/10.5518/136