Non-equilibrium relaxation in a stochastic lattice Lotka–Volterra model

Computing macroscopic reaction rates in reaction-diffusion systems using Monte Carlo simulations

Stochastic reaction-diffusion models are employed to represent many complex physical, biological, societal, and ecological systems. The macroscopic reaction rates describing the large-scale, long-time kinetics in such systems are effective, scale-dependent renormalized parameters that need to be either measured experimentally or computed by means of a microscopic model. In a Monte Carlo simulation of stochastic reaction-diffusion systems, microscopic probabilities for specific events to happen serve as the input control parameters. To match the results of any computer simulation to observations or experiments carried out on the macroscale, a mapping is required between the microscopic probabilities that define the Monte Carlo algorithm and the macroscopic reaction rates that are experimentally measured. Finding the functional dependence of emergent macroscopic rates on the microscopic probabilities (subject to specific rules of interaction) is a very difficult problem, and there is currently no systematic, accurate analytical way to achieve this goal. Therefore, we introduce a straightforward numerical method of using lattice Monte Carlo simulations to evaluate the macroscopic reaction rates by directly obtaining the count statistics of how many events occur per simulation time step. Our technique is first tested on well-understood fundamental examples, namely, restricted birth processes, diffusion-limited two-particle coagulation, and two-species pair annihilation kinetics. Next we utilize the thus gained experience to investigate how the microscopic algorithmic probabilities become coarse-grained into effective macroscopic rates in more complex model systems such as the Lotka-Volterra model for predator-prey competition and coexistence, as well as the rock-paper-scissors or cyclic Lotka-Volterra model and its May-Leonard variant that capture population dynamics with cyclic dominance motifs. Thereby we achieve a more thorough and deeper understanding of coarse graining in spatially extended stochastic reaction-diffusion systems and the nontrivial relationships between the associated microscopic and macroscopic model parameters, with a focus on ecological systems. The proposed technique should generally provide a useful means to better fit Monte Carlo simulation results to experimental or observational data.

Coupled two-species model for the pair contact process with diffusion

The contact process with diffusion (PCPD) defined by the binary reactions B+B→B+B+B, B+B→∅ and diffusive particle spreading exhibits an unusual active to absorbing phase transition whose universality class has long been disputed. Multiple studies have indicated that an explicit account of particle pair degrees of freedom may be required to properly capture this system's effective long-time, large-scale behavior. We introduce a two-species representation for the PCPD in which single particles B and particle pairs A are dynamically coupled according to the stochastic reaction processes B+B→A, A→A+B, A→∅, and A→B+B, with each particle type diffusing independently. Mean-field analysis reveals that the phase transition of this model is driven by competition and balance between the two species. We employ Monte Carlo simulations in one, two, and three dimensions to demonstrate that this model consistently captures the pertinent features of the PCPD. In the inactive phase, A particles rapidly go extinct, effectively leaving the B species to undergo pure diffusion-limited pair annihilation kinetics B+B→∅. At criticality, both A and B densities decay with the same exponents (within numerical errors) as the corresponding order parameters of the original PCPD, and display mean-field scaling above the upper critical dimension d_{c}=2. In one dimension, the critical exponents for the B species obtained from seed simulations also agree well with previously reported exponent value ranges. We demonstrate that the scaling properties of consecutive particle pairs in the PCPD are identical with that of the A species in the coupled model. This two-species picture resolves the conceptual difficulty for seed simulations in the original PCPD and naturally introduces multiple length scales and timescales to the system, which are also the origin of strong corrections to scaling. The extracted moment ratios from our simulations indicate that our model displays the same temporal crossover behavior as the PCPD, which further corroborates its full dynamical equivalence with our coupled model.

Aging and equilibration in bistable contagion dynamics

We analyze the late-time relaxation dynamics for a general contagion model. In this model, nodes are either active or failed. Active nodes can fail either "spontaneously" at any time or "externally" if their neighborhoods are sufficiently damaged. Failed nodes may always recover spontaneously. At late times, the breaking of time-translation invariance is a necessary condition for physical aging. We observe that time-translational invariance is lost for initial conditions that lie between the basins of attraction of the model's two stable stationary states. Based on corresponding mean-field predictions, we characterize the observed model behavior in terms of a phase diagram spanned by the fractions of spontaneously and externally failed nodes. For the square lattice, the phases in which the dynamics approaches one of the two stable stationary states are not linearly separable due to spatial correlation effects. Our results provide insights into aging and relaxation phenomena that are observable in a model of social contagion processes.

Critical dynamics of anisotropic antiferromagnets in an external field

We numerically investigate the nonequilibrium critical dynamics in three-dimensional anisotropic antiferromagnets in the presence of an external magnetic field. The phase diagram of this system exhibits two critical lines that meet at a bicritical point. The nonconserved components of the staggered magnetization order parameter couple dynamically to the conserved component of the magnetization density along the direction of the external field. Employing a hybrid computational algorithm that combines reversible spin precession with relaxational Monte Carlo updates, we study the aging scaling dynamics for the model C critical line, identifying the critical initial slip, autocorrelation, and aging exponents for both the order parameter and the conserved field, thus also verifying the dynamic critical exponent. We further probe the model F critical line by investigating the system size dependence of the characteristic spin wave frequencies near criticality and measure the dynamic critical exponents for the order parameter including its aging scaling at the bicritical point.

Analytical treatment for cyclic three-state dynamics on static networks

Whenever a dynamical process unfolds on static networks, the dynamical state of any focal individual will be exclusively influenced by directly connected neighbors, rather than by those unconnected ones, hence the arising of the dynamical correlation problem, where mean-field-based methods fail to capture the scenario. The dynamic correlation coupling problem has always been an important and difficult problem in the theoretical field of physics. The explicit analytical expressions and the decoupling methods often play a key role in the development of corresponding field. In this paper, we study the cyclic three-state dynamics on static networks, which include a wide class of dynamical processes, for example, the cyclic Lotka-Volterra model, the directed migration model, the susceptible-infected-recovered-susceptible epidemic model, and the predator-prey with empty sites model. We derive the explicit analytical solutions of the propagating size and the threshold curve surface for the four different dynamics. We compare the results on static networks with those on annealed networks and made an interesting discovery: for the symmetrical dynamical model (the cyclic Lotka-Volterra model and the directed migration model, where the three states are of rotational symmetry), the macroscopic behaviors of the dynamical processes on static networks are the same as those on annealed networks; while the outcomes of the dynamical processes on static networks are different with, and more complicated than, those on annealed networks for asymmetric dynamical model (the susceptible-infected-recovered-susceptible epidemic model and the predator-prey with empty sites model). We also compare the results forecasted by our theoretical method with those by Monte Carlo simulations and find good agreement between the results obtained by the two methods.

Lineage space and the propensity of bacterial cells to undergo growth transitions

The molecular makeup of the offspring of a dividing cell gradually becomes phenotypically decorrelated from the parent cell by noise and regulatory mechanisms that amplify phenotypic heterogeneity. Such regulatory mechanisms form networks that contain thresholds between phenotypes. Populations of cells can be poised near the threshold so that a subset of the population probabilistically undergoes the phenotypic transition. We sought to characterize the diversity of bacterial populations around a growth-modulating threshold via analysis of the effect of non-genetic inheritance, similar to conditions that create antibiotic-tolerant persister cells and other examples of bet hedging. Using simulations and experimental lineage data in Escherichia coli, we present evidence that regulation of growth amplifies the dependence of growth arrest on cellular lineage, causing clusters of related cells undergo growth arrest in certain conditions. Our simulations predict that lineage correlations and the sensitivity of growth to changes in toxin levels coincide in a critical regime. Below the critical regime, the sizes of related growth arrested clusters are distributed exponentially, while in the critical regime clusters sizes are more likely to become large. Furthermore, phenotypic diversity can be nearly as high as possible near the critical regime, but for most parameter values it falls far below the theoretical limit. We conclude that lineage information is indispensable for understanding regulation of cellular growth.