Universality class of absorbing transitions with continuously varying critical exponents

Transition to period-3 synchronized state in coupled gauss maps

We study coupled Gauss maps in one dimension with nearest-neighbor interactions. We observe transitions from spatiotemporal chaos to period-3 states in a coarse-grained sense and synchronized period-3 states. Synchronized fixed points are frequently observed in one dimension. However, synchronized periodic states are rare. The obvious reason is that it is very easy to create defects in one dimension. We characterize all transitions using the following order parameter. Let x∗ be the fixed point of the map. The values above (below) x∗ are classified as +1 (-1) spins. We expect all sites to return to the same band after three time steps for a coarse-grained periodic or three-period state. We define the flip rate F(t) as the fraction of sites i such that si(3t-3)≠si(t). It is zero in the coarse-grained periodic state. This state may or may not be synchronized. We observe three different transitions. (a) If different sites reach different bands, the transition is in the directed-percolation universality class. (b) If all sites reach the same band, we find an Ising-type transition. (c) A synchronized period-3 state where a new exponent is observed. We also study the finite-size scaling at critical points. The exponents obtained indicate that the synchronized period-3 transition is in a new universality class.

Subdiffusive Activity Spreading in the Diffusive Epidemic Process

The diffusive epidemic process is a paradigmatic example of an absorbing state phase transition in which healthy and infected individuals spread with different diffusion constants. Using stochastic activity spreading simulations in combination with finite-size scaling analyses we reveal two qualitatively different processes that characterize the critical dynamics: subdiffusive propagation of infection clusters and diffusive fluctuations in the healthy population. This suggests the presence of a strong-coupling regime and sheds new light on a long-standing debate about the theoretical classification of the system.

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.

Full nonuniversality of the symmetric 16-vertex model on the square lattice

We consider the symmetric two-state 16-vertex model on the square lattice whose vertex weights are invariant under any permutation of adjacent edge states. The vertex-weight parameters are restricted to a critical manifold which is self-dual under the gauge transformation. The critical properties of the model are studied numerically with the Corner Transfer Matrix Renormalization Group method. Accuracy of the method is tested on two exactly solvable cases: the Ising model and a specific version of the Baxter eight-vertex model in a zero field that belong to different universality classes. Numerical results show that the two exactly solvable cases are connected by a line of critical points with the polarization as the order parameter. There are numerical indications that critical exponents vary continuously along this line in such a way that the weak universality hypothesis is violated.

Universality of the local persistence exponent for models in the directed Ising class in one dimension

We investigate local persistence in five different models and their variants in the directed Ising (DI) universality class in one dimension. These models have right-left symmetry. We study Grassberger's models A and B. We also study branching and annihilating random walks with two offspring: the nonequilibrium kinetic Ising model and the interacting monomer-dimer model. Grassberger's models are updated in parallel. This is not the case in other models. We find that the local persistence exponent in all these models is unity or very close to it. A change in the mode of the update does not change the exponent unless the universality class changes. In general, persistence exponents are not universal. Thus it is of interest that the persistence exponent in a range of models in the DI class is the same. Excellent scaling behavior of finite-size scaling is obtained using exponents in the DI class in all models. We also study off-critical scaling in some models and DI exponents yield excellent scaling behavior. We further study graded persistence, which shows similar behavior. However, for a logistic map with delay, which also has the transition in the DI class, there is no transition from nonzero to zero persistence at the critical point. Thus the accompanying transition in persistence and universality of the persistence exponent hold when the underlying model has right-left symmetry.

Universality classes of absorbing phase transitions in generic branching-annihilating particle systems with nearest-neighbor bias

We study absorbing phase transitions in systems of branching annihilating random walkers and pair contact process with diffusion on a one-dimensional ring, where the walkers hop to their nearest neighbor with a bias ε. For ε=0, three universality classes-directed percolation (DP), parity-conserving (PC), and pair contact process with diffusion (PCPD)-are typically observed in such systems. We find that the introduction of ε does not change the DP universality class but alters the other two universality classes. For nonzero ε, the PCPD class crosses over to DP, and the PC class changes to a new universality class.

Clearing out a maze: A model of chemotactic motion in porous media

We study the anomalous dynamics of a biased "hungry" (or "greedy") random walk on a percolating cluster. The model mimics chemotaxis in a porous medium: In close resemblance to the 1980s arcade game PAC-MAN^®, the hungry random walker consumes food, which is initially distributed in the maze, and biases its movement towards food-filled sites. We observe that the mean-squared displacement of the process follows a power law with an exponent that is different from previously known exponents describing passive or active microswimmer dynamics. The change in dynamics is well described by a dynamical exponent that depends continuously on the propensity to move towards food. It results in slower differential growth when compared to the unbiased random walk.

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.

Continuously Varying Critical Exponents Beyond Weak Universality

Renormalization group theory does not restrict the form of continuous variation of critical exponents which occurs in presence of a marginal operator. However, the continuous variation of critical exponents, observed in different contexts, usually follows a weak universality scenario where some of the exponents (e.g., β, γ, ν) vary keeping others (e.g., δ, η) fixed. Here we report ferromagnetic phase transition in (Sm_1−yNd_y)_0.52Sr_0.48MnO₃ (0.5 ≤ y ≤ 1) single crystals where all three exponents β, γ, δ vary with Nd concentration y. Such a variation clearly violates both universality and weak universality hypothesis. We propose a new scaling theory that explains the present experimental results, reduces to the weak universality as a special case, and provides a generic route leading to continuous variation of critical exponents and multi-criticality.

Stationary and dynamic critical behavior of the contact process on the Sierpinski carpet

We investigate the critical behavior of a stochastic lattice model describing a contact process in the Sierpinski carpet with fractal dimension d=log8/log3. We determine the threshold of the absorbing phase transition related to the transition between a statistically stationary active and the absorbing states. Finite-size scaling analysis is used to calculate the order parameter, order parameter fluctuations, correlation length, and their critical exponents. We report that all static critical exponents interpolate between the line of the regular Euclidean lattices values and are consistent with the hyperscaling relation. However, a short-time dynamics scaling analysis shows that the dynamical critical exponent Z governing the size dependence of the critical relaxation time is found to be larger then the literature values in Euclidean d=1 and d=2, suggesting a slower critical relaxation in scale-free lattices.

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.

Universality classes of the absorbing state transition in a system with interacting static and diffusive populations

In this work, we study the critical behavior of a one-dimensional model that mimics the propagation of an epidemic process mediated by a density of diffusive individuals which can infect a static population upon contact. We simulate the above model on linear chains to determine the critical density of the diffusive population, above which the system achieves a statistically stationary active state, as a function of two relevant parameters related to the average lifetimes of the diffusive and nondiffusive populations. A finite-size scaling analysis is employed to determine the order parameter and correlation length critical exponents. For high-recovery rates, the critical exponents are compatible with the usual directed percolation universality class. However, in the opposite regime of low-recovery rates, the diffusion is a relevant mechanism responsible for the propagation of the disease and the absorbing state phase transition is governed by a distinct set of critical exponents.

Discontinuous phase transitions of conserved threshold transfer process with deterministic hopping

The deterministic conserved threshold transfer process, which is a variant of the conserved threshold transfer process modified in a way as that hopping of a particle is to be deterministic, is proposed. The critical behavior of the model is investigated in one, two, and four dimensions. It is found that the order parameter yields a discontinuous transition; i.e., the transition appears to be first ordered in all dimensions considered. The origin of such a discontinuous transition is investigated, considering clustering of active sites and accumulation of critical sites just before the steady state is reached.

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.

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.

Universality split in absorbing phase transition with conserved field on fractal lattices

The universality split in absorbing phase transition between the conserved lattice gas (CLG) model and the conserved threshold transfer process (CTTP) is investigated on a checkerboard fractal and on a Sierpinski gasket. The critical exponents theta, beta, nu||, and z, which are associated with, respectively, the density of active particles in time, the order parameter, the temporal correlation length, and the dynamics of active particles, are elaborately measured for two models on selected fractal lattices. The exponents for the CLG model are found to be distinctly different from those of the CTTP model on a checkerboard fractal, whereas the two models exhibit the same critical behavior on a Sierpinski gasket, indicating that the universality split between the two models occurs only on a checkerboard fractal. Such a universality split is attributed from the dominant hopping mechanisms caused by the intrinsic properties of the underlying fractal lattice.

Absorbing phase transition in a conserved lattice gas model in one dimension

We investigated the nonequilibrium phase transition of the conserved lattice gas model in one dimension using two update methods: i.e., the sequential update and the parallel update. We measured the critical indices of theta, beta, nu(parallel), and nu(perpendicular), and found that, for a parallel update, the exponents were delicately influenced by the hopping rule of active particles. When the hopping rule was designed to be symmetric, the results were found to be consistent with those of the sequential update. The exponents we obtained were precisely the same as the corresponding results of a recently presented lattice model of two species of particles with a conserved field in one dimension, in contrast with the authors' claim. We also found that one of the scaling relations known for absorbing phase transition is violated.

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

At nonequilibrium phase transitions into absorbing (trapped) states, it is well known that the directed percolation (DP) critical scaling is shared by two classes of models with a single (S) absorbing state and with infinitely many (IM) absorbing states. We study the crossover behavior in one dimension, arising from a considerable reduction of the number of absorbing states (typically from the IM-type to the S -type DP models) by following two different (excitatory or inhibitory) routes which make the auxiliary field density abruptly jump at the crossover. Along the excitatory route, the system becomes overly activated even for an infinitesimal perturbation and its crossover becomes discontinuous. Along the inhibitory route, we find a continuous crossover with universal crossover exponent phi approximately=1.78(6), which is argued to be equal to nu||, the relaxation time exponent of the DP universality class on a general footing. This conjecture is also confirmed in the case of the directed Ising (parity-conserving) class. Finally, we discuss the effect of diffusion on the IM-type models and suggest an argument why diffusive models with some hybrid-type reactions should belong to the DP class.

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.

Moment ratios for the pair-contact process with diffusion

We study the continuous absorbing-state phase transition in the one-dimensional pair contact process with diffusion (PCPD). In previous studies [Dickman and de Menezes, Phys. Rev. E 66, 045101(R) (2002)], the critical point moment ratios of the order parameter showed anomalous behavior, growing with system size rather than taking universal values, as expected. Using the quasistationary simulation method we determine the moments of the order parameter up to fourth order at the critical point, in systems of up to 40,960 sites. Due to strong finite-size effects, the ratios converge only for large system sizes. Moment ratios and associated order-parameter histograms are compared with those of directed percolation. We also report an improved estimate [pc=0.077092(1)] for the location of the critical point in the nondiffusive pair contact process.

Kinetic roughening and pinning of two coupled interfaces in disordered media

We studied the kinetic roughening dynamic of two coupled interfaces formed in paper wetting experiments at low evaporation rate. We observed three different regimes of impregnation in which kinetic roughening dynamics of coupled precursor and main fronts belong to different universality classes; nevertheless both interfaces are pinned in the same configuration. Reported experimental observations provide a novel insight into the nature of kinetic roughening phenomena occurring in the vast variety of systems far from equilibrium.

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.

Anomalous roughness with system-size-dependent local roughness exponent

We note that in a system far from equilibrium the interface roughening may depend on the system size which plays the role of control parameter. To detect the size effect on the interface roughness, we study the scaling properties of rough interfaces formed in paper combustion experiments. Using paper sheets of different width lambda L0, we found that the turbulent flame fronts display anomalous multiscaling characterized by nonuniversal global roughness exponent alpha and the system-size-dependent spectrum of local roughness exponents zeta(q) (lambda) = zeta(1) (1) q(-omega) lambda(phi) <alpha, whereas the burning fronts possess conventional multiaffine scaling. The structure factor of turbulent flame fronts also exhibit unconventional scaling dependence on lambda. These results are expected to apply to a broad range of far from equilibrium systems, when the kinetic energy fluctuations exceed a certain critical value.

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.

Critical behavior of the two-dimensional 2A-->3A, 4A--> phi binary system

The phase transitions of the recently introduced 2A-->3A, 4A--> phi reaction-diffusion model [G. Odor, Phys. Rev. E 69, 036112 (2004)]] are explored in two dimensions. This model exhibits site-occupation restriction and explicit diffusion of isolated particles. A reentrant phase diagram in the diffusion-creation rate space is confirmed, in agreement with cluster mean-field and one-dimensional results. For strong diffusion, a mean-field transition can be observed at zero branching rate characterized by an alpha=1/3 density decay exponent. In contrast, for weak diffusion the effective 2A-->3A-->4A--> phi reaction becomes relevant and the mean-field transition of the 2A-->3A, 2A--> phi model characterized by alpha=1/2 also appears for nonzero branching rates.

Phase transitions of the binary production 2A-->3A, 4A--> X model

Phase transitions of the 2A-->3A, 4A--> X reaction-diffusion model is explored by dynamical, N-cluster approximations and by simulations. The model exhibits site occupation restriction and explicit diffusion of isolated particles. While the site mean-field approximation shows a single transition at zero branching rate introduced by Odor [G. Odor, Phys. Rev. E 67, 056114 (2003)], N>2 cluster approximations predict the appearance of another transition line for weak diffusion (D) as well. The latter phase transition is continuous, occurs at finite branching rate, and exhibits different scaling behavior. I show that the universal behavior of these transitions is in agreement with that of the diffusive pair contact process model both on the mean-field level and in one dimension. Therefore this model exhibiting annihilation by quadruplets does not fit in the recently suggested classification of universality classes of absorbing state transitions in one dimension [J. Kockelkoren and H. Chaté, Phys. Rev. Lett. 90, 125701 (2003)]. For high diffusion rates the effective 2A-->3A-->4A--> X reaction becomes irrelevant and the model exhibits a mean-field transition only. The two regions are separated by a nontrivial critical end point at D*.