We introduce branching annihilating attracting walk (BAAW) in one dimension. The attracting walk is implemented by a biased hopping in such a way that a particle prefers hopping to a nearest neighbor located on the side where the nearest particle is found within the range of attraction. We study the BAAW with four offspring by extensive Monte Carlo simulation. At first, we find the critical exponents of the BAAW with infinite range of attraction, which are different from those of the directed Ising (DI) universality class. Our results are consistent with the recent observation [Daga and Ray, Phys. Rev. E 99, 032104 (2019)2470-004510.1103/PhysRevE.99.032104]. Then, by studying crossover behaviors, we show that as far as the range of attraction is finite the BAAW belongs to the DI class. We conclude that the origin of non-DI critical behavior of the BAAW with infinite range of attraction is the long-range nature of the attraction.
We introduce branching annihilating attracting walk (BAAW) in one dimension. The attracting walk is implemented by a biased hopping in such a way that a particle prefers hopping to a nearest neighbor located on the side where the nearest particle is found within the range of attraction. We study the BAAW with four offspring by extensive Monte Carlo simulation. At first, we find the critical exponents of the BAAW with infinite range of attraction, which are different from those of the directed Ising (DI) universality class. Our results are consistent with the recent observation [Daga and Ray, Phys. Rev. E 99, 032104 (2019)2470-004510.1103/PhysRevE.99.032104]. Then, by studying crossover behaviors, we show that as far as the range of attraction is finite the BAAW belongs to the DI class. We conclude that the origin of non-DI critical behavior of the BAAW with infinite range of attraction is the long-range nature of the attraction.
We introduce and numerically study the branching annihilating random walks with long-range attraction (BAWL). The long-range attraction makes hopping biased in such a manner that particle's hopping along the direction to the nearest particle has larger transition rate than hopping against the direction. Still, unlike the Lévy flight, a particle only hops to one of its nearest-neighbor sites. The strength of bias takes the form x^{-σ} with non-negative σ, where x is the distance to the nearest particle from a particle to hop. By extensive Monte Carlo simulations, we show that the critical decay exponent δ varies continuously with σ up to σ=1 and δ is the same as the critical decay exponent of the directed Ising (DI) universality class for σ≥1. Investigating the behavior of the density in the absorbing phase, we argue that σ=1 is indeed the threshold that separates the DI and non-DI critical behavior. We also show by Monte Carlo simulations that branching bias with symmetric hopping exhibits the same critical behavior as the BAWL.
Recently, Takayasu and Tretyakov [Phys. Rev. Lett. 68, 3060 (1992)] studied branching annihilating random walks with n=1-5 offsprings. These models exhibit a continuous phase transition to an absorbing state. Steady-state simulations yielded an estimate for the order parameter critical exponent \ensuremath{\beta} different from that of directed percolation. This result is quite surprising, as the universality class of directed percolation is known to be very robust. I have studied the critical behavior of the one-dimensional model with n=1 and 3 using time-dependent Monte Carlo simulations, and determined three critical exponents, all of which are in agreement with directed percolation.
We study the effects of hard core (HC) interactions between different species of particles on two-species branching annihilating random walks with one offspring(BAW2(1)). The single-species model belongs to the directed percolation (DP) universality class. In the BAW2(1) model, a particle creates one particle of the same species in its neighborhood with the probability σ(1−p) and of the different species with (1−σ)(1−p), where p is the hopping probability. Without HC interactions, this model always exhibits the DP-type absorbing transition for all σ. Even with HC interactions, the nature of the phase transitions does not change except near σ=0, where the HC interaction destabilizes and completely wipes away the absorbing phase. The model is always active except at the annihilation fixed point of zero branching rate (p=1). Critical behavior near the annihilation fixed point is characterized by exponents β=ν⊥=1/2 and ν∣∣=1.
Recently, Takayasu and Tretyakov [Phys. Rev. Lett. 68, 3060 (1992)] studied branching annihilating random walks (BAWs) with n=1--5 offspring. These models exhibit a continuous phase transition to an absorbing state. For odd n the models belong to the universality class of directed percolation. For even n the particle number is conserved modulo 2 and the critical behavior is not compatible with directed percolation. In this article I study the BAW with n=4 using time-dependent simulations and finite-size scaling, obtaining precise estimates for various critical exponents. The results are consistent with the conjecture: \ensuremath{\beta}/${\ensuremath{\nu}}_{\mathrm{\ensuremath{\perp}}}$=1/2, ${\ensuremath{\nu}}_{\mathrm{\ensuremath{\parallel}}}$/${\ensuremath{\nu}}_{\mathrm{\ensuremath{\perp}}}$=7/4, \ensuremath{\gamma}=0, \ensuremath{\delta}=2/7, \ensuremath{\eta}=0, and ${\mathrm{\ensuremath{\delta}}}_{\mathit{h}}$=9/2. These critical exponents characterize, respectively, the dependence of the order parameter (\ensuremath{\beta}/${\ensuremath{\nu}}_{\mathrm{\ensuremath{\perp}}}$) and relaxation time (${\ensuremath{\nu}}_{\mathrm{\ensuremath{\parallel}}}$/${\ensuremath{\nu}}_{\mathrm{\ensuremath{\perp}}}$) on system size, the growth of fluctuations (\ensuremath{\gamma}) close to the critical point, the long-time behavior of the probability of survival (\ensuremath{\delta}) and average number of particles (\ensuremath{\eta}) when starting at time zero with just two particles, and finally the decay of the order parameter (${\mathrm{\ensuremath{\delta}}}_{\mathit{h}}$) at the critical point in the presence of an external source.
With extensive Monte Carlo simulations, we present high-precision estimates of the critical exponents of branching annihilating random walks with two offspring, a prototypical model of the directed Ising universality class in one dimension. To estimate the exponents accurately, we propose a systematic method to find corrections to scaling whose leading behavior is supposed to take the form t −χ in the long-time limit at the critical point. Our study shows that χ ≈ 0.75 for the number of particles in defect simulations and χ ≈ 0.5 for other measured quantities, which should be compared with the widely used value of χ = 1. Using χ so obtained, we analyze the effective exponents to find that β/ν ‖ = 0.2872(2), z = 1.7415(5), η = 0.0000(2), and accordingly, β/ν ⊥ = 0.5000(6). Our numerical results for β/ν ‖ and z are clearly different from the conjectured rational numbers β/ν ‖ = $\tfrac{2} {7} $ ≈ 0.2857, z = $\tfrac{7} {4} $ = 1.75 by Jensen [Phys. Rev. E, 50, 3623 (1994)]. Our result for β/ν ⊥, however, is consistent with $\tfrac{1} {2} $ , which is believed to be exact.
The effect of blocking between different species occurring in one dimension is investigated here numerically in the case of particles following branching and annihilating random walk. It is shown that two-dimensional simulations confirm the field theoretical results with logarithmic corrections. In one dimension, however, if particles exhibit hard core interaction I confirm the very recent predictions of Kwon et al. [Phys. Rev. Lett. 85, 1682 (2000)] that there are two different universality classes depending on the spatial symmetry of the offspring production characterized by beta(S)=0.5 and beta(A)=2. Elaborate analysis of simulation data shows that the order parameter exponent beta does not depend on initial conditions or on diffusion rates of species but strong correction to scaling is observed. By systematic numerical simulations the critical point properties have been explored and initial condition dependence of the dynamical exponents Z and alpha is shown. In the case of a random initial state the particle-density decay at the critical point follows the t(-1/4) law with logarithmic corrections with two offsprings.