Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics

Critical behavior of branching annihilating random walks with an odd number of offsprings

Iwan Jensen

Published Jan 1, 1993• 66 citations
Publisher

Abstract

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.

DOI
10.1103/physreve.47.r1

Related papers

Critical exponents for branching annihilating random walks with an even number of offspring

Iwan Jensen

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.

arXiv:cond-mat/9405006•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Nov 1, 1994•140 citations

High-precision estimate of the critical exponents for the directed Ising universality class

Su-Chan Park

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.

arXiv:1406.7070•Journal of the Korean Physical Society•Feb 1, 2013•16 citations

Critical behavior in a model of random walkers in processes with conservation of parity, branching, multiple occupation and annihilation.

Fabiana Carvalho dos Anjos

In this dissertation we address the problem of second order phase transitions with absorbing states. We analyze the critical behavior of a one dimensional model simulated on one dimensional chains. The sites on the linear chain can be multiple occupied by particles that annhilates each other with a finite probability. It also presents parity-conserving branching with an even number of offsprings. The number of particles are constant modulo 2. We determine the threshold of the phase transition between the statistically stationary active state and the absorbing state. From steady-state simulations and a finite-size scaling analysis, we calculate the order-parameter, order-parameter fluctuations, and special correlation length critical exponents. We also followed the short-time critical relaxation to determine some dynamical critical exponents. The hyperscaling relation is checked for both sets of stationary and dynamical critical exponents which we found to be consistent with the BARW-PC universality class.

•Jul 27, 2016•0 citations

Critical branching-annihilating random walk of two species

Géza Ódor

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.

arXiv:cond-mat/0008381•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jan 26, 2001•16 citations

Pair contact process with diffusion: A new type of nonequilibrium critical behavior?

Haye Hinrichsen B

In the preceding article Carlon et al. investigate the critical behavior of the pair contact process with diffusion. Using density matrix renormalization group methods, they estimate the critical exponents, raising the possibility that the transition might belong to the same universality class as branching annihilating random walks with even numbers of offspring. This is surprising since the model does not have an explicit parity-conserving symmetry. In order to understand this contradiction, we estimate the critical exponents by Monte Carlo simulations. The results suggest that the transition might belong to a different universality class that has not been investigated before.

arXiv:cond-mat/0001177•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Feb 15, 2001•63 citations

DIRECTED PERCOLATION AND OTHER SYSTEMS WITH ABSORBING STATES: IMPACT OF BOUNDARIES

PER FRÖJDH, MARTIN HOWARD, K. B. Lauritsen

We review the critical behavior of nonequilibrium systems, such as directed percolation (DP) and branching-annihilating random walks (BARW), which possess phase transitions into absorbing states. After reviewing the bulk scaling behavior of these models, we devote the main part of this review to analyzing the impact of walls on their critical behavior. We discuss the possible boundary universality classes for the DP and BARW models, which can be described by a general scaling theory which allows for two independent surface exponents in addition to the bulk critical exponents. Above the upper critical dimension d c , we review the use of mean field theories, whereas in the regime d<d c , where fluctuations are important, we examine the application of field theoretic methods. Of particular interest is the situation in d=1, which has been extensively investigated using numerical simulations and series expansions. Although DP and BARW fit into the same scaling theory, they can still show very different surface behavior: for DP some exponents are degenerate, a property not shared with the BARW model. Moreover, a "hidden" duality symmetry of BARW in d=1 is broken by the boundary and this relates exponents and boundary conditions in an intricate way.

arXiv:cond-mat/0101387•International Journal of Modern Physics B•May 20, 2001•33 citations
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.
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.
In this dissertation we address the problem of second order phase transitions with absorbing states. We analyze the critical behavior of a one dimensional model simulated on one dimensional chains. The sites on the linear chain can be multiple occupied by particles that annhilates each other with a finite probability. It also presents parity-conserving branching with an even number of offsprings. The number of particles are constant modulo 2. We determine the threshold of the phase transition between the statistically stationary active state and the absorbing state. From steady-state simulations and a finite-size scaling analysis, we calculate the order-parameter, order-parameter fluctuations, and special correlation length critical exponents. We also followed the short-time critical relaxation to determine some dynamical critical exponents. The hyperscaling relation is checked for both sets of stationary and dynamical critical exponents which we found to be consistent with the BARW-PC universality class.
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.
In the preceding article Carlon et al. investigate the critical behavior of the pair contact process with diffusion. Using density matrix renormalization group methods, they estimate the critical exponents, raising the possibility that the transition might belong to the same universality class as branching annihilating random walks with even numbers of offspring. This is surprising since the model does not have an explicit parity-conserving symmetry. In order to understand this contradiction, we estimate the critical exponents by Monte Carlo simulations. The results suggest that the transition might belong to a different universality class that has not been investigated before.
We review the critical behavior of nonequilibrium systems, such as directed percolation (DP) and branching-annihilating random walks (BARW), which possess phase transitions into absorbing states. After reviewing the bulk scaling behavior of these models, we devote the main part of this review to analyzing the impact of walls on their critical behavior. We discuss the possible boundary universality classes for the DP and BARW models, which can be described by a general scaling theory which allows for two independent surface exponents in addition to the bulk critical exponents. Above the upper critical dimension d c , we review the use of mean field theories, whereas in the regime d<d c , where fluctuations are important, we examine the application of field theoretic methods. Of particular interest is the situation in d=1, which has been extensively investigated using numerical simulations and series expansions. Although DP and BARW fit into the same scaling theory, they can still show very different surface behavior: for DP some exponents are degenerate, a property not shared with the BARW model. Moreover, a "hidden" duality symmetry of BARW in d=1 is broken by the boundary and this relates exponents and boundary conditions in an intricate way.