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

arXiv:cond-mat/0412749 · Physical Review E

Monte Carlo simulations of bosonic reaction-diffusion systems

Su‐Chan Park

Published Sep 9, 2005• 8 citations
arXivHTMLPDF

Abstract

An efficient Monte Carlo simulation method for bosonic reaction-diffusion systems which are mainly used in the renormalization group (RG) study is proposed. Using this method, one-dimensional bosonic single species annihilation model is studied and, in turn, the results are compared with RG calculations. The numerical data are consistent with RG predictions. As a second application, a bosonic variant of the pair contact process with diffusion (PCPD) is simulated and shown to share the critical behavior with the PCPD. The invariance under the Galilean transformation of this boson model is also checked and discussion about the invariance in conjunction with other models are in order.

DOI
10.1103/physreve.72.036111

Related papers

Driven Pair Contact Process with Diffusion

Su‐Chan Park, Hyunggyu Park

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.

arXiv:cond-mat/0406606•Physical Review Letters•Feb 17, 2005•27 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

Real-space renormalization for reaction-diffusion systems

Jef Hooyberghs, Carlo Vanderzande

The stationary state of stochastic processes such as reaction-diffusion systems can be related to the ground state of a suitably defined quantum Hamiltonian. Using this analogy, we investigate the applicability of a real-space renormalization group approach, originally developed for quantum spin systems, to interacting particle systems. We apply the technique to an exactly solvable reaction-diffusion system and to the contact process (both in d = 1). In the former case, several exact results are recovered. For the contact process, surprisingly good estimates of critical parameters are obtained from a small-cell renormalization.

Journal of Physics A Mathematical and General•Feb 1, 2000•27 citations

Universality class of the pair contact process with diffusion

Frank Smallenburg, G. T. Barkema

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.

arXiv:0710.3692•Physical Review E•Sep 23, 2008•16 citations

Critical properties of the reaction-diffusion model 2A→3A, 2A→0

Enrico Carlon, Malte Henkel, Ulrich Schollwöck

The steady-state phase diagram of the one-dimensional reaction-diffusion model $2\stackrel{\ensuremath{\rightarrow}}{A}3A,$ $2\stackrel{\ensuremath{\rightarrow}}{A}0$ is studied through the non-Hermitian density matrix renormalization group. In the absence of single-particle diffusion the model reduces to the pair-contact process, which has a phase transition in the universality class of directed percolation (DP) and an infinite number of absorbing steady states. When single-particle diffusion is added, the number of absorbing steady states is reduced to 2 and the model no longer shows DP critical behavior. The exponents $\ensuremath{\theta}={\ensuremath{\nu}}_{\ensuremath{\Vert}}/{\ensuremath{\nu}}_{\ensuremath{\perp}}$ and $\ensuremath{\beta}/{\ensuremath{\nu}}_{\ensuremath{\perp}}$ are calculated numerically. The value of $\ensuremath{\beta}/{\ensuremath{\nu}}_{\ensuremath{\perp}}$ is close to the value of the parity conserving universality class, in spite of the absence of local conservation laws.

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

Universality in the pair contact process with diffusion

G. T. Barkema, Enrico Carlon

The pair contact process with diffusion is studied by means of multispin Monte Carlo simulations and density matrix renormalization group calculations. Effective critical exponents are found to behave nonmonotonically as functions of time or of system length and extrapolate asymptotically towards values consistent with the directed percolation universality class. We argue that an intermediate regime exists where the effective critical dynamics resembles that of a parity conserving process.

arXiv:cond-mat/0302151•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Sep 18, 2003•28 citations
An efficient Monte Carlo simulation method for bosonic reaction-diffusion systems which are mainly used in the renormalization group (RG) study is proposed. Using this method, one-dimensional bosonic single species annihilation model is studied and, in turn, the results are compared with RG calculations. The numerical data are consistent with RG predictions. As a second application, a bosonic variant of the pair contact process with diffusion (PCPD) is simulated and shown to share the critical behavior with the PCPD. The invariance under the Galilean transformation of this boson model is also checked and discussion about the invariance in conjunction with other models are in order.
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.
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.
The stationary state of stochastic processes such as reaction-diffusion systems can be related to the ground state of a suitably defined quantum Hamiltonian. Using this analogy, we investigate the applicability of a real-space renormalization group approach, originally developed for quantum spin systems, to interacting particle systems. We apply the technique to an exactly solvable reaction-diffusion system and to the contact process (both in d = 1). In the former case, several exact results are recovered. For the contact process, surprisingly good estimates of critical parameters are obtained from a small-cell renormalization.
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.
The steady-state phase diagram of the one-dimensional reaction-diffusion model $2\stackrel{\ensuremath{\rightarrow}}{A}3A,$ $2\stackrel{\ensuremath{\rightarrow}}{A}0$ is studied through the non-Hermitian density matrix renormalization group. In the absence of single-particle diffusion the model reduces to the pair-contact process, which has a phase transition in the universality class of directed percolation (DP) and an infinite number of absorbing steady states. When single-particle diffusion is added, the number of absorbing steady states is reduced to 2 and the model no longer shows DP critical behavior. The exponents $\ensuremath{\theta}={\ensuremath{\nu}}_{\ensuremath{\Vert}}/{\ensuremath{\nu}}_{\ensuremath{\perp}}$ and $\ensuremath{\beta}/{\ensuremath{\nu}}_{\ensuremath{\perp}}$ are calculated numerically. The value of $\ensuremath{\beta}/{\ensuremath{\nu}}_{\ensuremath{\perp}}$ is close to the value of the parity conserving universality class, in spite of the absence of local conservation laws.