Moment ratios for the pair-contact process with diffusion

Machine learning of pair-contact process with diffusion

The pair-contact process with diffusion (PCPD), a generalized model of the ordinary pair-contact process (PCP) without diffusion, exhibits a continuous absorbing phase transition. Unlike the PCP, whose nature of phase transition is clearly classified into the directed percolation (DP) universality class, the model of PCPD has been controversially discussed since its infancy. To our best knowledge, there is so far no consensus on whether the phase transition of the PCPD falls into the unknown university classes or else conveys a new kind of non-equilibrium phase transition. In this paper, both unsupervised and supervised learning are employed to study the PCPD with scrutiny. Firstly, two unsupervised learning methods, principal component analysis (PCA) and autoencoder, are taken. Our results show that both methods can cluster the original configurations of the model and provide reasonable estimates of thresholds. Therefore, no matter whether the non-equilibrium lattice model is a random process of unitary (for instance the DP) or binary (for instance the PCP), or whether it contains the diffusion motion of particles, unsupervised learning can capture the essential, hidden information. Beyond that, supervised learning is also applied to learning the PCPD at different diffusion rates. We proposed a more accurate numerical method to determine the spatial correlation exponent [Formula: see text], which, to a large degree, avoids the uncertainty of data collapses through naked eyes.

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.

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.

Critical properties of the Ziff, Gulari, and Barshad (ZGB) model with inert sites

We study the model proposed by Ziff, Gulari, and Barshad to mimic the oxidation of carbon monoxide (CO) in the presence of fixed impurities distributed over the catalytic surface. Our focus is on the continuous phase transition between the active phase, where occurs the production of carbon dioxide (CO2), and the inactive phase, where all the non inert sites become filled with oxygen molecules. We employ Monte Carlo simulations to calculate the different ratios between moments of the order parameter at the critical point, as well as, we determine the critical exponents β and ν⊥ as a function of the concentration of impurities. We show that the presence of impurities over the catalytic surface changes the critical behavior of the system. The critical exponents depend on the concentration of impurities and the model does not belong to the directed percolation universality class.

Critical short-time dynamics in a system with interacting static and diffusive populations

We study the critical short-time dynamical behavior of a one-dimensional model where diffusive individuals can infect a static population upon contact. The model presents an absorbing phase transition from an active to an inactive state. Previous calculations of the critical exponents based on quasistationary quantities have indicated an unusual crossover from the directed percolation to the diffusive contact process universality classes. Here we show that the critical exponents governing the slow short-time dynamic evolution of several relevant quantities, including the order parameter, its relative fluctuations, and correlation function, reinforce the lack of universality in this model. Accurate estimates show that the critical exponents are distinct in the regimes of low and high recovery rates.

Moment ratios and dynamic critical behavior of a reactive system with several absorbing configurations

We determine the critical behavior of a reactive model with many absorbing configurations. Monomers A and B land on the sites of a linear lattice and can react depending on the state of their nearest-neighbor sites. The probability of a reaction depends on temperature of the catalyst as well as on the energy coupling between pairs of nearest-neighbor monomers. We employ Monte Carlo simulations to calculate the moments of the order parameter of the model as a function of temperature. Some ratios between pairs of moments are independent of temperature and are in the same universality class of the contact process. We also find the dynamical critical exponents of the model and we show that they are in the directed percolation universality class whatever the values of temperature.

Influence of zero range process interaction on diffusion

We study the aspects of diffusion for the case of zero range process interaction on scale-free networks, through statistical quantities such as the mean first passage time, coverage, mean square displacement etc., and pay attention to how the interaction, especially the resulted condensation, influences the diffusion. By mean-field theory we show that the statistical quantities of diffusion can be significantly reduced by the condensation and can be figured out by the waiting time of a particle staying at a node. Numerical simulations have confirmed the theoretical predictions.

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.

Contact process on a Voronoi triangulation

We study the continuous absorbing-state phase transition in the contact process on the Voronoi-Delaunay lattice. The Voronoi construction is a natural way to introduce quenched coordination disorder in lattice models. We simulate the disordered system using the quasistationary simulation method and determine its critical exponents and moment ratios. Our results suggest that the critical behavior of the disordered system is unchanged with respect to that on a regular lattice, i.e., that of directed percolation.

Absorbing-state phase transitions: exact solutions of small systems

I derive precise results for absorbing-state phase transitions using exact (numerically determined) quasistationary (QS) probability distributions for small systems. Analysis of the contact process on rings of 23 or fewer sites yields critical properties (control parameter, order-parameter ratios, and critical exponents z and betanu_(perpendicular)) with an accuracy of 0.06% or better; precise results are also obtained for the pair contact process. The QS distribution yields insights on the statistical entropy of these models. Preliminary application to a model in the stochastic sandpile class is also described.

Contact process with long-range interactions: a study in the ensemble of constant particle number

We analyze the properties of the contact process with long-range interactions by the use of a kinetic ensemble in which the total number of particles is strictly conserved. In this ensemble, both annihilation and creation processes are replaced by a unique process in which a particle of the system chosen at random leaves its place and jumps to an active site. The present approach is particularly useful for determining the transition point and the nature of the transition, whether continuous or discontinuous, by evaluating the fractal dimension of the cluster at the emergence of the phase transition. We also present another criterion appropriate to identify the phase transition that consists of studying the system in the supercritical regime, where the presence of a "loop" characterizes the first-order transition. All results obtained by the present approach are in full agreement with those obtained by using the constant rate ensemble, supporting that, in the thermodynamic limit the results from distinct ensembles are equivalent.