Universality in the pair contact process with diffusion

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.

Stochastic dynamics and logistic population growth

The Verhulst model is probably the best known macroscopic rate equation in population ecology. It depends on two parameters, the intrinsic growth rate and the carrying capacity. These parameters can be estimated for different populations and are related to the reproductive fitness and the competition for limited resources, respectively. We investigate analytically and numerically the simplest possible microscopic scenarios that give rise to the logistic equation in the deterministic mean-field limit. We provide a definition of the two parameters of the Verhulst equation in terms of microscopic parameters. In addition, we derive the conditions for extinction or persistence of the population by employing either the momentum-space spectral theory or the real-space Wentzel-Kramers-Brillouin approximation to determine the probability distribution function and the mean time to extinction of the population. Our analytical results agree well with numerical simulations.

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.

Finite-scale singularity in the renormalization group flow of a reaction-diffusion system

We study the nonequilibrium critical behavior of the pair contact process with diffusion (PCPD) by means of nonperturbative functional renormalization group techniques. We show that usual perturbation theory fails because the effective potential develops a nonanalyticity at a finite length scale: Perturbatively forbidden terms are dynamically generated and the flow can be continued once they are taken into account. Our results suggest that the critical behavior of PCPD can be either in the directed percolation or in a different (conjugated) universality class.

Percolation transition at growing spatiotemporal fractal patterns in models of mesoscopic neural networks

Spike packet propagation is modeled in mesoscopic-scale networks, composed of locally and recurrently coupled neural pools, and embedded in a two-dimensional lattice. Site dynamics is governed by three key parameters--pool connectedness probability, synaptic strength (following the steady-state distribution of some realizations of spike-timing-dependent plasticity learning rule), and the neuron refractoriness. Formation of spatiotemporal patterns in our model, synfire chains, exhibits critical behavior, with the emerging percolation phase transition controlled by the probability for nonzero synaptic strength value. Applying the finite-size scaling method, we infer the critical probability dependence on synaptic strength and refractoriness and determine the effects of connection topology on the pertaining percolation clusters fractal dimensions. We find that the directed percolation and the pair contact process with diffusion constitute the relevant universality classes of phase transitions observed in a class of mesoscopic-scale network models, which may be related to recently reported data on in vitro cultures.

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.

Universality class of the pair contact process with diffusion

The pair contact process with diffusion (PCPD) is studied with a standard Monte Carlo approach and with simulations at fixed densities. A standard analysis of the simulation results, based on the particle densities or on the pair densities, yields inconsistent estimates for the critical exponents. However, if a well-chosen linear combination of the particle and pair densities is used, leading corrections can be suppressed, and consistent estimates for the independent critical exponents delta=0.16(2) , beta=0.28(2) , and z=1.58 are obtained. Since these estimates are also consistent with their values in directed percolation (DP), we conclude that the PCPD falls in the same universality class as DP.

Classification of phase transitions in reaction-diffusion models

Equilibrium phase transitions are associated with rearrangements of minima of a (Lagrangian) potential. Treatment of nonequilibrium systems requires doubling of degrees of freedom, which may be often interpreted as a transition from the "coordinate"- to the "phase"-space representation. As a result, one has to deal with the Hamiltonian formulation of the field theory instead of the Lagrangian one. We suggest a classification scheme of phase transitions in reaction-diffusion models based on the topology of the phase portraits of corresponding Hamiltonians. In models with an absorbing state such a topology is fully determined by intersecting curves of zero "energy." We identify four families of topologically distinct classes of phase portraits stable upon renormalization group transformations.

Crossover from the pair contact process with diffusion to directed percolation

Crossover behaviors from the pair contact process with diffusion (PCPD) and the driven PCPD (DPCPD) to the directed percolation (DP) are studied in one dimension by introducing a single particle annihilation and/or branching dynamics. The crossover exponents phi are estimated numerically as 1/phi approximately 0.58 +/- 0.03 for the PCPD and 1/phi approximately 0.49+/-0.02 for the DPCPD. Nontriviality of the PCPD crossover exponent strongly supports the non-DP nature of the PCPD critical scaling, which is further evidenced by the anomalous critical amplitude scaling near the PCPD point. In addition, we find that the DPCPD crossover is consistent with the mean field prediction of the tricritical DP class as expected.

Driven pair contact process with diffusion

The pair contact process with diffusion (PCPD) has been recently investigated extensively, but its critical behavior is not yet clearly established. By introducing biased diffusion, we show that the external driving is relevant and the driven PCPD exhibits a mean-field-type critical behavior even in one dimension. In systems which can be described by a single-species bosonic field theory, the Galilean invariance guarantees that the driving is irrelevant. The well-established directed percolation (DP) and parity-conserving (PC) classes are such examples. This leads us to conclude that the PCPD universality class should be distinct from the DP or the PC class. Moreover, it implies that the PCPD is generically a multispecies model and a field theory of two species is suitable for proper description.

Cluster mean-field approximations with the coherent-anomaly-method analysis for the driven pair contact process with diffusion

The cluster mean-field approximations are performed, up to 13 cluster sizes, to study the critical behavior of the driven pair contact process with diffusion (DPCPD) and its precedent, the PCPD in one dimension. Critical points are estimated by extrapolating our data to the infinite cluster size limit, which are in good accordance with recent simulation results. Within the cluster mean-field approximation scheme, the PCPD and the DPCPD share the same mean-field critical behavior. The application of the coherent anomaly method, however, shows that the two models develop different coherent anomalies, which lead to different true critical scaling. The values of the critical exponents for the particle density, the pair density, the correlation length, and the relaxation time are fairly well estimated for the DPCPD. These results support and complement our recent simulation results for the DPCPD.

Pair contact process with diffusion: failure of master equation field theory

We demonstrate that the "microscopic" field theory representation, directly derived from the corresponding master equation, fails to adequately capture the continuous nonequilibrium phase transition of the pair contact process with diffusion (PCPD). The ensuing renormalization group (RG) flow equations do not allow for a stable fixed point in the parameter region that is accessible by the physical initial conditions. There exists a stable RG fixed point outside this regime, but the resulting scaling exponents, in conjunction with the predicted particle anticorrelations at the critical point, would be in contradiction with the positivity of the equal-time mean-square particle number fluctuations. We conclude that a more coarse-grained effective field theory approach is required to elucidate the critical properties of the PCPD.

Universality class of absorbing transitions with continuously varying critical exponents

The well-established universality classes of absorbing critical phenomena are directed percolation (DP) and directed Ising (DI) classes. Recently, the pair contact process with diffusion (PCPD) has been investigated extensively and claimed to exhibit a different type of critical phenomenon distinct from both DP and DI classes. Noticing that the PCPD possesses a long-term memory effect, we introduce a generalized version of the PCPD (GPCPD) with a parameter controlling the memory strength. The GPCPD connects the DP fixed point to the PCPD point continuously. Monte Carlo simulations strongly suggest that the GPCPD displays, to our knowledge, novel critical phenomena which are characterized by continuously varying critical exponents. The same critical behaviors are also observed in models where two species of particles are coupled cyclically. We present one possible scenario that the long-term memory may serve as a marginal perturbation to the ordinary DP fixed point.