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:1305.6711 · Physical Review E

Phase transition in an exactly solvable reaction-diffusion process

Somayeh Zeraati, Farhad H Jafarpour, Haye Hinrichsen B

Published Jun 14, 2013• 9 citations
arXivHTMLPDF

Abstract

We study a nonconserved one-dimensional stochastic process which involves two species of particles $A$ and $B$. The particles diffuse asymmetrically and react in pairs as $A\ensuremath{\emptyset}\ensuremath{\leftrightarrow}AA\ensuremath{\leftrightarrow}BA\ensuremath{\leftrightarrow}A\ensuremath{\emptyset}$ and $B\ensuremath{\emptyset}\ensuremath{\leftrightarrow}BB\ensuremath{\leftrightarrow}AB\ensuremath{\leftrightarrow}B\ensuremath{\emptyset}$. We show that the stationary state of the model can be calculated exactly by using matrix product techniques. The model exhibits a phase transition at a particular point in the phase diagram which can be related to a condensation transition in a particular zero-range process. We determine the corresponding critical exponents and provide a heuristic explanation for the unusually strong corrections to scaling seen in the vicinity of the critical point.

DOI
10.1103/physreve.87.062120

Related papers

Phase separation transition in a nonconserved two-species model

Urna Basu

A one-dimensional stochastic exclusion process with two species of particles, + and -, is studied where density of each species can fluctuate but the total particle density is conserved. From the exact stationary state weights we show that, in the limiting case where density of negative particles vanishes, the system undergoes a phase separation transition where a macroscopic domain of vacancies form in front of a single surviving negative particle. We also show that the phase-separated state is associated with a diverging correlation length for any density and that the critical exponents characterizing the behavior in this region are different from those at the transition line. The static and the dynamical critical exponents are obtained from the exact solution and numerical simulations, respectively.

arXiv:1610.08898•Physical review. E•Dec 27, 2016•2 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

On matrix product ground states for reaction - diffusion models

Haye Hinrichsen B, Sven Sandow, Ingo Peschel

We discuss a new mechanism leading to a matrix product form for the stationary state of one-dimensional stochastic models. The corresponding algebra is quadratic and involves four different matrices. For the example of a coagulation - decoagulation model explicit four-dimensional representations are given and exact expressions for various physical quantities are recovered. We also find the general structure of n-point correlation functions at the phase transition.

arXiv:cond-mat/9510104•Journal of Physics A Mathematical and General•Jun 7, 1996•77 citations

Exact solution of a reaction-diffusion process with three-site interactions

Malte Henkel, Haye Hinrichsen B

The one-dimensional reaction diffusion process AA → A and A ∅ A → AAA is exactly solvable through the empty interval method if the diffusion rate equals the coagulation rate. Independently of the particle production rate, the model is always in the universality class of diffusion-annihilation. This allows us to check analytically the universality of finite-size scaling in a non-equilibrium critical point.

arXiv:cond-mat/0010062•Journal of Physics A Mathematical and General•Feb 19, 2001•48 citations

Coupled two-species model for the pair contact process with diffusion

Shengfeng Deng, Wei Li, Uwe Claus Täuber

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.

arXiv:2007.09496•Physical review. E•Oct 22, 2020•5 citations

Exact Analysis of Two-Species Totally Asymmetric Exclusion Process with Open Boundary Condition

Chikashi Arita

In recent years, one-dimensional driven diffusion models consisting of particles on a lattice have been studied. The most famous topics of these models are solvability and phase transitions. For example, the exact stationary-state solution of the asymmetric exclusion process (ASEP) with an open boundary condition can be obtained as a matrix product form and its physical quantities exhibit boundary-induced phase transitions in the thermodynamic limit. There are three regimes, namely, maximal-current (MC) phase, low-density (LD) phase and high-density (HD) phase in the phase diagram. The two-species models or three-state models are also interesting. The stationary-state solutions of the two-species totally asymmetric exclusion process (TASEP) with a periodic boundary condition and an open boundary condition can be obtained as matrix product forms. The model in ref. 6 exhibits two phases in which symmetry breaks spontaneously. Furthermore, the model in ref. 7 exhibits phase separation. In this note, the two-species TASEP with a new open boundary condition different from the condition in ref. 6 is introduced and we obtain the stationary-state solution exactly as a matrix product form and calculate the physical quantities in the thermodynamic limit. Let us consider a stochastic process on an L-site chain with two-species particles. We note j 1⁄4 0 if site j is empty, j 1⁄4 1 if site j is occupied by the first-class particle and j 1⁄4 2 if site j is occupied by the second-class particle. The exchange rule in the bulk of the chain is expressed as 10 ! 01 with a rate 1; ð1Þ 20 ! 02 with a rate 1; ð2Þ 12 ! 21 with a rate 1: ð3Þ

Journal of the Physical Society of Japan•Jun 8, 2006•16 citations
We study a nonconserved one-dimensional stochastic process which involves two species of particles AAA and BBB. The particles diffuse asymmetrically and react in pairs as A\ensuremath∅\ensuremath↔AA\ensuremath↔BA\ensuremath↔A\ensuremath∅A\ensuremath{\emptyset}\ensuremath{\leftrightarrow}AA\ensuremath{\leftrightarrow}BA\ensuremath{\leftrightarrow}A\ensuremath{\emptyset}A\ensuremath∅\ensuremath↔AA\ensuremath↔BA\ensuremath↔A\ensuremath∅ and B\ensuremath∅\ensuremath↔BB\ensuremath↔AB\ensuremath↔B\ensuremath∅B\ensuremath{\emptyset}\ensuremath{\leftrightarrow}BB\ensuremath{\leftrightarrow}AB\ensuremath{\leftrightarrow}B\ensuremath{\emptyset}B\ensuremath∅\ensuremath↔BB\ensuremath↔AB\ensuremath↔B\ensuremath∅. We show that the stationary state of the model can be calculated exactly by using matrix product techniques. The model exhibits a phase transition at a particular point in the phase diagram which can be related to a condensation transition in a particular zero-range process. We determine the corresponding critical exponents and provide a heuristic explanation for the unusually strong corrections to scaling seen in the vicinity of the critical point.
A one-dimensional stochastic exclusion process with two species of particles, + and -, is studied where density of each species can fluctuate but the total particle density is conserved. From the exact stationary state weights we show that, in the limiting case where density of negative particles vanishes, the system undergoes a phase separation transition where a macroscopic domain of vacancies form in front of a single surviving negative particle. We also show that the phase-separated state is associated with a diverging correlation length for any density and that the critical exponents characterizing the behavior in this region are different from those at the transition line. The static and the dynamical critical exponents are obtained from the exact solution and numerical simulations, respectively.
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.
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.
In recent years, one-dimensional driven diffusion models consisting of particles on a lattice have been studied. The most famous topics of these models are solvability and phase transitions. For example, the exact stationary-state solution of the asymmetric exclusion process (ASEP) with an open boundary condition can be obtained as a matrix product form and its physical quantities exhibit boundary-induced phase transitions in the thermodynamic limit. There are three regimes, namely, maximal-current (MC) phase, low-density (LD) phase and high-density (HD) phase in the phase diagram. The two-species models or three-state models are also interesting. The stationary-state solutions of the two-species totally asymmetric exclusion process (TASEP) with a periodic boundary condition and an open boundary condition can be obtained as matrix product forms. The model in ref. 6 exhibits two phases in which symmetry breaks spontaneously. Furthermore, the model in ref. 7 exhibits phase separation. In this note, the two-species TASEP with a new open boundary condition different from the condition in ref. 6 is introduced and we obtain the stationary-state solution exactly as a matrix product form and calculate the physical quantities in the thermodynamic limit. Let us consider a stochastic process on an L-site chain with two-species particles. We note j 1⁄4 0 if site j is empty, j 1⁄4 1 if site j is occupied by the first-class particle and j 1⁄4 2 if site j is occupied by the second-class particle. The exchange rule in the bulk of the chain is expressed as 10 ! 01 with a rate 1; ð1Þ 20 ! 02 with a rate 1; ð2Þ 12 ! 21 with a rate 1: ð3Þ