Nontrivial critical crossover between directed percolation models: effect of infinitely many absorbing states

Universality-class crossover by a nonorder field introduced to the pair contact process with diffusion

The one-dimensional pair contact process with diffusion (PCPD), an interacting particle system with diffusion, pair annihilation, and creation by pairs, has defied consensus about the universality class to which it belongs. An argument by Hinrichsen [Physica A 361, 457 (2006)PHYADX0378-437110.1016/j.physa.2005.06.101] claims that freely diffusing particles in the PCPD should play the same role as frozen particles when it comes to the critical behavior. Therefore, the PCPD is claimed to have the same critical phenomena as a model with infinitely many absorbing states that belongs to the directed percolation (DP) universality class. To investigate if diffusing particles are really indistinguishable from frozen particles in the sense of the renormalization group, we study numerically a variation of the PCPD by introducing a nonorder field associated with infinitely many absorbing states. We find that a crossover from the PCPD to DP occurs due to the nonorder field. By studying a similar model, we exclude the possibility that the mere introduction of a nonorder field to one model can entail a nontrivial crossover to another model in the same universality class, thus we attribute the observed crossover to the difference of the universality class of the PCPD from the DP class.

Critical decay exponent of the pair contact process with diffusion

We investigate the one-dimensional pair contact process with diffusion (PCPD) by extensive Monte Carlo simulations, mainly focusing on the critical density decay exponent δ. To obtain an accurate estimate of δ, we first find the strength of corrections to scaling using the recently introduced method [S.-C. Park. J. Korean Phys. Soc. 62, 469 (2013)KPSJAS0374-488410.3938/jkps.62.469]. For small diffusion rate (d≤0.5), the leading corrections-to-scaling term is found to be ∼t^{-0.15}, whereas for large diffusion rate (d=0.95) it is found to be ∼t^{-0.5}. After finding the strength of corrections to scaling, effective exponents are systematically analyzed to conclude that the value of critical decay exponent δ is 0.173(3) irrespective of d. This value should be compared with the critical decay exponent of the directed percolation, 0.1595. In addition, we study two types of crossover. At d=0, the phase boundary is discontinuous and the crossover from the pair contact process to the PCPD is found to be described by the crossover exponent ϕ=2.6(1). We claim that the discontinuity of the phase boundary cannot be consistent with the theoretical argument supporting the hypothesis that the PCPD should belong to the DP. At d=1, the crossover from the mean field PCPD to the PCPD is described by ϕ=2 which is argued to be exact.

Tuning spreading and avalanche-size exponents in directed percolation with modified activation probabilities

We consider the directed percolation process as a prototype of systems displaying a nonequilibrium phase transition into an absorbing state. The model is in a critical state when the activation probability is adjusted at some precise value p(c). Criticality is lost as soon as the probability to activate sites at the first attempt, p(1), is changed. We show here that criticality can be restored by "compensating" the change in p(1) by an appropriate change of the second time activation probability p(2) in the opposite direction. At compensation, we observe that the bulk exponents of the process coincide with those of the normal directed percolation process. However, the spreading exponents are changed and take values that depend continuously on the pair (p(1),p(2)). We interpret this situation by acknowledging that the model with modified initial probabilities has an infinite number of absorbing states.

Coexistence and critical behavior in a lattice model of competing species

In the present paper we study a lattice model of two species competing for the same resources. Monte Carlo simulations for d = 1,2, and 3 show that when resources are easily available both species coexist. However, when the supply of resources is on an intermediate level, the species with slower metabolism becomes extinct. On the other hand, when resources are scarce it is the species with faster metabolism that becomes extinct. The range of coexistence of the two species increases with dimension. We suggest that our model might describe some aspects of the competition between normal and tumor cells. With such an interpretation, examples of tumor remission, recurrence, and different morphologies are presented. In the d = 1 and d = 2 models, we analyze the nature of phase transitions: they are either discontinuous or belong to the directed-percolation universality class, and in some cases they have an active subcritical phase. In the d = 2 case, one of the transitions seems to be characterized by critical exponents that differ from directed-percolation ones, but this transition could be also weakly discontinuous. In the d = 3 version, Monte Carlo simulations are in a good agreement with the solution of the mean-field approximation. This approximation predicts that oscillatory behavior occurs in the present model but only for d ≳ 2. For d ≥ 2, a steady state depends on the initial configuration in some cases.

Generalized contact process with two symmetric absorbing states in two dimensions

We explore the two-dimensional generalized contact process with two absorbing states by means of large-scale Monte-Carlo simulations. In part of the phase diagram, an infinitesimal creation rate of active sites between inactive domains is sufficient to take the system from the inactive phase to the active phase. The system, therefore, displays two different nonequilibrium phase transitions. The critical behavior of the generic transition is compatible with the generalized voter universality class, implying that the symmetry-breaking and absorbing transitions coincide. In contrast, the transition at zero domain-boundary activation rate is not critical.

Phase transitions of the generalized contact process with two absorbing states

We investigate the generalized contact process with two absorbing states in one space dimension by means of large-scale Monte Carlo simulations. Treating the creation rate of active sites between inactive domains as an independent parameter leads to a rich phase diagram. In addition to the conventional active and inactive phases we find a parameter region where the simple contact process is inactive, but an infinitesimal creation rate at the boundary between inactive domains is sufficient to take the system into the active phase. Thus, the generalized contact process has two different phase transition lines. The point separating them shares some characteristics with a multicritical point. We also study in detail the critical behaviors of these transitions and their universality.

Absence of the discontinuous transition in the one-dimensional triplet creation model

Although Hinrichsen in his unpublished work theoretically rebutted the possibility of the discontinuous transition in one-dimensional nonequilibrium systems unless there are additional conservation laws, long-range interactions, macroscopic currents, or special boundary conditions, we have recently observed the resurrection of the claim that the triplet creation (TC) model introduced by Dickman and Tomé [Phys. Rev. A 44, 4833 (1991)] would show the discontinuous transition. By extensive simulations, however, we find that the one-dimensional TC does belong to the directed percolation universality class even for larger diffusion constant than the suggested tricritical point in the literature. Furthermore, we find that the phase boundary is well described by the crossover from the mean field to the directed percolation, which supports the claim that the one-dimensional TC does not exhibit a discontinuous transition.

Crossover from the parity-conserving pair contact process with diffusion to other universality classes

The pair contact process with diffusion (PCPD) with modulo 2 conservation (PCPD2) [ 2A-->4A , 2A-->0 ] is studied in one dimension, focused on the crossover to other well established universality classes: the directed Ising (DI) and the directed percolation (DP). First, we show that the PCPD2 shares the critical behaviors with the PCPD, both with and without directional bias. Second, the crossover from the PCPD2 to the DI is studied by including a parity-conserving single-particle process (A-->3A) . We find the crossover exponent 1/varphi_{1}=0.57(3) , which is argued to be identical to that of the PCPD-to-DP crossover by adding A-->2A . This suggests that the PCPD universality class has a well-defined fixed point distinct from the DP. Third, we study the crossover from a hybrid-type reaction-diffusion process belonging to the DP [ 3A-->5A , 2A-->0 ] to the DI by adding A-->3A . We find 1/varphi_{2}=0.73(4) for the DP-to-DI crossover. The inequality of varphi_{1} and varphi_{2} further supports the non-DP nature of the PCPD scaling. Finally, we introduce a symmetry-breaking field in the dual spin language to study the crossover from the PCPD2 to the DP. We find 1/varphi_{3}=1.23(10) , which is associated with a new independent route from the PCPD to the DP.

Crossovers from parity conserving to directed percolation universality

The crossover behavior of various models exhibiting phase transition to absorbing phase with parity conserving class has been investigated by numerical simulations and cluster mean-field method. In case of models exhibiting Z_2 symmetric absorbing phases (the cellular automaton version of the nonequilibrium kinetic Ising model (NEKIMCA) and a stochastic cellular automaton invented by Grassberger, Krause, and von der Twer [J. Phys. A 17, L105 (1984)]) the introduction of an external symmetry breaking field causes a crossover to kink parity conserving models characterized by dynamical scaling of the directed percolation (DP) and the crossover exponent: 1/phi approximately equal to 0.53(2) . In the case of a branching and annihilating random walk model with an even number of offspring (dual to NEKIMCA) the introduction of spontaneous particle decay destroys the parity conservation and results in a crossover to the DP class characterized by the crossover exponent: 1/phi approximately equal to 0.205(5) . The two different kinds of crossover operators cannot be mapped onto each other and the resulting models show a diversity within the DP universality class in one dimension. These subclasses differ in cluster scaling exponents.

Three different routes from the directed Ising to the directed percolation class

The scaling nature of absorbing critical phenomena is well understood for the directed percolation (DP) and the directed Ising (DI) systems. However, a full analysis of the crossover behavior is still lacking, which is of our interest in this study. In one dimension, we find three different routes from the DI to the DP classes by introducing a symmetry-breaking field (SB), breaking a modulo 2 conservation (CB), or making channels connecting two equivalent absorbing states (CC). Each route can be characterized by a crossover exponent, which is found numerically as phi=2.1+/-0.1 (SB), 4.6+/-0.2 (CB), and 2.9+/-0.1 (CC), respectively. The difference between the SB and CB crossover can be understood easily in the domain wall language, while the CC crossover involves an additional critical singularity in the auxiliary field density with the memory effect to identify itself independent.

Generic absorbing transition in coevolution dynamics

We study a coevolution voter model on a complex network. A mean-field approximation reveals an absorbing transition from an active to a frozen phase at a critical value [see text for formula] that only depends on the average degree micro of the network. In finite-size systems, the active and frozen phases correspond to a connected and a fragmented network, respectively. The transition can be seen as the sudden change in the trajectory of an equivalent random walk at the critical point, resulting in an approach to the final frozen state whose time scale diverges as tau approximately |p(c) - p|(-)} near p(c).